ANS is committed to advancing, fostering, and promoting the development and application of nuclear sciences and technologies to benefit society.
Explore the many uses for nuclear science and its impact on energy, the environment, healthcare, food, and more.
Explore membership for yourself or for your organization.
Conference Spotlight
2026 Nuclear Energy Conference & Expo (NECX)
August 24–27, 2026
Dallas, TX|Hilton Anatole
Latest Magazine Issues
Jun 2026
Jan 2026
2026
Latest Journal Issues
Nuclear Science and Engineering
July 2026
Nuclear Technology
June 2026
Fusion Science and Technology
May 2026
Latest News
Spent fuel recycling and conditioning topic of U.S.-Japan meeting
Officials with the Department of Energy’s Office of Environmental Management discussed spent nuclear fuel recycling and conditioning with counterparts from Japan during the 13th U.S.-Japan Technical Meeting of the Civil Nuclear Energy Research and Development Working Group, held recently in Santa Fe, N.M.
Hans R. Hammer, Jim E. Morel, Yaqi Wang
Nuclear Science and Engineering | Volume 193 | Number 4 | April 2019 | Pages 388-403
Technical Paper | doi.org/10.1080/00295639.2018.1525977
Articles are hosted by Taylor and Francis Online.
Second-order forms of the transport equation allow the use of continuous finite elements (CFEMs). This can be desired in multiphysics calculations where other physics require CFEM discretizations. Second-order transport operators are generally self-adjoint, yielding symmetric positive-definite (SPD) matrices, which allow the use of efficient linear algebra solvers with an enormous advantage in memory usage.
Least-squares (LS) forms of the transport equation can circumvent the void problems of other second-order forms but are almost always nonconservative. Additionally, the standard LS form is not compatible with discrete ordinates method (SN) iterative solution techniques such as source iteration. A new form of the LS transport equation has recently been developed that is compatible with voids and standard SN iterative solution techniques. Performing nonlinear diffusion acceleration (NDA) using an independently differenced low-order equation enforces conservation for the whole system and makes this equation suitable for reactor physics calculations. In this context, “independent” means that both the transport and low-order solutions converge to the same scalar flux and current as the spatial mesh is refined, but for a given mesh, the solutions are not necessarily equal.
In this paper we show that introducing a weight function into this LS equation improves issues with causality and can render our equation equal to the self-adjoint angular flux (SAAF) equation. Causality is a principle of the transport equation that states that information travels only downstream along characteristics. This principle can be violated numerically. We show how to limit the weight function in voids and demonstrate the effect of this limit on accuracy. Using the C5G7 benchmark, we compare our method to the SAAF formulation with a void treatment (SAAFτ) that is not self-adjoint and has a nonsymmetric coefficient matrix. We show that the weighted LS equation with NDA gives acceptable accuracy relative to the SAAFτ equation while maintaining a SPD system matrix.
