Expansion waves: Difference between revisions
From Flowpedia
No edit summary |
No edit summary |
||
| Line 1: | Line 1: | ||
[[Category:Compressible flow]] | <!-- | ||
[[Category:Two-dimensional flow]] | -->[[Category:Compressible flow]]<!-- | ||
[[Category:Inviscid flow]] | -->[[Category:Two-dimensional flow]]<!-- | ||
<noinclude> | -->[[Category:Inviscid flow]]<!-- | ||
[[Category:Compressible flow:Topic]] | --><noinclude><!-- | ||
</noinclude> | -->[[Category:Compressible flow:Topic]]<!-- | ||
--></noinclude><!-- | |||
__TOC__ | --><nomobile><!-- | ||
-->__TOC__<!-- | |||
--></nomobile><!-- | |||
--><noinclude><!-- | |||
-->{{#vardefine:secno|4}}<!-- | |||
-->{{#vardefine:eqno|0}}<!-- | |||
--></noinclude><!-- | |||
--> | |||
=== Prandtl-Meyer Expansion Waves === | === Prandtl-Meyer Expansion Waves === | ||
| Line 22: | Line 30: | ||
A single Mach wave has a insignificant effect on the flow passing it but an expansion region constitutes an infinite number of Mach waves and the integrated effect is significant. The net turning of the flow by a single Mach wave is depicted schematically in Fig. \ref{fig:machwave}. It can be shown geometrically that | A single Mach wave has a insignificant effect on the flow passing it but an expansion region constitutes an infinite number of Mach waves and the integrated effect is significant. The net turning of the flow by a single Mach wave is depicted schematically in Fig. \ref{fig:machwave}. It can be shown geometrically that | ||
<math | {{NumEqn|<math> | ||
d\theta=\sqrt{M^2-1}\frac{dV}{V} | d\theta=\sqrt{M^2-1}\frac{dV}{V} | ||
</math> | </math>}} | ||
Since Eqn. \ref{eq:mach:turning} is derived from the flow turning geometry with the assumption that the net flow tuning is small, it is valid for all gas models. | Since Eqn. \ref{eq:mach:turning} is derived from the flow turning geometry with the assumption that the net flow tuning is small, it is valid for all gas models. | ||
| Line 30: | Line 38: | ||
To get the integrated effect of all Mach waves in the expansion region, we integrate Eqn. \ref{eq:mach:turning} over the expansion region | To get the integrated effect of all Mach waves in the expansion region, we integrate Eqn. \ref{eq:mach:turning} over the expansion region | ||
<math | {{NumEqn|<math> | ||
\int_{\theta_1}^{\theta_2}d\theta=\int_{M_1}^{M_2}\sqrt{M^2-1}\frac{dV}{V} | \int_{\theta_1}^{\theta_2}d\theta=\int_{M_1}^{M_2}\sqrt{M^2-1}\frac{dV}{V} | ||
</math> | </math>}} | ||
To be able to do the integration, we need to rewrite it | To be able to do the integration, we need to rewrite it | ||
<math | {{NumEqn|<math> | ||
V=Ma \Rightarrow \ln V=\ln M + \ln a | V=Ma \Rightarrow \ln V=\ln M + \ln a | ||
</math> | </math>}} | ||
Differentiate to get | Differentiate to get | ||
<math | {{NumEqn|<math> | ||
\frac{dV}{V}=\frac{dM}{M}+\frac{da}{a} | \frac{dV}{V}=\frac{dM}{M}+\frac{da}{a} | ||
</math> | </math>}} | ||
Each Mach wave is isentropic and thus the expansion is an isentropic process, which means that we can use the adiabatic energy equation | Each Mach wave is isentropic and thus the expansion is an isentropic process, which means that we can use the adiabatic energy equation | ||
<math | {{NumEqn|<math> | ||
\frac{T_o}{T}=1+\frac{\gamma-1}{2}M^2 | \frac{T_o}{T}=1+\frac{\gamma-1}{2}M^2 | ||
</math> | </math>}} | ||
For a calorically perfect gas <math>a=\sqrt{\gamma RT}</math> and <math>a_o=\sqrt{\gamma RT_o}</math> and thus | For a calorically perfect gas <math>a=\sqrt{\gamma RT}</math> and <math>a_o=\sqrt{\gamma RT_o}</math> and thus | ||
<math | {{NumEqn|<math> | ||
\frac{T_o}{T}=\left(\frac{a_o}{a}\right)^2\Rightarrow | \frac{T_o}{T}=\left(\frac{a_o}{a}\right)^2\Rightarrow | ||
</math> | </math>}} | ||
<math | {{NumEqn|<math> | ||
\left(\frac{a_o}{a}\right)^2=1+\frac{\gamma-1}{2}M^2 | \left(\frac{a_o}{a}\right)^2=1+\frac{\gamma-1}{2}M^2 | ||
</math> | </math>}} | ||
Solve for <math>a</math> gives | Solve for <math>a</math> gives | ||
<math | {{NumEqn|<math> | ||
a=a_o \left(1+\frac{\gamma-1}{2}M^2\right)^{-1/2} | a=a_o \left(1+\frac{\gamma-1}{2}M^2\right)^{-1/2} | ||
</math> | </math>}} | ||
Differentiate Eqn \ref{eq:adiatbatic:energy:b} to get | Differentiate Eqn \ref{eq:adiatbatic:energy:b} to get | ||
<math | {{NumEqn|<math> | ||
da=a_o\left(-\frac{1}{2}\right)\left(\frac{\gamma-1}{2}\right)2M\left(1+\frac{\gamma-1}{2}M^2\right)^{-3/2}dM | da=a_o\left(-\frac{1}{2}\right)\left(\frac{\gamma-1}{2}\right)2M\left(1+\frac{\gamma-1}{2}M^2\right)^{-3/2}dM | ||
</math> | </math>}} | ||
Eqn. \ref{eq:adiatbatic:energy:b} in Eqn. \ref{eq:adiatbatic:energy:c} gives | Eqn. \ref{eq:adiatbatic:energy:b} in Eqn. \ref{eq:adiatbatic:energy:c} gives | ||
<math | {{NumEqn|<math> | ||
\frac{da}{a}=-\left(\frac{\gamma-1}{2}\right)\left(1+\frac{\gamma-1}{2}M^2\right)^{-1}MdM | \frac{da}{a}=-\left(\frac{\gamma-1}{2}\right)\left(1+\frac{\gamma-1}{2}M^2\right)^{-1}MdM | ||
</math> | </math>}} | ||
From Eqn. \ref{eq:mach:turning:c}, we have | From Eqn. \ref{eq:mach:turning:c}, we have | ||
<math | {{NumEqn|<math> | ||
\frac{dV}{V}=\frac{dM}{M}+\frac{da}{a} | \frac{dV}{V}=\frac{dM}{M}+\frac{da}{a} | ||
</math> | </math>}} | ||
With <math>da/a</math> from Eqn. \ref{eq:adiatbatic:energy:d}, we get | With <math>da/a</math> from Eqn. \ref{eq:adiatbatic:energy:d}, we get | ||
<math | {{NumEqn|<math> | ||
\frac{dV}{V}=\frac{dM}{M}-\left(\frac{\gamma-1}{2}\right)\left(1+\frac{\gamma-1}{2}M^2\right)^{-1}MdM | \frac{dV}{V}=\frac{dM}{M}-\left(\frac{\gamma-1}{2}\right)\left(1+\frac{\gamma-1}{2}M^2\right)^{-1}MdM | ||
</math> | </math>}} | ||
<math | {{NumEqn|<math> | ||
\frac{dV}{V}=dM\left[\frac{\left(1+\dfrac{\gamma-1}{2}M^2\right)-\left(\dfrac{\gamma-1}{2}\right)M^2}{\left(1+\dfrac{\gamma-1}{2}M^2\right)M}\right]=\left(1+\dfrac{\gamma-1}{2}M^2\right)^{-1}\frac{dM}{M} | \frac{dV}{V}=dM\left[\frac{\left(1+\dfrac{\gamma-1}{2}M^2\right)-\left(\dfrac{\gamma-1}{2}\right)M^2}{\left(1+\dfrac{\gamma-1}{2}M^2\right)M}\right]=\left(1+\dfrac{\gamma-1}{2}M^2\right)^{-1}\frac{dM}{M} | ||
</math> | </math>}} | ||
Now, insert <math>dV/V</math> in Eqn. \ref{eq:mach:turning:b} to get | Now, insert <math>dV/V</math> in Eqn. \ref{eq:mach:turning:b} to get | ||
<math | {{NumEqn|<math> | ||
\int_{\theta_1}^{\theta_2}d\theta=\int_{M_1}^{M_2}\sqrt{M^2-1}\left(1+\dfrac{\gamma-1}{2}M^2\right)^{-1}\frac{dM}{M} | \int_{\theta_1}^{\theta_2}d\theta=\int_{M_1}^{M_2}\sqrt{M^2-1}\left(1+\dfrac{\gamma-1}{2}M^2\right)^{-1}\frac{dM}{M} | ||
</math> | </math>}} | ||
The integral on the right hand side of Eqn. \ref{eq:mach:turning:c} is the Prandtl-Meyer function, which is usually denoted <math>\nu</math>. The Prandtl-Meyer function evaluated for Mach number <math>M</math> becomes | The integral on the right hand side of Eqn. \ref{eq:mach:turning:c} is the Prandtl-Meyer function, which is usually denoted <math>\nu</math>. The Prandtl-Meyer function evaluated for Mach number <math>M</math> becomes | ||
<math | {{NumEqn|<math> | ||
\nu(M)=\sqrt{\frac{\gamma+1}{\gamma-1}}\tan^{-1}\sqrt{\frac{\gamma-1}{\gamma+1}(M^2-1)}-\tan^{-1}\sqrt{M^2-1} | \nu(M)=\sqrt{\frac{\gamma+1}{\gamma-1}}\tan^{-1}\sqrt{\frac{\gamma-1}{\gamma+1}(M^2-1)}-\tan^{-1}\sqrt{M^2-1} | ||
</math> | </math>}} | ||
and thus the net turning of the flow can be calculated as | and thus the net turning of the flow can be calculated as | ||
<math | {{NumEqn|<math> | ||
\theta_2-\theta_1=\nu(M_2)-\nu(M_1) | \theta_2-\theta_1=\nu(M_2)-\nu(M_1) | ||
</math> | </math>}} | ||
==== Solving Problems using the Prandtl Meyer Function ==== | ==== Solving Problems using the Prandtl Meyer Function ==== | ||
| Line 138: | Line 146: | ||
The aim is to derive relations of temperature, pressure and density over an expansion wave. Total temperature upstream and downstream of the expansion wave is calculated as | The aim is to derive relations of temperature, pressure and density over an expansion wave. Total temperature upstream and downstream of the expansion wave is calculated as | ||
<math | {{NumEqn|<math> | ||
\frac{T_{o_1}}{T_1}=1+\frac{\gamma-1}{2}M_1^2 | \frac{T_{o_1}}{T_1}=1+\frac{\gamma-1}{2}M_1^2 | ||
</math> | </math>}} | ||
<math | {{NumEqn|<math> | ||
\frac{T_{o_2}}{T_2}=1+\frac{\gamma-1}{2}M_2^2 | \frac{T_{o_2}}{T_2}=1+\frac{\gamma-1}{2}M_2^2 | ||
</math> | </math>}} | ||
The temperature ratio over the expansion wave may now be calculated as | The temperature ratio over the expansion wave may now be calculated as | ||
<math | {{NumEqn|<math> | ||
\frac{T_2}{T_1}=\frac{T_2}{T_{o_2}}\frac{T_{o_1}}{T_1}=\frac{1+\dfrac{\gamma-1}{2}M_1^2}{1+\dfrac{\gamma-1}{2}M_2^2}=\frac{2+(\gamma-1)M_1^2}{2+(\gamma-1)M_2^2} | \frac{T_2}{T_1}=\frac{T_2}{T_{o_2}}\frac{T_{o_1}}{T_1}=\frac{1+\dfrac{\gamma-1}{2}M_1^2}{1+\dfrac{\gamma-1}{2}M_2^2}=\frac{2+(\gamma-1)M_1^2}{2+(\gamma-1)M_2^2} | ||
</math> | </math>}} | ||
The expansion is isentropic and thus total temperature and both total pressure are unaffected by the expansion. Therefore, <math>T_{o_1}=T_{o_2}</math> and thus | The expansion is isentropic and thus total temperature and both total pressure are unaffected by the expansion. Therefore, <math>T_{o_1}=T_{o_2}</math> and thus | ||
<math | {{NumEqn|<math> | ||
\frac{T_2}{T_1}=\frac{2+(\gamma-1)M_1^2}{2+(\gamma-1)M_2^2} | \frac{T_2}{T_1}=\frac{2+(\gamma-1)M_1^2}{2+(\gamma-1)M_2^2} | ||
</math> | </math>}} | ||
The pressure and density ratios can be obtained from Eqn. \ref{eq:tr} using the isentropic relations | The pressure and density ratios can be obtained from Eqn. \ref{eq:tr} using the isentropic relations | ||
<math | {{NumEqn|<math> | ||
\frac{p_2}{p_1}=\left[\frac{2+(\gamma-1)M_1^2}{2+(\gamma-1)M_2^2}\right]^{\gamma/(\gamma-1)} | \frac{p_2}{p_1}=\left[\frac{2+(\gamma-1)M_1^2}{2+(\gamma-1)M_2^2}\right]^{\gamma/(\gamma-1)} | ||
</math> | </math>}} | ||
<math | {{NumEqn|<math> | ||
\frac{\rho_2}{\rho_1}=\left[\frac{2+(\gamma-1)M_1^2}{2+(\gamma-1)M_2^2}\right]^{1/(\gamma-1)} | \frac{\rho_2}{\rho_1}=\left[\frac{2+(\gamma-1)M_1^2}{2+(\gamma-1)M_2^2}\right]^{1/(\gamma-1)} | ||
</math> | </math>}} | ||
<!-- | <!-- | ||
