PLUTO Test Problems  4.4-patch2
Functions
init.c File Reference

Disk-Planet interaction problem. More...

Functions

void InitDomain (Data *d, Grid *grid)
 

Detailed Description

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 $\propto T^{-1}$ 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:

The test can be carried out in polar (2D or 3D) or spherical (3D) coordinates and the following parameters determine the initial configuration:

  1. g_inputParam[Mstar]: controls the star mass (in solar masses)
  2. g_inputParam[Mdisk]: controls the disk mass (in solar masses)
  3. g_inputParam[Mplanet]: sets the planet mass (in earth masses)
  4. g_inputParam[Viscosity]: sets the amount of viscosity

Computation can be carried in the rotating or in the observer's frame of reference (ROTATING_FRAME to YES or NO, respectively). In particular:

hd_disk_planet.08.png
Density map for configuration #08 using AMR
Author
A. Mignone (migno.nosp@m.ne@t.nosp@m.o.inf.nosp@m.n.it)
Date
Jul 19, 2019

References:

Function Documentation

◆ InitDomain()

void InitDomain ( Data *  d,
Grid *  grid 
)

Assign initial condition by looping over the computational domain. Called after the usual Init() function to assign initial conditions on primitive variables. Value assigned here will overwrite those prescribed during Init().