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 \]

An exponential 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 region.
  2. g_inputParam[RHO_IN]: density inside the initial region.
  3. g_inputParam[PRS_OUT]: ambient pressure
  4. g_inputParam[RHO_OUT]: ambient density
  5. g_inputParam[BMAG]: magnetic field strength
  6. g_inputParam[RADIUS]: radius of the initial over-pressurized region.
  7. g_inputParam[ETA]: resistivity.

We use the configuration from [MAR18] who employ Bx = 0.1 although similar configurations (with Bx=0.05) can be found also [Pal_etal09], [Dion_etal13], [Miz13]. 3D configurations can be found in [Kom07] (larger resistivity) and in [Dion_etal13] (smaller field).

Here we employ the following configurations:

Conf.DIM divB eta Res Notes Ref
#01 2 GLM 1.e-6 400^2 [MAR18, Kom07]
#02 2 CT 1.e-6 400^2 [MAR18, Kom07]
#03 3 GLM 1.e-6 128^3 [1] -
#04 3 CT 1.e-6 128^3 -

Note that oscillations appear in 2D when CFL > 0.4

  1. Only runs with CHARGE_SCHEME = 1 with large overshoots in the charge
Authors
A. Mignone (migno.nosp@m.ne@p.nosp@m.h.uni.nosp@m.to.i.nosp@m.t)
Date
Oct 29, 2018

References: