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

Hydrostatic atmosphere in a smoothed gravitational potential. More...

Functions

void InitDomain (Data *d, Grid *grid)
 

Detailed Description

Set initial conditions for an hydrostatic atmosphere in cylindrical coordinates. Gravity is given as

\[ g = \left\{\begin{array}{ll} -1/R^2 & \quad{\rm for} \quad R > 1 \\ \noalign{\medskip} aR + bR^2 + cR^3 & \quad{\rm for} \quad R < 1 \end{array}\right.\,, \]

where $R = \sqrt{r^2 + z^2}$ is the spherical radius. The coefficients $a$, $b$ and $c$ are chosen to guarantee continuity of $g$, its first and second derivative (optionally). Density and pressure are tied by the isothermal condition $P = \rho/a$ so that the hydrostatic condition is

\[ \frac{1}{\rho} \frac{d\rho}{dr} = ag\,, \]

with the normalization $\rho = 1$ at $R = 1$.

The runtime parameters that are read from pluto.ini are

Configurations:

References:

Author
A. Mignone (migno.nosp@m.ne@t.nosp@m.o.inf.nosp@m.n.it)
Date
Jul 26, 2019

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().