Hypersonic Small-Disturbance Equations

• Hypersonic small-disturbance equations are reduced models that capture key high-Mach flow physics by retaining only leading-order terms from the full Euler or Navier–Stokes systems.
• They employ asymptotic similarity scalings for steady flows and linearization techniques for stability and receptivity analyses, revealing shock interactions and flow instabilities.
• Rigorous convergence rates and measure-theoretic formulations validate these models, enabling practical predictions in aerodynamic design and high-enthalpy flow analysis.

Hypersonic small-disturbance equations are reduced descriptions of high-Mach-number flow obtained by exploiting either hypersonic similarity scalings or infinitesimal perturbations about a computed base state. In the cited literature, the term appears in two related senses: the classical steady asymptotic equations for slender wedges, cones, and other bodies in the hypersonic limit, and the linearized evolution equations used for modal, non-modal, and receptivity analyses of high-enthalpy hypersonic flows (Kang et al., 23 Aug 2025, Antón-Álvarez et al., 4 Apr 2026). In both senses, the objective is to retain the leading flow physics while simplifying the full compressible governing equations enough to expose similarity structure, shock interaction, instability mechanisms, or receptivity pathways.

1. Terminology and conceptual scope

The literature uses “hypersonic small-disturbance equations” for two distinct but connected reductions. One is an asymptotic reduction of the steady compressible Euler or potential-flow equations under a hypersonic similarity limit, usually with a fixed similarity parameter such as $K:=M_\infty \tau$, $K=M_0\theta$, or $a_0=\tau M_0$. The other is a linearization of the compressible Navier–Stokes equations about a steady base flow $\bm{q}_0$, with disturbances $\bm{q}'$ treated as infinitesimal (Kang et al., 23 Aug 2025, Kuang et al., 2021, Antón-Álvarez et al., 4 Apr 2026).

Usage in the literature	Starting equations	Reduced form
Asymptotic steady formulation	Steady compressible Euler or steady potential flow	Limit system as $\tau \to 0$ with fixed $K$ or $a_0$
Linearized disturbance formulation	Compressible Navier–Stokes about $\bm{q}_0$	$$\partial \bm{q}'/\partial t = \bm{L}(\bm{q}_0)\bm{q}'$$

In the asymptotic setting, the equations formalize the Mach-number independence principle, also called the hypersonic similarity law or Van Dyke’s similarity theory. For a fixed similarity parameter, the flow solution after scaling depends on that parameter and on the gas constant $K=M_0\theta$0, rather than on Mach number and slenderness separately (Chen et al., 2023, Kang et al., 23 Aug 2025). In the linearized setting, the same small-disturbance idea is applied to base-flow stability and receptivity, especially when shocks, boundary layers, and real-gas effects must be retained in the operator (Antón-Álvarez et al., 4 Apr 2026).

A recurrent interpretive issue is that these two usages are not interchangeable. The asymptotic systems are nonlinear steady conservation laws; the stability systems are linearized evolution equations. This distinction is explicit in the cited work and is essential for reading the modern literature correctly (Kang et al., 23 Aug 2025, Antón-Álvarez et al., 4 Apr 2026).

2. Asymptotic derivation from Euler and potential-flow models

A standard starting point is the two-dimensional steady compressible Euler system for a polytropic gas,

$K=M_0\theta$1

with $K=M_0\theta$2 in one formulation (Qu et al., 2019). For steady non-isentropic Euler flow past slender bodies, one scaling used in the literature is

$K=M_0\theta$3

with fixed $K=M_0\theta$4. After dropping $K=M_0\theta$5 terms, the limiting two-dimensional hypersonic small-disturbance system becomes (Kang et al., 23 Aug 2025)

$K=M_0\theta$6

with

$K=M_0\theta$7

For an axisymmetric cone, the corresponding three-dimensional equations acquire radial weights: $K=M_0\theta$8 again with $K=M_0\theta$9 and $a_0=\tau M_0$0 on the cone surface (Kang et al., 23 Aug 2025).

Potential-flow specializations give simpler systems. In one scaled irrotational formulation, neglecting all terms involving $a_0=\tau M_0$1 yields

$a_0=\tau M_0$2

with wall condition $a_0=\tau M_0$3 (Kuang et al., 2021). In the isothermal hypersonic small-disturbance limit studied for large data, this further appears as

$a_0=\tau M_0$4

with $a_0=\tau M_0$5 on the wedge boundary (Chen et al., 2024).

The mathematical character of these asymptotic systems remains nontrivial. In the non-isentropic Euler setting, the derived small-disturbance system is non-strictly hyperbolic with two genuinely nonlinear and two linearly degenerate characteristic fields (Kang et al., 23 Aug 2025).

3. Boundary conditions, shocks, and singular limit structure

Hypersonic small-disturbance theory is inseparable from boundary conditions and shock relations. For steady wedge problems, the slip condition appears as $a_0=\tau M_0$6 on $a_0=\tau M_0$7 in one formulation (Qu et al., 2019), or as $a_0=\tau M_0$8 on the scaled boundary in similarity variables (Chen et al., 2023). In free-boundary potential-flow formulations over wedges, the shock location is unknown and Rankine–Hugoniot conditions determine it through relations such as

$a_0=\tau M_0$9

with $\bm{q}_0$0 across the shock (Hu, 5 Aug 2025).

A central feature of the hypersonic limit is collapse of the post-shock layer. For two-dimensional steady compressible Euler flow past a straight wedge, the Mach number of the upcoming uniform supersonic flow increases to infinite may be taken as the adiabatic exponent $\bm{q}_0$1 decreases to $\bm{q}_0$2, with fixed total energy $\bm{q}_0$3 (Qu et al., 2019). In that limit, the pressure on the surface obeys the Newtonian sine-square law,

$\bm{q}_0$4

the density grows like $\bm{q}_0$5 with $\bm{q}_0$6, and the post-shock region becomes vanishingly thin (Qu et al., 2019).

Because the support of the post-shock state collapses toward the body while density diverges, ordinary function-valued weak solutions are not adequate in general. Several papers therefore formulate the limit problem in terms of Radon measure solutions, with singular concentrations supported on the body boundary. One representative form is

$\bm{q}_0$7

for the mass flux measure on the wedge, with analogous expressions for momentum and energy (Qu et al., 2019). For slender-body similarity problems, the scaled density measure takes the form

$\bm{q}_0$8

again combining an absolutely continuous part with a Dirac part on the body surface (Kang et al., 23 Aug 2025).

This measure-theoretic viewpoint also recovers force laws that are classical in hypersonic aerodynamics. The straight-wedge limit justifies Newton’s theory through $\bm{q}_0$9 (Qu et al., 2019), while ramp and slender-body analyses recover Newton-Busemann-type laws for pressure force on the boundary (Jin et al., 2019, Kang et al., 23 Aug 2025). A common misconception is that hypersonic small-disturbance theory is purely regular. The cited mathematical literature shows that, in the hypersonic limit, the natural solution class may contain Dirac measures supported on walls or free layers (Jin et al., 2019).

4. Similarity law, entropy solutions, and convergence theory

The asymptotic equations are closely tied to the hypersonic similarity law, also called the Mach-number independence principle or Van Dyke’s similarity theory. In the slender-body setting, if the parameter $\bm{q}'$0 is fixed and $\bm{q}'$1, the flow field structures after scaling no longer depend on the body’s shape and the Mach number $\bm{q}'$2 independently, but only on $\bm{q}'$3 and adiabatic index $\bm{q}'$4 (Kang et al., 23 Aug 2025). Parallel statements are formulated for wedges with $\bm{q}'$5 or $\bm{q}'$6 fixed (Kuang et al., 2019, Kuang et al., 2021).

A substantial modern development is the rigorous validation of this similarity law in $\bm{q}'$7 for entropy solutions. For two-dimensional steady potential flow with large data, a modified Glimm scheme is used to construct global entropy solutions for fixed $\bm{q}'$8 and sufficiently large Mach number $\bm{q}'$9, and as $\tau \to 0$0, the scaled solutions converge in $\tau \to 0$1 to the solution of the corresponding initial-boundary value problem of the hypersonic small-disturbance equations (Kuang et al., 2019). This is described there as the first rigorous mathematical global result on the validation of the hypersonic similarity for the two dimensional steady potential flow.

Quantitative convergence rates are also available. For inviscid steady irrotational Euler flow over a two-dimensional Lipschitz slender wedge, the $\tau \to 0$2 difference between the solution to the scaled equations and the solution to the hypersonic small-disturbance equations is of order $\tau \to 0$3, provided the total variations of the initial data and the tangential derivative of the boundary are sufficiently small (Kuang et al., 2021). For compressible full Euler flows over two-dimensional slender Lipschitz wedges, the corresponding estimate is

$\tau \to 0$4

for all $\tau \to 0$5, with $\tau \to 0$6 independent of $\tau \to 0$7 and $\tau \to 0$8 (Chen et al., 2023). That work states that this is the first mathematical result on the comparison of two solutions of the compressible Euler equations with characteristic boundary conditions.

For steady potential flows with large data past a straight wedge, the convergence rate can be written as

$\tau \to 0$9

and the rate is shown to be optimal (Chen et al., 2024). The proofs use the Standard Riemann Semigroup, local $K$0 comparison of Riemann solutions, and front-tracking constructions (Chen et al., 2024). This suggests that the asymptotic reduction is not merely formal; it supports a sharp, global comparison theory between full and reduced models.

5. Linearized hypersonic small-disturbance equations in high-enthalpy stability theory

In high-enthalpy stability analysis, the small-disturbance formulation is obtained by linearizing the compressible Navier–Stokes equations about a steady base flow. The total flow field is decomposed as

$K$1

and retaining only terms linear in $K$2 yields

$K$3

where $K$4 is the Jacobian of the non-linear equation discretization evaluated at the base state $K$5 (Antón-Álvarez et al., 4 Apr 2026). For two-dimensional analysis,

$K$6

and with shock-fitting the state also includes shock displacement perturbations $K$7 (Antón-Álvarez et al., 4 Apr 2026).

This linearized framework is implemented in HYMOR, which provides global modal, non-modal, and freestream receptivity analyses for high-enthalpy hypersonic flows (Antón-Álvarez et al., 4 Apr 2026). Several thermochemical models are available for treatment of real-gas effects in high-enthalpy regimes: calorically perfect gas (CPG), Frozen-RTV, Chemical-RTV, Chemical-RTVE, and NonEq-RTVE (Antón-Álvarez et al., 4 Apr 2026). The base and disturbance equations are closed using effective thermodynamic and transport properties, represented via fast surrogate fits for computational efficiency (Antón-Álvarez et al., 4 Apr 2026).

A defining methodological feature is shock-fitting. The bow shock is treated as a zero-thickness, sharp discontinuity using a spline representation and Rankine–Hugoniot jump conditions, rather than a numerically smeared shock (Antón-Álvarez et al., 4 Apr 2026). In the linearized setting, perturbations of the jump conditions couple upstream and downstream disturbances through transmission matrices, and the shock-front displacement vector $K$8 is incorporated into the operator so that post-shock flow and shock motion are consistently coupled (Antón-Álvarez et al., 4 Apr 2026). The stated purpose is to ensure that interaction of infinitesimal disturbances with the shock reproduces the exact response predicted by linear interaction analysis.

The same linear operator supports three standard analyses. Modal stability seeks eigenpairs

$K$9

with $a_0$0 indicating exponential growth. Non-modal analysis maximizes finite-time growth in a Chu-energy-based norm through a generalized eigenvalue problem for

$a_0$1

Freestream receptivity analysis uses an extended upstream–downstream state,

$a_0$2

to compute the freestream disturbance that is most efficiently transmitted and amplified in the post-shock region (Antón-Álvarez et al., 4 Apr 2026).

The numerical implementation uses finite volume, second-order discretization on curvilinear structured grids, automatic linearization of the discrete operators, Arnoldi iteration for modal analysis, Lanczos or Arnoldi iteration for non-modal and receptivity analysis, and GPU acceleration for large sparse problems (Antón-Álvarez et al., 4 Apr 2026). In this sense, the modern small-disturbance framework is not a replacement for the asymptotic steady theory but a different reduction aimed at instability, amplification, and shock-mediated receptivity.

6. Extensions, applications, and interpretive issues

The scope of hypersonic small-disturbance ideas extends beyond slender straight wedges. For hypersonic potential flow onto a large curved wedge, a global shock wave attached to the tip and a smooth post-shock flow field can be constructed when the incoming Mach number is sufficiently large, without smallness assumptions on the height of the wedge or on $a_0$3-perturbation of a line (Hu et al., 2024). If the slope of the wedge tends to a positive limit, the shock slope tends to that of the self-similar straight-wedge case with the same limiting slope; if the slope of the wedge is parallel to the incoming flow at infinity, the strength of the shock diminishes to zero at infinity (Hu et al., 2024). This places classical similarity arguments in a larger geometric setting.

A separate extension treats hypersonic potential flow past a large curved convex wedge through a hodograph transformation. For $a_0$4, the free-boundary quasilinear problem is reduced to a linear equation in the hodograph plane, an approximate boundary is constructed from the asymptotic state, and uniform weighted Schauder estimates are used to obtain an error estimate for the free boundary (Hu, 5 Aug 2025). The shock polar in the hypersonic limit reduces to

$a_0$5

which is a circle in the $a_0$6-plane (Hu, 5 Aug 2025).

Related work on transition and vehicle dynamics uses small-disturbance linearization in ways that are adjacent to, but not identical with, the classical similarity equations. Direct numerical simulations on a Mach 6 ogive-cylinder consider linear and nonlinear disturbance evolution about a steady laminar base flow, with linear perturbations at $a_0$7 to $a_0$8 and nonlinear forcing at $a_0$9; Mack modes are reported as the dominant primary instabilities, with fundamental resonance as the dominant secondary instability mechanism (Goparaju et al., 2022). In guidance and control, nonlinear hypersonic vehicle dynamics can be linearized about an open-loop reference trajectory to obtain a linear time-varying system for a zero-sum Linear Quadratic Differential Game (Lee et al., 2021). In lateral flight dynamics, introducing time-variant pressure fluctuations into a simplified small-disturbance model leads to

$\bm{q}_0$0

a Mathieu differential equation (Wei et al., 2012).

These adjacent uses clarify a final interpretive point. “Hypersonic small-disturbance equations” do not denote a single immutable PDE system. In the cited literature, the phrase names a family of reductions: steady asymptotic similarity equations for compressible Euler or potential flow, measure-theoretic limit equations with singular concentrations, and linearized disturbance equations for stability, receptivity, and vehicle dynamics. What unifies them is the same organizing principle: in a hypersonic regime, the dominant flow behavior can often be isolated by retaining leading-order disturbance structure while discarding higher-order terms (Kang et al., 23 Aug 2025, Antón-Álvarez et al., 4 Apr 2026).