Unit 5: Brief introduction to the methods of characteristics

I. Foundations of Compressible-Flow Analysis

Compressible-flow methods describe how density changes influence aerodynamic disturbances. The governing classification depends primarily on the free-stream Mach number (M\infty=U\infty/a\infty), where (U\infty) is free-stream velocity and (a_\infty) is free-stream sound speed.

• Flow regimes:
  • Subsonic: (M_\infty<1); disturbances can propagate upstream and downstream.
  • Transonic: (M_\infty\approx1); subsonic and supersonic regions coexist, often with shocks.
  • Supersonic: (M_\infty>1); disturbances are confined by Mach waves.
  • Hypersonic: conventionally (M_\infty\gtrsim5); strong shocks and high-temperature effects may arise.
• Basic assumptions: Introductory theories commonly assume steady, two-dimensional, inviscid, adiabatic, irrotational flow of a perfect gas.
• Velocity potential: For irrotational flow, velocity is written as (\mathbf{V}=\nabla\Phi), where (\Phi) is the total velocity potential.
• Small-disturbance representation:

Φ = U∞x + φ,     |φx|, |φy| ≪ U∞

Here (x,y) are Cartesian coordinates, (\phi) is the perturbation potential, and subscripts denote partial differentiation.

• Equation type: The potential equation is elliptic in subsonic flow, approximately parabolic near sonic conditions, and hyperbolic in supersonic flow.

II. Method of Characteristics — Supersonic-Flow Construction

A. Method of characteristics

The method of characteristics converts a hyperbolic partial differential equation into compatibility relations integrated along special curves called characteristics.

• Purpose: It determines velocity, pressure, Mach number, and flow direction in supersonic regions such as nozzles, expansion corners, and shock-free compression systems.
• Mach angle:

μ = sin⁻¹(1/M)

Here (\mu) is the Mach angle and (M>1) is the local Mach number. Mach lines make angles (\theta+\mu) and (\theta-\mu) with the reference axis, where (\theta) is the local flow angle.

• Characteristic directions:

dy/dx = tan(θ + μ),     dy/dx = tan(θ − μ)

Information propagates along these two characteristic families rather than freely in every direction.

• Prandtl–Meyer function:

ν(M) = √[(γ+1)/(γ−1)] tan⁻¹√{[(γ−1)/(γ+1)](M²−1)}
       − tan⁻¹√(M²−1)

Here (\nu) is the Prandtl–Meyer angle and (\gamma) is the ratio of specific heats.

• Compatibility relations:

  1. Along a characteristic of slope (\tan(\theta+\mu)), (\theta-\nu) remains constant.
  2. Along a characteristic of slope (\tan(\theta-\mu)), (\theta+\nu) remains constant.

  Thus, at the intersection of characteristics carrying constants (K-) and (K+),

θ = (K+ + K−)/2,     ν = (K+ − K−)/2

• Expansion fan: At a convex corner, infinitely many Mach waves form a centered Prandtl–Meyer fan; the flow turns smoothly, (M) rises, and pressure falls.
• Applications and limitations: The method is effective for two-dimensional or axisymmetric supersonic design, but shocks require separate Rankine–Hugoniot relations, and viscous boundary layers are not represented.

III. Compressibility Transformations — Relating Compressible and Incompressible Flow

A. Prandtl-Glauert and Goethert rules

The Prandtl–Glauert and Goethert rules transform linear subsonic compressible flow into an equivalent incompressible-flow problem.

• Linear equation:

β²φxx + φyy = 0,     β = √(1 − M∞²)

Here (\beta) is the subsonic compressibility factor and (M_\infty<1).

• Prandtl–Glauert rule: For the same thin airfoil and angle of attack,

Cp = Cp,0/β,     Cl = Cl,0/β,     Cm = Cm,0/β

Here (C_p), (C_l), and (C_m) are compressible pressure, lift, and pitching-moment coefficients; subscript (0) denotes corresponding incompressible values.

• Physical meaning: As (M_\infty) increases, a given geometric disturbance produces larger pressure and lift coefficients according to linear theory.
• Goethert rule: Through the coordinate transformation (Y=\beta y), compressible flow around a given thin profile can be related to incompressible flow around a geometrically transformed profile. An equivalent incompressible ordinate may be expressed as

y0 = y/β

where (y) is the compressible-flow profile ordinate and (y_0) is the equivalent incompressible-profile ordinate.

• Contrast:
  1. Prandtl–Glauert: retains the airfoil geometry and corrects aerodynamic coefficients.
  2. Goethert: interprets compressibility through a geometry or coordinate transformation.
• Limitations: Both rely on small disturbances and attached flow. Their singular prediction as (M_\infty\to1) is nonphysical because nonlinear transonic effects and shocks then become important.

IV. Linearized Supersonic Airfoil Analysis

A. Ackeret’s supersonic airfoil theory

Ackeret’s theory relates the pressure on a thin supersonic airfoil directly to its small local turning angle.

• Pressure relation:

Cp = 2θ/√(M∞² − 1)

Here (Cp=(p-p\infty)/q\infty), (p) is local pressure, (p\infty) is free-stream pressure, (q\infty=\tfrac12\rho\infty U_\infty^2) is dynamic pressure, and (\theta) is the signed flow-deflection angle in radians.

• Compression and expansion:
  1. A surface turning the flow toward itself produces (\theta>0), compression, and (C_p>0).
  2. A surface turning the flow away produces (\theta<0), expansion, and (C_p<0).
• Flat-plate airfoil: For a thin flat plate at small angle of attack (\alpha),

Cl = 4α/√(M∞² − 1)
Cd,w = 4α²/√(M∞² − 1)

Here (Cl) is lift coefficient and (C{d,w}) is wave-drag coefficient.

• Wave drag: Unlike ideal incompressible potential flow, supersonic flow produces drag through pressure forces associated with compression and expansion waves.
• Validity: The theory requires (M\infty>1), small surface slopes, small (\alpha), thin profiles, and attached weak waves. It becomes unreliable near (M\infty=1) and across strong shocks.

V. Perturbation Theory Across Mach-Number Regimes

A. Small perturbation equations for subsonic, transonic, supersonic and hypersonic flow

Small-perturbation theory simplifies compressible potential flow by expanding about a uniform free stream and retaining terms appropriate to each Mach-number regime.

• Common linear equation:

(1 − M∞²)φxx + φyy = 0

The sign of (1-M_\infty^2) determines the mathematical and physical character of the solution.

1. Subsonic flow

   • Equation type: For (M_\infty<1), both second-derivative coefficients have the same sign, so the equation is elliptic.
   • Consequence: Boundary conditions over the complete body influence the entire flow field; disturbances can travel upstream.
   • Solution basis: The transformation (Y=\sqrt{1-M_\infty^2}\,y) reduces the equation to Laplace’s equation.

2. Transonic flow

   • Nonlinearity: Near (M\infty=1), the coefficient of (\phi{xx}) becomes small, so neglected nonlinear terms become comparable to linear terms.
   • Transonic small-disturbance equation:

[1 − M∞² − (γ+1)M∞²φx/U∞]φxx + φyy = 0

• Mixed type: The bracketed coefficient may change sign within the field, producing adjacent subsonic and supersonic regions.
• Physical result: Local supersonic pockets may terminate in shocks even when (M_\infty<1).

3. Supersonic flow
   • Hyperbolic form:

(M∞² − 1)φxx − φyy = 0

• Consequence: Solutions depend on data propagated along two characteristic families.
• Domain of dependence: A point is influenced only by disturbances inside its upstream Mach cone, explaining why downstream geometry cannot affect an upstream supersonic region.

4. Hypersonic flow
   • Scaling: Although the formal linear equation remains hyperbolic, linearization is not uniformly valid when (M_\infty\theta=O(1)), where (\theta) is body slope.
   • Hypersonic similarity: Bodies having comparable values of the similarity parameter

K = M∞τ

 can exhibit related pressure fields; \(\tau\) is a representative thickness ratio or surface slope.

• Dominant effects: Strong shocks, entropy variation, viscous interaction, and aerodynamic heating increasingly limit potential-flow descriptions.

VI. Measured Airfoil Behaviour

A. Experimental characteristics of airfoils in incompressible flow

Wind-tunnel experiments establish real airfoil performance by measuring forces, moments, and surface pressures at low Mach number.

• Incompressible condition: Density variation is negligible, commonly approximated when (M_\infty\lesssim0.3).
• Measured coefficients:

Cl = L/(q∞S),     Cd = D/(q∞S),     Cm = M/(q∞Sc)

Here (L), (D), and (M) are lift, drag, and pitching moment; (S) is reference area and (c) is chord.

• Lift curve: Before separation, a two-dimensional thin airfoil approximately follows

Cl = 2π(α − αL=0)

Angles are in radians, and (\alpha{L=0}) is the zero-lift angle. A symmetric airfoil has (\alpha{L=0}\approx0), while a cambered airfoil generally has a negative zero-lift angle.

• Stall: Beyond a critical angle, extensive upper-surface separation causes (Cl) to reach (C{l,\max}) and then decrease; drag rises sharply.
• Drag behaviour: Profile drag combines skin-friction and pressure drag. A useful experimental approximation is (Cd=C{d0}+kCl^2), where (C{d0}) is minimum-profile drag and (k) is an empirical constant.
• Pitching moment: For conventional subsonic airfoils, the aerodynamic center lies near the quarter-chord, and (C_{m,c/4}) is approximately constant over the attached-flow range.
• Pressure distribution: Pressure taps determine (Cp=(p-p\infty)/q_\infty); integrating (C_p) over the surface gives pressure contributions to lift, drag, and moment.
• Reynolds-number effect:

Re = ρ∞U∞c/μ∞

Here (\rho\infty) is density and (\mu\infty) is dynamic viscosity. Reynolds number controls transition, boundary-layer thickness, separation, drag, and maximum lift.

• Experimental limitations: Wall interference, support forces, turbulence level, surface roughness, and finite-span effects must be corrected before measurements represent a two-dimensional airfoil accurately.