Prandtl–Meyer expansion fan: Difference between revisions

Content deleted Content added

Inline

Latest revision as of 16:47, 25 September 2023

A supersonic expansion fan, technically known as Prandtl–Meyer expansion fan, a two-dimensional simple wave, is a centered expansion process that occurs when a supersonic flow turns around a convex corner. The fan consists of an infinite number of Mach waves, diverging from a sharp corner. When a flow turns around a smooth and circular corner, these waves can be extended backwards to meet at a point.

Each wave in the expansion fan turns the flow gradually (in small steps). It is physically impossible for the flow to turn through a single "shock" wave because this would violate the second law of thermodynamics.^[1]

Across the expansion fan, the flow accelerates (velocity increases) and the Mach number increases, while the static pressure, temperature and density decrease. Since the process is isentropic, the stagnation properties (e.g. the total pressure and total temperature) remain constant across the fan.

The theory was described by Theodor Meyer on his thesis dissertation in 1908, along with his advisor Ludwig Prandtl, who had already discussed the problem a year before.^[2]^[3]

Flow properties

The expansion fan consists of an infinite number of expansion waves or Mach lines.^[4] The first Mach line is at an angle $\mu _{1}=\arcsin \left({\frac {1}{M_{1}}}\right)$ with respect to the flow direction, and the last Mach line is at an angle $\mu _{2}=\arcsin \left({\frac {1}{M_{2}}}\right)$ with respect to final flow direction. Since the flow turns in small angles and the changes across each expansion wave are small, the whole process is isentropic.^[1] This simplifies the calculations of the flow properties significantly. Since the flow is isentropic, the stagnation properties like stagnation pressure ( $p_{0}$ ), stagnation temperature ( $T_{0}$ ) and stagnation density ( $\rho _{0}$ ) remain constant. The final static properties are a function of the final flow Mach number ( $M_{2}$ ) and can be related to the initial flow conditions as follows, where $\gamma$ is the heat capacity ratio of the gas (1.4 for air):

{\begin{aligned}{\frac {T_{2}}{T_{1}}}&=\left({\frac {1+{\frac {\gamma -1}{2}}M_{1}^{2}}{1+{\frac {\gamma -1}{2}}M_{2}^{2}}}\right)\\[3pt]{\frac {p_{2}}{p_{1}}}&=\left({\frac {1+{\frac {\gamma -1}{2}}M_{1}^{2}}{1+{\frac {\gamma -1}{2}}M_{2}^{2}}}\right)^{\frac {\gamma }{\gamma -1}}\\[3pt]{\frac {\rho _{2}}{\rho _{1}}}&=\left({\frac {1+{\frac {\gamma -1}{2}}M_{1}^{2}}{1+{\frac {\gamma -1}{2}}M_{2}^{2}}}\right)^{\frac {1}{\gamma -1}}.\end{aligned}}

The Mach number after the turn ( $M_{2}$ ) is related to the initial Mach number ( $M_{1}$ ) and the turn angle ( $\theta$ ) by,

\theta =\nu (M_{2})-\nu (M_{1})\,

where, $\nu (M)\,$ is the Prandtl–Meyer function. This function determines the angle through which a sonic flow (M = 1) must turn to reach a particular Mach number (M). Mathematically,

{\begin{aligned}\nu (M)&=\int {\frac {\sqrt {M^{2}-1}}{1+{\frac {\gamma -1}{2}}M^{2}}}{\frac {\,dM}{M}}\\&={\sqrt {\frac {\gamma +1}{\gamma -1}}}\arctan {\sqrt {{\frac {\gamma -1}{\gamma +1}}\left(M^{2}-1\right)}}-\arctan {\sqrt {M^{2}-1}}.\\\end{aligned}}

By convention, $\nu (1)=0.\,$

Thus, given the initial Mach number ( $M_{1}$ ), one can calculate $\nu (M_{1})\,$ and using the turn angle find $\nu (M_{2})\,$ . From the value of $\nu (M_{2})\,$ one can obtain the final Mach number ( $M_{2}$ ) and the other flow properties. The velocity field in the expansion fan, expressed in polar coordinates $(r,\phi )$ are given by^[5]

v_{r}={\sqrt {2(h_{0}-h)-c^{2}}},\quad v_{\phi }=c,\quad {\text{where}}\quad \phi =-\int {\frac {d(\rho c)}{\rho {\sqrt {2(h_{0}-h)-c^{2}}}}},

$h$ is the specific enthalpy and $h_{0}$ is the stagnation specific enthalpy.

Maximum turn angle

As Mach number varies from 1 to $\infty$ , $\nu \,$ takes values from 0 to $\nu _{\text{max}}\,$ , where

\nu _{\text{max}}={\frac {\pi }{2}}\left({\sqrt {\frac {\gamma +1}{\gamma -1}}}-1\right).

This places a limit on how much a supersonic flow can turn through, with the maximum turn angle given by,

\theta _{\text{max}}=\nu _{\text{max}}-\nu (M_{1}).\,

One can also look at it as follows. A flow has to turn so that it can satisfy the boundary conditions. In an ideal flow, there are two kinds of boundary condition that the flow has to satisfy,

Velocity boundary condition, which dictates that the component of the flow velocity normal to the wall be zero. It is also known as no-penetration boundary condition.
Pressure boundary condition, which states that there cannot be a discontinuity in the static pressure inside the flow (since there are no shocks in the flow).

If the flow turns enough so that it becomes parallel to the wall, we do not need to worry about pressure boundary condition. However, as the flow turns, its static pressure decreases (as described earlier). If there is not enough pressure to start with, the flow won't be able to complete the turn and will not be parallel to the wall. This shows up as the maximum angle through which a flow can turn. The lower the Mach number is to start with (i.e. small $M_{1}$ ), the greater the maximum angle through which the flow can turn.

The streamline which separates the final flow direction and the wall is known as a slipstream (shown as the dashed line in the figure). Across this line there is a jump in the temperature, density and tangential component of the velocity (normal component being zero). Beyond the slipstream the flow is stagnant (which automatically satisfies the velocity boundary condition at the wall). In case of real flow, a shear layer is observed instead of a slipstream, because of the additional no-slip boundary condition.

Notes

^ ^a ^b
An expansion process through a single "shock" is impossible, because it will violate the second law of thermodynamics.

Impossibility of expanding a flow through a single "shock" wave: Consider the scenario shown in the adjacent figure. As a supersonic flow turns, the normal component of the velocity increases ( $w_{2}>w_{1}$ ), while the tangential component remains constant ( $v_{2}=v_{1}$ ). The corresponding change is the entropy ( $\Delta s=s_{2}-s_{1}$ ) can be expressed as follows,

${\begin{aligned}{\frac {\Delta s}{R}}&=\ln \left[\left({\frac {p_{2}}{p_{1}}}\right)^{\frac {1}{\gamma -1}}\left({\frac {\rho _{2}}{\rho _{1}}}\right)^{-{\frac {\gamma }{\gamma -1}}}\right]\\&\approx {\frac {\gamma +1}{12\gamma ^{2}}}\left({\frac {p_{2}-p_{1}}{p_{1}}}\right)^{3}\\&\approx {\frac {\gamma +1}{12\gamma ^{2}}}\left[{\frac {\rho _{1}w_{1}^{2}}{p_{1}}}\left(1-{\frac {w_{2}}{w_{1}}}\right)\right]^{3}\end{aligned}}$

where, $R$ is the universal gas constant, $\gamma$ is the ratio of specific heat capacities, $\rho$ is the static density, $p$ is the static pressure, $s$ is the entropy, and $w$ is the component of flow velocity normal to the "shock". The suffix "1" and "2" refer to the initial and final conditions respectively.
Since $w_{2}>w_{1}$ , this would mean that $\Delta s<0$ . Since this is not possible, it means that it is impossible to turn a flow through a single shock wave. The argument may be further extended to show that such an expansion process can occur only if we consider a turn through infinite number of expansion waves in the limit $\Delta s\rightarrow 0$ . Accordingly, an expansion process is an isentropic process.
^ Meyer, T. (1908). Über zweidimensionale Bewegungsvorgänge in einem Gas, das mit Überschallgeschwindigkeit strömt (Doctoral dissertation) (in German). Georg-August Universität, Göttingen. OCLC 77709738.
^ Prandtl, L. (1907). "Neue Untersuchungen über die strömende Bewegung der Gase und Dämpfe". Physikalische Zeitschrift (in German). 8: 23–30. Reprinted in Riegels, F. W., ed. (1961). Ludwig Prandtl Gesammelte Abhandlungen. Berlin: Springer. doi:10.1007/978-3-662-11836-8_78.
^
For an object moving at supersonic speeds ( $u>c$ ) as it moves from point A to B (distance u·t), the disturbances originating from point A travel a distance c·t. The corresponding angle is known as a Mach angle and the lines enclosing the disturbed region are known as Mach lines (in 2-D case) or Mach cone (in 3-D).
Mach lines (cone) and Mach angle:
Mach lines are a concept usually encountered in 2-D supersonic flows (i.e. $M\geq 1$ ). They are a pair of bounding lines which separate the region of disturbed flow from the undisturbed part of the flow. These lines occur in pairs and are oriented at an angle

$\mu =\arcsin \left({\frac {c}{u}}\right)=\arcsin \left({\frac {1}{M}}\right)$

with respect to the direction of motion (also known as the Mach angle). In case of 3-D flow field, these lines form a surface known as Mach cone, with Mach angle as the half angle of the cone.
To understand the concept better, consider the case sketched in the figure. We know that when an object moves in a flow, it causes pressure disturbances (which travel at the speed of sound, also known as Mach waves). The figure shows an object moving from point A to B along the line AB at supersonic speeds ( $u>c$ ). By the time the object reaches point B, the pressure disturbances from point A have travelled a distance c·t and are now at circumference of the circle (with centre at point A). There are infinite such circles with their centre on the line AB, each representing the location of the disturbances due to the motion of the object. The lines propagating outwards from point B and tangent to all these circles are known as Mach lines.
Note: These concepts have a physical meaning only for supersonic flows ( $u\geq c$ ). In case of subsonic flows, the disturbances will travel faster than the source and the argument of the $\arcsin()$ function will be greater than one.
^ Landau, L. D., & Lifshitz, E. M. (2013). Fluid Mechanics: Landau and Lifshitz: Course of Theoretical Physics, Volume 6 (Vol. 6). Elsevier.

References

Liepmann, Hans W.; Roshko, A. (2001) [1957]. Elements of Gasdynamics. Dover Publications. ISBN 0-486-41963-0.
Von Mises, Richard (2004) [1958]. Mathematical theory of compressible fluid flow. Dover Publications. ISBN 0-486-43941-0.
Courant, Richard; Friedrichs, K. O. (1999) [1948]. Supersonic flow and shock waves. Springer Science+Business Media. ISBN 0387902325.
Anderson, John D. Jr. (January 2001) [1984]. Fundamentals of Aerodynamics (3rd ed.). McGraw-Hill Science/Engineering/Math. ISBN 0-07-237335-0.
Shapiro, Ascher H. (1953). The Dynamics and Thermodynamics of Compressible Fluid Flow, Volume 1. Ronald Press. ISBN 978-0-471-06691-0.

External links

[shockExpansion-1] 
An expansion process through a single "shock" is impossible, because it will violate the second law of thermodynamics.

Impossibility of expanding a flow through a single "shock" wave: Consider the scenario shown in the adjacent figure. As a supersonic flow turns, the normal component of the velocity increases ( $w_{2}>w_{1}$ ), while the tangential component remains constant ( $v_{2}=v_{1}$ ). The corresponding change is the entropy ( $\Delta s=s_{2}-s_{1}$ ) can be expressed as follows,

${\begin{aligned}{\frac {\Delta s}{R}}&=\ln \left[\left({\frac {p_{2}}{p_{1}}}\right)^{\frac {1}{\gamma -1}}\left({\frac {\rho _{2}}{\rho _{1}}}\right)^{-{\frac {\gamma }{\gamma -1}}}\right]\\&\approx {\frac {\gamma +1}{12\gamma ^{2}}}\left({\frac {p_{2}-p_{1}}{p_{1}}}\right)^{3}\\&\approx {\frac {\gamma +1}{12\gamma ^{2}}}\left[{\frac {\rho _{1}w_{1}^{2}}{p_{1}}}\left(1-{\frac {w_{2}}{w_{1}}}\right)\right]^{3}\end{aligned}}$

where, $R$ is the universal gas constant, $\gamma$ is the ratio of specific heat capacities, $\rho$ is the static density, $p$ is the static pressure, $s$ is the entropy, and $w$ is the component of flow velocity normal to the "shock". The suffix "1" and "2" refer to the initial and final conditions respectively.
Since $w_{2}>w_{1}$ , this would mean that $\Delta s<0$ . Since this is not possible, it means that it is impossible to turn a flow through a single shock wave. The argument may be further extended to show that such an expansion process can occur only if we consider a turn through infinite number of expansion waves in the limit $\Delta s\rightarrow 0$ . Accordingly, an expansion process is an isentropic process.

[2] Meyer, T. (1908). Über zweidimensionale Bewegungsvorgänge in einem Gas, das mit Überschallgeschwindigkeit strömt (Doctoral dissertation) (in German). Georg-August Universität, Göttingen. OCLC 77709738.

[3] Prandtl, L. (1907). "Neue Untersuchungen über die strömende Bewegung der Gase und Dämpfe". Physikalische Zeitschrift (in German). 8: 23–30. Reprinted in Riegels, F. W., ed. (1961). Ludwig Prandtl Gesammelte Abhandlungen. Berlin: Springer. doi:10.1007/978-3-662-11836-8_78.

[4] 
For an object moving at supersonic speeds ( $u>c$ ) as it moves from point A to B (distance u·t), the disturbances originating from point A travel a distance c·t. The corresponding angle is known as a Mach angle and the lines enclosing the disturbed region are known as Mach lines (in 2-D case) or Mach cone (in 3-D).
Mach lines (cone) and Mach angle:
Mach lines are a concept usually encountered in 2-D supersonic flows (i.e. $M\geq 1$ ). They are a pair of bounding lines which separate the region of disturbed flow from the undisturbed part of the flow. These lines occur in pairs and are oriented at an angle

$\mu =\arcsin \left({\frac {c}{u}}\right)=\arcsin \left({\frac {1}{M}}\right)$

with respect to the direction of motion (also known as the Mach angle). In case of 3-D flow field, these lines form a surface known as Mach cone, with Mach angle as the half angle of the cone.
To understand the concept better, consider the case sketched in the figure. We know that when an object moves in a flow, it causes pressure disturbances (which travel at the speed of sound, also known as Mach waves). The figure shows an object moving from point A to B along the line AB at supersonic speeds ( $u>c$ ). By the time the object reaches point B, the pressure disturbances from point A have travelled a distance c·t and are now at circumference of the circle (with centre at point A). There are infinite such circles with their centre on the line AB, each representing the location of the disturbances due to the motion of the object. The lines propagating outwards from point B and tangent to all these circles are known as Mach lines.
Note: These concepts have a physical meaning only for supersonic flows ( $u\geq c$ ). In case of subsonic flows, the disturbances will travel faster than the source and the argument of the $\arcsin()$ function will be greater than one.

[5] Landau, L. D., & Lifshitz, E. M. (2013). Fluid Mechanics: Landau and Lifshitz: Course of Theoretical Physics, Volume 6 (Vol. 6). Elsevier.

[1]

[2]

[3]

[4]

[5]

@@ Line 1: / Line 1: @@
+{{Short description|Phenomenon in fluid dynamics}}
-[[Image:Expansion fan.png|thumb|300px|When a supersonic flow encounters a convex corner, it forms an expansion fan, which consists of an infinite number of expansion waves centered at the corner. The figure shows one such ideal expansion fan.]]
+{{No footnotes|date=November 2020}}
+[[Image:Expansion fan.png|thumb|300px|When a supersonic flow encounters a convex corner, it forms an expansion fan, which consists of an infinite number of expansion waves centred at the corner. The figure shows one such ideal expansion fan.]]
-A '''Prandtl–Meyer expansion fan''' is a centered expansion process, which turns a supersonic flow around a [[Convex set|convex]] corner. The fan consists of an infinite number of [[Mach wave]]s, diverging from a sharp corner. In case of a smooth corner, these waves can be extended backwards to meet at a point. Each wave in the expansion fan turns the flow gradually (in small steps). It is physically impossible to turn the flow away from itself through a single "shock" wave because it will violate the [[second law of thermodynamics]].<ref name="shockExpansion" >[[Image:Shock expansion.png|thumb|300px|An expansion process through a single "shock" is impossible, because it will violate the second law of thermodynamics.]] '''Impossibility of expanding a flow through a single "shock" wave:'''
+A supersonic expansion fan, technically known as '''Prandtl–Meyer expansion fan''', a two-dimensional [[simple wave]], is a centered expansion process that occurs when a supersonic flow turns around a [[Convex set|convex]] corner. The fan consists of an infinite number of [[Mach wave]]s, diverging from a sharp corner. When a flow turns around a smooth and circular corner, these waves can be extended backwards to meet at a point.
+Each wave in the expansion fan turns the flow gradually (in small steps). It is physically impossible for the flow to turn through a single "shock" wave because this would violate the [[second law of thermodynamics]].<ref name="shockExpansion">
+[[Image:Shock expansion.png|thumb|300px|An expansion process through a single "shock" is impossible, because it will violate the second law of thermodynamics.]]
+Impossibility of expanding a flow through a single "shock" wave:
-Consider the scenario shown in the adjacent figure. As a supersonic flow turns, the normal component of the velocity increases (<math> w_2 > w_1 </math>), while the tangential component remains constant (<math> v_2 = v_1 </math>). The corresponding change is the entropy (<math>\Delta s = s_2-s_1</math>) can be expressed as follows,
+Consider the scenario shown in the adjacent figure. As a supersonic flow turns, the normal component of the velocity increases (<math> w_2 > w_1 </math>), while the tangential component remains constant (<math> v_2 = v_1 </math>). The corresponding change is the entropy (<math>\Delta s = s_2 - s_1</math>) can be expressed as follows,
-:<math>\begin{align} \frac{\Delta s}{R}
+:<math>\begin{align}
+  \frac{\Delta s}{R}
-& = ln \bigg[ \bigg( \frac{p_2}{p_1} \bigg)^{1/(\gamma-1)} \bigg( \frac{\rho_2}{\rho_1} \bigg)^{-\gamma/(\gamma-1)} \bigg] \\
+    & = \ln \left[ \left( \frac{p_2}{p_1} \right)^\frac{1}{\gamma - 1} \left( \frac{\rho_2}{\rho_1} \right)^{-\frac{\gamma}{\gamma - 1}} \right] \\
-& \approx \frac{\gamma+1}{12\gamma^2} \bigg( \frac{p_2 - p_1}{p_1} \bigg)^3 \\
+    & \approx \frac{\gamma + 1}{12\gamma^2} \left( \frac{p_2 - p_1}{p_1} \right)^3 \\
-& \approx \frac{\gamma+1}{12\gamma^2} \bigg[ \frac{\rho_1 w_1^2}{p_1} \bigg( 1-\frac{w_2}{w_1}\bigg) \bigg]^3
+    & \approx \frac{\gamma + 1}{12\gamma^2} \left[ \frac{\rho_1 w_1^2}{p_1} \left(1 - \frac{w_2}{w_1}\right) \right]^3
 \end{align}</math>
-where, <math>R</math> is the universal gas constant, <math>\gamma</math> is the ratio of specific heat capacities, <math>\rho</math> is the static density, <math>p</math> is the static pressure, <math>s</math> is the static entropy, and  <math>w</math> is the component of flow velocity normal to the "shock". The suffix "1" and "2" refer to the initial and final conditions respectively.
+where, <math>R</math> is the universal gas constant, <math>\gamma</math> is the ratio of specific heat capacities, <math>\rho</math> is the static density, <math>p</math> is the static pressure, <math>s</math> is the entropy, and <math>w</math> is the component of flow velocity normal to the "shock". The suffix "1" and "2" refer to the initial and final conditions respectively.
-Since <math> w_2 > w_1 </math>, this would mean that <math>\Delta s < 0</math>. Since this is not possible it means that it is impossible to turn a flow through a single shock wave. The argument may be further extended to show that such an expansion process can occur only if we consider a turn through infinite number of expansion waves in the limit <math> \Delta s \rightarrow 0</math>. Accordingly an expansion process is an [[isentropic process]].
-</ref> Across the expansion fan, the flow accelerates (velocity increases) and the [[Mach number]] increases, while the [[static pressure]], [[temperature]] and [[density]] decrease. Since the process is [[isentropic]], the stagnation properties remain constant across the fan.
+Since <math> w_2 > w_1 </math>, this would mean that <math>\Delta s < 0</math>. Since this is not possible, it means that it is impossible to turn a flow through a single shock wave. The argument may be further extended to show that such an expansion process can occur only if we consider a turn through infinite number of expansion waves in the limit <math> \Delta s \rightarrow 0</math>. Accordingly, an expansion process is an [[isentropic process]].
+</ref>
+Across the expansion fan, the flow accelerates (velocity increases) and the [[Mach number]] increases, while the [[static pressure]], [[temperature]] and [[density]] decrease. Since the process is [[isentropic]], the [[stagnation point|stagnation]] properties (e.g. the total pressure and total temperature) remain constant across the fan.
+The theory was described by [[Theodor Meyer]] on his thesis dissertation in 1908, along with his advisor [[Ludwig Prandtl]], who had already discussed the problem a year before.<ref>{{cite thesis |last=Meyer |first=T. |year=1908 |title=Über zweidimensionale Bewegungsvorgänge in einem Gas, das mit Überschallgeschwindigkeit strömt |type=Doctoral dissertation |publisher=Georg-August Universität, Göttingen |oclc=77709738 |language=de }}</ref><ref>{{cite journal |last=Prandtl |first=L. |year=1907 |title=Neue Untersuchungen über die strömende Bewegung der Gase und Dämpfe |journal=[[Physikalische Zeitschrift]] |volume=8 |pages=23–30 |language=de }} Reprinted in {{cite book |editor-last=Riegels |editor-first=F. W. |title=Ludwig Prandtl Gesammelte Abhandlungen |publisher=Springer |location=Berlin |year=1961 |doi=10.1007/978-3-662-11836-8_78 }}</ref>
 == Flow properties ==
-The expansion fan consists of infinite number of expansion waves or [[Mach line]]s.<ref>
+The expansion fan consists of an infinite number of expansion waves or [[Mach line]]s.<ref>
-[[Image:Mach angle.jpg|thumb|300px|For an object moving at supersonic speeds (<math>u > c</math>) as it moves from point A to B (distance u·t), the disturbances originating from point A travel a distance c·t. The corresponding angle is known as Mach angle and the lines enclosing the disturbed region are known as Mach lines (in 2-D case) or Mach cone (in 3-D).]] '''Mach lines (cone) and Mach angle:'''
+[[Image:Mach angle.jpg|thumb|300px|For an object moving at supersonic speeds (<math>u > c</math>) as it moves from point A to B (distance u·t), the disturbances originating from point A travel a distance c·t. The corresponding angle is known as a Mach angle and the lines enclosing the disturbed region are known as Mach lines (in 2-D case) or Mach cone (in 3-D).]] '''Mach lines (cone) and Mach angle:'''
-'''Mach lines''' are a concept usually encountered in 2-D supersonic flows (i.e. <math>M \ge 1</math>). They are a pair of bounding lines which separate the region of disturbed flow from the undisturbed part of the flow. These lines occur in pair and are oriented at an angle
+'''Mach lines''' are a concept usually encountered in 2-D supersonic flows (i.e. <math>M \ge 1</math>). They are a pair of bounding lines which separate the region of disturbed flow from the undisturbed part of the flow. These lines occur in pairs and are oriented at an angle
-:<math> \textstyle \mu = \arcsin \left( \frac{c}{u} \right) = \arcsin \left( \frac{1}{M} \right)</math>
+:<math>\mu = \arcsin \left( \frac{c}{u} \right) = \arcsin \left( \frac{1}{M} \right)</math>
-with respect to the direction of motion (also known as the '''Mach angle'''). In case of 3-D flow field these lines form a surface known as '''Mach cone''', with Mach angle as the half angle of the cone.
+with respect to the direction of motion (also known as the '''Mach angle'''). In case of 3-D flow field, these lines form a surface known as '''Mach cone''', with Mach angle as the half angle of the cone.
-To understand the concept better, consider the case sketched in the figure. We know that when an object moves in a flow, it causes pressure disturbances (which travel at the speed of sound, also known as [[Mach wave]]s). The figure shows an object moving from point A to B along the line AB at supersonic speeds (<math>u > c</math>). By the time the object reaches point B, the pressure disturbances from point A have traveled a distance c·t and are now at circumference of the circle (with center at point A). There are infinite such circles with their center on the line AB, each representing the location of the disturbances due to the motion of the object. The lines propagating outwards from point B and tangent to all these circles are known as Mach lines.
+To understand the concept better, consider the case sketched in the figure. We know that when an object moves in a flow, it causes pressure disturbances (which travel at the speed of sound, also known as [[Mach wave]]s). The figure shows an object moving from point A to B along the line AB at supersonic speeds (<math>u > c</math>). By the time the object reaches point B, the pressure disturbances from point A have travelled a distance c·t and are now at circumference of the circle (with centre at point A). There are infinite such circles with their centre on the line AB, each representing the location of the disturbances due to the motion of the object. The lines propagating outwards from point B and tangent to all these circles are known as Mach lines.
-''Note:'' These concepts have a physical meaning only for supersonic flows (<math>u \ge c</math>). In case of subsonic flows the disturbances will travel faster than the source and the argument of the <math>\arcsin( )</math> function will be greater than one.
+''Note:'' These concepts have a physical meaning only for supersonic flows (<math>u \ge c</math>). In case of subsonic flows, the disturbances will travel faster than the source and the argument of the <math>\arcsin( )</math> function will be greater than one.
-</ref> The first Mach line is at an angle <math>\mu_1 = \arcsin \left( \frac{1}{M_1} \right)</math> with respect to the flow direction, and the last Mach line is at an angle <math>\mu_2 = \arcsin \left( \frac{1}{M_2} \right)</math> with respect to final flow direction. Since the flow turns in small angles and the changes across each expansion wave are small, the whole process is isentropic.<ref name="shockExpansion" /> This simplifies the calculations of the flow properties significantly. Since the flow is isentropic, the stagnation properties like stagnation pressure (<math>p_0</math>), [[stagnation temperature]] (<math>T_0</math>) and stagnation density (<math>\rho_0</math>) remain constant. The final static properties are a function of the final flow Mach number (<math>M_2</math>) and can be related to the initial flow conditions as follows,
+</ref> The first Mach line is at an angle <math>\mu_1 = \arcsin \left( \frac{1}{M_1} \right)</math> with respect to the flow direction, and the last Mach line is at an angle <math>\mu_2 = \arcsin \left( \frac{1}{M_2} \right)</math> with respect to final flow direction. Since the flow turns in small angles and the changes across each expansion wave are small, the whole process is isentropic.<ref name="shockExpansion" /> This simplifies the calculations of the flow properties significantly. Since the flow is isentropic, the [[stagnation point|stagnation]] properties like [[stagnation pressure]] (<math>p_0</math>), [[stagnation temperature]] (<math>T_0</math>) and stagnation density (<math>\rho_0</math>) remain constant. The final static properties are a function of the final flow Mach number (<math>M_2</math>) and can be related to the initial flow conditions as follows, where <math>\gamma</math> is the [[heat capacity ratio]] of the gas (1.4 for air):
-:<math>\begin{array}{lcl}
-\frac{T_2}{T_1} &= & \bigg( \frac{1+\frac{\gamma -1}{2}M_1^2}{1+\frac{\gamma -1}{2}M_2^2} \bigg) \\
+:<math>\begin{align}
-\frac{p_2}{p_1} &= & \bigg( \frac{1+\frac{\gamma -1}{2}M_1^2}{1+\frac{\gamma -1}{2}M_2^2} \bigg)^{\gamma/(\gamma-1)} \\
-\frac{\rho_2}{\rho_1} &= &\bigg( \frac{1+\frac{\gamma -1}{2}M_1^2}{1+\frac{\gamma -1}{2}M_2^2} \bigg)^{1/(\gamma-1)}.
+  \frac{T_2}{T_1} &= \left( \frac{1 + \frac{\gamma - 1}{2}M_1^2}{1 + \frac{\gamma - 1}{2}M_2^2} \right) \\[3pt]
+  \frac{p_2}{p_1} &= \left( \frac{1 + \frac{\gamma - 1}{2}M_1^2}{1 + \frac{\gamma - 1}{2}M_2^2} \right)^\frac{\gamma}{\gamma - 1} \\[3pt]
-\end{array}</math>
+  \frac{\rho_2}{\rho_1} &= \left( \frac{1 + \frac{\gamma - 1}{2}M_1^2}{1 + \frac{\gamma - 1}{2}M_2^2} \right)^\frac{1}{\gamma - 1}.
+\end{align}</math>
 The Mach number after the turn (<math>M_2</math>) is related to the initial Mach number (<math>M_1</math>) and the turn angle (<math>\theta</math>) by,
@@ Line 39: / Line 52: @@
 where, <math>\nu(M) \,</math> is the [[Prandtl–Meyer function]]. This function determines the angle through which a sonic flow ([[Mach number|M]] = 1) must turn to reach a particular Mach number (M). Mathematically,
-: <math>\begin{align} \nu(M)
+: <math>\begin{align}
+  \nu(M)
-& = \int \frac{\sqrt{M^2-1}}{1+\frac{\gamma -1}{2}M^2}\frac{\,dM}{M} \\
+    &= \int \frac{\sqrt{M^2 - 1}}{1 + \frac{\gamma - 1}{2}M^2}\frac{\,dM}{M} \\
-& = \sqrt{\frac{\gamma + 1}{\gamma -1}} \cdot \arctan \sqrt{\frac{\gamma -1}{\gamma +1} (M^2 -1)} - \arctan \sqrt{M^2 -1}. \\
+    &= \sqrt{\frac{\gamma + 1}{\gamma - 1}} \arctan \sqrt{\frac{\gamma - 1}{\gamma + 1} \left(M^2 - 1\right)} - \arctan\sqrt{M^2 - 1}. \\
 \end{align} </math>
 By convention, <math>\nu(1) = 0. \,</math>
-Thus, given the initial Mach number (<math> M_1 </math>), one can calculate <math>\nu(M_1) \,</math> and using the turn angle find <math>\nu(M_2) \,</math>. From the value of <math>\nu(M_2) \,</math> one can obtain the final Mach number (<math> M_2 </math>) and the other flow properties.
+Thus, given the initial Mach number (<math> M_1 </math>), one can calculate <math>\nu(M_1) \,</math> and using the turn angle find <math>\nu(M_2) \,</math>. From the value of <math>\nu(M_2) \,</math> one can obtain the final Mach number (<math> M_2 </math>) and the other flow properties. The velocity field in the expansion fan, expressed in polar coordinates <math>(r,\phi)</math> are given by<ref>Landau, L. D., & Lifshitz, E. M. (2013). Fluid Mechanics: Landau and Lifshitz: Course of Theoretical Physics, Volume 6 (Vol. 6). Elsevier.</ref>
+:<math>v_r = \sqrt{2(h_0-h)-c^2}, \quad v_\phi=c, \quad \text{where} \quad \phi = - \int \frac{d(\rho c)}{\rho \sqrt{2(h_0-h)-c^2}},</math>
+<math>h</math> is the specific enthalpy and <math>h_0</math> is the stagnation specific enthalpy.
 == Maximum turn angle ==
-[[Image:Maximum turning angle.png|thumb|300px|There is a limit on the maximum angle (<math>\theta_{max}</math>) through which a supersonic flow can turn.]]
+[[Image:Maximum turning angle.png|thumb|300px|There is a limit on the maximum angle (<math>\theta_\text{max}</math>) through which a supersonic flow can turn.]]
 As Mach number varies from 1 to <math>\infty</math>, <math>\nu \,</math> takes values from 0 to <math>\nu_\text{max} \,</math>, where
-: <math>\nu_\text{max} = \frac{\pi}{2} \bigg( \sqrt{\frac{\gamma+1}{\gamma-1}} -1 \bigg).</math>
+: <math>\nu_\text{max} = \frac{\pi}{2} \left( \sqrt{\frac{\gamma + 1}{\gamma - 1}} - 1 \right).</math>
 This places a limit on how much a supersonic flow can turn through, with the maximum turn angle given by,
@@ Line 59: / Line 78: @@
 One can also look at it as follows. A flow has to turn so that it can satisfy the boundary conditions. In an ideal flow, there are two kinds of boundary condition that the flow has to satisfy,
 # Velocity boundary condition, which dictates that the component of the flow velocity [[normal component|normal]] to the wall be zero. It is also known as no-penetration boundary condition.
 # Pressure boundary condition, which states that there cannot be a discontinuity in the static pressure inside the flow (since there are no shocks in the flow).
 If the flow turns enough so that it becomes parallel to the wall, we do not need to worry about pressure boundary condition. However, as the flow turns, its static pressure decreases (as described earlier). If there is not enough pressure to start with, the flow won't be able to complete the turn and will not be parallel to the wall. This shows up as the maximum angle through which a flow can turn. The lower the Mach number is to start with (i.e. small <math>M_1</math>), the greater the maximum angle through which the flow can turn.
@@ Line 80: / Line 99: @@
 * {{cite book
   | last = Liepmann | first = Hans W.  |author2=Roshko, A.
-  | title = Elements of Gasdynamics    | origyear = 1957
+  | title = Elements of Gasdynamics    | orig-year = 1957
   | publisher = [[Dover Publications]] | year = 2001
   | isbn = 0-486-41963-0 }}
+* {{cite book
+  | last = Von Mises | first = Richard
+  | title = Mathematical theory of compressible fluid flow    | orig-year = 1958
+  | publisher = [[Dover Publications]] | year = 2004
+  | isbn = 0-486-43941-0 }}
+* {{cite book | last = Courant | first = Richard | author2 = Friedrichs, K. O. | title = Supersonic flow and shock waves | orig-year = 1948 | publisher = [[Springer Science+Business Media]] | year = 1999 | isbn = 0387902325 | url-access = registration | url = https://archive.org/details/supersonicflowsh0000cour }}
 * {{cite book
   | last = Anderson | first = John D. Jr.
-  | title = Fundamentals of Aerodynamics | origyear = 1984 | edition = 3rd Edition
+  | title = Fundamentals of Aerodynamics | orig-year = 1984 | edition = 3rd
   | publisher = [[McGraw-Hill|McGraw-Hill Science/Engineering/Math]] |date=January 2001
   | isbn = 0-07-237335-0 }}
@@ Line 95: / Line 120: @@
 == External links ==
-* [http://exploration.grc.nasa.gov/education/rocket/expans.html Expansion fan] ([[NASA]])
+* [https://web.archive.org/web/20130216144918/http://exploration.grc.nasa.gov/education/rocket/expans.html Expansion fan] ([[NASA]])
-* [http://meweb.ecn.purdue.edu/~meapplet/java/comp_calculator/expansionfan_calculator.html Prandtl- Meyer expansion fan calculator] ([[Java applet]]).
+* [https://web.archive.org/web/20070205180438/http://meweb.ecn.purdue.edu/~meapplet/java/comp_calculator/expansionfan_calculator.html Prandtl- Meyer expansion fan calculator] ([[Java applet]]).
 {{DEFAULTSORT:Prandtl-Meyer expansion fan}}
 [[Category:Aerodynamics]]
+[[Category:Conservation equations]]
 [[Category:Fluid dynamics]]