Efficient, Nonlinear Second Moment Methods for Multigroup Thermal Radiative Transfer
This addresses computational bottlenecks in high-fidelity thermal radiative transfer simulations for applications like astrophysics or nuclear engineering, representing a methodological advance rather than a paradigm shift.
The authors tackled the computational challenges of multigroup thermal radiative transfer by developing an efficient nonlinear acceleration framework based on the Second Moment method, which achieved scalable linear solves via algebraic multigrid-preconditioned conjugate gradient by creating a symmetric positive definite low-order system.
Thermal radiative transfer (TRT) presents significant computational challenges due to the stiff, nonlinear coupling between radiation and material energy, particularly in multigroup, high-fidelity transport models. In this work, we develop an efficient nonlinear acceleration framework for TRT based on the Second Moment (SM) method. Our approach couples high-order discrete ordinates transport to a gray, diffusion-based low-order system that implicitly resolves the stiff absorption-emission physics, isolating this stiffness from the high-order system. The resulting algorithm alternates between transport sweeps and a Newton-type solution of the coupled low-order and material energy balance equations, utilizing nonlinear temperature elimination for improved robustness. Crucially, our approach is the first moment-based TRT algorithm with a symmetric and positive definite (SPD) low-order system enabling scalable linear solves via algebraic multigrid-preconditioned conjugate gradient. We investigate both consistent and independent low-order discretizations within a discontinuous Galerkin framework and assess their performance on one and two-dimensional gray and multigroup benchmark problems.