Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock physics”

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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Studying the Richtmyer–Meshkov instability in convergent geometry under high energy density conditions using the Decel platform

The “Decel” platform at Sandia National Laboratories investigates the Richtmyer–Meshkov instability (RMI) in converging geometry under high energy density conditions [Knapp et al., Phys. Plasmas 27, 092707 (2020)]. In Decel, the Z machine magnetically implodes a cylindrical beryllium liner filled with liquid deuterium, launching a converging shock toward an on-axis beryllium rod machined with sinusoidal perturbations. The passage of the shock deposits vorticity along the Be/D2 interface, causing the perturbations to grow. Here, we present platform improvements along with recent experimental results. To improve the stability of the imploding liner to the magneto Rayleigh–Taylor instability, we modified its acceleration history by shortening the Z electrical current pulse. Next, we introduce a “split rod” configuration that allows two axial modes to be fielded simultaneously in different axial locations along the rod, doubling our data per experiment. We then demonstrate that asymmetric slots in the return current structure modify the magnetic drive pressure on the surface of the liner, advancing the evolution on one side of the rod by multiple ns compared to its 180° counterpart. This effectively enables two snapshots of the instability at different stages of evolution per radiograph with small deviations of the cross-sectional profile of the rod from the circular. Using this platform, we acquired RMI data at 272 and 157 μm wavelengths during the single shock stage. Finally, we demonstrate the utility of these data for benchmarking simulations by comparing calculations using ALEGRA MHD and RageRunner.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Wave-particle interactions on the FAST satellite

NASA's Fast Auroral Snapshot, or 'FAST' satellite, scheduled for launch in 1993, will investigate the plasma physics of the low altitude auroral zone from a 3500-km apogee polar orbit. FAST will give attention to wave, double-layer, and soliton production processes due to electrons and ions, as well as to wave-wave interactions, and the acceleration of electrons and ions by waves and electric fields. FAST will employ an intelligent data-handling system capacle of data acquisition at rates of up to 1 Mb/sec, in addition to a 1-Gbit solid-state memory. The data need be gathered for only a few minutes during passes through the auroral zone, since the most interesting auroral phenomena occur in such narrow regions as auroral arcs, electrostatic shocks, and superthermal electron bursts.

Temerin, M. A.↗
NASA NTRS2003

Coordination of Lower Limb Joints During Locomotion: The Effects of Vestibulo-Ocular Reflex Adaptation

Controlling locomotion while maintaining a stable gaze requires precise coordination between several, interdependent full-body sensorimotor subsystems (Bloomberg and Mulavara, 2003; McDonald, et al., 1997). The overall goal of this study is to determine how this full-body gaze stabilization system responds to adaptive changes in vestibuloocular reflex (VOR) function. Locomotion involves cyclical physical interactions (impacts) with the environment. Hence, focusing on a target and maintaining visual acuity during this activity may require mechanisms to manage the energy flow, so it does not disrupt the visual and vestibular sensory information processing that stabilizes gaze. It has been shown that increasing the difficulty of a gaze task (reading numbers on a screen as opposed to simply focusing on a central dot pattern) resulted in an increase in the amount of knee flexion movement during the critical phase immediately following the heel strike event (Mulavara and Bloomberg, 2003). The increase in knee flexion during the stance phase of the gait cycle has been suggested to function as a shock absorbing mechanism associated with the rapid weight transfer from the trailing to the leading leg during walking. To understand this full-body coordination, the relative contributions of each component and the resulting effects should be assessed. In this study, we hypothesized that VOR adaptation would result in a reorganization of the lower limb joint coordination during treadmill walking in a manner to facilitate the gaze stabilization task and preserve locomotor function.

Miller, Chris↗
NASA NTRS2009

Hypersonic Viscous Flow Over Large Roughness Elements

Viscous flow over discrete or distributed surface roughness has great implications for hypersonic flight due to aerothermodynamic considerations related to laminar-turbulent transition. Current prediction capability is greatly hampered by the limited knowledge base for such flows. To help fill that gap, numerical computations are used to investigate the intricate flow physics involved. An unstructured mesh, compressible Navier-Stokes code based on the space-time conservation element, solution element (CESE) method is used to perform time-accurate Navier-Stokes calculations for two roughness shapes investigated in wind tunnel experiments at NASA Langley Research Center. It was found through 2D parametric study that at subcritical Reynolds numbers, spontaneous absolute instability accompanying by sustained vortex shedding downstream of the roughness is likely to take place at subsonic free-stream conditions. On the other hand, convective instability may be the dominant mechanism for supersonic boundary layers. Three-dimensional calculations for both a rectangular and a cylindrical roughness element at post-shock Mach numbers of 4.1 and 6.5 also confirm that no self-sustained vortex generation from the top face of the roughness is observed, despite the presence of flow unsteadiness for the smaller post-shock Mach number case.

Chang, Chau-Lyan↗
DOE OSTI2022

Fast ignition inertial fusion energy using laser-driven ion beams

Ion fast ignition (IFI), or fusion fast ignition initiated by a laser-driven ion beam, is a promising path to high-gain inertial fusion energy (IFE). In IFI, cold, dense deuterium-tritium (DT) fuel is first assembled using lasers or pulsed power drivers. Then, a high-power ion beam is focused onto a small volume within the fuel (the hot spot), heating the fuel rapidly to conditions where fusion ignition takes place. Fusion burn in this hot spot propagates to the fuel surrounding the hot spot, leading to burnup of a significant fraction of this fuel and the possibility of high gain (G~100), as needed for inertial fusion energy. IFI uses separate drivers for the two basic elements, fuel compression and ignition, allowing maximum control and optimization of each. On the other hand, conventional laser fusion uses multiple beams of the same driver to compress the fuel and shock-heat its very center to ignite a burn wave. Despite impressive progress in conventional laser fusion, the precise spatial symmetry, temporal pulse shaping and timing required for high gain and IFE remain a serious unmet challenge.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Transportation Fatigue Testing of the pHB650 Power Coupler Antenna for PIP-II at FNAL

The PIP-II Project will see international shipment of cryomodules Proton Improvement Plan-II from Europe to the United States, and as such, the shocks which can occur during shipment pose a risk to the internal components. Of particular concern is the coupler ceramic window and surrounding brazes, which will see stresses during an excitation event. Since the antenna design is new, and because of the setback failure would create, a cyclic stress test was devised for the antenna. This paper presents the experimental methods, setup, and results of the test.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1976

Interplanetary disturbances in the solar wind produced by density, temperature, or velocity pulses at 0.08 AU

Time-dependent solutions of a one-fluid model of the interplanetary medium are investigated. This set of unsteady hydrodynamic equations has been written in conserved form in order to apply the Lax-Wendroff (1960) method for the solution of this problem. The initial disturbance is specified by a pulse at 0.08 AU. Physically, this pulse can be interpreted as having been caused by a solar flare, surge, or any other solar disturbance. The equilibrium condition is determined to be the steady solution of the governing equations and represents the quiet solar wind. Results are presented in terms of density, temperature, and velocity profiles of the interplanetary gas flow at heliocentric distances up to 6 AU at several times. Also, the trajectories of disturbances for various initial pulses are shown. Some June 1972 interplanetary observational data are compared with these theoretical calculations. On the basis of these results, the effects of solar disturbances on the interplanetary environment (such as the generation of large nonlinear wave trains in the shock wakes) can be inferred.

Wu, S. T.↗
DOE OSTI2025

Thermal and kinematic properties of ejecta in SN1987A revealed by XRISM

We present an analysis of high-resolution spectra from the shock-heated plasmas in SN 1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM). The 1.7–10 keV Resolve spectra are accurately represented by a single-component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58}$ × $10^{11}\,\,{\rm s\,\, cm}^{-3}$. The Resolve spectra are also well reproduced by the 3D magneto-hydrodynamic simulation presented by Orlando et al. (2020, A&A, 636, A22) suggesting substantial contribution from the ejecta. The metal abundances obtained with Resolve align with the Large Magellanic Cloud value, indicating that the X-rays in 2024 originate from “non-metal-rich” shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of ejecta. Doppler widths of the atomic lines from Si, S, and Fe correspond to velocities of 1500–1700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a $90\%$ upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV range, aligning with NuSTAR findings by Greco et al. (2022, ApJ, 931, 132). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4}\, M_{\odot }$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015, Science, 348, 670; Grebenev et al. 2012, Nature, 490, 373; Leising 2006, ApJ, 651, 1019).

ISM: supernova remnants↗
DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗
DOE OSTI2023

VIPIR: A High-Throughput Drop-Weight Impact Instrument for Imaging the Initiation and Propagation of Reactions in Energetic Materials

Characterizing the handling safety and sensitivity of explosives has been a challenging area of study for over 60 years. Historically one of the most accessible and widely utilized experiments has been the drop-weight impact test, which involves dropping a weight on a small sample sandwiched between two anvils. Because this experiment generally only utilizes sound thresholds to determine whether or not a sample reacted, the physical and chemical properties governing sensitivity remain convolved. Better understanding of chemical and material characteristics is needed to give the chemistry and engineering communities a predictive tool to determine the handling sensitivity of explosives prior to pursuing expensive and potentially hazardous synthesis and formulation operations. Here, we are developing a high throughput drop tower instrument capable of imaging the deformation and flow of energetic materials during impact and the resulting thermal ignition and propagation events. This instrument is based on previous designs but has been modified for higher throughput and tailorable modifications in the future. Herein, we present key design features that enable high-speed visible and thermal imaging of explosive initiation by sub-shock impacts, as well as preliminary results in which ignition sites were observed in an erythritol tetranitrate sample.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Sierra/SD – Verification Test Manual – 5.22

Verification and validation (V&V) of scientific computing programs are important at Sandia National Labs due to the expanding role of computational simulation in managing the United States nuclear stockpile. The complexities of structural response calculations used to analyze physical problems, the varieties of codes applied to the calculations, and the importance of accurate predictions when assessing field conditions demand confidence in the consistency and accuracy of computer codes. Confidence in the accuracy of the predictions arising from computer simulations must ultimately be gained through verification and validation. The Sierra salinas structural dynamics analysis code, Sierra/SD, is used at the DOE Laboratories, and in several DOD projects. The roles of Sierra/SD in the qualification of weapon systems and components for normal and hostile environments throughout the Stockpile-to-Target Sequence include to, • Redesign weapon components. • Certify weapon components and systems for target environments such as hypersonic vehicles. • Certify that components will survive the thermal mechanical shock loads associated with hostile environments. • Evaluate current stockpile issues, including issues associated with uncertainty quantification. • Address many other problems that are encountered in stockpile management. The Sierra/SD verification plan is described, and an evolving set of key verification tests are described in detail. The verification tests ensure the correctness of the mathematics and numerical algorithms associated with functionality describing engineering phenomena. Development is in accordance with a set of tailored Software Quality Engineering (SQE) practices. SQE practices guide the overall verification and validation effort.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Ion-scale Current Structures in Short Large-amplitude Magnetic Structures

Magnetic reconnection is an explosive phenomenon widely happen in plasma environments. Recent studies show that reconnection could occur in the shock transition region. We study thin current layers in the region in front of the shock and find evidence of reconnection. The analysis also shows the possible generation mechanism of the thin current sheet that further enables reconnection. A large-scale wave is further generated, and the thin current sheets are generated through compression or secondary waves. The thin current structures also lead to electron heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Velocity-space Origins of the Pressure–Strain Interaction in Multipopulation Distributions and Its Application to Magnetic Reconnection

A forefront research question is how energy evolves in weakly collisional plasmas for which departures from local thermodynamic equilibrium (LTE) are significant. The standard approach is studying the terms in the non-LTE energy evolution equation derived by taking the second moment of the Boltzmann equation, but the resultant fluid metrics do not retain information about which particles at which velocities drive energy evolution. A widely studied channel for internal energy density evolution is the pressure–strain interaction. Here, we employ the kinetic pressure–strain, a phase-space diagnostic whose velocity-space integral recovers the pressure–strain interaction to disambiguate the contributions to the pressure–strain interaction from disparate particle populations in composite phase-space densities. We develop phase-space analogs of the pressure–strain interaction decompositions to provide the phase-space origins of normal versus sheared flow. We introduce the “kinetic strain-rate” tensor, the phase-space analog of the strain-rate tensor, which we argue is needed to interpret the phase-space origins of the pressure–strain interaction. To demonstrate the utility of these quantities, we investigate them for composite electron distributions near the electron diffusion region in two-dimensional particle-in-cell simulations of antiparallel symmetric magnetic reconnection. We find that the phase-space-based diagnostics isolate the roles of distinct populations. These results contribute to a growing body of work providing new methods for quantifying phase-space energy evolution for a broad array of processes, from magnetic reconnection to collisionless shocks and turbulence, opening new pathways for answering longstanding problems of particle energization in weakly collisional plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1999

What We Don't Understand About Ion Acceleration Flares

There are now strong associations between the (3)He-rich, Fe-rich ions in "impulsive" solar energetic particle (SEP) events and the similar abundances derived from gamma-ray lines from flares. Compact flares, where wave energy can predominate, are ideal sites for the study of wave-particle physics. Yet there are nagging questions about the magnetic geometry, the relation between ions that escape and those that interact, and the relative roles of cascading Alfven waves and the EMIC waves required to enhance He-3. There are also questions about the relative timing of ion and electron acceleration and of heating; these relate to the variation of ionization states before and during acceleration and during transport out of the corona. We can construct a model that addresses many of these issues, but problems do remain. Our greatest lack is realistic theoretical simulations of element abundances, spectra, and their variations. By contrast, we now have a much better idea of the acceleration at CME-driven shock waves in the rare but large "gradual" SEP events, largely because of their slow temporal evolution and great spatial extent.

Reames, Donald V.↗
NASA NTRS2002

A Comparative Study of Three Methodologies for Modeling Dynamic Stall

During the past two decades, there has been an increased reliance on the use of computational fluid dynamics methods for modeling rotors in high speed forward flight. Computational methods are being developed for modeling the shock induced loads on the advancing side, first-principles based modeling of the trailing wake evolution, and for retreating blade stall. The retreating blade dynamic stall problem has received particular attention, because the large variations in lift and pitching moments encountered in dynamic stall can lead to blade vibrations and pitch link fatigue. Restricting to aerodynamics, the numerical prediction of dynamic stall is still a complex and challenging CFD problem, that, even in two dimensions at low speed, gathers the major difficulties of aerodynamics, such as the grid resolution requirements for the viscous phenomena at leading-edge bubbles or in mixing-layers, the bias of the numerical viscosity, and the major difficulties of the physical modeling, such as the turbulence models, the transition models, whose both determinant influences, already present in static maximal-lift or stall computations, are emphasized by the dynamic aspect of the phenomena.

Sankar, L.↗
NASA NTRS2010

The First Detection of (O IV) from an Ultraluminous X-Ray Source with Spitzer. 2. Evidence for High Luminosity in Holmberg II ULX

This is the second of two papers examining Spitzer Infrared Spectrograph observations of the ultraluminous X-ray source (ULX) in Holmberg II. Here, we perform detailed photoionization modeling of the infrared lines. Our analysis suggests that the luminosity and morphology of the [Oiv] 25.89 micronmeters emission line is consistent with photoionization by the soft X-ray and far ultraviolet (FUV) radiation from the accretion disk of the binary system and inconsistent with narrow beaming. We show that the emission nebula is matter bounded both in the line-of-sight direction and to the east, and probably radiation-bounded to the west. A bolometric luminosity in excess of 10(exp 40) erg s(exp -1) would be needed to produce the measured [O iv] flux. We use modeling and previously published studies to conclude that shocks likely contribute very little, if at all,to the high-ionization line fluxes observed in the Holmberg II ULX. Additionally, we find that the spectral type of the companion star has a surprisingly strong effect on the predicted strength of the [O iv] emission. This finding could explain the origin of [O iv] in some starburst systems containing black hole binaries.

GALAXIES↗
DOE OSTI2026

A point-particle-based hydride shell-shedding model for ejecta particle transport in reactive environments

A shock wave passing over a rough or perturbed metal surface will induce a limiting case of Richtmyer–Meshkov instability and will cause small particles to eject from the surface and transport into the surrounding medium. These particles are known as ejecta and can be either solid or liquid in nature. Recent experiments have shown that liquid cerium ejecta clouds exhibit unexpected non-monotonic acceleration behaviors as well as temperature plateaus after a brief temperature rise if they are transporting in a chemically reactive, hydrogen-based medium while they act as expected in an inert medium. This work details a point-particle model developed for reactive cerium ejecta transport, which attempts to account for these new physics through the behavior of a developing solid hydride shell, which is believed to form as a product of the reaction. The overall model incorporates the effects of the reaction on the particle properties as well as the effects of potential shedding of the shell into sub-micrometer scale flakes and potential phase change of the hydride if the ejecta particles reach the melt point of the hydride layer. The model is tested by performing simulations of the original motivating experiments and comparing quantities, such as ejected mass, velocimetry, and temperature profiles, against the experimental data. While the model is able to capture many general features of the observed anomalies, some inaccuracies still exist. These point to both missing physics in the model (such as a deuterium adsorption mechanism on the hydride layer) as well as a lack of knowledge of certain material properties (such as the strength of cerium hydride to determine dynamic fracture thicknesses) needed to fully reduce the uncertainties in the model by up to an order of magnitude and perform a true attempt at model validation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 2: Homogeneous medium at Mach 3 with gamma = 5/3

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing spheres of varying sizes (from 10 down to 0.01 accretion radii) move at Mach 3 relative to a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3). To accommodate the long-range gravitational forces, the extent of the computational volume is 32(exp 3) accretion radii. We examine the influence of numerical procedure on physical behavior. The hydrodynamics is modeled by the 'piecewise parabolic method.' No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (5-10) grids around the sphere, each finer grid being a factor of 2 smaller in zone dimension that the next coarser grid. The largest dynamic range (ratio of size of the largest grid to size of the finest zone) is 16,384. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid, while at the same time evolving the gas on the coarser grids. Initially (at time t = 0-10), a shock front is set up, a Mach cone develops, and the accretion column is observable. Eventually the flow becomes unstable, destroying axisymmetry. This happens approximately when the mass accretion rate reaches the values (+/- 10%) predicted by the Bondi-Hoyle accretion formula (factor of 2 included). However, our three-dimensional models do not show the highly dynamic flip-flop flow so prominent in two-dimensional calculations performed by other authors. The flow, and thus the accretion rate of all quantities, shows quasi-periodic (P approximately equals 5) cycles between quiescent and active states. The interpolation formula proposed in an accompanying paper is found to follow the collected numerical data to within approximately 30%. The specific angular momentum accreted is of the same order of magnitude as the values previously found for two-dimensional flows.

Ruffert, Maximilian↗