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 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.
Compressible wall modeled large-eddy simulations of a Mach eight turbulent boundary layer over a flat plate were carried out for the conditions of the Hypersonic Wind Tunnel at Sandia National Laboratories. Overall good agreement of the velocity and temperature profiles is obtained with reference data from a direct numerical simulation and a theoretical relationship. Profiles of the resolved root-mean-square velocity fluctuations are in adequate agreement with the reference data. The refractive index is calculated from the density field and integrated along an expected beam path to calculate the optical path length. Then, by subtracting a bilinear fit of the instantaneous optical path length, the optical path difference is obtained. The computed aero-optical path difference shows a similar dependence on the aperture size as in the literature. The normalized root-mean-square optical path difference from the present wall-modeled large-eddy simulations and a reference direct numerical simulation and experiment are in good agreement. The optical path distortion is slightly above the value predicted by a semi-analytical relationship from the literature. Finally, instantaneous snapshots of the flow are analyzed via proper orthogonal decomposition and the optical path distortion is computed from subsets of the modes. The optical path distortion converges quickly with increasing number of modes which suggests that the main contribution comes from large energetic flow structures.
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.
Four Direct Numerical Simulation (DNS) datasets covering effective freestream Mach numbers of 8 through 14 are used to investigate the behavior of turbulence-induced aero-optical distortions in hypersonic boundary layers. The datasets include two from simulations of flat plate boundary layers (Mach 8 and 14) and two from simulations of flow over a sharp cone (Mach 8 and 14). Instantaneous three-dimensional fields of density from each DNS are converted to refraction index and integrated to produce distributions of the Optical Path Differences (OPD) caused by turbulence. These values are then compared to experimental data from the literature and to an existing model for the root-mean-square of the OPD. Although the model was originally developed for flows with Mach ≤ 5, it provides a basis to which we compare the hypersonic data.
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.
The character of aero-optical distortions produced by turbulence is investigated for subsonic, supersonic, and hypersonic boundary layers. Data from four Direct Numerical Simulations (DNS) of boundary layers with nominal Mach numbers ranging from 0.5 to 8 are used. The DNS data for the subsonic and supersonic boundary layers are of flow over flat plates. Two hypersonic boundary layers are both from flows with a Mach 8 inlet condition, one of which is flow over a flat plate while the other is a boundary layer on a sharp cone. Density fields from these datasets are converted to index-of-refraction fields which are integrated along an expected beam path to determine the effective Optical Path Lengths that a beam would experience while passing through the refractions of the turbulent field. By then accounting for the mean path length and tip/tilt issues related to bulk boundary layer effects, the distribution of Optical Path Differences (OPD s) is determined. Comparisons of the root-mean-squares of the OPDs are made to an existing model. The OPDr m s values determined from the subsonic and supersonic data were found to match the existing model well. As could be expected, the hypersonic data does not match as well due to assumptions like the Strong Reynold Analogy that were made in the derivation of the model. Until now, the model has never been compared to flows with Mach numbers as high as included herein or to flow over a sharp cone geometry.
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.
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.
Measurements are presented of the aero-optic distortion produced by a Mach 8 turbulent boundary layer in the Sandia Hypersonic Wind Tunnel. Flat optical inserts installed in the test section walls enabled a double-pass arrangement of a collimated laser beam. The distortion of this beam was imaged by a high-speed Shack-Hartmann sensor at a sampling rate of up to 1 MHz. Analysis is performed using two processing methods to extract the aero-optic distortion from the data. A novel de-aliasing algorithm is proposed to extract convective-only spectra and is demonstrated to correctly quantify the physical spectra even in case of relatively low sampling rates. The results are compared with an existing theoretical model, and it is shown that this model under-predicts the experimentally measured distortions regardless of the processing method used. Possible explanations for this discrepancy are presented. The presented results represent to-date the highest Mach number for which aero-optic boundary layer distortion measurements are available.
An experimental characterization of the flow environment for the Sandia Axisymmetric Transonic Hump is presented. This is an axisymmetric model with a circular hump tested at a transonic Mach number, similar to the classic Bachalo-Johnson configuration. The flow is turbulent approaching the hump and becomes locally supersonic at the apex. This leads to a shock-wave/boundary-layer interaction, an unsteady separation bubble, and flow reattachment downstream. The characterization focuses on the quantities required to set proper boundary conditions for computational efforts described in the companion paper, including: 1) stagnation and test section pressure and temperature; 2) turbulence intensity; and 3) tunnel wall boundary layer profiles. Model characterization upstream of the hump includes: 1) surface shear stress; and 2) boundary layer profiles. Note: Numerical values characterizing the experiment have been redacted from this version of the paper. Model geometry and boundary conditions will be withheld until the official start of the Validation Challenge, at which time a revised version of this paper will become available. Data surrounding the hump are considered final results and will be withheld until completion of the Validation Challenge.