Mach: Multifaceted Applications in Science

• Mach is a dimensionless measure comparing flow speed to a characteristic wave speed, crucial in compressible fluid dynamics and shock analysis.
• The term extends to specialized applications such as interferometry, thermal clocks, and innovative computing systems, highlighting its multidisciplinary relevance.
• Diverse definitions of Mach across fluid dynamics, plasma physics, and astrophysics drive unique numerical challenges and experimental insights in various regimes.

In the literature represented here, mach is not a single concept but a cluster of technical usages. In fluid dynamics and plasma physics it most often denotes a Mach number, that is, a ratio of flow speed to a characteristic propagation speed, and it appears in related constructions such as low-Mach limits, Mach reflection, Mach stems, and Mach cones. The same string also appears in other specialized contexts, including the Mach-Zehnder-Fano interferometer, the Mach thermal clock, and several recent software systems named MACH or mach (Hu et al., 2012, Sulaiman et al., 2015, Xu et al., 2012, Milburn et al., 2017, Spring et al., 29 Apr 2025, Essendelft et al., 18 Jun 2025, Guan et al., 6 Apr 2026).

1. Mach number as a control parameter across continua

Across the cited literature, the most common usage is the dimensionless comparison of flow speed with a characteristic wave speed. In compressible Euler flow, the Mach number is written

$M = \frac{|u|}{c},$

and the low-Mach regime is the regime $M \ll 1$ (Hu et al., 2012). In collisionless plasma shocks, the emphasis can shift from sonic to Alfvénic scaling, with

$M_A = \frac{\sqrt{\mu_0 P_{\mathrm{dyn}}}}{B_u},$

while sonic and fast magnetosonic Mach numbers are also distinguished (Sulaiman et al., 2015). In interstellar MHD turbulence, the paper on Tsallis statistics uses the volume-averaged forms

$$\mathcal{M}_s \equiv \left\langle \frac{|{\bf v}|}{C_s} \right\rangle,\qquad \mathcal{M}_A \equiv \left\langle \frac{|{\bf v}|}{v_A} \right\rangle,$$

with $v_A = |{\bf B}|/\sqrt{\rho}$ (Tofflemire et al., 2011). In radiative turbulent mixing layers, the relevant definition is

$\mathcal M \equiv \frac{v_{\rm hot}}{c_{\rm s,h}},$

and the regime transition occurs near $\mathcal M \sim 1$ (Yang et al., 2022).

Context	Definition	Representative source
Compressible Euler flow	$M=|u|/c$	(Hu et al., 2012)
Euler–Korteweg scaling	pressure and capillarity enter as $1/M^2$ terms	(Giesselmann, 2013)
Collisionless plasma shocks	$M_A=\sqrt{\mu_0 P_{\mathrm{dyn}}}/B_u$	(Sulaiman et al., 2015)
ISM turbulence	$M \ll 1$0, $M \ll 1$1	(Tofflemire et al., 2011)
Radiative TMLs	$M \ll 1$2	(Yang et al., 2022)

The low-Mach limit is not merely definitional; it changes the asymptotic structure of the equations. In the Euler–Korteweg model, the pressure and capillarity forces appear with a singular prefactor $M \ll 1$3, and the low-Mach limit yields a frozen leading-order density satisfying

$M \ll 1$4

together with a divergence constraint

$M \ll 1$5

The semi-implicit finite-difference scheme in that paper is described as asymptotic preserving, because for fixed discretization parameters it converges to a stable discretization of the incompressible limit as $M \ll 1$6 (Giesselmann, 2013).

A related computational problem appears in stellar hydrodynamics. Conventional compressible Godunov-type solvers were developed primarily for transonic and supersonic problems and “show difficulties in representing already moderately low Mach numbers,” typically below about $M \ll 1$7, depending on scheme details and resolution. The proposed remedy is a modified Roe solver whose numerical dissipation has the correct low-Mach scaling while retaining the compressible Euler equations (Miczek et al., 2014).

2. Reflection, stems, and low-Mach numerical pathologies

In shock-capturing numerics, Mach-related phenomena include both physical reflection structures and numerical artifacts. Hu and Adams connect two defects of Godunov-type schemes that are often treated separately: the inaccurate pressure profile in low-Mach flow and the carbuncle phenomenon in high-Mach shock calculations. For Roe fluxes, the dissipative part contains an advection dissipation proportional to $M \ll 1$8 and an acoustic dissipation proportional to $M \ll 1$9. When $M_A = \frac{\sqrt{\mu_0 P_{\mathrm{dyn}}}}{B_u},$0, the ratio is effectively of order $M_A = \frac{\sqrt{\mu_0 P_{\mathrm{dyn}}}}{B_u},$1, and the scheme produces pressure fluctuations of order $M_A = \frac{\sqrt{\mu_0 P_{\mathrm{dyn}}}}{B_u},$2 even when physically they should be $M_A = \frac{\sqrt{\mu_0 P_{\mathrm{dyn}}}}{B_u},$3. In a Mach-10 double-Mach-reflection test, two opposite modifications of the Roe eigenvalues—Roe-M1, which decreases acoustic dissipation, and Roe-M2, which increases advection dissipation—both eliminate the kinked Mach stem, suggesting that the carbuncle is strongly related to the non-comparability of acoustic and advection contributions rather than to insufficient dissipation alone (Hu et al., 2012).

The classical physical problem of Mach reflection is studied directly for inviscid flow over a wedge with a higher-order discontinuous Galerkin method and overset grids. The paper distinguishes regular reflection (RR) from Mach reflection (MR) and resolves both the detachment criterion and the Von Neumann condition. For $M_A = \frac{\sqrt{\mu_0 P_{\mathrm{dyn}}}}{B_u},$4, impulsive start gives the RR-to-MR transition between $M_A = \frac{\sqrt{\mu_0 P_{\mathrm{dyn}}}}{B_u},$5 and $M_A = \frac{\sqrt{\mu_0 P_{\mathrm{dyn}}}}{B_u},$6, in essentially perfect agreement with the reported three-shock-theory detachment angle $M_A = \frac{\sqrt{\mu_0 P_{\mathrm{dyn}}}}{B_u},$7. For the same Mach number, continuation from an established MR state preserves a nonzero Mach stem down to between $M_A = \frac{\sqrt{\mu_0 P_{\mathrm{dyn}}}}{B_u},$8 and $M_A = \frac{\sqrt{\mu_0 P_{\mathrm{dyn}}}}{B_u},$9, close to the reported Von Neumann angle $$\mathcal{M}_s \equiv \left\langle \frac{|{\bf v}|}{C_s} \right\rangle,\qquad \mathcal{M}_A \equiv \left\langle \frac{|{\bf v}|}{v_A} \right\rangle,$$0. For $$\mathcal{M}_s \equiv \left\langle \frac{|{\bf v}|}{C_s} \right\rangle,\qquad \mathcal{M}_A \equiv \left\langle \frac{|{\bf v}|}{v_A} \right\rangle,$$1, the corresponding thresholds are $$\mathcal{M}_s \equiv \left\langle \frac{|{\bf v}|}{C_s} \right\rangle,\qquad \mathcal{M}_A \equiv \left\langle \frac{|{\bf v}|}{v_A} \right\rangle,$$2–$$\mathcal{M}_s \equiv \left\langle \frac{|{\bf v}|}{C_s} \right\rangle,\qquad \mathcal{M}_A \equiv \left\langle \frac{|{\bf v}|}{v_A} \right\rangle,$$3 and $$\mathcal{M}_s \equiv \left\langle \frac{|{\bf v}|}{C_s} \right\rangle,\qquad \mathcal{M}_A \equiv \left\langle \frac{|{\bf v}|}{v_A} \right\rangle,$$4–$$\mathcal{M}_s \equiv \left\langle \frac{|{\bf v}|}{C_s} \right\rangle,\qquad \mathcal{M}_A \equiv \left\langle \frac{|{\bf v}|}{v_A} \right\rangle,$$5, matching theoretical values $$\mathcal{M}_s \equiv \left\langle \frac{|{\bf v}|}{C_s} \right\rangle,\qquad \mathcal{M}_A \equiv \left\langle \frac{|{\bf v}|}{v_A} \right\rangle,$$6 and $$\mathcal{M}_s \equiv \left\langle \frac{|{\bf v}|}{C_s} \right\rangle,\qquad \mathcal{M}_A \equiv \left\langle \frac{|{\bf v}|}{v_A} \right\rangle,$$7. The finite Mach stem height $$\mathcal{M}_s \equiv \left\langle \frac{|{\bf v}|}{C_s} \right\rangle,\qquad \mathcal{M}_A \equiv \left\langle \frac{|{\bf v}|}{v_A} \right\rangle,$$8 functions as the order parameter of the RR/MR transition and was resolved down to $$\mathcal{M}_s \equiv \left\langle \frac{|{\bf v}|}{C_s} \right\rangle,\qquad \mathcal{M}_A \equiv \left\langle \frac{|{\bf v}|}{v_A} \right\rangle,$$9 (Kochi et al., 2023).

A related astrophysical use of Mach stem appears in intersecting bow shocks. In two-dimensional simulations of dense clumps in a supersonic wind with $v_A = |{\bf B}|/\sqrt{\rho}$0, the critical included angle for Mach stem formation is reported as approximately

$v_A = |{\bf B}|/\sqrt{\rho}$1

with numerical values $v_A = |{\bf B}|/\sqrt{\rho}$2 for $v_A = |{\bf B}|/\sqrt{\rho}$3, $v_A = |{\bf B}|/\sqrt{\rho}$4 for $v_A = |{\bf B}|/\sqrt{\rho}$5, $v_A = |{\bf B}|/\sqrt{\rho}$6 for $v_A = |{\bf B}|/\sqrt{\rho}$7, and $v_A = |{\bf B}|/\sqrt{\rho}$8 for $v_A = |{\bf B}|/\sqrt{\rho}$9. The inferred critical clump separations are $\mathcal M \equiv \frac{v_{\rm hot}}{c_{\rm s,h}},$0, $\mathcal M \equiv \frac{v_{\rm hot}}{c_{\rm s,h}},$1, and $\mathcal M \equiv \frac{v_{\rm hot}}{c_{\rm s,h}},$2 for $\mathcal M \equiv \frac{v_{\rm hot}}{c_{\rm s,h}},$3, $\mathcal M \equiv \frac{v_{\rm hot}}{c_{\rm s,h}},$4, and $\mathcal M \equiv \frac{v_{\rm hot}}{c_{\rm s,h}},$5, respectively; the $\mathcal M \equiv \frac{v_{\rm hot}}{c_{\rm s,h}},$6 case is unstable (Hansen et al., 2014).

Mach number also controls instability around spiked blunt bodies. For a round-tip aerospike with $\mathcal M \equiv \frac{v_{\rm hot}}{c_{\rm s,h}},$7 and $\mathcal M \equiv \frac{v_{\rm hot}}{c_{\rm s,h}},$8, the paper reports pulsation mode at Mach 2 and 3, an almost stable state at Mach 4, and oscillatory mode at Mach 5, 6, and 7. The dominant frequency rises with Mach number, while the reported Strouhal-number boundary separating pulsation from oscillation for this configuration is about $\mathcal M \equiv \frac{v_{\rm hot}}{c_{\rm s,h}},$9 (Vashishtha et al., 2021).

3. Collisionless shocks, dense matter, and astrophysical high-Mach regimes

In collisionless plasma physics, Mach number organizes both structure and nonstationarity. Cassini observations of Saturn’s quasi-perpendicular bow shock probe a regime spanning two orders of magnitude in $\mathcal M \sim 1$0. The paper reports evidence for cyclic reformation at a timescale $\mathcal M \sim 1$1, equivalently $\mathcal M \sim 1$2–$\mathcal M \sim 1$3, consistent with specularly reflected ions. Reformation becomes common in the highest-$\mathcal M \sim 1$4 subset, especially for $\mathcal M \sim 1$5, but high $\mathcal M \sim 1$6 is described as necessary rather than sufficient. At a given $\mathcal M \sim 1$7, reforming shocks also exhibit systematically larger magnetic amplification $\mathcal M \sim 1$8 than non-reforming shocks (Sulaiman et al., 2015).

In dense-matter kinetic plasma physics, high Mach number is necessary but not monotonically beneficial. For head-on collisions of two quantum-degenerate deuterium jets, the scanned collision-speed range $\mathcal M \sim 1$9 to $M=|u|/c$0 corresponds to $M=|u|/c$1 to $M=|u|/c$2. The hydrodynamic model predicts monotonic approach to the strong-shock compression limit $M=|u|/c$3, but first-principles kinetic simulations show that the compression peaks at about

$M=|u|/c$4

around $M=|u|/c$5, $M=|u|/c$6, and then decreases at higher Mach number because the shock thickness scales roughly as $M=|u|/c$7. For $M=|u|/c$8, $M=|u|/c$9, and $1/M^2$0, the paper estimates $1/M^2$1 and a shock thickness $1/M^2$2, comparable to the collision region size (Zhang et al., 2023).

High-Mach radiative turbulent mixing layers also change character near transonic conditions. For $1/M^2$3, the TML behaves as a single cooling-and-mixing structure and the surface brightness and intermediate-temperature ion columns scale approximately as $1/M^2$4. For $1/M^2$5, the layer separates into a Mach-independent mixing zone and an expanding turbulent zone; the cooling is then predominantly balanced by turbulent dissipation rather than enthalpy consumption, and both surface brightness and ion columns saturate at $1/M^2$6. In this high-Mach regime, inflow velocities and hot-gas entrainment are suppressed and cold gas evaporates rather than growing by condensation (Yang et al., 2022).

Astrophysical Mach cones appear in two very different settings. In viscous heavy-ion collisions, Mach cones generated by jets in a full $1/M^2$7-dimensional expanding medium are visible for small viscosity, especially $1/M^2$8, but are progressively smeared out as $1/M^2$9 rises to $M_A=\sqrt{\mu_0 P_{\mathrm{dyn}}}/B_u$0 and $M_A=\sqrt{\mu_0 P_{\mathrm{dyn}}}/B_u$1. The paper argues that a double peak in azimuthal correlations is not by itself a reliable Mach-cone signature, because averaging over many jet paths can generate a double peak from deflected head shocks and diffusion wakes rather than from a clean Mach angle (Bouras et al., 2014). In the Virgo cluster, diffuse X-ray halos and tails of NGC 4569, NGC 4388, and NGC 4501 are interpreted with simple Mach-cone geometries inferred from three-dimensional galaxy velocities; for NGC 4569, the derived opening angle is difficult to reconcile with a Mach number based on sound speed alone, and the paper therefore argues that a magnetosonic Mach number is more appropriate, implying magnetic fields of order a few $M_A=\sqrt{\mu_0 P_{\mathrm{dyn}}}/B_u$2G for $M_A=\sqrt{\mu_0 P_{\mathrm{dyn}}}/B_u$3 (Wezgowiec et al., 2011).

4. Mach-like wakes and cone geometries in dispersive media

In dispersive wave media, Mach-type geometry need not arise from ordinary acoustic cones. In a two-dimensional hydrodynamic Fermi sea, the plasmon dispersion

$M_A=\sqrt{\mu_0 P_{\mathrm{dyn}}}/B_u$4

interpolates between $M_A=\sqrt{\mu_0 P_{\mathrm{dyn}}}/B_u$5 and $M_A=\sqrt{\mu_0 P_{\mathrm{dyn}}}/B_u$6. The corresponding wake geometry is governed by

$M_A=\sqrt{\mu_0 P_{\mathrm{dyn}}}/B_u$7

with sharp changes at $M_A=\sqrt{\mu_0 P_{\mathrm{dyn}}}/B_u$8 and $M_A=\sqrt{\mu_0 P_{\mathrm{dyn}}}/B_u$9. For $M \ll 1$00, no stationary wake occurs. For $M \ll 1$01, the wake is confined to the standard Mach sector of angle

$M \ll 1$02

and consists only of transverse wavefronts. For $M \ll 1$03, an additional Kelvin-like outer region appears, and the full wake angle becomes

$M \ll 1$04

The paper therefore describes the resulting structure as a Kelvin-Mach wake (Kolomeisky et al., 2017).

A formally analogous but physically distinct transition appears in surface-gravity wakes of ships. Kelvin’s classical prediction gives a constant half-angle

$M \ll 1$05

independent of speed, but the cited analysis of airborne images shows that the visible wake angle narrows at high speed in a Mach-like way. The control parameter is the hull Froude number

$M \ll 1$06

and the transition occurs at approximately

$M \ll 1$07

Above this value, the finite hull size excludes the Kelvin-cusp wavelength, and the asymptotic wake angle becomes

$M \ll 1$08

that is, $M \ll 1$09 for fixed $M \ll 1$10 (Rabaud et al., 2013).

These two papers jointly show that “Mach-like” wake narrowing in the literature is not restricted to non-dispersive acoustics. A plausible implication is that the term can denote a geometric scaling law even when the underlying propagation medium is dispersive, provided a minimum or effective signal speed still organizes the far-field pattern.

5. Mach in interferometry and thermal clocks

In photonics, the term appears in the nonlinear Mach-Zehnder-Fano interferometer (MZFI). This device is a Mach-Zehnder interferometer whose arm contains a side-coupled nonlinear Fano defect. The interaction between ordinary loop resonances and the defect’s Fano resonance produces hybrid resonant states with stronger field localization and sharper transmission features than either a conventional Mach-Zehnder interferometer or a standalone Fano structure. The paper reports enhanced bistability, “low threshold 100% switching operation,” and a figure of merit

$M \ll 1$11

that can be enhanced by more than 60 times at the hybrid resonance $M \ll 1$12 relative to the conventional Fano resonance $M \ll 1$13 (Xu et al., 2012).

A very different usage appears in quantum thermodynamics. The paper on a quantum optomechanical Mach clock explicitly invokes the thermal clock “first introduced by Ernst Mach,” in which a temperature difference between cooling bodies functions as a relational time variable. The proposed realization uses two optical cavities and a mechanical resonator coupled to a heat bath. After adiabatic elimination of the mechanics, the effective cavity dynamics becomes

$M \ll 1$14

so photons are transferred irreversibly between cavities through thermally biased jump processes. The clock variable is the cavity photon-number difference, which in a semiclassical limit obeys

$M \ll 1$15

The paper contrasts this ensemble-average relaxation law with the stochastic record of a single continuously measured system and studies how quantum fluctuations modify the ideal thermal-clock picture (Milburn et al., 2017).

6. MACH and mach as names of contemporary technical systems

Recent literature also uses MACH or mach as the proper name of several systems unrelated to compressible-flow Mach numbers.

Name	Expansion or description	Representative source
MACH	Multi-Agent Coordination for RSU-centric Handovers	(Spring et al., 29 Apr 2025)
MACH	Multiple-Architecture Compiler for Advanced Computing Hardware	(Essendelft et al., 18 Jun 2025)
mach	open-source, GPU-accelerated delay-and-sum ultrasound beamformer	(Guan et al., 6 Apr 2026)

In vehicular edge computing, MACH is a decentralized RSU-centric handover method in which Road Side Units, rather than a centralized coordinator or vehicles, decide where and when tasks should be offloaded. The system model uses vehicle state $M \ll 1$16, RSU state $M \ll 1$17, and a product QoS model combining distance-based and load-based terms. The evaluation uses a Python Mesa simulator driven by the Créteil Roundabout Dataset, with 4-RSU and 9-RSU configurations and compute capacities modeled after NVIDIA Tesla T4 GPUs. The paper notes a naming inconsistency in which some text refers to ARHC instead of MACH (Spring et al., 29 Apr 2025).

In compiler design, MACH denotes the Multiple-Architecture Compiler for Advanced Computing Hardware, intended for massively parallel, spatial, dataflow architectures such as the Wafer Scale Engine. Its central abstraction is a hardware-agnostic Virtual Machine composed of a controller and one or more workers. The compiler stack includes object-oriented data structures, an intermediate language, an intermediate representation graph, and a backend that lowers NumPy-like tensor programs to Cerebras’ Tungsten and Paint. The same conceptual program can also run on traditional unified-memory hardware (Essendelft et al., 18 Jun 2025).

In ultrasound imaging, mach is an open-source Python/CUDA beamformer for ultrafast ultrasound. The paper reports 1.1 trillion beamforming points per second on a consumer NVIDIA GeForce RTX 5090, a runtime of 0.23 ms on the PyMUST rotating-disk benchmark, and numerical agreement with other beamformers at below $M \ll 1$18 dB for Power Doppler and below $M \ll 1$19 dB for B-mode. The key implementation idea is a hybrid delay strategy: transmit wavefront arrival times are precomputed and stored once, while receive delays are computed on the fly within CUDA thread blocks and reused across frames via shared memory (Guan et al., 6 Apr 2026).

Taken together, these usages show that mach functions in current technical writing as a genuinely polysemous marker. In some fields it denotes a dimensionless control parameter for compressibility and wave propagation; in others it labels geometric shock structures, interferometric architectures, thermodynamic clocks, or modern software systems. The commonality is therefore lexical rather than conceptual, and the meaning is determined almost entirely by disciplinary context.