PLUTO Test Problems
4.4-patch2
|
Accretion disk in 3D spherical coordinates. More...
Functions | |
void | InitDomain (Data *d, Grid *grid) |
Setup a magnetized disk in 3D spherical coordinates ). The disk model is taken form [Flo11],[Mig12] and it consists of an initial equilibrium in a point mass gravity
with equilibrium profiles given by
where is the cylindrical radius while
and
is the sound speed. Here
H
is the scale height and it is proportional to the cylindrical radius. The constant of proportionality is the ratio defined by the input parameter
g_inputParam[H_R]
. With the current setting, one rotation period of the inner orbit (r=1
) is .
Differently from [Flo11], a polodial magnetic field is used here, with vector potential prescribed as follows:
where A0
is a constant chosen to prescribe a given plasma beta. The exponential terms on the right hand side confine the magnetic field in the midplane around R=6
. We also make the vector potential vanish identically for or
.
The boundary conditions are all set to userdef:
For the sake of simplicity only a quarter of a disk is considered. This test problems can be run normally or using the FARGO module, giving a saving factor of (roughly) 5. The figure below was produced by running the current setup at twice the resolution and for 90 orbits using FARGO.
References