Oblique Shock Wave Calculator (θ-β-Mach Relation)

Solve the supersonic compressible θ-β-M relation, weak and strong shock angles, downstream Mach numbers, pressure jumps, and detachment boundaries.

📁 Physics🔒 100% Free & Private⚡ Instant Output
Quick Overview & Direct Answer

Oblique Shock Wave Calculator (θ-β-M Relation): Free Oblique Shock Wave Calculator tool: calculate supersonic oblique shock wave angles, downstream Mach number, pressure ratio, temperature jump.

Primary Use Physical formula evaluation, mechanics, thermodynamics calculations
Input Formats Physical values, mass, velocity, acceleration, energy variables
Output Format Standardized SI units, converted physics quantities
Processing Mode 100% Client-Side Mathematical Engine
Access Tier Free, unlimited, no sign-up needed

Upstream Flow Parameters

Supersonic freestream ($M_1 > 1.0$)
Flow deflection in degrees (°)
Air = 1.40, Hot gas = 1.33, Mono = 1.67

Freestream Thermodynamic State

kPa (absolute)
Kelvin (K)

Downstream Aerodynamic State

36.95°
Shock Wave Angle ($beta$)
1.87
Downstream Mach ($M_2$)
2.47
Pressure Ratio ($p_2 / p_1$)
0.963
Total Pressure Ratio ($p_{02} / p_{01}$)
Max Deflection ($theta_{max}$): 29.80°
Mach Angle ($mu$): 23.58°
Normal Mach ($M_{n1} to M_{n2}$): 1.503 → 0.701
Temperature Ratio ($T_2 / T_1$): 1.320 (380.4 K)
Density Ratio ($rho_2 / rho_1$): 1.871
Shock Regime: Attached Oblique Shock

Compressible Aerodynamics & Oblique Shock Wave Physics

In supersonic flow ($M > 1$), fluid particles travel faster than the speed at which acoustic pressure disturbances can propagate upstream. When a supersonic stream encounters a compressive geometrical deflection—such as a sharp wedge, the leading edge of a supersonic diamond airfoil, or an internal engine inlet ramp—the fluid cannot adjust gradually through smooth subsonic streamlines. Instead, it undergoes an abrupt, virtually instantaneous discontinuity in velocity, pressure, temperature, and entropy termed an Oblique Shock Wave.

Unlike a normal shock wave which stands perpendicular to the upstream velocity vector ($beta = 90^circ$), an oblique shock wave is inclined at an angle $beta$ relative to the incoming freestream. As the flow traverses the shock sheet, it is deflected toward the wedge surface by angle $theta$.

1. The Fundamental $theta$-$beta$-$text{Mach}$ Relation

By decomposing upstream velocity $V_1$ into components normal ($V_{n1} = V_1 sin beta$) and tangential ($V_{t1} = V_1 cos beta$) to the shock wave, and applying the conservation laws of mass, normal momentum, tangential momentum, and energy across the control volume, one derives the famous $theta$-$beta$-$M$ relation:

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

Where:

  • $M_1$: Freestream Mach number upstream of the shock ($M_1 ge 1.0$).
  • $theta$: Wedge deflection angle of the flow ($0^circ le theta le theta_{max}$).
  • $beta$: Oblique shock wave angle relative to the upstream flow direction.
  • $gamma$: Ratio of specific heats ($c_p / c_v$, dimensionless; $1.40$ for diatomic air).

2. Weak vs. Strong Shock Solutions and Shock Detachment

For any supersonic Mach number $M_1$ and a given deflection angle $theta$ below a critical upper limit $theta_{max}$, the $theta$-$beta$-$M$ relation yields two mathematically valid positive roots for $beta$:

  1. Weak Shock Solution ($beta_{text{weak}}$): Characterized by a lower shock wave angle ($beta < beta_{text{sonic}}$). The flow deflection occurs with moderate compression and minimal stagnation pressure loss. Crucially, the downstream Mach number $M_2$ remains supersonic ($M_2 > 1$) across almost the entire operational range (except for a microscopic sliver immediately adjacent to $theta_{max}$). Because the downstream pressure does not feed upstream perturbations, the weak shock is the physical solution observed on isolated wedges in unconfined free flow.
  2. Strong Shock Solution ($beta_{text{strong}}$): Characterized by a much steeper wave angle ($beta > beta_{text{sonic}}$, approaching $90^circ$). The normal Mach component $M_{n1}$ is large, producing severe compression, high thermal jumps, and substantial entropy increase ($p_{02} ll p_{01}$). The downstream flow is strictly subsonic ($M_2 < 1$). The strong shock can only exist when forced by external downstream boundary conditions, such as high backpressure in an engine diffuser or wind tunnel diffuser throttle.

As the deflection angle $theta$ is incrementally increased, the weak and strong roots converge toward each other until they merge at the maximum deflection angle ($theta_{max}$):

sin^2 beta_{theta_{max}} = frac{1}{gamma M_1^2} left[ frac{(gamma + 1) M_1^2 - 4}{4} + sqrt{(gamma + 1)left( 1 + frac{gamma - 1}{2} M_1^2 + frac{gamma + 1}{16} M_1^4 right)} right]

If the physical wedge angle exceeds $theta_{max}$, a straight oblique shock wave is physically impossible. The shock detaches from the wedge tip, moving upstream to form a detached, curved bow shock with a subsonic stagnation bubble between the shock apex and the blunt body.

3. Exact Thermodynamic Jump Relations

Because tangential velocity is perfectly preserved across an inviscid oblique shock ($V_{t2} = V_{t1}$), all changes in thermodynamic state variables are governed strictly by the normal component of the upstream Mach number:

M_{n1} = M_1 sin beta

The normal component of the downstream Mach number $M_{n2}$ satisfies standard normal shock relations:

M_{n2}^2 = frac{2 + (gamma - 1) M_{n1}^2}{2gamma M_{n1}^2 - (gamma - 1)}

The actual downstream Mach number $M_2$ is recovered via trigonometry:

M_2 = frac{M_{n2}}{sin(beta - theta)}

Static pressure, density, temperature, and total pressure jump ratios are calculated as:

frac{p_2}{p_1} = 1 + frac{2gamma}{gamma + 1}(M_{n1}^2 - 1) = frac{2gamma M_{n1}^2 - (gamma - 1)}{gamma + 1}
frac{rho_2}{rho_1} = frac{(gamma + 1) M_{n1}^2}{2 + (gamma - 1) M_{n1}^2}
frac{T_2}{T_1} = left(frac{p_2}{p_1}right) cdot left(frac{rho_1}{rho_2}right) = left[ frac{2gamma M_{n1}^2 - (gamma - 1)}{gamma + 1} right] left[ frac{2 + (gamma - 1) M_{n1}^2}{(gamma + 1) M_{n1}^2} right]
frac{p_{02}}{p_{01}} = left[ frac{(gamma + 1) M_{n1}^2}{2 + (gamma - 1) M_{n1}^2} right]^{frac{gamma}{gamma - 1}} left[ frac{gamma + 1}{2gamma M_{n1}^2 - (gamma - 1)} right]^{frac{1}{gamma - 1}}

Because the shock process is adiabatic without external work or heat transfer, the stagnation (total) temperature is conserved: $T_{02} = T_{01}$.

Step-by-Step Worked Aerodynamic Example

Consider a supersonic flight vehicle traveling at $M_1 = 2.50$ through standard sea-level air ($gamma = 1.40$, $p_1 = 101.325 text{kPa}$, $T_1 = 288.15 text{K}$). The nose cone possesses a half-wedge angle of $theta = 15.0^circ$.

Step 1: Check for Shock Attachment ($theta_{max}$)

Evaluating $beta_{theta_{max}}$ for $M_1 = 2.50$ and $gamma = 1.40$ yields $theta_{max} approx 29.80^circ$. Because $theta = 15.0^circ < 29.80^circ$, an attached straight oblique shock exists!

Step 2: Solve the $theta$-$beta$-$text{Mach}$ Equation for $beta$

Using the weak shock root:

tan(15.0^circ) = 0.26795

Iterating the $theta$-$beta$-$M$ equation yields:

beta_{text{weak}} = 36.95^circ quad (text{Note: } beta_{text{strong}} = 84.62^circ)

Step 3: Evaluate Normal Mach Components

M_{n1} = M_1 sin beta = 2.50 times sin(36.95^circ) = 2.50 times 0.6011 = 1.5028
M_{n2}^2 = frac{2 + 0.4 times (1.5028)^2}{2.8 times (1.5028)^2 - 0.4} = frac{2 + 0.9034}{6.323 - 0.4} = frac{2.9034}{5.923} approx 0.4902 implies M_{n2} = 0.7001

Step 4: Compute Downstream Mach Number $M_2$

beta - theta = 36.95^circ - 15.0^circ = 21.95^circ
M_2 = frac{M_{n2}}{sin(beta - theta)} = frac{0.7001}{sin(21.95^circ)} = frac{0.7001}{0.3738} approx 1.873

Notice that $M_2 = 1.873 > 1.0$; the downstream flow remains supersonic, as expected for a weak shock solution.

Step 5: Compute Thermodynamic Jumps

frac{p_2}{p_1} = 1 + frac{2 times 1.4}{2.4} (1.5028^2 - 1) = 1 + 1.1667 times (2.2584 - 1) = 1 + 1.468 = 2.468
p_2 = 101.325 times 2.468 = 250.1 text{kPa}
frac{rho_2}{rho_1} = frac{2.4 times 2.2584}{2 + 0.4 times 2.2584} = frac{5.420}{2.9034} approx 1.867
frac{T_2}{T_1} = frac{p_2 / p_1}{rho_2 / rho_1} = frac{2.468}{1.867} approx 1.322 implies T_2 = 288.15 times 1.322 = 381.0 text{K} (+107.8^circtext{C})
frac{p_{02}}{p_{01}} approx 0.9298 (7.02% text{loss in total pressure})

Comparison of Oblique vs. Normal Shock Waves

Aerodynamic Characteristic Oblique Shock (Weak Solution) Oblique Shock (Strong Solution) Normal Shock ($beta = 90^circ$)
Wave Angle ($beta$) Mach angle $mu le beta < beta_{text{sonic}}$ $beta_{text{sonic}} < beta < 90^circ$ Strictly $beta = 90.0^circ$
Downstream Mach ($M_2$) Supersonic ($M_2 > 1$) Subsonic ($M_2 < 1$) Subsonic ($M_2 < 1$)
Stagnation Pressure Loss Small to moderate ($p_{02}/p_{01} approx 0.90 - 0.99$) Severe ($p_{02}/p_{01} approx 0.30 - 0.70$) Maximum possible loss for given $M_1$
Physical Occurrence Sharp wedges, supersonic wings, aircraft noses High-backpressure internal ducts, wind tunnels Blunt body stagnation line, convergent-divergent nozzles

Frequently Asked Questions

The θ-β-M relation is the governing implicit equation relating the upstream supersonic Mach number (M1), the flow deflection or wedge angle (θ), and the resulting oblique shock wave angle (β). Mathematically: tan(θ) = 2 * cot(β) * (M1² * sin²(β) - 1) / (M1² * (γ + cos(2β)) + 2).

For any given upstream Mach number and deflection angle below the detachment limit (θ < θ_max), the θ-β-M equation yields two distinct physical shock angles. The 'weak shock' has a smaller wave angle (β), lower pressure jump, and the downstream flow remains supersonic (M2 > 1, except for a tiny zone near θ_max). This is the naturally occurring shock observed on unconfined supersonic wedges. The 'strong shock' has a steeper angle (closer to a normal shock at 90°), much higher entropy increase, and downstream flow is strictly subsonic (M2 < 1), requiring high downstream backpressure to sustain.

When the flow deflection angle exceeds the maximum theoretical deflection angle θ_max, straight oblique shock solutions cease to exist. The shock detaches from the wedge apex and stands upstream as a curved detached bow shock. Along the stagnation streamline directly ahead of the nose, the bow shock is normal (subsonic downstream), transitioning smoothly into a weak oblique shock and eventually an acoustic Mach wave far away.

Because no heat is added to or extracted from the gas as it traverses the shock layer, the flow is strictly adiabatic, meaning total enthalpy and stagnation temperature remain constant (T02 = T01). However, the shock is an irreversible thermodynamic process driven by viscous dissipation and thermal conduction inside the thin shock structure, generating entropy (Δs > 0). This entropy generation causes a mandatory drop in stagnation total pressure (p02 < p01).