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

Relativistic shock tube problems. More...

Functions

void InitDomain (Data *d, Grid *grid)
 

Detailed Description

Shock tube problems are widely used for code benchmarking since they directly probe the solver's ability in resolving elementary waves and their propagation. In addition, numerical results can be compared with analytical solutions.

Here we provide a simple suite of shock tube for the relativistic hydro module, by following Mignone & Bodo (2005).

\[ \left(\rho,\, v_x,\, p\right) = \left\{\begin{array}{ll} \left(\rho,\, v_x,\, p\right)_L & \quad\mathrm{for}\quad x < 0.5 \\ \noalign{\medskip} \left(\rho,\, v_x,\, p\right)_R & \quad\mathrm{otherwise} \end{array}\right. \]

The first four configurations correspond to Problem 1-4 of Mignone & Bodo (2005):

Test rho vx p Gamma
1L 1 0.9 1 4/3
1R 1 0.0 10 -
2L 1 -0.6 10 5/3
2R 10 0.5 20 -
3L 10 -0.6 40/35/3
3R 1 0.0 0 -
4L 1 0.0 1.e35/3
4R 1 0.0 1.e-2-

Conf #05 and #06 are equivalent to the first two, using the entropy equation. Results are shown in the four figures below:

rhd_shock_tube.01.jpg
Flow profiles for the 1st shock tube problem
rhd_shock_tube.02.jpg
Flow profiles for the 2nd shock tube problem
rhd_shock_tube.03.jpg
Flow profiles for the 3rd shock tube problem
rhd_shock_tube.04.jpg
Flow profiles for the 4th shock tube problem
Author
A. Mignone (migno.nosp@m.ne@t.nosp@m.o.inf.nosp@m.n.it)
Date
May 03, 2020

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