When a shock wave reflects obliquely, the three-shock (Mach) reflection configuration contains a small interaction zone — a few tens of molecular mean free paths wide — in which the classical inviscid Rankine–Hugoniot jump relations and three-shock theory fail: a finite viscous/kinetic 'non-Rankine–Hugoniot' transition zone replaces the singular triple point, and the local pressure overshoots the inviscid prediction by roughly 10%. A 2011 combined Navier–Stokes (NS) and direct simulation Monte Carlo (DSMC) study established this zone for a strong-shock case (M∞=4M_\infty=4, wedge angle Θw=25∘\Theta_w=25^\circ) and a weak-shock von Neumann paradox case (M∞=1.7M_\infty=1.7, Θw=13.5∘\Theta_w=13.5^\circ) at Kn‘≈‘`\approx`0.001, but left its quantitative theory open. The open problem is to determine the complete scaling laws of the zone (width, peak pressure excess, and internal field profiles as functions of free-stream Mach number, wedge angle, Knudsen/Reynolds number, and molecular gas properties), the general parametric boundary within which the continuum NS equations quantitatively reproduce kinetic (DSMC) internal structure at the shock intersection, the role of the viscous connector in unsteady transitions between regular and irregular reflection including its stability, and direct experimental confirmation of the zone's internal structure and pressure excess in planar reflection configurations.

Background

When an oblique shock wave meets a wedge or a plane of symmetry, the reflection takes one of two canonical configurations: regular reflection, in which the incident and reflected shocks meet on the surface, or irregular (Mach) reflection, in which the incident and reflected shocks are joined to a third shock — the Mach stem — at a nominal triple point. Inviscid theory predicts the geometry and post-shock states by shock polar analysis: each oblique shock maps the upstream state onto a curve of possible pressure/flow-deflection pairs using the one-dimensional Rankine–Hugoniot jump relations (conservation of mass, momentum, and energy across a discontinuity), and a three-shock solution exists only where the polars of the incident and reflected shocks intersect consistently. For weak incident shocks in air-like gases there is a parameter range — the von Neumann paradox regime — in which experiments and simulations show a Mach-reflection-like pattern although the inviscid three-shock construction admits no solution. Experiments in the early 1990s (Olim and Dewey, 1992) had already shown that weak Mach reflection requires a pressure excess over the classical three-shock prediction, motivating 'revised' three-shock solutions. The physical resolution is that the triple point is not a mathematical discontinuity: the region where the three shocks meet has a characteristic size of a few tens of mean free paths λ\lambda (the average distance between molecular collisions), so viscosity and molecular (kinetic) non-equilibrium smear the intersection into a finite transition zone. Here λ\lambda enters through the Knudsen number Kn, the ratio of mean free path to a macroscopic flow scale, while the Reynolds number Re measures inertia relative to viscous stresses; NS denotes the Navier–Stokes continuum equations and DSMC the direct simulation Monte Carlo method, a stochastic particle scheme that solves the Boltzmann kinetic equation and remains valid at finite Kn. In 2011 Shoev, Khotyanovsky, Bondar, Kudryavtsev, and Ivanov computed the triple-shock internal structure with high-order NS and with DSMC for steady reflection of a shock generated by a wedge: for a strong case (M∞=4M_\infty=4, Θw=25∘\Theta_w=25^\circ) and a weak von Neumann case (M∞=1.7M_\infty=1.7, Θw=13.5∘\Theta_w=13.5^\circ), both near Kn‘≈‘`\approx`0.001, the two methods agreed, and the numerical fields showed that away from the intersection the shocks follow their shock polars, while across a finite connector region the local (θ,p)(\theta, p) pairs (flow deflection angle and pressure) move continuously from the incident to the reflected polar, with a maximum pressure about 10% above the inviscid three-shock value (p/p∞=19.57p/p_\infty=19.57 at the triple point for the strong case). This continuous non-Rankine–Hugoniot zone explains mechanistically how an irregular reflection can exist where inviscid theory says it cannot. Subsequent work has added pieces: a 2019 DSMC study (Liu, Chen, Zhang, and Liu) found the zone size at Mach numbers 5–8 is set by shock thickness and independent of Mach number while the pressure hot spot grows with Mach number, with a boundary-layer integral predicting its magnitude; a 2023 Journal of Fluid Mechanics study (Shoev, Kudryavtsev, Khotyanovsky, Bondar) showed for a monatomic gas that in strong reflection the zone is self-similar in mean-free-path-normalized coordinates at all Reynolds numbers examined (physical size shrinking with Re, peak pressure excess constant near 10%), whereas in weak reflection the structure is not self-similar and converges to the inviscid three-shock prediction as Re grows; a 2024 JFM study (Bondar, Shoev, Timokhin) showed that at the regular-reflection point the Navier–Stokes–Fourier equations fail qualitatively and the regularized 13-moment kinetic approximation tracks DSMC much better; a 2025 Shock Waves study (Shoev, Nazarov, Timokhin) found that viscosity delays the steady regular-to-Mach-reflection transition at M=2M=2, Re=1000=1000; and a 2025 experimental/computational study (Higa, Matsunaga, Ogawa, Ohtani, Handa) imaged a non-Rankine–Hugoniot zone behind the Mach disk of an axisymmetric reflection at M∞=2.5M_\infty=2.5 using planar laser-induced fluorescence, attributing it to heat conduction across the slip line and cooling along the Rayleigh line. None of these, however, supplies the general quantitative theory of the zone that the 2011 study called for.

Problem Statement

Establish a validated quantitative theory of the viscous/kinetic non-Rankine–Hugoniot transition zone that replaces the singular triple point in irregular shock reflection, jointly determining: (i) the scaling laws of the zone — its width, peak pressure excess, and internal field profiles — as functions of free-stream Mach number, wedge angle, Knudsen/Reynolds number, and molecular gas properties (specific-heat ratio γ\gamma and transport coefficients), unified across the strong-shock, weak-shock, and von Neumann paradox regimes; (ii) the general parametric boundary in this same parameter space within which Navier–Stokes computations quantitatively reproduce kinetic (Boltzmann/DSMC) predictions for the triple-shock internal structure — the comparison quantities (peak pressure excess, zone width, and internal field profiles) and the numerical tolerance applied to each being stated explicitly, so that the boundary is a determinate level set rather than a qualitative judgement — a threshold law beyond the single examined point Kn‘≈‘`\approx`0.001 — together with the identification of which continuum approximation (Navier–Stokes–Fourier or higher-order moment systems) is adequate in each region; (iii) the role of the viscous/kinetic connector in unsteady reflection phenomena, i.e. how it affects time-dependent transitions between regular and irregular reflection and the stability of the transition zone itself; and (iv) direct experimental confirmation of the zone's internal structure and pressure excess in the immediate triple-point region of planar reflection configurations. These four deliverables jointly support one objective — replacing the inviscid singular triple point by a predictive, experimentally anchored, multi-scale physical description of the three-shock interaction — and none of them is separately resolved by the existing literature. Deliverable (iv) requires measurements that resolve the zone interior in a planar reflection configuration; comparisons against measurements that do not resolve the zone interior constrain at most the pressure excess.

References

  1. G. V. Shoev, D. V. Khotyanovsky, Ye. A. Bondar, A. N. Kudryavtsev, M. S. Ivanov, 'Numerical Study of Triple-Shock-Wave Structure in Steady Irregular Reflection', AIP Conference Proceedings (27th International Symposium on Rarefied Gas Dynamics), pp. 325–330, 2011. DOI: 10.1063/1.3562669, https://doi.org/10.1063/1.3562669
  2. G. V. Shoev, A. N. Kudryavtsev, D. V. Khotyanovsky, Ye. A. Bondar, 'Viscous effects in Mach reflection of shock waves and passage to the inviscid limit', Journal of Fluid Mechanics 959, A2, 2023. DOI: 10.1017/jfm.2023.89, https://doi.org/10.1017/jfm.2023.89
  3. Y. A. Bondar, G. V. Shoev, M. Y. Timokhin, 'Regular reflection of shock waves in steady flows: viscous and non-equilibrium effects', Journal of Fluid Mechanics 984, A10, 2024. DOI: 10.1017/jfm.2024.112, https://doi.org/10.1017/jfm.2024.112
  4. Haochen Liu, Hao Chen, Bin Zhang, Hong Liu, 'Effects of Mach Number on Non-Rankine–Hugoniot Shock Zone of Mach Reflection', Journal of Spacecraft and Rockets 56, pp. 761–770, 2019. DOI: 10.2514/1.A34251, https://doi.org/10.2514/1.A34251
  5. Yoshitaka Higa, Masanobu Matsunaga, Hideaki Ogawa, Kiyonobu Ohtani, Taro Handa, 'Study on non-Rankine-Hugoniot zone in axisymmetric supersonic flow', Physics of Fluids 37, 056110, 2025. DOI: 10.1063/5.0267059, https://doi.org/10.1063/5.0267059
  6. G. Shoev, V. Nazarov, M. Timokhin, 'Numerical study of viscous effects on the formation of the Mach reflection in a transition from regular reflection', Shock Waves 35, pp. 109–124, 2025. DOI: 10.1007/s00193-024-01213-1, https://doi.org/10.1007/s00193-024-01213-1
  7. M. Olim, J. M. Dewey, 'A revised three-shock solution for the Mach reflection of weak shocks (1.1 < Mi < 1.5)', Shock Waves 2, pp. 167–176, 1992. DOI: 10.1007/BF01414639, https://doi.org/10.1007/BF01414639
  8. Yazhong Jiang, Xuxu Sun, Jie Niu, Jun Zhang, 'Study of shock–shock interactions in rarefied flows using direct simulation Monte Carlo method', Aerospace Science and Technology 173, 111768, 2026. DOI: 10.1016/j.ast.2026.111768, https://doi.org/10.1016/j.ast.2026.111768