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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Stability of radiative shocks with time-dependent cooling

A full evolutionary calculation of ion abundances and radiative cooling has been incorporated into an accurate one-dimensional gasdynamics calculation in order to investigate the dynamics of radiative shock instability and to determine the extent to which power-law cooling models describe the situation. Radiative shocks are shown to be unstable for velocities of greater than about 140 km/s. The temperature dependence of cooling functions behind steady shocks suggests that the stability limit may depend upon the amplitude of the perturbation. Oscillation is noted in the fundamental oscillation mode.

Gaetz, T. J.↗
DOE OSTI2022

Adaptive Space-Time Methods for Large Scale Optimal Design

When modeling complex physical systems with advanced dynamics, such as shocks and singularities, many classic methods for solving partial differential equations can return inaccurate or unusable results. One way to resolve these complex dynamics is through r-adaptive refinement methods, in which a fixed number of mesh points are shifted to areas of high interest. The mesh refinement map can be found through the solution of the Monge-Ampére equation, a highly nonlinear partial differential equation. Due to its nonlinearity, the numerical solution of the Monge-Ampére equation is nontrivial and has previously required computationally expensive methods. In this report, we detail our novel optimization-based, multigrid-enabled solver for a low-order finite element approximation of the Monge-Ampére equation. This fast and scalable solver makes r-adaptive meshing more readily available for problems related to large-scale optimal design. Beyond mesh adaptivity, our report discusses additional applications where our fast solver for the Monge-Ampére equation could be easily applied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A description of structured waves in shock compressed particulate composites

Dynamic compression of composite materials is of scientific interest because the mechanical mismatch between internal phases challenges continuum theories. Typical assumptions about steady wave propagation and quasi-instantaneous state changes require reexamination along with the need for time-dependent models. To that end, data and models are presented here for the shock compression of an idealized particulate composite. To serve as a generic representative of this material class, a polymer matrix was filled with tungsten particles, ranging from 1 to 50 vol. %. This creates a simple microstructure containing randomly scattered particles with an extreme impedance mismatch to the binding matrix. These materials were parallel plate impact loaded by Al flyers traveling at 1.8–5.0 km/s. Velocimetry provided records of the equilibrium state and the compression wave structure for each case with trends quantified by an empirical fit. The same quantities were also studied as a function of the wave's propagation distance. A homogenized viscoelastic model then made it possible to progress from cataloging trends to making predictions. Starting from a Mie–Gruneisen equation of state, additional time varying terms were added to capture the transient response. After calibration, accurate predictions of the steady wave structure were possible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2000

A Shocking New Pump

Hydro Dynamics, Inc. received a technical helping hand from NASA that made their Hydrosonic Pump (HPump) a reality. Marshall engineers resolved a bearing problem in the rotor of the pump and recommended new bearings, housings and mounting hardware as a solution. The resulting HPump is able to heat liquids with greater energy efficiency using shock waves to generate heat.

Source record↗
DOE OSTI2024

High-throughput quantification of quasistatic, dynamic and spall strength of materials across 10 orders of strain rates

Abstract The response of metals and their microstructures under extreme dynamic conditions can be markedly different from that under quasistatic conditions. Traditionally, high strain rates and shock stresses are achieved using cumbersome and expensive methods such as the Kolsky bar or large spall experiments. These methods are low throughput and do not facilitate high-fidelity microstructure–property linkages. In this work, we combine two powerful small-scale testing methods, custom nanoindentation, and laser-driven microflyer (LDMF) shock, to measure the dynamic and spall strength of metals. The nanoindentation system is configured to test samples from quasistatic to dynamic strain-rate regimes. The LDMF shock system can test samples through impact loading, triggering spall failure. The model material used for testing is magnesium alloys, which are lightweight, possess high-specific strengths, and have historically been challenging to design and strengthen due to their mechanical anisotropy. We adopt two distinct microstructures, solutionized (no precipitates) and peak-aged (with precipitates) to demonstrate interesting upticks in strain-rate sensitivity and evolution of dynamic strength. At high shock-loading rates, we unravel an interesting paradigm where the spall strength vs. strain rate of these materials converges, but the failure mechanisms are markedly different. Peak aging, considered to be a standard method to strengthen metallic alloys, causes catastrophic failure, faring much worse than solutionized alloys. Our high-throughput testing framework not only quantifies strength but also teases out unexplored failure mechanisms at extreme strain rates, providing valuable insights for the rapid design and improvement of materials for extreme environments.

Eswarappa Prameela, Suhas (ORCID:0000000334530184)↗
NASA NTRS1985

Evolution and interaction of interplanetary shocks

Whang's unsteady, one-dimensional, one-fluid MHD model is used to carry out a simulation study of the evolution of the solar wind based on two shock events over a distance of the order of 10 AU in the outer heliosphere. The study is based on the observation of two events, each observed by two or more spacecraft. The results show that the shock process, including the formation, collision, and merging of shocks, dominates the dynamical evolution of large-scale solar wind structures. In the outer heliosphere, the large-scale solar wind and magnetic field evolve into a much simpler structure, and MHD shocks are present as a principal component of the solar wind. The simulation results shed new light on the interaction and evolution of large interplanetary streams.

Whang, Y. C.↗
NASA NTRS1991

The dynamic instability of adiabatic blast waves

Adiabatic blastwaves, which have a total energy injected from the center E varies as t(sup q) and propagate through a preshock medium with a density rho(sub E) varies as r(sup -omega) are described by a family of similarity solutions. Previous work has shown that adiabatic blastwaves with increasing or constant postshock entropy behind the shock front are susceptible to an oscillatory instability, caused by the difference between the nature of the forces on the two sides of the dense shell behind the shock front. This instability sets in if the dense postshock layer is sufficiently thin. The stability of adiabatic blastwaves with a decreasing postshock entropy is considered. Such blastwaves, if they are decelerating, always have a region behind the shock front which is subject to convection. Some accelerating blastwaves also have such region, depending on the values of q, omega, and gamma where gamma is the adiabatic index. However, since the shock interface stabilizes dynamically induced perturbations, blastwaves become convectively unstable only if the convective zone is localized around the origin or a contact discontinuity far from the shock front. On the other hand, the contact discontinuity of accelerating blastwaves is subject to a strong Rayleigh-Taylor instability. The frequency spectra of the nonradial, normal modes of adiabatic blastwaves have been calculated. The results have been applied to the shocks propagating through supernovae envelopes. It is shown that the metal/He and He/H interfaces are strongly unstable against the Rayleigh-Taylor instability. This instability will induce mixing in supernovae envelopes. In addition the implications of this work for the evolution of planetary nebulae is discussed.

Ryu, Dongsu↗
NASA NTRS1995

Kinetic structure of intermediate shocks: Implications for the magnetopause

A general study of the structure and stability of intermediate shocks (IS) in an isotropic plasma is presented using a hybrid as well as a resistive Hall MHD code. Special emphasis is put on the question of whether the rotational layers observed at the magnetopause can be intermediate shocks. The shocks are formed dynamically by the interaction between a flowing plasma and a stationary piston. Coplanar ISs (both strong and weak) are found to be stable in a collisionless plasma. The existence of slow shocks in a high beta plasma is also established for the first time. Noncoplanar ISs are found to be time-dependent, evolving toward a rotational discontinuity (RD) after some characteristic time tau which can be quite long (1000 Omega(exp -1), where Omega is the ion gyrofrequency). The value tau is larger the closer the rotation angle is to 180 deg. Rotations larger than 180 deg are found to be unstable, decaying into a state of minimum shear (i.e., rotation angle less than 180 deg). There are various length scales associated with an IS in the kinetic regime. The shortest scale is found to be the length scale over which rotation of the transverse component of the magnetic field takes place. This scale can have a half width as small as one ion inertial length (c/omega(sub p)) for electron sense rotations and 3c/omega(sub p) for ion sense rotations, for an upstream ion beta of unity. Both of these scales are consistent with the observed thickness at the magnetopause and identical to the corresponding RD scales. A detailed study of the mode conversion of the Alfven ion cyclotron waves (A/IC) waves across both slow and intermediate shocks and the resulting downstream wave spectrum are presented. The possibility that the large number of relfected ions observed at the magnetopause may be due to the presence of strong ISs is considered. The identification of strong ISs and their distinction from RDs should be possible in observations due to significant differences that exist between jump conditions and overall structure of the two discontinuities. The jumps in the plasma parameters across a weak IS are typically small. This together with the fact that the weak ISs and RDs have very similar thickness and other overall properties makes the distinction between weak ISs and RDs in the observations largely inconsequential. However, at large noncoplanarity angles the weak IS approaches the RD limit in a relatively short time (approximately less than 100 Omega (exp -1)). Thus, magnetopause rotations with large noncoplanarity angles are most likely either RDs or strong ISs. Finally, direct comparisons between fluid (resistive Hall MHD) and kinetic simulations show that fluid theory is not applicable to study of ISs in a collisionless plasma.

Karimabadi, H.↗
NASA NTRS1994

Structure, Shape and Movement of the Termination Shock

In this paper we review current ideas concerning the morphology and dynamics of the heliospheric termination shock. The mean distance of the shock is expected to lie somewhere in the range 50-200 AU, and should vary with latitude, depending on mean latitude variation of the dynamical pressure of the solar wind. The shock is expected to move in response to variation in upstream solar wind conditions, so that at any particular instant the termination shock will resemble a distorted asymmetric balloon with some parts moving inward and others moving outward. The simplest model of the shock is that of an infinitely thin gasdynamic or magnetohydrodynamic shock wave, which moves inward or outward at speeds of approximately 100 kilometers per second. The structure and motion of the shock will be strongly modified if an appreciable fraction of the incident solar wind energy goes into acceleration of energetic particles, specifically the anomalous cosmic rays.

Barnes, Aaron↗
NASA NTRS1999

Gas Evolution Dynamics in Godunov-Type Schemes and Analysis of Numerical Shock Instability

In this paper we are going to study the gas evolution dynamics of the exact and approximate Riemann solvers, e.g., the Flux Vector Splitting (FVS) and the Flux Difference Splitting (FDS) schemes. Since the FVS scheme and the Kinetic Flux Vector Splitting (KFVS) scheme have the same physical mechanism and similar flux function, based on the analysis of the discretized KFVS scheme the weakness and advantage of the FVS scheme are closely observed. The subtle dissipative mechanism of the Godunov method in the 2D case is also analyzed, and the physical reason for shock instability, i.e., carbuncle phenomena and odd-even decoupling, is presented.

Xu, Kun↗
DOE OSTI2023

Observation of subcritical shocks in a collisional laboratory plasma: scale dependence near the resistive length

We present a study of subcritical shocks in a highly collisional laboratory plasma with a dynamically significant magnetic field. Shocks were produced by placing cylindrical obstacles into the supermagnetosonic ( $M_{{\rm MS}} \sim 1.9$ ) outflow from an inverse wire array z-pinch at the MAGPIE pulsed power facility ( $n_{e} \sim 8.5 \times 10^{17}\,{\rm cm}^{-3}$ , $v \sim 45\,{\rm km}\,{\rm s}^{-1}$ ). We demonstrate the existence of subcritical shocks in this regime and find that secondary stagnation shocks form in the downstream which we infer from interferometry and optical Thomson scattering measurements are hydrodynamic in nature. The subcritical shock width is found to be approximately equal to the resistive diffusion length and we demonstrate the absence of a jump in hydrodynamic parameters. Temperature measurements by collective optical Thomson scattering showed little temperature change across the subcritical shock ( ${<}10\,\%$ of the ion kinetic energy) which is consistent with a balance between adiabatic and Ohmic heating and radiative cooling. We demonstrate the absence of subcritical shocks when the obstacle diameter is less than the resistive diffusion length due to decoupling of the magnetic field from the plasma. These findings are supported by magnetohydrodynamic simulations using the Gorgon and AstroBEAR codes and discrepancies between the simulations and experiment are discussed.

Physics↗
NASA NTRS2009

Shock Layer Radiation Measurements and Analysis for Mars Entry

NASA's In-Space Propulsion program is supporting the development of shock radiation transport models for aerocapture missions to Mars. A comprehensive test series in the NASA Antes Electric Arc Shock Tube facility at a representative flight condition was recently completed. The facility optical instrumentation enabled spectral measurements of shocked gas radiation from the vacuum ultraviolet to the near infrared. The instrumentation captured the nonequilibrium post-shock excitation and relaxation dynamics of dispersed spectral features. A description of the shock tube facility, optical instrumentation, and examples of the test data are presented. Comparisons of measured spectra with model predictions are also made.

Bose, Deepak↗
DOE OSTI2024

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗
NASA NTRS1971

Dynamic compression of enstatite

Shock wave data for Bamle enstatite in 60-480 kb range, considering Hugoniot elastic limit and phase transition produced shock front

Ahrens, T. J.↗
NASA NTRS1988

The 58th Shock and Vibration Symposium

Dynamic analysis, dynamic testing, space shuttle main engine vibration, isolation and damping and analytical methods are among the topics discussed.

Pilkey, Walter D.↗
NASA NTRS1991

Numerical studies of electron dynamics in oblique quasi-perpendicular collisionless shock waves

Linear and nonlinear electron damping of the whistler precursor wave train to low Mach number quasi-perpendicular oblique shocks is studied using a one-dimensional electromagnetic plasma simulation code with particle electrons and ions. In some parameter regimes, electrons are observed to trap along the magnetic field lines in the potential of the whistler precursor wave train. This trapping can lead to significant electron heating in front of the shock for low beta(e). Use of a 64-processor hypercube concurrent computer has enabled long runs using realistic mass ratios in the full particle in-cell code and thus simulate shock parameter regimes and phenomena not previously studied numerically.

Liewer, P. C.↗