PLUTO Test Problems
4.4-patch2
|
Circularly polarized Alfven waves. More...
Functions | |
void | InitDomain (Data *d, Grid *grid) |
This is a rotated version of a 1D setup where
where is the translation velocity in the
direction,
is the phase (
in 1D and (
) in 2D), and
is the wave amplitude (
implies right going waves;
implies left going waves).
With this normalization, the Alfven velocity is always unity.
The configuration is rotated by specifying and
which express the ratios between the
- and
- components of the wave vector with the
component. In order to apply periodic boundary conditions everywhere, an integer number of wavelegnths must be contained in each direction, that is,
,
,
.
We use the tools in rotate.c which requires to specify the four integer shifts such that
The final time step is one period and is found from
–>
The runtime parameters that are read from pluto.ini
are
g_inputParam[EPS]
: sets the wave amplitude g_inputParam[VEL0]
: sets g_inputParam[PR0]
: sets the pressure of the 1D solution;g_inputParam[ALPHA_GLM]
: ;Configurations:
References: