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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Shock induced magnetic effects in fine particle iron dispersions

Magnetic effects associated with shock induced transformation of fcc antiferromagnetic iron precipitates in polycrystalline copper disks at levels up to 5 GPa in weak magnetic fields (H not greater than 0.5 Oe) were investigated. The demagnetization and anisotropy associated with second order transition, the effects of plastic deformation in imparting magnetic anisotropy and magnetic hardening, and the influence of post shock thermal transients on magnetization associated with recovery, recrystallization and grain growth were studied. It was found that on the microsecond time scale of the shock induced first order transformation, the field sense is recorded in the transformed iron particles. For a given particle size the degree of transformation of fcc iron depends on the level of the shock. For a given shock level the resultant magnetic properties depend on the particle size distribution, with maximum effects noted in specimens with 400 to 600 A particles.

Wasilewski, P. J.↗
DOE OSTI2024

Bayesian blacksmithing: discovering thermomechanical properties and deformation mechanisms in high-entropy refractory alloys

Finding alloys with specific design properties is challenging due to the large number of possible compositions and the complex interactions between elements. This study introduces a multi-objective Bayesian optimization approach guiding molecular dynamics simulations for discovering high-performance refractory alloys with both targeted intrinsic static thermomechanical properties and also deformation mechanisms occurring during dynamic loading. The objective functions are aiming for excellent thermomechanical stability via a high bulk modulus, a low thermal expansion, a high heat capacity, and for a resilient deformation mechanism maximizing the retention of the BCC phase after shock loading. Contrasting two optimization procedures, we show that the Pareto-optimal solutions are confined to a small performance space when the property objectives display a cooperative relationship. Conversely, the Pareto front is much broader in the performance space when these properties have antagonistic relationships. Density functional theory simulations validate these findings and unveil underlying atomic-bond changes driving property improvements.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Teleseismic detection of a slow precursor to the great 1989 Macquarie Ridge earthquake

Low-frequency spectra for the 1989 Macquarie Ridge earthquake (magnitude 8.2) show an amplitude increase and a phase-delay decrease below 6 millihertz that require a short-term slow precursor. This earthquake can be modeled as a compound event in which a fast-rupturing, ordinary earthquake was initiated by an episode of slow, smooth deformation that began more than 100 seconds before the main shock. The moment released in the slow precursor was large, about 3 x 10 exp 20 newton-meters, equivalent to an event of magnitude 7.6. The data are consistent with the precursor being generated in a region of the oceanic upper mantle below the main rupture.

Ihmle, Pierre F.↗
NASA NTRS1995

Control of unsteady separated flow associated with the dynamic stall of airfoils

An effort to understand and control the unsteady separated flow associated with the dynamic stall of airfoils was funded for three years through the NASA cooperative agreement program. As part of this effort a substantial data base was compiled detailing the effects various parameters have on the development of the dynamic stall flow field. Parameters studied include Mach number, pitch rate, and pitch history, as well as Reynolds number (through two different model chord lengths) and the condition of the boundary layer at the leading edge of the airfoil (through application of surface roughness). It was found for free stream Mach numbers as low as 0.4 that a region of supersonic flow forms on the leading edge of the suction surface of the airfoil at moderate angles of attack. The shocks which form in this supersonic region induce boundary-layer separation and advance the dynamic stall process. Under such conditions a supercritical airfoil profile is called for to produce a flow field having a weaker leading-edge pressure gradient and no leading-edge shocks. An airfoil having an adaptive-geometry, or dynamically deformable leading edge (DDLE), is under development as a unique active flow-control device. The DDLE, formed of carbon-fiber composite and fiberglass, can be flexed between a NACA 0012 profile and a supercritical profile in a controllable fashion while the airfoil is executing an angle-of-attack pitch-up maneuver. The dynamic stall data were recorded using point diffraction interferometry (PDI), a noninvasive measurement technique. A new high-speed cinematography system was developed for recording interferometric images. The system is capable of phase-locking with the pitching airfoil motion for real-time documentation of the development of the dynamic stall flow field. Computer-aided image analysis algorithms were developed for fast and accurate reduction of the images, improving interpretation of the results.

Wilder, M. C.↗
NASA NTRS1992

Simulated meteorite impacts and volcanic explosions: Ejecta analyses and planetary implications

Past cratering studies have focused primarily on crater morphology. However, important questions remain about the nature of crater deposits. Phenomena that need to be studied include the distribution of shock effects in crater deposits and crater walls; the origin of mono- and polymict breccia; differences between local and distal ejecta; deformation induced by explosive volcanism; and the production of unshocked, high-speed ejecta that could form the lunar and martian meteorites found on the Earth. To study these phenomena, one must characterize ejecta and crater wall materials from impacts produced under controlled conditions. New efforts at LLNL simulate impacts and volcanism and study resultant deformation. All experiments use the two-stage light-gas gun facility at LLNL to accelerate projectiles to velocities of 0.2 to 4.3 km/s, including shock pressures of 0.9 to 50 GPa. We use granite targets and novel experimental geometries to unravel cratering processes in crystalline rocks. We have thus far conducted three types of simulations: soft recovery of ejecta, 'frozen crater' experiments, and an 'artificial volcano. Our ejecta recovery experiments produced a useful separation of impactites. Material originally below the projectile remained trapped there, embedded in the soft metal of the flyer plate. In contrast, material directly adjacent to the projectile was jetted away from the impact, producing an ejecta cone that was trapped in the foam recovery fixture. We find that a significant component of crater ejecta shows no signs of strong shock; this material comes from the near-surface 'interference zone' surrounding the impact site. This phenomenon explains the existence of unshocked meteorites on the Earth of lunar and martian origin. Impact of a large bolide on neighboring planets will produce high-speed, weakly shocked ejecta, which may be trapped by the Earth's gravitational field. 'Frozen crater' experiments show that the interference zone is highly localized; indeed, disaggregation does not extend beyond approx. 1.5 crater radii. A cone-shaped region extending downward from the impact site is completely disaggregated, including powdered rock that escaped into the projectile tube. Petrographic analysis of crater ejecta and wall material will be presented. Finally, study of ejecta from 0.9- and 1.3-GPa simulations of volcanic explosions reveal a complete lack of shock metamorphism. The ejecta shows no evidence of PDF's, amorphization, high-pressure phases, or mosaicism. Instead, all deformation was brittle, with fractures irregular (not planar) and most intergranular. The extent of fracturing was remarkable, with the entire sample reduced to fragments of gravel size and smaller.

Gratz, A. J.↗
NASA NTRS2014

Impact-Induced Chondrule Deformation and Aqueous Alteration of CM2 Murchison

Deformed chondrules in CM2 Murchison have been found to define a prominent foliation [1,2] and lineation [3] in 3D using X-ray computed tomography (XCT). It has been hypothesized that chondrules in foliated chondrites deform by "squeezing" into surrounding pore space [4,5], a process that also likely removes primary porosity [6]. However, shock stage classification based on olivine extinction in Murchison is consistently low (S1-S2) [4-5,7] implying that significant intracrystalline plastic deformation of olivine has not occurred. One objective of our study is therefore to determine the microstructural mechanisms and phases that are accommodating the impact stress and resulting in relative displacements within the chondrules. Another question regarding impact deformation in Murchison is whether it facilitated aqueous alteration as has been proposed for the CMs which generally show a positive correlation between degree of alteration and petrofabric strength [7,2]. As pointed out by [2], CM Murchison represents a unique counterpoint to this correlation: it has a strong petrofabric but a relatively low degree of aqueous alteration. However, Murchison may not represent an inconsistency to the proposed causal relationship between impact and alteration, if it can be established that the incipient aqueous alteration post-dated chondrule deformation. Methods: Two thin sections from Murchison sample USNM 5487 were cut approximately perpendicular to the foliation and parallel to lineation determined by XCT [1,3] and one section was additionally polished for EBSD. Using a combination of optical petrography, SEM, EDS, and EBSD several chondrules were characterized in detail to: determine phases, find microstructures indicative of strain, document the geometric relationships between grain-scale microstructures and the foliation and lineation direction, and look for textural relationships of alteration minerals (tochilinite and Mg-Fe serpentine) that indicate timing of their formation relative to deformation event(s). Preliminary Results: Deformed chondrules are dominated by forsterite and clinoenstatite with lesser amounts of Fe-Mg serpentine, sulfides, and low calcium pyroxene. Olivine grains are commonly fractured but generally show sharp optical extinction. The pyroxene, in contrast, is not only fractured but also often displays undulose extinction. In addition, the clinoenstatite is frequently twinned but it is unclear whether the twins are the result of mechanical deformation or inversion from protoenstatite [8]. EBSD work is currently ongoing to determine if areas of higher crystallographic strain can be imaged and mapped, and to determine the pyroxene twin orientations. In regards to alteration, we have found evidence for post-deformation formation of tochilinite and Mg-Fe serpentine indicating that aqueous alteration has indeed post-dated the deformation of the chondrules.

Hanna, R. D.↗
DOE OSTI2023

Shock response of cobalt-based cemented tungsten carbides at pressures up to 100 GPa

In the present study, shock wave experiments are conducted on General Carbide cemented tungsten carbides with 3.7 wt.% and 6.0 wt.% cobalt binder to determine their shock compression response up to 100 GPa. 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 post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relationship, and the longitudinal stress vs. specific volume Hugoniot relationship for the two samples. Here, the HEL for the WC 3.7wt% Co and WC 6.0wt% Co samples was determined to be ~4.45 GPa and 3.72 GPa, respectively. The U s – u p relation was determined to be U s = 4.97(0.006)+1. 446 (0.018) u p for WC 3.7 wt.% Co and U s = 4.93(0.006)+1. 454(0.017)u p for the WC 6.0 wt.% Co sample at peak particle velocities>0.75 km/s. For both WC grades, in the particle velocity regime less than ~0.75 km/s, the measured shock wave velocities were found to be larger than those predicted by the linear U s –u p Hugoniot relationship, indicating the two WC samples to preserve substantial strength in the post-yield deformation regime. Both WC grades show a catastrophic drop in shear stress carrying capacity when shocked to longitudinal stresses greater than ~70 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗
NASA NTRS1981

Unsteady environments and responses of the shuttle combined loads orbiter test

Both separate and combined wind tunnel and vibration shaker tests were conducted on two structural panels representative of the Shuttle orbiter in the NASA LaRC 8-foot transonic pressure tunnel (TPT) to determine the effects of combined loads on the thermal protection system (TPS). The primary objective of this test was to provide a combined full-scale load environment and realistic time history of the dynamic pressures, Mach numbers (through transonic), and dynamic structural responses of these panels. The panels were selected from orbiter locations where interactive load sources such as aerodynamic shock waves, turbulent boundary layers, strut-induced vorticity, and substrate deformation combined to provide high bonding loads between the TPS and the orbiter structure. The test panels for this program were highly instrumented with static and dynamic pressure gages, accelerometers, deflectometers, strain gages, Schileren and high-speed photography, and special instrumentation necessary to determine TPS/structure interface loads and tile motions. Two test specimens of each orbiter panel were utilized. Both were high-fidelity representations of the selected orbiter location. The first panel (calibration panel) for each test was utilized to calibrate the tunnel and shaker system as well as provide a data-gathering source for applied static and dynamic loads. The second panel (test panel) for each test was subjected to flight time exposures only so that TPS characteristics as a function of flight mission exposure could be determined. These tests were the last in a series that was suce6ssfully completed and necessary to certify the TPS prior to the first launching of the Shuttle.

P. H. Schuetz↗
DOE OSTI2026

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

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

42 ENGINEERING↗
DOE OSTI2022

Effect of Topology on Transient Dynamic and Shock Response of Polymeric Lattice Structures

Architected cellular materials, such as lattice structures, offer potential for tunable mechanical properties for dynamic applications of energy absorption and impact mitigation. In this work, the static and dynamic behavior of polymeric lattice structures was investigated through experiments on octet-truss, Kelvin, and cubic topologies with relative densities around 8%. Here, dynamic testing was conducted via direct impact experiments (25–70 m/s) with high-speed imaging coupled with digital image correlation and a polycarbonate Hopkinson pressure bar. Mechanical properties such as elastic wave speed, deformation modes, failure properties, particle velocities, and stress histories were extracted from experimental results. At low impact velocities, a transient dynamic response was observed which was composed of a compaction front initiating at the impact surface and additional deformation bands whose characteristics matched low strain-rate behavior. For higher impact velocities, shock analysis was carried out using compaction wave velocity and Eulerian Rankine–Hugoniot jump conditions with parameters determined from full-field measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2024

New dynamic diamond anvil cell for time-resolved radial x-ray diffraction

The dynamic diamond anvil cell (dDAC) is a recently developed experimental platform that has shown promise for studying the behavior of materials at strain rates ranging from intermediate to quasi-static and shock compression regimes. Combining dDAC with time-resolved x-ray diffraction (XRD) in the radial geometry (i.e., with incident x-rays perpendicular to the axis of compression) enables the study of material properties such as strength, texture evolution, and deformation mechanisms. This work describes a radial XRD dDAC setup at beamline P02.2 (Extreme Conditions Beamline) at DESY’s PETRA III synchrotron. Time-resolved radial XRD data are collected for titanium, zirconium, and zircon samples, demonstrating the ability to study the strength and texture of materials at compression rates above 300 GPa/s. In addition, the simultaneous optical imaging of the DAC sample chamber is demonstrated. The ability to conduct simultaneous radial XRD and optical imaging provides the opportunity to characterize plastic strain and deviatoric strain rates in the DAC at intermediate rates, exploring the strength and deformation mechanisms of materials in this regime.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗
NASA NTRS1974

Evaluate ERTS imagery for mapping and detection of changes of snowcover on land and on glaciers

The author has identified the following significant results. The area of snowcover on 10 individual drainage basins in the North Cascades, Washington, has been determined by use of a semi-automatic radiance threshold technique. The result is a unique record of the changing water storage as snow in these important hydrologic units, the runoff of which is utilized for hydroelectric power, dilution of wastes and heat, support of salmon migration, and irrigation. These data allow a new type of hydrologic modelling to proceed which should permit more accurate forecasts of streamflow. A new technique has been developed for measuring snow-covered area or snowline altitude semi-automatically. This variable contour overlay method permits the snowcover to be matched efficiently to the best fit contour of altitude. The motion of the Yentna Glacier during the concluding phase of its surge was successfully measured by a flicker technique using images of two dates. It appears that displacements as small as 100 to 200 m can be measured. Motion of the Tweedsmuir Glacier in Alaska was measured using ERTS-1 images enlarged to 1:50,000. Changes detected included a shock wave moving down the glacier, the margin expanding, the moraine pattern deforming, and the marginal valley deepening.

Meier, M. F.↗
NASA NTRS1993

The impact of supernova fragments on the evolution of multisupernova remnants

Analytical approximations and 2D hydrodynamical simulations are used to examine the interaction of supernova fragments with the internal structure of large multisupernova remnants (MSRs). The fragments are thermalized by reverse shocks generated in the interaction with the MSR interior, which is assumed to be hot and rarefied. The evolution is divided into two stages: before and after reaching a reference distance, R(E), from the explosion site. As the density of the expanding fragment drops, the reverse shock accelerates, and, when the distance R(E) is reached, it begins to effectively erode the fragment. At some selected evolutionary times, the X-ray emission from the shocked fragment is also calculated. The direct bombardment of the MRS shell by the shocked fragment has a series of important consequences: it excites, punctures, and deforms the expanding shell.

Franco, J.↗
DOE OSTI2020

High-speed x-ray phase contrast imaging and digital image correlation analysis of microscale shock response of an additively manufactured energetic material simulant

The performance of energetic materials subjected to dynamic loading significantly depends on their micro- and meso-scale structural morphology. The geometric versatility offered by additive manufacturing opens new pathways to tailor the performance of these materials. Additively manufactured energetic materials (AMEMs) have a wide range of structural characteristics with a hierarchy of length scales and process-inherent heterogeneities, which are hitherto difficult to precisely control. It is important to understand how these features affect AMEMs’ response under dynamic/shock loading. Therefore, temporally and spatially resolved measurements of both macroscopic behavior and micro- and meso-level processes influencing macroscopic behavior are required. In this paper, we analyze the shock compression response of an AMEM simulant loaded under several impact conditions and orientations. Furthermore, x-ray phase contrast imaging (PCI) is used to track features across the observed shock front and determine the linear shock velocity vs particle velocity equation of state, as well as to quantify the interior deformation fields via digital image correlation (DIC) analyses. Photon Doppler velocimetry is simultaneously used to measure the particle velocities of the specimens, which are consistent with those obtained from x-ray PCI. The DIC analyses provide an assessment of the average strain fields inside the material, showing that the average axial strain depends on the loading intensity and reaches as high as 0.23 for impact velocities up to 1.5 km/s. The overall results demonstrate the utility of x-ray PCI for probing “in-material” equation of state and interior strains associated with dynamic shock compression behavior of the AMEM simulant.

3D printing↗