Moving expansion waves: Difference between revisions

Latest revision as of 13:36, 1 April 2026

Moving Expansion Waves

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

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.

J_{a}^{+} = J_{b}^{+}

(Eq. 6.135)

$J^{+}$ invariants constant along $C^{+}$ characteristics

J_{a}^{+} = J_{c}^{+} = J_{e}^{+}

(Eq. 6.136)

J_{b}^{+} = J_{d}^{+} = J_{f}^{+}

(Eq. 6.137)

Since $J_{a}^{+} = J_{b}^{+}$ this also implies $J_{e}^{+} = J_{f}^{+}$ . In fact, since the flow properties ahead of the expansion are constant, all $C^{+}$ lines will have the same $J^{+}$ value.

$J^{-}$ invariants constant along $C^{-}$ characteristics

J_{c}^{-} = J_{d}^{-}

(Eq. 6.138)

J_{e}^{-} = J_{f}^{-}

(Eq. 6.139)

\begin{matrix} 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{matrix}} \Rightarrow u_{e} = u_{f} \Rightarrow a_{e} = a_{f}

(Eq. 6.140)

Due to the fact the $J^{+}$ is constant in the entire expansion region, $u$ and $a$ will be constant along each $C^{-}$ line.

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

u + \frac{2 a}{γ - 1} = u_{4} + \frac{2 a_{4}}{γ - 1} = 0 + \frac{2 a_{4}}{γ - 1}

(Eq. 6.141)

and thus

\frac{a}{a_{4}} = 1 - \frac{γ - 1}{2} (\frac{u}{a_{4}})

(Eq. 6.142)

Eqn. \ref{eq:expansion:a} and $a = \sqrt{γ R T}$ gives

\frac{T}{T_{4}} = {[1 - \frac{γ - 1}{2} (\frac{u}{a_{4}})]}^{2}

(Eq. 6.143)

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

\frac{p}{p_{4}} = {[1 - \frac{γ - 1}{2} (\frac{u}{a_{4}})]}^{2 γ / (γ - 1)}

(Eq. 6.144)

\frac{ρ}{ρ_{4}} = {[1 - \frac{γ - 1}{2} (\frac{u}{a_{4}})]}^{2 / (γ - 1)}

(Eq. 6.145)

Moving expansion waves: Difference between revisions

Latest revision as of 13:36, 1 April 2026

Contents

Moving Expansion Waves

Navigation menu

@@ Line 1: / Line 1: @@
-[[Category:Compressible flow]]
+<!--
-[[Category:Unsteady waves]]
+-->[[Category:Compressible flow]]<!--
-[[Category:Inviscid flow]]
+-->[[Category:Unsteady waves]]<!--
+-->[[Category:Inviscid flow]]<!--
+--><noinclude><!--
+-->[[Category:Compressible flow:Topic]]<!--
+--></noinclude><!--
-__TOC__
+--><nomobile><!--
+-->__TOC__<!--
+--></nomobile><!--
-\section{Moving Normal Shock Waves}
+--><noinclude><!--
+-->{{#vardefine:secno|6}}<!--
+-->{{#vardefine:eqno|134}}<!--
+--></noinclude><!--
+-->
+==== Moving Expansion Waves ====
-\noindent The starting point is the governing equations for stationary normal shocks (repeated here for convenience).
+The expansion wave propagation into the driver section in a shock tube can be described using characteristic lines.
-\begin{equation}
+<!--
-\rho_1 u_1 = \rho_2 u_2
+\begin{figure}[ht!]
-\label{eq:stationary:cont}
+\begin{center}
-\end{equation}
+\includegraphics[]{figures/standalone-figures/Chapter07/pdf/expansion-region.pdf}
+\caption{Expansion fan centered at $(x,t)=(0.0,0.0)$}
+\label{fig:characteristics}
+\end{center}
+\end{figure}
+-->
-\begin{equation}
+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.
-\rho_1 u_1^2+p_1 = \rho_2 u_2^2 + p_2
-\label{eq:stationary:mom}
-\end{equation}
-\begin{equation}
+{{NumEqn|<math>
-h_1 + \frac{1}{2}u_1^2 = h_2 + \frac{1}{2}u_2^2
+J^+_a=J^+_b
-\label{eq:stationary:energy}
+</math>}}
-\end{equation}
-\noindent Shock moving to the right with the constant speed $W$ into a gas that is standing still. Moving with the shock, we would see a gas velocity ahead of the shock $u_1=W$, and the gas behind the shock moves to the right with the velocity $u_2=W-u_p$. Now, let's insert $u_1$ and $u_2$ in the stationary shock relations \ref{eq:stationary:cont} - \ref{eq:stationary:energy}.
+<math>J^+</math> invariants constant along <math>C^+</math> characteristics
-\begin{equation}
+{{NumEqn|<math>
-\rho_1 W = \rho_2 (W-u_p)
+J^+_a=J^+_c=J^+_e
-\label{eq:unsteady:cont}
+</math>}}
-\end{equation}
-\begin{equation}
+{{NumEqn|<math>
-\rho_1 W^2+p_1 = \rho_2 (W-u_p)^2 + p_2
+J^+_b=J^+_d=J^+_f
-\label{eq:unsteady:mom}
+</math>}}
-\end{equation}
-\begin{equation}
+Since <math>J^+_a=J^+_b</math> this also implies <math>J^+_e=J^+_f</math>. In fact, since the flow properties ahead of the expansion are constant, all <math>C^+</math> lines will have the same <math>J^+</math> value.
-h_1 + \frac{1}{2}W^2 = h_2 + \frac{1}{2}(W-u_p)^2
-\label{eq:unsteady:energy}
-\end{equation}
-%\newpage
+<math>J^-</math> invariants constant along <math>C^-</math> characteristics
-\noindent Rewriting Eqn. \ref{eq:unsteady:cont}
+{{NumEqn|<math>
+J^-_c=J^-_d
+</math>}}
-\begin{equation}
+{{NumEqn|<math>
-(W-u_p) = W \frac{\rho_1}{\rho_2}
+J^-_e=J^-_f
-\label{eq:unsteady:cont:mod}
+</math>}}
-\end{equation}\\
-\noindent Inserting Eqn. \ref{eq:unsteady:cont:mod} in Eqn. \ref{eq:unsteady:mom} gives\\
+{{NumEqn|<math>
+\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
+</math>}}
-\begin{equation*}
+Due to the fact the <math>J^+</math> is constant in the entire expansion region, <math>u</math> and <math>a</math> will be constant along each <math>C^-</math> line.
-p_1+\rho_1 W^2 = p_2+\rho_2 W^2\left(\frac{\rho_1}{\rho_2}\right)^2 \Rightarrow p_2-p_1 = \rho_1W^2\left(1-\frac{\rho_1}{\rho_2}\right)
-%\label{eq:unsteady:mom:mod}
-\end{equation*}\\
-\begin{equation}
+The constant <math>J^+</math> value can be used to obtain relations for the variation of flow properties through the expansion region. Evaluation of the <math>J^+</math> invariant at any position within the expansion region should give the same value as in region 4.
-W^2=\frac{p_2-p_1}{\rho_2-\rho_1}\left(\frac{\rho_2}{\rho_1}\right)
-\label{eq:unsteady:mom:mod}
-\end{equation}\\
-\noindent From the continuity equation \ref{eq:unsteady:cont}, we get \\
+{{NumEqn|<math>
+u+\frac{2a}{\gamma-1}=u_4+\frac{2a_4}{\gamma-1}=0+\frac{2a_4}{\gamma-1}
+</math>}}
-\begin{equation}
+and thus
-W = (W-u_p) \left(\frac{\rho_2}{\rho_1}\right)
-\label{eq:unsteady:cont:modb}
-\end{equation}\\
-\noindent Inserting Eqn. \ref{eq:unsteady:cont:modb} in Eqn. \ref{eq:unsteady:mom:mod} gives\\
+{{NumEqn|<math>
+\frac{a}{a_4}=1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)
+</math>}}
-\begin{equation}
+Eqn. \ref{eq:expansion:a} and <math>a=\sqrt{\gamma RT}</math> gives
-(W-u_p)^2=\frac{p_2-p_1}{\rho_2-\rho_1}\left(\frac{\rho_1}{\rho_2}\right)
-\label{eq:unsteady:mom:modb}
-\end{equation}\\
-\noindent Now, let's insert Eqns. \ref{eq:unsteady:mom:mod} and \ref{eq:unsteady:mom:modb} in the energy equation (Eqn. \ref{eq:unsteady:energy}).\\
+{{NumEqn|<math>
+\frac{T}{T_4}=\left[1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)\right]^2
+</math>}}
-\begin{equation}
+Using isentropic relations, we can get pressure ratio and density ratio
-h_1 + \frac{1}{2}\left[\frac{p_2-p_1}{\rho_2-\rho_1}\left(\frac{\rho_2}{\rho_1}\right)\right] = h_2 + \frac{1}{2}\left[\frac{p_2-p_1}{\rho_2-\rho_1}\left(\frac{\rho_1}{\rho_2}\right)\right]
-\label{eq:unsteady:energy:mod}
-\end{equation}
-\begin{equation}
+{{NumEqn|<math>
-h=e+\frac{p}{\rho}
+\frac{p}{p_4}=\left[1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)\right]^{2\gamma/(\gamma-1)}
-\label{eq:enthalpy}
+</math>}}
-\end{equation}
-\begin{equation}
+{{NumEqn|<math>
-e_1 + \frac{p_1}{\rho_1} +  \frac{1}{2}\left[\frac{p_2-p_1}{\rho_2-\rho_1}\left(\frac{\rho_2}{\rho_1}\right)\right] = e_2 + \frac{p_2}{\rho_2} + \frac{1}{2}\left[\frac{p_2-p_1}{\rho_2-\rho_1}\left(\frac{\rho_1}{\rho_2}\right)\right]
+\frac{\rho}{\rho_4}=\left[1-\frac{\gamma-1}{2}\left(\frac{u}{a_4}\right)\right]^{2/(\gamma-1)}
-\label{eq:unsteady:energy:mod:b}
+</math>}}
-\end{equation}\\
-\noindent which can be rewritten as\\
-\begin{equation}
-e_2-e_1=\frac{p_1+p_2}{2}\left(\frac{1}{\rho_1}-\frac{1}{\rho_2}\right)
-\label{eq:unsteady:hugonoit}
-\end{equation}\\
-\noindent Eqn \ref{eq:unsteady:hugonoit} is the same Hugoniot equation as we get for a stationary normal shock. The Hugoniot equation is a relation of thermodynamic properties over a shock. As the shock in the unsteady case is moving with a constant velocity, the frame of reference moving with the shock is an inertial frame and thus the same physical relations apply in the moving shock case as in the stationary shock case. The fact that the Hugoniot relation does not include any velocities or Mach numbers but only thermodynamic properties, the relation will be unchanged for a moving shock.
-\subsection{Moving Shock Relations}
-\noindent For a calorically perfect gas we have $e=C_v T$. Inserted in the Hugoniot relation above this gives\\
-\begin{equation}
-C_v(T_2-T_1)=\frac{p_1+p_2}{2}\left(\nu_1-\nu_2\right)
-\label{eq:unsteady:hugonoit:b}
-\end{equation}\\
-where $\nu=1/\rho$\\
-\noindent Now, using the ideal gas law $T=p\nu/R$ and $C_v/R=1/(\gamma-1)$ gives\\
-\begin{equation*}
-\left(\frac{1}{\gamma-1}\right)(p_2\nu_2-p_1\nu_1)=\frac{p_1+p_2}{2}\left(\nu_1-\nu_2\right)
-\end{equation*}
-\begin{equation*}
-\Leftrightarrow
-\end{equation*}
-\begin{equation*}
-p_2\left(\frac{\nu_2}{\gamma-1}-\frac{\nu_1-\nu_2}{2}\right)=p_1\left(\frac{\nu_1}{\gamma-1}+\frac{\nu_1-\nu_2}{2}\right)
-\end{equation*}\\
-\noindent From this result, we can derive a relation for the pressure ratio over the shock as a function of density ratio\\
-\begin{equation}
-\frac{p_2}{p_1}=\frac{\left(\dfrac{\gamma+1}{\gamma-1}\right)\left(\dfrac{\nu_1}{\nu_2}\right)-1}{\left(\dfrac{\gamma+1}{\gamma-1}\right)-\left(\dfrac{\nu_1}{\nu_2}\right)}
-\label{eq:unsteady:hugonoit:c}
-\end{equation}\\
-\noindent $\nu=RT/p$ and thus
-\begin{equation}
-\frac{\nu_1}{\nu_2}=\frac{T_1}{T_2}\frac{p_2}{p_1}
-\label{eq:unsteady:density:ratio}
-\end{equation}\\
-\noindent Eqn. \ref{eq:unsteady:density:ratio} in Eqn. \ref{eq:unsteady:hugonoit:c} gives\\
-\begin{equation}
-\frac{p_2}{p_1}=\frac{\left(\dfrac{\gamma+1}{\gamma-1}\right)\left(\dfrac{T_1}{T_2}\dfrac{p_2}{p_1}\right)-1}{\left(\dfrac{\gamma+1}{\gamma-1}\right)-\left(\dfrac{T_1}{T_2}\dfrac{p_2}{p_1}\right)}
-\label{eq:unsteady:hugonoit:c}
-\end{equation}\\
-\noindent Now, we can get a relation for calculation of the temperature ratio over the moving shock as function of the shock pressure ratio\\
-\begin{equation}
-\frac{T_2}{T_1}=\frac{p_2}{p_1}\left[\frac{\left(\dfrac{\gamma+1}{\gamma-1}\right)+\left(\dfrac{p_2}{p_1}\right)}{1+\left(\dfrac{\gamma+1}{\gamma-1}\right)\left(\dfrac{p_2}{p_1}\right)}\right]
-\label{eq:unsteady:temperature:ratio}
-\end{equation}\\
-\noindent Once again using the ideal gas law\\
-\begin{equation}
-\frac{\rho_2}{\rho_1}=\frac{\left(\dfrac{\gamma+1}{\gamma-1}\right)+\left(\dfrac{p_2}{p_1}\right)}{1+\left(\dfrac{\gamma+1}{\gamma-1}\right)\left(\dfrac{p_2}{p_1}\right)}
-\label{eq:unsteady:density:ratio}
-\end{equation}\\
-\noindent Going back to the momentum equation\\
-\begin{equation*}
-p_2-p_1 = \rho_1W^2\left(1-\frac{\rho_1}{\rho_2}\right)=\left\{W=M_s a_1\right\}=\rho_1M_s^2a_1^2\left(1-\frac{\rho_1}{\rho_2}\right)
-\end{equation*}\\
-\noindent with $a_1^2=\gamma p_1/\rho_1$, we get\\
-\begin{equation}
-\frac{p_2}{p_1} = \gamma M_s^2\left(1-\frac{\rho_1}{\rho_2}\right)+1
-\label{eq:unsteady:Mach:a}
-\end{equation}\\
-\noindent From the normal shock relations, we have\\
-\begin{equation}
-\frac{\rho_1}{\rho_2} = \frac{2+(\gamma-1)M_s^2}{(\gamma+1)M_s^2}
-\label{eq:unsteady:Mach:b}
-\end{equation}\\
-\noindent Eqn. \ref{eq:unsteady:Mach:b} in \ref{eq:unsteady:Mach:a} gives\\
-\begin{equation}
-\frac{p_2}{p_1} = 1 + \left(\frac{2\gamma}{\gamma+1}\right)(M_s^2-1)
-\label{eq:unsteady:Mach:c}
-\end{equation}\\
-or\\
-\begin{equation}
-M_s=\sqrt{\left(\frac{\gamma+1}{2\gamma}\right)\left(\frac{p_2}{p_1}-1\right)+1}
-\label{eq:unsteady:Mach}
-\end{equation}\\
-\noindent Eqn. \ref{eq:unsteady:Mach} with $M_s=W/a_1$\\
-\begin{equation}
-W=a_1\sqrt{\left(\frac{\gamma+1}{2\gamma}\right)\left(\frac{p_2}{p_1}-1\right)+1}
-\label{eq:unsteady:W}
-\end{equation}\\
-\subsection{Induced Flow Behind Moving Shock}
-\noindent Let's try to find a relation for calculation of the induced velocity behind the moving shock. Once again, the starting point is the continuity equation for moving shocks (Eqn. \ref{eq:unsteady:cont}) repeated here for convenience\\
-\begin{equation*}
-\rho_1 W = \rho_2 (W-u_p)
-\end{equation*}\\
-\noindent The induced velocity appears on the right side of the continuity equation\\
-\begin{equation*}
-W (\rho_1-\rho_2) = -\rho_2 u_p
-\end{equation*}\\
-\begin{equation}
-u_p = W \left(1-\frac{\rho_1}{\rho_2}\right)
-\label{eq:unsteady:up:a}
-\end{equation}\\
-\noindent From before we have a relation for $W$ as a function of pressure ratio and one for $\rho_1/\rho_2$, also as a function of pressure ratio.\\
-Eqn. \ref{eq:unsteady:up:a} togheter with Eqns. \ref{eq:unsteady:W} and \ref{eq:unsteady:density:ratio} gives\\
-\begin{equation}
-u_p=a_1\underbrace{\sqrt{\left(\frac{\gamma+1}{2\gamma}\right)\left(\frac{p_2}{p_1}-1\right)+1}}_{I}\underbrace{\left[1-\dfrac{\left(\dfrac{\gamma+1}{\gamma-1}\right)+\left(\dfrac{p_2}{p_1}\right)}{1+\left(\dfrac{\gamma+1}{\gamma-1}\right)\left(\dfrac{p_2}{p_1}\right)}\right]}_{II}
-\label{eq:unsteady:up:b}
-\end{equation}\\
-\noindent The equation subsets I and II can be rewritten as:\\
-Term I:
-\begin{equation*}
-\sqrt{\left(\frac{\gamma+1}{2\gamma}\right)\left(\frac{p_2}{p_1}-1\right)+1}=\sqrt{\frac{\gamma+1}{2\gamma}\left[\left(\frac{p_2}{p_1}\right)+\left(\frac{\gamma-1}{\gamma+1}\right)\right]}
-\end{equation*}\\
-Term II:
-\begin{equation*}
-\left[1-\dfrac{\left(\dfrac{\gamma+1}{\gamma-1}\right)+\left(\dfrac{p_2}{p_1}\right)}{1+\left(\dfrac{\gamma+1}{\gamma-1}\right)\left(\dfrac{p_2}{p_1}\right)}\right]=\frac{1}{\gamma}\left(\frac{p_2}{p_1}-1\right)\frac{\left(\dfrac{2\gamma}{\gamma+1}\right)}{\left(\dfrac{\gamma-1}{\gamma+1}\right)+\left(\dfrac{p_2}{p_1}\right)}
-\end{equation*}\\
-\noindent With the rewritten terms I and II implemented, Eqn. \ref{eq:unsteady:up:b} becomes\\
-\begin{equation}
-u_p=\frac{a_1}{\gamma}\left(\frac{p_2}{p_1}-1\right)\sqrt{\frac{\left(\dfrac{2\gamma}{\gamma+1}\right)}{\left(\dfrac{\gamma-1}{\gamma+1}\right)+\left(\dfrac{p_2}{p_1}\right)}}
-\label{eq:unsteady:up}
-\end{equation}\\
-\noindent Since the region behind the moving shock is region 2, the induced flow Mach number is obtained as\\
-\begin{equation*}
-M_p=\frac{u_p}{a_2}=\frac{u_p}{a_1}\frac{a_1}{a_2}=\frac{u_p}{a_1}\sqrt{\frac{\gamma R T_1}{\gamma R T_2}}=\frac{u_p}{a_1}\sqrt{\frac{T_1}{T_2}}
-\end{equation*}\\
-\noindent With $up/a_1$ from Eqn. \ref{eq:unsteady:up}  and $T_1/T_2$ from Eqn. \ref{eq:unsteady:temperature:ratio}\\
-\begin{equation}
-M_p=\frac{1}{\gamma}\left(\frac{p_2}{p_1}-1\right)\left(\frac{\left(\dfrac{2\gamma}{\gamma+1}\right)}{\left(\dfrac{\gamma-1}{\gamma+1}\right)+\left(\dfrac{p_2}{p_1}\right)}\right)^{1/2}
-\left(\frac{1+\left(\dfrac{\gamma+1}{\gamma-1}\right)\left(\dfrac{p_2}{p_1}\right)}{\left(\dfrac{\gamma+1}{\gamma-1}\right)\left(\dfrac{p_2}{p_1}\right)+\left(\dfrac{p_2}{p_1}\right)^2}\right)^{1/2}
-\label{eq:unsteady:Mp}
-\end{equation}\\
-\noindent There is a theoretical upper limit for the induced Mach number $M_p$\\
-\begin{equation*}
-\lim_{p_2/p_1\rightarrow\infty} M_p\left(\frac{p_2}{p_1}\right)=\sqrt{\frac{2}{\gamma(\gamma-1)}}
-\end{equation*}\\
-\noindent As can be seen, at the upper limit the induced Mach number is a function of $\gamma$ and for air ($\gamma=1.4$) we get\\
-\begin{equation*}
-\lim_{p_2/p_1\rightarrow\infty} M_p\left(\frac{p_2}{p_1}\right)\simeq 1.89
-\end{equation*}\\

Moving expansion waves: Difference between revisions

Latest revision as of 13:36, 1 April 2026

Moving Expansion Waves

Navigation menu

Search