PLUTO Test Problems
4.4-patch2
|
Relativistic magnetized blast wave. More...
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.
An exponential smoothing is applied in the region . The input parameters used in this problem are:
g_inputParam[PRS_IN]
: pressure inside the initial region.g_inputParam[RHO_IN]
: density inside the initial region.g_inputParam[PRS_OUT]
: ambient pressureg_inputParam[RHO_OUT]
: ambient densityg_inputParam[BMAG]
: magnetic field strengthg_inputParam[RADIUS]
: radius of the initial over-pressurized region.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
CHARGE_SCHEME
= 1 with large overshoots in the chargeReferences: