Compressible Flow and Gas Dynamics

Introduction and Context

Compressible flow is a flow in which density changes are sufficiently important to influence the pressure, velocity, temperature, or mass-flow distribution. In liquids, density variations are often negligible; in gases, they become central when the flow speed approaches the speed of sound, when pressure ratios are large, or when heat transfer and chemical reactions are significant. Gas dynamics is the study of such flows, especially the interaction among pressure, density, temperature, velocity, and area.

Compressible-flow analysis is fundamental to turbines, compressors, nozzles, diffusers, rockets, aircraft, internal-combustion engines, supersonic inlets, high-speed wind tunnels, and safety systems involving gas release. The central conceptual transition is from subsonic behaviour, in which pressure disturbances can propagate upstream, to supersonic behaviour, in which they cannot. This transition produces choking, shocks, expansion fans, and strong geometric constraints.

Learning Objectives

By the end of this chapter, the reader should be able to:

1. Explain compressibility, Mach number, stagnation properties, and acoustic wave propagation.
2. Apply the conservation equations to one-dimensional compressible flow.
3. Derive and use the isentropic area–velocity and mass-flow relations.
4. Analyse choked flow through nozzles and restrictions.
5. Explain normal shocks, oblique shocks, and Prandtl–Meyer expansion waves.
6. Distinguish Fanno flow, Rayleigh flow, and isentropic flow.
7. Evaluate the assumptions and limitations of ideal-gas gas-dynamic models.
8. Interpret compressible-flow results in engineering design and safety applications.

1. Fundamental Concepts and Terminology

1.1 Compressibility and acoustic speed

The speed of sound is the speed at which an infinitesimal pressure disturbance travels through a medium. For a calorically perfect gas undergoing an isentropic disturbance,

where a is the speed of sound (m/s), γ = c_p/c_v is the ratio of specific heats, R is the specific gas constant (J/(kg·K)), and T is absolute temperature (K). For air near ordinary conditions, γ ≈ 1.4 and R ≈ 287 J/(kg·K).

The Mach number is

where V is local flow speed. The usual classification is:

The threshold M = 0.3 is a practical approximation, not a fundamental discontinuity. Density variation depends on the required accuracy, pressure field, heat transfer, and gas properties.

1.2 Thermodynamic model

For an ideal gas,

where p is pressure (Pa), ρ is density (kg/m³), and T is temperature (K). For a calorically perfect gas, c_p, c_v, γ, and R are treated as constants, with

At high temperatures, vibrational excitation, dissociation, ionisation, and variable specific heats invalidate the calorically perfect approximation. Real-gas equations of state may then be necessary.

1.3 Stagnation properties

The stagnation, or total, state is the state obtained when a moving fluid is brought to rest adiabatically and reversibly. For steady adiabatic flow with no shaft work and negligible potential-energy change, stagnation enthalpy is conserved:

For a calorically perfect gas,

Thus, for an isentropic process,

The subscript 0 denotes stagnation conditions. These equations require a perfect gas and reversible adiabatic deceleration. Across a shock, stagnation temperature remains constant for an adiabatic flow without work, but stagnation pressure decreases because entropy increases.

2. Governing Equations for One-Dimensional Flow

Consider steady flow through a duct of area A, with mean velocity V. The one-dimensional continuity equation is

where ṁ is mass-flow rate (kg/s).

The steady-flow energy equation is

where q is heat added per unit mass, w_s is shaft work delivered by the fluid per unit mass, and z is elevation. For a nozzle, with negligible heat transfer, shaft work, and elevation change, h + V²/2 is constant.

For inviscid flow along a streamline, the momentum equation can be expressed as

For a general flow, entropy and transport effects must be included. The differential equation of state and the definition of the speed of sound give

where the derivative is evaluated at constant entropy. This relation is more general than a = √(γRT).

3. Isentropic Flow in Variable-Area Ducts

An isentropic flow is adiabatic and reversible. Combining continuity, momentum, and the thermodynamic relation dh = dp/ρ gives the important area–velocity relation:

This equation explains the fundamentally different geometry required for subsonic and supersonic acceleration. If M < 1, acceleration dV > 0 requires dA < 0; a converging duct accelerates subsonic flow. If M > 1, acceleration requires dA > 0; a diverging duct accelerates supersonic flow. At M = 1, a smooth extremum of area is required, normally a throat.

The area–Mach-number relation for a perfect gas is

Here A* is the area at which the flow would have M = 1 under the same stagnation conditions. For a given area ratio greater than one, two solutions may exist: one subsonic and one supersonic. The physically realised branch depends on back pressure and the complete flow configuration.

3.1 Nozzle mass-flow relation and choking

Using the isentropic relations in continuity gives

The mass-flow function reaches its maximum at M = 1. At the throat,

Choking occurs when the throat reaches sonic conditions. Further reduction of downstream pressure does not increase the mass flow, provided the upstream stagnation state and throat area remain fixed and the assumptions remain valid. The critical pressure ratio is

For air with γ = 1.4, this ratio is approximately 0.528. Choking is not synonymous with supersonic exit flow: a converging nozzle can be choked at its throat while its exit remains sonic, whereas a converging–diverging nozzle can develop a supersonic divergent section.

Worked Example 1: Choked air nozzle

Air enters a converging nozzle from a reservoir at p₀ = 500 kPa and T₀ = 300 K. The throat area is A* = 1.0 × 10⁻⁴ m². Take γ = 1.4 and R = 287 J/(kg·K). Find the maximum mass-flow rate.

Using the choked-flow equation,

The factors are approximately 0.02887, 0.06982, and 0.57870, respectively. Therefore,

The result is dimensionally consistent: pressure times area divided by the square root of temperature, multiplied by √(1/R), produces kg/s. The calculation assumes steady, one-dimensional, adiabatic, inviscid flow and a uniform stagnation state upstream.

4. Shock Waves

A shock wave is a thin region across which pressure, temperature, density, and entropy change rapidly. It is not isentropic. In an adiabatic shock without shaft work, total temperature remains constant, while total pressure decreases.

4.1 Normal shocks

For a stationary normal shock, upstream Mach number is M₁ > 1. For a perfect gas,

The downstream Mach number is

For every finite-strength normal shock, M₂ < 1. The entropy increase is positive, and the stagnation-pressure ratio is less than one. The stagnation-pressure ratio may be obtained from the isentropic-to-stagnation relations on each side:

A normal shock causes large total-pressure loss, so designs should avoid it when efficient compression is required.

4.2 Oblique shocks

When supersonic flow turns into itself, the compression wave coalesces into an oblique shock. The upstream Mach number normal to the shock is

where β is the shock angle. The normal-shock relations apply using M_n1. The shock angle, turning angle θ, and Mach number satisfy the θ–β–M relation:

Weak oblique shocks usually produce supersonic downstream flow and lower losses than a normal shock. Strong solutions may produce subsonic downstream flow. A detached bow shock forms when the required turning angle exceeds the maximum attached-shock turning angle or when the geometry is blunt.

4.3 Expansion waves

When supersonic flow turns away from itself, it undergoes a Prandtl–Meyer expansion. The turning is continuous through a centred fan and is approximately isentropic. The Prandtl–Meyer function is

The turning angle is θ = ν(M₂) − ν(M₁). Expansion increases Mach number and decreases pressure and temperature. Unlike a shock, an ideal expansion fan has no stagnation-pressure loss.

5. Duct Flows with Friction and Heat Transfer

5.1 Fanno flow

Fanno flow is steady, one-dimensional, constant-area flow of a perfect gas with adiabatic wall friction and no shaft work. Friction drives both subsonic and supersonic flow toward M = 1, but from opposite directions: subsonic flow accelerates, while supersonic flow decelerates.

For a constant-friction-factor model, the Fanno parameter is

where f is the Darcy friction factor, L* is the duct length required to reach sonic conditions, and D is hydraulic diameter. Tables or numerical integration are normally used for the Fanno function. Wall friction causes stagnation-pressure loss, even though stagnation temperature is constant.

5.2 Rayleigh flow

Rayleigh flow is steady, one-dimensional, constant-area flow with heat transfer and negligible friction. Heat addition drives subsonic flow toward sonic conditions by increasing velocity, whereas it drives supersonic flow toward sonic conditions by decreasing velocity. Heat addition increases entropy and generally reduces stagnation pressure. The model is useful for combustors and heated ducts, but it neglects wall friction and multidimensional mixing.

Fanno and Rayleigh models are limiting idealisations. Real devices often involve simultaneous friction, heat transfer, area change, chemical reaction, and turbulence; numerical compressible-flow solvers are then required.

6. Practical Design and Measurement Considerations

A compressible-flow design should begin by identifying the controlling state: reservoir or stagnation conditions, static conditions, back pressure, gas composition, and expected Mach-number range. The following practices are important:

• Use stagnation quantities for reservoirs and upstream supply systems, and static quantities for local flow measurements.
• Check choking before applying a pressure-driven mass-flow equation.
• Select the correct subsonic or supersonic branch of the area–Mach relation.
• Account for boundary layers, discharge coefficients, separation, and manufacturing tolerances.
• Treat shocks as sources of entropy and total-pressure loss.
• Verify whether variable specific heats or real-gas effects are important.
• Validate analytical models against experiments or computational-fluid-dynamics results.

Pressure measurements are especially delicate. A pressure tap aligned with the local flow may approximate static pressure, while a properly designed pitot probe measures stagnation pressure in subsonic flow. In supersonic flow, a pitot probe creates a normal shock ahead of the probe; the measured pressure is therefore related to the upstream Mach number through normal-shock theory rather than the subsonic Bernoulli equation.

7. Illustrative Case: Supersonic Nozzle with Excessive Back Pressure

Background

Consider an illustrative converging–diverging nozzle supplied from a large reservoir. The nozzle is designed to deliver a supersonic jet, but its back pressure is adjustable. The reservoir stagnation pressure and temperature remain nearly constant.

Problem

Determine qualitatively what happens as back pressure is reduced from a high value to a very low value.

Analysis

At high back pressure, the flow remains subsonic throughout the nozzle, and the divergent portion behaves as a subsonic diffuser: velocity decreases and static pressure rises. As back pressure is reduced to the critical value, the throat reaches M = 1, and mass flow becomes choked. Further reduction creates a supersonic region downstream of the throat, but a shock may stand in the divergent section if the back pressure is still too high for shock-free expansion.

As back pressure is reduced further, the shock moves downstream. Across the shock, Mach number and stagnation pressure decrease, producing a loss in thrust potential. At a sufficiently low back pressure, the shock is expelled from the nozzle and the flow expands supersonically to the exit. If the exit pressure differs from ambient pressure, external compression or expansion waves adjust the jet.

Approach and outcome

A robust engineering analysis uses the isentropic area–Mach relation in the smooth regions and normal-shock relations if an internal shock exists. The back pressure determines the operating regime; geometry alone does not determine the complete solution.

Learning points

1. Choking fixes the maximum mass flow for fixed upstream conditions and throat area.
2. A converging–diverging nozzle can have subsonic, transonic, shock-containing, or fully supersonic operation.
3. Internal shocks reduce stagnation pressure and therefore performance.
4. Nozzle design must include the operating back-pressure range, not only the design point.

8. Limitations, Safety, and Contemporary Practice

Ideal one-dimensional theory neglects boundary-layer displacement, three-dimensionality, turbulence, unsteady shock motion, heat losses, and chemical nonequilibrium. These effects can cause separation in overexpanded nozzles, side loads, vibration, acoustic emissions, and structural damage. In high-temperature propulsion, dissociation and ionisation alter γ, thermodynamic properties, and reaction rates.

Safety analysis must consider accidental gas release, blast waves, pressure-vessel rupture, toxic or flammable mixtures, and safe vent sizing. Standards and regulatory requirements should govern pressure-vessel design, relief systems, hazardous-area classification, and experimental facilities. Models should be conservative where uncertainty affects human safety.

Modern computational fluid dynamics solves the compressible Navier–Stokes equations using finite-volume, finite-difference, or finite-element methods. Shock-capturing schemes require numerical dissipation or limiters to avoid nonphysical oscillations. Turbulence models, adaptive meshing, uncertainty quantification, and high-speed experimental diagnostics remain active research areas. Machine learning can assist reduced-order modelling, surrogate construction, anomaly detection, and design optimisation, but it should not replace conservation laws, validation, or engineering judgement. Training data may inadequately represent shock interactions or rare operating conditions; therefore, extrapolative AI predictions require particular caution.

Chapter Summary

Compressible flow couples fluid mechanics and thermodynamics through density variation and acoustic propagation. Mach number classifies the influence of pressure disturbances. For ideal-gas, steady, adiabatic, reversible flow, stagnation properties and isentropic relations provide the foundation for nozzle and diffuser analysis. The area–velocity relation shows why subsonic acceleration requires convergence while supersonic acceleration requires divergence. Choking occurs when a minimum area reaches sonic conditions and establishes a maximum mass-flow rate.

Normal and oblique shocks are irreversible compression processes with entropy increase and stagnation-pressure loss. Prandtl–Meyer fans provide approximately isentropic supersonic expansion. Fanno and Rayleigh models isolate the effects of friction and heat transfer, respectively. Practical design requires attention to back pressure, boundary layers, separation, measurement methods, real-gas behaviour, safety, and computational validation.

Key Terms

Mach number: Ratio of flow speed to local speed of sound.

Speed of sound: Propagation speed of an infinitesimal pressure disturbance.

Compressible flow: Flow in which density variation materially affects the solution.

Ideal gas: Gas model satisfying p = ρRT.

Stagnation state: State obtained by reversible adiabatic deceleration to zero velocity.

Static property: Local thermodynamic property of the moving fluid.

Isentropic flow: Reversible adiabatic flow with constant entropy.

Choking: Condition in which a controlling section reaches M = 1 and mass flow is maximised.

Critical state: Sonic state associated with a specified stagnation condition.

Converging nozzle: Nozzle whose area decreases in the flow direction.

Diverging nozzle: Nozzle whose area increases in the flow direction.

Normal shock: Shock whose wave front is normal to the upstream flow.

Oblique shock: Shock inclined to the upstream flow direction.

Prandtl–Meyer expansion: Isentropic supersonic turning away from itself.

Fanno flow: Constant-area adiabatic flow with friction.

Rayleigh flow: Constant-area flow with heat transfer and negligible friction.

Stagnation pressure: Pressure obtained by isentropically bringing the flow to rest.

Back pressure: Pressure imposed downstream of a nozzle or duct.

Calorically perfect gas: Perfect gas with constant specific heats.

Review Questions

Short Answer Questions

1. Define Mach number and explain its physical significance.
2. What is meant by a stagnation state?
3. Why does a converging nozzle accelerate subsonic flow?
4. Define choking and state its effect on mass flow.
5. Why does a normal shock cause stagnation-pressure loss?

Long Answer Questions

1. Derive the area–velocity relation for steady, quasi-one-dimensional, isentropic flow.
2. Explain the operation of a converging–diverging nozzle as back pressure is reduced.
3. Derive and discuss the normal-shock pressure, density, temperature, and Mach-number relations.
4. Compare Fanno flow, Rayleigh flow, and isentropic flow in terms of assumptions and physical effects.
5. Discuss the limitations of ideal-gas gas dynamics in hypersonic and high-temperature applications.

Analytical / Application Questions

1. Air enters a choked nozzle with specified reservoir pressure, reservoir temperature, and throat area. Calculate the mass-flow rate and identify every modelling assumption.
2. A supersonic nozzle develops an internal normal shock. Explain how the shock location changes when back pressure is reduced, and identify the properties that are discontinuous or conserved.
3. Design a measurement and validation strategy for a high-speed duct experiment, including pressure instrumentation, uncertainty sources, shock detection, and CFD comparison.

Multiple Choice Questions with Answers

Question 1

For a calorically perfect gas, the local speed of sound is given by which expression?

A. a = V/M

B. a = √(γRT)

C. a = γRT

D. a = p/(ρV)

Correct Answer: B. a = √(γRT)

Question 2

In a steady, isentropic, variable-area duct, subsonic flow accelerates when the area:

A. Increases

B. Remains constant

C. Decreases

D. Changes discontinuously only

Correct Answer: C. Decreases

Question 3

For air with γ = 1.4, the approximate critical pressure ratio p*/p₀ is:

A. 0.528

B. 0.750

C. 1.000

D. 0.250

Correct Answer: A. 0.528

Question 4

At a finite-strength adiabatic normal shock, which quantity remains constant when shaft work and heat transfer are negligible?

A. Stagnation pressure

B. Entropy

C. Stagnation temperature

D. Static temperature

Correct Answer: C. Stagnation temperature

Question 5

A normal shock in a supersonic perfect-gas flow produces downstream flow that is:

A. Always supersonic

B. Always sonic

C. Subsonic for a finite-strength shock

D. At zero velocity

Correct Answer: C. Subsonic for a finite-strength shock

Question 6

Which process is approximately isentropic for an ideal Prandtl–Meyer expansion fan?

A. Supersonic turning away from itself

B. Supersonic turning into itself

C. Flow across a normal shock

D. Constant-area frictional flow

Correct Answer: A. Supersonic turning away from itself

Question 7

For choked flow through a nozzle, further reduction of back pressure generally:

A. Reduces the mass flow immediately

B. Does not increase mass flow if upstream stagnation state and throat area are fixed

C. Makes throat Mach number subsonic

D. Makes stagnation temperature zero

Correct Answer: B. Does not increase mass flow if upstream stagnation state and throat area are fixed

Question 8

Fanno flow is conventionally modelled as:

A. Constant-area adiabatic flow with friction

B. Isentropic flow with changing area

C. Constant-area frictionless flow with heat transfer

D. Flow through a normal shock only

Correct Answer: A. Constant-area adiabatic flow with friction

Question 9

In Rayleigh flow, heat addition drives both subsonic and supersonic states toward:

A. M = 0

B. M = 1

C. M = ∞

D. The same static pressure

Correct Answer: B. M = 1

Question 10

The principal performance penalty associated with a normal shock is a decrease in:

A. Stagnation pressure

B. Stagnation temperature in an adiabatic no-work flow

C. Mass conservation

D. Gas constant

Correct Answer: A. Stagnation pressure