|
Input files
|
ShockTube.KIN
Download ShockTube.zip
|
|
Description
|
In the one-dimensional calculation, the shock tube has a high-pressure region (-10.0 m ≤ z ≤ 0.0 m) and a low pressure region (0.0 m ≤ z ≤ 10.0 m). A diaphragm at z = 0.0 m ruptures at t = 0.0 s. A shock wave propagates in the +z direction and a rarefaction wave in the -z-direction. The material is an ideal gamma-law gas with gamma = 1.4. The high pressure region has initial conditions are p = 10^5 Pa and rho = 1.0 kG/m3. The initial conditions in the low pressure region are p = 10^4 Pa and rho = 0.125 kG/m3. For a gamma-law gas, the internal energy is given by u = p/(Gamma-1)*rho.
The internal energy in the high-pressure region is u = 2.5E5 J/kg and is equal to u = 2.0E5 J/kg in the low-pressure region.
|
|
Results
|
The figure shows a plot of material density as a function of z at t = 10.05 ms along with a inset of the Comsol result. The symbols in the KB1V plot mark the centers of an initially uniform set of elements. In the Lagrangian description, the spacing between symbols indicates compression and expansion. The symbol color shows the material identity of the elements. The term "contact wave" designates the expansion of the high pressure material into the low pressure region.
|
|
Comments
|
For a comparison, the example uses the same geometry and parameters as a Comsol benchmark test : https://www.comsol.com/model/shock-tube-43591. The results from the Lagrangian KB calculation closely match those of the Comsol Eulerian model.
|
