Compressible flow (gas dynamics) is the branch of fluid mechanics that deals with flows having significant changes in fluid density. Gases, but not liquids, display such behaviour.^{[1]} To distinguish between compressible and incompressible flow in air, the Mach number (the ratio of the speed of the flow to the speed of sound) must be greater than about 0.3 (since the density change is greater than 5% in that case) before significant compressibility occurs. The study of compressible flow is relevant to highspeed aircraft, jet engines, rocket motors, hyperloops, highspeed entry into a planetary atmosphere, gas pipelines, commercial applications such as abrasive blasting, and many other fields.
Contents

History 1

Introductory concepts 2

Mach number, wave motion, and sonic speed 3

Onedimensional flow 4

Convergingdiverging Laval nozzles 4.1

Maximum achievable velocity of a gas 4.2

Isentropic flow Mach number relationships 4.3

Achieving supersonic flow 4.4

Nonisentropic 1D channel flow of a gas  normal shock waves 4.5

Twodimensional flow 5

Oblique shock waves 5.1

Shock polar diagram 5.1.1

Oblique shock reflection 5.1.2

Solid boundary 5.1.2.1

Irregular reflection 5.1.2.2

PrandtlMeyer fans 5.2

PrandtlMeyer expansion fans 5.2.1

PrandtlMeyer compression fans 5.2.2

Applications 6

Supersonic wind tunnels 6.1

Supersonic aircraft inlets 6.2

Natural gas pipeline 6.3

See also 7

References 8

External links 9
History
The study of gas dynamics is often associated with the flight of modern highspeed aircraft and atmospheric reentry of spaceexploration vehicles; however, its origins lie with simpler machines. At the beginning of the 19th century, investigation into the behaviour of fired bullets led to improvement in the accuracy and capabilities of guns and artillery.^{[2]} As the century progressed, inventors such as Gustaf de Laval advanced the field, while researchers such as Ernst Mach sought to understand the physical phenomenon involved through experimentation.
At the beginning of the 20th century, the focus of gas dynamics research shifted to what would eventually become the aerospace industry. Ludwig Prandtl and his students proposed important concepts ranging from the boundary layer to supersonic shock waves, supersonic wind tunnels, and supersonic nozzle design.^{[2]} Theodore von Kármán, a student of Prandtl, continued to improve the understanding of supersonic flow. Other notable figures (Meyer, Luigi Crocco, and Shapiro) also contributed significantly to the principles considered fundamental to the study of modern gas dynamics. Many others also contributed to this field.
Accompanying the improved conceptual understanding of gas dynamics in the early 20th century was a public misconception that there existed a barrier to the attainable speed of aircraft, commonly referred to as the "sound barrier." In truth, the barrier to supersonic flight was merely a technological one, although it was a stubborn barrier to overcome. Amongst other factors, conventional aerofoils saw a dramatic increase in drag coefficient when the flow approached the speed of sound. Overcoming the larger drag proved difficult with contemporary designs, thus the perception of a sound barrier. However, aircraft design progressed sufficiently to produce the Bell X1. Piloted by Chuck Yeager, the X1 officially achieved supersonic speed in October 1947.^{[3]} This achievement paved the way to the future of modern aircraft, missiles, and spacecraft.
Historically, two parallel paths of research have been followed in order to further gas dynamics knowledge. Experimental gas dynamics undertakes wind tunnel model experiments and experiments in shock tubes and ballistic ranges with the use of optical techniques to document the findings. Theoretical gas dynamics considers the equations of motion applied to a variabledensity gas, and their solutions. Much of basic gas dynamics is analytical, but in the modern era Computational fluid dynamics applies computing power to solve the otherwiseintractable nonlinear partial differential equations of compressible flow for specific geometries and flow characteristics.
Introductory concepts
Breakdown of fluid mechanics chart
There are several important assumptions involved in the underlying theory of compressible flow. All fluids are composed of molecules, but tracking a huge number of individual molecules in a flow (for example at atmospheric pressure) is unnecessary. Instead, the continuum assumption allows us to consider a flowing gas as a continuous substance except at low densities. This assumption provides a huge simplification which is accurate for most gasdynamic problems. Only in the lowdensity realm of rarefied gas dynamics does the motion of individual molecules become important.
A related assumption is the noslip condition where the flow velocity at a solid surface is presumed equal to the velocity of the surface itself, which is a direct consequence of assuming continuum flow. The noslip condition implies that the flow is viscous, and as a result a boundary layer forms on bodies traveling through the air at high speeds, much as it does in lowspeed flow.
Most problems in incompressible flow involve only two unknowns: pressure and velocity, which are typically found by solving the two equations that describe conservation of mass and of linear momentum, with the fluid density presumed constant. In compressible flow, however, the gas density and temperature also become variables. This requires two more equations in order to solve compressibleflow problems: an equation of state for the gas and a conservation of energy equation. For the majority of gasdynamic problems, the simple Ideal gas law is the appropriate state equation.
Fluid dynamics problems have two overall types of references frames, called Lagrangian and Eulerian (see JosephLouis Lagrange and Leonhard Euler). The Lagrangian approach follows a fluid mass of fixed identity as it moves through a flowfield. The Eulerian reference frame, in contrast, does not move with the fluid. Rather it is a fixed frame or control volume that fluid flows through. The Eulerian frame is most useful in a majority of compressible flow problems, but requires that the equations of motion be written in a compatible format.
Finally, although space is known to have 3 dimensions, an important simplification can be had in describing gas dynamics mathematically if only one spatial dimension is of primary importance, hence 1dimensional flow is assumed. This works well in duct, nozzle, and diffuser flows where the flow properties change mainly in the flow direction rather than perpendicular to the flow. However, an important class of compressible flows, including the external flow over bodies traveling at high speed, requires at least a 2dimensional treatment. When all 3 spatial dimensions and perhaps the time dimension as well are important, we often resort to computerized solutions of the governing equations.
Mach number, wave motion, and sonic speed
The Mach number (M) is defined as the ratio of the speed of an object (or of a flow) to the speed of sound. For instance, in air at room temperature, the speed of sound is about 340 m/s (1,100 ft/s). M can range from 0 to ∞, but this broad range falls naturally into several flow regimes. These regimes are subsonic, transonic, supersonic, hypersonic, and hypervelocity flow. The figure below illustrates the Mach number "spectrum" of these flow regimes.
Mach number flow regimes spectrum
These flow regimes are not chosen arbitrarily, but rather arise naturally from the strong mathematical background that underlies compressible flow (see the cited reference textbooks). At very slow flow speeds the speed of sound is so much faster that it is mathematically ignored, and the Mach number is irrelevant. Once the speed of the flow approaches the speed of sound, however, the Mach number becomes allimportant, and shock waves begin to appear. Thus the transonic regime is described by a different (and much more difficult) mathematical treatment. In the supersonic regime the flow is dominated by wave motion at oblique angles similar to the Mach angle. Above about Mach 5, these wave angles grow so small that a different mathematical approach is required, defining the Hypersonic speed regime. Finally, at speeds comparable to that of planetary atmospheric entry from orbit, in the range of several km/s, the speed of sound is now comparatively so slow that it is once again mathematically ignored in the Hypervelocity regime.
As an object accelerates from subsonic toward supersonic speed in a gas, different types of wave phenomena occur. To illustrate these changes, the next figure shows a stationary point (M = 0) that emits symmetric sound waves. The speed of sound is the same in all directions in a uniform fluid, so these waves are simply concentric spheres. As the soundgenerating point begins to accelerate, the sound waves "bunch up" in the direction of motion and "stretch out" in the opposite direction. When the point reaches sonic speed (M = 1), it travels at the same speed as the sound waves it creates. Therefore, an infinite number of these sound waves "pile up" ahead of the point, forming a Shock wave. Upon achieving supersonic flow, the particle is moving so fast that it continuously leaves its sound waves behind. When this occurs, the locus of these waves trailing behind the point creates an angle known as the Mach wave angle or Mach angle, µ:

\mu = \arcsin(a/V) = \arcsin(1/M)
where a represents the speed of sound in the gas and V represents the velocity of the object. Although named for Austrian physicist Ernst Mach, these oblique waves were actually first discovered by Christian Doppler.^{[4]}
Wave motion and the speed of sound
Onedimensional flow
Onedimensional (1D) flow refers to the flow of gas through a duct or channel in which the flow parameters are assumed to change significantly along only one spatial dimension, namely, the duct length. In analysing the 1D channel flow, a number of assumptions are made:
Convergingdiverging Laval nozzles
As the speed of a flow accelerates from the subsonic to the supersonic regime, the physics of nozzle and diffuser flows is altered. Using the conservation laws of fluid dynamics and thermodynamics, the following relationship for channel flow is developed (combined mass and momentum conservation):

dP(1{M}^2) = \rho{V}^2\left(\frac {dA} {A}\right) ,
where dP is the differential change in pressure, M is the Mach number, ρ is the density of the gas, V is the velocity of the flow, A is the area of the duct, and dA is the change in area of the duct. This equation states that, for subsonic flow, a converging duct (dA<0) increases the velocity of the flow and a diverging duct (dA>0) decreases velocity of the flow. For supersonic flow, the opposite occurs due to the change of sign of (1M^{2}). A converging duct (dA<0) now decreases the velocity of the flow and a diverging duct (dA>0) increases the velocity of the flow. At Mach = 1, a special case occurs in which the duct area must be either a maximum or minimum. For practical purposes, only a minimum area can accelerate flows to Mach 1 and beyond. See table of subsupersonic diffusers and nozzles.
Table showing the reversal in the physics of nozzles and diffusers with changing Mach numbers
Therefore, to accelerate a flow to Mach 1, a nozzle must be designed to converge to a minimum crosssectional area and then expand. This type of nozzle – the convergingdiverging nozzle – is called a de Laval nozzle after Gustaf de Laval, who invented it. As subsonic flow enters the converging duct and the area decreases, the flow accelerates. Upon reaching the minimum area of the duct, also known as the throat of the nozzle, the flow can reach Mach 1. If the speed of the flow is to continue to increase, its density must decrease in order to obey conservation of mass. To achieve this decrease in density, the flow must expand, and to do so, the flow must pass through a diverging duct. See image of de Laval Nozzle.
Nozzle de Laval diagram
Maximum achievable velocity of a gas
Ultimately, because of the energy conservation law, a gas is limited to a certain maximum velocity based on its energy content. The maximum velocity, V_{max}, that a gas can attain is:

V_{max} = \sqrt {2c_{p} T_{t}}
where c_{p} is the specific heat of the gas and T_{t} is the stagnation temperature of the flow.
Isentropic flow Mach number relationships
Using conservations laws and thermodynamics, a number of relationships of the form

\frac=f(M,\gamma)
can be obtained, where M is the Mach number and γ is the ratio of specific heats (1.4 for air). See table of isentropic flow Mach number relationships.
Isentropic flow relationship table. Equations to relate the field properties in isentropic flow.
Achieving supersonic flow
As previously mentioned, in order for a flow to become supersonic, it must pass through a duct with a minimum area, or sonic throat. Additionally, an overall pressure ratio, P_{b}/P_{t}, of approximately 2 is needed to attain Mach 1. Once it has reached Mach 1, the flow at the throat is said to be "choked." Because changes downstream can only move upstream at sonic speed, the mass flow through the nozzle cannot be affected by changes in downstream conditions after the flow is choked.
Nonisentropic 1D channel flow of a gas  normal shock waves
Normal shock waves are shock waves that are perpendicular to the local flow direction. These shock waves occur when pressure waves build up and coalesce into an extremely thin shockwave that converts useful energy into heat. The waves thus overtake and reinforce one another, forming a finite shock wave from an infinite series of infinitesimal sound waves. Because a loss of energy occurs over the thin shock wave, the shock is considered nonisentropic and enthalpy increases across the shock. When analysing a normal shock wave, onedimensional, steady, and adiabatic (stagnation temperature does not change across the shock wave) flow of a perfect gas is assumed.
The RankineHugoniot equations relate conditions before and after a normal shock wave.
Normal shock waves can occur in two reference frames: the standing normal shock and the moving shock. The flow before a normal shock wave must be supersonic, and the flow after a normal shock must be subsonic. The RankineHugoniot equations are used to solve for the flow conditions.
Twodimensional flow
Although onedimensional flow can be directly analysed, it is merely a specialized case of twodimensional flow. It follows that one of the defining phenomena of onedimensional flow, a normal shock, is likewise only a special case of a larger class of oblique shocks. Further, the name "normal" is with respect to geometry rather than frequency of occurrence. Oblique shocks are much more common in applications such as: aircraft inlet design, objects in supersonic flight, and (at a more fundamental level) supersonic nozzles and diffusers. Depending on the flow conditions, an oblique shock can either be attached to the flow or detached from the flow in the form of a bow shock.
Attached shock wave shown on a X15 Model in a supersonic wind tunnel


Bowshock example for a blunt body

Oblique shock waves
Diagram of obstruction
Oblique shock waves are similar to normal shock waves, but they occur at angles less than 90° with the direction of flow. When a disturbance is introduced to the flow at a nonzero angle (δ), the flow must respond to the changing boundary conditions. Thus an oblique shock is formed, resulting in a change in the direction of the flow.
Shock polar diagram
Shock polar diagram
Based on the level of flow deflection (δ), oblique shocks are characterized as either strong or weak. Strong shocks are characterized by larger deflection and more entropy loss across the shock, with weak shocks as the opposite. In order to gain cursory insight into the differences in these shocks, a shock polar diagram can be used. With the static temperature after the shock, T*, known the speed of sound after the shock is defined as,

{A}^*= \sqrt{\gamma R{T}^*}
with R as the gas constant and γ as the specific heat ratio. The Mach number can be broken into Cartesian coordinates

{M}^*_{2x} = \frac

{M}^*_{2y} = \frac
with V_{x} and V_{y} as the x and ycomponents of the fluid velocity V. With the Mach number before the shock given, a locus of conditions can be specified. At some δ _{max} the flow transitions from a strong to weak oblique shock. With δ = 0°, a normal shock is produced at the limit of the strong oblique shock and the Mach wave is produced at the limit of the weak shock wave.
Oblique shock reflection
Due to the inclination of the shock, after an oblique shock is created, it can interact with a boundary in three different manners, two which are explained below.
Solid boundary
Incoming flow is first turned by angle δ with respect to the flow. This shockwave is reflected off the solid boundary, and the flow is turned by – δ to again be parallel with the boundary. It is important to note that each progressive shock wave is weaker and the wave angle is increased.
Irregular reflection
An irregular reflection is much like the case described above, with the caveat that δ is larger than the maximum allowable turning angle. Thus a detached shock is formed and a more complicated reflection occurs.
PrandtlMeyer fans
PrandtlMeyer fans can be expressed as both compression and expansion fans. PrandtlMeyer fans also cross a boundary layer (i.e. flowing and solid) which reacts in different changes as well. When a shock wave hits a solid surface the resulting fan returns as one from the opposite family while when one hits a free boundary the fan returns as a fan of opposite type.
PrandtlMeyer expansion fans
PrandtlMeyer expansion fan diagram
To this point, the only flow phenomena that have been discussed are shock waves, which slow the flow and increase its entropy. It is possible to accelerate supersonic flow in what has been termed a PrandtlMeyer expansion fan, after Ludwig Prandtl and Theodore Meyer. The mechanism for the expansion is shown in the figure below.
As opposed to the flow encountering an inclined obstruction and forming an oblique shock, the flow expands around a convex corner and forms an expansion fan through a series of isentropic Mach waves. The expansion "fan" is composed of Mach waves that span from the initial Mach angle to the final Mach angle. Flow can expand around either a sharp or rounded corner equally, as the increase in Mach number is proportional to only the convex angle of the passage (δ). The expansion corner that produces the PrandtlMeyer fan can be sharp (as illustrated in the figure) or rounded. If the total turning angle is the same, then the PM flow solution is also the same.
The PrandtlMeyer expansion can be seen as the physical explanation of the operation of the Laval nozzle. The contour of the nozzle creates a smooth and continuous series of PrandtlMeyer expansion waves.
PrandtlMeyer compression fans
Basic PM compression diagram
A PrandtlMeyer compression is the opposite phenomenon to a PrandtlMeyer expansion. If the flow is gradually turned through an angle of δ, a compression fan can be formed. This fan is a series of Mach waves that eventually coalesce into an oblique shock. Because the flow is defined by an isentropic region (flow that travels through the fan) and an anisentropic region (flow that travels through the oblique shock), a slip line results between the two flow regions.
Applications
Supersonic wind tunnels
Supersonic wind tunnels are used for testing and research in supersonic flows, approximately over the Mach number range of 1.2 to 5. The operating principle behind the wind tunnel is that a large pressure difference is maintained upstream to downstream, driving the flow.
Supersonic wind tunnel classification list
Wind tunnels can be divided into two categories: continuousoperating and intermittentoperating wind tunnels. Continuous operating supersonic wind tunnels require an independent electrical power source that drastically increases with the size of the test section. Intermittent supersonic wind tunnels are less expensive in that they store electrical energy over an extended period of time, then discharge the energy over a series of brief tests. The difference between these two is analogous to the comparison between a battery and a capacitor.
Blowdown supersonic wind tunnel schematic

Langley indraft supersonic wind tunnel vacuum sphere

Blowdown type supersonic wind tunnels offer high Reynolds number, a small storage tank, and readily available dry air. However, they cause a high pressure hazard, result in difficulty holding a constant stagnation pressure, and are noisy during operation.
Indraft supersonic wind tunnels are not associated with a pressure hazard, allow a constant stagnation pressure, and are relatively quiet. Unfortunately, they have a limited range for the Reynolds number of the flow and require a large vacuum tank. There is no dispute that knowledge is gained through research and testing in supersonic wind tunnels; however, the facilities often require vast amounts of power to maintain the large pressure ratios needed for testing conditions. For example, Arnold Engineering Development Complex has the largest supersonic wind tunnel in the world and requires the power required to light a small city for operation. For this reason, large wind tunnels are becoming less common at universities.
Supersonic aircraft inlets
Perhaps the most common application for oblique shocks is in highspeed aircraft inlets. The purpose of the inlet is to slow incoming supersonic flow to the subsonic regime before it enters the turbojet engine, with the caveat of minimizing losses across the shock. Knowledge of normal and oblique shocks suggests that this be accomplished with a series of weakening oblique shocks followed by a very weak normal shock, usually less than M = 1.4. This may sound relatively straight forward, but there is one rather large issue to be dealt with when designing a supersonic aircraft inlet: acceleration. Between taking off, manoeuvring, and cruising, an aircraft travels at a range of Mach numbers. In order to ensure efficient flight, the aircraft intake must be capable of variable geometry. If it is not, the shock waves will not reflect properly through the inlet and negatively affect performance.
McDonnell Douglas F15C

Concorde

Lockheed SR71 Blackbird

Although variable geometry is a universally recognized approach to improve aircraft efficiency and performance over a range of Mach numbers, there is no one method to achieve variable geometry. The F15 Eagle employs wedge inlets with adjustable flaps to control the flow. For subsonic flow, the flaps are completely closed and for supersonic flow, the flaps are open. The Concorde employed an externalcompression inlet, using a series of oblique shocks followed by a normal shock to slow the flow sufficiently for the turbojet engine. Perhaps the most recognizable supersonic aircraft, the SR71, used a hydraulically actuated cone to reduce the speed of the supersonic flow through the aircraft inlet.
Natural gas pipeline
Constructing natural gas line in winter, Finland
Natural gas pipelines are used to transport natural gas from extraction sites to refinement or chemical processing facilities. In the United States there are more than 210 natural gas pipeline systems with more than 305,000 miles (491,000 km) of intrastate transmission pipelines.^{[5]} Two compressible flow phenomenon characterize the flow through these pipelines: friction (Fanno flow) and (Rayleigh flow) and heat transfer. Natural gas pipelines are buried in the ground at a constant temperature of 15 °C. However, the friction generated by the flow offsets the heat loss to the Earth, thus resulting in an isothermal flow.
Friction factor as a function of Mach number in Fanno flow
The relationship between fL_{max}/D and Mach number for Fanno flow suggests that only subsonic flow can be used in the long pipes used to transport natural gas (even these pipes must be broken into shorter segments with compressor stations at the discontinuities in the pipeline). Additionally using conservation, an equation can be derived to describe the flow.

\frac=\frac \frac
This equation describes flow that chokes at M^{*} = 0.87 for natural gas γ = 1.32; however choking requires an infinite heat flux. Therefore, a combination of intuition and mathematics explains why it is most economically feasible that subsonic natural gas is pumped through long sections of pipe to reach its intended destination.
See also
References

^ But see compressibility which lists compressibilities for water and some other liquids

^ ^{a} ^{b} [4]

^ [5]

^ P. M. Schuster:Moving the Stars: Christian Doppler  His Life, His Works and Principle and the World After, Pollauberg, Austria:Living Edition Publishers, 2005

^ http://www.eia.gov/pub/oil_gas/natural_gas/analysis_publications/ngpipeline/index.html

Liepmann, Hans W.; Roshko, A. (1957) [1957]. Elements of Gasdynamics.

Anderson, John D. Jr. (2003) [1982]. Modern Compressible Flow (3rd ed.).

John, James E.; Keith, T. G. (2006) [1969]. Gas Dynamics (3rd ed.).

Oosthuizen, Patrick H.; Carscallen, W. E. (2013) [1997]. Introduction to Compressible Flow (2nd ed.).

Zucker, Robert D.; Biblarz, O. (2002) [1977]. Fundamentals of Gas Dynamics (2nd ed.).

Shapiro, Ascher H. (1953). The Dynamics and Thermodynamics of Compressible Fluid Flow, Volume 1.

Anderson, John D. Jr. (2000) [1989]. Hypersonic and High Temperature Gas Dynamics.
External links

NASA Beginner's Guide to Compressible Aerodynamics

Virginia Tech Compressible Flow Calculators

[6]
This article was sourced from Creative Commons AttributionShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and USA.gov, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for USA.gov and content contributors is made possible from the U.S. Congress, EGovernment Act of 2002.
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a nonprofit organization.