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.
A new reflected shock tunnel has been commissioned at Sandia capable of generating hypersonic environments at realistic flight enthalpies. The tunnel uses an existing free-piston driver and shock tube coupled to a conical nozzle to accelerate the flow to approximately Mach 9. The facility design process is outlined and compared to other ground test facilities. A representative flight enthalpy condition is designed using an in-house state-to-state solver and piston dynamics model and evaluated using quasi-1D modeling with the University of Queensland L1d code. This condition is demonstrated using canonical models and a calibration rake. A 25 cm core flow with 4.6 MJ/kg total enthalpy is achieved over an approximately 1 millisecond test time. Analysis shows that increasing piston mass should extend test time by a factor of 2-3.
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.
Here we present results from experiments within Sandia National Labs’ multiphase shock tube on the shock-induced dispersal of dense particle curtains. This study builds on previous work by examining the effect of particle density on the dynamics of a shock-particle interaction in a dense volume fraction regime. We present results gathered from high-speed schlieren images used to track the propagation of the upstream and downstream fronts of the particle curtain. The effect of particle density on the curtain spread rate was examined by comparing curtains comprised of soda lime, stainless steel, and tungsten particles at two distinct volume fractions ϕp = 9% and ϕp ≈ 20%, and various incident shock strengths. Time scales of the spreading process were non-dimensionalized using two scaling methods from literature; one defined by the pressure ratio across a reflected shock and the other related to the incompressible drag through a grid. Both scaling methods successfully collapsed the spreading rate of curtains with different particle densities, while only the drag based scaling could account for variation in volume fraction. In addition, a new scaling based on a simple force balance that uses the pressure ratio across the curtain was found to achieve the tightest collapse of all methods tested.
A high-speed thermometry diagnostic based on spontaneous Raman scattering (SRS) was demonstrated using a pulse-burst laser at a 3-kHz data acquisition rate, with a pulse duration of 200 ns and wavelength of 532 nm. The technique was evaluated in a challenging free-piston shock tube environment operated at conditions up to 1653 K and 112 bar following an incident shock Mach number of 3.5 and a reflected shock Mach number of 2.2. The SRS thermometry resolved the temperature in post-incident and post-reflected shock flows and clearly showed cooling associated with driver expansion waves. A detailed spectral physics model inferred temperatures within 1% of the predicted post-shock temperatures, when SNR was greater than 2.0. This was a significant advancement of spontaneous Raman vibrational thermometry.
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.