Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hypersonic flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effects of non-thermal termolecular reactions on wedge-induced oblique detonation waves

The shock-induced combustion ramjet (Shcramjet) based on oblique detonation waves (ODWs) is among the promising choices for hypersonic propulsion systems. An understanding of the ignition, propagation, and stability of ODWs is critical to harnessing their propulsive potential. In such high speed reacting flows, there is a high probability of occurrence of non-thermal reactions due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 as demonstrated recently [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078–1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. The present work focuses on examining the initiation, propagation and structure of oblique detonation waves in stoichiometric H 2 -air mixtures through numerical simulations with and without non-thermal reactivity on a two-dimensional adaptive grid. Non-thermal reactions were included in the macroscopic kinetic model as chemically termolecular reactions facilitated by the H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Since, the non-thermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions, an additional simulation was performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. Results show that, under ODWE relevant conditions, non-thermal reactivity fundamentally alters the induction length, intensity as well as the structure of the ODW. Specifically, it is found that non-thermal reactivity leads to a noticeable reduction in initiation length and a simultaneous increase in instantaneous peak heat release rate and the degree of unsteadiness of the ODW. Finally, statistical analysis of key thermo-chemical variables is performed to elucidate the important species as well as reactions responsible for the observed variations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Impact of Vibrational Nonequilibrium on the Supersonic Mode in Hypersonic Boundary Layers

The supersonic mode in hypersonic boundary layers has recently been shown to achieve a larger peak growth rate than the traditional second mode in particular circumstances relevant to free flight. In this environment, high levels of aerodynamic heating are encountered such that incorporation of chemical and vibrational nonequilibrium effects can become critical. The impact of thermal nonequilibrium (TNE) on the supersonic mode in particular has not been thoroughly and systematically investigated. This work uses thermochemical nonequilibrium (TCNE) direct numerical simulation (DNS) and TCNE linear stability theory (LST) to obtain a more complete investigation of the supersonic mode using both nonequilibrium and frozen gas models. The simulation is Mach 5 flow over a 1-mm-nose-radius axisymmetric cone 1 m in length, representative of experimental conditions in ground test facilities. Both LST and DNS results indicate that TNE effects are destabilizing compared with frozen flow and can lead to an N-factor difference of approximately +1 at the end of the cone. LST results also indicate that TNE effects are stabilizing to the second and supersonic modes in comparison with thermal equilibrium. The supersonic mode appeared to be slightly more affected by TNE than the second mode. In conclusion, collectively, the findings in this work suggest that even in ground test facilities where TCNE effects are relatively small compared with free flight, incorporating nonequilibrium models in both the mean flow and stability solvers can have a practical effect on hypersonic boundary-layer transition prediction analyses.

42 ENGINEERING↗
DOE OSTI2020

Coupling of Laminar-Turbulent Transition with RANS Computational Fluid Dynamics

This project combines several new concepts to create a boundary layer transition prediction capability that is suitable for analyzing modern hypersonic flight vehicles. The first new concept is the use of ''optimization'' methods to detect the hydrodynamic instabilities that cause boundary layer transition; the use of this method removes the need for many limiting assumptions of other methods and enables quantification of the interactions between boundary layer instabilities and the flow field imperfections that generate them. The second new concept is the execution of transition analysis within a conventional hypersonics CFD code, using the same mesh and numerical schemes for the transition analysis and the laminar flow simulation. This feature enables rapid execution of transition analysis with less user oversight required and no interpolation steps needed.

42 ENGINEERING↗
DOE OSTI2023

Arrays of glass wedges for multi-dimensional optical diagnostics

There is a common need in the advancement of optical diagnostic techniques to increase the dimensionality of measurements. For example, point measurements could be improved to multi-point, line, planar, volumetric, or time-resolved volumetric measurements. In this work, a unique optical element is presented to enable multi-dimensional measurements, namely, an array of glass wedges. A light source is passed through the wedges, and different portions of the illumination are refracted by different amounts depending on the glass wedge angle. Subsequent optics can be used to focus the light to multiple points, lines, or planes. Basic characterization of a glass wedge array is presented. Additional wedge-array configurations are discussed, including the use of a periodic intensity mask for multi-planar measurements via structured illumination. The utility of this optical element is briefly demonstrated in (a) multi-planar flame particulate measurements, (b) multi-point femtosecond-laser electronic excitation tagging for flow velocimetry, and (c) multi-line nitric oxide molecular tagging velocimetry in a hypersonic shock-tunnel. One significant advantage of this optical component is its compatibility with high-energy laser sources, which may be a limiting factor with other beam-splitting or beam-forming elements such as some diffractive optics. Additionally, an array of glass wedges is simple and easily customizable compared to other methods for forming multiple closely spaced illumination patterns. In conclusion, suggestions for further development and applications are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Assessment of the Amplification Factor Transport Transition Model for High-Mach Number Flows

In its present state, the Amplification Factor Transport (AFT) model provides an efficient alternative for high-fidelity stability methods for complex geometries subject to low-speed conditions considered in the industrial-production CFD environment. Its performance has not yet been documented for the typical transition scenarios in hypersonic conditions. The goal of the present study is to apply the AFT model to a yawed- and flared-cone configuration, for which the laminar flows respectively support the crossflow and second-mode instability mechanisms, and to compare the amplification results against high-fidelity results computed with Linear Parabolized Stability Equations (LPSE). High-resolution viscous-flow solutions are provided with OVERFLOW and COFFE. For the yawed-cone case, the AFT model yields the largest amplification in a region of the flow where the expected crossflow instability is not amplified (in the leeward symmetry plane) according to LPSE. For the flared-cone case with a wall-temperature slightly lower than the adiabatic value, the AFT model yields extremely large amplification factors (in excess of 100), exceeding the computed LPSE amplification of second-mode instability results by a factor 14. Upon considering a cold-wall condition, the model shows much smaller amplification factors, which is the opposite with respect to the expected behavior of the second mode. The AFT model does not adequately capture the crossflow and second-mode instability mechanisms in its present state.

Groot, Koen J.↗
DOE OSTI2022

Comprehensive uncertainty quantification (UQ) for full engineering models by solving probability density function (PDF) equation

This report details a new method for propagating parameter uncertainty (forward uncertainty quantification) in partial differential equations (PDE) based computational mechanics applications. The method provides full-field quantities of interest by solving for the joint probability density function (PDF) equations which are implied by the PDEs with uncertain parameters. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known apriori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the driving (dynamics) term of the PDEs based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the dynamics term of the governing PDEs. The efficacy of the method is demonstrated by applying it to two target problems: solid mechanics explicit dynamics with uncertain material model parameters, and reacting hypersonic fluid mechanics with uncertain chemical kinetic rate parameters. A minimally invasive implementation of the method for representative codes SPARC (reacting hypersonics) and NimbleSM (finite- element solid mechanics) and associated software details are described. For solid mechanics demonstration problems the method shows order of magnitudes improvement in accuracy over traditional stochastic collocation. For the reacting hypersonics problem, the method is implemented as a streamline integration and results show very good accuracy for the approximate sample solutions of re-entry flow past the Apollo capsule geometry at Mach 30.

42 ENGINEERING↗
DOE OSTI2024

Surface Pressure Fluctuations Induced by a Hypersonic Turbulent Boundary Layer on a Sharp Cone at Angle of Attack

High-fidelity simulations are performed to characterize the turbulence-induced wall pressure fluctuations on a sharp cone at a 5.5° angle-of-attack in a Mach 8 flow. Wall-resolved large-eddy simulation (LES) and wall-modeled large-eddy simulation (WMLES) results are compared to measurements at several locations on the cone body, where comparisons with high-frequency PCB sensors are good, while comparison with low-frequency kulite sensors varies depending on the location. For a given streamwise location, simulation results show significant azimuthal variation in the wall pressure fluctuation statistics. Comparisons between LES and WMLES results indicate that WMLES can accurately reproduce the low-frequency component of the autospectra quite well, while slight sensitivity is apparent in the coherence. A grid sensitivity study for WMLES further demonstrated sensitivity in the wall pressure fluctuation coherence, but not in the autospectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Expanded verification and validation studies of hypersonic aerodynamics with multiple physics-fidelity models

Hypersonic aerothermodynamics is an important domain of modern multiphysics simulation. The Multi-Fidelity Toolkit is a simulation tool being developed at Sandia National Laboratories to predict aerodynamic properties for compressible flows from a range of physics fidelities and computational speeds. These models include the Reynolds-averaged Navier–Stokes (RANS) equations, the Euler equations with momentum-energy integral technique (MEIT), and modified Newtonian aerodynamics with flat-plate boundary layer (MNA+FPBL) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities, as well as model validation with uncertainty considerations. Code verification activities on the MNA+FPBL model build on previous work by focusing on the viscous portion of the model. Viscous quantities of interest are compared against those from an analytical solution for flat-plate, inclined-plate, and cone geometries. The code verification methodology for the MEIT model is also presented. Test setup and results of code verification tests on the laminar and turbulent models within MEIT are shown. Solution-verification activities include grid-refinement studies on simulations that model the HIFiRE-1 wind tunnel experiments. These experiments are used for validation of all model fidelities. A thorough validation comparison with prediction error and uncertainty is also presented. Three additional HIFiRE-1 experimental runs are simulated in this study, and the solution verification and validation work examines the effects of the associated parameter changes on model performance. Finally, a study is presented that compares the computational costs and fidelities from each of the different models.

42 ENGINEERING↗
DOE OSTI2021

HTR-1.3 solver: Predicting electrified combustion using the hypersonic task-based research solver

Here this manuscript presents an updated open-source version of the Hypersonics Task-based Research (HTR) solver. The solver, whose main features are presented in Di Renzo et al. (2020) and Di Renzo & Pirozzoli (2021), is designed for direct numerical simulation of reacting flows at high Reynolds numbers. This new version extends the applications of the HTR solver to turbulent combustion in the presence of external electric fields. In particular, a new distributed Poisson solver compatible with heterogeneous architectures has been incorporated in the algorithm to compute the electric potential distribution in bi-periodic configurations. The drift fluxes of the electrically charged species are now included in the transport equations using a targeted essentially non-oscillatory scheme. A verification of these new features of the solver is provided using one-dimensional burner stabilized flames, whereas a three dimensional turbulent flame is utilized to discuss the scalability of the proposed numerical tool.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Time-resolved unsteady shock motion and spectral analysis of a hypersonic shock-wave/boundary-layer interaction over a double-cone

Direct experimental measurements of double-cone shock-wave/boundary-layer interaction dynamics at Mach 7.2 were collected using a schlieren-based shock-detection algorithm at the University of Texas at San Antonio. The points of flow separation and reattachment were spatially and spectrally evaluated at 300 kHz for time-resolved characterization of the unsteadiness of this flow field. The separation bubble was found to rapidly expand and contract about a large mean length. Spectral analysis of the velocity fluctuations of this motion revealed this behavior to be a turbulent phenomenon, with spectral energy trends consistent with the Kolmogorov theory for the turbulent energy cascade.

Glasby, Ryan [ORNL]↗
DOE OSTI2024

Hypersonic Wake Velocity Measurements Using Acetone Molecular Tagging Velocimetry

A non-intrusive laser diagnostic known as molecular tagging velocimetry was used to find quantitative off-body velocity measurements in the wake of a sphere in the Mach 7 Ludwieg Tube Wind Tunnel located at The University of Texas at San Antonio. Acetone gas seeded in the flow was excited using the 4th harmonic of a pulse-burst Nd:YAG laser. The experimental results were used to validate and compare to simulations using a continuous Galerkin flow solver. Both the experimental and simulated results agreed on negative velocities in the viscous shear layer. Near the centerline of the sphere, average velocities of -200 m/s were observed due to reverse flow in the recirculation region. Outside of the shear layer, velocities of 800 m/s were observed in both the experimental and simulated results. The overall average uncertainty for the strut-mounted and free-flight case was estimated to be ±27 m/s. The flow features of the wake were visualized using schlieren imaging, and the experimental results agreed well with the simulated results.

Shoemake, Lawton↗
DOE OSTI2024

Thermochemical Nonequilibrium Modeling in a Continuous-Galerkin, Finite-Element Framework

The presented work discusses the implementation, verification, and validation of Park's two-temperature model in a scalable, computational fluid dynamics (CFD) code developed at the US Department of Energy's Oak Ridge National Laboratory (ORNL). The implementation of Park's two-temperature model was verified through 0D test cases involving an adiabatic reactor and a nitrogen thermal bath. The implementation was then validated through comparisons with other validated CFD codes and experimental data on a hypersonic cylinder and double cones. These are standard benchmark test cases for thermochemical non-equilibrium (TCNE) modeling, and all data are shared publicly. The verification and validation results showed that ORNL's in-house CFD code could model complex, high-speed flow problems with and without TCNE modeling. This work is essential for future research involving 3D shock wave/boundary layer interactions (SBLIs).

Nutter, Nicole↗
DOE OSTI2020

Laser-Plasma Spatiotemporal Cyanide Spectroscopy and Applications

This article reports new measurements of laser-induced plasma hypersonic expansion measurements of diatomic molecular cyanide (CN). Focused, high-peak-power 1064 nm Q-switched radiation of the order of 1 TW/cm 2 generated optical breakdown plasma in a cell containing a 1:1 molar gas mixture of N 2 and CO 2 at a fixed pressure of 1.1 × 10 5 Pascal and in a 100 mL/min flow of the mixture. Line-of-sight (LOS) analysis of recorded molecular spectra indicated the outgoing shockwave at expansion speeds well in excess of Mach 5. Spectra of atomic carbon confirmed increased electron density near the shockwave, and, equally, molecular CN spectra revealed higher excitation temperature near the shockwave. Results were consistent with corresponding high-speed shadowgraphs obtained by visualization with an effective shutter speed of 5 nanoseconds. In addition, LOS analysis and the application of integral inversion techniques allow inferences about the spatiotemporal plasma distribution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Surface Properties of Double-Fin Generated Shock-Wave/Boundary-Layer Interactions

Shock-Wave/Boundary-Layer Interaction (SBLI) is ubiquitous in high-speed air vehicles. SBLI surface flowfield properties of crossing shocks induced by symmetric Double-Fins (DFs) at fin angles of 8° and 10° at Mach 2 are experimentally investigated; a 10° Single-Fin (SF) is also explored as a comparison. Surface visualization is used to capture distinct flowfield structures such as separation and upstream influence. Steady and unsteady Pressure-Sensitive Paints (PSPs) are used to obtain global surface pressure fields. Characteristics of DF-SBLI surface topology are discussed by comparing the corresponding SF-SBLI. Here, the effect of fin angle on DF-SBLI is characterized, and unsteady dynamics of the pressure field are examined. Spectral Proper Orthogonal Decomposition (SPOD) is performed on unsteady PSP, and traveling surface pressure waves and their frequency-dependent behavior are discussed. Dispersion relation of surface pressure waves is revealed by performing two-dimensional space–time Fourier transform on unsteady PSP, and the results are compared with the phase velocity distributions obtained from SPOD. Accompanying group velocity behavior is examined for distinct SBLI regions wherein wave propagation is found dependent on its frequency. Steady surface pressures along the centerline are compared with previous experimental data. Fin-tip spacing is found to be an essential parameter for scaling of upstream influence length.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Hypersonic Jets of Detonation Products in the Hydrodynamic Collapse of Macroscopic Voids

Localizing the energetic output from detonation waves has been a long-standing challenge in applied detonation physics. Here, energy localization is achieved via machined millimeter scale voids in pressed samples of PBX 9501, an HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane)-based plastic bonded explosive. A main mechanism of energy localization in these systems, the formation of hydrodynamic jets of dense product gases, is characterized experimentally using a semicylindrical geometry in witness plate impact experiments and streak imaging of the jet propagating into the air. The distance at which the jet is optimally developed is identified and the supersonic flow structure in the vicinity of this feature is explored using hydrocode simulations. This analysis found that most of the kinetic energy of the hydrodynamic jet arises from pressure gradients induced by geometrically mediated squeeze flow lateral to the direction of detonation propagation. This work presents a new development in the control of energetic output from detonation waves and applications to detonation wave shaping are discussed.

42 ENGINEERING↗
DOE OSTI2021

Hybrid Large Eddy Simulation/Reynolds-Averaged Navier–Stokes Analysis of a Premixed Ethylene-Fueled Dual-Mode Scramjet Combustor

Hydrocarbon fuels offer optimal high energy per volume for scramjet applications for sustained hypersonic flight but require additional residence time due to slower ignition delays (compared with hydrogen fuel). The injection of ethylene at the start of the isolator of a dual-mode scramjet combustor, operating in ramjet mode, allows sufficient mixing to achieve efficient premixed turbulent combustion. A cavity flameholder anchors the flame, supplying sufficient radicals to sustain a stable flame in the high-speed environment. This work investigates flame structure and stabilization limits of a new configuration with a scaled-down cavity embedded in the flow path of the combustor with a strut and insert. The cavity is reduced in size by one-third to enable eventual direct numerical simulations of the flame stabilization process. This work, however, focuses on modeling the full isolator/combustor geometry using a hybrid large eddy simulation/Reynolds-averaged Navier–Stokes simulation strategy. Additionally, particle image velocimetry, planar laser-induced fluorescence, coherent anti-Stokes Raman spectroscopy, and pressure measurements are compared with numerical predictions to analyze and characterize the conditions within the combustor, including flame structure, flow velocities, species composition, and wall pressure. An adjustable air throttle in the extender is able to control the placement of the shock train in the isolator and maintain a stable flame at various equivalence ratios. The simulations show reasonably good agreement with the experimental scalar and velocity data and predict a flame angle consistent with premixed turbulent flame-speed correlations.

42 ENGINEERING↗
DOE OSTI2024

Inverse-Velocity Transformation Wall Model for Reacting Turbulent Hypersonic Boundary Layers

The present study builds on prior work by taking advantage of the novel framework proposed by Griffin et al. (hereafter referred to as the GFM) as a baseline. The model is progressively extended to multicomponent reacting mixtures, accounting for differential diffusion and finite-rate chemistry in a similar fashion to Di Renzo & Urzay and Di Renzo et al. The accuracy of the present approach, as well the prior model of Di Renzo & Urzay, is assessed in an a priori sense for the first time in a turbulent reacting boundary layer, using the boundary-layer data of Williams et al. Five species are included in the present analysis, i.e., N s = 5, namely N 2 , O 2 , NO, N and O, a neutral mixture most representative of dissociation/recombination phenomena for temperatures below 6000 K. A schematic of a flow over a wedge representative of the described configuration is presented in Figure 1. The brief is organized as follows: In Section 2, the wall-model equations and the computational framework are presented. In Section 3, the a priori results of the proposed model are described and compared to the extended EWM. Finally, in Section 4, some conclusions are offered.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The Computational Analysis of a Generic Hypersonic Store Separation

Many current hypersonic vehicles involve a main vehicle from which components separate during flight and presents simulation challenges such as shock-shock interactions, highly nonlinear interactions, and shock wave boundary layer interactions. These aerodynamic phenomena influence the attitude and trajectory of each object, which needs to be understood before large-scale experiments can be run. To gain a generic understanding of the process, high-fidelity Reynolds-averaged Navier-Stokes solutions are performed on a conical vehicle geometry passing through an oblique shock wave that is representative of different configurations and scenarios. Trajectory and applied forces are tracked and show the vehicle’s dynamics are predominantly the result of differential flow incidence angles causing a strong shock; this leads to a large pressure increase over a fraction of the vehicle, which influences pitch. The vehicle appears to follow conventional stability theory with detached eddy simulation and small disturbances in initial attitude shown to have minimal influence on the scenario.

Savery, Ryan↗