PLUTO Test Problems
4.4-patch2
|
Relativistic magnetized Kelvin-Helmholtz problem. More...
Functions | |
void | InitDomain (Data *d, Grid *grid) |
Implements the initial condition for a relativistic KH flow in 2D Cartesian coordinates as in section 6.6 of Mignone et al. (ApJS, 2012). Density and pressure are initially constant and equal to and
where the latter value is retrieved from the definition of the Mach number
. The sound speed is
for the
IDEAL
equation of state and for the
TAUB
EoS (see Eq. [17], [31] and [32] of Mignone & McKinney [MNRAS, 2007]), where is the gas enthalpy.
The velocity is parallel to the x
axis with a hyperbolic tangent shear profile. Magnetic field has the form
where and
control the strength of the poloidal and toroidal component, respectively.
Periodic boundary conditions are assumed in the horizontal direction while zero-gradient (outflow) conditions hold at the lower and upper boundaries in the y
direction.
The presence of a magnetic field gives rise to a filamentary structure and to the formation of a number of vortices arranged symmetrically with respect to the central point. The figure below shows density for configuration #03 using AMR with 4 levels of refinement at t
= 5
The control parameters for this problem are:
g_inputParam[SIGMA_TOR]
: control the strength of the toroidal (z
) component of magnetic field;g_inputParam[SIGMA_POL]
: control the strength of the poloidal (x
) component of magnetic field;g_inputParam[VEL0]
: the flow velocity (in units of the speed of light);g_inputParam[MACH]
: the flow Mach number used to recover the pressure value;.
References