PLUTO Test Problems
4.4-patch2
|
Disk-Planet interaction problem. More...
Functions | |
void | InitDomain (Data *d, Grid *grid) |
Simulate the interaction of a planet embedded in a disk as described in section 3.4 of Mignone et al., A&A (2012) 545, A152. This test is a nice benchmark for the FARGO module and the ROTATING_FRAME
switch. For testing-purposes no viscosity is used here. The initial condition consists of a locally isothermal configuration with temperature profile yielding a disk vertical height to radius of
0.05
. The gravitational potential due to the presence of the star and the planet is defined in BodyForcePotential() function.
The conventions used throught the implementation are the following:
r
= spherical radiusR
= cylindrical radiusz
= cylindrical heightth
= meridional angleThe test can be carried out in polar (2D or 3D) or spherical (3D) coordinates and the following parameters determine the initial configuration:
g_inputParam[Mstar]
: controls the star mass (in solar masses)g_inputParam[Mdisk]
: controls the disk mass (in solar masses)g_inputParam[Mplanet]
: sets the planet mass (in earth masses)g_inputParam[Viscosity]
: sets the amount of viscosityComputation can be carried in the rotating or in the observer's frame of reference (ROTATING_FRAME
to YES
or NO
, respectively). In particular:
References: