A radiation-hydrodynamics scheme valid from the transport to the diffusion limit

E. Audit,P. Charrier,J.-P. Chièze,B. Dubroca

DOI: https://doi.org/10.48550/arXiv.astro-ph/0206281

2002-06-17

Abstract:We present in this paper the numerical treatment of the coupling between hydrodynamics and radiative transfer. The fluid is modeled by classical conservation laws (mass, momentum and energy) and the radiation by the grey moment $M_1$ system. The scheme introduced is able to compute accurate numerical solution over a broad class of regimes from the transport to the diffusive limits. We propose an asymptotic preserving modification of the HLLE scheme in order to treat correctly the diffusion limit. Several numerical results are presented, which show that this approach is robust and have the correct behavior in both the diffusive and free-streaming limits. In the last numerical example we test this approach on a complex physical case by considering the collapse of a gas cloud leading to a proto-stellar structure which, among other features, exhibits very steep opacity gradients.

Astrophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to accurately perform numerical calculations under different physical conditions from the transport limit to the diffusion limit in radiation hydrodynamics. Specifically, the authors focus on how to develop a mathematical model and numerical method that can handle these different physical conditions when the material opacity values vary greatly. These problems are very common in fields such as astrophysics, laser fusion, and plasma physics, but traditional numerical methods often struggle to maintain high precision and stability under different physical conditions.

### Main contributions of the paper:
1. **Model and method**: The authors propose a new model and numerical scheme that can describe radiation hydrodynamics under a wide range of physical conditions. This model describes the fluid based on the classical conservation laws (mass, momentum, and energy) and uses the gray - moment M1 system to describe radiation.
2. **Numerical scheme**: To handle the diffusion limit, the authors made an asymptotically - preserving modification to the HLLE scheme. This modification enables the numerical scheme to exhibit correct behavior in both the transport limit and the diffusion limit.
3. **Numerical experiments**: The paper presents the results of multiple numerical experiments, verifying the robustness and accuracy of the proposed model and numerical scheme under different physical conditions. In particular, the authors tested the performance of this method in complex physical situations by simulating the process of gas cloud collapse to form a protostellar structure.

### Key equations and concepts:
- **Fluid conservation equations**:
\[
\begin{aligned}
&\frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{u}) = 0\\
&\frac{\partial (\rho\mathbf{u})}{\partial t}+\nabla\cdot(\rho\mathbf{u}\otimes\mathbf{u}+P\mathbf{I})=\kappa\rho\mathbf{F}_r+\mathbf{F}\\
&\frac{\partial E}{\partial t}+\nabla\cdot(\mathbf{u}(E + P))=-\kappa\rho c(a_rT^4 - E_r)+(\kappa\rho\mathbf{F}_r+\mathbf{F})\cdot\mathbf{u}
\end{aligned}
\]

- **Radiation transport equations**:
\[
\begin{aligned}
&\frac{\partial E_r}{\partial t}+\nabla\cdot(E_r\mathbf{u})+\nabla\cdot\mathbf{F}_r+P_r:\nabla\mathbf{u}=\kappa\rho c(a_rT^4 - E_r)\\
&\frac{\partial \mathbf{F}_r}{\partial t}+\nabla\cdot(\mathbf{F}_r\mathbf{u})+c^2\nabla\cdot\mathbf{P}_r+(\mathbf{F}_r\cdot\nabla)\mathbf{u}=-\kappa\rho c\mathbf{F}_r
\end{aligned}
\]

- **M1 closure relation**:
\[
\mathbf{P}_r = D E_r
\]
where the Eddington tensor \(D\) is defined as:
\[
D=\frac{1-\chi}{2}\mathbf{I}+\frac{3\chi - 1}{2}\mathbf{n}\otimes\mathbf{n}
\]