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 829 records · Page 46

Intergranular Hotspots: A Molecular Dynamics Study on the Influence of Compressive and Shear Work

Numerous crystal- and microstructural-level mechanisms are at play in the formation of hotspots, which are known to govern high explosives initiation behavior. Most of these mechanisms, including pore collapse, interfacial friction, and shear banding, involve both compressive and shear work done within the material and have thus far remained difficult to separate. Here, we assess hotspots formed at shocked crystal–crystal interfaces using quasi-1D molecular dynamics simulations that isolate effects due to compression and shear. Two high explosive materials are considered (TATB and PETN) that exhibit distinctly different levels of molecular conformational flexibility and crystal packing anisotropy. Temperature and intramolecular strain energy localization in the hotspot are assessed through parametric variation of the crystal orientation and two velocity components that respectively modulate compression and shear work. The resulting hotspots are found to be highly localized to a region within 5–20 nm of the crystal–crystal interface. Compressive work plays a considerably larger role in localizing temperature and intramolecular strain energy for both materials and all crystal orientations considered. Shear induces a moderate increase in energy localization relative to unsheared cases only for relatively weak compressive shock pressures of approximately 10 GPa. These results help isolate and rank the relative importance of hotspot generation mechanisms and are anticipated to guide the treatment of crystal–crystal interfaces in coarse-grained models of polycrystalline high-explosive materials.

36 MATERIALS SCIENCE↗

Time resolved radiographs and fuel density of the stagnation phase of Inertial Confinement Fusion implosions (NSA report)

This report contains a figure and caption. Left: Compton radiographs of a 3-shock HDC THD layered implosion recorded 10ps before and 160ps after bang-time (BT) at the NIF, using >50keV point projection backlighters generated by irradiating a 25μm-diameter Au wire with 30ps-long ARC pulses. Right: Reconstructed density distributions of fuel and remaining ablator showing significant 3D structure and top-down asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of Hydrogen Fuel Cell Technology for Freight Rail Propulsion and Review of Relevant Industry Standards.

Alternatives to conventional diesel electric propulsion are currently of interest to rail operators. In the U.S., smaller railroads have implemented natural gas and other railroads are exploring hydrogen technology as a cleaner alternative to diesel. Diesel, battery, hydrogen fuel cell, or track electrification all have trade-offs for operations, economics, safety, and public acceptability. A framework to compare different technologies for specific applications is useful to optimize the desired results. Standards from the Association of American Railroads (AAR) and other industry best practices were reviewed for applicability with hydrogen fuel cell technology. Some technical gaps relate to the physical properties of hydrogen, such as embrittlement of metals, invisible flames, and low liquid temperatures. A reassessment of material selection, leak/flame detection, and thermal insulation methods is required. Hydrogen is less dense and diffuses more easily than natural gas, and liquid hydrogen is colder than liquefied natural gas. Different densities between natural gas and hydrogen require modifications to tank designs and flow rates. Leaked hydrogen will rise rather than pool on the ground like diesel, requiring a modification to the location of hydrogen tanks on rolling stock. Finally, the vibration and shock experienced in the rail environment is higher than light-duty vehicles and stationary applications for which current fuel cell technology has been developed, requiring a modification in tank design requirements and testing.

08 HYDROGEN↗

3D xRAGE simulation of inertial confinement fusion implosion with imposed mode 2 laser drive asymmetry

Low-mode asymmetries represent an important obstacle to achieving high-gain inertial confinement fusion implosions. As a step in learning how to control such effects, an OMEGA experiment with imposed mode 2 laser drive asymmetries was done to study the expected signatures of this type of asymmetry [M. Gatu Johnson et al., PRE 2018]. In the present work, a 3D xRAGE simulation including the stalk mount has been brought to bear on the data from that experiment. Furthermore, comprehensive comparisons between simulated and measured observables are made. Good agreement between simulated and measured x-ray image-inferred shell trajectories, bang times and neutron emission widths are seen, showing that the hydrodynamics are well captured in the simulation. Asymmetries seen in simulated and measured time-resolved and time-integrated x-ray images and areal densities also compare well, showing impact of both stalk and mode 2. On the other hand, important differences in measured and simulated neutron emission histories, yield, and ion temperature (T ion ) asymmetries are seen, suggesting that the simulation is overestimating shock yield. Overall, the results clearly demonstrate the importance of considering all asymmetry sources when interpreting measured signatures of asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fuel–shell mix and yield degradation in kinetic shock-driven inertial confinement fusion implosions

Fuel–shell mix in kinetic plasma conditions is probed using nuclear and x-ray self-emission in shock-driven, D 3 He-gas-filled inertial confinement fusion implosions. As initial gas fill density decreases, measured nuclear yields and ion temperatures are lower than expected as compared to radiation-hydrodynamic simulations. Spatially and temporally resolved x-ray emissions indicate significant mixing at the fuel–shell interface in implosions with low initial gas fill density. In conclusion, this observed fuel–shell mix introduces a substantial amount of shell ions into the center of the implosion prior to and during shock flash and is the key mechanism needed in the kinetic-ion simulations to match experimental nuclear yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shock formation in flowing plasmas by temporally and spatially smoothed laser beams

The cumulative impact of multiple laser speckles on a supersonic plasma flow across optically smoothed laser beams is investigated. The bending of laser beams caused by ponderomotive laser–plasma coupling, together with flow, leads to plasma a momentum-conserving response that results in a deceleration of the flow. Once the flow velocity decreases to a subsonic level, the action of the laser beams can generate a shock within the plasma. This scenario has been predicted theoretically and confirmed by hydrodynamic simulations. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size, and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging only in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser–plasma coupling in the inertial confinement fusion experiments are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evidence of non-Maxwellian ion velocity distributions in spherical shock-driven implosions

The ion velocity distribution functions of thermonuclear plasmas generated by spherical laser direct drive implosions are studied using deuterium-tritium (DT) and deuterium-deuterium (DD) fusion neutron energy spectrum measurements. A hydrodynamic Maxwellian plasma model accurately describes measurements made from lower temperature (< 10 keV), hydrodynamiclike plasmas, but is insufficient to describe measurements made from higher temperature more kineticlike plasmas. The high temperature measurements are more consistent with Vlasov-Fokker-Planck (VFP) simulation results which predict the presence of a bimodal plasma ion velocity distribution near peak neutron production. Furthermore, these measurements provide direct experimental evidence of non-Maxwellian ion velocity distributions in spherical shock driven implosions and provide useful data for benchmarking kinetic VFP simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fermi-LAT Gamma-Ray Emission Discovered from the Composite Supernova Remnant B0453-685 in the Large Magellanic Cloud

Abstract We report the second extragalactic pulsar wind nebula (PWN) to be detected in the megaelectronvolt–gigaelectronvolt band by the Fermi-LAT, located within the Large Magellanic Cloud. The only other known PWN to emit in the Fermi band outside of the Milky Way is N157B, which lies to the west of the newly detected gamma-ray emission at an angular distance of 4°. Faint, pointlike gamma-ray emission is discovered at the location of the composite supernova remnant (SNR) B0453-685 with a ∼4 σ significance with energies ranging from 300 MeV–2 TeV. We present the Fermi-LAT data analysis of the new gamma-ray source, coupled with a detailed multiwavelength investigation to understand the nature of the observed emission. Combining the observed characteristics of the SNR and the physical implications from broadband modeling, we argue it is unlikely that the SNR is responsible for the gamma-ray emission. While the gamma-ray emission is too faint for a pulsation search, we try to distinguish between any pulsar and PWN component of SNR B0453-685 that could be responsible for the observed gamma-ray emission using semi-analytic models. We determine the most likely scenario is that the old PWN ( τ ∼ 14,000 yr) within B0453-685 has been impacted by the return of the SNR reverse shock with a possible substantial pulsar component below 5 GeV.

79 ASTRONOMY AND ASTROPHYSICS↗

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems Using MOOSE-based Multiphysics Simulations

As an increasing number of companies look toward commercial Inertial Fusion Energy (IFE) designs, there is a pressing need to understand the physics of thick-wall chamber gas dynamics. The thick liquid wall approach implements a renewable wall to mitigate the fusion target emissions, thereby reducing the radiation damage rate and significantly extending the lifetime of chamber structures, leading to increased plant availability and reduced waste streams in comparison to dry wall chamber designs. It is necessary, however, to assess the critical performance and safety aspects of these systems. For example, it is crucial to predict (1) where the mass ablated from the liquid walls will vent, which determines the placement of condensing surfaces; (2) debris propagation up the beam lines, which provides essential information for design and protection requirements; (3) peak pressures and impulse on chamber walls, which affect chamber structural design; and (4) momentum transfer to the liquid jets, which constrains the shape and positioning of the jets. In turn, the chamber design and its liquid walls affect shielding requirements, material activation, and tritium fuel cycle. Currently available simulation tools, however, are unable to accurately capture key thick-wall chamber dynamics. Significant assumptions are often made to simplify the system and reduce computational cost and modeling capability needs, but the impact of these assumptions on simulation predictions has not been evaluated. For example, no three-dimensional simulations can be found in the open literature to evaluate gas venting and momentum transfer to the jets with simulations using two-dimensional domains to represent complex three-dimensional geometries. Moreover, limited studies have been dedicated to jet breakup due to both turbulence and neutron heating, and no studies have been found that evaluate how jet breakup can impact shock-jet interaction. Furthermore, effects of radiative heat transfer have rarely been included for the hydrodynamic phase of shock propagation, and integration of proper equations of state in shock dynamics codes has been mostly exploratory. In this study, we use the flexible, high-fidelity Multiphysics Object-Oriented Simulation Environment (MOOSE) to model these complex phenomena and inform design and safety studies. Capabilities to model thick-wall chamber gas dynamics are being developed, and the impact of the assumptions listed above (i.e., two-dimensional vs three-dimensional, absence of jet breakout, no radiative heat transfer, and ideal gas behavior) are being quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Steady detonation propagation in thin channels with strong confinement

We examine asymptotically the dynamics of two-dimensional, steady detonation wave propagation and failure for a strongly confined high explosive (HE), in which the width of the explosive is small relative to the reaction zone length. An energy balance equation is derived, which shows how the longitudinal acceleration of subsonic flow behind the detonation shock is influenced both by chemical reaction and by the effects of HE boundary streamline deflection, specifically via the induced rate of change of mass flux through the detonation wave. The latter serves to either counteract or reinforce the acceleration of longitudinal flow, depending on the sign of the gradient of the boundary streamline deflection at the detonation shock. The analysis is valid for general equations of state and chemical reaction rates in the HE. The asymptotically derived form of the energy equation represents an eigenvalue problem for the determination of the steady detonation propagation speed, solved via a shooting method. We explore specific results for ideal and stiffened equations of state, along with a pressure-dependent reaction rate for which changes in the pressure exponent and reaction order are also studied. We consider the influences of both straight and curved HE boundary streamline shapes. The asymptotic analysis reveals significant physical insights into how detonation propagation and failure are affected by strong confinement.

42 ENGINEERING↗

Radiation Damage Studies on Titanium Alloys as High Intensity Proton Accelerator Beam Window Materials

A high-strength dual a+ß phase titanium alloy Ti–6Al–4V is utilized as a material for beam windows in several accelerator target facilities. However, relatively little is known about how material properties of this alloy are affected by high-intensity proton beam irradiation. With plans to upgrade neutrino facilities at J-PARC and Fermilab to over 1 MW beam power, the radiation damage in the window material will reach a few displacements per atom (dpa) per year, significantly above the ~0.3 dpa level of existing data. The international RaDIATE collaboration, Radiation Damage In Accelerator Target Environments, has conducted a high intensity proton beam irradiation of various target and window material specimens at Brookhaven Linac Isotope Producer (BLIP) facility, including a variety of titanium alloys. Post-Irradiation Examination (PIE) of the specimens in the 1st capsule, irradiated at up to 0.25 dpa, is in progress. Tensile tests in a hot cell at Pacific Northwest National Laboratory (PNNL) exhibited a clear signature of radiation hardening and loss of ductility for Ti–6Al–4V, while Ti–3Al–2.5V, with less ß phase, exhibited less severe hardening. Microstructural investigations will follow to study the cause of the difference in tensile behavior between these alloys. High-cycle fatigue (HCF) performance is critical to the lifetime estimation of beam windows exposed to a periodic thermal stress from a pulsed proton beam. The first HCF data on irradiated titanium alloys are to be obtained by a conventional bend fatigue test at Fermilab and by an ultrasonic mesoscale fatigue test at Culham Laboratory. Specimens in the 2nd capsule, irradiated at up to ~1 dpa, cover typical titanium alloy grades, including possible radiation-resistant candidates. These systematic studies on the effects of radiation damage of titanium alloy materials are intended to enable us not only to predict realistic lifetimes of current beam windows made of Ti–6Al–4V, but also to extend the lifetime by choosing a more radiation and thermal shock tolerant alloy with a preferable heat treatment, or even by developing new materials.

Titanium Alloy, Beam Window, Proton Beam, Radiatio↗

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems using MOOSE

As an increasing number of companies look toward commercial Inertial Fusion Energy (IFE) designs, there is a pressing need to understand the physics of thick-wall chamber gas dynamics. The thick liquid wall approach implements a renewable wall to mitigate the fusion target emissions, thereby reducing the radiation damage rate and significantly extending the lifetime of chamber structures, leading to increased plant availability and reduced waste streams in comparison to dry wall chamber designs. It is necessary, however, to assess the critical performance and safety aspects of these systems. For example, it is crucial to predict (1) where the mass ablated from the liquid walls will vent, which determines the placement of condensing surfaces; (2) debris propagation up the beam lines, which provides essential information for design and protection requirements; (3) peak pressures and impulse on chamber walls, which affect chamber structural design; and (4) momentum transfer to the liquid jets, which constrains the shape and positioning of the jets. In turn, the chamber design and its liquid walls affect shielding requirements, material activation, and tritium fuel cycle. Currently available simulation tools, however, are unable to accurately capture key thick-wall chamber dynamics. Significant assumptions are often made to simplify the system and reduce computational cost and modeling capability needs, but the impact of these assumptions on simulation predictions has not been evaluated. For example, no three-dimensional simulations can be found in the open literature to evaluate gas venting and momentum transfer to the jets with simulations using two-dimensional domains to represent complex three-dimensional geometries. Moreover, limited studies have been dedicated to jet breakup due to both turbulence and neutron heating, and no studies have been found that evaluate how jet breakup can impact shock-jet interaction. Furthermore, effects of radiative heat transfer have rarely been included for the hydrodynamic phase of shock propagation, and integration of proper equations of state in shock dynamics codes has been mostly exploratory. In this study, we use the flexible, high-fidelity Multiphysics Object-Oriented Simulation Environment (MOOSE) to model these complex phenomena and inform design and safety studies. Capabilities to model thick-wall chamber gas dynamics are being developed, and the impact of the assumptions listed above (i.e., two-dimensional vs three-dimensional, absence of jet breakout, no radiative heat transfer, and ideal gas behavior) are being quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dijet impact factor in DIS at next-to-leading order in the Color Glass Condensate

In this work, we compute the next-to-leading order impact factor for inclusive dijet production in deeply inelastic electron-nucleus scattering at small x Bj . Our computation, performed in the framework of the Color Glass Condensate effective field theory, includes all real and virtual contributions in the gluon shock wave background of all-twist lightlike Wilson line correlators. We demonstrate explicitly that the rapidity evolution of these correlators, to leading logarithmic accuracy, is described by the JIMWLK Hamiltonian. When combined with the next-to-leading order JIMWLK Hamiltonian, our results for the impact factor improve the accuracy of the inclusive dijet cross-section to \( \mathcal{O} \) ( \( {\alpha}_s^2 \) ln( x f /x Bj )), where x f is a rapidity factorization scale. These results are an essential ingredient in assessing the discovery potential of inclusive dijets to uncover the physics of gluon saturation at the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Next-generation anion exchange membrane water electrolyzers operating for commercially relevant lifetimes

Alkaline anion exchange membrane (AEM) water electrolysis has gained increasing attention due to its potential to achieve low-cost, high performance hydrogen production. However, most existing membranes are not durable in industrial settings. In this work, we demonstrate good performance relative to industrial parameters using Sustainion® anion exchange membranes. Long-duration tests showed stable performance of 1 A/cm 2 at 1.85 V with a degradation rate of less than 1 μV/h over 10,000 h. The projected lifetime is thus over 20 years. Daily on/off cycling performance over the course of 30 years was simulated experimentally through accelerated voltage shock tests, resulting in a performance loss of only 0.15 μV/cycle over 11,000 cycles. As shown through impact and crossover testing, an improvement in performance is achieved by the addition of zirconia to the polymer matrix and mechanically reinforcing the membrane.

08 HYDROGEN↗

Peening Techniques for Mitigating Chlorine-Induced Stress Corrosion Cracking of Dry Storage Canisters for Nuclear Applications

Fusion-welded austenitic stainless steel (ASS) was predominantly employed to manufacture dry storage canisters (DSCs) for the storage applications of spent nuclear fuel (SNF). However, the ASS weld joints are prone to chloride-induced stress corrosion cracking (CISCC), a critical safety issue in the nuclear industry. DSCs were exposed to a chloride-rich environment during storage, creating CISCC precursors. The CISCC failure leads to nuclear radiation leakage. Therefore, there is a critical need to enhance the CISCC resistance of DSC weld joints using promising repair techniques. This review article encapsulates the current state-of-the-art of peening techniques for mitigating the CISCC in DSCs. More specifically, conventional shot peening (CSP), ultrasonic impact peening (UIP), and laser shock peening (LSP) were elucidated with a focus on CISCC mitigation. The underlying mechanism of CISCC mitigation in each process was summarized. Finally, this review provides recent advances in surface modification techniques, repair techniques, and developments in welding techniques for CISCC mitigation in DSCs.

Chemistry↗

Double cylinder implosion experiments at the National Ignition Facility

Here, cylindrical implosion experiments are used to directly measure instability growth in a convergent geometry, providing a wealth of data for model validation. Double cylinders are a natural extension of the platform and enable measurements at a classically unstable interface, the outer surface of the inner cylinder, which experiences no ablative stabilization from the laser drive. However, the utility of this platform relies upon maintaining adequate axial uniformity of the inner cylinder during the implosion. Although previous smaller-scale double cylinder experiments exhibited acceptable levels of axial uniformity, radiation-hydrodynamics simulations of larger-scale double cylinders predict more axial non-uniformity induced by the impedance mismatch as the shock wraps around the axial ends of the inner cylinder. A mechanism to reduce axial non-uniformity in these larger double cylinder implosions is presented, and preliminary experimental data confirms the efficacy of the selected mitigation approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Re-Examination of the N 2 O + O Reaction

The reaction of N 2 O with O is a key step in consumption of nitrous oxide in thermal processes. It has two product channels, NO + NO (R2) and N 2 + O 2 (R3). The rate constant for R2 has been measured both in the forward and the reverse direction at elevated temperature and is well established. However, the rate constant for the N 2 + O 2 channel (R3) has been difficult to quantify and has significant error limits. Here, the direct reaction on the triplet surface has a barrier of around 40 kcal mol –1 , and it is too slow for the N 2 + O 2 channel to have any practical significance. Recently, Pham and Lin (2022) suggested an alternative low activation energy reaction path that involves intersystem crossing and reaction on the singlet surface. In the present work, we re-examined a wide range of experiments relevant for the N 2 O + O reaction through kinetic modeling, paying attention to the impact of artifacts such as impurities and surface reactions. Experimental results from shock tubes and batch reactors on the final NO yield in N 2 O decomposition, covering temperatures of 973–2200 K and pressures of 0.013–11.5 atm, support k 3 ~ 0, consistent with the high activation energy for reaction on the triplet surface and a low probability of ISC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗