Jurnal Solusi Persamaan Difusi Neutron Satu Grup Satu Dimensi

𝕸𝖆𝖘 𝕬𝖓𝖉𝖞 𝕯𝖎𝖌𝖎𝖙𝖆𝖑 𝕸𝖊𝖉𝖎𝖆

3 min read

·

Apr 27, 2025

61

Share

A Complete Recipe for Solving the One-Group, One-Dimensional Neutron Diffusion Equation

The one-group, one-dimensional neutron diffusion equation is a fundamental equation in reactor physics, used to model the neutron flux distribution within a nuclear reactor. Solving this equation is crucial for reactor design and safety analysis. This article provides a comprehensive recipe for solving this equation, encompassing both analytical and numerical methods.

Understanding the Equation

The time-independent, one-group, one-dimensional neutron diffusion equation is expressed as:

∇²Φ(x) + (k_∞ - 1)/M² * Φ(x) = 0

where:

• Φ(x) is the neutron flux as a function of position x.
• k_∞ is the infinite multiplication factor, representing the neutron production rate.
• M² is the migration area, characterizing the average distance a neutron travels before absorption or leakage. The migration area incorporates diffusion and slowing-down properties of the reactor core. It is usually expressed as M² = L² + τ, where L² is the diffusion length and τ is the slowing down length.

The specific form of the Laplacian (∇²) depends on the geometry:

• Slab Geometry: ∂²Φ(x)/∂x²
• Cylindrical Geometry: (1/r) * ∂/∂r (r * ∂Φ(r)/∂r)
• Spherical Geometry: (1/r²) * ∂/∂r (r² * ∂Φ(r)/∂r)

Boundary Conditions

Accurate solutions require appropriate boundary conditions. Common boundary conditions include:

• Zero Flux Boundary Condition: Φ(x) = 0 at the boundaries. This represents a highly absorbing boundary, such as a strong neutron absorber material.
• Extrapolated Boundary Condition: Φ(x) = 0 at a distance of approximately 0.71 * λ_tr beyond the physical boundary, where λ_tr is the transport mean free path. This condition accounts for the fact that neutrons don't abruptly stop at a material boundary.

Analytical Solution Methods

For simple geometries and material properties, analytical solutions are possible. These often involve solving ordinary differential equations. The specific solution method depends on the geometry and boundary conditions. Techniques may include separation of variables and application of appropriate eigenfunctions.

For instance, in slab geometry with zero flux boundary conditions, the solution involves trigonometric functions and determining the eigenvalues that satisfy the boundary conditions. This often leads to transcendental equations that may require numerical methods for solving.

Numerical Solution Methods

For complex geometries or heterogeneous materials, numerical methods are essential. Popular numerical techniques include:

• Finite Difference Method (FDM): This method discretizes the spatial domain and approximates the derivatives using difference quotients. This leads to a system of linear algebraic equations which can be solved using matrix methods.
• Finite Element Method (FEM): This method divides the domain into smaller elements and uses basis functions to approximate the neutron flux within each element. This method is particularly suitable for complex geometries.

Selection of the appropriate numerical method depends on the geometry, material properties, and desired accuracy. Each method offers trade-offs in terms of computational cost and accuracy.

Iterative Methods

Solving the resulting linear system, especially in large problems, often requires iterative methods like:

• Gauss-Seidel Iteration: An iterative method that updates the solution iteratively based on the previous iterations.
• Successive Over-Relaxation (SOR): Improves the convergence rate of Gauss-Seidel iteration.
• Other advanced iterative solvers.

The choice of an iterative method is important to ensure efficient convergence.

Post-Processing and Interpretation

Once the neutron flux distribution is obtained (either analytically or numerically), further analysis can be performed, including:

• Calculating the effective multiplication factor (k_eff): This indicates whether the reactor is critical (k_eff = 1), subcritical (k_eff < 1), or supercritical (k_eff > 1).
• Determining the power distribution: This informs the design and operational safety of the reactor.
• Analyzing the neutron flux shape: Understanding this helps in optimizing fuel loading and control rod placement.

This detailed recipe provides a comprehensive framework for solving the one-group, one-dimensional neutron diffusion equation. Choosing the appropriate method depends heavily on the specifics of the problem at hand. Remember to carefully consider the boundary conditions and select suitable numerical techniques for efficient and accurate solutions. The results of these calculations are critical for reactor design, analysis, and safety.