A High-Order Finite Volume GENO Scheme with Implicit Time Integration for Three-Temperature Radiation Diffusion Equations
This work addresses computational efficiency and accuracy issues in radiation hydrodynamics simulations, which is incremental as it builds on existing finite volume and ENO methods.
The authors tackled the challenge of solving three-temperature radiation diffusion equations with large gradients and stiff terms by developing a high-order finite volume scheme with a central generalized ENO reconstruction and implicit time integration, achieving high accuracy, robustness, and large time-step capability in numerical tests.
This study presents a high-order finite volume scheme capable of large time-step integration for three-temperature radiation diffusion (3TRD) equations, where conservation is naturally achieved through energy update. To handle local large gradients and discontinuities in temperature, a central generalized ENO (GENO) reconstruction is developed for diffusion systems, which achieves essentially non-oscillatory reconstruction for discontinuous solutions. Compared to conventional nonlinear reconstruction methods, its most distinctive feature is the central-type symmetric sub-stencils, which ensure consistency between the numerics and the isotropic nature of thermal diffusion. Additionally, the central GENO method provides smooth states of temperature and temperature gradient at interfaces, facilitating the evaluation of numerical fluxes. Furthermore, interface flux evaluation for cases with discontinuous physical property parameters is modeled. To address the extremely small time-step issue caused by stiff diffusion and source terms, a dual-time-stepping method based on implicit time discretization is developed for the first time in 3TRD systems, with the advantage of decoupling temporal discretization from complex nonlinear spatial discretization. A series of numerical examples validates the high accuracy, physical property preservation, strong robustness, and large time-step integration capability of the present high-order central GENO scheme.