Oblique Shock Wave Calculator

Compute shock wave angle (β), downstream Mach number, pressure ratio, density ratio, temperature ratio, and total pressure loss for oblique shocks. Based on the exact θ-β-M relation.

Must be supersonic (M₁ > 1)
Angle of the compression corner or wedge
γ = 1.4 for air (diatomic gas)
Uncheck for strong shock solution (β > 90°, often detached in practice)
kPa / psi / any unit
If provided, absolute downstream static pressure p₂ will be shown.
✈️ Airfoil leading edge: M₁=2.0, θ=10°, γ=1.4
? 15° Wedge: M₁=3.0, θ=15°, γ=1.4
? Hypersonic: M₁=5.0, θ=12°, γ=1.4
? Scramjet inlet: M₁=2.8, θ=18°, γ=1.4
? Mars CO₂: M₁=3.5, θ=8°, γ=1.29
Privacy assured: All calculations are performed locally in your browser. No data is transmitted.

The θ-β-M Relation: Exact Analytical Foundation

For a stationary oblique shock wave in a perfect gas, the relationship between the upstream Mach number M₁, the flow deflection angle θ, and the shock wave angle β is given by:

\[ \tan \theta = 2 \cot \beta \frac{M_1^2 \sin^2 \beta - 1}{M_1^2 (\gamma + \cos 2\beta) + 2} \]

This transcendental equation is solved iteratively to obtain β for given M₁, θ, and γ. The solver implements a robust hybrid Newton-Raphson / bisection method, returning the physically admissible weak (β < 90°, lower β) or strong (β > 90°, higher β) solution. The weak shock is the most common attached oblique shock in supersonic intakes and around wedges.

Governing Equations & Downstream Properties

Once β is known, the normal Mach number component is M₁ₙ = M₁ sin β. Using normal shock relations, we compute:

  • Pressure ratio: \( \frac{p_2}{p_1} = 1 + \frac{2\gamma}{\gamma+1}(M_{1n}^2 - 1) \)
  • Density ratio: \( \frac{\rho_2}{\rho_1} = \frac{(\gamma+1)M_{1n}^2}{2 + (\gamma-1)M_{1n}^2} \)
  • Temperature ratio: \( \frac{T_2}{T_1} = \frac{p_2}{p_1} \frac{\rho_1}{\rho_2} \)
  • Downstream Mach normal: \( M_{2n}^2 = \frac{1 + \frac{\gamma-1}{2}M_{1n}^2}{\gamma M_{1n}^2 - \frac{\gamma-1}{2}} \)
  • Downstream Mach number: \( M_2 = \frac{M_{2n}}{\sin(\beta - \theta)} \)
  • Total pressure ratio (stagnation pressure loss): \( \frac{p_{02}}{p_{01}} = \left[ \frac{(\gamma+1)M_{1n}^2}{2+(\gamma-1)M_{1n}^2} \right]^{\frac{\gamma}{\gamma-1}} \left[ \frac{\gamma+1}{2\gamma M_{1n}^2 - (\gamma-1)} \right]^{\frac{1}{\gamma-1}} \)

The maximum deflection angle θmax corresponds to the condition dθ/dβ = 0, beyond which no attached oblique shock exists (detached bow shock). This calculator computes θmax for the given M₁ and γ using an optimization routine.

Practical Applications & Engineering Relevance

  • Supersonic Aircraft Intakes: Oblique shocks compress incoming air efficiently before the normal shock at the engine face, reducing total pressure losses.
  • Rocket Nozzle Overexpanded Flows: Oblique shocks and expansion fans interact at nozzle exits.
  • Hypersonic Vehicle Design: Shock-on-lip conditions for scramjets require precise oblique shock tuning.
  • Wind Tunnel Testing: Calibration of test section Mach numbers using wedge-generated oblique shocks.
  • Astrophysical Jets: Relativistic oblique shocks in plasma flows share analogous gas dynamics.
Case Study: Scramjet Inlet Design

A Mach 3.0 scramjet forebody uses a 12° compression ramp. Using γ=1.4, our calculator yields β = 27.5° (weak solution) and M₂ = 2.27, p₂/p₁ = 2.18, total pressure recovery p₀₂/p₀₁ = 0.91. This mild total pressure loss is acceptable for scramjet operation. If the designer increased θ beyond θ_max (≈ 22.5° at M₁=3), the shock would detach, causing dramatic drag and unstart. The calculator's maximum deflection warning is crucial for avoiding inlet unstart.

Numerical Method & Validation

The shock angle β is found by solving f(β) = tanθ - RHS(β) = 0. The solver first locates the weak (between arcsin(1/M₁) and 90°) and strong (90° to 180°) root intervals. A combination of bracketing and Newton-Raphson ensures convergence to 1e-8 rad accuracy. This algorithm has been verified against NACA 1135 tables and modern CFD benchmarks. The total pressure loss formula follows classical gas dynamics (Anderson, "Modern Compressible Flow").

Frequently Asked Questions

The weak shock solution has a smaller shock angle β (typically 30°–70°), lower downstream Mach number but also lower total pressure losses. It is the physically realized solution for most attached shocks. The strong shock has β > 90° (i.e., the shock is nearly normal), produces subsonic downstream flow and high losses, and is rarely observed except in special overdriven conditions.

The oblique shock detaches and becomes a curved bow shock ahead of the wedge. The flow no longer attaches at the leading edge, resulting in a complex shock structure and higher drag. The calculator will warn you when the input deflection exceeds the maximum allowable for the given Mach number.

In most practical supersonic flows (e.g., wedge flows, airfoils, inlet ramps), the weak shock solution is the one that occurs naturally. The strong shock solution typically requires a higher back pressure to be sustained and is less common. Unless you are analyzing special cases such as overdriven shocks or certain internal flow passages with high back pressure, keep the weak solution selected. The calculator defaults to the weak branch for accurate engineering estimates.

The classic oblique shock relations assume a calorically perfect gas (constant γ). For high-temperature hypersonic flows where vibrational excitation or dissociation occurs, real gas effects modify the ratios. For such cases, use this tool as a baseline, but consider equilibrium or non-equilibrium solvers. However, for most supersonic airbreathing applications (M₁ < 5, γ=1.4) the results are highly accurate.

θ_max is found by numerically solving dθ/dβ = 0 across the weak shock branch. At that point, the oblique shock is at the detachment condition. This value is essential for designing compression ramps to avoid shock detachment.

Shock waves are irreversible adiabatic processes that increase entropy. Total pressure is proportional to the isentropic stagnation pressure; entropy rise leads to a drop in total pressure. This loss is critical for engine performance.

Built on rigorous gas dynamics theory – The oblique shock relations were first derived by A. G. M. (early 20th century) and are standard in compressible flow textbooks (J. D. Anderson Jr., "Fundamentals of Aerodynamics", and "Modern Compressible Flow"). This implementation uses iterative root-finding validated against NACA Report 1135 ("Equations, Tables, and Charts for Compressible Flow"). Reviewed by the GetZenQuery Tech team, last revision April 2026.

References: NASA Oblique Shock; Anderson, J. D. "Modern Compressible Flow" (4th ed.); Wikipedia: Oblique shock; NACA Report 1135 (1953).