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

Radial & vertical drift test problem. More...

Functions

void Init (double *v, double x1, double x2, double x3)
 
void InitDomain (Data *d, Grid *grid)
 

Detailed Description

Implement the equilibrium solution of Nelson, Gressel & Umurhan (2013). Units are chosen such that GM = 1, R0 = 1.

The indices p and q are user-supplied. The parameters rho0 and T0 are the density and temperature at the fiducial radius R0.

In the program output, particles positions and velocity are written to disk together with the analytical solution, obtained by integrating the equations of motion with an ODE solver.

Config. Geometry Grid Particles Stopping time
#01 Cyl Uniform Rad osc. Userdef
#02 Cyl Uniform Vert osc. Userdef
#03 SphericalUniform Rad osc. Userdef
#04 SphericalUniform Vert osc. Userdef
Author
A. Mignone (migno.nosp@m.ne@t.nosp@m.o.inf.nosp@m.n.it) M. Flock
Date
March 30, 2019

References

Function Documentation

◆ Init()

void Init ( double *  v,
double  x1,
double  x2,
double  x3 
)

Here we use

  • r = spherical radius
  • R = cylindrical radius

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