PLUTO Test Problems
4.4-patch2
|
Relativistic magnetized jet in axisymmetric coordinates. More...
Functions | |
void | InitDomain (Data *d, Grid *grid) |
The jet problem is set up in axisymmetric cylindrical coordinates as described in section 4.2.3 of Mignone, Ugliano & Bodo (2009, [MUB09] hereafter). We label ambient and jet values with the suffix "a" and "j", respectively. The ambient medium is at rest and is characterized by uniform density and pressure,
and
. The beam enters from the lower-z boundary from a circular nozzle of radius 1, carries a constant poloidal field
and a (radially-varying) toroidal component
. Flow variables are prescribed as
Here a
is the magnetization radius while . The pressure distribution is found by solving the radial momentum balance between thermal, centrifugal and magnetic forces. Neglecting rotation and assuming
Bz
to be constant the solution is given by
Here is the jet/ambient pressure at
r=1
and its value depends on the Mach number, see [MUB09].
The parameters controlling the problem is
g_inputParam[MACH]
: the jet Mach number g_inputParam[LORENTZ]
: the jet Lorentz factor;g_inpurParam[RHOJ]
: the jet density;g_inpurParam[SIGMA_POL]
: magnetization strength for poloidal magnetic field component (see [MUB09], Eq [65]);g_inpurParam[SIGMA_TOR]
: magnetization strength for toroidal magnetic field component (see [MUB09], Eq [65]);The different configurations are:
The following image show the (log) density map at the end of simulation for setup #01.
References