PLUTO Test Problems
4.4-patch2
|
RMHD rotor test problem in 2 or 3 dimensions. More...
Functions | |
void | InitDomain (Data *d, Grid *grid) |
The 2D rotor problem [dZBL03][Mig12] consists of a rapidly spinning disk embedded in a uniform background medium treaded by a costant magnetic field. The initial conditin reads as
where is the angular frequency of rotation. The computational domain is the unit square and outflow boundary conditions are imposed everywhere.
As the disk rotates, strong torsional Alfven waves form and propagate outward carrying angular momentum from the disk to the ambient. The emerging flow structure is enclosed by a circular fast forward shock traveling into the surrounding medium. An inward fast shock bounds the innermost oval region where density has been depleted to lower values. The presence of the magnetic field slows down the rotor, and the maximum Lorentz factor decreases from the nominal value of 10 to 2.2 (approx).
A list of tested configurations is given in the following table:
Conf. | GEOMETRY | DIM | T. STEPPING | RECONSTRUCTION | divB | AMR |
---|---|---|---|---|---|---|
#01 | CARTESIAN | 2 | RK3 | LINEAR | CT | NO |
#02 | CARTESIAN | 2 | RK2 | LINEAR | CT | NO |
#03 | CARTESIAN | 2 | HANCOCK | LINEAR | 8W | NO |
#04 | CARTESIAN | 2 | HANCOCK | LINEAR | CT | NO |
#05 | CARTESIAN | 3 | RK2 | LINEAR | GLM | NO |
#06 | CARTESIAN | 3 | RK2 | LINEAR | GLM | NO |
#07 | CARTESIAN | 3 | HANCOCK | LINEAR | CT | NO |
#08 | CARTESIAN | 3 | RK3 | LimO3 | GLM | NO |
The final solutin at t = 0.4 on a grid of 400x400 grid points is shown below.
The three dimensional version extends the current configuration to a spinning sphere and it is described in [MUB09].
References: