11.4.3.1: The pressure Mach number relationship
- Page ID
- 800
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Before going further in the mathematical derivations it is worth looking at the physical meaning of equation (42). The term \({\rho\, U^{2} / A}\) is always positive (because all the three terms can be only positive). Now, it can be observed that \(dP\) can be positive or negative depending on the \(dA\) and Mach number. The meaning of the sign change for the pressure differential is that the pressure can increase or decrease. It can be observed that the critical Mach number is one. If the Mach number is larger than one than \(dP\) has opposite sign of \(dA\). If Mach number is smaller than one \(dP\) and \(dA\) have the same sign. For the subsonic branch \(M<1\) the term \(1 /(1 -M^2)\) is positive hence
\[
dA > 0 \Longrightarrow dP> 0 \\
dA < 0 \Longrightarrow dP< 0
\]
\begin{align*}
dA > 0 \Rightarrow dP< 0 \\
dA < 0 \Rightarrow dP> 0
\end{align*}
This behavior is opposite to incompressible flow behavior. For the special case of \(M=1\) (sonic flow) the value of the term \(1 - M^{2}= 0\) thus mathematically \(dP \rightarrow \infty\) or \(dA = 0\). Since physically \(dP\) can increase only in a finite amount it must that \(dA = 0\). It must also be noted that when \(M=1\) occurs only when \(dA = 0\). However, the opposite, not necessarily means that when \(dA = 0\) that \(M=1\). In that case, it is possible that \(dM=0\) thus the diverging side is in the subsonic branch and the flow isn't choked. The relationship between the velocity and the pressure can be observed from equation (36) by solving it for \(dU\).
\[
dU = - \dfrac{ dP }{ P \,U}
\label{gd:iso:eq:relationshipU-rho}
\]
From equation (43) it is obvious that \(dU\) has an opposite sign to \(dP\) (since the term \(PU\) is positive). Hence the pressure increases when the velocity decreases and vice versa. From the speed of sound, one can observe that the density, \(\rho\), increases with pressure and vice versa (see equation (44)).
\[
d\rho = \dfrac{ 1 }{ c^2} dP
\label{gd:iso:eq:relationshipP-rho}
\]
It can be noted that in the derivations of the above equations (??) - (??), the equation of state was not used. Thus, the equations are applicable for any gas (perfect or imperfect gas). The second law (isentropic relationship) dictates that \(ds=0\) and from thermodynamics \(ds = 0 = C_p \,\dfrac{dT }{ T} - R \,\dfrac{dP }{ P}\) and for perfect gas
\[
\dfrac{dT }{ T} = \dfrac{k - 1 }{ k } \; \dfrac{dP }{ P}
\label{gd:iso:eq:isentropicDiff}
\]
Thus, the temperature varies in the same way that pressure does. The relationship between the Mach number and the temperature can be obtained by utilizing the fact that the process is assumed to be adiabatic \(dT_0 = 0\). Differentiation of equation (27), the relationship between the temperature and the stagnation temperature becomes
\[
dT_0 = 0 = dT \left( 1 + \dfrac{ k - 1 }{ 2 } M^2 \right) + T (k -1) M dM
\label{gd:iso:eq:temp-M1}
\]
and simplifying equation (46) yields
\[
\dfrac{ dT }{ T } = - \dfrac{(k -1) \,M\, dM }{ 1 + \dfrac{ k - 1 }{ 2 }\, M^{2} }
\label{gd:iso:eq:temp-M}
\]
Contributors and Attributions
Dr. Genick Bar-Meir. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or later or Potto license.