PLUTO Test Problems
4.4-patch2
|
Radial & vertical drift test problem. More...
Functions | |
void | Init (double *v, double x1, double x2, double x3) |
void | InitDomain (Data *d, Grid *grid) |
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 | Spherical | Uniform | Rad osc. | Userdef |
#04 | Spherical | Uniform | Vert osc. | Userdef |
References
void Init | ( | double * | v, |
double | x1, | ||
double | x2, | ||
double | x3 | ||
) |
Here we use