Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

7.6 Moving Expansion Waves

Chalmers University of Technology
Department of Mechanical Engineering
Division of Fluid Dynamics

The expansion wave propagation into the driver section in a shock tube can be described using characteristic lines.

expansion fan

Figure 7.3:Expansion fan centered at (x,t)=(0.0,0.0)(x,t)=(0.0,0.0)

The expansion is propagating into stagnant fluid in region four (the driver section), which means that the flow properties ahead of the expansion wave are constant.

Ja+=Jb+J^+_a=J^+_b

J+J^+ invariants constant along C+C^+ characteristics

Ja+=Jc+=Je+J^+_a=J^+_c=J^+_e
Jb+=Jd+=Jf+J^+_b=J^+_d=J^+_f

Since Ja+=Jb+J^+_a=J^+_b this also implies Je+=Jf+J^+_e=J^+_f. In fact, since the flow properties ahead of the expansion are constant, all C+C^+ lines will have the same J+J^+ value.

JJ^- invariants constant along CC^- characteristics

Jc=JdJ^-_c=J^-_d
Je=JfJ^-_e=J^-_f
ue=12(Je++Je)uf=12(Jf++Jf)Je=JfJe+=Jf+}ue=ufae=af\left. \begin{aligned} &u_e=\frac{1}{2}(J^+_e+J^-_e)\\ &u_f=\frac{1}{2}(J^+_f+J^-_f)\\ &J^-_e=J^-_f\\ &J^+_e=J^+_f \end{aligned} \right\}\Rightarrow u_e=u_f \Rightarrow a_e=a_f

Due to the fact the J+J^+ is constant in the entire expansion region, uu and aa will be constant along each CC^- line.

The constant J+J^+ value can be used to obtain relations for the variation of flow properties through the expansion region. Evaluation of the J+J^+ invariant at any position within the expansion region should give the same value as in region 4.

u+2aγ1=u4+2a4γ1=0+2a4γ1u+\frac{2a}{\gamma-1}=u_4+\frac{2a_4}{\gamma-1}=0+\frac{2a_4}{\gamma-1}

and thus

aa4=1γ12(ua4)\frac{a}{a_4}=1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)

Eqn. (7.145) and a=γRTa=\sqrt{\gamma RT} gives

TT4=[1γ12(ua4)]2\frac{T}{T_4}=\left[1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)\right]^2

Using isentropic relations, we can get pressure ratio and density ratio

pp4=[1γ12(ua4)]2γ/(γ1)\frac{p}{p_4}=\left[1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)\right]^{2\gamma/(\gamma-1)}
ρρ4=[1γ12(ua4)]2/(γ1)\frac{\rho}{\rho_4}=\left[1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)\right]^{2/(\gamma-1)}