Search NASA⌕ Search

SEARCH · Search NASA

Results for “shock deformation”

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 73 records · Page 4

Numerical investigation of laser-driven shock interaction with a deformable particle

A laser-driven shock propagating through an isolated particle embedded in a plastic (CH) target was studied using the radiation-hydrodynamic code FLASH. Preliminary simulations using IONMIX equations of state (EOS) showed significant differences in the shock Hugoniot of aluminum compared to experimental data in the low-pressure regime [O(10) GPa], resulting in higher streamwise compression and deformation of an aluminum particle. Hence, a simple modification to the ideal gas EOS was developed and employed to describe the target materials and examine the particle dynamics. The evolution of the pressure field demonstrated a complex wave interaction, resulting in a highly unsteady particle drag which featured two drag minima due to shock focusing at the rear end of the particle and rarefaction stretching due to laser shut-off. Although ~30% lateral expansion and ~25% streamwise compression were observed, the aluminum particle maintained considerable integrity without significant distortion. Additional simulations examined the particle response for a range of particle densities, sizes, and acoustic impedances. The data revealed that lighter particles such as aluminum gained significant momentum, reaching up to ~96% of the shocked CH's speed, compared to ~29% for the heavier tungsten particles. Despite the differences seen in the early stage of shock interaction, particles with varying acoustic impedances ultimately reached the same peak velocity. This identified particle-to-host density ratio is an important factor in determining the inviscid terminal velocity of the particle. In addition, the modified EOS model presented in this study could be used to approximate solid materials in hydrocodes that lack material strength models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗

Electron petrography of Apollo 14 and 15 breccias and shock-produced analogs

Six lunar breccias were studied in detail by optical and transmission electron microscopy and separated into two classes on the basis of presence or absence of recovery and recrystallization textures: Class A breccias with no recrystallization include samples 14063, 14301, and 15498; Class B breccias with recrystallization include samples 14321, 15455, and 15418. Samples of Apollo 16 soil (65010) experimentally shocked to 50, 100, and 250 kilobars pressure were sectioned and studied for comparison with the breccias. It is concluded that the Class A breccias were formed from unconsolidated lunar regolith by deformation due to shock waves of a few tens of kilobars. The Class B breccias could also have been produced from unconsolidated material at higher shock stresses, the annealing being achieved by the heating resulting from the shock deformation.

Christie, J. M.↗

Mineralogical Analysis of Calcium-Aluminum-Rich Inclusions Provides Insight Into Post-Formation Processes

Introduction: Calcium-aluminum inclusions (CAIs) are cm- to mm-sized intergrowths of refractory phases found in chondritic meteorites [1]. Their mineral compositions closely match the compositions of the solids thought to condense from an extremely hot (>1500K) gas with a bulk solar composition [2-4], suggesting that CAIs were the earliest solids within the Solar System [2-4]. These inclusions provide invaluable insights into the conditions and dynamics of the early Solar System. CAIs must be transported from their formation region near the protosun to the chondrite parent-body accretion region. The nature of their journey could affect how, when, and where the secondary processes recorded in these CAIs occurred [4]. Did these secondary processes occur in the solar nebula or during accretion with the parent body, or both? How were the CAIs affected by varying thermochemical processes? To better answer these questions, we have undertaken an in-depth, textural study of the secondary alteration of select CAI samples. Methods: We chose three CAIs that experienced a range of post-formation processing to gain better understanding of secondary alteration based on previous preliminary examination [5]. Representative CAIs were chosen from NWA 5508 (CV3), NWA 12772 (CV3), and Coolidge (CL4) carbonaceous chondrites. We used scanning electron microscopy (SEM) and electron backscatter diffraction (EBSD) techniques for detailed chemical, mineralogical, and textural analysis of the CAIs. Backscattered electron (BSE) images and X-ray elemental maps were taken of Coolidge and NWA 5508 using the Lunar and Planetary Institute (LPI) Phenom SEM. The JEOL 7900F SEM at NASA Johnson Space Center (JSC) was used to obtain energy dispersive spectroscopy (EDS) chemical maps of all samples. High resolution EBSD analysis identified the mineral phases and textures, providing key information on their nature. Based on the SEM and EBSD data, minerals of interest were selected for quantitative chemical electron probe micro-analysis (EPMA) using the JEOL JXA-8530F at NASA JSC. Results and Discussion: The EDS maps show that the NWA 5508 and NWA 12772 CAIs designated “Saguaro” and “Hoopoe” respectively [5] are enriched in calcium while the Coolidge CAI designated “Cottonwood” is aluminum and magnesium rich. Saguaro. A ~1.5cm diameter igneous Type B CAI with a rounded shape. The dominant phases are melilite, spinel, and Al-Ti pyroxene with minor anorthite. Spinel is subhedral and occurs in clusters. Some of these clusters have a circular geometry which encloses other minerals, known as a palisade structure [6] (Fig. 1). Some palisades form near perfect circles while others are more irregular in shape. Melilite in Saguaro ranges in size from coarse (>250m) to fine-grained (5-7m) and forms intergrown laths. A third of the melilite grains exhibit simple twinning about their <001> axis. The spinel palisades and twinned melilite in Saguaro suggest an igneous history, and the lack of secondary minerals suggests minimal aqueous alteration. At some point in time after the initial condensation of the minerals and formation of the inclusion, the sample was remelted and quickly solidified. The formation of the palisades is still heavily debated. One hypothesis is that the palisades are the rims of smaller CAIs that accreted early on, essentially acting as xenoliths within the larger CAIs [7]. Another hypothesis suggests an igneous origin for palisade structures [6-8] wherein the melt traps gas bubbles, and the spinel nucleates on the surface of this bubble. Based on the WDS spot analyses of 41 melilite grains using EPMA, the data suggest that the composition inside and outside the palisades is nearly identical. This finding indicates that these palisades are likely not exogenous but rather formed from melt-vapor reactions. Our results are consistent with studies by Simon and Grossman,1997 [6] and Zhang et al. (2019) [8]. Hoopoe. A ~0.5cm compact Type A CAI with an irregular shape. The dominant mineral phases are melilite, spinel, and hibonite with minor amounts of anorthite, augite, and perovskite. The melilite ranges in size from ~500 to 50µm. The larger melilite grains have simple twinning along the <001> axis like melilite in Saguaro. Melilite in Hoopoe exhibits crystal-plastic strain with misorientation dominantly about the <010> and <110> axes. Spinel shows subhedral to euhedral morphology and appears in clusters. Hibonite grains are similar to spinel in habit and size but show more plastic strain. The two minerals are often found together with one appearing to replace the other. Perovskite appears in fine grained recrystallized regions alongside fine augite and spinel. The abundant strain and deformation features in the melilite and hibonite suggest that Hoopoe experienced shock. This shock could have occurred in the nebula [9] or from an impact of another body on the parent body asteroid. The appearance of fine-grained (<10m) areas of augite, perovskite, and spinel in the dominantly coarse-grained inclusion suggest recrystallization, possibly due to the sudden increase in pressure and temperature. Cottonwood. This ~0.5cm CAI exhibits distinct mineralogy and textures suggesting a high degree of alteration. It is irregular in shape. The dominant mineral phases are spinel and anorthite with minor amounts of augite and rutile. The two main texture types can be seen in Fig. 2. The first type consists of coarse euhedral to subhedral spinel and anorthite. The second includes fine grained spinel, anorthite, rutile, and iron sulfides. Within these fine-grained regions, the anorthite grains are clustered into domains exhibiting the same crystallographic orientation. Rutile exclusively occurs with fine anorthite indicating a potential relationship between the two. Cottonwood has a clear and unbroken Wark-Lovering [10] rim (Fig. 2) on one side that consists of a sequence of spinel followed by anorthite and an outer layer of augite. The abundance of the fine-grained regions containing iron oxides and iron sulfides, secondary phases such as rutile, and oriented anorthite grains is evidence for recrystallization associated with a high degree of thermal metamorphism. The EPMA analyses on both coarse- and fine-grained spinel show that the fine spinel grains are more enriched in Cr (1-2 wt. %) than their coarse counterparts (0.1-0.5 wt. %). The data suggest thermal metamorphism drove chemical exchange of the previously refractory inclusion, introducing chromium and iron as well as sulfur, which is moderately volatile. Conclusions: The three CAIs analyzed record distinct secondary nebular and parent body processes. Saguaro melted in the nebula as seen by the twinned melilite and spinel palisades. Hoopoe experienced intense shock which deformed its melilite and recrystallized perovskite, augite, and spinel in fine-grained regions. Cottonwood shows evidence of recrystallization and chemical changes consistent with thermal metamorphism occurring after accretion into the parent-body.

V E Burnette↗

Geometry and dynamics of passive scalar structures in compressible turbulent mixing

We present a structure-based numerical analysis of passive scalar mixing in decaying homogeneous isotropic turbulence (DHIT) and shock-turbulence interaction canonical configurations. The analysis focuses on the temporal evolution of ensembles of passive scalar structures, initialized as spheres of different sizes relative to the Taylor microscale. An algorithm is introduced to track the evolution of each individual structure and the interactions with other structures in the ensemble, relating changes in the surface geometry and the underlying physical processes (turbulent transport, scalar dissipation, and shock compression). The tracking algorithm is applied to datasets from shock-capturing direct numerical simulations of DHIT, with Taylor microscale Reynolds number Re λ =40 and turbulence Mach number M t =0.2 ⁠, and STI cases in which the turbulence is processed by a shock wave at Mach numbers M = 1.5 and 3.0. Temporal surface convolution increases for initially larger structures, resulting in a higher probability of locally hyperbolic geometries where breakup into smaller structures occurs. Shock-induced deformation of the structures amplifies breakup processes, enhancing mixing, particularly for larger structures. Mixing enhancement by the shock is manifested as an amplification of the surface-averaged scalar gradient, which increases for initially larger structures. The alignment between the scalar gradient and the most extensional strain-rate eigendirection on the scalar isosurfaces also increases across the shock. Larger magnitudes of the scalar gradient and its alignment with the most compressive strain-rate eigendirection correlate with flatter surface regions. Shock-induced structure compression increases the area coverage of flat regions, where the amplification of scalar gradient is localized.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Inelastic deformation of diamond single crystals shock compressed to multimegabar stresses: Wave profile calculations

As the archetypal strong solid, the response of diamond shock compressed to multimegabar stresses is important for fundamental science and for numerical simulations of wave profiles for applications in high energy density physics experiments. Previous experiments and analysis have shown that the commonly used hydrodynamic assumption is invalid for diamond shock compressed to stresses below melt and an elastic–inelastic description is needed. Here, we present a phenomenological material model for calculating wave profiles in shock compressed diamond single crystals that incorporates this description. Also, to support the modeling effort, we carried out wave profile measurements on shock compressed diamond single crystals at the Sandia Z facility to augment previous measurements. Wave profiles for [100] and [111] diamond calculated using the material model provide a good match to the elastic–inelastic response (observed two-wave structure) measured at ∼325 and ∼360 GPa. Furthermore, the calculated peak stresses for single (overdriven) waves provide a good match to the measured Hugoniot states for stresses reaching ∼700 GPa, which is near melting conditions. The present results show that the diamond single crystal response at multimegabar shock stresses is characteristic of a brittle solid—pressure-dependent strength and strength loss due to inelastic deformation.

Deformation↗

Shock effects in meteorites

The impacts that can occur between objects on intersecting solar system orbits can generate shock-induced deformations and transformations, creating new mineral phases or melting old ones. These shock-metamorphic effects affect not only the petrography but the chemical and isotopic properties and the ages of primordial meteoritic materials. A fuller understanding of shock metamorphism and breccia formation in meteorites will be essential not only in the study of early accretion, differentiation, and regolith-evolution processes, but in the characterization of the primordial composition of the accreted material itself.

Stoeffler, D.↗

Virtual diffraction simulations using the quasi-coarse-grained dynamics method to understand and interpret plasticity contributions during in situ shock experiments

The dynamic deformation response of metallic materials has contributions from dislocations, deformation twinning, and plastic deformation. The current state-of-art techniques can detail the complex mechanistic history of deformation modes under shock loading in real-time using in situ X-ray diffraction (XRD). However, the capability of these experiments to unravel plasticity contributions is challenging due to limitations in interpreting results and the lack of validation from atomistic simulations. Molecular dynamics (MD) simulations can successfully capture various deformation modes in metals and complement experiments using simulated diffractograms at various stages of evolution. However, the difference in length and time scales of MD simulations and experiments is a substantial obstacle in completing and interpreting in situ diffractograms. Therefore, the existing modeling methods require various approximations to model defect evolution and interaction at the mesoscales. However, while using approximations to correlate the peak broadening behavior to the density of dislocations or the shifts/splitting due to the presence of twins, the interpretations of the plasticity contributions from diffractograms are non-trivial, especially when multiple modes of deformation may be operating. This viewpoint discusses combining a mesoscale modeling method called quasi-coarse-grained dynamics and virtual XRD to characterize the plasticity contributions in BCC metals from slip, twinning, and phase transformation behavior. Furthermore, the combined approach shows promise in bridging the mesoscale gap between the capabilities of atomic-scale simulations and in situ experiments to characterize the dynamic deformation of materials.

36 MATERIALS SCIENCE↗

Structure of Shock Waves and Inelasticity in Shock-Compressed Cemented Tungsten Carbides

In this study, shock wave experiments are conducted on General Carbide cemented tungsten carbide with 3.7wt.% cobalt binder to determine its shock-induced compression behavior up to 100 GPa. The measured wave profiles indicate the cemented tungsten carbide to undergo elastic-plastic deformation during shock compression. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot Elastic Limit (HEL), an elastic-plastic ramp indicating substantial post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, the shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relation, and the longitudinal stress (σ x ) vs. specific volume (V) curve for the samples. The HEL of the material was determined to lie between 4.41 GPa and 4.58 GPa. The U s - u p relation was determined to be U s = 4.97 + 1. 457u p for particle velocities greater than 0.75 km/s. The measured plastic shock velocities for particle velocities less than 0.7 km/s were found to be larger than those predicted using the linear U s - u p Hugoniot relationship, indicating the cemented WC samples to preserve substantial shear strength in the post-yield deformation region. No phase transformation was observed up to 100 GPa.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The History of Allan Hills 84001 Revised: Multiple Shock Events

The geologic history of Martian meteorite Allan Hills (ALH) 84001 is more complex than previously recognized, with evidence for four or five crater-forming impacts onto Mars. This history of repeated deformation and shock metamorphism appears to weaken some arguments that have been offered for and against the hypothesis of ancient Martian life in ALH 84001. Allan Hills 84001 formed originally from basaltic magma. Its first impact event (I1) is inferred from the deformation (D1) that produced the granular-textured bands ("crush zones") that transect the original igneous fabric. Deformation D1 is characterized by intense shear and may represent excavation or rebound flow of rock beneath a large impact crater. An intense thermal metamorphism followed D1 and may be related to it. The next impact (I2) produced fractures, (Fr2) in which carbonate "pancakes" were deposited and produced feldspathic glass from some of the igneous feldspars and silica. After I2, carbonate pancakes and globules were deposited in Fr2 fractures and replaced feldspathic glass and possibly crystalline silicates. Next, feldspars, feldspathic glass, and possibly some carbonates were mobilized and melted in the third impact (I3). Microfaulting, intense fracturing, and shear are also associated with 13. In the fourth impact (I4), the rock was fractured and deformed without significant heating, which permitted remnant magnetization directions to vary across fracture surfaces. Finally, ALH 84001 was ejected from Mars in event I5, which could be identical to I4. This history of multiple impacts is consistent with the photogeology of the Martian highlands and may help resolve some apparent contradictions among recent results on ALH 84001. For example, the submicron rounded magnetite grains in the carbonate globules could be contemporaneous with carbonate deposition, whereas the elongate magnetite grains, epitaxial on carbonates, could be ascribed to vapor-phase deposition during I3.

Treiman, Allan H.↗

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↗

Dynamic Strength of Iron at High Pressures and Strain Rates

Accurate modeling of meteorite impacts, and deformation of planetary cores require characterization of the flow strength and in-elasticity of iron in its different phases. In this letter, we investigate the flow strength of both the ambient α-phase and high-pressure ϵ-phase of iron at strain-rates of 1 × 10 5 s –1 and pressures up to 42 GPa using high pressure - pressure shear plate impact experiments. We report the strength of the ϵ-iron to be significantly higher than α-phase but consequently one order smaller than the previously reported dynamic strength at high pressures. As a result, the complete stress-strain response of the ϵ-phase is reported for the first time.

36 MATERIALS SCIENCE↗

Vanishing dynamic strength measured in a transient high-pressure phase of tin

This work reports two independent experimental estimations of the dynamic strength of the body-centered tetragonal (bct) 𝛾 phase of tin, which occurs at pressures above about 9 GPa and transforms promptly back to the ambient 𝛽 phase upon pressure release. Measuring strength in such a transient high-pressure phase is challenging. One measurement used free-surface Richtmyer-Meshkov instabilities generated with gas gun impact. Another set of measurements used ramp-release loading in a pulsed-power facility. Strength estimations came from comparison to simulations using a comprehensive multiphase modeling framework that includes equation of state, shear moduli, and strength separately for each phase, and treats mixed-phase regions. Both experimental methods found the 𝛾-phase strength to be very low, within experimental uncertainty of zero. The existing literature on other materials, by contrast, almost universally reports higher strength in transient high-pressure phases compared with ambient phases. Recent Molecular-Dynamics simulations on tin in the literature showed almost zero deviatoric strength during a deformation-induced bct → bct transformation in which one of the 𝛾-phase 𝑎 axes “flips” to the 𝑐 axis. This reorientation currently provides the most plausible explanation for the low observed strength in 𝛾-phase tin.

36 MATERIALS SCIENCE↗