PLUTO Test Problems  4.4-patch2
init.c File Reference

Relativistic magnetized blast wave. More...

Detailed Description

Set the initial condition for the relativistic magnetized blast wave problem in 2D or 3D. It consists of a highly pressurized region inside a circle (in 2D) or a sphere (in 3D) embeddd in a static uniform medium with lower pressure. The magnetic field is constant and threads the whole computational domain.

\[ (\rho,\, p) = \left\{\begin{array}{lcl} (\rho_{\rm in},\, p_{\rm in}) & \quad\mathrm{or}\quad & r < r_c \\ \noalign{\medskip} (\rho_{\rm out},\, p_{\rm out}) & \quad\mathrm{or}\quad & r \ge r_c \end{array}\right. \,,\qquad |\vec{B}| = B_0 \]

In 3D, a linear smoothing is applied in the region $ r_c<r<1$. The input parameters used in this problem are:

  1. g_inputParam[PRS_IN]: pressure inside the initial circular (2D) or spherical (3D) region.
  2. g_inputParam[PRS_OUT]: ambient pressure
  3. g_inputParam[RHO_OUT]: ambient density
  4. g_inputParam[BMAG]: magnetic field intensity
  5. g_inputParam[THETA]: angle between mag. field and z-axis (Cartesian only )
  6. g_inputParam[PHI]: angle between mag. field and xy-plane (Cartesian only)
  7. g_inputParam[RADIUS]: radius of the initial over-pressurized region.

Note that a given choice of parameters can be re-scaled by an arbitrary factor $\eta$ by letting $ \{\rho,\, p\} \to \eta^2\{\rho,\,p\},\, B \to \eta B $.

The different configurations are given in the following table

Conf.rho, p (in) rho, p (out) B Dim Ref
#01 1, 10^3 1.0, 0.01 4.0 2 [dZBL2003], Sec 4.2
#02 1, 10^3 1.0, 0.01 4.0 2 [dZBL2003], Sec 4.2
#03 0.01, 1 1.e-4, 3.e-5 1.0 3 [Mig_etal2007], sec 5.7
#04 0.01, 1 1.e-4, 3.e-5 1.0 2(ax)[Mig_etal2007, MB2006]
#05 0.01, 1 1.e-4, 5.e-3 1.0 2 [BS2011], Sec 4.6
#06 0.01, 1 1.e-4, 5.e-3 1.0 2 [BS2011]
#07 0.01, 1 1.e-4, 5.e-3 1.0 2 [BS2011]
#08 0.01, 1 1.e-4, 5.e-4 0.1 2 [Leis_etal2005, dZZBL2007]

Some of the strongly magnetized configurations can pass this test only by taking some precautions (e.g. correcting total energy with staggered magnetic field).

rmhd_blast.02.jpg
Density map (in log scale) for configuration #02

References

Authors
A. Mignone (migno.nosp@m.ne@t.nosp@m.o.inf.nosp@m.n.it)
Date
Jan 03, 2020