In front of a blunt body moving supersonically, a detached bow shock forms whose nose curvature controls both the shock standoff distance and the entire near-axis flow. The analytic stagnation-line (hodograph-like) treatment of Keshet & Naor reproduces measured and simulated flow for a sphere once the normalized shock curvature is known, but the curvature itself is only interpolated between weak- and strong-shock limits: the power-law exponent β≈1/2 in the successful fit is empirical, its dependence on the adiabatic index and on body shape is undetermined, and the conjectured extension to non-axisymmetric bodies via the stagnation streamline is unproven. The open problem is to derive the shock-curvature function from the steady compressible Euler equations, including the origin and body-dependence of the exponent, and to establish or refute the stagnation-streamline generalization to non-axisymmetric bodies.

Background

When a blunt body moves through a compressible inviscid gas faster than the sound speed, a detached bow shock forms ahead of it at a standoff distance from the body's nose. Predicting this shock — its standoff distance and its curvature at the nose — is a classical problem of gas dynamics with applications from hypersonic vehicles to planets in the solar wind and gas clumps in galaxy clusters. Most practical treatments use empirical shock-shape correlations (e.g. Billig 1967) or case-by-case numerical simulation; the classic engineering fits give the shock shape for spheres and cylinders as functions of Mach number, but their parameters are fitted to data rather than derived.

Keshet & Naor (2014, 2016) developed an analytic approximation for the flow along the symmetry axis in front of an axisymmetric blunt body. Writing the flow velocity near the axis in cylindrical coordinates as v≃(−u, q(u)θ, 0)\mathbf{v}\simeq(-u,\,q(u)\theta,\,0) — a hodograph-like ansatz (a hodograph is a formulation in which the velocity components, rather than the spatial coordinates, serve as the independent variables) in which the perpendicular velocity gradient is expanded in powers of the parallel velocity uu — the steady Euler equations for a polytropic gas with adiabatic index γ\gamma reduce to a single integral relation. Here velocities are normalized so that u=1u=1 far upstream and u=0u=0 at the stagnation point on the body, and M0≡v/cˉM_0\equiv v/\bar c is the Mach number based on the stagnation sound speed, whose sonic (transonic) and strong-shock limiting values are W=2/(γ+1)W=\sqrt{2/(\gamma+1)} and S=2/(γ−1)S=\sqrt{2/(\gamma-1)}. In the supersonic regime, expanding the Euler equations about the point just downstream of the shock with the Rankine–Hugoniot jump conditions — the conservation laws relating the flow states on the two sides of a shock — as boundary conditions gives the expansion coefficients as explicit functions of a single unknown: the normalized shock curvature ξ−1≡(R/rs)θ=0\xi^{-1}\equiv(R/r_s)_{\theta=0}, the ratio of the shock's local radius of curvature RR to its nose distance rsr_s from the body center measured on the axis. Once ξ(M~0)\xi(\tilde M_0) is known (with M~0\tilde M_0 the far-upstream value of M0M_0), the entire axial flow, the standoff distance, and post-shock quantities such as the compression ratio gg follow algebraically; the approach and its astrophysical applications (planetary bow shocks, galaxy-cluster cold fronts, magnetically draped layers) are laid out in Keshet & Naor (2016).

The unsolved link in this program is the curvature input. Anchoring a second-order expansion of ξ\xi on the weak-shock constraint (smooth transonic transition) and the strong-shock limit (the shock curvature approaches the body's nose curvature, ξ→1\xi\to1 for a sphere) yields a parameter-free approximate formula for ξ(M~0)\xi(\tilde M_0), which works well but is not a derivation of the exact function. Alternatively, a one-parameter power law

ξ≃(M~0−WS−W)β\xi\simeq\left(\frac{\tilde M_0-W}{S-W}\right)^{\beta}

reproduces the measured shock curvature for a sphere over most of the supersonic range with β≃1/2\beta\simeq1/2 (best fits β=0.48\beta=0.48 for γ=7/5\gamma=7/5 and β=0.52\beta=0.52 for γ=5/3\gamma=5/3), but this exponent has no first-principles origin: it is fitted to classic wind-tunnel data. The same papers state that for an aspherical body, which modifies the weak- and strong-shock anchors themselves, higher-order terms or a re-tuned β\beta may be needed — leaving the body-dependence of the curvature law undetermined — and conjecture, without proof, that the whole stagnation-line analysis could be generalized to a non-axisymmetric body by replacing the symmetry axis with the stagnation streamline, provided the axial-velocity profile remains monotonic between the shock and the body and suitable perpendicular-gradient conditions close the problem. A follow-up within the program (Wallerstein & Keshet 2022) extended the subsonic Janzen–Rayleigh expansion — a perturbation series in the square of the upstream Mach number for potential flow — to arbitrary order, sharpened the sphere coefficients, and found an approximate universal scaling for prolate spheroids in subsonic flow, but left the supersonic curvature function, its exponent, and the non-axisymmetric extension untouched; adjacent communities still rely on empirical conic fits or simulation calibration (e.g. Mars bow-shock modeling, galaxy-cluster shocks).

Problem Statement

Derive, from the steady compressible Euler equations for an inviscid gas, the normalized nose curvature ξ(M~0)\xi(\tilde M_0) of the detached bow shock in front of a blunt body — the quantity that closes the Keshet–Naor stagnation-line expansion — and thereby place the empirical power-law fit ξ≃[(M~0−W)/(S−W)]β\xi\simeq[(\tilde M_0-W)/(S-W)]^{\beta} on a first-principles footing. Here M~0\tilde M_0 is the far-upstream value of the Mach number M0≡v/cˉM_0\equiv v/\bar c based on the stagnation sound speed cˉ\bar c, whose transonic and strong-shock limiting values are W=2/(γ+1)W=\sqrt{2/(\gamma+1)} and S=2/(γ−1)S=\sqrt{2/(\gamma-1)} for a polytropic gas of adiabatic index γ\gamma; ξ−1≡(R/rs)θ=0\xi^{-1}\equiv(R/r_s)_{\theta=0} is the ratio, evaluated on the symmetry axis, of the bow shock's local radius of curvature RR to its nose distance rsr_s from the body center; and β\beta is the exponent of the fit. Concretely, for the steady, axisymmetric supersonic regime: (i) for a sphere, derive the curvature function over the full supersonic range, and in particular the exponent β\beta and its dependence on the adiabatic index γ\gamma (fitted values are 0.48 for γ=7/5\gamma=7/5 and 0.52 for γ=5/3\gamma=5/3); (ii) determine how ξ(M~0)\xi(\tilde M_0) and the effective exponent are modified for an aspherical axisymmetric blunt body, whose nose shape changes the weak- and strong-shock limits themselves — deciding whether the correction appears as a re-tuned β\beta, higher-order terms, or a different functional form; and (iii) establish or refute the conjecture that the stagnation-line analysis generalizes to a non-axisymmetric body by replacing the symmetry axis with the stagnation streamline: formulate the conditions (monotonicity of the axial-velocity profile between shock and body, and the perpendicular-gradient constraints that close the expansion) under which the generalized expansion is well-posed and convergent, and prove or disprove that they hold. The curvature result should reproduce the measured and simulated shock curvature and standoff distance for the sphere and give testable predictions for aspherical and non-axisymmetric bodies.

References

  1. Keshet, U., & Naor, Y. (2016). Compressible flow in front of an axisymmetric blunt object: analytic approximation and astrophysical implications. The Astrophysical Journal 830(2), 147. DOI: 10.3847/0004-637X/830/2/147. https://arxiv.org/abs/1611.06234
  2. Keshet, U., & Naor, Y. (2014). Compressible flow in front of an axisymmetric blunt object. arXiv:1412.8406 (letter superseded by the 2016 ApJ paper). https://arxiv.org/abs/1412.8406
  3. Wallerstein, I. S., & Keshet, U. (2022). Compressible potential flows around round bodies: Janzen–Rayleigh expansion inferences. Journal of Fluid Mechanics 932, A6. DOI: 10.1017/jfm.2021.965. https://arxiv.org/abs/2011.13381
  4. Keshet, U., Raveh, I., & Naor, Y. (2021). Protruding bullet heads indicating dark matter pull. Monthly Notices of the Royal Astronomical Society 508(3), 3455–3462. DOI: 10.1093/mnras/stab2808. https://arxiv.org/abs/2101.06271
  5. Zhang, C., Churazov, E., Forman, W. R., & Jones, C. (2019). Standoff distance of bow shocks in galaxy clusters as proxy for Mach number. Monthly Notices of the Royal Astronomical Society 482(1), 20–29. DOI: 10.1093/mnras/sty2501. https://arxiv.org/abs/1808.02885
  6. Simon Wedlund, C., Volwerk, M., Beth, A., Mazelle, C., Möstl, C., Halekas, J., Gruesbeck, J., & Rojas-Castillo, D. (2022). A fast bow shock location predictor-estimator from 2D and 3D analytical models: Application to Mars and the MAVEN mission. Journal of Geophysical Research: Space Physics 127(1), e2021JA029942. DOI: 10.1029/2021JA029942. https://arxiv.org/abs/2109.04366
  7. Kleimann, J., Röken, C., & Fichtner, H. (2017). An improved analytical model of the local interstellar magnetic field: the extension to compressibility. The Astrophysical Journal 838(1), 75. DOI: 10.3847/1538-4357/aa5f15. https://arxiv.org/abs/1702.01619
  8. Amoloye, T. O., Oladimeji, L. T., Hayajnh, M. A., & Olayemi, O. A. (2026). Modeling the compressible flow field of an impulsively started circular cylinder with refined potential flow theory. Physics of Fluids 38(8), 086108. DOI: 10.1063/5.0332888. https://arxiv.org/abs/2608.19246
  9. Billig, F. S. (1967). Shock-wave shapes around spherical- and cylindrical-nosed bodies. Journal of Spacecraft and Rockets 4(6), 822–823. DOI: 10.2514/3.28969. https://doi.org/10.2514/3.28969