Scribd
Upload
55 views

Topic 4: Computational Fluid Dynamics: Lecture 4-2: The Shock Tube Problem

The document discusses Sod's shock tube problem, a test of computational fluid dynamics algorithms. It involves the simulation of gas in a tube where the pressure and density are initially discontinuous in the middle. This setup results in contact discontinuities and shock waves that propagate through the gas. The document outlines the governing equations of gas dynamics and describes Sod's initial conditions for the shock tube problem, which has since been used to evaluate CFD methods for solving problems involving shock waves.

Uploaded by

rdeepak01
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PDF, TXT or read online on Scribd
55 views

Topic 4: Computational Fluid Dynamics: Lecture 4-2: The Shock Tube Problem

The document discusses Sod's shock tube problem, a test of computational fluid dynamics algorithms. It involves the simulation of gas in a tube where the pressure and density are initially discontinuous in the middle. This setup results in contact discontinuities and shock waves that propagate through the gas. The document outlines the governing equations of gas dynamics and describes Sod's initial conditions for the shock tube problem, which has since been used to evaluate CFD methods for solving problems involving shock waves.

Uploaded by

rdeepak01
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

Topic 4: Computational Fluid Dynamics

Lecture 4-2: The Shock Tube Problem

Wednesday, March 17, 2010

Contents 1 Sods Shock Tube Problem 1.1 The equations of gas dynamics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1.2 The shock tube problem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 2 4

1 SODS SHOCK TUBE PROBLEM

Sods Shock Tube Problem

A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 27, 1 (1978), to test the ability of various algorithms in solving uid dynamics problems with shock wave behavior.

1.1

The equations of gas dynamics

The equations of uid dynamics are mathematical statements of three fundamental physical principles: Mass is conserved F = ma, i.e., Newtons second law Energy is conserved The one-dimensional equations for the uid dynamics of a gas can be written in conservation form as follows:

1 SODS SHOCK TUBE PROBLEM

1.1 The equations of gas dynamics

(u) + = 0 t x (u) + (u2 + p) = 0 t x e + (u(e + p)) = 0 t x where is the density of the uid, u is the uid velocity, e is the energy per unit volume (length), and p is the pressure. We need one more equation to close the system. This is the equation of state 1 p = ( 1) e u2 2 where is the adiabatic gas index. For an ideal gas = 1.4. These equations can be written in vector form U + F(U) = 0 , t x where U = u , e u F = u2 + p . u(e + p)
3

1 SODS SHOCK TUBE PROBLEM

1.2 The shock tube problem

1.2

The shock tube problem

Sod considered a one-dimensional tube of unit length 0 x 1 and the following initial conditions at t = 0: (x, 0) = p(x, 0) = 1.0 for x
1 2 1 2

0.125 for x > 1.0 for x

1 2 1 2

0.1 for x > u(x, 0) = 0 .

This initial state can be produced by having a diaphragm in the middle of the tube. The gas to the left and right of the diaphragm is initially at rest. The pressure and density are discontinuous across the diaphragm. At t = 0, the diaphragm is broken. Two types of singularities then propagate through the gas: Contact discontinuities: The pressure p and velocity u are continuous, but the density and energy per unit volume e are discontinuous. Shock waves: All quantities p, u, and e are in general discontinuous across the shock front. To simulate a closed tube, reection boundary conditions can be applied at x = 0, 1. The shock tube then exhibits interesting behavior with shock waves and contact discontinuities bouncing back and forth in the tube and interacting with one another.
4

REFERENCES

REFERENCES

References [Recipes-C19-5] W.H. Press, S.A. Teukolsky, W. Vetterling, and B.P. Flannery, Numerical Recipes in C, Chapter 19 5: Relaxation Methods for Boundary Value Problems, http://www.nrbook.com/a/bookcpdf/c19-5.pdf.

You might also like