Hypersonic Similarity Law in Compressible Flows

Updated 9 July 2026

Hypersonic Similarity Law is a set of asymptotic principles for high-speed compressible flows that achieve Mach-number independence by scaling shock and post-shock fields.
It uses fixed similarity parameters like K = M∞θ to reduce the full Euler and potential flow equations to hypersonic small-disturbance models with quantifiable convergence rates.
The framework extends to measure-valued solutions with Dirac masses, rigorously justifying classical Newton and Newton–Busemann pressure laws in varied geometrical settings.

Searching arXiv for recent and foundational papers on hypersonic similarity laws, Newton–Busemann/Newtonian limits, and related boundary-layer similitude. Hypersonic similarity law denotes a class of asymptotic principles for high-speed compressible flow in which, after a problem-specific scaling, the flow structure becomes independent of the Mach number as an independent parameter and is instead controlled by a reduced set of similarity variables. In the classical wedge and slender-body setting, this principle is formulated as the Mach-number independence principle or Van Dyke’s similarity theory: for a fixed hypersonic similarity parameter, the scaled shock and post-shock fields converge to a hypersonic small-disturbance description as $M_\infty \to \infty$ (Kuang et al., 2019). In a distinct but related line of work on steady Euler flows in the hypersonic limit, the limiting solution is not generally a regular function but a Radon measure with Dirac mass supported on the body surface, and classical Newtonian or Newton–Busemann pressure laws emerge from that singular limit (Jin et al., 2019). More recent work extends similarity beyond classical nondimensional-group matching by using Lie equivalence transformations of the similarity-reduced boundary-layer equations, while rarefied blunt-body studies delineate a regime in which single-scale similarity breaks down (Bader et al., 19 Jun 2026, Roohi et al., 16 May 2026).

1. Core definitions and parameters

The most common formulation in two-dimensional wedge theory uses the hypersonic similarity parameter

$K = M_\infty \theta,$

where $\theta$ is the wedge angle and $M_\infty$ is the upstream Mach number. In this form, the law states that for a fixed $K$ , when the Mach number is sufficiently large, the governing equations after scaling are approximated by the hypersonic small-disturbance equation, and the scaled flowfields become similar (Kuang et al., 2021). The same principle is described in the large-data potential-flow literature as the assertion that, if the hypersonic similarity parameter $K$ is fixed, the shock solution structures after scaling are consistent when $M_\infty$ is sufficiently large (Kuang et al., 2019).

For slender-body formulations, the similarity parameter is written as

$K \doteq M_\infty \tau,$

with $\tau$ the slenderness ratio. In this setting, as $\tau \to 0$ and $K = M_\infty \theta,$ 0 with $K = M_\infty \theta,$ 1 fixed, the scaled flow field depends only on $K = M_\infty \theta,$ 2 and the adiabatic index $K = M_\infty \theta,$ 3, not on the body shape and Mach number independently (Kang et al., 23 Aug 2025).

A complementary hypersonic-limit formulation arises in non-isentropic steady Euler flow past a straight wedge. There, for fixed total energy $K = M_\infty \theta,$ 4, the incoming Mach number and the adiabatic exponent are linked by

$K = M_\infty \theta,$ 5

so the hypersonic limit $K = M_\infty \theta,$ 6 may be taken mathematically as $K = M_\infty \theta,$ 7 (Qu et al., 2019). This formulation is not merely a reparameterization: it leads directly to measure-valued limits with concentration on the body.

These uses of the term are related but not identical. One concerns convergence of scaled entropy solutions to hypersonic small-disturbance equations for fixed similarity parameter; another concerns singular measure limits of the steady Euler equations in which the shock layer collapses onto the body. The literature treats both as rigorous expressions of hypersonic similarity.

2. Classical scaling and hypersonic small-disturbance equations

In the slender-wedge literature, the starting point is a scaled steady compressible flow system in which the transverse coordinate and disturbance variables are normalized by a small parameter $K = M_\infty \theta,$ 8. A representative scaling is

$K = M_\infty \theta,$ 9

with fixed similarity parameter $\theta$ 0 (Kuang et al., 2021). As $\theta$ 1, the full equations reduce to a simpler hypersonic small-disturbance system. In the potential-flow formulation, the limiting equations are

$\theta$ 2

which define the universal limit problem for fixed $\theta$ 3 (Kuang et al., 2019).

For full Euler flow past slender bodies, the reduced system is likewise explicit. In the two-dimensional curved-wedge setting, the scaled variables

$\theta$ 4

lead, as $\theta$ 5, to a hypersonic small-disturbance Euler system together with the state relation

$\theta$ 6

and slip boundary condition $\theta$ 7 on the scaled boundary (Kang et al., 23 Aug 2025).

The significance of this reduction is twofold. First, it explains the phrase “Mach-number independence principle”: the reduced equations do not retain $\theta$ 8 as an independent control parameter. Second, it identifies the precise limit problem to which rigorous $\theta$ 9 convergence results can be attached. In the cited work, the similarity law is therefore not a heuristic collapse alone, but a statement about convergence between entropy solutions of two distinct initial-boundary value problems.

3. Radon measure solutions and Newtonian pressure laws

A major development in the mathematical theory is the introduction of Radon measure solutions for steady Euler boundary-value problems in the hypersonic limit. Classical weak solutions are inadequate when the shock layer collapses to zero thickness and the density, momentum, and energy concentrate on the body surface. In the straight-wedge problem, the limiting density measure is written as

$M_\infty$ 0

with analogous decompositions for momentum and energy, and the explicit density weight

$M_\infty$ 1

along the wedge (Qu et al., 2019).

In this limit, the surface pressure obeys Newton’s sine-squared law,

$M_\infty$ 2

which the measure-theoretic construction justifies rigorously for two-dimensional straight wedges (Qu et al., 2019). For straight cones, the corresponding conical hypersonic limit yields

$M_\infty$ 3

again as the limit of weak entropy solutions converging vaguely as Radon measures to a singular boundary-supported solution (Li et al., 2024).

For two-dimensional ramps, the analogous result is the Newton–Busemann pressure law. In the ramp setting, the cited work proposes rigorous definitions of Radon measure solutions for steady compressible Euler boundary-value problems modeling hypersonic-limit inviscid flows passing two-dimensional ramps, and proves the Newton–Busemann pressure law of drags of body in hypersonic flow (Jin et al., 2019). The same framework constructs measure solutions with density containing Dirac measures supported on curves and treats interactions with still gas and pressureless jets, including examples of blow up of certain measure solutions (Jin et al., 2019).

This body of work alters a common interpretation of hypersonic similarity. In these Euler problems, similarity is not only a statement about regular asymptotic profiles. The limit may be singular, and the correct mathematical objects are measures with concentrated mass on the wall, wedge, ramp, or cone. The classical Newtonian and Newton–Busemann laws then emerge from the weak formulation of the conservation laws rather than from an external impact argument.

4. Rigorous convergence, optimal rates, and well-posedness

The potential-flow literature provides the first rigorous global validation of the Mach-number independence principle for large data in two-dimensional steady potential flow. For fixed $M_\infty$ 4 and sufficiently large $M_\infty$ 5, a modified Glimm scheme yields global entropy solutions for the scaled system, and as $M_\infty$ 6 these solutions converge to the entropy solution of the corresponding hypersonic small-disturbance problem (Kuang et al., 2019).

Subsequent work establishes quantitative rates. For inviscid steady irrotational Euler flow over a two-dimensional Lipschitz slender wedge in $M_\infty$ 7, the $M_\infty$ 8 difference between the scaled solution and the small-disturbance solution is of order $M_\infty$ 9,

$K$ 0

and this rate is identified as the same as that predicted by the Newtonian-Busemann law (Kuang et al., 2021).

For steady compressible full Euler flow over two-dimensional slender Lipschitz wedges, the optimal convergence rate is again

$K$ 1

under smallness assumptions on the total variation of the initial data and the tangential slope function of the wedge boundary. The same work establishes global BV entropy solutions and $K$ 2-stability, and describes the result as the first mathematical comparison of two solutions of the compressible Euler equations with characteristic boundary conditions (Chen et al., 2023).

Large-data potential-flow theory refines the rate further in a parameter-dependent form. For two-dimensional steady potential flows with large data past a straight wedge, one obtains

$K$ 3

with optimal order $K$ 4 (Chen et al., 2024).

These theorems show that the phrase “similarity law” has a precise analytical content: it can mean convergence in $K$ 5, sharp $K$ 6 error bounds, existence of Lipschitz continuous solution maps, and optimality results obtained through explicit constructions. The methods are correspondingly technical, including wave-front tracking, Standard Riemann Semigroups, detailed comparison of Riemann solvers, and stability theory for characteristic boundaries.

5. Extensions: curved bodies, friction, and group-theoretic similitude

The classical straight-wedge paradigm has been extended in several directions. For hypersonic potential flow past a curved smooth wedge with neither smallness assumption on the height of the wedge nor that it is a BV perturbation of a line, there exists a global shock wave attached to the tip together with a smooth flow field between it and the wedge when the incoming Mach number is sufficiently large. If the slope of the wedge has a positive limit as $K$ 7, the slope of the shock tends to that of the self-similar straight-wedge case with the same limiting slope; if 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 suggests that self-similar downstream structure can persist even beyond the small-disturbance, near-straight geometry usually associated with classical similarity arguments.

Skin friction provides another extension. For stationary non-isentropic compressible Euler equations on hypersonic-limit flows passing ramps with frictions on their boundaries, Radon measure solutions again contain Dirac measures on the ramp boundary, and the Newton–Busemann law is generalized to include drag and lift distributions under three friction models (Qu et al., 2023). In the full-stick limit $K$ 8, the normal and tangential force laws reduce to

$K$ 9

and the normal law matches the Newton “sine-squared” law for a straight ramp (Qu et al., 2023).

A different extension appears in hypersonic stagnation-point boundary layers. There, Lie equivalence symmetry analysis treats the thermochemical property laws $K$ 0 as transformable quantities, not fixed coefficients. The resulting invariants yield nonlinear maps from laboratory to flight solutions, most notably

$K$ 1

for the similarity-reduced non-dimensional temperature, with the map constructed from the actual lab- and flight-side property laws (Bader et al., 19 Jun 2026). In this framework, classical similarity appears as a particular case, while the broader equivalence group allows extrapolation even when standard similitude groups are not exactly matched.

Taken together, these developments show that hypersonic similarity law is no longer restricted to straight inviscid wedges with frictionless walls and fixed property models. The term now encompasses geometric generalization, wall-friction generalization, and transformation-based extrapolation in thermochemically complex boundary layers.

A closely related similarity result appears in hypersonic boundary-layer theory for blunt-nosed bodies. Using self-similar boundary-layer analysis, the ratio of skin-friction coefficient to heat-transfer coefficient is found to vary linearly with the local wall slope angle,

$K$ 2

with an explicit expression for circular cylinders and extensions to other nose shapes and chemically reacting gases (Chen et al., 2015). In the near-continuum regime, Burnett-based corrections lead to

$K$ 3

and, for circular cylinders in the transition regime, a bridge function is proposed: $K$ 4 This is not the same law as Mach-number independence for wedges, but it is a similarity relation in hypersonic aerothermodynamics built on self-similar structure and asymptotic scaling (Chen et al., 2015).

The limits of similarity are equally important. In rarefied hypersonic bow shocks over a circular cylinder, direct simulation Monte Carlo data show that as $K$ 5 increases, a sharp density-gradient ridge is replaced by a broad kinetic compression layer, and rarefied bow-shock inflation is a coupled compression–relaxation process rather than a single-scale rescaling of a continuum-like shock (Roohi et al., 16 May 2026). At low rarefaction, the continuum normal-shock density ratio remains a useful reference compression scale; at higher Knudsen number, the measured standoff growth is governed primarily by the kinetic mean free path, and density-based registration no longer collapses Mach-number and thermal profiles into a one-parameter structure (Roohi et al., 16 May 2026).

This limitation clarifies a persistent misconception. Hypersonic similarity is not a universal statement that all high-speed flowfields collapse under a single rescaling. The rigorous results are regime-specific: attached slender-body inviscid flows, measure-valued hypersonic Euler limits, similarity-reduced boundary layers with transformed property laws, and near-continuum blunt-body analogies each admit their own precise formulation. Outside those settings—particularly in rarefied blunt-body flow—the hypothesis of a single coherent similarity structure can fail.