High-enthalpy hypersonic flight represents an application space of significant concern within the current national-security landscape. The hypersonic environment is characterized by high-speed compressible fluid mechanics and complex reacting flow physics, which may present both thermal and chemical nonequilibrium effects. We report on the results of a three-year LDRD effort, funded by the Engineering Sciences Research Foundation (ESRF) investment area, which has been focused on the development and deployment of new high-speed thermochemical diagnostics capabilities for measurements in the high-enthalpy hypersonic environment posed by Sandia's free-piston shock tunnel. The project has additionally sponsored model development efforts, which have added thermal nonequilibrium modeling capabilities to Sandia codes for subsequent design of many of our shock-tunnel experiments. We have cultivated high-speed, chemically specific, laser-diagnostic approaches that are uniquely co-located with Sandia's high-enthalpy hypersonic test facilities. These tools include picosecond and nanosecond coherent anti-Stokes Raman scattering at 100-kHz rates for time-resolved thermometry, including thermal nonequilibrium conditions, and 100-kHz planar laser-induced fluorescence of nitric oxide for chemically specific imaging and velocimetry. Key results from this LDRD project have been documented in a number of journal submissions and conference proceedings, which are cited here. The body of this report is, therefore, concise and summarizes the key results of the project. The reader is directed toward these reference materials and appendices for more detailed discussions of the project results and findings.
X-ray stereo digital image correlation (DIC) measurements were performed at 10 kHz on the internal surface of a jointed structure in a shock tube at a shock Mach number of 1.42 and compared with optical stereo DIC measurements on the outer, visible surface of the structure. The shock tube environment introduces temperature and density gradients in the gas through which the structure was imaged, resulting in spatial and temporal index of refraction variations. These variations cause bias errors in optical DIC measurements due to beam-steering but have minimal influence on x-ray DIC measurements. These results demonstrate the utility of time-resolved x-ray DIC measurements in complicated environments where optical measurements suffer severe errors and/or are precluded by lack of optical access.
Aero-optics refers to optical distortions due to index-of-refraction gradients that are induced by aerodynamic density gradients. At hypersonic flow conditions, the bulk velocity is many times the speed of sound and density gradients may originate from shock waves, compressible turbulent structures, acoustic waves, thermal variations, etc. Due to the combination of these factors, aero-optic distortions are expected to differ from those common to sub-sonic and lower super-sonic speeds. This report summarizes the results from a 2019-2022 Laboratory Directed Research and Development (LDRD) project led by Sandia National Laboratories in collaboration with the University of Notre Dame, New Mexico State University, and the Georgia Institute of Technology. Efforts extended experimental and simulation methodologies for the study of turbulent hypersonic boundary layers. Notable experimental advancements include development of spectral de-aliasing techniques for highspeed wavefront measurements, a Spatially Selective Wavefront Sensor (SSWFS) technique, new experimental data at Mach 8 and 14, a Quadrature Fringe Imaging Interferometer (QFII) technique for time-resolved index-of-refraction measures, and application of QFII to shock-heated air. At the same time, model advancements include aero-optic analysis of several Direct Numerical Simulation (DNS) datasets from Mach 0.5 to 14 and development of wall-modeled Large Eddy Simulation (LES) techniques for aero-optic predictions. At Mach 8 measured and predicted root mean square Optical Path Differences agree within confidence bounds but are higher than semi-empirical trends extrapolated from lower Mach conditions. Overall, results show that aero-optic effects in the hypersonic flow regime are not simple extensions from prior knowledge at lower speeds and instead reflect the added complexity of compressible hypersonic flow physics.
This work presents a high-speed laser-absorption-spectroscopy diagnostic capable of measuring temperature, pressure, and nitric oxide (NO) mole fraction in shock-heated air at a measurement rate of 500 kHz. This diagnostic was demonstrated in the High-Temperature Shock Tube (HST) facility at Sandia National Laboratories. The diagnostic utilizes a quantum-cascade laser to measure the absorbance spectra of two rovibrational transitions near 5.06 µm in the fundamental vibration bands (v" = 0 and 1) of NO in its ground electronic state (X2 Π1/2 ). Gas properties were determined using scanned-wavelength direct absorption and a recently established fitting method that utilizes a modified form of the time-domain molecular free-induction-decay signal (m-FID). This diagnostic was applied to acquire measurements in shock-heated air in the HST at temperatures ranging from approximately 2500 to 5500 K and pressures of 3 to 12 atm behind both incident and reflected shocks. The measurements agree well with the temperature predicted by NASA CEA and the pressure measured simultaneously using PCB pressure sensors. The measurements presented demonstrate that this diagnostic is capable of resolving the formation of NO in shock-heated air and the associated temperature change at the conditions studied.
Coherent anti-Stokes Raman scattering of the N2 molecule is performed at rates up to 100 kHz for thermometry in the Sandia free-piston, high-temperature shock-tube facility (HST) for reflected-shock conditions in excess of T = 4000 K at pressures up to P = 10 atm. A pulse-burst laser architecture delivers picosecond-duration pulses to provide both the CARS pump and probe photons, and to pump a solid-state optical parametric generator (OPG)/optical parametric amplifier (OPA) source, which provides frequency tunable Stokes pulses with a bandwidth of 100-120 cm-1 . Single-laser-shot and averaged CARS spectra obtained in both the incident (P = 1.1 atm, T = 2090 K) and reflected (P ~ 8-10.5 atm, T > 4000 K) shock regions of HST are presented. The results indicate that burst-mode CARS is capable of resolving impulsive, high-temperature events in HST.
This work presents a high-speed laser-absorption-spectroscopy diagnostic capable of measuring temperature, pressure, and nitric oxide (NO) mole fraction in shock-heated air at a measurement rate of 500 kHz. This diagnostic was demonstrated in the High-Temperature Shock Tube (HST) facility at Sandia National Laboratories. The diagnostic utilizes a quantum-cascade laser to measure the absorbance spectra of two rovibrational transitions near 5.06 µm in the fundamental vibration bands (v" = 0 and 1) of NO in its ground electronic state (X2 Π1/2 ). Gas properties were determined using scanned-wavelength direct absorption and a recently established fitting method that utilizes a modified form of the time-domain molecular free-induction-decay signal (m-FID). This diagnostic was applied to acquire measurements in shock-heated air in the HST at temperatures ranging from approximately 2500 to 5500 K and pressures of 3 to 12 atm behind both incident and reflected shocks. The measurements agree well with the temperature predicted by NASA CEA and the pressure measured simultaneously using PCB pressure sensors. The measurements presented demonstrate that this diagnostic is capable of resolving the formation of NO in shock-heated air and the associated temperature change at the conditions studied.
Density fluctuations in compressible turbulent boundary layers cause aero-optical distortions that affect the performance of optical systems such as sensors and lasers. The development of models for predicting the aero-optical distortions relies on theory and reference data that can be obtained from experiments and time-resolved simulations. This paper reports on wall-modeled large-eddy simulations of turbulent boundary layers over a flat plate at Mach 3.5, 7.87, and 13.64. The conditions for the Mach 3.5 case match those for the DNS presented by Miller et al.1 The Mach 7.87 simulation match those inside the Hypersonic Wind Tunnel at Sandia National Laboratories. For the Mach 13.64, the conditions inside the Arnold Engineering Development Complex Hypervelocity Tunnel 9 are matched. Overall, adequate agreement of the velocity and temperature as well as Reynolds stress profiles with reference data from direct numerical simulations is obtained for the different Mach numbers. For all three cases, the normalized root-mean-square optical path difference was computed and compared with data obtained from the reference direct numerical simulations and experiments, as well as predictions obtained with a semi-analytical relationship by Notre Dame University. Above Mach five, the normalized path difference obtained from the simulations is above the model prediction. This provides motivation for future work aimed at evaluating the assumptions behind the Notre Dame model for hypersonic boundary layer flows.
A high-speed temperature diagnostic based on spontaneous Raman scattering (SRS) was demonstrated using a pulse-burst laser. The technique was first benchmarked in near-adiabatic H2-air flames at a data-acquisition rate of 5 kHz using an integrated pulse energy of 1.0 J per realization. Both the measurement precision and accuracy in the flame were within 3% of adiabatic predictions. This technique was then evaluated in a challenging free-piston shock tube environment operated at a shock Mach number of 3.5. SRS thermometry resolved the temperature in post-incident and post-reflected shock flows at a repetition rate of 3 kHz and clearly showed cooling associated with driver expansion waves. Collectively, this Letter represents a major advancement for SRS in impulsive facilities, which had previously been limited to steady state regions or single-shot acquisition.
Measurements of bifurcated reflected shocks over a wide range of incident shock Mach numbers, 2.9 < Ms < 9.4, are carried out in Sandia’s high temperature shock tube. The size of the non-uniform flow region associated with the bifurcation is measured using high speed schlieren imaging. Measurements of the bifurcation height are compared to historical data from the literature. A correlation for the bifurcation height from Petersen et al. [1] is examined and found to over estimate the bifurcation height for Ms > 6. An improved correlation is introduced that can predict the bifurcation height over the range 2.15 < Ms < 9.4. The time required for the non-uniform flow region to pass over a stationary sensor is also examined. A non-dimensional time related to the induced velocity behind the shock and the distance to the endwall is introduced. This non-dimensional time collapses the data and yields a new correlation that predicts the temporal duration of the bifurcation.
This paper validates the concept of a spatially filtered wavefront sensor, which uses a convergent-divergent beam to reduce sensitivity to aero-optical distortions near the focal point while retaining sensitivity at large beam diameters. This sensor was used to perform wavefront measurements in a cavity flow test section. The focal point was traversed to various spanwise locations across the test section, and the overall OPDRMS levels and aperture-averaged spectra of wavefronts were computed. It was demonstrated that the sensor was able to effectively suppress the stronger aero-optical signal from the cavity flow and recover the aero-optical signal from the boundary layer when the focal point was placed inside the shear region of the cavity flow. To model these measured quantities, additional collimated beam wavefronts were taken at various subsonic speeds in a wind tunnel test section with two turbulent boundary layers, and then in the cavity flow test section, where the signal from the cavity was dominant. The results from the experimental model agree with the measured convergent-divergent beam results, confirming that the spatial filtering properties of the proposed sensor are due to attenuating effects at small apertures.
