Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock dynamics”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Accuracy of the Burnett equations for hypersonic real gas flows

Burnett solutions for the structure of normal shock waves for Mach numbers ranging from 1.1 to 50 have been obtained by using the time-dependent equation set in conjunction with modern computational fluid dynamic techniques. The current work seeks to provide corrected and improved high-resolution direct simulation Monte Carlo (DSMC) results for hypersonic shocks in monatomic gases. These DSMC results will serve as a reliable test base for assessing the accuracy of various continuum equations more advanced than Navier-Stokes, such as the Burnett equations. The main conclusion of Fisko and Chapman (1988) that the Burnett equations are more accurate than the Navier-Stokes equations for all gases and at all Mach numbers in matching DSMC shock wave results is reassesed by comparing the improved and corrected DSMC results with continuum Burnett results. Various Burnett results for shocks in a diatomic gas, nitrogen, which include improvements to the work of Lumpkin and Chapman (1989) are presented. It is concluded that the Burnett equations represent a significant improvement over Navier-Stokes equations for flows exhibiting translational nonequilibrium.

Lumpkin, Forrest E., III↗
NASA NTRS1981

Shock waves, atmospheric structure and mass loss in Miras

Large amplitude shock waves are observed to be present in the atmospheres of the Mira variables: spectral line doubling with Delta v 30 km/s is present in infrared spectra. Even the visible spectra contain some evidence for such shocks. These shocks are sufficiently large to clearly dominate the energy balance of the atmosphere. Mira variables also show symptoms of substantial mass loss rates: they are strong maser and infrared continuum sources and have strong circumstellar absorption features. The pulsation induced shocks which are seen to be present are obvious suspects for causing or enhancing the mass loss from these stars. The Miras thus present an ideal case for the study of dynamical effects on atmospheric structure, since both the dynamics and the results are clearly observable. The results are given of calculations of the thermalization and cooling of the shock heated material passing through shock fronts whose properties were selected to be consistent with both the isothermal models and the spectroscopic observations.

Willson, L. A.↗
DOE OSTI2022

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗
NASA NTRS1993

Spectropolarimetry of the post-main-sequence bipolar nebulae GL 618, M2-56, and M1-92

New high-quality spectropolarimetry of the post-main-sequence bipolar nebulae GL 618, M2-56, and M1-92 is presented which permits accurate separation of the scattered and unscattered components. Shock emission dominates the optical line spectrum of the three nebulae and probably plays a significant role in their dynamical evolution. The central H II region spectrum for GL 6189 is isolated and T(stellar) of 36,000-40,000 and log (N/O) = 0.0 +/- 0.2 are derived. GL 618 and M2-56 have shock velocities of 40-60 km/s; that of M1-92 is 60-100 km/s. Log (N/O) is derived for different regions of the three nebulae. There is an apparent correlation of outflow velocity with chemical abundance.

Trammell, Susan R.↗
DOE OSTI2022

The role of pre-existing heterogeneities in materials under shock and spall

There has been a challenge for many decades to understand how heterogeneities influence the behavior of materials under shock loading, eventually leading to spall formation and failure. Experimental, analytical, and computational techniques have matured to the point where systematic studies of materials with complex microstructures under shock loading and the associated failure mechanisms are feasible. This is enabled by more accurate diagnostics as well as characterization methods. As interest in complex materials grows, understanding and predicting the role of heterogeneities in determining the dynamic behavior becomes crucial. Early computational studies, hydrocodes, in particular, historically preclude any irregularities in the form of defects and impurities in the material microstructure for the sake of simplification and to retain the hydrodynamic conservation equations. Contemporary computational methods, notably molecular dynamics simulations, can overcome this limitation by incorporating inhomogeneities albeit at a much lower length and time scale. This review discusses literature that has focused on investigating the role of various imperfections in the shock and spall behavior, emphasizing mainly heterogeneities such as second-phase particles, inclusions, and voids under both shock compression and release. Pre-existing defects are found in most engineering materials, ranging from thermodynamically necessary vacancies, to interstitial and dislocation, to microstructural features such as inclusions, second phase particles, voids, grain boundaries, and triple junctions. This literature review explores the interaction of these heterogeneities under shock loading during compression and release. Systematic characterization of material heterogeneities before and after shock loading, along with direct measurements of Hugoniot elastic limit and spall strength, allows for more generalized theories to be formulated. Further, continuous improvement toward time-resolved, in situ experimental data strengthens the ability to elucidate upon results gathered from simulations and analytical models, thus improving the overall ability to understand and predict how materials behave under dynamic loading.

36 MATERIALS SCIENCE↗
NASA NTRS1989

The magnetospheric response to 8-minute period strong-amplitude upstream pressure variations

This paper documents a series of brief, strong (delta p/p = 1), dynamic pressure oscillations that occurred in the region upstream of the earth's bow shock during a period of radial interplanetary magnetic field. The analyzed set of oscillations, which may be either intrinsic solar wind or bow shock-related phenomena, recur approximately every 8-10 min, and their magnetic field signatures occur nearly simultaneously over great distances transverse to the earth-sun line. The pressure oscillations appear to drive tailward-moving magnetopause surface wavelets. In turn, the surface wavelets can be identified as hydromagnetic waves with strong compressional components in the outer magnetosphere and as quasi-periodic variations in electron precipitation and high-latitude ground pulsations. Observations by spacecraft in the outer dayside magnetosphere are used to predict geosynchronous and subsolar magnetic field strengths, the location of the subsolar magnetopause, the solar wind dynamic pressure, and variations in the energetic magnetospheric ion flux.

Sibeck, D. G.↗
NASA NTRS2006

Notes on Earth Atmospheric Entry for Mars Sample Return Missions

The entry of sample return vehicles (SRVs) into the Earth's atmosphere is the subject of this document. The Earth entry environment for vehicles, or capsules, returning from the planet Mars is discussed along with the subjects of dynamics, aerodynamics, and heat transfer. The material presented is intended for engineers and scientists who do not have strong backgrounds in aerodynamics, aerothermodynamics and flight mechanics. The document is not intended to be comprehensive and some important topics are omitted. The topics considered in this document include basic principles of physics (fluid mechanics, dynamics and heat transfer), chemistry and engineering mechanics. These subjects include: a) fluid mechanics (aerodynamics, aerothermodynamics, compressible fluids, shock waves, boundary layers, and flow regimes from subsonic to hypervelocity; b) the Earth s atmosphere and gravity; c) thermal protection system design considerations; d) heat and mass transfer (convection, radiation, and ablation); e) flight mechanics (basic rigid body dynamics and stability); and f) flight- and ground-test requirements; and g) trajectory and flow simulation methods.

Rivell, Thomas↗
NASA NTRS2000

Screech Noise Generation From Supersonic Underexpanded Jets Investigated

Many supersonic military aircraft and some of the modern civilian aircraft (such as the Boeing 777) produce shock-associated noise. This noise is generated from the jet engine plume when the engine nozzle is operated beyond the subsonic operation limit to gain additional thrust. At these underexpanded conditions, a series of shock waves appear in the plume. The turbulent vortices present in the jet interact with the shock waves and produce the additional shock-associated noise. Screech belongs to this noise category, where sound is generated in single or multiple pure tones. The high dynamic load associated with screech can damage the tailplane. One purpose of this study at the NASA Glenn Research Center at Lewis Field was to provide an accurate data base for validating various computational fluid dynamics (CFD) codes. These codes will be used to predict the frequency and amplitude of screech tones. A second purpose was to advance the fundamental physical understanding of how shock-turbulence interactions generate sound. Previously, experiments on shock-turbulence interaction were impossible to perform because no suitable technique was available. As one part of this program, an optical Rayleigh-scattering measurement technique was devised to overcome this difficulty.

Panda, Jayanta↗
NASA NTRS2008

Estimating Shock Spectra: Extensions beyond GEVS

Shock response spectra (SRS) are the standard description of some vibration environments on spacecraft for equipment qualification. For shock events produced by pyrotechnic devices, SRS can have significant frequency content as high as 10 kHz. It is difficult to construct and analyze finite element models that can resolve dynamic behavior at such high frequencies. GEVS provides simple, empirically based methods for approximating the SRS for a wide variety of shock events. It begins with a base SRS according to the type of pyrotechnic device, and then provides attenuation relations to adjust this SRS according to distance from the shock source, the type of structural frame and the properties of any structural joints between the source and equipment. In our paper we extend GEVS to include more detailed information about the spacecraft structure. To retain the general framework of GEVS, we begin with a base SRS and adjust this SRS using attenuation relations. We use modal and traveling wave concepts to derive the attenuation relations for simple canonical structures. Then we show how these concepts can be used to analyze more complex structures using finite element mode shapes to explicitly calculate the attenuation factors. Since the low- to mid-frequency finite element modal information is extrapolated to obtain the low- to high-frequency attenuation relations, the resulting attenuated SRS is formulated as an upper bound rather than as mean predicted values. We illustrate the extended GEVS approach by analyzing the impact response of composite tubes and the shock response of the STEREO spacecraft.

Igusa, Takeru↗
DOE OSTI2023

Direct imaging of shock wave splitting in diamond at Mbar pressure

Understanding the behavior of matter at extreme pressures of the order of a megabar (Mbar) is essential to gain insight into various physical phenomena at macroscales—the formation of planets, young stars, and the cores of super-Earths, and at microscales—damage to ceramic materials and high-pressure plastic transformation and phase transitions in solids. Under dynamic compression of solids up to Mbar pressures, even a solid with high strength exhibits plastic properties, causing the induced shock wave to split in two: an elastic precursor and a plastic shock wave. This phenomenon is described by theoretical models based on indirect measurements of material response. The advent of x-ray free-electron lasers (XFELs) has made it possible to use their ultrashort pulses for direct observations of the propagation of shock waves in solid materials by the method of phase-contrast radiography. However, there is still a lack of comprehensive data for verification of theoretical models of different solids. Here, we present the results of an experiment in which the evolution of the coupled elastic–plastic wave structure in diamond was directly observed and studied with submicrometer spatial resolution, using the unique capabilities of the x-ray free-electron laser (XFEL). The direct measurements allowed, for the first time, the fitting and validation of the 2D failure model for diamond in the range of several Mbar. Our experimental approach opens new possibilities for the direct verification and construction of equations of state of matter in the ultra-high-stress range, which are relevant to solving a variety of problems in high-energy-density physics.

36 MATERIALS SCIENCE↗
NASA NTRS1978

Radiative cooling in shock-heated hydrogen-helium plasmas

Axial and off-axis radiative cooling of cylindrical shock-heated hydrogen-helium plasmas is investigated theoretically and experimentally. The coupled fluid dynamic-radiative transfer equations are solved by a combination of approximation techniques aimed at simplifying the computation of the flux divergence term, namely, the quasi-isothermal approximation and the exponential approximation developed for the solid angle integration. The accuracy of the approximation schemes has been assessed and found acceptable for applying the methods to the rapid computation of the radiatively coupled flow problem. Radiative cooling experiments were conducted in a 6-inch annular arc accelerator shock tube (ANAA) for an initial pressure of 1 torr and shock speeds from 35 to 45 Km/sec. The results indicate that the lateral cooling is small compared with the axial cooling, and that better agreement is achieved between the data and the theoretical results by inclusion of the lateral temperature gradient.

Poon, P. T. Y.↗
NASA NTRS1981

Shock-associated low-energy ion enhancements observed by Voyagers 1 and 2

Observations of shock-associated ion enhancements at energies of not less than 30 keV are presented from data gathered by Voyagers 1 and 2. Observations include examples of energetic storm particle events associated with flare-produced shocks and examples of corotating particle events (CPE) associated with forward and reverse shocks that bound corotating interaction regions in the outer heliosphere. Most of the CPE have recurrent double-peaked intensity enhancements showing little or no velocity dispersion at peak intensities, time durations of several days, and soft energy spectra extending up to 5 MeV/nucleon. The lowest energy ion enhancements are confined mainly downstream of the corotating interaction regions with the magnitude and duration of the upstream enhancements increasing with increasing ion energy. The observations are consistent with the dynamical and kinematical effects expected when low energy ions are accelerated at quasi-perpendicular shocks.

Decker, R. B.↗
NASA NTRS2010

X-Ray Measured Dynamics of Tycho's Supernova Remnant

We present X-ray proper-motion measurements of the forward shock and reverse-shocked ejecta in Tycho's supernova remnant, based on three sets of archival Chandra data taken in 2000, 2003, and 2007. We find that the proper motion of the edge of the remnant (i.e., the forward shock and protruding ejecta knots) varies from 0.''20 yr-1 (expansion index m = 0.33, where R = tm ) to 0.''40 yr-1 (m = 0.65) with azimuthal angle in 2000-2007 measurements, and 0.''14 yr-1 (m = 0.26) to 0.''40 yr-1 (m = 0.65) in 2003-2007 measurements. The azimuthal variation of the proper motion and the average expansion index of [approx]0.5 are consistent with those derived from radio observations. We also find proper motion and expansion index of the reverse-shocked ejecta to be 0.''21-0.''31 yr-1 and 0.43-0.64, respectively. From a comparison of the measured m-value with Type Ia supernova evolutionary models, we find a pre-shock ambient density around the remnant of [less, similar]0.2 cm-3.

Katsuda, Satoru↗
DOE OSTI2021

Controlled Pyroshock Test Transients Can Be Used To Better Match Operational Shock Transients and Reduce Artificial Pyroshock Test Failures [Slides]

Shock analysis is typically done via SRS (shock response spectrum). Oftentimes shock analysis is not performed because most loads specifications only supply SRS requirements with no representative shock transients that can be used for analysis. SRS analysis is too conservative. There is no structural dynamics involved as all modes are assumed to have peak responses simultaneously and in-phase with each other. It is possible that all other loads and combination of loads show the design to be good, whereas SRS analysis (1000 g’s are involved) results may require redesign. Multiple shock transients can satisfy a SRS, yet most loads documents don’t provide a representative shock transient that could be used for shock analysis. Shock testing is less common than harmonic and random vibration testing. Engineers typically go into shock tests with very little knowledge of what to expect. A pass/fail type of approach with very little engineering done to prepare for the test or to control the type of shock transient is used. Most engineers will not model the shock test fixture due to lack of knowledge about the details of the test setup. Shock analysis and testing is often an open loop type situation with no test setup modeling, no test correlation, and no control over the test shock transient. In this work, if we can create applied test force transient that comes close to matching calculated force transient then a more reasonable acceleration transient will be created at the test reference accelerometer that better matches operational shocks.

42 ENGINEERING↗
DOE OSTI2021

Equation of State Measurements on Iron Near the Melting Curve at Planetary Core Conditions by Shock and Ramp Compressions

The outer core of the Earth is composed primarily of liquid iron, and the inner core boundary is governed by the intersection of the melt line and the geotherm. While there are many studies on the thermodynamic equation of state for solid iron, the equation of state of liquid iron is relatively unexplored. In this work, we use dynamic compression to diagnose the high-pressure liquid equation of state of iron by utilizing the shock-ramp capability at Sandia National Laboratories’ Z-Machine. This technique enables measurements of material states off the Hugoniot by initially shocking samples and subsequently driving a further, shockless compression. Planetary studies benefit greatly from isentropic, off-Hugoniot experiments since they can cover pressure-temperature (P-T) conditions that are close to adiabatic profiles found in planetary interiors. We used this method to drive iron to P-T conditions similar to those of the Earth’s outer-inner core boundary, along an elevated-temperature isentrope in the liquid from 275 GPa to 400 GPa. We derive the equation of state using a hybrid backward integration – forward Lagrangian technique on particle velocity traces to determine the pressure-density history of the sample. Our results are in excellent agreement with SESAME 92141, a previously published equation of state table. With our data and previous experimental data on liquid iron we provide new information on the iron melting line and derive new parameters for a Vinet-based equation of state. The table and our parameterized equation of state are applied to provide an updated means of modeling the pressure, mass, and density of liquid iron cores in exoplanetary interiors.

58 GEOSCIENCES↗
NASA NTRS1989

Direct simulation of compressible turbulence

Several direct simulations of 3-D homogeneous, compressible turbulence are presented with emphasis on the differences with incompressible turbulent simulations. A fully spectral collocation algorithm, periodic in all directions coupled with a 3rd order Runge-Kutta time discretization scheme is sufficient to produce well-resolved flows at Taylor Reynolds numbers below 40 on grids of 128x128x128. A Helmholtz decomposition of velocity is useful to differentiate between the purely compressible effects and those effects solely due to vorticity production. In the context of homogeneous flows, this decomposition in unique. Time-dependent energy and dissipation spectra of the compressible and solenoidal velocity components indicate the presence of localized small scale structures. These structures are strongly a function of the initial conditions. Researchers concentrate on a regime characterized by very small fluctuating Mach numbers Ma (on the order of 0.03) and density and temperature fluctuations much greater than sq Ma. This leads to a state in which more than 70 percent of the kinetic energy is contained in the so-called compressible component of the velocity. Furthermore, these conditions lead to the formation of curved weak shocks (or shocklets) which travel at approximately the sound speed across the physical domain. Various terms in the vorticity and divergence of velocity production equations are plotted versus time to gain some understanding of how small scales are actually formed. Possible links with Burger turbulence are examined. To visualize better the dynamics of the flow, new graphic visualization techniques have been developed. The 3-D structure of the shocks are visualized with the help of volume rendering algorithms developed in-house. A combination of stereographic projection and animation greatly increase the number of visual cues necessary to properly interpret the complex flow.

Zang, T. A.↗
DOE OSTI2023

Shock compression of poly(methyl methacrylate) PMMA in the 1000 GPa regime: Z machine experiments

Hydrocarbon polymers are used in a wide variety of practical applications. In the field of dynamic compression at extreme pressures, these polymers are used at several high energy density (HED) experimental facilities. One of the most common polymers is poly(methyl methacrylate) or PMMA, also called Plexiglass® or Lucite®. Here, we present high-fidelity, hundreds of GPa range experimental shock compression data measured on Sandia’s Z machine. We extend the principal shock Hugoniot for PMMA to more than threefold compression up to 650 GPa and re-shock Hugoniot states up to 1020 GPa in an off-Hugoniot regime, where experimental data are even sparser. These data can be used to put additional constraints on tabular equation of state (EOS) models. The present results provide clear evidence for the need to re-examine the existing tabular EOS models for PMMA above ~120 GPa as well as perhaps revisit EOSs of similar hydrocarbon polymers commonly used in HED experiments investigating dynamic compression, hydrodynamics, or inertial confinement fusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1981

Unsteady environments and responses of the shuttle combined loads orbiter test

Both separate and combined wind tunnel and vibration shaker tests were conducted on two structural panels representative of the Shuttle orbiter in the NASA LaRC 8-foot transonic pressure tunnel (TPT) to determine the effects of combined loads on the thermal protection system (TPS). The primary objective of this test was to provide a combined full-scale load environment and realistic time history of the dynamic pressures, Mach numbers (through transonic), and dynamic structural responses of these panels. The panels were selected from orbiter locations where interactive load sources such as aerodynamic shock waves, turbulent boundary layers, strut-induced vorticity, and substrate deformation combined to provide high bonding loads between the TPS and the orbiter structure. The test panels for this program were highly instrumented with static and dynamic pressure gages, accelerometers, deflectometers, strain gages, Schileren and high-speed photography, and special instrumentation necessary to determine TPS/structure interface loads and tile motions. Two test specimens of each orbiter panel were utilized. Both were high-fidelity representations of the selected orbiter location. The first panel (calibration panel) for each test was utilized to calibrate the tunnel and shaker system as well as provide a data-gathering source for applied static and dynamic loads. The second panel (test panel) for each test was subjected to flight time exposures only so that TPS characteristics as a function of flight mission exposure could be determined. These tests were the last in a series that was suce6ssfully completed and necessary to certify the TPS prior to the first launching of the Shuttle.

P. H. Schuetz↗