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.

SourceResearch connectionFrom yearThrough year

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

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Observations of twinning microstructure in iron ramp-compressed through the α – ε phase transition

Single crystal iron was ramp-compressed along the [100] axis at 4 GPa/ns through the α–ϵ phase transformation at 13 GPa. Unlike shock experiments, in these ramp compression experiments the bcc lattice is observed to be nearly isotropically relaxed before the phase transition. There was a mixed α/ϵ phase region starting at ~13 GPa which was observed for 0.75 ns until the peak stress of 18 GPa was reached. In situ x-ray diffraction measurements show the formation of new hcp orientations not reported in shock or quasistatic compression experiments. The new hcp orientations appear to be caused by {10¯11} <10¯1 ¯2> → {¯1011} <¯101¯2> sequential twins that occur during the phase transition. This twinning mechanism relieves the shear strain caused by the bcc-hcp phase transition from the isotropically relaxed bcc phase while subject to uniaxial compression. Furthermore, these results show that in iron the induced microstructure through a phase transition and the phase transition mechanism depend on the loading history.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1981

Fully implicit shock tracking

A tracking method for solving hyperbolic free boundary problems arising in fluid dynamics is presented. The principal feature of the new approach is that the implicit form is used to simultaneously treat both interior points and boundary conditions. In particular, the location of the free boundary (shock) surface is treated implicitly and coupled with all other unknowns. The method is presented in the context of unsteady one-dimensional flow in a variable area duct with an internal shock wave. The fully implicit method and other strategies for advancing the shock are compared for computing a steady solution via a time asymptotic approach. Issues regarding extension of the method to multiple dimensions are also discussed.

Bell, J. B.↗

NASA NTRS2014

Increased Mach Number Capability for the NASA Glenn 10x10 Supersonic Wind Tunnel

Computational simulations and wind tunnel testing were conducted to explore the operation of the Abe Silverstein Supersonic Wind Tunnel at the NASA Glenn Research Center at test section Mach numbers above the current limit of Mach 3.5. An increased Mach number would enhance the capability for testing of supersonic and hypersonic propulsion systems. The focus of the explorations was on understanding the flow within the second throat of the tunnel, which is downstream of the test section and is where the supersonic flow decelerates to subsonic flow. Methods of computational fluid dynamics (CFD) were applied to provide details of the shock boundary layer structure and to estimate losses in total pressure. The CFD simulations indicated that the tunnel could be operated up to Mach 4.0 if the minimum width of the second throat was made smaller than that used for previous operation of the tunnel. Wind tunnel testing was able to confirm such operation of the tunnel at Mach 3.6 and 3.7 before a hydraulic failure caused a stop to the testing. CFD simulations performed after the wind tunnel testing showed good agreement with test data consisting of static pressures along the ceiling of the second throat. The CFD analyses showed increased shockwave boundary layer interactions, which was also observed as increased unsteadiness of dynamic pressures collected in the wind tunnel testing.

Wind tunnel↗

NASA NTRS1984

Impact ejection, spallation, and the origin of meteorites

A model for the ejection of material from an impact crater which links ejection velocity, fragment size, and shock pressure through a simplified stress-wave propagation and reflection scheme is presented. It is shown that a small amount of material (0.01 to 0.05 projectile mass) may be ejected at high velocity without suffering petrologically detectable shock pressures. The largest fragments ejected at any velocity are spalls that originate from the target planet's surface. The spall size is proportional to the radius of the primary impactor and the target tensile strength and inversely proportional to ejection velocity. The shock level in the spalls is low, typically half of the dynamic crushing strength of the rock. The model also predicts the aspect ratio of the spalled fragments, the angle of ejection, and the sizes and shock level of other fragments originating deeper in the target. Comparison with observational and experimental data shows generally good agreement.

Melosh, H. J.↗

NASA NTRS1990

Assessment of shock capturing schemes for resonant flows in nonlinear instability analysis

The paper presents computational assessment of advanced numerical schemes for nonlinear acoustic problems related to combustion instabilities in liquid rocket engines. Several time-accurate, shock capturing schemes have been evaluated on a benchmark, closed-end resonant pipe flow problem. It involves the numerical solution of inviscid, compressible gas dynamics equations to predict acoustic wave propagation, wave steepening, formation of shocks, acoustic energy dissipation and wave-wall reflection for several hundred wave cycles. It was demonstrated that high accuracy TVD type schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time. The selected scheme was then applied to analyze the acoustic responses of resonant pipe-resonator, radial acoustic modes and hub-baffle configurations. Interesting observations of wave shape and damping characteristics have been drawn from presented computational studies.

Przekwas, A. J.↗

NASA NTRS2023

On-Board Deployment Event Verification for GOES-R Spacecraft

As is common with many spacecraft designs, the GOES-R vehicles require a series of deployment events to transition from the launch configuration to the operational configuration. Rather than implementing additional sensors to verify various deployments, the GOES-R program developed an alternate approach that uses existing gyro rate sensing. The approach includes two pieces: the first is a new onboard shock detection capability to confirm initiation of individual deployment events, and the second is a ground-based dynamics verification step to confirm completion of deployment events. We first present the algorithm that detects shock events for various deployment devices along with the parameter tuning performed during ground tests. We then show inflight performance of the shock detection algorithm. While shock detection is useful for observing initiation of deployment events, completion of some deployment events cannot be determined by shock detection alone. For these events, such as solar panel latch-up and deployable boom extension, the program developed dynamics models for the deployment transient responses. High-rate gyro data were recorded for these events, which allow the ground team to verify that the appendages were fully deployed. We show the predictive models for these events and corresponding flight results.

GeoXO↗

DOE OSTI2022

Characterization of a new ultra-high pressure shock tube facility for combustion and propulsion studies

A new shock tube facility has been designed, constructed, and characterized at the University of Central Florida. This facility is capable of withstanding pressures of up to 1000 atm, allowing for combustion diagnostics of extreme conditions, such as in rocket combustion chambers or in novel power conversion cycles. For studies with toxic gas impurities, the high initial pressures required the development of a gas delivery system to ensure the longevity of the facility and the safety of the personnel. Data acquisition and experimental propagation were implemented with remote access to ensure safety, paired with a LabVIEW- and Python-based user interface. Thus far, test pressures of 270 atm, blast pressures of 730 atm, and temperatures approaching 10 000 K have been achieved. The extreme limitations of this facility allow for emission spectroscopy to be performed during the oxidation of fuel mixtures, e.g., alkanes diluted in argon and carbon dioxide. Ignition delay times were determined and compared to simulations using chemical kinetic mechanisms. The design, experimental procedures, processes of analysis, and uncertainty determination are outlined, and typical pressure profiles are compared with a new gas dynamics solver and empirical correlations developed across multiple shock tube facilities. Preliminary reactive mixture analyses are included with further investigation of the mixtures outlined.

47 OTHER INSTRUMENTATION↗

NASA NTRS1984

Dynamic compression and volatile release of carbonates

Particle velocity profiles upon shock compression and isentropic releases were measured for polycrystalline calcite. The Solenhofen limestone release paths lie, close to the Hugoniot. Calcite 3 to 2 transition, upon release, was observed, but rarefaction shocks were not detected. The equation of state is used to predict the fraction of material devolatilized upon isentropic release as a function of shock pressure. The effect of ambient partial pressure of CO2 on the calculations is demonstrated and considered in models of atmospheric evolution by impact induced mineral devolatilization. The radiative characteristics of shocked calcite indicate that localization of thermal energy occurs under shock compression. Shock entropy calculations result in a minimum estimate of 90% devolatilization upon complete release from 10 GPa. Isentropic release paths from calculated continuum Hugoniot temperatures cross into the CaO (solid) + CO2 (vapor) field at improbably low pressures. It is found that release paths from measured shock temperatures cross into the melt plus vapor field at pressures greater than .5 GPa, which suggests that devolatilization is initiated at the shear banding sites.

Tyburczy, J. A.↗

NASA NTRS1986

Dynamic compression and volatile release of carbonates

Particle velocity profiles upon shock compression and isentropic releases were measured for polycrystalline calcite. The Solenhofen limestone release paths lie close to the Hugoniot. Calcite 3 to 2 transition, upon release was observed, but rarefaction shocks were not detected. The equation of state is used to predict the fraction of material devolatilized upon isentropic release as a function of shock pressure. The effect of ambient partial pressure of CO2 on the calculations is demonstrated and considered in models of atmospheric evolution by impact induced mineral devolatilization. The radiative characteristics of shocked calcite indicate that localization of thermal energy occurs under shock compression. Shock entropy calculations result in a minimum estimate of 90 percent devolatilization upon complete release from 10 GPa. Isentropic release paths from calculated continuum Hugoniot temperatures cross into the CaO (solid) + CO2 (vapor) field at improbably low pressures. It is found that release paths from measured shock temperatures cross into the melt plus vapor field at pressures greater than .5 GPa, which suggests that devolatilization is initiated at the shear banding sites.

Tyburczy, J. A.↗

DOE OSTI2026

Rotation‐Controlled Diurnal Evolution of Uranus' Asymmetric Bow Shock at Equinox

Uranus possesses the most extreme magnetic and rotational geometry in the solar system, resulting in a uniquely dynamic and asymmetric interaction between its magnetosphere and the solar wind. Here we investigate the diurnal evolution of the Uranian bow shock (BS) at equinox using global multifluid magnetohydrodynamic simulations constrained by Voyager 2 observations. The BS structure and variability are quantified using the stand‐off distance, terminator distance, flaring parameter, and BS parameter, enabling a systematic assessment of its global asymmetry over one planetary rotation. Our results show that the BS exhibits pronounced diurnal asymmetry that is primarily controlled by planetary rotation. Even under steady upstream solar wind conditions, the BS undergoes periodic expansion and contraction, reflecting rotation‐driven reconfiguration of the planetary magnetospheric topology. This diurnal modulation represents a remarkably large asymmetry in planetary BS geometry within the solar system. These results establish planetary rotation as the intrinsic primary driver of Uranus' BS dynamics at equinox and provide essential constraints for the design and formulation of future space missions to the ice giants. These findings also offer a unique framework for understanding shock physics in the abundant population of ice‐giant exoplanets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

Impact-generated winds on Venus: Causes and effects

The pressure of the dense atmosphere of Venus significantly changes the appearance of ejecta deposits relative to craters on the Moon and Mercury. Conversely, specific styles and sequences of ejecta emplacement can be inferred to represent different intensities of atmospheric response winds acting over different timescales. Three characteristic timescales can be inferred from the geologic record: surface scouring and impactor-controlled (angle and direction) initiation of the long fluidized run-out flows; nonballistic emplacement of inner, radar-bright ejecta facies and radar-dark outer facies; and very late reworking of surface materials. These three timescales roughly correspond to processes observed in laboratory experiments that can be scaled to conditions on Venus (with appropriate assumptions): coupling between the atmosphere and earlytime vapor/melt (target and impactor) that produces an intense shock that subsequently evolves into blast/response winds; less energetic dynamic response of the atmosphere to the outward-moving ballistic ejecta curtain that generates nonthermal turbulent eddies; and late recovery of the atmosphere to impact-generated thermal and pressure gradients expressed as low-energy but long-lived winds. These different timescales and processes can be viewed as the atmosphere equivalent of shock melting, material motion, and far-field seismic response in the target. The three processes (early Processes, Atmospheric Processes, and Late Recovery Winds) are discussed at length.

Schultz, Pater H.↗

DOE OSTI2020

High Explosive Ignition through Chemically Activated Nanoscale Shear Bands

Shock initiation and detonation of high explosives is considered to be controlled through hot spots, which are local regions of elevated temperature that accelerate chemical reactions. Using classical molecular dynamics, we predict the formation of nanoscale shear bands through plastic failure in shocked TATB high explosive crystal. By scale-bridging with quantum-based molecular dynamics, we show that shear bands exhibit lower reaction barriers. While shear bands quickly cool,they remain chemically activated and support increased reaction rates without the local heating typically evoked by the hot spot paradigm. We describe this phenomenon as chemical activation through shear banding

Chemistry - Inorganic↗

NASA NTRS2015

Increased Mach Number Capability for the NASA Glenn 10x10 Supersonic Wind Tunnel

Diffuser↗

DOE OSTI2020

Femtosecond X-ray Diffraction of Laser-shocked Forsterite (Mg 2 SiO 4 ) to 122 GPa

The response of forsterite, Mg 2 SiO 4 , under dynamic compression is of fundamental importance for understanding its phase transformations and high-pressure behavior. Here, we have carried out an in situ X-ray diffraction study of laser-shocked polycrystalline and single-crystal forsterite from 19 to 122 GPa using the Matter in Extreme Conditions end-station of the Linac Coherent Light Source. Under laser-based shock loading, forsterite does not transform to the high-pressure equilibrium assemblage of MgSiO 3 bridgmanite and MgO periclase, as has been suggested previously. Instead, we observe forsterite and forsterite III, a metastable polymorph of Mg 2 SiO 4 , coexisting in a mixed-phase region from 33 to 75 GPa for both polycrystalline and single-crystal samples. Densities inferred from X-ray diffraction are consistent with earlier gas-gun shock data. At higher stress, the response is sample-dependent. Polycrystalline samples undergo amorphization above 79 GPa. For [010]- and [001]-oriented crystals, a mixture of crystalline and amorphous material is observed to 108 GPa, whereas the [100]-oriented forsterite adopts an unknown phase at 122 GPa. The first two sharp diffraction peaks of amorphous Mg 2 SiO 4 show a similar trend with compression as those observed for MgSiO 3 in both recent static- and laser-driven shock experiments. This study provides new insight into the transformation of forsterite under nanosecond-duration shock loading. This work emphasizes the importance of formation of metastable phases along the Hugoniot and adds to evidence that the 300-K single-crystal diamond anvil cell experiments have relevance for understanding structures formed under shock compression. In particular, the metastable phase forsterite III has now been shown to form under dynamic compression from 10s to 100s of nanoseconds as well as under 300-K static compression. Upon compression to higher pressures, Mg 2 SiO 4 transforms to an amorphous phase. These results have broad relevance for understanding the behavior of silicates under dynamic compression.

36 MATERIALS SCIENCE↗

NASA NTRS2002

Solar Flares, Type III Radio Bursts, Coronal Mass Ejections, and Energetic Particles

In this correlative study between greater than 20 MeV solar proton events, coronal mass ejections (CMEs), flares, and radio bursts it is found that essentially all of the proton events are preceded by groups of type III bursts and all are preceded by CMEs. These type III bursts (that are a flare phenomenon) usually are long-lasting, intense bursts seen in the low-frequency observations made from space. They are caused by streams of electrons traveling from close to the solar surface out to 1 AU. In most events the type III emissions extend into, or originate at, the time when type II and type IV bursts are reported (some 5 to 10 minutes after the start of the associated soft X-ray flare) and have starting frequencies in the 500 to approximately 100 MHz range that often get lower as a function of time. These later type III emissions are often not reported by ground-based observers, probably because of undue attention to type II bursts. It is suggested to call them type III-1. Type III-1 bursts have previously been called shock accelerated (SA) events, but an examination of radio dynamic spectra over an extended frequency range shows that the type III-1 bursts usually start at frequencies above any type II burst that may be present. The bursts sometimes continue beyond the time when type II emission is seen and, furthermore, sometimes occur in the absence of any type II emission. Thus the causative electrons are unlikely to be shock accelerated and probably originate in the reconnection regions below fast CMEs. A search did not find any type III-1 bursts that were not associated with CMEs. The existence of low-frequency type III bursts proves that open field lines extend from within 0.5 radius of the Sun into the interplanetary medium (the bursts start above 100 MHz, and such emission originates within 0.5 solar radius of the solar surface). Thus it is not valid to assume that only closed field lines exist in the flaring regions associated with CMEs and some interplanetary particles originating in such flare regions might be expected in all solar particle events.

Cane, Hilary V.↗

NASA NTRSDate not supplied

Application of Flux Conserved Modeling to an Unsteady Combustion Driven Pseudo-Shock

Quasi-steady correlation-based and flux conserved models of the pseudo-shock pressure profile are experimentally assessed in the Air Force Research Lab RC-18 Sub-Scale DirectConnect Isolator Combustion Facility. A bifurcated normal pseudo-shock is established in the 2.67 aspect ratio isolator test section by a downstream cavity-stabilized combustor. The Waltrup-Billig correlation, Smart flux conserved model, and Fiévet data driven flux conserved models are compared against measured pressure profiles. The literature models are found to perform poorly in a high-amplitude forced environment of RC-18, where shock speeds upwards of 20% the isolator Mach number occur. Two modifications to the quasi-steady Fiévet data driven flux conserved model are proposed including incorporating an estimate of pseudo-shock static pressure recovery and exercising the model in a shock stationary reference frame. It is shown that these modifications allow the quasi-steady flux conserved modeling approach to accurately predict pseudo-shock pressure profiles over the discrete phases of high-amplitude forced dynamics.

Shock Train↗

DOE OSTI2022

Shock revival in core-collapse supernovae assisted by heavy axionlike particles

We report axionlike particles (ALPs) are a class of hypothetical pseudoscalar particles which feebly interact with ordinary matter. The hot plasma of core-collapse supernovae is a possible laboratory to explore physics beyond the standard model including ALPs. Once produced, some of the ALPs can be absorbed by the supernova matter and affect energy transfer. In this study, we calculate the ALP emission in core-collapse supernovae and the backreaction on supernova dynamics consistently. It is found that the stalled bounce shock can be revived if the coupling between ALPs and photons is as high as $g_{aγ}$ ~ 10 -9 GeV -1 and the ALP mass is 40–400 MeV. Most of the models result in more energetic explosions than the average observed supernova. While this can be used to place constraints on those ALPs, long-term simulations across multiple progenitors need to be further investigated to place robust limits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

LANL Test Engineering’s Structural Dynamics Research and Development Efforts [Slides]

LANL's priority roles are serving as a nuclear weapons design agency and a nuclear weapons production agency; addressing nuclear threats; and performing national security science, technology, and engineering. LANL’s Test Engineering organization serves the national mission of maintaining a strategic nuclear deterrent by providing high-quality, empirical evidence through the execution and assessment of weapons system and component testing. Shock and vibration environmental testing provides one key piece of evidence for evaluation and qualification of these systems and components in the service environments incurred during a lifetime in the stockpile. LANL’s Test Engineering group continuously pursues research and development in structural dynamics to improve shock and vibration testing capabilities, processes, and analysis techniques to increase the fidelity of service environment recreations in laboratory testing settings. This talk will highlight several ongoing projects to represent these R&D efforts.

42 ENGINEERING↗