Oblique Shock Calculator

Pri Geens

Pri Geens

Oblique Shock Calculator

Weak Shock Solution Properties

Wave Angle (β)
Downstream Mach Number (M₂)
Static Pressure Ratio (p₂ / p₁)
Static Temperature Ratio (T₂ / T₁)
Density Ratio (ρ₂ / ρ₁)
Stagnation Pressure Ratio (p₀₂ / p₀₁)
Aerodynamic Analysis
Calculations utilize standard θ-β-M relations and normal shock tables aligned with NACA 1135 and modern compressible flow aerodynamic principles. Results assume a weak shock solution, calorically perfect gas, and steady adiabatic flow with no external work.

What Is an Oblique Shock Calculator?

An Oblique Shock Calculator solves the relationship between upstream Mach number, flow deflection angle, shock angle, and gas properties. It models a supersonic flow that turns through a compression corner and forms an angled shock wave.

The calculator uses the entered Mach number, deflection angle, and specific heat ratio to find the attached weak shock solution. It then applies normal-shock equations to the velocity component perpendicular to the wave. The result shows how Mach number, static pressure, temperature, density, and stagnation pressure change across the shock.

This tool is designed for quick calculations and technical study. It can identify a Mach wave at zero deflection, calculate an attached weak oblique shock, or report that the shock is detached when no attached solution is found.

How the Oblique Shock Calculator Formula Works

The calculator first computes the Mach angle, which is the lowest possible wave angle for the entered upstream Mach number.

μ=sin1(1M1)\mu = \sin^{-1}\left(\frac{1}{M_1}\right)

It then searches between the Mach angle and 90 degrees for roots of the coded theta-beta-M relationship:

θ=atan2(2cosβ(M12sin2β1),sinβ[M12(γ+cos2β)+2])\theta = \operatorname{atan2}\left(2\cos\beta\left(M_1^2\sin^2\beta-1\right),\;\sin\beta\left[M_1^2\left(\gamma+\cos 2\beta\right)+2\right]\right)

Here, M₁ is upstream Mach number, θ is flow deflection angle, β is wave angle, and γ is the specific heat ratio. The first root found is used as the weak shock solution. A second root, when found, is reported only in the aerodynamic notes as the strong solution.

M1n=M1sinβM_{1n}=M_1\sin\beta
M2n2=1+γ12M1n2γM1n2γ12M_{2n}^2=\frac{1+\frac{\gamma-1}{2}M_{1n}^2}{\gamma M_{1n}^2-\frac{\gamma-1}{2}}
M2=M2nsin(βθ)M_2=\frac{M_{2n}}{\sin\left(\beta-\theta\right)}

The static property ratios come from the normal component of upstream Mach number:

p2p1=1+2γγ+1(M1n21)\frac{p_2}{p_1}=1+\frac{2\gamma}{\gamma+1}\left(M_{1n}^2-1\right)
ρ2ρ1=(γ+1)M1n22+(γ1)M1n2\frac{\rho_2}{\rho_1}=\frac{\left(\gamma+1\right)M_{1n}^2}{2+\left(\gamma-1\right)M_{1n}^2}
T2T1=p2/p1ρ2/ρ1\frac{T_2}{T_1}=\frac{p_2/p_1}{\rho_2/\rho_1}
p02p01=[(γ+1)M1n22+(γ1)M1n2]γγ1[2γM1n2(γ1)γ+1]1γ1\frac{p_{02}}{p_{01}}=\left[\frac{\left(\gamma+1\right)M_{1n}^2}{2+\left(\gamma-1\right)M_{1n}^2}\right]^{\frac{\gamma}{\gamma-1}}\left[\frac{2\gamma M_{1n}^2-\left(\gamma-1\right)}{\gamma+1}\right]^{-\frac{1}{\gamma-1}}

For example, use air with γ = 1.4, M₁ = 2.5, and θ = 15 degrees. The calculator finds β = 36.94 degrees and M₁n = about 1.503. It displays M₂ = 1.874, p₂/p₁ = 2.4675, ρ₂/ρ₁ = 1.8665, T₂/T₁ = 1.3220, and p₀₂/p₀₁ = 0.9290.

How to Use the Oblique Shock Calculator: Step by Step

  1. Select a Gas Type. Available choices are air or another diatomic gas, carbon dioxide, helium or another monatomic gas, methane, and a custom specific heat ratio.
  2. If you select Custom γ, enter a specific heat ratio strictly greater than 1 and strictly less than 2.
  3. Enter the Upstream Mach Number (M₁). The calculation requires a value strictly greater than 1 because the modeled incoming flow must be supersonic.
  4. Enter the Deflection Angle (θ) in degrees. Use zero for an unturned flow and a positive value for a compression turn.
  5. Select Calculate. The tool searches for the weak attached shock solution and displays the calculated flow properties.
  6. Select Reset to restore air as the selected gas, clear the numeric fields, hide the custom gamma input, and remove the displayed results.

The wave angle is shown in degrees. Downstream Mach number is displayed to three decimal places. Static pressure, temperature, density, and stagnation pressure ratios are displayed to four decimal places. Each ratio compares the downstream value with the upstream value. For example, p₂/p₁ greater than 1 means static pressure increased across the shock.

How to Read Oblique Shock Calculator Results

Wave angle and downstream Mach number

The wave angle β measures the angle between the shock wave and the upstream flow direction. For an attached weak shock, β is greater than the Mach angle and less than 90 degrees. The downstream Mach number shows how much the flow slows after crossing the wave.

Static and stagnation property ratios

OutputMeaning
p₂/p₁Downstream static pressure divided by upstream static pressure
T₂/T₁Downstream static temperature divided by upstream static temperature
ρ₂/ρ₁Downstream density divided by upstream density
p₀₂/p₀₁Downstream stagnation pressure divided by upstream stagnation pressure

A stagnation pressure ratio below 1 represents total pressure loss. When that ratio is below 0.95, the calculator adds a note showing the loss as a percentage. That percentage is calculated as one minus the stagnation pressure ratio, multiplied by 100.

Zero deflection and detached shocks

At zero deflection, the code reports the Mach angle instead of a finite-strength shock. It keeps M₂ equal to M₁ and sets every property ratio to 1.0000. This represents unchanged properties across an infinitesimally weak Mach wave.

If the root search finds no attached solution, the wave angle displays as “Detached.” All other numerical outputs display N/A. The tool does not calculate the curved bow-shock structure or its downstream property field.

Assumptions and limitations

The calculation assumes steady, adiabatic flow with no external work and a calorically perfect gas. This means γ stays constant. Real gases can have temperature-dependent properties. The tool also focuses on the weak solution and does not model viscosity, heat transfer, three-dimensional effects, boundary layers, shock interactions, or chemical reactions.

The input field displays a deflection range up to 90 degrees, but the calculation script directly rejects only negative angles. Use physically meaningful compression angles. The numerical solver searches in fixed angular increments before applying bisection, so results near a solution limit may depend on the root search finding a sign change.

Frequently Asked Questions

What is an oblique shock wave?

An oblique shock wave is an angled compression wave in supersonic flow. It can form when the flow turns into itself at a wedge, ramp, inlet, or leading edge. Across the wave, static pressure, temperature, and density rise, while Mach number and stagnation pressure generally decrease.

How do I calculate the oblique shock angle?

Enter upstream Mach number, deflection angle, and the gas specific heat ratio. The calculator solves its theta-beta-M equation numerically between the Mach angle and 90 degrees. It uses the first root as the weak shock angle and may identify a second root as the strong solution.

What is the difference between weak and strong oblique shocks?

The weak solution has the smaller wave angle and is the solution used for the displayed property calculations. The strong solution has a larger wave angle. When the code detects that second root, it mentions the angle in the aerodynamic analysis but does not display separate strong-shock property ratios.

Why does the calculator say the shock is detached?

The calculator reports a detached shock when its numerical search cannot find an attached weak solution for the entered Mach number, deflection angle, and γ. It then displays N/A for downstream properties because the coded attached-shock equations do not model the curved and non-uniform flow behind a detached bow shock.

What specific heat ratio should I use?

Use the gas option that matches your problem. The calculator assigns γ = 1.400 for air or a diatomic gas, 1.289 for carbon dioxide, 1.667 for helium or a monatomic gas, and 1.299 for methane. A custom value must be strictly between 1 and 2.

What happens when the deflection angle is zero?

At zero deflection, the calculator reports an infinitesimally weak Mach wave rather than a finite shock. The wave angle equals arcsine of 1 divided by M₁. Downstream Mach number remains equal to upstream Mach number, and all static and stagnation property ratios remain exactly 1.0000.

How accurate is the Oblique Shock Calculator?

The calculator follows the equations and numerical search coded into the tool. Its results are suited to ideal compressible-flow calculations under the stated assumptions. Real aerodynamic results may differ because of variable gas properties, viscosity, boundary layers, geometry, heat transfer, shock interactions, measurement error, and other effects not included here.