Previous efforts determined a set of calibrated, optimal model parameter values for Reynolds-averaged Navier–Stokes (RANS) simulations of a compressible jet in crossflow (JIC) using a $k–ε$ turbulence model. These parameters were derived by comparing simulation results to particle image velocimetry (PIV) data of a complementary JIC experiment under a limited set of flow conditions. Here, a $k–ε$ model using both nominal and calibrated parameters is validated against PIV data acquired from a much wider variety of JIC cases, including a realistic flight vehicle. The results from the simulations using the calibrated model parameters showed considerable improvements over those using the nominal values, even for cases that were not used in the calibration procedure that defined the optimal parameters. This improvement is demonstrated using a number of quality metrics that test the spatial alignment of the jet core, the magnitudes of multiple flow variables, and the location and strengths of vortices in the counter-rotating vortex cores on the PIV planes. These results suggest that the calibrated parameters have applicability well outside the specific flow case used in defining them and that with the right model parameters, RANS solutions for the JIC can be improved significantly over those obtained from the nominal model.
Femtosecond laser electronic excitation tagging (FLEET) is a powerful unseeded velocimetry technique typically used to measure one component of velocity along a line, or two or three components from a dot. In this Letter, we demonstrate a dotted-line FLEET technique which combines the dense profile capability of a line with the ability to perform two-component velocimetry with a single camera on a dot. Our set-up uses a single beam path to create multiple simultaneous spots, more than previously achieved in other FLEET spot configurations. We perform dotted-line FLEET measurements downstream of a highly turbulent, supersonic nitrogen free jet. Dotted-line FLEET is created by focusing light transmitted by a periodic mask with rectangular slits of 1.6 × 40 mm2 and an edge-to-edge spacing of 0.5 mm, then focusing the imaged light at the measurement region. Up to seven symmetric dots spaced approximately 0.9 mm apart, with mean full-width at half maximum diameters between 150 and 350 µm, are simultaneously imaged. Both streamwise and radial velocities are computed and presented in this Letter.
Time-resolved particle image velocimetry (TR-PIV) has become widespread in fluid dynamics. Essentially a velocity field movie, the dynamic content provides temporal as well as spatial information, in contrast to conventional PIV offering only statistical ensembles of flow quantities. From these time series arise further analyses such as accelerometry, space-time correlations, frequency spectra of turbulence including spatial variability, and derivation of pressure fields and forces. The historical development of TR-PIV is chronicled, culminating in an assessment of the current state of technology in high-repetition-rate lasers and high-speed cameras. Commercialization of pulse-burst lasers has expanded TR-PIV into more flows, including the compressible regime, and has achieved MHz rates. Particle response times and peak locking during image interrogation require attention but generally are not impediments to success. Accuracy considerations are discussed, including the risks of noise and aliasing in spectral content. Oversampled TR-PIV measurements allow use of multi-frame image interrogation methods, which improve the precision of the correlation and raise the velocity dynamic range of PIV. In combination with volumetric methods and data assimilation, a full four-dimensional description of a flow is not only achievable but becoming standardized. A survey of exemplary applications is followed by a few predictions concerning the future of TR-PIV.
The development of new hypersonic flight vehicles is limited by the physical understanding that may be obtained from ground test facilities. This has motivated the present development of a temporally and spatially resolved velocimetry measurement for Sandia National Laboratories (SNL) Hypersonic Wind Tunnel (HWT) using Femtosecond Laser Electronic Excitation Tagging (FLEET). First, a multi-line FLEET technique has been created for the first time and tested in a supersonic jet, allowing simultaneous measurements of velocities along multiple profiles in a flow. Secondly, two different approaches have been demonstrated for generating dotted FLEET lines. One employs a slit mask pattern focused into points to yield a dotted line, allowing for two- or three-component velocity measurements free of contamination between components. The other dotted-line approach is based upon an optical wedge array and yields a grid of points rather than a dotted line. Two successful FLEET measurement campaigns have been conducted in SNL’s HWT. The first effort established optimal diagnostic configurations in the hypersonic environment based on earlier benchtop reproductions, including validation of the use of a 267 nm beam to boost the measurement signal-to-noise ratio (SNR) with minimal risk of perturbing the flow and greater simplicity than a comparable resonant technique at 202 nm. The same FLEET system subsequently was reconstituted to demonstrate the ability to make velocimetry measurements of hypersonic turbulence in a realistic flow field. Mean velocity profiles and turbulence intensity profiles of the shear layer in the wake of a hypersonic cone model were measured at several different downstream stations, proving the viability of FLEET as a hypersonic diagnostic.
This experimental study explores the fluid-structure interactions occurring between a control surface and the hypersonic flow deflected by it. The control surface is simplified for this work as a spanwise finite wedge placed on a longitudinally sliced part of the cone. The front surface of the wedge is a thin panel which is designed to respond to the unsteady fluid loading arising from the shock-wave/boundary layer interactions. Experiments have been conducted in the Sandia Hypersonic Wind Tunnel at Mach 5 and Mach 8 at wedge angles of 10◦, 20◦ and 30◦ . High-speed schlieren and backside panel accelerometer measurements capture the unsteady flow dynamics and structural response of the thin panel, respectively. For attached or small separation interactions, the transitional regime has the strongest panel fluctuations with convective shock undulations induced by the boundary layer disturbance shown to be associated with dominant panel vibrations. For large separated interactions, shear layer flapping can excite select panel modes. Heating of the panel causes a downward shift in natural mode frequencies.
This study seeks to simplify the optical requirements for multi-line FLEET (Femtosecond Laser Electronic Excitation Tagging) generation by focusing the image of a periodic slit-mask with a cylindrical and spherical lens. Geometry effects on the signal were analyzed over fifteen mask iterations. The signal for each mask was found to vary with mask standoff from the focusing optics, which was optimized based on maximizing the Signal-to-Noise Ratio (SNR) for each mask. The number of generated lines was found to decrease with slit spacing while the separation of the lines increased. FLEET line spacing was determined by a constant magnification value of the imaged masks’ slit spacing. From the geometry study, two masks that produced three to five lines spaced at 0.8–1 mm apart with SNR > 4 were chosen to demonstrate the multi-line technique in a supersonic free-jet. Velocity calculations from this data showed good agreement with schlieren imaging of compressible flow structures.
Multi-frame correlation algorithms for time-resolved PIV have been shown in previous studies to reduce noise and error levels in comparison with conventional two-frame correlations. However, none of these prior efforts tested the accuracy of the algorithms in spectral space. Even should a multi-frame algorithm reduce the error of vector computations summed over an entire data set, this does not imply that these improvements are observed at all frequencies. The present study examines the accuracy of velocity spectra in comparison with simultaneous hot-wire data. Results indicate that the high-frequency content of the spectrum is very sensitive to choice of the interrogation algorithm and may not return an accurate response. A top-hat-weighted sliding sum-of-correlation is contaminated by high-frequency ringing whereas Gaussian weighting is indistinguishable from a low-pass filtering effect. Some evidence suggests the pyramid correlation modestly increases bandwidth of the measurement at high frequencies. The apparent benefits of multi-frame interrogation algorithms may be limited in their ability to reveal additional spectral content of the flow.
The primary parameter of a standard k-ϵ model, Cμ, was calculated from stereoscopic particle image velocimetry (PIV) data for a supersonic jet exhausting into a transonic crossflow. This required the determination of turbulent kinetic energy, turbulent eddy viscosity, and turbulent energy dissipation rate. Image interrogation was optimized, with different procedures used for mean strain rates and Reynolds stresses, to produce useful turbulent eddy viscosity fields. The eddy viscosity was calculated by a least-squares fit to all components of the three-dimensional strain-rate tensor that were available from the PIV data. This eliminated artifacts and noise observed when using a single strain component. Local dissipation rates were determined via Kolmogorov’s similarity hypotheses and the second-order structure function. The eddy viscosity and dissipation rates were then combined to determine Cμ. Considerable spatial variation was observed in Cμ, with the highest values found in regions where turbulent kinetic energy was relatively ow but where turbulent mixing was important, e.g., along the high-strain jet edges and in the wake region. This suggests that use of a constant Cμ in modeling may lead to poor Reynolds stress predictions at mixing interfaces. A data-driven modeling approach that can predict this spatial variation of Cμ based on known state variables may lead to improved simulation results without the need for calibration.
A simple linear configuration for multi-line femtosecond laser electronic excitation tagging (FLEET) velocimetry is used for the first time, to the best of our knowledge, to image an overexpanded unsteady supersonic jet. The FLEET lines are spaced 0.5-1.0mmapart, and up to six lines can be used simultaneously to visualize the flowfield. These lines are created using periodic masks, despite the mask blocking 25%-30%of the 10 mJ incident beam.Maps of mean singlecomponent velocity in the direction along the principal flow axis, and turbulence intensity in that same direction, are created using multi-line FLEET, and computed velocities agree well with those obtained from single-line (traditional) FLEET. Compared to traditional FLEET, multi-line FLEET offers increased simultaneous spatial coverage and the ability to produce spatial correlations in the streamwise direction. This FLEET permutation is especially well suited for short-duration test facilities.
Two techniques have extended the effective frequency limits of postage-stamp PIV, in which a pulse-burst laser and very small fields of view combine to achieve high repetition rates. An interpolation scheme reduced measurement noise, raising the effective frequency response of previous 400-kHz measurements from about 120 kHz to 200 kHz. The other technique increased the PIV acquisition rate to very nearly MHz rates (990 kHz) by using a faster camera. Charge leaked through the camera shift register at these framing rates but this was shown not to bias the measurements. The increased framing rate provided oversampled data and enabled use of multi-frame correlation algorithms for a lower noise floor, increasing the effective frequency response to 240 kHz where the interrogation window size begins to spatially filter the data. Good agreement between the interpolation technique and the MHz-rate PIV measurements was established. The velocity spectra suggest turbulence power-law scaling in the inertial subrange steeper than the theoretical-5/3 scaling, attributed to an absence of isotropy.
A blind CFD validation challenge is being organized for the unsteady transonic shock motion induced by the Sandia Axisymmetric Transonic Hump, which echoes the Bachalo-Johnson configuration. The wind tunnel and model geometry will be released at the start of the validation challenge along with flow boundary conditions. Primary data concerning the unsteady separation region will be released at the conclusion of the challenge after computational entrants have been submitted. This paper details the organization of the challenge, its schedule, and the metrics of comparison by which the models will be assessed.
Bench-top tests are conducted to characterize Femtosecond Laser Electronic Excitation Tagging (FLEET) in static low pressure (35 mTorr-760 Torr) conditions, and to measure the acoustic disturbance caused by the resulting filament as a function of tagging wavelength and energy. The FLEET line thickness as a function of pressure and delay is described by a simple diffusion model. Initial FLEET measurements in a Mach 8 flow show that gate times of ≥ 1µs can produce visible smearing of the FLEET emission and challenge the traditional Gaussian fitting methods used to find the line center. To minimize flow perturbations and uncertainty of the final line position, several recommendations are offered: using third harmonic FLEET at 267 nm for superior signal levels with lower energy deposition than both 800 nm and 400 nm FLEET, and short camera delays and exposure times to reduce fitting uncertainty. This guidance is implemented in a Mach 8 test condition and results are presented.
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.
Here, experiments were performed within Sandia National Labs’ Multiphase Shock Tube to measure and quantify the shock-induced dispersal of a shock/dense particle curtain interaction. Following interaction with a planar travelling shock wave, schlieren imaging at 75 kHz was used to track the upstream and downstream edges of the curtain. Data were obtained for two particle diameter ranges ($d_{p}=106{-}125$,$300{-}355~\unicode[STIX]{x03BC}\text{m}$) across Mach numbers ranging from 1.24 to 2.02. Using these data, along with data compiled from the literature, the dispersion of a dense curtain was studied for multiple Mach numbers (1.2–2.6), particle sizes ($100{-}1000~\unicode[STIX]{x03BC}\text{m}$) and volume fractions (9–32 %). Data were non-dimensionalized according to two different scaling methods found within the literature, with time scales defined based on either particle propagation time or pressure ratio across a reflected shock. The data refelct that spreading of the particle curtain is a function of the volume fraction, with the effectiveness of each time scale based on the proximity of a given curtain’s volume fraction to the dilute mixture regime. It is observed that volume fraction corrections applied to a traditional particle propagation time scale result in the best collapse of the data between the two time scales tested here. In addition, a constant-thickness regime has been identified, which has not been noted within previous literature.
Time-resolved particle image velocimetry was conducted at 40 kHz using a pulse-burst laser in the supersonic wake of a wall-mounted hemisphere. Velocity fields suggest a recirculation region with two lobes, in which flow moves away from the wall near the centerline and recirculates back toward the hemisphere off the centerline, contrary to transonic configurations. Spatio-temporal cross-correlations and conditional ensemble averages relate the characteristic behavior of the unsteady shock motion to the flapping of the shear layer. At Mach 1.5, oblique shocks develop, associated with vortical structures in the shear layer and convect downstream in tandem; a weak periodicity is observed. Shock motion at Mach 2.0 appears somewhat different, wherein multiple weak disturbances propagate from shear-layer turbulent structures to form an oblique shock that ripples as these vortices pass by. Bifurcated shock feet coalesce and break apart without evident periodicity. Power spectra show a preferred frequency of shear-layer flapping and shock motion for Mach 1.5, but at Mach 2.0, a weak preferred frequency at the same Strouhal number of 0.32 is found only for oblique shock motion and not shear-layer unsteadiness.
A new wind tunnel experiment is underway to provide a comprehensive CFD validation dataset of an unsteady, transonic flow. The experiment is based on the work of Bachalo and Johnson; an axisymmetric model with a spherical hump is tested at a transonic Mach number. 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. A suite of diagnostics characterizes the flow: oil-flow surface visualization for shock and reattachment locations, particle image velocimetry for mean flow and turbulence properties, fast pressure-sensitive paint for model pressure distributions and unsteadiness, high-speed Schlieren for shock position and motion, and oil-film interferometry for surface shear stress. This will provide a new level of detail for validation studies; therefore, a blind comparison, or ‘CFD Challenge’ is proposed to the community. Participants are to be provided the geometry, incoming boundary layer, and boundary conditions, and are free to simulate with their method of choice and submit their results. A blind comparison will be made to the new experimental data, with the goal of evaluating the state of various CFD methods for use in unsteady, transonic flows.
Femtosecond Laser Electronic Excitation Tagging (FLEET) is used to measure velocity flowfields in the wake of a sharp 7◦ half-angle cone in nitrogen at Mach 8, over freestream Reynolds numbers from 4.3∗106 /m to 13.8∗106 /m. Flow tagging reveals expected wake features such as the separation shear layer and two-dimensional velocity components. Frequency-tripled FLEET has a longer lifetime and is more energy efficient by tenfold compared to 800 nm FLEET. Additionally, FLEET lines written with 267 nm are three times longer and 25% thinner than that written with 800 nm at a 1 µs delay. Two gated detection systems are compared. While the PIMAX 3 ICCD offers variable gating and fewer imaging artifacts than a LaVision IRO coupled to a Photron SA-Z, its slow readout speed renders it ineffective for capturing hypersonic velocity fluctuations. FLEET can be detected to 25 µs following excitation within 10 mm downstream of the model base, but delays greater than 4 µs have deteriorated signal-to-noise and line fit uncertainties greater than 10%. In a hypersonic nitrogen flow, exposures of just several hundred nanoseconds are long enough to produce saturated signals and/or increase the line thickness, thereby adding to measurement uncertainty. Velocity calculated between the first two delays offer the lowest uncertainty (less than 3% of the mean velocity).
The mechanism by which aerodynamic effects of jet/fin interaction arise from the flow structure of a jet in crossflow is explored using particle image velocimetry measurements of the crossplane velocity field as it impinges on a downstream fin instrumented with high-frequency pressure sensors. A Mach 3.7 jet issues into a Mach 0.8 crossflow from either a normal or inclined nozzle, and three lateral fin locations are tested. Conditional ensemble-averaged velocity fields are generated based upon the simultaneous pressure condition. Additional analysis relates instantaneous velocity vectors to pressure fluctuations. The pressure differential across the fin is driven by variations in the spanwise velocity component, which substitutes for the induced angle of attack on the fin. Pressure changes at the fin tip are strongly related to fluctuations in the streamwise velocity deficit, wherein lower pressure is associated with higher velocity and vice versa. The normal nozzle produces a counter-rotating vortex pair that passes above the fin, and pressure fluctuations are principally driven by the wall horseshoe vortex and the jet wake deficit. In conclusion, the inclined nozzle produces a vortex pair that impinges the fin and yields stronger pressure fluctuations driven more directly by turbulence originating from the jet mixing.
Fluid-structure interactions were studies on a 7° half-angle cone in the Sandia Hypersonic Wind Tunnel at Mach 5 and 8 and in the Purdue Boeing/AFOSR Mach 6 Quiet Tunnel. A thin composite panel was integrated into the cone and the response to boundary-layer disturbances was characterized by accelerometers on the backside of the panel. Here, under quiet-flow conditions at Mach 6, the cone boundary layer remained laminar. Artificially generated turbulent spots excited a directionally dependent panel response which would last much longer than the spot duration.
The mechanism by which aerodynamic effects of jet/fin interaction arise from the flow structure of a jet in crossflow is explored using particle image velocimetry measurements of the crossplane velocity field as it impinges on a downstream fin instrumented with high-frequency pressure sensors. A Mach 3.7 jet issues into a Mach 0.8 crossflow from either a normal or inclined nozzle, and three lateral fin locations are tested. Conditional ensemble-averaged velocity fields are generated based upon the simultaneous pressure condition. Additional analysis relates instantaneous velocity vectors to pressure fluctuations. The pressure differential across the fin is driven by variations in the spanwise velocity component, which substitutes for the induced angle of attack on the fin. Pressure changes at the fin tip are strongly related to fluctuations in the streamwise velocity deficit, wherein lower pressure is associated with higher velocity and vice versa. The normal nozzle produces a counter-rotating vortex pair that passes above the fin, and pressure fluctuations are principally driven by the wall horseshoe vortex and the jet wake deficit. The inclined nozzle produces a vortex pair that impinges the fin and yields stronger pressure fluctuations driven more directly by turbulence originating from the jet mixing.