PLUTO Test Problems
4.4-patch2
|
MHD Rotor test problem. More...
Functions | |
void | InitDomain (Data *d, Grid *grid) |
void | UserDefBoundary (const Data *d, RBox *box, int side, Grid *grid) |
The rotor problem consists of a rapidly spinning cylinder embedded in a static background medium with uniform density and pressure and a constant magnetic field along the x direction . The cylinder rotates uniformly with constant angular velocity
and has larger density:
Here and
is a taper function. The ideal equation of state with
is used. As the disk rotates, strong torsional Alfven waves form and propagate outward carrying angular momentum from the disk to the ambient.
A list of tested configurations is given in the following table:
Conf. | GEOMETRY | divB | BCK_FIELD | AMR |
---|---|---|---|---|
#01 | CARTESIAN | CT | NO | NO |
#02 | POLAR | CT | NO | NO |
#03 | POLAR | 8W | NO | NO |
#04 | POLAR | CT | YES | NO |
#05 | POLAR | GLM | NO | NO |
#06 | CARTESIAN | GLM | NO | NO |
#07 | CARTESIAN | GLM | NO | YES |
#08 | CARTESIAN | 8W | NO | YES |
#09 | POLAR | CT | YES | NO |
#10 | POLAR | CT | YES | NO |
#11 | POLAR | GLM | YES | YES |
A snapshot of the solution using static and AMR grid is given below.
Reference:
void InitDomain | ( | Data * | d, |
Grid * | grid | ||
) |
void UserDefBoundary | ( | const Data * | d, |
RBox * | box, | ||
int | side, | ||
Grid * | grid | ||
) |
Provide inner radial boundary condition in polar geometry. Zero gradient is prescribed on density, pressure and magnetic field. For the velocity, zero gradient is imposed on v/r (v = vr, vphi).