PLUTO Test Problems  4.4-patch2
init.c File Reference

Radiative magnetized blast wave. More...

Detailed Description

Sets up a 2D radiative magnetized blast wave. Initial conditions are set such that $p_g\sim\mathbf{B}^2/2<E_r/3$, so that matter dynamics is magnetically dominated when the opacities are low, and radiation-dominated in the opposite case. The latter case also serves to investigate the high-absorption regime in which both the diffusion approximation and LTE are valid.

Gas pressure and density are initially set as follows:

\[ \left(\begin{array}{c} p \\ \rho \end{array}\right) = \left(\begin{array}{c} p_1 \\ \rho_1 \end{array}\right) \delta + \left(\begin{array}{c} p_0 \\ \rho_0 \end{array}\right) (1-\delta) \]

where $p_0$ and $\rho_0$ are the ambient values while $p_1$ and $\rho_1$ identify the over-pressurized region. Here $R=\sqrt{x^2 + y^2}$ is the cylindrical radius while $\delta\equiv\delta(R/R_0)$ is a taper function that decreases linearly for $R_0<R\le1$ (we use $R_0=0.8$). The radiation field is initially introduced in equilibrium with the gas, with $\mathbf{F}_r=\mathbf{0}$. A null scattering opacity is set for every configuration.

The following setups are implemented: