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

Sod shock tube problem. More...

Functions

void InitDomain (Data *d, Grid *grid)
 

Detailed Description

The Sod shock tube problem is one of the most used benchmark for shock-capturing schemes. It is a one-dimensional problem with initial condition given by a discontinuity separating two constant states:

\[ \begin{array}{lcll} \left(\rho,\, v_x,\, p\right)_L &=& \left(1, 0, 1\right) & \qquad\mathrm{for}\quad x < 0.5 \\ \noalign{\medskip} \left(\rho,\, v_x,\, p\right)_R &=& \left(\frac{1}{8}, 0, \frac{1}{10}\right) & \qquad\mathrm{for}\quad x > 0.5 \end{array} \]

The evolved structured at t=0.2 is shown in the panels below and consists of a left-going rarefaction wave, a right-going contact discontinutity and a right-going shock wave. The results shown here were carried with PARABOLIC interpolation, CHARACTERISIC_TRACING time stepping and the two_shock Riemann solver on 400 zones (configuration #04).

hd_sod.04.jpg
Flow profiles for the Sod shock tube at t = 0.2 using configuration #04

.

Author
A. Mignone (migno.nosp@m.ne@p.nosp@m.h.uni.nosp@m.to.i.nosp@m.t)
Date
June 08, 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().