VM-EXD-MECH-014

VM-EXD-MECH-014
Propagation of Shock and Rarefaction Waves in a Shock Tube

Overview

Reference: Harlow, Francis H. et al., Fluid Dynamics – A LASL Monograph, LA-4700, June 1971.
Analysis Type(s): Explicit Dynamics 3-D
Elements: Multi-Material Euler Quad
Boundary Conditions: Default Wall
Structural Interactions: No
Fluid-Structure Interactions: No
Bonds: No
Materials: Ideal Gas

Test Case

A 1000 mm (1000 element) shock tube is divided into two sections by a diaphragm at its midpoint. Initially, gas is at rest on both sides. To the left of the diaphragm, the gas is initially at a higher density and pressure than on the right . Both sides of the diaphragm have the same internal energy .

At t = 0, the diaphragm is removed and the gas in the tube is observed. A shock is observed moving from left to right, as well as a contact discontinuity moving from the left to the right, and a rarefaction wave moving from the right to the left.

There is no significant length to the system. The appearance of the system at a later time is a magnification of an earlier appearance (similarity solution).

Analysis Assumptions and Modeling Notes

Two interacting parts are created, one part is assigned a material with Euler inflow properties and the other ideal gas properties. The parts are then made to interact The initial conditions for the ideal gas in both shock tubes are:

γ = 1.4

ρ- = 0.01 gm/ccρ+ = 0.001 gm/cc
p- = 1000 MPap+ = 100 MPa
I0 = 2.5E5 J/Kgλ = ρ-/ρ+ = p-/p+ = 10

The calculation is run to t = 0.8 ms, the time at which the shock front nears the right boundary.

Results Comparison

Applying the theoretical solution to the initial conditions used for the simulation, the following values are obtained:

Final pressure in tube: MPa
Density of gas on left of diaphragm: g/mL
Density of gas on right of diaphragm: g/mL
Continuous pressure: MPa

Pressure and density profiles along the length of the shock tube at the end of the calculation show that the theoretical values are accurately predicted. The graph is produced by exporting the results from the user-defined function PRESSURE and is plotted against node number to show the pressure distribution at the end time:

Figure 249: Pressure Along the Shock Tube at t = 0.8 ms

Pressure Along the Shock Tube at t = 0.8 ms

Figure 250: Density of Gas Along the Shock Tube at t = 0.8 ms

Density of Gas Along the Shock Tube at t = 0.8 ms

Figure 251: Pressure Distribution in Shock Tube at t = 0.8 ms

Pressure Distribution in Shock Tube at t = 0.8 ms