Normal Shock Wave Relations lab report

1. Title and aim

Normal Shock Wave Pressure and Total Pressure Loss in Supersonic Ducts. Measure static pressure jump P_2 / P_1 and total pressure ratio P_02 / P_01 across a stationary normal shock wave in a supersonic duct as upstream Mach M_1 varies.

2. Abstract

Six controlled aerodynamic test points were executed across standardized inputs M_1^2 - 1 = 0.8 dimensionless to M_1^2 - 1 = 4.29 dimensionless for Normal Shock Wave Pressure and Total Pressure Loss in Supersonic Ducts. An ordinary least-squares linear regression yielded an empirical Rankine-Hugoniot Slope 2*gamma/(gamma+1) of 1.1667 dimensionless (difference from theoretical reference = 0.0023%). The empirical telemetry conclusively validates the fundamental aerospace engineering formulation with R² = 1.000000. The investigation demonstrates stable flow behavior and mechanical response across all test conditions with negligible residual variance. Data recorded through the virtual wind tunnel and telemetry bus exhibit strict adherence to established boundary layer, shock wave, flight mechanics, or structural mechanics governing equations, confirming the dimensional and similitude assertions of the curriculum.

3. Apparatus and fixed settings

• Supersonic Wind Tunnel Test Section: Continuously variable supersonic duct operating from Mach M = 1.2 to 2.5.
• Variable Back-Pressure Exit Throttle Valve: Regulates downstream back-pressure to position a stable normal shock wave in the test section.
• Multi-Channel Fast-Response Pressure Scanner: Samples wall static pressures upstream and downstream of the shock wave at 1 kHz.
• Supersonic Pitot-Static Total Pressure Rake: Measures total pressure upstream P_01 and downstream P_02 across the shock discontinuity.
• High-Speed Schlieren Imaging Camera: Captures optical shadowgraph and Schlieren images of the normal shock wave thickness and location.
• Digital Environmental DAQ Computer: Logs stagnation conditions and computes experimental Rankine-Hugoniot pressure ratios.
• Test gas: Air (gamma = 1.400)
• Nozzle design Mach range: 1.20 to 2.50
• Stagnation pressure P_01: 300.0 kPa
• Stagnation temperature T_01: 293.15 K

4. Theory and graph choice

A normal shock wave is an extremely thin discontinuous compression front (thickness on the order of the molecular mean free path, ~10^-7 m) oriented perpendicular to oncoming supersonic flow. As supersonic gas (Mach M_1 > 1) traverses this microscopic zone, intense viscous dissipation and thermal conduction abruptly decelerate the flow to subsonic speed (M_2 < 1), causing violent irreversible thermodynamic changes. By applying the fundamental conservation equations of mass, linear momentum, and energy across the normal shock wave, the classical Rankine-Hugoniot jump relations are derived. For a calorically perfect gas with specific heat ratio gamma = 1.4, the static pressure ratio across the shock scales strictly with the upstream Mach number squared according to P_2 / P_1 = 1 + [2*gamma / (gamma+1)] * (M_1^2 - 1) = (7/6)*M_1^2 - 1/6. Crucially, the entropy generation across the irreversible shock causes a severe loss of total stagnation pressure (P_02 < P_01), while total temperature remains constant (T_02 = T_01). In this laboratory experiment, a supersonic wind tunnel test section equipped with a variable back-pressure throttle positions a stationary normal shock over a static pressure tap array. Upstream Pitot-static measurements determine M_1, while downstream taps record the static pressure jump and stagnation pressure deficit across Mach numbers from M_1 = 1.34 to 2.30, validating the Rankine-Hugoniot gradient.

The governing physical relationship was linearized into the standard regression form y = mx + b. Here, the independent coordinate x represents Mach Factor (M_1^2 - 1), and the dependent coordinate y represents Static Pressure Jump Ratio (P_2/P_1 - 1). The slope m directly extracts Rankine-Hugoniot Slope 2*gamma/(gamma+1), while the intercept b accounts for zero-lift offsets, tare baselines, or residual reference shifts. Linear transformation allows robust parameter estimation via ordinary least squares without solver instability or local minima traps. Plotting transformed observation pairs highlights flow similitude coefficients, lift curve slopes, pressure ratios, or stiffness parameters directly as linear gradients. This linear graphical formulation provides an immediate visual check of homoscedasticity, linearity, and potential flow separation or structural non-linearities.

5. Method

1. Start the supersonic wind tunnel, establishing stable supersonic flow at nominal Mach M_1 = 1.34 in the test section.
2. Gradually close the downstream throttle valve until Schlieren optics confirm a stationary normal shock centered over tap station 12.
3. Record upstream static pressure P_1, downstream static pressure P_2, and stagnation pressures P_01 and P_02.
4. Adjust the supersonic nozzle block to step upstream Mach M_1 through values corresponding to (M_1^2 - 1) = 0.80, 1.25, 1.89, 2.61, 3.41, and 4.29.
5. Compute static pressure jump ratio (P_2/P_1 - 1) and total pressure ratio P_02 / P_01 at each Mach operating point.
6. Plot (P_2/P_1 - 1) versus (M_1^2 - 1), perform linear regression to determine slope 2*gamma/(gamma+1), and verify gamma = 1.4 for air.

6. Observations and calculation table

All aerodynamic observations use strictly controlled 64-bit IEEE 754 precision buffers. Residuals are computed relative to the ordinary least-squares regression line.

Download the sample observations (CSV) ↓

7. Worked calculation for one observation

8. Graph, gradient and final calculation

Circular markers represent empirical aerodynamic observations. The solid line represents the ordinary least-squares linear fit.

9. Uncertainty and comparison

Percentage difference = |estimate − theoretical value| / |theoretical value| × 100 = 0.0023%.

For n = 6 observations, Σresidual² = 1.316e-09; s² = Σresidual²/(n − 2). Gradient standard error = √[s²/Σ(x − x̄)²] = 6.122e-06. Intercept standard error = √[s²(1/n + x̄²/Σ(x − x̄)²)] = 1.632e-05.

Propagating the fit parameters yields a fit standard error of 6.122e-06 dimensionless, with an ideal coefficient of determination R² = 1.000000.

In experimental aeronautics and flight vehicle testing, empirical uncertainties originate from dynamic pressure transducer calibration drift, wind tunnel boundary layer blockage corrections, tare balance interaction tare loads, free-stream turbulence intensity, and numerical discretization (64-bit IEEE 754 precision). The ordinary least-squares regression across n = 6 standardized test points yields unbiased standard errors for both gradient and intercept. Because the governing aerodynamic and structural equations adhere strictly to conservation of mass, momentum, and energy, residual deviations are tightly bounded, yielding a coefficient of determination R² ≥ 0.999999 and confirming the theoretical scaling model. Furthermore, systematic errors—including test section solid blockage, wake blockage, buoyancy pressure gradients along the test section axis, and pitch balance sting flexure deflections—are corrected in accordance with standard aerodynamic testing protocols. Standard deviations are propagated using first-order Taylor series expansion (delta method) across transformed dimensionless parameters. Residual diagnostics demonstrate homoscedasticity across the entire operational envelope without systematic bias.

10. Result and discussion

The empirical wind tunnel and flight test observations provide rigorous quantitative verification of Normal Shock Wave Pressure and Total Pressure Loss in Supersonic Ducts. In aeronautical engineering, gas dynamics, and aerospace structures, physical systems are governed by Navier-Stokes conservation laws, isentropic compressible relations, rigid-body flight dynamics equations, and continuum elasticity. Across all six observations, the empirical slope m = 1.1667 confirms the theoretical model with near-zero residual variance. Furthermore, comparison with classical aeronautical benchmarks underscores critical design principles including aerodynamic efficiency, static stability margins, control authority, thermodynamic cycle limits, and structural mass optimization. Residual analysis reveals homoscedastic dispersion without systematic non-linear patterning, validating the underlying model assumptions. Confidence intervals calculated via asymptotic standard errors align with theoretical quantiles, indicating that finite-sample distortions remain negligible under controlled test conditions. The experimental protocol confirms the robustness of the methodology across varying flight conditions.

11. Precautions and limitations

• Maintain strict throttle valve stability to prevent high-amplitude normal shock oscillation (shock buzz) inside the test section.
• Verify that wall boundary layer separation does not induce a complex lambda-shock bifurcation near the static pressure taps.
• Ensure the supersonic Pitot tube is aligned with the tunnel axis to within +/- 0.5 degrees.
• Apply the Rayleigh supersonic Pitot formula when evaluating total pressure behind the detached bow shock of the probe.

Shock-wave boundary layer interaction (SWBLI) creates a lambda-foot near duct walls that slightly smears the idealized 1D jump. Linear regression of the pressure jump ratio against (M_1^2 - 1) yields a slope of 1.1667 matching the Rankine-Hugoniot relation with R^2 > 0.999999.

12. Conclusion

This laboratory investigation quantitatively verified the mathematical formulation and physical behavior of Normal Shock Wave Pressure and Total Pressure Loss in Supersonic Ducts. Through rigorous tracking across six observation points, we confirmed an empirical Rankine-Hugoniot Slope 2*gamma/(gamma+1) of 1.1667 dimensionless matching the theoretical reference 1.1667 dimensionless within 0.0023% error. The investigation demonstrates why this aerospace method serves as a fundamental cornerstone in aircraft and spacecraft design. The verified scaling laws and parameter estimates demonstrate the practical efficacy of aerodynamic modeling and aerospace engineering principles.