What is the Prandtl-Glauert correction?

The Prandtl-Glauert rule is a simple compressibility correction that converts incompressible pressure coefficient C_p and lift coefficient C_L into their subsonic compressible values using the factor 1/beta = 1/sqrt(1-M^2). For example at M = 0.6 the factor is 1.25, so C_p and C_L are roughly 25 percent larger than the incompressible values. The rule is accurate for thin airfoils up to about M = 0.7 and is widely used in preliminary design and wind-tunnel data correction.

Why does the pressure coefficient grow with Mach number?

In subsonic compressible flow the local density and temperature change around an airfoil, so the same shape produces larger pressure differences than in incompressible flow. Linearized analysis shows that the correction takes the simple form 1/sqrt(1-M^2). Physically, as the freestream Mach number approaches one, disturbances propagate more slowly so the suction peaks deepen and the lower-surface pressure rises.

How far is the Prandtl-Glauert rule valid?

The correction is based on linearized subsonic potential flow, so accuracy degrades quickly above M = 0.7 where local supersonic patches and shocks appear. As a rule of thumb the rule is used for M up to about 0.6 to 0.7. Beyond that the Karman-Tsien or Laitone corrections give better agreement, and full transonic CFD using Euler or Navier-Stokes solvers is required for design-quality answers.

How is the stagnation-to-static pressure ratio computed?

For an isentropic perfect gas the stagnation-to-static pressure ratio is P_0/P = (1 + (gamma-1)/2 * M^2)^(gamma/(gamma-1)). At gamma = 1.4 and M = 0.6 this gives P_0/P about 1.276. This formula is used to extract the local Mach number from pitot-static measurements and is independent of the Prandtl-Glauert factor: 1/beta describes the linearized surface pressure correction, while P_0/P describes the freestream state.

Prandtl-Glauert Simulator — Subsonic Compressibility Correction

Parameters

Incompressible C_p

—

Mach number M

—

Incompressible C_L

—

Specific heat ratio γ

—

Assumes thin-airfoil, small-disturbance, subsonic potential flow. Accurate up to roughly M = 0.7; the correction factor diverges as M approaches one.

Results

—

Compressible C_p

—

Compressible C_L

—

Factor 1/β

—

Stagnation/static P_0/P

Subsonic Airfoil Pressure Field

Upper surface shown in blue (suction), lower surface in red (pressure). Higher Mach number widens the pressure difference and raises the lift.

Correction Factor 1/β vs Mach

Horizontal axis = Mach number M, vertical axis = factor 1/β. The yellow marker shows the current M; 1/β diverges as M → 1.

Theory & Key Formulas

The Prandtl-Glauert correction comes from linearized small-disturbance subsonic potential flow around a thin airfoil and rescales the incompressible coefficients:

$$C_{p,\text{compr}} = \frac{C_{p,\text{inc}}}{\sqrt{1 - M^2}}, \qquad C_{L,\text{compr}} = \frac{C_{L,\text{inc}}}{\sqrt{1 - M^2}}$$

The correction factor 1/β is defined as:

$$\frac{1}{\beta} = \frac{1}{\sqrt{1 - M^2}}$$

For an isentropic perfect gas the stagnation-to-static pressure ratio is:

$$\frac{P_0}{P} = \left(1 + \frac{\gamma - 1}{2}\,M^2\right)^{\gamma/(\gamma - 1)}$$

As M → 1 the factor 1/β diverges and the rule breaks down in the transonic regime. The practical range is M ≤ 0.7.

What is the Prandtl-Glauert Compressibility Correction Simulator?

🙋

Does the pressure field around an aircraft change with Mach number, or is incompressible flow still a good model below the speed of sound?

🎓

Incompressible flow is fine up to about Mach 0.3, but between 0.5 and 0.7 the density changes start to matter and you need a compressibility correction. The simplest one is Prandtl-Glauert. Set C_p,inc = −0.50, M = 0.60 and γ = 1.40 in the simulator and you should see the factor 1/β = 1.250 and C_p,compr = −0.625 — about 25 percent larger.

🙋

Twenty-five percent is a lot. So lift coefficient grows as well?

🎓

Yes. The same factor multiplies C_L, so with C_L,inc = 0.50 the compressible value becomes 0.625. At cruise the wing produces more lift than the low-speed wind-tunnel data would suggest, which directly changes trim, drag and pitching moment. That is why airliner pre-design always carries compressibility through the lift, drag and moment estimates.

🙋

The curve looks like it blows up near Mach 0.9. How far should I trust it?

🎓

Good catch. The factor 1/β diverges as M → 1, but the rule itself is only based on linearized small-disturbance theory. Beyond about Mach 0.7 local supersonic patches and shocks appear, and the linear theory breaks down well before the divergence. For practical work treat the left half of the curve, M ≤ 0.7, as reliable; above that switch to Karman-Tsien, Laitone or a full transonic CFD solver.

🙋

What about the stagnation-to-static pressure ratio? Is that the same correction?

🎓

No — it is a separate isentropic relation that tells you the freestream state, not the surface correction. P_0/P = (1 + (γ−1)/2 · M²)^(γ/(γ−1)) gives 1.276 at γ = 1.4, M = 0.6. It is the basis for pitot-static airspeed indication. The simulator shows both 1/β and P_0/P so you can see two different Mach-number effects side by side without confusing them.

Frequently Asked Questions

Prandtl-Glauert is the simplest linearized correction, C_p,compr = C_p,inc / sqrt(1−M²). Karman-Tsien folds part of the isentropic gas relations back into the correction to extend usefulness toward transonic Mach numbers, while Laitone adds a local-Mach refinement and works best for slightly thicker airfoils. As rules of thumb, Prandtl-Glauert is reliable up to M = 0.6, Karman-Tsien up to 0.85, and beyond that full-potential or Euler/Navier-Stokes solvers are needed. The simulator covers the most basic rule, which is the right starting point for understanding the physical limits.

In linearized subsonic aerodynamics every linearly-scaled aerodynamic coefficient is multiplied by 1/β, so the pitching-moment coefficient C_m grows the same way as C_p and C_L. Induced drag C_Di scales with lift squared, so it goes roughly as 1/β². Viscous drag (skin friction and pressure drag from separation) does not follow the linear theory and must be computed separately from boundary-layer methods. Wave drag, which appears above M = 0.7, is outside the Prandtl-Glauert framework entirely and requires shock-capturing solvers.

The Prandtl-Glauert correction is based on thin-airfoil, small-disturbance and small angle-of-attack assumptions. In practice it gives within about 5 percent error for thickness-to-chord ratios up to 12 percent, angles of attack under five degrees and camber up to about 4 percent, all at M ≤ 0.6. Beyond those bounds non-linear effects matter: thick airfoils or high alpha drive local Mach to one and trigger shocks, and the linearized theory breaks. The rule is excellent for preliminary screening but final aerodynamic data must come from wind-tunnel tests or CFD.

Yes, it is the correct physical limit: in the incompressible limit there is no correction. So for hang gliders, UAVs at low speeds and cars (Mach below 0.3) incompressible Navier-Stokes or incompressible potential-flow methods are enough. By contrast at airliner cruise Mach numbers, 0.78 to 0.85, the correction factor reaches 1.6 to 1.9; ignoring it would underestimate design lift by 40 to 90 percent. The simulator lets you sweep M to see how quickly the correction grows.

Physical Model & Key Equations

The Prandtl-Glauert rule is derived from the linearized small-disturbance subsonic potential equation $(1-M_\infty^2)\phi_{xx} + \phi_{yy} = 0$. The Prandtl coordinate stretching $x = X$, $y = Y\sqrt{1-M_\infty^2}$ converts the compressible equation into the incompressible Laplace equation $\phi_{XX} + \phi_{YY} = 0$. Together with linearized boundary conditions and the surface pressure relation, this gives the correction:

$$C_{p,\text{compr}}(x, y) = \frac{C_{p,\text{inc}}(x, y/\beta)}{\sqrt{1 - M_\infty^2}}$$

where $\beta = \sqrt{1 - M_\infty^2}$ is the Prandtl-Glauert factor. The lift and pitching-moment coefficients scale with the same factor:

$$C_L^{\text{compr}} = \frac{C_L^{\text{inc}}}{\sqrt{1 - M_\infty^2}}, \qquad C_m^{\text{compr}} = \frac{C_m^{\text{inc}}}{\sqrt{1 - M_\infty^2}}$$

From isentropic perfect-gas relations the stagnation-to-static pressure ratio is:

$$\frac{P_0}{P} = \left(1 + \frac{\gamma - 1}{2}M^2\right)^{\gamma/(\gamma - 1)}$$

This is a separate identity used in pitot-static instrumentation. With $T_0/T = 1 + (\gamma-1)/2 \cdot M^2$ it provides the freestream state for any subsonic compressible flow. The simulator displays both 1/β and P_0/P so the surface correction and the freestream state can be studied at the same time.

Real-World Applications

Subsonic jet pre-design: In the early sizing of cruise-Mach 0.7 to 0.85 airliners and business jets, incompressible airfoil databases (NACA and supercritical SC-series) are scaled with the Prandtl-Glauert rule to estimate compressible C_L and C_p. Final aerodynamics are confirmed by CFD and wind-tunnel tests, but the simplicity of Prandtl-Glauert is still valuable for trim, surface area and stability iterations. At a typical cruise of M = 0.78 the factor is about 1.6, so ignoring it produces large errors in initial lift sizing.

Compressibility correction of wind-tunnel data: Data taken in low-speed tunnels at M ≤ 0.3 are routinely scaled with the Prandtl-Glauert rule to the cruise Mach number, and vice versa. The procedure is standard practice in JAXA, AIAA and NASA test reports and lets the same airfoil data set serve multiple Mach numbers without having to repeat costly transonic tunnel runs.

Propeller-tip Mach corrections: Large propeller tips can run at local Mach 0.7 to 0.9 while the inner blade sees much lower Mach. Applying Prandtl-Glauert section-by-section gives a quick lift distribution, blade efficiency and noise estimate. The same approach is used for helicopter rotor blades and UAV propellers in preliminary design.

High-altitude UAV operation: Stratospheric platforms cruise at low air density and high true airspeed, sometimes reaching Mach 0.5 to 0.7. Combining the local static-pressure variation with altitude and Prandtl-Glauert gives a fast initial estimate of the pressure load on the wing and helps drive structural sizing before more detailed CFD is run.

Common Misconceptions & Caveats

The most common mistake is to extend Prandtl-Glauert into the transonic regime. The factor 1/β formally produces a number for any M < 1, but its physical validity ends near M = 0.7 where local supersonic regions and shocks appear; the linear theory cannot describe wave drag or shock-induced separation. In the simulator setting M = 0.9 gives 1/β ≈ 2.29, but in reality the surface flow is dominated by shock dynamics and the linear C_p,compr value is meaningless.

A second confusion is between the correction factor 1/β and the stagnation-to-static pressure ratio P_0/P. Both depend on M but they mean different things. 1/β scales the local surface pressure coefficient (a linearized theory result), while P_0/P is a thermodynamic identity describing the freestream state of the gas. At the default values M = 0.6 and γ = 1.4 the two happen to be close (1.250 vs 1.276) but they diverge as Mach grows: 1/β blows up as M → 1, while P_0/P only rises from 1.28 to 1.89 between M = 0.6 and 1.0.

Finally there is the misconception that Prandtl-Glauert is only used to scale wind-tunnel data. In fact it is also used to apply compressibility to incompressible CFD output. Tools like XFOIL and several panel-method codes are formulated as incompressible solvers and combine Karman-Tsien (or Prandtl-Glauert as a fallback) to deliver compressible C_p distributions. Understanding how 1/β changes with Mach number, as the simulator shows, develops the intuition needed to use these mixed methods correctly.