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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Oxidation of methylamine

A detailed chemical kinetic model for oxidation of methylamine has been developed, based on theoretical work and a critical evaluation of data from the literature. In this work, the rate coefficients for the reactions of CH 3 NH 2 + O 2 → CH 2 NH 2 / CH 3 NH + HO 2 , CH 3 NH 2 + H → CH 3 + NH 3 , CH 3 NH → CH 2 NH 2 , and CH 3 NH + O 2 → CH 2 NH + HO 2 were calculated from ab initio theory. The mechanism was validated against experimental results from batch reactors, flow reactors, shock tubes, and premixed flames. The model predicts satisfactorily explosion limits for CH 3 NH 2 and its oxidation in a flow reactor. However, oxidation in the presence of nitric oxide, which strongly promotes reaction at lower temperatures, is only described qualitatively. Furthermore, calculated flame speeds are higher than reported experimental values; the model does not capture the inhibiting effect of the NH 2 group in CH 3 NH 2 compared to CH 4 . More work is desirable to confirm the products of the CH 3 NH + NO reaction and to look into possible pathways to NH 3 in methylamine oxidation.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Modified Winterbottom Construction Including Boundaries

There has been extensive work on the equilibrium shape of isolated nanoparticles with internal boundaries, and also single crystals on substrates. Surprisingly, almost shockingly, there has been very little work on the equilibrium shape of particles with internal boundaries on substrates. Here in this paper, the general solution is given for the configuration of particles which contain twin and other grain boundaries on a flat substrate, which can be applied to any polycrystalline or multiphase nanoparticle configuration. The solution is based upon combining the established modified-Wulff construction that has been extensively validated for twinned particles with the Winterbottom construction for single particles on a substrate. The solution is illustrated for the specific case of five-fold multiply twinned particles. Good agreement is observed between both existing experimental data in the literature as well as some experimental data included within this work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Development and modeling for a small-scale, rapidly heated high explosives initiation time (HEIT) experiment

Most small-scale assessments of explosive sensitivity, including the popular drop-weight impact test, convolute thermal and mechanical phenomena to the extent that it has been extremely challenging to decipher how an explosive ignites and propagates reactions. For instance, an impact generates heat through a range of dissipation mechanisms, which can in turn, depending on the reaction rates of the explosive, lead to chemical decomposition. To deconvolute the various contributions to the sub-shock initiation and propagation of explosive reactions, we describe the development and modeling of the High Explosives Initiation Time (HEIT) test - a new, small-scale, high-throughput experiment designed to rapidly heat milligram quantities of energetic materials confined within small diameter steel needles. Specifically, we have modeled and designed a 250 joule pulsed power system capable of rapidly delivering electrical current to the needles, resulting in rapid heat delivery to the sample. In conclusion, the energy deposition rate into the sample is controlled by different transmission line topologies. Modeling in COMSOL is performed to understand the energy required to heat up the explosive sample.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Pump depletion and hot-electron generation in long-density-scale-length plasma with shock-ignition high-intensity laser

Two-dimensional particle-in-cell simulations for laser plasma interaction with laser intensity of 10 16 W/cm 2 , plasma density range of 0.01–0.28$n_c$, and scale length of 230–330 $\textit{μ}$m showed significant pump depletion of the laser energy due to stimulated Raman scattering (SRS) and stimulated Brillouin scattering (SBS) in the low-density region ($n_e$ = 0.01–0.2$n_c$). The simulations identified hot electrons generated by SRS in the low-density region with moderate energy and by two-plasmon-decay near $n_e$ = 0.25$n_c$ with higher energy. The overall hot electron temperature (46 keV) and conversion efficiency (3%) were consistent with the experiment's measurements. Further, the simulations also showed artificially reducing SBS would lead to stronger SRS and a softer hot-electron spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Cross-correlation of the extragalactic gamma-ray background with the thermal Sunyaev-Zel’dovich effect in the cosmic microwave background

Cosmic rays in galaxy clusters are unique probes of energetic processes operating with large-scale structures in the Universe. Precise measurements of cosmic rays in galaxy clusters are essential for improving our understanding of nonthermal components in the intracluster medium (ICM) as well as the accuracy of cluster mass estimates in cosmological analyses. In this paper, we perform a cross-correlation analysis with the extragalactic gamma-ray background and the thermal Sunyaev-Zeldovich (tSZ) effect in the cosmic microwave background. The expected cross-correlation signal would contain rich information about the cosmic-ray-induced gamma-ray emission in the most massive galaxy clusters at $\textit{z}$ ~ 0.1 – 0.2 . We analyze the gamma-ray background map with 8 years of data taken by the Large Area Telescope onboard Fermi satellite and the publicly available tSZ map by Planck. We confirm that the measured cross-correlation is consistent with a null detection, and thus it enables us to put the tightest constraint on the acceleration efficiency of cosmic-ray protons at shocks in and around galaxy clusters. In this work, we find the acceleration efficiency must be below 5% with a 2σ confidence level when the hydrostatic mass bias of clusters is assumed to be 30%, and our result is not significantly affected by the assumed value of the hydrostatic mass bias. Our constraint implies that the nonthermal cosmic-ray pressure in the ICM can introduce only a≲3% level of the hydrostatic mass bias, highlighting that cosmic rays alone do not account for the mass bias inferred by the Planck analyses. Finally, we discuss future detectability prospects of cosmic-ray-induced gamma rays from the Perseus cluster for the Cherenkov Telescope Array.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Reduced mixing in inertial confinement fusion with early-time interface acceleration

In inertial confinement fusion (ICF) implosions, the interface between the cryogenic DT fuel and the ablator is unstable to shock acceleration (the Richtmyer-Meshkov instability, RM) and constant acceleration (Rayleigh-Taylor instability, RT). Instability growth at this interface can reduce the final compression, limiting fusion burnup. If the constant acceleration is in the direction of the lighter material (negative Atwood number), the RT instability produces oscillatory motion that can stabilize against RM growth. Theory and simulations suggest this scenario occurred at early times in some ICF experiments on the National Ignition Facility, possibly explaining their favorable performance compared to one-dimensional simulations. In conclusion, this characteristic is being included in newer, lower adiabat designs, seeking to improve compression while minimizing ablator mixing into the fuel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Three-dimensional signatures of self-similarity in a high-energy-density plasma shear-driven mixing layer

A hydrodynamic shear mixing layer experiment at the National Ignition Facility had previously demonstrated Eulerian scaling of integrated, late-time quantities, including turbulent kinetic energy. In this manuscript, the experiment is repeated with new materials. In this work, using the new dataset, we demonstrate that Euler-number scalings hold not just for late time, but dynamically throughout the experiment, for measurements in all three spatial dimensions. Incorporating the dynamic scaling leads to an enhanced calculation that the heavier of the two scaled experiments has approached three generations of mergers of its primary instability's structures and a consistent observation of such a merger in action in the lighter of the two scaled experiments. Furthermore, the improved scrutiny of the time evolution of instability structures leads to sharper estimates of turbulent kinetic energy, including a demonstration of different behaviors correlating with surface roughness (quantitatively consistent with transitions between laminar and turbulent initial states), as predicted by a Reynolds-averaged turbulent model, which evidently correctly handles the differing shock-roughness interactions to drive its internal state of the model into different regimes. Altogether, a picture arises of the analytical improvements in treating these variations (of times, densities, and roughnesses) as a unified whole and of multiple ways by which deviations from the scaling could indicate an onset of non-hydrodynamic behavior. Such deviations were not expected for these experiments (which models correctly indicated would remain hydrodynamic) but could be introduced by, for example, imposing external fields or increasing drive energy to test conditions relevant to inertial confinement fusion or other high-energy-density experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The Coupling of an EUV Coronal Wave and Ion Acceleration in a Fermi-LAT Behind-the-Limb Solar Flare

We present the Fermi-LAT observations of the behind-the-limb (BTL) flare of 2021 July 17 and the joint detection of this flare by STIX on board the Solar Orbiter. The separation between Earth and the Solar Orbiter was $99^{\circ}_{\centerdot}2$ at 05:00 UT, allowing STIX to have a front view of the flare. The location of the flare was S20E140 in Stonyhurst heliographic coordinates, making this the most distant behind-the-limb flare ever detected in >100 MeV gamma-rays. The LAT detection lasted for ~16 minutes, the peak flux was 3.6 ± 0.8 (10 –5 ) ph cm –2 s –1 with a significance >15σ. A coronal wave was observed from both STEREO-A and SDO in extreme ultraviolet (EUV), with an onset on the visible disk in coincidence with the LAT onset. A complex type II radio burst was observed by GLOSS also in coincidence with the onset of the LAT emission, indicating the presence of a shock wave. We discuss the relation between the time derivative of the EUV wave intensity profile at 193 Å as observed by STEREO-A and the LAT flux to show that the appearance of the coronal wave at the visible disk and the acceleration of protons as traced by the observed >100 MeV gamma-ray emission are coupled. We also report how this coupling is present in the data from three other BTL flares detected by Fermi-LAT, suggesting that the protons driving the gamma-ray emission of BTL solar flares and the coronal wave share a common origin.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Effects of Turbulence on Heliosheath Ions and Implications for Energetic Neutral Atoms

The distribution of ions in the heliosheath—the region between the heliospheric termination shock and the heliopause—is important for understanding remote observations of energetic neutral atoms (ENAs). The ion distributions were estimated previously based on hybrid simulations of the heating and evolution of solar wind and interstellar pickup ions across the solar wind termination shock, but these estimates only provide the distributions near the shock. In this work, we use self-consistent hybrid kinetic simulations to investigate the effects of turbulence on ion distributions in the heliosheath. The simulations are compared against Voyager observations, constraining the feasible amplitude and compressibility of turbulence. We find that the heating due to turbulent dissipation can lead to a significant increase in the temperature of thermal solar wind ions. Both turbulent velocity fluctuations and the heating of solar wind ions increase the charge-exchange source for ENAs at low energies (around 100 eV), where current ENA models underpredict observations by more than an order of magnitude. However, the effects of turbulence are likely not strong enough to fully explain these discrepancies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Design of a shipping fixture for a compact cryomodule hermetic assembly

Two conduction-cooled 915 MHz superconducting radio frequency hermetic assemblies must be safely transported from the Jefferson Lab in Newport News, VA to General Atomics in San Diego, CA for perfor-mance testing in a custom horizontal test cryostat. One hermetic assembly consists of a 2-cell 915 MHz cavity, a coaxial fundamental power coupler, and the warm-to-cold transition beam tubes. The second hermetic assembly consists of a 2-cell 915 MHz cavity only. The assemblies will be transported on a flatbed air-ride trailer over the approximate 4000 km distance. Design requirements included adequate attenuation of 4g vertical axis, 5g beamline axis, and 1.5g lateral axis shock events. The isolation system was designed using helical wire-rope isolators with modal and transient finite element analysis performed in Ansys. Results show shock attenuation of a 10 ms half-sine pulse input to < 1g in the vertical axis, < 1.5g in the beam-line axis, and < 0.5g in the lateral axis for both assem-blies at the specified design loads and all structural stresses are kept below the material yield limits. Addi-tionally, the natural frequencies of both isolation sys-tems adequately attenuate the fundamental modes of the critical structures.

Accelerator Physics↗
DOE OSTI2022

Collision of localized shocks in AdS 5 as a series expansion in transverse gradients

We introduce a computational framework to more efficiently calculate the collision of localized shocks in five dimensional asymptotically Anti-de Sitter space. We expand the Einstein equations in transverse gradients and find that our numerical results agree well with exact solutions already at first order in the expansion. Moreover, the Einstein equations at first order in transverse gradients can be decoupled into two sets of differential equations. Here, the bulk fields of one of these sets has only a negligible contribution to boundary observables, such that the computation on each time slice can be simplified to the solution of several planar shockwave equations plus four further differential equations for each transverse plane ‘pixel’. At the cost of errors of ≲ 10% at the hydrodynamization time and for low to mid rapidities, useful numerical solutions can be sped up by roughly one order of magnitude.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Evaluation of RMI Ejecta Sourcing from Reconstructed Interfaces in 2D FLAG Flyer-Plate Simulations

This study investigates the performance of RMI ejecta sourcing from a reconstructed interface in two-dimensional FLAG ejecta flyer-plate simulations. Enabling the RMI ejecta model to operate on reconstructed interfaces has been a longstanding challenge in FLAG, limiting its applicability to experiments in which interface reconstruction is not required. Addressing this limitation is therefore an important step toward expanding the model’s predictive capability. The first phase of this effort focused on identifying a flyer-plate/target configuration that produces large ejecta clouds from two distinct shocks using the standard RMI ejecta model. A series of simulations examined the effects of flyer composition, layer thickness, and mesh resolution on shock structure and ejecta formation. VISAR traces were used to evaluate shock timing and identify configurations that produced a clear double-shock response. A copper middle-layer flyer with a doubled tin layer yielded the most favorable shock characteristics and was selected as the baseline configuration. After confirming this configuration performed as expected when running a reference FLAG version, it became the basis for three simulations: (1) a standard two-dimensional ejecta flyer-plate calculation using a reference FLAG version, (2) the same calculation with the reconstructed interface option enabled, and (3) a calculation employing ALE strategies to force the ejecta model to operate on the reconstructed interface. While the baseline simulation produced the expected ejecta response, the reconstructed interface approach did not result in measurable ejecta, indicating that this capability remains under development. These results highlight both the progress made and the remaining challenges in extending the RMI ejecta model to reconstructed interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Energy Flow in Thin Shell Implosions and Explosions

Energy flow and balance in convergent systems beyond petapascal energy densities controls the fate of late-stage stars and the potential for controlling thermonuclear inertial fusion ignition. Timeresolved x-ray self-emission imaging combined with a Bayesian inference analysis is used to describe the energy flow and the potential information stored in the rebounding spherical shock at 0.22 petaPascal (2.2 Gbar or billions of atmospheres pressure). This analysis, together with a simple mechanical model, describes the trajectory of the shell and the time history of the pressure at the fuel-shell interface, ablation pressure, and energy partitioning including kinetic energy of the shell and internal energy of the fuel. The techniques used here provide a fully self-consistent uncertainty analysis of integrated implosion data, a thermodynamic-path independent measurement of pressure in the petaPascal range, and can be used to deduce the energy flow in a wide variety of implosion systems to petapascal energy densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

First-principles equation of state of CHON resin for inertial confinement fusion applications

A wide-range (0 to 1044.0 g/ cm 3 and 0 to 10 9 K) equation-of-state (EOS) table for a CH 1.72 O 0.37 N 0.086 quaternary compound has been constructed based on density-functional theory (DFT) molecular-dynamics (MD) calculations using a combination of Kohn-Sham DFT MD, orbital-free DFT MD, and numerical extrapolation. The first-principles EOS data are compared with predictions of simple models, including the fully ionized ideal gas and the Fermi-degenerate electron gas models, to chart their temperature-density conditions of applicability. The shock Hugoniot, thermodynamic properties, and bulk sound velocities are predicted based on the EOS table and compared to those of C-H compounds. The Hugoniot results show the maximum compression ratio of the C-H-O-N resin is larger than that of CH polystyrene due to the existence of oxygen and nitrogen; while the other properties are similar between CHON and CH. Radiation hydrodynamic simulations have been performed using the table for inertial confinement fusion targets with a CHON ablator and compared with a similar design with CH. Here, the simulations show CHON outperforms CH as the ablator for laser-direct-drive target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Theoretical predictive screening of noble-metal-containing M 3 AuC 2 (M = Ti, V, and Cr) MAX phases

Using first-principles density functional theory (DFT), the ground state physical properties of the newly synthesized noble-metal-containing Ti 3 AuC 2 MAX phase have been investigated. The effect of transition-element Ti replacement (V and Cr) on physical properties, including structural, elastic, electronic, thermal, and optical are presented. The optimized lattice parameters of Ti 3 AuC 2 MAX phase are in good agreement with the experimental values and decrease when we replace Ti with V and then Cr. The magnetic properties of M 3 AuC 2 are predicted by using generalized gradient approximation (GGA) incorporating onsite electron correlation parameter (Hubbard U parameter) within DFT. Out of three studied MAX phases, V 3 AuC 2 and Cr 3 AuC 2 phases exhibit magnetism. The thermodynamical stability of M 3 AuC 2 (M = Ti, V, and Cr) phases is discussed by the formation enthalpy with respect to their most competing phases. The electronic structure reveals that these phases have metallic nature and are electrically anisotropic. The M-element replacement has an effect on the bonding properties of M 3 AuC 2 . The bonding between M-C is in the order of Cr 3 AuC 2 > V 3 AuC 2 > Ti 3 AuC 2 according to their peak positions and heights of density of states in the occupied site, which is also confirmed from the charge density distribution. The elastic properties indicate that M 3 AuC 2 phases are ductile, machinable, less stiff, and better resistant to thermal shock and are elastically anisotropic. To understand the properties of M 3 AuC 2 for the extreme environment, we employed a quasi-harmonic Debye model at the pressure and temperature range of 0–50 GPa and 0–1600 K, respectively. Additionally, we evaluated thermal conductivity and melting temperature to further understand the potential of these MAX phases in coating applications at elevated temperature. In relation to optical properties, these MAX phases have a reasonable absorption coefficient in the visible as well as in the ultra-violet regions. The reflectivity of M 3 AuC 2 phases is polarization dependent and of 45 % against the visible region, which reveals its potential as coating material to minimize solar heating.

312 MAX phases↗
DOE OSTI2020

Role of opacity at the 9 keV back lighter energy used in measuring the equation of state of boron at pressures up to a Gbar

Many experiments have been conducted at the National Ignition Facility to measure the equation of state Hugoniot of plastic, boron, and diamond at extreme pressures up to a Gbar. The “Gbar” design employs a strong spherically converging shock launched through a solid ball of material using a hohlraum radiation drive. The shock front conditions are characterized using X-ray radiography, typically at energies near 9 keV. Here we examine how the opacity of boron at 9 keV changes at high pressures and temperatures. Understanding this is vital to unfolding the density in the shock front as pressures exceed 100 Mbar. We compare opacity calculations from a number of methods including the legacy XSN opacity tables, super transition array (STA), average atom (AVAT), and detailed configuration accounting (DCA) methods. We examine how the changing opacity is correlated with the K-shell occupation calculated using the opacity methods as well as other electronic structure methods such as MECCA, Purgatorio, and an all-electron Kohn-Sham density functional theory calculation that uses optimized all-electron norm-conserving Vanderbilt pseudopotentials (ONCV). We also examine the relative contribution of free-free opacity to the overall opacity, which is dominated by bound-free absorption.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Radiation Belt Response to Fast Reverse Shock at Geosynchronous Orbit

Fast reverse shocks (FRSs) cause the magnetosphere to expand, by contrast to the well-known compressions caused by the impact of fast forward shocks (FFS). Usually, FFSs are more geoeffective than FRSs, and consequently the inner magnetosphere dynamic responses to both shock types can be quite different. In this study, we investigate for the first time the radiation belt response to an FRS impact using multi-satellite observations and numerical simulations. Spacecraft on the dayside observed decreases in magnetic field strength and energetic (~40–475 keV) particle fluxes. Timing analysis shows that the magnetic field signature propagated from the dayside to the nightside magnetosphere. Particles with different energies vary simultaneously at each spacecraft, implying a non-dispersive particle response to the shock. Spacecraft located at lower L-shells did not record any significant signatures. The observations indicate a local time dependence of the response associated with the shock inclination, with the clearest signatures being observed in the dusk–midnight sector. Simulations underestimate the amplitude of the magnetic field variations observed on the nightside. The observed decreases in the electron intensities result from a combination of radial gradient and adiabatic effects. The radial gradients in the spectral index appear to be the dominant contributor to the observed variations of electrons seen on the dayside (near noon and dusk) and on the nightside (near midnight). This work shows that even an FRS can affect the radiation belts significantly and provides an opportunity to understand their dynamic response to a sudden expansion of the magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Subnanosecond phase transition dynamics in laser-shocked iron

Iron is one of the most studied chemical elements due to its sociotechnological and planetary importance; hence, understanding its structural transition dynamics is of vital interest. By combining a short pulse optical laser and an ultrashort free electron laser pulse, we have observed the subnanosecond structural dynamics of iron from high-quality x-ray diffraction data measured at 50-ps intervals up to 2500 ps. We unequivocally identify a three-wave structure during the initial compression and a two-wave structure during the decaying shock, involving all of the known structural types of iron (α-, γ-, and ε-phase). In the final stage, negative lattice pressures are generated by the propagation of rarefaction waves, leading to the formation of expanded phases and the recovery of γ-phase. Our observations demonstrate the unique capability of measuring the atomistic evolution during the entire lattice compression and release processes at unprecedented time and strain rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗