A direct Eulerian GRP scheme for radiation hydrodynamical equations in diffusion limit
This work provides a numerical method for solving radiation hydrodynamical equations, which is important for astrophysical and high-energy density physics simulations.
The paper develops a second-order accurate direct Eulerian GRP scheme for radiation hydrodynamical equations in the zero diffusion limit, achieving second-order accuracy and high resolution of strong discontinuities as demonstrated by numerical examples.
The paper proposes a second-order accurate direct Eulerian generalized Riemann problem (GRP) scheme for the radiation hydrodynamical equations (RHE) in the zero diffusion limit. The difficulty comes from no explicit expression of the flux in terms of the conservative vector. The characteristic fields and the relations between the left and right states across the elementary-waves are first studied, and then the solution of the one-dimensional Riemann problem is analyzed and given. Based on those, the direct Eulerian GRP scheme is derived by directly using the generalized Riemann invariants and the Runkine-Hugoniot jump conditions to analytically resolve the left and right nonlinear waves of the local GRP in the Eulerian formulation. Several numerical examples show that the GRP scheme can achieve second-order accuracy and high resolution of strong discontinuity.