PLUTO Test Problems  4.4-patch2
Functions
init.c File Reference

Relativistic magnetized blast wave. More...

Functions

void InitDomain (Data *d, Grid *grid)
 

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:

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
Authors
A. Mignone (migno.nosp@m.ne@t.nosp@m.o.inf.nosp@m.n.it)
Date
Aug 09, 2019

Function Documentation

◆ InitDomain()

void InitDomain ( Data *  d,
Grid *  grid 
)

Assign initial condition by looping over the computational domain. Called after the usual Init() function to assign initial conditions on primitive variables. Value assigned here will overwrite those prescribed during Init().