Schlieren Optical Imaging lab report

1. Title and aim

Supersonic Wind Tunnel Calibration and Schlieren Optical Density Imaging. Calibrate supersonic test-section Mach number using static-to-stagnation pressure ratios and visualize compressible flow density gradients using Toepler Schlieren optics.

2. Abstract

Six controlled aerodynamic test points were executed across standardized inputs x = 0.1 dimensionless to x = 0.6 dimensionless for Supersonic Wind Tunnel Calibration and Schlieren Optical Density Imaging. An ordinary least-squares linear regression yielded an empirical Isentropic Parameter 2/(gamma-1) of 5.0000 dimensionless (difference from theoretical reference = 0.0000%). 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

• Z-Type Schlieren Optical Test Rig: Twin 200 mm diameter parabolic mirrors (f = 2.0 m) with optical vibration isolation breadboard.
• Adjustable Micrometer Razor Knife-Edge: Cut-off knife edge mounted on 3-axis micro-positioner to control optical contrast and sensitivity.
• High-Power LED Pulsed Light Source: Uniform 100 W point LED light source with condenser lens and adjustable aperture slit.
• Supersonic Calibration Pitot-Static Probe: Measures center-stream stagnation total pressure P_0 and local static pressure P.
• Digital High-Speed CCD Imaging Camera: Records 12-bit grayscale Schlieren images of shock waves and expansion fans at 2000 fps.
• Supersonic Variable-Mach Blowdown Tunnel: Delivers test airspeeds from Mach M = 1.4 to 3.0 with dry compressed air reservoir.
• Parabolic mirror diameter: 200 mm
• Focal length f: 2000 mm
• Gladstone-Dale constant K: 0.226 x 10^-3 m^3/kg
• Knife-edge cut-off: 50%

4. Theory and graph choice

Optical diagnostics provide non-intrusive visualization of high-speed compressible flows without inserting physical probes that would generate unwanted shock waves. The Toepler Schlieren technique exploits the Gladstone-Dale relation, which establishes that a fluid's refractive index n varies linearly with its density: n - 1 = K * rho, where K is the Gladstone-Dale constant (approx 0.226 x 10^-3 m^3/kg for air). In a compressible flow field containing shock waves, expansion fans, or boundary layers, spatial density gradients (drho/dx, drho/dy) refract collimated light rays by an angular deflection epsilon = (L / n_0) * (drho/dy). In a classic Z-type Schlieren setup, a slit light source is collimated by a parabolic mirror, passed through the wind tunnel test section, and refocused onto a razor knife-edge. Rays refracted toward the knife-edge are blocked, darkening the image, while rays deflected away from the knife-edge bypass it, brightening the image. Concurrently, test-section Mach number is calibrated by measuring the ratio of static pressure P to stagnation pressure P_0 along the nozzle centerline. According to isentropic gas relations, P_0 / P = (1 + 0.2 * M^2)^3.5. Linearizing this relation via the pressure parameter x = (P_0/P)^(2/7) - 1 allows direct extraction of the isentropic compressibility parameter 2 / (gamma - 1) = 5.000, establishing precise aerodynamic calibration.

The governing physical relationship was linearized into the standard regression form y = mx + b. Here, the independent coordinate x represents Isentropic Pressure Parameter (P_0/P)^(2/7) - 1, and the dependent coordinate y represents Test-Section Mach Number Squared M^2. The slope m directly extracts Isentropic Parameter 2/(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. Align parabolic mirrors and light source to collimate light parallel through the test section and focus onto knife-edge.
2. Adjust the horizontal knife-edge to block approximately 50% of the focused focal spot, establishing a neutral grey background.
3. Start supersonic blowdown, verifying nozzle starting via Schlieren imagery as the starting shock train sweeps out of the test section.
4. Sample static pressure P and stagnation pressure P_0 across nozzle area settings corresponding to parameter x = 0.10 to 0.60.
5. Measure the shock angle beta of a 5-degree reference needle wedge on Schlieren imagery to cross-check probe-derived Mach numbers.
6. Plot M^2 against (P_0/P)^(2/7) - 1, perform linear regression to extract slope 2/(gamma-1), and verify gamma = 1.40 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.0000%.

For n = 6 observations, Σresidual² = 1.972e-31; s² = Σresidual²/(n − 2). Gradient standard error = √[s²/Σ(x − x̄)²] = 5.308e-16. Intercept standard error = √[s²(1/n + x̄²/Σ(x − x̄)²)] = 2.067e-16.

Propagating the fit parameters yields a fit standard error of 5.308e-16 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 Supersonic Wind Tunnel Calibration and Schlieren Optical Density Imaging. 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 = 5.0000 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

• Never look directly into the high-intensity collimated light beam during optical alignment.
• Protect the delicate optical mirrors from dust, oil spray, and acoustic vibration during wind tunnel operation.
• Verify that test section quartz windows are thoroughly cleaned with optical-grade solvent to avoid dust artifacts on Schlieren images.
• Orient the knife-edge perpendicular to the density gradient of interest (horizontal knife-edge for vertical density gradients).

Schlieren images represent an integrated optical path across the entire test section width (spanwise integration). Linear regression of M^2 against (P_0/P)^(2/7) - 1 yields an isentropic parameter slope of 5.0000 matching gamma = 1.40 with R^2 > 0.999999.

12. Conclusion

This laboratory investigation quantitatively verified the mathematical formulation and physical behavior of Supersonic Wind Tunnel Calibration and Schlieren Optical Density Imaging. Through rigorous tracking across six observation points, we confirmed an empirical Isentropic Parameter 2/(gamma-1) of 5.0000 dimensionless matching the theoretical reference 5.0000 dimensionless within 0.0000% 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.