Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock loading”

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

Directional amorphization of covalently-bonded solids: A generalized deformation mechanism in extreme loading

Shock compression subjects materials to a unique regime of high quasi-hydrostatic pressure and coupled shear stresses for durations on the order of 1–10 nanoseconds for laser-driven loading of samples. There is, additionally, an attendant temperature increase due to the shock and the mechanisms of plastic deformation in metals whereby dislocations, twins, and phase transitions nucleate and propagate at velocities near the sound speed. Covalently bonded materials have, by virtue of the directionality of their bonds, great difficulty in responding by conventional plastic deformation to this extreme regime of shock compression. Here, we propose that the shear from shock compression induces amorphization, as observed in Si, Ge, B4C, SiC, and olivine ((Mg, Fe) 2 SO 4 ) and that this is a general deformation mechanism in a broad class of covalently bonded materials. The crystalline structure transforms to amorphous along regions of maximum shear stress, forming nanoscale bands, and thereby relaxing the shear component of the imposed shock stress. This process is usually preceded by the emission and propagation of a critical concentration of dislocations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The role of pre-existing heterogeneities in materials under shock and spall

There has been a challenge for many decades to understand how heterogeneities influence the behavior of materials under shock loading, eventually leading to spall formation and failure. Experimental, analytical, and computational techniques have matured to the point where systematic studies of materials with complex microstructures under shock loading and the associated failure mechanisms are feasible. This is enabled by more accurate diagnostics as well as characterization methods. As interest in complex materials grows, understanding and predicting the role of heterogeneities in determining the dynamic behavior becomes crucial. Early computational studies, hydrocodes, in particular, historically preclude any irregularities in the form of defects and impurities in the material microstructure for the sake of simplification and to retain the hydrodynamic conservation equations. Contemporary computational methods, notably molecular dynamics simulations, can overcome this limitation by incorporating inhomogeneities albeit at a much lower length and time scale. This review discusses literature that has focused on investigating the role of various imperfections in the shock and spall behavior, emphasizing mainly heterogeneities such as second-phase particles, inclusions, and voids under both shock compression and release. Pre-existing defects are found in most engineering materials, ranging from thermodynamically necessary vacancies, to interstitial and dislocation, to microstructural features such as inclusions, second phase particles, voids, grain boundaries, and triple junctions. This literature review explores the interaction of these heterogeneities under shock loading during compression and release. Systematic characterization of material heterogeneities before and after shock loading, along with direct measurements of Hugoniot elastic limit and spall strength, allows for more generalized theories to be formulated. Further, continuous improvement toward time-resolved, in situ experimental data strengthens the ability to elucidate upon results gathered from simulations and analytical models, thus improving the overall ability to understand and predict how materials behave under dynamic loading.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dynamic experiments to study the α - ε phase transition in cerium

The ability to understand and predict the response of matter at extreme conditions requires knowledge of a material’s equation-of-state including the location of phase boundaries, transition kinetics, and the evolution of material strength. Cerium is a material with a complex phase diagram that continues to attract significant scientific interest. Recent dynamic experiments have provided information on the low-pressure γ–α phase transition, sound speed, and Hugoniot data for the higher-pressure α phase, as well as the incipient shock melt transition. Despite these efforts, there are still regions of the phase diagram that are largely unexplored dynamically, including the high-pressure region below the melt boundary. Along a room temperature isotherm, diamond anvil cell data report a transition to the ϵ phase between 13 and 17 GPa. At higher temperatures, similar diamond anvil cell data show significant disagreement regarding the existence, location, and slope of the ϵ-phase boundary. In this work, double-shock loading was used to access the α–ϵ region of the phase diagram to obtain equation-of-state information and to determine the location of the ϵ-phase boundary for shock loading.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Effect of local structure and stoichiometry on the dynamic behavior of bi-metal interfaces

Molecular dynamics (MD) simulations are performed to investigate the failure behavior of Al/Al, Al/Fe, and Al/AlFe bi-layers under high strain-rate shock loading. Additional simulations are performed to invesigate solely the tensile response of these systems and to isolate the effects of loading history on damage and failure. Specifically, the compression stage of the shock is excluded from the tensile simulations. This study shows that local structure variation, including the introduction of serrations, plays different roles in controlling damage nucleation depending on the specific loading condition. Under shock loading, the stress for void nucleation and the resulting void distribution is insensitive to an interface structure, whereas under pure tensile loading, the opposite is true, even though the strain rate is comparable. The above difference can be explained based on the corresponding change in the total deformation prior to damage nucleation under shock loading, which is missing from the pure tensile loading. Therefore, it is concluded that whether local interface structure matters under high strain-rate loading depends on the specific loading history and the stress state: the local structure variation matters only when it alters the deformation behavior, and subsequently, damage nucleation and evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The role of micro-inertia on the shock structure in porous metals

The behavior of porous materials under shock loading is a multi-scale problem bridging orders of magnitude across the macroscale geometry and the microscale pores. Under static loading, this problem is well understood, relating mechanisms of pore closure and crushing to the equivalent macroscale models. The dynamic response of porous solids under shock loading is related to the effects of viscoplasticity and micro-acceleration fields around the void boundaries. The significance of the micro-inertia effects in modeling the dynamic behavior of porous materials remains an open question. In this work, an experimental investigation on closed-cell porous aluminum with small porosity provides the evidence for the first time of micro-inertia’s fundamental role in describing the shock structure in these materials. Materials with different levels of porosity were manufactured using a modified process of additive manufacturing to achieve a mean pore size below 50$\mu$m. Plate impact experiments on porous aluminum samples were conducted at pressures in the range of 2 to 11 GPa. The structure of the steady shock was characterized as a function of porosity and shown to validate behavior revealed by an analytical approach (Czarnota et al. [J. Mech. Phys. Solids 107 (2017)]), highlighting the fundamental role of micro-inertia effects in such cases.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Viscoelastic Modeling of Polymers in ALEGRA with the GAP Model

The Glassy Amorphous Polymer (GAP) model is a viscoelastic/plastic model developed at Los Alamos National Laboratory to accurately model a variety of polymers across a wide range of conditions and loading rates, including shock loading. In the present report we introduce and assess this model, newly implemented in the ALEGRA shock and multiphysics code, using a series of verification and application-related validation problems. We describe the mathematical and theoretical formulation of the model, as well as its implementation in ALEGRA, in detail. We provide verification results that assess the model implementation against published computational results, as well as validation results which we compare to existing experimental results when possible. These comparisons instill confidence in the implementation and indicate that the addition of the GAP model to the ALEGRA code provides users with a high-fidelity polymer modeling capability that is capable of recreating complex polymer phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dislocation generation in diamond under extreme loading

Diamond is, by virtue of its high bonding forces and Peierls-Nabarro barrier stresses, among the hardest materials on earth due to the difficulty of generating and moving dislocations. We demonstrate, using molecular dynamics and analytical calculations, that the generation of defects is dependent on loading orientation. Shock-loading single-crystal diamond along [001] and [011] to a stress of 137 GPa did not reveal dislocations, whereas loading along [111] generated profuse dislocations. The introduction of a void at the latter orientation generated dislocations at 72.6 GPa. Two slip systems were identified: <011>{100} and <112>{111}. The threshold for plastic deformation in diamond is orientation dependent and significantly lower than reported in previous studies: 0.14 G (shear modulus). The void collapse generated localized amorphization for [001]. In conclusion, these results on the anisotropy of plasticity and void effects have relevance for the improvement of symmetry in the collapse of diamond capsules in inertial fusion experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

High-pressure and temperature neural network reactive force field for energetic materials

Reactive force fields for molecular dynamics have enabled a wide range of studies in numerous material classes. These force fields are computationally inexpensive compared with electronic structure calculations and allow for simulations of millions of atoms. However, the accuracy of traditional force fields is limited by their functional forms, preventing continual refinement and improvement. Therefore, we develop a neural network-based reactive interatomic potential for the prediction of the mechanical, thermal, and chemical responses of energetic materials at extreme conditions. The training set is expanded in an automatic iterative approach and consists of various CHNO materials and their reactions under ambient and shock-loading conditions. Further, this new potential shows improved accuracy over the current state-of-the-art force fields for a wide range of properties such as detonation performance, decomposition product formation, and vibrational spectra under ambient and shock-loading conditions.

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

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase Transitions in (Mg,Fe) 2 SiO 4 Olivine under Shock Compression

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

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spatial description of dislocation nucleation in the shock response of single-crystal aluminum

Nonequilibrium molecular dynamics simulations of shock loaded single-crystal Al in the $\langle$100$\rangle$, $\langle$110$\rangle$, $\langle$111$\rangle$, and $\langle$123$\rangle$ orientations are conducted to study elastic and plastic shockwave formation and details associated with dislocation activity. A computer vision-based approach is implemented to capture the presence of dislocations and describe their spatial characteristics in the zone of nucleation behind the propagating shockwave. The methodology developed relies on the sequences of images extracted during shock loading that show dislocation activity within a cross section of the sample. Results reveal that the spacing between activated slip systems is orientation dependent and exhibits a modest reduction for the $\langle$100$\rangle$ and $\langle$111$\rangle$ orientations as shock pressure increases. Comparisons are made to existing theoretical models. Such relationships between shock pressure and dislocation activity, extracted from molecular dynamics simulations, can be used to inform higher length scale simulations or modeling of dislocation-based plasticity during shock.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Investigation of novel twinning in additive body centered cubic titanium (Fiscal Year 2024 Final Report)

This report summarizes the result of a one year seedling project to investigate unusual twinning behavior in shock loaded additive Ti5552. The twinning behavior only occurs when the β phase of Ti5553 is metastable, and it appears to be a type of double twin involving two different twin variants, first a {332}⟨113⟩ twin forms before being consumed by a specific {112}⟨111⟩ twin variant to create a 20 o . This behavior has only been detected during shock loading around incipient spall damage. The twinning is investigated by performing postmortem EBSD and PED analysis of gas-gun loaded specimens and preliminary molecular dynamics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Probabilistic Predictions for Fastener Failure in the Sandia Mechanics Challenge Using the Discrete-Direct Uncertainty Quantification Approach

This paper documents the blind and post-blind analysis predictions for the 2023 Sandia Mechanics Challenge (SMC), which involved predicting the behavior of a threaded fastener joint structure subjected to shock loading. Utilizing repeat sets of fastener calibration data from various experimental configurations including tension, double shear, and joint tension, we developed a library of calibrated models which were propagated through the application model using the Discrete-Direct (DD) uncertainty quantification (UQ) approach. Although the initial blind predictions did not incorporate spare-sample processing to quantify fastener failure probabilities, the analyses yielded reasonable conclusions aligned with experimental results. In the post-blind analysis phase, we focused on enhancing the fidelity of the aluminum constitutive model and innovating the DD approach to obtain probabilistic predictions for fastener failure, particularly when quantities of interest (QoIs) approach their bounds. The improved aluminum model captures the behavior of the cantilever under shock loading more accurately, predicting both partial and complete cracks, although it tends to underpredict failure propagation. The enhanced DD approach facilitates probabilistic predictions that reflect the interdependent failure mechanisms of the fasteners and the cantilever, revealing that while certain fasteners are more likely to fail, the failure does not necessarily follow a progressive pattern. Overall, the post-blind analyses significantly improved the predictive capabilities of the model, providing valuable insights into the SMC application and establishing a robust foundation for informed engineering decisions. The methodology demonstrates a cost-effective and extensible approach suitable for a wide range of applications, highlighting the importance of uncertainty quantification to provide context for engineering decision making.

42 ENGINEERING↗
DOE OSTI2022

Role of α → ε → α phase transformation on the spall behavior of iron at atomic scales

Shock compression of iron microstructures above a threshold stress results in a α (BCC) → ε (HCP) transformation and the propagation of the release wave brings the metal back to the α phase following the ε→α transformation. Predicting failure behavior under shock loading conditions (spallation) relies on understanding the evolution of defects in the microstructure as it undergoes the α→ε→α phase transformation. This study uses molecular dynamics (MD) simulations to investigate the role of defect evolution during α→ε→α phase transformation on the spall strength values of single-crystal (sc) Fe microstructures. The MD simulations aim to characterize the ε phase fraction formed during shock compression and the defects during shock release for variations in loading orientations and shock stresses. The simulations are carried out for loading along the [100], [110], [111], and [112] orientations and for impact velocities ranging from 600 m/s to 1 km/s. The ε phase fractions during compression and defects (dislocations, twins) characterized during spall failure show an orientation dependence that affects the spall strength values. Here, the lowest value for spall strength is observed for the <110> loading orientation that shows a high density of twinning at the spall plane, whereas the highest value is observed for the <100> orientation and is associated with a α (BCC) → γ (FCC) transformation at the spall plane. The correlations of the spall strength values with the strain rates and with the ε phase fractions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computational Design of a Simple Flyer Plate Launcher

The response of materials to shock loading is important to understand for a variety of applications. When shock physics emerged during and after WWII, direct explosive loading or explosively driven plate impact was the primary tool for these studies. Subsequent decades have seen the widespread use of large caliber guns for plate impact studies, laser-shock facilities and pulsed power facilities. INL currently lacks a gun suitable for plate impact or explosives casting and machining facilities; however, it does possess explosives use and handling capabilities. An option for performing plate impact experiments was needed, therefore continuum scale models were utilized to explore a few simple donor-acceptor explosive plane wave lens designs, one of which could be hand packed with plastic explosives to launch flyer plates. 2D simulations were performed to study different geometries in an effort to minimize the difference in shock arrival across the central portion of a small copper flyer plate. A shock wave arrival time difference under 50 ns across 50 percent of the center of the flyer was achieved with a few designs. This work summarizes the computational models and results.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bounds on the Rate-Dependent Plastic Flow of Tantalum up to 75 GPa

In this study we report on the continued development of thermodynamics-based analysis of shock waves propagation with the objective of extracting information related to materials strength at high strain rates and pressures. Building on previous results reported for peak stresses of 10 GPa and 25 GPa, we present a series of three-step gas-gun shock experiments designed to explore the pressure and strain rate dependence of plastic flow in polycrystalline tantalum. These experiments at nominal peak stresses of 50 GPa and 75 GPa show the irreversible deformation before pullback to be almost entirely confined to the shock loading, with negligible plastic relaxation on the post-shock plateau. We also add a reverse-ballistics shot at 25 GPa, which was designed to reveal the pullback response with negligible interference from free-surface effects. General thermodynamic considerations allow us to place bounds on the plastic behavior even for parts of the curve that change far too rapidly for the velocimetric time resolution of (conservatively) ~ 5 ns. To analyze the data, we found it necessary to substantially improve the interpolation/extrapolation scheme in order to improve its robustness, flexibility and range of applicability. We describe the new scheme based on splines, as well as an extension of free-surface corrections to the post-shock rarefaction waves. Reanalysis with the new scheme produces results essentially within the error bars previously reported, showing that the known systematic errors associated with free-surface effects are relatively inconsequential for determining thermodynamic paths.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Microphysics of Relativistic Collisionless Electron-ion-positron Shocks

Abstract We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizations σ ≲ 10 −4 and pair-loading factors Z ± ≲ 10 are studied, where Z ± is the number of loaded electrons and positrons per ion. We find the following: (1) The shock becomes mediated by the ion Larmor gyration in the mean field when σ exceeds a critical value σ L that decreases with Z ± . At σ ≲ σ L the shock is mediated by particle scattering in the self-generated microturbulent fields, the strength and scale of which decrease with Z ± , leading to lower σ L . (2) The energy fraction carried by the post-shock pairs is robustly in the range between 20% and 50% of the upstream ion energy. The mean energy per post-shock electron scales as E ¯ e ∝ Z ± + 1 − 1 . (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low as σ ≈ 5 × 10 −6 . The ions then become essentially thermal with mean energy E ¯ i , while electrons form a nonthermal tail, extending from E ∼ Z ± + 1 − 1 E ¯ i to E ¯ i . When σ = 0, particle acceleration is enhanced by the formation of intense magnetic cavities that populate the precursor during the late stages of shock evolution. Here, the maximum energy of the nonthermal ions and electrons keeps growing over the duration of the simulation. Alongside the simulations, we develop theoretical estimates consistent with the numerical results. Our findings have important implications for models of early gamma-ray burst afterglows.

79 ASTRONOMY AND ASTROPHYSICS↗