|
– author: –
ŠKORPÍK, Jiří (LinkedIn.com/in/jiri-skorpik)
– issue date: –
January 2006, February 2023, February 2026 (3rd ed.)
– title: –
Mach number and high velocity flow effects
– web: –
– provenance: – Brno (Czech Republic)
– email: – skorpik.jiri@email.cz
Copyright©Jiří Škorpík, 2006-2026 |
Basic concepts in high velocity flowWhen fluids flow at speeds V significantly lower than the speed of sound a, we expect, as observers, a certain behavior of this fluid when it flows around bodies. A typical manifestation of subsonic flow is a smooth change in the flow trajectories already in front of the body being flowed around, as in Figure 520. This behavior is due to sound as a pressure disturbance propagating in a continuum from the source-S of pressure disturbance, which is the surface of the bodies. This means that in supersonic flow it is impossible for changes in pressure and other state quantities to propagate against the direction of flow. For this reason, we observe fundamental differences between the propagation of a pressure disturbance in a subsonic, sonic or supersonic flow in channels or when flowing around bodies and the formation of special phenomena (shock and expansion waves) that do not arise at subsonic speeds. The ratio of V to a is called the Mach number. – 520: – Subsonic flow around ![]() S-source of pressure disturbance.
– 337: – ![]() a [m·s-1] speed of sound in continuum under investigation; M [Mach] Mach number; p [Pa] pressure; r [J·kg-1·K-1] specific gas constant; T [K] absolute gas temperature (static temperature); V [m·s-1] velocity of body or flow; ρ [kg·m-3] density; κ [1] heat capacity ratio. The derivation of the equation for the speed of sound is shown in Appendix 337.
|
Effects of finite magnitude of speed of sound on stream continuityUsing sound as a pressure disturbance, information about the size of bodies in the investigated volume is propagated through the continuum. This allows, in subsonic flow, the streamlines to react to the flowed around bodies with an advance as in Figure 520, p. 3.3. On the contrary, in supersonic mutual motion of the fluid and the body, such an advance cannot occur and this body is forced to displace the surrounding fluid with its volume by a stepwise compression and thus change its direction through a effect called shock wave. The consequence of the finite size of the speed of sound is also a "mirror" change in properties of the channel flows when the speed of sound is exceeded.
– 772: – Propagation of sound waves when pressure disturbance source moves ![]() τ [s] time. The circles 0, 1, 2, 3 represent the boundary of sound waves in the environment at time τ=0...3. At time 0 the source is just at coordinate 0 at time 1 at coordinate 1 etc. This means that at point 0 the source will cause a pressure disturbance that propagates at the speed of sound in the spherical plane, after the source travels the distance 0-S it will have the radius of the sound wave indicated by the symbol 0 in the figure. The same procedure applies to the pressure disturbance induced by the source at point 1, etc. |
– 339: – ![]() (a) source is moving at supersonic speed; (b) streamlines around oblique shock wave (compare with Figure 520, p. 3.3); (c) source is moving at speed of sound. SW-shock wave. βSW [°] shock wave angle (μ<βSW); μ [°] Mach angle. The figure does not deal with the situation of the shock waves at the time before τ=0, nor with the situation behind the shock wave and behind the body; this problem is discussed in the next section of the paper.
|
– 518: – Hugoniot's theorem ![]() A [m] flow area. This equation is referred to as Hugoniot's theorem or the Area-Mach number relation. The flow tube can be formed by solid walls or a sharp boundary between two media with very different states or properties (liquid versus gas; rarefied gas versus shock wave at edge, etc.). The derivation of Hugoniot's theorem is given in Appendix 518.
|
Normal shock waveIn normal shock wave, the velocity decreases almost stepwise (wave thickness is approximately 10-7 m [Hloušek, 1992]) to subsonic without changing direction, as shown in Figure 519. Changes of state quantities in normal shock wave can be calculated using the Rankine-Hugoniot equations. Normal shock waves are generated in channels and around isolated bodies at the speed of sound. – 519: – Passage of gas through normal shock wave ![]() 1-gas state before shock wave; 2-gas state behind shock wave. p [Pa] pressure; V* [m·s-1] critical flow velocity. The derivation of the equations for normal shock wave is carried out, for example, in [Macur, 2010, p. 372].
|
|
– 338: – ![]() h [kJ·kg-1] gas enthalpy; s [J·kg-1·K-1] gas entropy; Lh [J·kg-1] shock wave loss; p* [Pa] critical pressure (pressure when flow hits speed of sound during expansion from point 1s); Vt [m·s-1] theoretical gas velocity during isentropic expansion from pressure p2s to pressure p1; ξ [1] shock wave loss coefficient. Index s denotes total state.
– 333: – Rankine-Hugoniot equations ![]() The equations are derived for the stable normal shock wave and ideal gas. The derivation of the equations is shown in Appendix 333. – Problem 896: –
A normal shock wave was generated in the Laval nozzle. Calculate the loss as the gas passes through this wave. The measured pressure and temperature in front of and behind the wave are shown in the attached figure. The calculated velocity before the wave from the nozzle cross section and mass flow rate is 583,72 m·s-1. Dry air flows through the nozzle. The solution of this problem is shown in Appendix 896.
![]() p [MPa]; t [°C]. |
Oblique shock wavesThe oblique shock wave is slope relative to the flow direction in front of it by the shock wave angle βSW. In front of the oblique shock wave, velocity must be supersonic, but behind, it can be subsonic or supersonic. As the flow traverses the oblique shock wave, its direction is deflected by angle δ (Figure 107). Normal velocity components (V1n, V2n) share properties with those passing through a normal shock wave (see Problem 1007). The equality of tangential velocity components (V1t=V2t) can be proved (e.g., [Kadrnožka, 2004, p. 126-127]). Oblique shock waves are formed, for example, around surfaces when gases flow at supersonic speeds. – 107: – Passage of gas through oblique shock wave ![]() δ [°] velocity deflection post-shock from original direction. The index n denotes the normal velocity components, the index t denotes tangential velocity components.
– Problem 1007: –
What is the angle of the shock wave formed by a missile at M=2,5 Mach? Find the velocity, temperature and pressure in the stream behind the wave? The geometry of the missile is shown in the figure. The other parameters are: κ=1,4, t1=20 °C, p1=101 325,25 Pa, r=287 J·kg-1·K-1. The solution of this problem is shown in Appendix 1007.
![]() |
– 808: – Formation of oblique shock wave at root of stepwise rising contact surface ![]() δS [°] angle of contact surface.
Unreachable compression wavesThe compression wave is an equivalent to the shock wave. It is smooth isentropic compression of supersonic flow in a narrowing space as described by Hugoniot's theorem. In practice, however, this process is not possible because the reduction of the flow area would have to be infinitesimal and the surfaces perfectly smooth [Dejč, 1967, p. 405]. In real supersonic flow, effect called crossed shock wave can be used to split shock wave. |
– 481: – ![]() (a) crossing of oblique shock waves on a stepped surface; (b) crossing of oblique shock waves on slowly rising surface (here the crossing is similar to compression wave). CW-set of compression waves.
– 905: – Project Quiet Spike ![]() The project successfully explored the possibility of reducing the intensity of sound effects by using a stepped extension of the nose of aircraft. Here, testing of the telescoping nose of F-15B aircraft. NASA Photo by:Carla Thomas
|
λ-shock waveThe λ-shock wave (Figure 865) forms when bodies are flown around at transonic speed (approximately between the critical Mach number and M=1..1,3) with a laminar boundary layer. In this layer, pressure gradually rises at the expense of velocity due to subsonic flow, increasing thickness. This creates a wedge, leading to oblique shock waves crossing, as shown in the simplified diagram (Figure 867). The resulting shock wave is often slightly inclined forward [Hošek, 1949]. In turbulent flow, the wedge is small (turbulent flow is less sensitive to pressure change), generating a directly normal shock wave at the boundary layer interface. – 865: – ![]() (a) overall view of λ-shock wave; (b) pressure change in λ-shock wave and in boundary layer. LBL-laminar boundary layer; i-pressure distribution in core of flow before and behind shock wave; ii-pressure distribution in laminar boundary layer; OSW-oblique shock waves due to increase in thickness of boundary layer. t [m] thickness of the boundary layer; x [m] distance.
– 867: – Principle of Boundary layer separation behind λ-shock wave ![]() CFD predictions of flow separation behind a shock wave are presented in the article [Marbona et al., 2026]. |
Expansion wavesIncreasing flow area also creates obtuse angles on bodies, like the trailing edge of projectiles, the beginning of fuselage taper, in nozzles etc., as seen in Figure 340, illustrating a typical supersonic obtuse angle flow. When wrapping around obtuse angles at supersonic velocity, gas must expand from pressure p1 to p2, flow velocity increases from V1 to V2, and the flow direction is deflected by angle δ from the original. In the expansion wave, there's a gradual change in state variables with very low losses (isentropic expansion). – 340: – Supersonic flow near obtuse angle ![]() ML-Mach line; Δ [m] distance traveled by fluid particle at velocity V before the pressure disturbance from source S reaches it.
– 521: – ![]() Prandtl-Meyer function for flow deflection: ν(M) [°] Prandtl-Meyer function
|
Effect of high velocities on airfoil aerodynamicsThe difference between subsonic and supersonic flow around airfoils is not only in the velocity distribution but mainly in the changes in values of aerodynamic quantities describing the force effects of the flow on the airfoil.
– 800: – ![]() (a) subsonic flow around airfoil; (b) transonic flow around airfoil, [Kneubuehl, 2004, p. 78]. EF-expansion fan; λW-λ-shock wave.
– 522: – ![]() (a) sound speed in front of airfoil; (b) supersonic speed in front of airfoil [Kneubuehl, 2004, p. 78].
|
|
– 897: – ![]() Photo source [O’Farrell and Rieckhoff, 2011].
– 894: – Profile types suitable for high speeds in compressible flow ![]() (a) transonic airfoil; (b) supersonic (lenticular shape); (c) supersonic (rhombic shape); (d) supersonic (trapezoidal shape); (e) hypersonic.
– 906: – Glauert-Prandtl rule ![]() (a) Glauert-Prandtl rule for pressure coefficient; (b) Glauert-Prandtl rule for lift coefficient; (c) effect of increasing velocity on drag coefficient. CD [1] drag coefficient of airfoil; CL [1] lift coefficient of airfoil; M [Mach] Mach number (before airfoil); CP [1] pressure coefficient of airfoil. The index i indicates incompressible flow, the index c compressible flow. The derivation is shown in [Hošek, 1949, p. 49]. |
– 893: – Changes in lift coefficient for rhombic airfoil ![]() Changes in lift coefficient for rhombic airfoil: i [°] angle of attack.
– 907: – Practical uses of Glauert-Prandtl rule ![]() (a) airfoil in incompressible flow; (b) airfoil in compressible flow. x [m] airfoil coordinates in direction of inflow velocity; yi, c [m] airfoil thickness.
|
Aerodynamics of profile cascades in compressible flowCompressible flow effects at high velocity also manifest in profile cascade, which has an impact on the resulting velocity triangles of turbomachines. The formation of these effects can be predicted by both numerical and analytical methods. However, there are supersonic wind tunnels with profile casacdes for measuring aerodynamic quantities and visualizing supersonic flow.
– 636: – ![]() left-scheme of situation recorded on interferogram; right-interferogram of supersonic flow through turbine profile cascade. Taken with Mach-Zehnder interferometer. Taken and images provided by Aerodynamic Laboratory in Nový Knín at Institute of Thermomechanics of AVČR, v.v.i.
|
|
– 635: – Examples of numerical modelling of compressible flow in blade cascade ![]() (a) water vapour M1=0,42 (in front of cascade), M2=0,7 (behind cascade), created at the Energy Institute of Brno University of Technology; (b) turbine blade cascade, working gas air [Tajč et al., 2007]. M [Mach]; ρ [kg-m-3] density.
ReferencesŠKORPÍK, Jiří, 2022, Essential equations of turbomachines, turbomachinery.education, Brno, https://turbomachinery.education/essential-equations-of-turbomachines.html.
ŠKORPÍK, Jiří, 2023, Technická matematika, engineering-sciences.education, Brno, https://engineering-sciences.education/technicka-matematika.html.
ŠKORPÍK, Jiří, 2024, Technická termomechanika, engineering-sciences.education, Brno, https://engineering-sciences.education/technicka-termomechanika.html.
ABBOTT, Ira, DOENHOFF, Albert, 1959, Theory of wing sections, including a summary of airfoil data, Dover publications, inc., New York, ISBN-10:0-486-60586-8.
ANON., 2010, Expansion fan – Isentropic flow, National Aeronautics and Space Administration – NASA, Washington, D.C., http://www.grc.nasa.gov.
DEJČ, Michail, 1967, Technická dynamika plynů, SNTL, Praha.
HLOUŠEK, Jiří, 1992, Termomechanika, VUT v Brně, Brno, ISBN 80-214-0387-X.
HOŠEK, Josef, 1949, Aerodynamika vysokých rychlostí, Naše vojsko, Praha.
HOŠEK, Josef, 1962, Dnešní letectví: nové vývojové směry v letectví, Práce, Praha
KADRNOŽKA, Jaroslav, 2004, Tepelné turbíny a turbokompresory, Akademické nakladatelství CERM, s.r.o., Brno, ISBN 80–7204–346–3.
KNEUBUEHL, Beat, 2004, Balistika střely, přesnost střelby, účinek, Naše vojsko, Praha, ISBN 80-206-0749-8.
MACUR, Milan, 2010, Úvod do analytické mechaniky a mechaniky kontinua, Vutium, Brno, ISBN 978-80-214-3944-3.
O’FARRELL, J.M., RIECKHOFF, T.J., 2011, Direct Visualization of Shock Waves in Supersonic Space Shuttle Flight, George C. Marshall Space Flight Center, AL 35812 – Technical Memorandum.
STEVER, Guyford, HAGGERTY James, 1966, Flight, Time Inc.
TAJČ, Ladislav, BEDNÁŘ, Lukáš, POLANSKÝ, Jiři, Šťastný, Miroslav, 2007, Radial Control Stage with Partial Steam Admission, Proceedings of the 8th International Symposium on Experimental and Computational Aerothermodynamics of Internal Flows, Lyon.
MARBONA, Himpu, RODRÍGUEZ, Daniel, MARTÍNEZ-CAVA AGUILAR, Alejandro, VALERO RAY VANDENHOECK, Eusebio, PURI, Kunal, PANIAGUA, Guillermo, 2026, Large-Eddy Simulations of A Turbulent Boundary Layer Separation Over a Bump With Strong Pressure Gradients and Free-Stream Turbulence, Journal of Fluids Engineering, ASME, New York, ISSN 0098-2202, doi: https://doi.org/10.1115/1.4071255.
©Jiří Škorpík, LICENCE
|