PLUTO Test Problems
4.4-patch2
|
Bell instability setup. 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) |
Initialize fluid variables using Eqns. of [MVBM18], sect. 4.4. using a 1D configuration. The configuration is then rotated around the y- and z-axis in a way similar to Mignone, Tzeferacos & Bodo, JCP (2010). The initial phase is computed as (which is invariant under rotations), where
and
are the wavevector and coordinate vector in the rotated system, while
and
are the original (unrotated) 1D vectors.
The amount of rotation depends on the dimensionality of the problem and is uniquely specified by the domain size in the three directions Lx, Ly, Lz:
where and
. In such a way the modulus of the wavevector is
In order to have one wavelength, Ly and Lz must be chosen so that . Thus in 1D, 2D and 3D we choose:
Dim | Lx | Ly | Lz |
---|---|---|---|
1 | 1 | - | - |
2 | ![]() | ![]() | - |
3 | 3 | 3/2 | 3/2 |
The 6 configurations come in pairs (serial/parallel) and test the instability for eps = 0.1 in 1D (Conf. #1/#2), in 2D (Conf. #3/#4) and 3D (Conf. #5/#6).
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 |