PLUTO Test Problems
4.4-patch2
|
Planar shock wave problem with particles. More...
Functions | |
void | InitDomain (Data *d, Grid *grid) |
void | Analysis (const Data *d, Grid *grid) |
void | BackgroundField (double x1, double x2, double x3, double *B0) |
void | UserDefBoundary (const Data *d, RBox *box, int side, Grid *grid) |
Set the initial condition for a Planar shocks wave problem. Ambient unshocked density and pressure are constant everywhere and equal to rho0
and Pg
. While rho0
is an input parameter, Pg=1.e-4
(for MHD runs) and Pg
= 1.0e-8 (for RMHD runs). Also for RMHD runs, a transverse velocity of 0.01 is set initially in the Ambient medium.
The shocked gas (x < 1.0) is allowed to move from left to right and the fluid quantities in this shocked region are given as input Parameters. The fluid quantities in the shocked region are estimated by solving the Rankine-Hugoniot conditions for specified shock velocities.
The input parameters read from pluto.ini are labeled as:
g_inputParam[RHOB]
: density in shocked gasg_inputParam[PGB]
: pressure in shocked gas;g_inputParam[VXB]/tt>: X-component velocity in shocked gas.
g_inputParam[VYB]/tt>: Y-component velocity in shocked gas.
g_inputParam[BYB]/tt>: Y-component mag. field in shocked gas.
g_inputParam[BETA]/tt>: Plasma beta.
g_inputParam[THETA_B]/tt>: Angle between BVec and X-axis in degrees.
The available configuration refer to: -#01 Cartesian (2.5D, MHD, Ideal, THETA_B = 30.0, BETA = 100.) -#02 Cartesian (2.5D, RMHD, Taub, THETA_B = 3.0, BETA = 0.01) -#03 Cartesian (2.5D, RMHD, Taub, THETA_B = 83.0, BETA = 0.01) -#04 Same as #02, but with a perturbation in the upstream.
References:
void Analysis | ( | const Data * | d, |
Grid * | grid | ||
) |
Perform runtime data analysis.
[in] | d | the PLUTO Data structure |
[in] | grid | pointer to array of Grid structures |
void BackgroundField | ( | double | x1, |
double | x2, | ||
double | x3, | ||
double * | B0 | ||
) |
Define the component of a static, curl-free background magnetic field.
[in] | x1 | position in the 1st coordinate direction ![]() |
[in] | x2 | position in the 2nd coordinate direction ![]() |
[in] | x3 | position in the 3rd coordinate direction ![]() |
[out] | B0 | array containing the vector componens of the background magnetic field |
void InitDomain | ( | Data * | d, |
Grid * | grid | ||
) |
void UserDefBoundary | ( | const Data * | d, |
RBox * | box, | ||
int | side, | ||
Grid * | grid | ||
) |
Assign user-defined boundary conditions.
[in,out] | d | pointer to the PLUTO data structure containing cell-centered primitive quantities (d->Vc) and staggered magnetic fields (d->Vs, when used) to be filled. |
[in] | box | pointer to a RBox structure containing the lower and upper indices of the ghost zone-centers/nodes or edges at which data values should be assigned. |
[in] | side | specifies the boundary side where ghost zones need to be filled. It can assume the following pre-definite values: X1_BEG, X1_END, X2_BEG, X2_END, X3_BEG, X3_END. The special value side == 0 is used to control a region inside the computational domain. |
[in] | grid | pointer to an array of Grid structures. |