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

Isentropic vortex problem. More...

Functions

void InitDomain (Data *d, Grid *grid)
 

Detailed Description

The isentropic vortex is a smooth exact solution of the 2D Euler equations. It consists of a single vortex centered at (5,5) in pressure equilibrium:

\[ \Big(\delta v_x\,, \delta v_y\Big) = - (y_c,\,x_c)\frac{\epsilon}{2\pi} \exp\left(\frac{1-r_c^2}{2}\right) \,,\quad T = \frac{p}{\rho} = 1 - \frac{(\Gamma-1)\epsilon^2}{8\Gamma\pi^2}\exp(1-r_c^2)\,,\quad s = \frac{p}{\rho^\Gamma} \]

where $ (x_c,\,y_c) = (x-5,\,y-5),\, r_c=\sqrt{x_c^2+y_c^2},\,\epsilon=5$ while $\Gamma = 1.4$ is the adiabatic index. A constant entropy is used with value s=1. The vortex shifts along the main diagonal of the computational domain with uniform velocity (1,1) and returns in its original position after t=10. The domain is assumed to be periodic and its size must be taken large enough to ensure there is no interaction between off-domain vortices.

The vortex problem is often used as computational benchmark to test the accuracy and dissipation of a numerical scheme in reproducing the vortex structure after several revolutions.

There're no input parameters for this problem.

hd_isentropic_vortex.13.jpg
Density cut at y=5 after 10 revolutions using PPM and Finite difference WENOZ and PPM (conf #01 and #03).
Author
A. Mignone (migno.nosp@m.ne@p.nosp@m.h.uni.nosp@m.to.i.nosp@m.t)
Date
Oct 3, 2014

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