One-dimensional inviscid flow: Difference between revisions

Latest revision as of 09:02, 14 March 2026

Acoustic Waves

In Fig. \ref{fig:soundwave}, station 1 represents the flow state ahead of the sound wave and station 2 the flow state behind the sound wave. Set up the continuity equation for one-dimensional flows between 1 and 2. If we could change frame of reference and follow the sound wave, we would see fluid approaching the wave with the propagation speed of the wave, $a$ , and behind the wave, the fluid would have a slightly modified speed, $a + d a$ . There would also be a slight in all other flow properties. Let's apply the one-dimensional continuity equation between station 1 and station 2.

$ρ_{1} u_{1} = ρ_{2} u_{2}$

$ρ a = (ρ + d ρ) (a + d a)$

$ρ a = ρ a + ρ d a + a d ρ + \underset{\sim 0}{\underset{⏟}{d ρ d a}} \Rightarrow$

$a = - ρ \frac{d a}{d ρ}$

The one-dimensional momentum equation between station 1 and station 2 gives

$ρ_{1} u_{1}^{2} + p_{1} = ρ_{2} u_{2}^{2} + p_{2}$

$ρ a^{2} + p = (ρ + d ρ) (a + d a)^{2} + (p + d p)$

$ρ a^{2} + p = ρ a^{2} + 2 ρ a d a + \underset{\sim 0}{\underset{⏟}{ρ d a^{2}}} + d ρ a^{2} + \underset{\sim 0}{\underset{⏟}{2 d ρ a d a}} + \underset{\sim 0}{\underset{⏟}{d ρ d a^{2}}} + p + d p \Rightarrow$

$d p = - 2 ρ a d a - d ρ a^{2} \Rightarrow$

$d a = - \frac{d p + d ρ a^{2}}{2 ρ a} = - \frac{d ρ}{2 a ρ} (\frac{d p}{d ρ} + a^{2}) \Rightarrow$

$\frac{d a}{d ρ} = - \frac{1}{2 a ρ} (\frac{d p}{d ρ} + a^{2})$

Eqn. \ref{eq:speedofsound:b} in \ref{eq:speedofsound:a} gives

$a = \frac{1}{2 a} (\frac{d p}{d ρ} + a^{2}) \Rightarrow$

$a^{2} = \frac{d p}{d ρ}$

Sound wave:

gradients are small
irreversible (dissipative effects are negligible)
no heat addition

Thus, the change of flow properties as the sound wave passes can be assumed to be an isentropic process

$a^{2} = {(\frac{d p}{d ρ})}_{s}$

$a = \sqrt{{(\frac{d p}{d ρ})}_{s}} = \sqrt{\frac{1}{ρ τ_{s}}}$

where $τ_{s}$ is the compressibility of the gas. Eqn. \ref{eq:speedofsound:d} is valid for all gases. It can be seen from the equation, that truly incompressible flow ( $τ_{s} = 0$ ) would imply infinite speed of sound.

Since the process is isentropic, we can use the isentropic relations if we also assume the gas to be calorically perfect\\

$\frac{p_{2}}{p_{1}} = {(\frac{ρ_{2}}{ρ_{1}})}^{γ} \Rightarrow p = C ρ^{γ}$

$a^{2} = {(\frac{d p}{d ρ})}_{s} = γ C ρ^{γ - 1} = γ \underset{= p}{\underset{⏟}{[C ρ^{γ}]}} ρ^{- 1} = \frac{γ p}{ρ} \Rightarrow$

$a = \sqrt{\frac{γ p}{ρ}}$

or

$a = \sqrt{γ R T}$

From the relation above, it is obvious that the local speed of sound is related to the temperature of the flow, which in turn is a measure of the motion of elementary particles (atoms and/or molecules) of the fluid at a specific location. This stems from the fact that sound waves are propagated via interaction of these elementary particles. Since information in a flow is propagated via molecular interaction the relation between the speed at which this information is conveyed and the speed of the flow has important physical implications. Figure~\ref{fig:speed:of:sound} compares three sound wave patterns generated by a a beacon. In the left picture, the sound transmitter is stationary and thus the acoustic waves are centered around the transmitter. In the middle image, the transmitter is moving to the left at a speed less than the speed of sound and thus the transmitter will always be within all sound wave circles but it will be off-centered with a bias in the direction that the transmitter is moving. In the right image the transmitter is moving faster than the speed of sound and thus it will always be located outside of all acoustic waves. In a supersonic flow, no information can travel upstream and therefore there is no way for the flow to adjust to downstream obstacles. This is compensated for by the introduction of shocks in the flow. Over a shock flow properties changes discontinuity. An example is given in Figure~\ref{fig:supersonic:flow}.

Shock Waves

The starting point is to set up the governing equations for one-dimensional steady compressible flow over a control volume enclosing the normal shock (Fig. \ref{fig:shock:cv}).

continuity:

$ρ_{1} u_{1} = ρ_{2} u_{2}$

momentum:

$ρ_{1} u_{1}^{2} + p_{1} = ρ_{2} u_{2}^{2} + p_{2}$

energy:

$h_{1} + \frac{1}{2} u_{1}^{2} = h_{2} + \frac{1}{2} u_{2}^{2}$

Divide the momentum equation by $ρ_{1} u_{1}$

$\frac{1}{ρ_{1} u_{1}} (ρ_{1} u_{1}^{2} + p_{1}) = \frac{1}{ρ_{1} u_{1}} (ρ_{2} u_{2}^{2} + p_{2}) = {ρ_{1} u_{1} = ρ_{2} u_{2}} = \frac{1}{ρ_{2} u_{2}} (ρ_{2} u_{2}^{2} + p_{2}) \Rightarrow$

$\frac{p_{1}}{ρ_{1} u_{1}} - \frac{p_{2}}{ρ_{2} u_{2}} = u_{2} - u_{1}$

For a calorically perfect gas $a = \sqrt{γ p / ρ}$ , which if implemented in Eqn. \ref{eq:governing:mom:b} gives

$\frac{a_{1}^{2}}{γ u_{1}} - \frac{a_{2}^{2}}{γ u_{2}} = u_{2} - u_{1}$

The energy equation (Eqn. \ref{eq:governing:energy}) with $h = C_{p} T$

$C_{p} T_{1} + \frac{1}{2} u_{1}^{2} = C_{p} T_{2} + \frac{1}{2} u_{2}^{2}$

Replacing $C_{p}$ with $γ R / (γ - 1)$ gives

$\frac{γ R T_{1}}{γ - 1} + \frac{1}{2} u_{1}^{2} = \frac{γ R T_{2}}{γ - 1} + \frac{1}{2} u_{2}^{2}$

With $a = \sqrt{γ R T}$ this becomes

$\frac{a_{1}^{2}}{γ - 1} + \frac{1}{2} u_{1}^{2} = \frac{a_{2}^{2}}{γ - 1} + \frac{1}{2} u_{2}^{2}$

Eqn. \ref{eq:governing:energy:d} can be set up between any two points in the flow. Specifically, we can use the relation to relate the flow velocity, $u$ , and speed of sound, $a$ , in any point to the corresponding flow properties at sonic conditions ( $u = a = a^{*}$ ).

$\frac{a^{2}}{γ - 1} + \frac{1}{2} u^{2} = \frac{γ + 1}{2 (γ - 1)} {a^{*}}^{2}$

If Eqn. \ref{eq:governing:energy:e} is evaluated in locations 1 and 2, we get

$\begin{matrix} a_{1}^{2} = \frac{γ + 1}{2} {a^{*}}^{2} - \frac{γ - 1}{2} u_{1}^{2} \\ a_{2}^{2} = \frac{γ + 1}{2} {a^{*}}^{2} - \frac{γ - 1}{2} u_{2}^{2} \end{matrix}$

Since the change in flow conditions over the shock is adiabatic (no heat is added inside the shock), critical properties will be constant over the shock. Especially $a^{*}$ will be constant.

Eqn. \ref{eq:governing:energy:f} inserted in \ref{eq:governing:mom:c} gives\\

$\frac{1}{γ u_{1}} (\frac{γ + 1}{2} {a^{*}}^{2} - \frac{γ - 1}{2} u_{1}^{2}) - \frac{1}{γ u_{2}} (\frac{γ + 1}{2} {a^{*}}^{2} - \frac{γ - 1}{2} u_{2}^{2}) = u_{2} - u_{1} \Rightarrow$

$(\frac{γ + 1}{2 γ}) {a^{*}}^{2} (\frac{1}{u_{1}} - \frac{1}{u_{2}}) = (\frac{γ + 1}{2 γ}) (u_{2} - u_{1}) \Rightarrow$

${a^{*}}^{2} (\frac{1}{u_{1}} - \frac{1}{u_{2}}) = (u_{2} - u_{1}) \Rightarrow$

${a^{*}}^{2} (\frac{u_{2}}{u_{1} u_{2}} - \frac{u_{1}}{u_{1} u_{2}}) = (u_{2} - u_{1}) \Rightarrow$

$\frac{1}{u_{1} u_{2}} {a^{*}}^{2} (u_{2} - u_{1}) = (u_{2} - u_{1}) \Rightarrow$

${a^{*}}^{2} = u_{1} u_{2}$

Eqn. \ref{eq:prandtl} is sometimes referred to as the Prandtl relation. Divide the Prandtl relation by ${a^{*}}^{2}$ on both sides gives

$1 = \frac{u_{1}}{a^{*}} \frac{u_{2}}{a^{*}} = M_{1}^{*} M_{2}^{*}$

or

$M_{2}^{*} = \frac{1}{M_{1}^{*}}$

The relation between $M^{*}$ and $M$ is given by

${M^{*}}^{2} = \frac{(γ + 1) M^{2}}{2 + (γ - 1) M^{2}}$

from which is can be seen that $M^{*}$ will follow the Mach number $M$ in the sense that

$M = 1 \Rightarrow M^{*} = 1$
$M < 1 \Rightarrow M^{*} < 1$
$M > 1 \Rightarrow M^{*} > 1$

Eqn. \ref{eq:MachStar} inserted in Eqn. \ref{eq:NormalMach} gives

$\frac{(γ + 1) M_{1}^{2}}{2 + (γ - 1) M_{1}^{2}} = \frac{2 + (γ - 1) M_{2}^{2}}{(γ + 1) M_{2}^{2}}$

$M_{2}^{2} = \frac{1 + [(γ - 1) / 2] M_{1}^{2}}{γ M_{1}^{2} - (γ - 1) / 2}$

The Mach number relations above effectively show that if the Mach number upstream of the shock is greater than one, the downstream Mach number must be less than one and vice versa. We can also see that a sonic upstream flow $M_{1} = 1.0$ gives sonic flow downstream of the shock. So, apparently the relation as such holds for both supersonic and subsonic upstream flow mathematically. The question is if it is also physically correct. For a supersonic upstream flow we will get a discontinuous compression and if the flow upstream of the control volume is subsonic we will instead get a discontinuous expansion inside the control volume but, again, is this physically correct? We will get the answer by analyzing the entropy change over the control volume.

Analyzing the energy equation and the second law of thermodynamics shows that there is a direct relation between entropy increase and total pressure drop.

$s_{2} - s_{1} = C_{p} \ln \frac{T_{2}}{T_{1}} - R \ln \frac{p_{2}}{p_{1}}$

$s_{2} - s_{1} = C_{p} \ln \frac{T_{2}}{T_{o_{2}}} \frac{T_{o_{1}}}{T_{1}} \frac{T_{o_{2}}}{T_{o_{1}}} - R \ln \frac{p_{2}}{p_{o_{2}}} \frac{p_{o_{1}}}{p_{1}} \frac{p_{o_{2}}}{p_{o_{1}}}$

using the isentropic relations we get

$s_{2} - s_{1} = C_{p} \ln \frac{T_{o_{2}}}{T_{o_{1}}} - R \ln \frac{p_{o_{2}}}{p_{o_{1}}}$

and since the process is adiabatic and thus $T_{o_{2}} = T_{o_{1}}$ the change in entropy is directly related to the change in total pressure as

$s_{2} - s_{1} = - R \ln \frac{p_{o_{2}}}{p_{o_{1}}}$

or

$\frac{p_{o_{2}}}{p_{o_{1}}} = e^{- (s_{2} - s_{1}) / R}$

Figure~\ref{fig:shock:entropy} shows the entropy change over a normal shock. As can be seen in the figure, a subsonic upstream Mach number leads to a reduction of entropy, which once and for all rules out all such solutions as non-physical and thus the question about the upstream conditions can now be considered answered. This in turn implies that the Mach number downstream of a normal shock will always be subsonic, which can be seen in Fig~\ref{fig:shock:downstream:Mach} below.

By rewriting the right-hand side of Eqn.\ref{eq:NormalMach:b}, it is easy to realize that the downstream Mach number $M_{2}$ approaches a finite value for large values of the upstream Mach number, $M_{1}$ .

${M_{2}^{2} |}_{M_{1} \to \infty} = {\frac{2 / M_{1}^{2} + (γ - 1)}{2 γ - (γ - 1) / M_{1}^{2}} |}_{M_{1} \to \infty} = \frac{γ - 1}{2 γ}$

Normal Shock Relations

Rewriting the continuity equation (Eqn. \ref{eq:governing:cont})

$\frac{ρ_{2}}{ρ_{1}} = \frac{u_{1}}{u_{2}} = \frac{u_{1}^{2}}{u_{1} u_{2}} = {{a^{*}}^{2} = u_{1} u_{2}} = \frac{u_{1}^{2}}{{a^{*}}^{2}} = {M_{1}^{*}}^{2}$

Eqn. \ref{eq:MachStar} in Eqn. \ref{eq:Normal:density:a} gives

$\frac{ρ_{2}}{ρ_{1}} = \frac{(γ + 1) M_{1}^{2}}{2 + (γ - 1) M_{1}^{2}}$

To get a corresponding relation for the pressure ratio over the shock, we go back to the momentum equation (Eqn. \ref{eq:governing:mom})

$p_{2} - p_{1} = ρ_{1} u_{1}^{2} - ρ_{2} u_{2}^{2} = {ρ_{1} u_{1} = ρ_{2} u_{1}} = ρ_{1} u_{1} (u_{1} - u_{2}) = ρ_{1} u_{1}^{2} (1 - \frac{u_{2}}{u_{1}})$

$\frac{p_{2} - p_{1}}{p_{1}} = \frac{ρ_{1} u_{1}^{2}}{p_{1}} (1 - \frac{u_{2}}{u_{1}}) = {a_{1} = \sqrt{\frac{γ p_{1}}{ρ_{1}}}} = γ \frac{u_{1}^{2}}{a_{1}^{2}} (1 - \frac{u_{2}}{u_{1}}) = γ M_{1}^{2} (1 - \frac{u_{2}}{u_{1}})$

$\frac{p_{2}}{p_{1}} - 1 = γ M_{1}^{2} (1 - \frac{u_{2}}{u_{1}}) = {\frac{u_{2}}{u_{1}} = \frac{ρ_{1}}{ρ_{2}}} = γ M_{1}^{2} (1 - \frac{2 + (γ - 1) M_{1}^{2}}{(γ + 1) M_{1}^{2}})$

$\frac{p_{2}}{p_{1}} = 1 + \frac{2 γ}{γ + 1} (M_{1}^{2} - 1)$

Figure~\ref{fig:shock:pressure:ratio} shows that the pressure must increase over the shock due to the fact that, based on the discussion above, the upstream Mach number must be greater than one and thus the shock is a discontinuous compression process.

The temperature ratio over the shock can be obtained using the already derived relations for pressure ratio and density ratio together with the equation of state $p = ρ R T$

$\frac{T_{2}}{T_{1}} = (\frac{p_{2}}{p_{1}}) (\frac{ρ_{1}}{ρ_{2}})$

$\frac{T_{2}}{T_{1}} = [1 + \frac{2 γ}{γ + 1} (M_{1}^{2} - 1)] [\frac{(γ + 1) M_{1}^{2}}{2 + (γ - 1) M_{1}^{2}}]$

Figure~\ref{fig:normal:shock:relations} below shows how different flow properties change over a normal shock as a function of upstream Mach number.

Now, one question remains. How come that we by analyzing the control volume using the upstream and downstream states get the normal shock relations. There is no way that the governing equations could have known about the fact that we assumed that there would be a shock inside of the control volume, or is it? The answer is that we have assumed that there will be a change in flow properties from upstream to downstream. We have further assumed that the flow is adiabatic (we are using the adiabatic energy equation) so there is no heat exchange. We are, however, allowing for irreversibilities in the flow. The only way to accomplish a change in flow properties under those constraints is a formation of a normal shock (a discontinuity in flow properties - a sudden flow compression) between station 1 and station 2.

One-dimensional inviscid flow: Difference between revisions

Latest revision as of 09:02, 14 March 2026

Contents

Acoustic Waves

Shock Waves

Normal Shock Relations

Navigation menu

@@ Line 1: / Line 1: @@
+[[Category:Compressible flow]]
+[[Category:One-dimensional flow]]
+[[Category:inviscid flow]]
 __TOC__
-== Acoustic Wave ==
+==Acoustic Waves==
 <!--
@@ Line 14: / Line 17: @@
 -->
-In Fig. \ref{fig:soundwave}, station 1 represents the flow state ahead of the sound wave and station 2 the flow state behind the sound wave. Set up the continuity equation for one-dimensional flows between 1 and 2. If we could change frame of reference and follow the sound wave, we would see fluid approaching the wave with the propagation speed of the wave, $a$, and behind the wave, the fluid would have a slightly modified speed, $a+da$. There would also be a slight in all other flow properties. Let's apply the one-dimensional continuity equation between station 1 and station 2.
+In Fig. \ref{fig:soundwave}, station 1 represents the flow state ahead of the sound wave and station 2 the flow state behind the sound wave. Set up the continuity equation for one-dimensional flows between 1 and 2. If we could change frame of reference and follow the sound wave, we would see fluid approaching the wave with the propagation speed of the wave, <math>a</math>, and behind the wave, the fluid would have a slightly modified speed, <math>a+da</math>. There would also be a slight in all other flow properties. Let's apply the one-dimensional continuity equation between station 1 and station 2.
-\[\rho_1 u_1=\rho_2 u_2\]\\
+<math display="block"> \rho_1 u_1=\rho_2 u_2</math>
-\[\rho a=(\rho+d\rho)(a+da)\]\\
+<math display="block">\rho a=(\rho+d\rho)(a+da)</math>
-\[\cancel{\rho a}=\cancel{\rho a} + \rho da + ad\rho +\underbrace{d\rho da}_{\sim 0} \Rightarrow\]\\
+<math display="block">\cancel{\rho a}=\cancel{\rho a} + \rho da + ad\rho +\underbrace{d\rho da}_{\sim 0} \Rightarrow</math>
-\begin{equation}
+<math display="block">
 a=-\rho\frac{da}{d\rho}
-\label{eq:speedofsound:a}
+</math>
-\end{equation}\\
-\noindent The one-dimensional momentum equation between station 1 and station 2 gives\\
+The one-dimensional momentum equation between station 1 and station 2 gives
-\[\rho_1 u_1^2+p_1=\rho_2 u_2^2+p_2\]\\
+<math display="block">\rho_1 u_1^2+p_1=\rho_2 u_2^2+p_2</math>
-\[\rho a^2+p=(\rho+d\rho)(a+da)^2+(p+dp)\]\\
+<math display="block">\rho a^2+p=(\rho+d\rho)(a+da)^2+(p+dp)</math>
-\[\cancel{\rho a^2}+\cancel{p}=\cancel{\rho a^2}+2\rho ada+\underbrace{\rho da^2}_{\sim 0}+d\rho a^2+\underbrace{2d\rho ada}_{\sim 0}+\underbrace{d\rho da^2}_{\sim 0}+\cancel{p}+dp \Rightarrow\]\\
+<math display="block">\cancel{\rho a^2}+\cancel{p}=\cancel{\rho a^2}+2\rho ada+\underbrace{\rho da^2}_{\sim 0}+d\rho a^2+\underbrace{2d\rho ada}_{\sim 0}+\underbrace{d\rho da^2}_{\sim 0}+\cancel{p}+dp \Rightarrow</math>
-\[dp=-2\rho ada-d\rho a^2 \Rightarrow\]\\
+<math display="block">dp=-2\rho ada-d\rho a^2 \Rightarrow</math>
-\[da=-\frac{dp+d\rho a^2}{2\rho a}=-\frac{d\rho}{2a\rho}\left(\frac{dp}{d\rho}+a^2\right) \Rightarrow\]\\
+<math display="block">da=-\frac{dp+d\rho a^2}{2\rho a}=-\frac{d\rho}{2a\rho}\left(\frac{dp}{d\rho}+a^2\right) \Rightarrow</math>
-\begin{equation}
+<math display="block">
 \frac{da}{d\rho}=-\frac{1}{2a\rho}\left(\frac{dp}{d\rho}+a^2\right)
-\label{eq:speedofsound:b}
+</math>
-\end{equation}\\
-\noindent Eqn. \ref{eq:speedofsound:b} in \ref{eq:speedofsound:a} gives\\
+Eqn. \ref{eq:speedofsound:b} in \ref{eq:speedofsound:a} gives
-\[a=\frac{1}{2a}\left(\frac{dp}{d\rho}+a^2\right) \Rightarrow \]\\
+<math display="block">a=\frac{1}{2a}\left(\frac{dp}{d\rho}+a^2\right) \Rightarrow</math>
-\begin{equation}
+<math display="block">
 a^2=\frac{dp}{d\rho}
-\label{eq:speedofsound:c}
+</math>
-\end{equation}\\
-\noindent Sound wave:
+Sound wave:
-\begin{itemize}
+* gradients are small
-\item gradients are small
+* irreversible (dissipative effects are negligible)
-\item irreversible (dissipative effects are negligible)
+* no heat addition
-\item no heat addition
-\end{itemize}
-\noindent Thus, the change of flow properties as the sound wave passes can be assumed to be an isentropic process\\
-\begin{equation}
+Thus, the change of flow properties as the sound wave passes can be assumed to be an isentropic process
+<math display="block">
 a^2=\left(\frac{dp}{d\rho}\right)_s
-\label{eq:speedofsound:c}
+</math>
-\end{equation}\\
-\begin{equation}
+<math display="block">
 a=\sqrt{\left(\frac{dp}{d\rho}\right)_s}=\sqrt{\frac{1}{\rho \tau_s}}
-\label{eq:speedofsound:d}
+</math>
-\end{equation}\\
-\noindent where $\tau_s$ is the compressibility of the gas. Eqn. \ref{eq:speedofsound:d} is valid for all gases. It can be seen from the equation, that truly incompressible flow ($\tau_s=0$) would imply infinite speed of sound. \\
+where <math>\tau_s</math> is the compressibility of the gas. Eqn. \ref{eq:speedofsound:d} is valid for all gases. It can be seen from the equation, that truly incompressible flow (<math>\tau_s=0</math>) would imply infinite speed of sound.
-\noindent Since the process is isentropic, we can use the isentropic relations if we also assume the gas to be calorically perfect\\
+Since the process is isentropic, we can use the isentropic relations if we also assume the gas to be calorically perfect\\
-\[\frac{p_2}{p_1}=\left(\frac{\rho_2}{\rho_1}\right)^\gamma \Rightarrow p=C\rho^\gamma\]\\
+<math display="block">\frac{p_2}{p_1}=\left(\frac{\rho_2}{\rho_1}\right)^\gamma \Rightarrow p=C\rho^\gamma</math>
-\[a^2=\left(\frac{dp}{d\rho}\right)_s=\gamma C\rho^{\gamma-1}=\gamma \underbrace{\left[C\rho^\gamma\right]}_{=p}\rho^{-1}=\frac{\gamma p}{\rho}\Rightarrow\]\\
+<math display="block">a^2=\left(\frac{dp}{d\rho}\right)_s=\gamma C\rho^{\gamma-1}=\gamma \underbrace{\left[C\rho^\gamma\right]}_{=p}\rho^{-1}=\frac{\gamma p}{\rho}\Rightarrow</math>
-\begin{equation}
+<math display="block">
 a=\sqrt{\frac{\gamma p}{\rho}}
-\label{eq:speedofsound:e}
+</math>
-\end{equation}\\
 or
-\begin{equation}
+<math display="block">
 a=\sqrt{\gamma RT}
-\label{eq:speedofsound:e}
+</math>
-\end{equation}\\
-\noindent From the relation above, it is obvious that the local speed of sound is related to the temperature of the flow, which in turn is a measure of the motion of elementary particles (atoms and/or molecules) of the fluid at a specific location. This stems from the fact that sound waves are propagated via interaction of these elementary particles. Since information in a flow is propagated via molecular interaction the relation between the speed at which this information is conveyed and the speed of the flow has important physical implications. Figure~\ref{fig:speed:of:sound} compares three sound wave patterns generated by a a beacon. In the left picture, the sound transmitter is stationary and thus the acoustic waves are centered around the transmitter. In the middle image, the transmitter is moving to the left at a speed less than the speed of sound and thus the transmitter will always be within all sound wave circles but it will be off-centered with a bias in the direction that the transmitter is moving. In the right image the transmitter is moving faster than the speed of sound and thus it will always be located outside of all acoustic waves. In a supersonic flow, no information can travel upstream and therefore there is no way for the flow to adjust to downstream obstacles. This is compensated for by the introduction of shocks in the flow. Over a shock flow properties changes discontinuity. An example is given in Figure~\ref{fig:supersonic:flow}.
+From the relation above, it is obvious that the local speed of sound is related to the temperature of the flow, which in turn is a measure of the motion of elementary particles (atoms and/or molecules) of the fluid at a specific location. This stems from the fact that sound waves are propagated via interaction of these elementary particles. Since information in a flow is propagated via molecular interaction the relation between the speed at which this information is conveyed and the speed of the flow has important physical implications. Figure~\ref{fig:speed:of:sound} compares three sound wave patterns generated by a a beacon. In the left picture, the sound transmitter is stationary and thus the acoustic waves are centered around the transmitter. In the middle image, the transmitter is moving to the left at a speed less than the speed of sound and thus the transmitter will always be within all sound wave circles but it will be off-centered with a bias in the direction that the transmitter is moving. In the right image the transmitter is moving faster than the speed of sound and thus it will always be located outside of all acoustic waves. In a supersonic flow, no information can travel upstream and therefore there is no way for the flow to adjust to downstream obstacles. This is compensated for by the introduction of shocks in the flow. Over a shock flow properties changes discontinuity. An example is given in Figure~\ref{fig:supersonic:flow}.
+<!--
 \begin{figure}[ht!]
 \begin{center}
@@ Line 111: / Line 107: @@
 \end{center}
 \end{figure}
+-->
+==Shock Waves==
+<!--
+\begin{figure}[ht!]
+\begin{center}
+\includegraphics[]{figures/standalone-figures/Chapter03/pdf/shock-wave.pdf}
+\caption{Stationary normal shock}
+\label{fig:shock}
+\end{center}
+\end{figure}
+-->
+The starting point is to set up the governing equations for one-dimensional steady compressible flow over a control volume enclosing the normal shock (Fig. \ref{fig:shock:cv}).
+<!--
+\begin{figure}[ht!]
+\begin{center}
+\includegraphics[]{figures/standalone-figures/Chapter03/pdf/shock-cv.pdf}
+\caption{Normal shock control volume}
+\label{fig:shock:cv}
+\end{center}
+\end{figure}
+-->
+continuity:
+<math display="block">
+\rho_1 u_1=\rho_2 u_2
+</math>
+momentum:
+<math display="block">
+\rho_1 u_1^2+p_1=\rho_2 u_2^2+p_2
+</math>
+energy:
+<math display="block">
+h_1 + \frac{1}{2}u_1^2=h_2 + \frac{1}{2}u_2^2
+</math>
+Divide the momentum equation by <math>\rho_1 u_1</math>
+<math display="block">
+\frac{1}{\rho_1 u_1}\left(\rho_1 u_1^2+p_1\right)=\frac{1}{\rho_1 u_1}\left(\rho_2 u_2^2+p_2\right)=\left\{\rho_1 u_1=\rho_2 u_2\right\}=\frac{1}{\rho_2 u_2}\left(\rho_2 u_2^2+p_2\right) \Rightarrow
+</math>
+<math display="block">
+\frac{p_1}{\rho_1 u_1}-\frac{p_2}{\rho_2 u_2}=u_2-u_1
+</math>
+For a calorically perfect gas <math>a=\sqrt{\gamma p/\rho}</math>, which if implemented in Eqn. \ref{eq:governing:mom:b} gives
+<math display="block">
+\frac{a_1^2}{\gamma u_1}-\frac{a_2^2}{\gamma u_2}=u_2-u_1
+</math>
+The energy equation (Eqn. \ref{eq:governing:energy}) with <math>h=C_p T</math>
+<math display="block">
+C_p T_1 + \frac{1}{2}u_1^2=C_p T_2 + \frac{1}{2}u_2^2
+</math>
+Replacing <math>C_p</math> with <math>\gamma R/(\gamma-1)</math> gives
+<math display="block">
+\frac{\gamma RT_1}{\gamma-1} + \frac{1}{2}u_1^2=\frac{\gamma RT_2}{\gamma-1} + \frac{1}{2}u_2^2
+</math>
+With <math>a=\sqrt{\gamma RT}</math> this becomes
+<math display="block">
+\frac{a_1^2}{\gamma-1} + \frac{1}{2}u_1^2=\frac{a_2^2}{\gamma-1} + \frac{1}{2}u_2^2
+</math>
+Eqn. \ref{eq:governing:energy:d} can be set up between any two points in the flow. Specifically, we can use the relation to relate the flow velocity, <math>u</math>, and speed of sound, <math>a</math>, in any point to the corresponding flow properties at sonic conditions (<math>u=a=a^*</math>).
+<math display="block">
+\frac{a^2}{\gamma-1} + \frac{1}{2}u^2=\frac{\gamma+1}{2(\gamma-1)}{a^*}^2
+</math>
+If Eqn. \ref{eq:governing:energy:e} is evaluated in locations 1 and 2, we get
+<math display="block">
+\begin{aligned}
+&a_1^2 = \frac{\gamma+1}{2}{a^*}^2 - \frac{\gamma-1}{2}u_1^2\\
+&a_2^2 = \frac{\gamma+1}{2}{a^*}^2 - \frac{\gamma-1}{2}u_2^2
+\end{aligned}
+</math>
+Since the change in flow conditions over the shock is adiabatic (no heat is added inside the shock), critical properties will be constant over the shock. Especially <math>a^*</math> will be constant.
+Eqn. \ref{eq:governing:energy:f} inserted in \ref{eq:governing:mom:c} gives\\
+<math display="block">
+\frac{1}{\gamma u_1}\left(\frac{\gamma+1}{2}{a^*}^2 - \frac{\gamma-1}{2}u_1^2\right)-\frac{1}{\gamma u_2}\left(\frac{\gamma+1}{2}{a^*}^2 - \frac{\gamma-1}{2}u_2^2\right)=u_2-u_1 \Rightarrow</math>
+<math display="block">
+\left(\frac{\gamma+1}{2\gamma}\right){a^*}^2\left(\frac{1}{u_1}-\frac{1}{u_2}\right)=\left(\frac{\gamma+1}{2\gamma}\right)\left(u_2-u_1\right) \Rightarrow
+</math>
+<math display="block">
+{a^*}^2\left(\frac{1}{u_1}-\frac{1}{u_2}\right)=\left(u_2-u_1\right) \Rightarrow
+</math>
+<math display="block">
+{a^*}^2\left(\frac{u_2}{u_1 u_2}-\frac{u_1}{u_1 u_2}\right)=\left(u_2-u_1\right) \Rightarrow
+</math>
+<math display="block">
+\frac{1}{u_1 u_2}{a^*}^2\left(u_2-u_1\right)=\left(u_2-u_1\right) \Rightarrow
+</math>
+<math display="block">
+{a^*}^2=u_1 u_2
+</math>
+Eqn. \ref{eq:prandtl} is sometimes referred to as the Prandtl relation. Divide the Prandtl relation by <math>{a^*}^2</math> on both sides gives
+<math display="block">
+=\frac{u_1}{a^*}\frac{u_2}{a^*}=M^*_1M^*_2
+</math>
+or
+<math display="block">
+M^*_2=\frac{1}{M^*_1}
+</math>
+The relation between <math>M^*</math> and <math>M</math> is given by
+<math display="block">
+{M^*}^2=\frac{(\gamma+1)M^2}{2+(\gamma-1)M^2}
+</math>
+from which is can be seen that <math>M^*</math> will follow the Mach number <math>M</math> in the sense that
+* <math>M=1\Rightarrow M^*=1</math>
+* <math>M<1\Rightarrow M^*<1</math>
+* <math>M>1\Rightarrow M^*>1</math>
+Eqn. \ref{eq:MachStar} inserted in Eqn. \ref{eq:NormalMach} gives
+<math display="block">
+\frac{(\gamma+1)M_1^2}{2+(\gamma-1)M_1^2}=\frac{2+(\gamma-1)M_2^2}{(\gamma+1)M_2^2}
+</math>
+<math display="block">
+M^2_2=\frac{1+\left[(\gamma-1)/2\right]M^2_1}{\gamma M^2_1-(\gamma-1)/2}
+</math>
+The Mach number relations above effectively show that if the Mach number upstream of the shock is greater than one, the downstream Mach number must be less than one and vice versa. We can also see that a sonic upstream flow <math>M_1=1.0</math> gives sonic flow downstream of the shock. So, apparently the relation as such holds for both supersonic and subsonic upstream flow mathematically. The question is if it is also physically correct. For a supersonic upstream flow we will get a discontinuous compression and if the flow upstream of the control volume is subsonic we will instead get a discontinuous expansion inside the control volume but, again, is this physically correct? We will get the answer by analyzing the entropy change over the control volume.
+Analyzing the energy equation and the second law of thermodynamics shows that there is a direct relation between entropy increase and total pressure drop.
+<math display="block">
+s_2-s_1=C_p\ln\dfrac{T_2}{T_1}-R\ln\dfrac{p_2}{p_1}
+</math>
+<math display="block">
+s_2-s_1=C_p\ln\dfrac{T_2}{T_{o_2}}\dfrac{T_{o_1}}{T_1}\dfrac{T_{o_2}}{T_{o_1}}-R\ln\dfrac{p_2}{p_{o_2}}\dfrac{p_{o_1}}{p_1}\dfrac{p_{o_2}}{p_{o_1}}
+</math>
+using the isentropic relations we get
+<math display="block">
+s_2-s_1=C_p\ln\dfrac{T_{o_2}}{T_{o_1}}-R\ln\dfrac{p_{o_2}}{p_{o_1}}
+</math>
+and since the process is adiabatic and thus <math>T_{o_2}=T_{o_1}</math> the change in entropy is directly related to the change in total pressure as
+<math display="block">
+s_2-s_1=-R\ln\dfrac{p_{o_2}}{p_{o_1}}
+</math>
+or
+<math display="block">
+\dfrac{p_{o_2}}{p_{o_1}}=e^{-(s_2-s_1)/R}
+</math>
+Figure~\ref{fig:shock:entropy} shows the entropy change over a normal shock. As can be seen in the figure, a subsonic upstream Mach number leads to a reduction of entropy, which once and for all rules out all such solutions as non-physical and thus the question about the upstream conditions can now be considered answered. This in turn implies that the Mach number downstream of a normal shock will always be subsonic, which can be seen in Fig~\ref{fig:shock:downstream:Mach} below.
+<!--
+\begin{figure}[ht!]
+\begin{center}
+\includegraphics[]{figures/standalone-figures/Chapter03/pdf/shock-entropy.pdf}
+\caption{Entropy change over a normal shock ($\Delta s$) as function of upstream Mach number ($M_1$)}
+\label{fig:shock:entropy}
+\end{center}
+\end{figure}
+\begin{figure}[ht!]
+\begin{center}
+\includegraphics[]{figures/standalone-figures/Chapter03/pdf/shock-downstream-Mach.pdf}
+\caption{Downstream Mach number ($M_2$) as function of upstream Mach number ($M_1$)}
+\label{fig:shock:downstream:Mach}
+\end{center}
+\end{figure}
+-->
+By rewriting the right-hand side of Eqn.\ref{eq:NormalMach:b}, it is easy to realize that the downstream Mach number <math>M_2</math> approaches a finite value for large values of the upstream Mach number, <math>M_1</math>.
+<math display="block">
+\left.M_2^2\right|_{M_1\rightarrow\infty}=\left.\dfrac{2/M_1^2+(\gamma-1)}{2\gamma-(\gamma-1)/M_1^2}\right|_{M_1\rightarrow\infty}=\dfrac{\gamma-1}{2\gamma}
+</math>
+===Normal Shock Relations===
+Rewriting the continuity equation (Eqn. \ref{eq:governing:cont})
+<math display="block">
+\frac{\rho_2}{\rho_1}=\frac{u_1}{u_2}=\frac{u_1^2}{u_1 u_2}=\left\{{a^*}^2=u_1u_2\right\}=\frac{u_1^2}{{a^*}^2}={M^*_1}^2
+</math>
+Eqn. \ref{eq:MachStar} in Eqn. \ref{eq:Normal:density:a} gives
+<math display="block">
+\frac{\rho_2}{\rho_1}=\frac{(\gamma+1)M_1^2}{2+(\gamma-1)M_1^2}
+</math>
+To get a corresponding relation for the pressure ratio  over the shock, we go back to the momentum equation (Eqn. \ref{eq:governing:mom})
+<math display="block">
+p_2-p_1=\rho_1 u^2_1 - \rho_2 u^2_2=\left\{\rho_1 u_1=\rho_2 u_1\right\}=\rho_1 u_1(u_1-u_2)=\rho_1 u^2_1\left(1-\frac{u_2}{u_1}\right)
+</math>
+<math display="block">
+\frac{p_2-p_1}{p_1}=\frac{\rho_1 u^2_1}{p_1}\left(1-\frac{u_2}{u_1}\right)=\left\{a_1=\sqrt{\frac{\gamma p_1}{\rho_1}}\right\}=\gamma\frac{u^2_1}{a^2_1}\left(1-\frac{u_2}{u_1}\right)=\gamma M^2_1\left(1-\frac{u_2}{u_1}\right)
+</math>
+<math display="block">
+\frac{p_2}{p_1}-1=\gamma M^2_1\left(1-\frac{u_2}{u_1}\right)=\left\{\frac{u_2}{u_1}=\frac{\rho_1}{\rho_2}\right\}=\gamma M^2_1\left(1-\frac{2+(\gamma-1)M_1^2}{(\gamma+1)M_1^2}\right)
+</math>
+<math display="block">
+\frac{p_2}{p_1}=1+\frac{2\gamma}{\gamma+1}(M^2_1-1)
+</math>
+Figure~\ref{fig:shock:pressure:ratio} shows that the pressure must increase over the shock due to the fact that, based on the discussion above, the upstream Mach number must be greater than one and thus the shock is a discontinuous compression process.
+<!--
+\begin{figure}[ht!]
+\begin{center}
+\includegraphics[]{figures/standalone-figures/Chapter03/pdf/shock-pressure-ratio.pdf}
+\caption{Pressure ratio over a normal shock ($p_2/p_1$) as function of upstream Mach number ($M_1$)}
+\label{fig:shock:pressure:ratio}
+\end{center}
+\end{figure}
+-->
+The temperature ratio over the shock can be obtained using the already derived relations for pressure ratio and density ratio together with the equation of state <math>p=\rho RT</math>
+<math display="block">
+\frac{T_2}{T_1}=\left(\frac{p_2}{p_1}\right)\left(\frac{\rho_1}{\rho_2}\right)
+</math>
+<math display="block">
+\frac{T_2}{T_1}=\left[1+\frac{2\gamma}{\gamma+1}(M^2_1-1)\right]\left[\frac{(\gamma+1)M_1^2}{2+(\gamma-1)M_1^2}\right]
+</math>
+Figure~\ref{fig:normal:shock:relations} below shows how different flow properties change over a normal shock as a function of upstream Mach number.
+<!--
+\begin{figure}[ht!]
+\begin{center}
+\includegraphics[]{figures/standalone-figures/Chapter03/pdf/normal-shock-relations.pdf}
+\caption{Changes of flow properties over a normal shock as a function of the upstream Mach number ($M_1$)}
+\label{fig:normal:shock:relations}
+\end{center}
+\end{figure}
+-->
+Now, one question remains. How come that we by analyzing the control volume using the upstream and downstream states get the normal shock relations. There is no way that the governing equations could have known about the fact that we assumed that there would be a shock inside of the control volume, or is it? The answer is that we have assumed that there will be a change in flow properties from upstream to downstream. We have further assumed that the flow is  adiabatic (we are using the adiabatic energy equation) so there is no heat exchange. We are, however, allowing for irreversibilities in the flow. The only way to accomplish a change in flow properties under those constraints is a formation of a normal shock (a discontinuity in flow properties - a sudden flow compression) between station 1 and station 2.

One-dimensional inviscid flow: Difference between revisions

Latest revision as of 09:02, 14 March 2026

Acoustic Waves

Shock Waves

Normal Shock Relations

Navigation menu

Search