On Energy-Dependent Neutron Diffusion

While the energy-dependent neutron diffusion equation is widely employed in nuclear engineering, its status as an approximation to the transport equation is not yet completely understood, and several different approximations are in use to determine the diffusion coefficients. Past work on the theory underlying the diffusion approximation has often made use of asymptotic arguments; in the energy-dependent case, however, papers have appeared that differ substantially in their findings. Here we present a formal asymptotic derivation of the multigroup diffusion equation which addresses these differences, along with the varying and sometimes physically stringent assumptions employed in these works. Further, we show a way to exactly invert the relationship between flux and current in the P1 approximation, giving a matricial expression for the multigroup diffusion coefficient which is formally exact, has clear physical meaning, and which can be easily computed to arbitrary precision on the basis of cross-section data already produced by lattice calculations. The resulting 2-group diffusion coefficient for an infinite medium of hydrogen is calculated with Monte Carlo, and compared to the those deriving from the Cumulative Migration Method and from the out-scatter approximation.

💡 Research Summary

The paper addresses a long‑standing ambiguity in the theoretical foundation of the energy‑dependent multigroup neutron diffusion equation, which is widely used in nuclear engineering as an approximation to the full neutron transport equation. While many authors have derived diffusion equations from transport theory using asymptotic or variational methods, the energy‑dependent case has produced conflicting results, especially concerning the assumptions required for the diffusion coefficients.

The authors begin by revisiting the classic asymptotic scaling used in diffusion theory. They introduce a small spatial parameter ϵ that measures the ratio of a characteristic length L to the mean free path, and a time scaling based on the mean free lifetime T_g of each energy group. By multiplying the multigroup transport equations by the diagonal matrix Λ = diag(σ_t,g⁻¹) they obtain a dimensionless system in which the spatial derivative is O(ϵ) and the reaction terms are O(1). No explicit scaling of the time derivative to ϵ or ϵ² is performed, allowing the authors to discuss the necessity (or lack thereof) of the traditional infinite‑medium criticality condition k_∞ = 1.

An asymptotic expansion ψ = ψ⁰ + ϵψ¹ + ϵ²ψ² + … is then inserted. At zeroth order the angular flux is isotropic and satisfies a balance equation that is independent of space; this reproduces the familiar “infinite‑medium spectrum” solution and would require k_∞ = 1 if the time derivative were also scaled. At first order the P₀ and P₁ moments appear. The P₁ moment yields the familiar relation between the current J_g and the spatial gradient of the scalar flux Φ_g, but the authors deliberately drop the time derivative of J_g, arguing that for prompt‑subcritical reactors the current changes negligibly over a mean free time (T_g < 10⁻⁴ s). This leads to the standard approximate P₁ condition

  σ_t,g J_g – Σ_{g′} σ_s¹,g′→g J_{g′} = –(1/3)∇Φ_g .

Traditionally, this equation is inverted by invoking the “out‑scatter approximation,” which replaces the group‑coupled scattering term with a diagonal one, yielding the familiar diffusion coefficient D_g =