Compressible Flow Equations

Why This Matters

Compressible flow equations form the mathematical backbone of high-speed aerodynamics, propulsion systems, and gas dynamics. You're being tested on your ability to apply conservation principles when density can no longer be treated as constant—a fundamental shift from the incompressible assumptions you've seen before. These equations govern everything from shock wave behavior to nozzle design to duct flow with friction and heat transfer, and understanding when and how to apply each one separates surface-level memorization from genuine problem-solving ability.

The key insight is that compressible flow introduces coupling between thermodynamic and fluid dynamic variables. Pressure, density, temperature, and velocity all become interdependent through the equation of state and energy conservation. Master the underlying physics—isentropic processes, shock discontinuities, friction effects, and heat addition mechanisms—and you'll be equipped to tackle any problem the exam throws at you. Don't just memorize formulas; know which physical scenario each equation set describes and what assumptions it requires.

---

Conservation Foundations

Every compressible flow analysis begins with the fundamental conservation laws. These equations express mass, momentum, and energy balance in forms that account for variable density—the defining characteristic of compressible flow.

Continuity Equation

• Conservation of mass—expressed as $\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{u}) = 0$, or for quasi-1D flow: $\frac{\partial (\rho A)}{\partial t} + \frac{\partial (\rho u A)}{\partial x} = 0$
• Density variations are explicitly tracked, unlike incompressible flow where $\nabla \cdot \mathbf{u} = 0$ suffices
• Steady-state simplification yields $\rho_1 u_1 A_1 = \rho_2 u_2 A_2$, essential for nozzle and duct analysis

Momentum Equation

• Newton's second law for fluids—balances inertial forces against pressure gradients, viscous stresses, and body forces
• Navier-Stokes form for compressible flow: $\rho \frac{D\mathbf{u}}{Dt} = -\nabla p + \nabla \cdot \boldsymbol{\tau} + \rho \mathbf{g}$, where $\boldsymbol{\tau}$ is the viscous stress tensor
• Euler equations result when viscosity is neglected, appropriate for high Reynolds number flows away from boundaries

Energy Equation

• First law of thermodynamics—accounts for internal energy, kinetic energy, pressure work, heat transfer, and viscous dissipation
• Total enthalpy $h_0 = h + \frac{u^2}{2}$ remains constant along streamlines for adiabatic, inviscid flow
• Stagnation temperature $T_0 = T\left(1 + \frac{\gamma - 1}{2}M^2\right)$ connects static and total conditions through Mach number

---

Thermodynamic Closure

The conservation equations alone don't close the system—you need relationships connecting pressure, density, and temperature. These constitutive relations complete the mathematical framework.

Equation of State

• Ideal gas law $p = \rho R T$ relates thermodynamic variables, where $R$ is the specific gas constant
• Closes the equation system—without it, you have more unknowns than equations
• Specific heat ratio $\gamma = c_p/c_v$ characterizes the gas and appears throughout compressible flow relations

Speed of Sound Equation

• Wave propagation speed defined as $a = \sqrt{\gamma R T}$ for an ideal gas, derived from isentropic compressibility
• Mach number $M = u/a$ is the fundamental dimensionless parameter distinguishing flow regimes
• Compressibility criterion—flows with $M > 0.3$ typically require compressible analysis due to density variations exceeding ~5%

---

Isentropic Flow Analysis

When flow is both adiabatic (no heat transfer) and reversible (no shocks or friction), entropy remains constant. These idealized conditions yield elegant closed-form relationships widely used in preliminary design.

Isentropic Flow Relations

• Pressure-density-temperature coupling: $\frac{p}{p_0} = \left(\frac{\rho}{\rho_0}\right)^\gamma = \left(\frac{T}{T_0}\right)^{\gamma/(\gamma-1)}$
• Mach number dependence: $\frac{T}{T_0} = \left(1 + \frac{\gamma-1}{2}M^2\right)^{-1}$ relates local conditions to stagnation values
• Nozzle design foundation—area-Mach relation $\frac{A}{A^*} = \frac{1}{M}\left[\frac{2}{\gamma+1}\left(1 + \frac{\gamma-1}{2}M^2\right)\right]^{(\gamma+1)/2(\gamma-1)}$ governs converging-diverging nozzle geometry

Prandtl-Meyer Expansion Equations

• Centered expansion waves—flow accelerates isentropically around convex corners, unlike shocks which are irreversible
• Prandtl-Meyer function $\nu(M) = \sqrt{\frac{\gamma+1}{\gamma-1}}\tan^{-1}\sqrt{\frac{\gamma-1}{\gamma+1}(M^2-1)} - \tan^{-1}\sqrt{M^2-1}$
• Turning angle relation: $\theta = \nu(M_2) - \nu(M_1)$ connects deflection angle to Mach number change for supersonic flow

---

Shock Wave Discontinuities

When supersonic flow decelerates abruptly, conservation laws must still hold across the discontinuity—but entropy increases irreversibly. Shock relations connect upstream and downstream states through algebraic jumps.

Normal Shock Equations

• Rankine-Hugoniot relations express conservation across the shock: $\frac{p_2}{p_1} = \frac{2\gamma M_1^2 - (\gamma-1)}{\gamma+1}$ for pressure ratio
• Mach number jump: $M_2^2 = \frac{M_1^2 + \frac{2}{\gamma-1}}{\frac{2\gamma}{\gamma-1}M_1^2 - 1}$, always yielding $M_2 < 1$ downstream of a normal shock
• Entropy increase $\Delta s > 0$ makes shocks irreversible—stagnation pressure drops while stagnation temperature remains constant

---

Duct Flow with Irreversibilities

Real internal flows involve friction and heat transfer, which modify the isentropic picture. Fanno and Rayleigh flows isolate these effects to reveal their individual influences on flow properties.

Fanno Flow Equations

• Adiabatic flow with friction—stagnation temperature $T_0$ remains constant while stagnation pressure $p_0$ decreases
• Choking behavior: friction drives flow toward $M = 1$ regardless of initial conditions—subsonic flow accelerates, supersonic flow decelerates
• Fanno line on $h$-$s$ diagram shows the locus of states; maximum entropy occurs at sonic condition

Rayleigh Flow Equations

• Frictionless flow with heat transfer—stagnation temperature changes while friction is neglected
• Heating effects: subsonic flow accelerates toward $M = 1$; supersonic flow decelerates toward $M = 1$—both choke with sufficient heat addition
• Maximum entropy at $M = 1$, but maximum stagnation temperature occurs at $M = 1/\sqrt{\gamma}$ for heating processes

---

Quick Reference Table

---

Self-Check Questions

1. Which two equation sets both predict that flow will approach $M = 1$, and what physical mechanism drives this choking in each case?

2. Compare and contrast normal shock waves and Prandtl-Meyer expansions: which conserves stagnation pressure, and why does the other not?

3. If you're analyzing a converging-diverging nozzle operating at design conditions with no shocks, which equation set provides the area-Mach relationship, and what assumptions does it require?

4. The energy equation introduces stagnation temperature $T_0$. Under what two duct flow conditions does $T_0$ remain constant, and under what condition does it change?

5. An FRQ describes supersonic flow entering a constant-area duct with significant wall friction. Which flow model applies, what happens to the Mach number, and what happens to the stagnation pressure?