Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock wave”

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 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

The Origin of Inlet Buzz in a Mach 1.7 Low Boom Inlet Design

Supersonic inlets with external compression, having a good level performance at the critical operating point, exhibit a marked instability of the flow in some subcritical operation below a critical value of the capture mass flow ratio. This takes the form of severe oscillations of the shock system, commonly known as "buzz". The underlying purpose of this study is to indicate how Detached Eddy Simulation (DES) analysis of supersonic inlets will alter how we envision unsteady inlet aerodynamics, particularly inlet buzz. Presented in this paper is a discussion regarding the physical explanation underlying inlet buzz as indicated by DES analysis. It is the normal shock wave boundary layer separation along the spike surface which reduces the capture mass flow that is the controlling mechanism which determines the onset of inlet buzz, and it is the aerodynamic characteristics of a choked nozzle that provide the feedback mechanism that sustains the buzz cycle by imposing a fixed mean corrected inlet weight flow. Comparisons between the DES analysis of the Lockheed Martin Corporation (LMCO) N+2 inlet and schlieren photographs taken during the test of the Gulfstream Large Scale Low Boom (LSLB) inlet in the NASA 8x6 ft. Supersonic Wind Tunnel (SWT) show a strong similarity both in turbulent flow field structure and shock wave formation during the buzz cycle. This demonstrates the value of DES analysis for the design and understanding of supersonic inlets.

Mach 1.7↗
NASA NTRSDate not supplied

Development and Assessment of a New Particle Image Velocimetry System in the NASA GRC 225 cm2 Wind Tunnel

As part of an ongoing campaign to provide detailed and thorough measurements in a two-dimensional impinging shock-wave/boundary-layer interaction (SWBLI), a custom implementation of the Particle Image Velocimetry (PIV) measurement technique in the NASA Glenn 225 cm2 Wind Tunnel was constructed. The flow seeding apparatus was devised to locally seed the region of interest in the overall flowfield to avoid coating the viewing window with seed material. The ability of the PIV system to provide small particles which faithfully track the underlying gas-phase flowfield was of key importance for the success of this system to make quality measurements in the SWBLI flowfield. In order to optimize and assess the performance of the particle seeding system, PIV measurements were made across a planar oblique shock wave and particle response was analyzed to quantify the particle lag effects due to the finite inertia of the seed droplets. Video recordings from a camera placed in the wind tunnel plenum tank were also used to discern the near-field behavior of the seeded plume and guide the design evolution of the seeding apparatus. The resulting optimized seeding configuration was applied to make PIV measurements in the wind tunnel test section boundary layer to assess the seed particle spatial coverage and concentration. This sample boundary layer dataset provides an analog to the eventual boundary layer measurements to be made in the axisymmetric test section under much more challenged viewing conditions.

PIV↗
NASA NTRSDate not supplied

Development and Assessment of a New Particle Image Velocimetry System in the NASA GRC 225 cm2 Wind Tunnel

As part of an ongoing campaign to provide detailed and thorough measurements in a two-dimensional impinging shock-wave/boundary-layer interaction (SWBLI), a custom implementation of the Particle Image Velocimetry (PIV) measurement technique in the NASA Glenn 225 cm2 Wind Tunnel was constructed. The flow seeding apparatus was devised to locally seed the region of interest in the overall flowfield to avoid coating the viewing window with seed material. The ability of the PIV system to provide small particles which faithfully track the underlying gas-phase flowfield was of key importance for the success of this system to make quality measurements in the SWBLI flowfield. In order to optimize and assess the performance of the particle seeding system, PIV measurements were made across a planar oblique shock wave and particle response was analyzed to quantify the particle lag effects due to the finite inertia of the seed droplets. Video recordings from a camera placed in the wind tunnel plenum tank were also used to discern the near-field behavior of the seeded plume and guide the design evolution of the seeding apparatus. The resulting optimized seeding configuration was applied to make PIV measurements in the wind tunnel test section boundary layer to assess the seed particle spatial coverage and concentration. This sample boundary layer dataset provides an analog to the eventual boundary layer measurements to be made in the axisymmetric test section under much more challenged viewing conditions.

PIV↗
DOE OSTI2022

Explosion behavior investigation and safety assessment of large-format lithium-ion pouch cells

We report large-format lithium-ion (Li-ion) batteries with high energy density for electric vehicles are prone to thermal runaway (or even explosion) under abusive conditions. In this study, overcharge induced explosion behaviors of large-format Li-ion pouch cells with Li[Ni 0.8 Co 0.1 Mn 0.1 ]O 2 cathode at different current rates (C-rates) (0.5C, 1C, 2C) were investigated. The explosion characteristics of the cells were elucidated by discussing the evolution of the cell voltage, the surface temperature and the shock wave pressure. Generally, the whole overcharge process could be divided into four stages according to the evolution of several key parameters and the overcharge behaviors; the overcharge C-rate has a great influence on cells’ thermal behaviors. The experimental results showed that the thermal runaway process of Li-ion cells caused by overcharging consisted of two kinds of explosions, physical explosion and chemical explosion. The existence of observable negative pressure zone in the pressure curves indicated that the Li-ion cells are not a self-supplying oxygen system during the explosion. Further, the explosion dynamics parameters were matched. An explosion TNT-equivalent conversion strategy that depended on the pressure of the shock wave was utilized to evaluate the released energy and its hazards. In addition, with respect to the overcharge of Li-ion pouch cells, a safety assessment method and a safety management method were proposed based on the explosion behaviors. From the perspective of battery safety, this study is of great significance for the safety design of Li-ion cells and can provide guidance for engineers to optimize the safety function of battery packs.

25 ENERGY STORAGE↗
NASA NTRS1985

Shock Compression of Liquid Helium to 56 GPa (560) Kbar

Shock-wave data are presented for liquid helium which has been compressed to densities up to five times greater than the normal liquid. The helium was heated to temperatures up to 21,000 K, while the maximum pressure attained was 56 GPa. The properties of helium and hydrogen are important for modeling the giant planets Saturn and Jupiter where these elements are the major constituents. Conditions on Saturn are of particular interest because studies have suggested that this planet has an internal energy source which is associated with unmixing and gravitational separation the hydrogen-helium fluid at pressures below 1 TPa. The existence of this phase transition depends very sensitively on the hydrogen and helium equation of state. In the experiments, strong shock waves were generated by the impact of planar projectiles into cryogenic specimen holders.

Nellis, W. J.↗
NASA NTRS1994

The gas-phase iron abundance in Herbig-Haro objects

The gas-phase abundance ratios Fe/S and Fe/O have been determined for the Herbig-Haro objects HH 1, HH 7, HH 11, HH 43A, and 'Burnham's Nebula' (HH 255). It is the purpose of this study to decide whether a sizeable fraction of the Fe in these HH objects is still bound in dust grains or whether the observed matter has gone through sufficiently fast shock waves so that the dust grains have been essentially destroyed and most of the iron has gone back into the gas phase. We have determined the abundance ratios using statistical equilibrium calculations for the ions Fe(+), S(+), and O(+). (These are the most abundant ions of the elements in question.) Abundance determinations have been made using homogeneous models of the HH objects for which electron temperatures and densities have been determined observationally from forbidden line ratios. The results show that the Fe/S ratio in the objects HH 1, HH 7, HH 11, and HH 43A agrees very well with the Population I abundance ratio. Only Burnham's Nebula (HH255) shows an Fe/S ratio which is about three times lower indicating a shock-wave history which is quite different from that of the other HH objects.

Beck-Winchatz, B.↗
NASA NTRS1984

Shock compression of liquid helium to 56 GPa (560 kbar)

Shock-wave data are presented for liquid helium which has been compressed to densities up to five times greater than the normal liquid. The helium was heated to temperatures up to 21,000 K, while the maximum pressure attained was 56 GPa. The properties of helium and hydrogen are important for modeling the giant planets Saturn and Jupiter where these elements are the major constituents. Conditions on Saturn are of particular interest because studies have suggested that this planet has an internal energy source which is associated with unmixing and gravitational separation of the hydrogen-helium fluid at pressures below 1 TPa. The existence of this phase transition depends very sensitively on the hydrogen and helium equation of state. In the experiments, strong shock waves were generated by the impact of planar projectiles into cryogenic specimen holders.

Nellis, W. J.↗
NASA NTRS1990

Global convergence of inexact Newton methods for transonic flow

In computational fluid dynamics, nonlinear differential equations are essential to represent important effects such as shock waves in transonic flow. Discretized versions of these nonlinear equations are solved using iterative methods. In this paper an inexact Newton method using the GMRES algorithm of Saad and Schultz is examined in the context of the full potential equation of aerodynamics. In this setting, reliable and efficient convergence of Newton methods is difficult to achieve. A poor initial solution guess often leads to divergence or very slow convergence. This paper examines several possible solutions to these problems, including a standard local damping strategy for Newton's method and two continuation methods, one of which utilizes interpolation from a coarse grid solution to obtain the initial guess on a finer grid. It is shown that the continuation methods can be used to augment the local damping strategy to achieve convergence for difficult transonic flow problems. These include simple wings with shock waves as well as problems involving engine power effects. These latter cases are modeled using the assumption that each exhaust plume is isentropic but has a different total pressure and/or temperature than the freestream.

Young, David P.↗
NASA NTRS1998

Low Dissipative High Order Shock-Capturing Methods using Characteristic-Based Filters

An approach which closely maintains the non-dissipative nature of classical fourth or higher- order spatial differencing away from shock waves and steep gradient regions while being capable of accurately capturing discontinuities, steep gradient and fine scale turbulent structures in a stable and efficient manner is described. The approach is a generalization of the method of Gustafsson and Olsson and the artificial compression method (ACM) of Harten. Spatially non-dissipative fourth or higher-order compact and non-compact spatial differencings are used as the base schemes. Instead of applying a scalar filter as in Gustafsson and Olsson, an ACM like term is used to signal the appropriate amount of second or third-order TVD or ENO types of characteristic based numerical dissipation. This term acts as a characteristic filter to minimize numerical dissipation for the overall scheme. For time-accurate computations, time discretizations with low dissipation are used. Numerical experiments on 2-D vortical flows, vortex-shock interactions and compressible spatially and temporally evolving mixing layers showed that the proposed schemes have the desired property with only a 10% increase in operations count over standard second-order TVD schemes. Aside from the ability to accurately capture shock-turbulence interaction flows, this approach is also capable of accurately preserving vortex convection. Higher accuracy is achieved with fewer grid points when compared to that of standard second-order TVD or ENO schemes. To demonstrate the applicability of these schemes in sustaining turbulence where shock waves are absent, a simulation of 3-D compressible turbulent channel flow in a small domain is conducted.

Yee, H. C.↗
NASA NTRS1998

Low Dissipative High Order Shock-Capturing Methods Using Characteristic-Based Filters

An approach which closely maintains the non-dissipative nature of classical fourth or higher- order spatial differencing away from shock waves and steep gradient regions while being capable of accurately capturing discontinuities, steep gradient and fine scale turbulent structures in a stable and efficient manner is described. The approach is a generalization of the method of Gustafsson and Oisson and the artificial compression method (ACM) of Harten. Spatially non-dissipative fourth or higher-order compact and non-compact spatial differencings are used as the base schemes. Instead of applying a scalar filter as in Gustafsson and Olsson, an ACM like term is used to signal the appropriate amount of second or third-order TVD or ENO types of characteristic based numerical dissipation. This term acts as a characteristic filter to minimize numerical dissipation for the overall scheme. For time-accurate computations, time discretizations with low dissipation are used. Numerical experiments on 2-D vortical flows, vortex-shock interactions and compressible spatially and temporally evolving mixing layers showed that the proposed schemes have the desired property with only a 10% increase in operations count over standard second-order TVD schemes. Aside from the ability to accurately capture shock-turbulence interaction flows, this approach is also capable of accurately preserving vortex convection. Higher accuracy is achieved with fewer grid points when compared to that of standard second-order TVD or ENO schemes. To demonstrate the applicability of these schemes in sustaining turbulence where shock waves are absent, a simulation of 3-D compressible turbulent channel flow in a small domain is conducted.

Yee, H. C.↗
DOE OSTI2020

An Eulerian multimaterial framework for simulating high-explosive aquarium tests

Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations. Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.

42 ENGINEERING↗
NASA NTRS2016

Background-Oriented Schlieren used in a hypersonic inlet test at NASA GRC

Background Oriented Schlieren (BOS) is a derivative of the classical schlieren technology, which is used to visualize density gradients, such as shock wave structures in a wind tunnel. Changes in refractive index resulting from density gradients cause light rays to bend, resulting in apparent motion of a random background pattern. The apparent motion of the pattern is determined using cross-correlation algorithms (between no-flow and with-flow image pairs) producing a schlieren-like image. One advantage of BOS is its simplified setup which enables a larger field-of-view (FOV) than traditional schlieren systems. In the present study, BOS was implemented into the Combined Cycle Engine Large-Scale Inlet Mode Transition Experiment (CCE LIMX) in the 10x10 Supersonic Wind Tunnel at NASA Glenn Research Center. The model hardware for the CCE LIMX accommodates a fully integrated turbine based combined cycle propulsion system. To date, inlet mode transition between turbine and ramjet operation has been successfully demonstrated. High-speed BOS was used to visualize the behavior of the flow structures shock waves during unsteady inlet unstarts, a phenomenon known as buzz. Transient video images of inlet buzz were recorded for both the ramjet flow path (high speed inlet) and turbine flow path (low speed inlet). To understand the stability limits of the inlet, operation was pushed to the point of unstart and buzz. BOS was implemented in order to view both inlets simultaneously, since the required FOV was beyond the capability of the current traditional schlieren system. An example of BOS data (Images 1-6) capturing inlet buzz are presented.

Background-oriented schlieren↗
NASA NTRS1991

A numerical study of shock wave/boundary layer interaction in nonequilibrium chemically reacting air - The effects of catalytic walls

This paper presents a numerical study that investigate the effects of nonequilibrium chemistry, and in particular, wall catalycity on the separated flow region generated by an oblique shock wave impinging upon a flat plate boundary layer. To obtain a solution to this problem, the full two dimensional Navier-Stokes equations were solved using MacCormack's predictor-corrector time dependent technique on a rectangular grid. Nonequilibrium chemistry was included by utilizing the 5 species, 17 reaction modified Dunn-Kang chemical kinetics model. Separate results were obtained for: calorically perfect, chemically reacting - noncatalytic wall, and chemically reacting - fully catalytic wall cases, for a given set of flow conditions. A direct comparison of all three cases revealed a slight decrease in the peak heat transfer for the noncatalytic wall case, as compared to the calorically perfect case. On the other hand, the fully catalytic wall case had a tremendous increase in the peak surface heat transfer. It is concluded that, for the particular conditons treated here (nearly frozen flow in the free stream), the effects of the nonequilibrium chemically reacting flow on the shock-wave/boundary-layer interaction depend critically on the catalycity of the wall, having virtually no effect for the case of a noncatalytic wall, and exerting a tremendous effect for a fully catalytic wall.

Grumet, Adam A.↗
NASA NTRS1985

Numerical simulation of compressible, turbulent, two-phase flow

A computer program for numerically simulating compressible, turbulent, two-phase flows is described and applied. Special attention is given to flows in which dust is ingested into the turbulent boundary layer behind shock waves moving over the earth's surface. it is assumed that the two phases are interpenetrating continua which are coupled by drag forces and heat transfer. The particle phase is assumed to be dilute, and turbulent effects are modeled by zero- and two-equation eddy viscosity models. An important feature of the turbulence modeling is the treatment of surface boundary conditions which control the ingestion of particles into the boundary layer by turbulent friction and diffusion. The numerical method uses second-order implicit upwind differencing of the inviscid terms of the equations and second-order central differencing of the viscous terms. A diagonal form of the implicit algorithm is used to improve efficiency, and the transformation to a curvilinear coordinate system is accomplished by the finite volume techniques. Applications to a series of representative flows include a two-phase nozzle flow, the steady flow of air over a sand bed, and the air flow behind a normal shock wave in uniform motion over a sand bed. Results of the latter two applications are compared with experimental results.

Coakley, t. J.↗
DOE OSTI2024

Helicity evolution at small x : quark to gluon and gluon to quark transition operators

We include the quark to gluon and gluon to quark shock-wave transition operators into the small Bjorken-x evolution equations for helicity in the flavor-singlet channel derived earlier in [1,2,3]. While such transitions do not affect the large-N c version of the evolution equations for helicity, the large-N c & N f equations are affected. (N c and N f are the numbers of quark colors and flavors, respectively.) We derive the corresponding corrected large-N c & N f equations for the polarized dipole amplitudes contributing to the flavor-singlet quark and gluon helicity distributions in the double-logarithmic approximation (DLA), resumming powers of α s ln 2 (1/x) with αs the strong coupling constant. We solve these equations iteratively and extract the polarized splitting functions up to four loops. We show that our splitting functions agree with the fixed-order perturbative calculations up to and including the existing three-loops results [4,5,6,7]. Similar to the large-Nc helicity evolution in the shock-wave approach [8], our large-N c & N f small-x splitting functions agree with those obtained in the infrared evolution equations framework from [9, 10] up to three loops, but appear to slightly disagree at four loops.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗