Search NASA⌕ Search

SEARCH · Search NASA

Results for “ELASTIC PROPERTY”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Anomalous elastic softening in ferroelectric hafnia under pressure

his study employs first-principles density-functional theory (DFT) calculations to explore the elastic and mechanical properties of ferroelectric hafnia (HfO 2 ) in its polar orthorhombic 𝑃⁢𝑐⁢𝑎⁢2 1 phase under varying hydrostatic pressure conditions up to 30 GPa. Utilizing a plane-wave basis set and Perdew-Burke-Ernzerhof generalized-gradient approximation for solids in our DFT calculations, we investigate both pure and yttrium-substituted HfO 2 . Our findings reveal an anomalous reduction in the 𝐶 33 component of the elastic tensor with increasing pressure, which becomes significant above 15 GPa and signals a potential pressure-driven structural phase transition at higher pressure. The analysis of atomic displacements under pressure sheds light on the unusual mechanical behavior and phase stability of this material. Additionally, we observe a transition from an indirect band gap to a direct band gap with increasing pressure, which could have significant implications for optical applications. Here, the effects of yttrium substitution on the mechanical and electronic properties are further examined, revealing that yttrium substitution softens the elastic response of this material and reduces the electronic band gap. These results enhance our understanding of elastic and mechanical responses of ferroelectric hafnia and its potential for applications in microelectronics, piezoelectric devices, and nonvolatile ferroelectric random-access memories. Further experimental validation is recommended to confirm our predictions and explore the practical implications of the observed phase transitions and electronic behavior of the ferroelectric hafnia under high-pressure conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Developing reliable machine learning interatomic potential for Fe–Cr–Ni austenitic alloys

Gaining atomistic understanding of mechanical behavior of heat-resistant structural materials such as Fe–Cr–Ni-based alloys requires an approach with an accuracy close to density functional theory (DFT) that considers the intrinsic properties of the bulk lattice and important defects such as stacking faults, grain boundaries, and surfaces. This work aims to develop reliable machine learning interatomic potential (MLIAP) at cross-scale for Fe–Cr–Ni ternary alloys with a focus on the face-centered-cubic (fcc) solid solution structure. Leveraging the advantages of moment tensor potentials, which typically necessitate a relatively small training dataset and enable rapid calculations using the large-scale atomic/molecular massively parallel simulator package, we ensure the stability and accuracy of the trained potentials. Important defects such as stacking faults, grain boundaries, and surfaces for wide-range compositions are investigated. Structural, thermal, elastic, and defect properties are determined from molecular dynamics simulations comprising several thousand atoms, generated via canonical Monte Carlo simulations guided by the trained potential. The trained potential allows efficient atomic simulations of structural, thermal, and mechanical properties of fcc Fe–Cr–Ni solid solution alloys as a function of composition and temperature. Therefore, the MLIAP approach represents a major advancement from DFT calculations that are limited to small simulation sizes and traditional molecular dynamics simulations using relatively low accuracy potentials. Furthermore, this work outlines a practical foundation for further investigating the structural evolution and mechanical behavior of austenitic stainless steel and nickel-based alloys in a wide array of applications in extreme environments.

Crystal structure↗
DOE OSTI2025

A Green’s function fast multipole method for computation of micromechanical fields in heterogeneous materials

Computation of micromechanical fields in heterogeneous materials is usually performed using either the finite element method or the Green’s function method based on FFTs. The finite element method allows for accurate discretization and for non-periodic boundary conditions but is computationally expensive. On the other hand, the FFT-based method is computationally efficient but requires discretization on a regular grid of hexahedral voxels. In this paper, a Green’s function method allowing for accurate discretization using tetrahedral elements and for non-periodic boundary conditions is proposed. The convolution is computed using the fast multipole method, which provides good accuracy even for low-order expansion due to the fast decay of interactions between elements. The proposed Green’s function fast multipole method is verified by comparison with analytical and FFT-based solutions. Furthermore, the computational time is analyzed and compared to the FFT-based method for non-periodic convolution. Finally, effective properties of an elastic polycrystalline microstructure containing thin intergranular cracks are computed and analyzed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Non-invasive Pipe Pressure Monitoring for Safeguards

Our project explores a new solution to complement current safeguards toolbox by independently, and noninvasively, determining the pressure inside uranium hexafluoride (UF 6 )-carrying Gas Centrifuge Enrichment Plant (GCEP) pipes. Proof-of-principle, non-invasive pipe pressure monitoring tests were conducted through measurements of hoop and axial strain using an optical-interference dilatometry technique based on fiber Bragg gratings (FBGs) inscribed in the core of telecommunication-type 125 µm diameter optical fibers. The fibers were affixed to commercial stainless steel and aluminum tubes and pipes of the type expected to be found at GCEP’s, each measuring four inches in diameter. When the internal air pressure in the pipes changes their external diameter and length follow, and the pitch in the FBG shrinks and/or expands accordingly. In turn, these changes result in a shift of the light wavelength reflected by the FBGs. Several FBGs physically attached to the pipe, as well as control FBGs not attached, were interrogated during the process using state-of-the-art commercial interrogators that operate on swept-wavelength laser technology. Strains measured were contrasted against calculations based on the known elastic and mechanical properties of the pipes. Good agreement is found, and strategies are proposed for further improvement of the experimental resolution as well as future implementation of this novel approach for GCEP safeguards.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling

Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.

Dynamic simulation↗
DOE OSTI2026

Elastic strain engineering of lattice thermal conductivity of silicon: An ab-initio study

Silicon (Si) is the most essential material in the semiconductor industry. It is important to manage the thermal properties of crystalline Si. Elastic strain engineering (ESE) has proven to be an effective tool in controlling the electrical conductivity of Si in strained-silicon technology; its effects on the thermal conductivity of silicon, therefore, warrants careful investigation. The ESE effect is much more pronounced for nanostructured materials due to the ultralarge elastic strains (on the order of 10%) achievable at the nanoscale. In this work, the lattice thermal conductivity (κ L ) of Si under hydrostatic, biaxial, and uniaxial strain states is studied with ab-initio simulations, and the values of strain-dependent κ L compare well with experimental results and existing molecular dynamics simulations. To understand the mechanisms of strain-modulated κ L , the phonon bands, scattering rate, and Grüneisen parameters of phonon modes are computed. It is shown that strain can significantly change the anharmonicity of the crystal system, thus changing phonon scattering rates and κ L . Our results demonstrate that ESE can reduce silicon κ L by up to approximately 90%. Furthermore, uniaxial and biaxial strains can induce highly anisotropic thermal conductivity in Si, with relative variations up to 58.5% and 14.5%, respectively.

Anisotropy in Thermal Conductivity↗
DOE OSTI2025

Development of a thermal creep model for aluminum alloy 6061 cladding in U-10Mo monolithic fuel plates

Plate-type fuel elements consisting of a high-density, low-enriched uranium (LEU) U–10Mo-based fuel foil encapsulated in an aluminum alloy (AA) cladding are fabricated using the hot isostatic pressing (HIP) technique. During the HIP process, the fuel plate system is heated to 560 °C, then cooled to room temperature. This heat cycle significantly affects the mechanical properties of the aluminum cladding, and experimental investigations have shown that, post-HIP bonding, the mechanical properties of the aluminum cladding transition from those of AA 6061-T6 to something closer to the O temper. More specifically, the ultimate strength of the cladding decreases while its ductility increases, making it challenging to capture the changes in mechanical behavior and material properties. Understanding the residual stresses generated during the HIP process is critical for assessing the fuel plate’s integrity under various temperature, pressure, and irradiation. To simulate the HIP bonding process, the elastic, plastic, and thermal properties of the cladding are assumed to be similar to those of AA 6061-O temper. However, the primary challenge lies in the lack of available data for the creep model of the AA 6061 cladding during this transient process of HIP. The present study focuses on developing a computational model that predicts the creep behavior of the aluminum cladding in the fuel plates during the HIP process, as cladding creep significantly influences the residual stresses generated in U-10Mo fuel plates during HIP fabrication. Furthermore, as HIP bonding occurs at high temperatures that are nearing the melting point of aluminum, the present work considered a temperature-dependent Arrhenius-type creep model. In particular, a hyperbolic sine creep model is employed to estimate the creep properties of the as-fabricated aluminum cladding. In conclusion, the residual stresses predicted in the U-10Mo fuel when using the newly calibrated creep model closely align with the experimental measurements, validating the model’s accuracy.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Reversible Polymer–Metal Mechanical Transitions Enabled by Electrochemical Modulation

The mechanical properties of materials, e.g., elastic modulus and hardness, are important for many engineering applications. Here, we introduce a hierarchically ordered structure (HOS) polymer system that exhibits exceptionally large, reversible changes in mechanical behavior upon electrochemical lithiation and delithiation. Full lithiation of the HOS polymer induces a substantial mechanical transition from polymer-like to metal-like attributes, yielding a 10-fold increase in elastic modulus and a 3-fold increase in hardness, with values comparable to those of aluminum. Upon removal of Li+ ions, the elastic modulus and hardness return to nearly pristine polymer levels, and this transformation remains highly repeatable over many electrochemical cycles. Reversible transitions between polymer-like and metal-like mechanical behaviors offer a new pathway for engineering materials for applications that require tunable mechanical properties, such as soft robotics and stimuli-responsive systems.

Thapa, Santosh↗
DOE OSTI2026

Effects of PLA/PHB Blend Ratio and Wood Flour Loading on the Melt Rheology and Thermomechanical Properties of Biobased Polymer Composites

The development of sustainable, biobased polymer composites is crucial for reducing supply chain dependence on fossil fuels. A major challenge lies in understanding and predicting the processability of these sustainable material alternatives, which directly impacts large-scale manufacturing. Here, this study systematically investigates the effect of poly(lactic acid) (PLA)/polyhydroxybutyrate (PHB) blend ratios and wood flour (WF) loadings on the thermal, rheological, and mechanical properties of the composites. Composites were prepared via melt compounding and characterized by using differential scanning calorimetry (DSC), thermogravimetric analysis (TGA), oscillatory shear melt rheology, tensile testing, and scanning electron microscopy. Results showed that increasing the PHB content lowered the glass transition temperature, while higher PHB and WF loadings decreased thermal stability. Complex viscosity decreased with increasing PHB in 0 and 10% WF, whereas it increased in 20% and 30% WF for PHB-rich blends (80/20 and 70/30 PLA/PHB) due to stronger filler interactions with flexible PHB chains. At 30% WF, PHB-containing composites also showed elastic dominance with G′ higher than G″. This behavior was confirmed by van Gurp–Palmen plots, where the phase angle decreased with increasing PHB and WF loadings. WF addition also shifted Cole–Cole plots away from semicircular terminal relaxation, while time–temperature superposition validity was maintained, confirming that polymer chain dynamics continue to dominate melt behavior. Young’s modulus increased with WF loading, while tensile strength and toughness decreased due to weaker interfacial adhesion between the matrix and the filler. Interestingly, 10% WF systems containing PHB exhibited the highest toughness among all of the filled systems, indicating synergistic reinforcement from WF fillers and PHB-induced ductility. Overall, melt rheology effectively captured the internal structure and stiffness of the composites as the formulation changed, reinforcing the concept that the PLA/PHB ratio and WF content can be tuned to adjust melt elasticity and balance mechanical property trade-offs for targeted applications.

melt rheology↗
DOE OSTI2024

Optimizing Deep Geothermal Drilling for Energy Sustainability in the Appalachian Basin

This study investigates the geological and geomechanical characteristics of the MIP 1S geothermal well in the Appalachian Basin to optimize drilling and address the wellbore stability issues encountered. Data from well logs, sidewall core analysis, and injection tests were used to derive elastic and rock strength properties, as well as stress and pore pressure profiles. A robust 1D-geomechanical model was developed and validated, correlating strongly with wellbore instability observations. This revealed significant wellbore breakout, widening the diameter from 12 ¼ inches to over 16 inches. Advanced technologies like Cerebro Force™ In-Bit Sensing were used to monitor drilling performance with high accuracy. This technology tracks critical metrics such as bit acceleration, vibration in the x, y, and z directions, Gyro RPM, stick-slip indicators, and bending on the bit. Cerebro Force™ readings identified hole drag caused by poor hole conditions, including friction between the drill string and wellbore walls and the presence of cuttings or debris. This led to higher torque and weight on bit (WOB) readings at the surface compared to downhole measurements, affecting drilling efficiency and wellbore stability. Optimal drilling parameters for future deep geothermal wells were determined based on these findings.

Environmental Sciences & Ecology↗
DOE OSTI2024

Deep-Learning Interatomic Potential Connects Molecular Structural Ordering to the Macroscale Properties of Polyacrylonitrile

Polyacrylonitrile (PAN) is an important commercial polymer, bearing atactic stereochemistry resulting from nonselective radical polymerization. As such, an accurate, fundamental understanding of governing interactions among PAN molecular units is indispensable for advancing the design principles of final products at reduced processability costs. While ab initio molecular dynamics (AIMD) simulations can provide the necessary accuracy for treating key interactions in polar polymers, such as dipole–dipole interactions and hydrogen bonding, and analyzing their influence on the molecular orientation, their implementation is limited to small molecules only. Herein, we show that the neural network interatomic potentials (NNIPs) that are trained on the small-scale AIMD data (acquired for oligomers) can be efficiently employed to examine the structures and properties at large scales (polymers). NNIP provides critical insight into intra- and interchain hydrogen-bonding and dipolar correlations and accurately predicts the amorphous bulk PAN structure validated by modeling the experimental X-ray structure factor. Furthermore, the NNIP-predicted PAN properties, such as density and elastic modulus, are in good agreement with their experimental values. Overall, the trend in the elastic modulus is found to correlate strongly with the PAN structural orientations encoded in the Hermans orientation factor. In conclusion, this study enables the ability to predict the structure–property relations for PAN and analogues with sustainable ab initio accuracy across scales.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effects of HCP/BCC element ratios on the room-temperature tensile properties of Ti-Zr-Hf-Nb-Ta refractory high-entropy alloys

Equiatomic and non-equiatomic Ti-Zr-Hf-Nb-Ta refractory high-entropy alloys (RHEAs) were arc melted, homogenized, cold rolled, and recrystallized to produce single-phase, body-centered cubic (BCC), microstructures with weak texture and equiaxed grains 76–199 μm in size. Here, the non-equiatomic alloys had either a 60:40 or 80:20 atomic ratio of hexagonal close-packed (HCP) elements (Ti + Zr + Hf) to BCC elements (Nb + Ta). Alloy compositions were measured after thermomechanical processing to determine the concentrations of the major (substitutional) and minor (interstitial) elements. We investigated how elastic constants and uniaxial tensile properties were affected by changes in the relative concentrations of the constituent elements at fixed HCP:BCC ratios. Yield strengths ranged from 801 to 922 MPa and ultimate tensile strengths from 815 to 933 MPa. Good agreement is obtained between the experimental yield strengths and those predicted by a strength theory based on edge dislocations indicating that the observed compositional effects are due to their effects on shear modulus and volume misfit. Fracture occurred by dimpled rupture with fracture strains of 19.4%–25.7%, but uniform strains were an order of magnitude lower at 1.1%–3.2%, calling into question the useable ductility (prior to necking) of RHEAs considered to be ductile based on their fracture strain. Contrary to predictions in the literature that HCP elements promote ductility, our present results show that increasing the HCP:BCC ratio decreases both the total strain and the uniform strain. Similar trends were not evident in the yield or ultimate strengths; rather, strengths were affected mainly by shear modulus and volume misfit.

BCC high-entropy alloys↗
DOE OSTI2024

Ductile Glassy Polymer Networks Capable of Large Plastic Deformation and Heat-Induced Elastic Recovery

Many thermoplastic polymers are ductile by combining strength and large deformations. These deformations are irreversible - known as plastic deformation. Elastomers can deform reversibly but have low strength. To this end, we developed glassy and ductile polyamide networks capable of large plastic deformation (>200% strain) and high strength (~50 MPa tensile strength and ~1500 MPa Young’s modulus), similar to those of polyolefins and Nylon-66. We discovered that hydrogen bonding between meta-phthalamide groups was essential to the ductility. Since these polyamide networks are covalently bonded, we demonstrated their unique durability by repeatable elastic recovery at elevated temperatures, exhibiting indifferent tensile properties in each cycle. Furthermore, when we fixed the strain during the elastic recovery, these polyamide networks actuated stresses of 9–18 MPa, among the highest reported in shape-memory polymer actuators. Here, we envision these ductile, glassy polymer networks as promising alternatives to ductile thermoplastics, given the combined benefit of ductility and durability.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Additively manufacturing high-temperature aluminum alloy for automotive engine pistons

The purpose of the project was to evaluate DuAlumin-3D, an ORNL designed Al alloy, for additive manufacturing of high-performance automotive components. The project consisted of additively manufacturing a cylindrical block of representative material which was then heat treated before machining into test coupons. Testing included tensile tests as a function of temperature, high-cycle fatigue testing as a function of temperature, density and elastic modulus measurements, and thermophysical property measurements. Fractography was also performed on failed high-cycle fatigue specimens. The results were compared against benchmark AM Al alloys and current high-performance conventionally manufactured alloys. Overall, the DuAlumin-3D material outperformed the incumbent materials in mechanical properties, especially at elevated temperatures, while maintaining broadly similar thermophysical properties. The results of this study support further maturation of the material for additive manufacturing of high-performance automotive components.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Prospects for exploring non-standard neutrino properties with argon-based CEvNS experiments

Coherent elastic neutrino-nucleus scattering (CEvNS) provides a powerful frame-work for testing the Standard Model (SM) and searching for new physics at low energies. In this work, we examine the prospects for argon-based CEvNS experiments at stopped-pion sources to perform precision measurements of weak interactions and probe non-standard neutrino properties. Our study focused on the CENNS-10 and CENNS-750 detectors at the Spallation Neutron Source at Oak Ridge National Laboratory, the Coherent Captain Mills (CCM) detector at Los Alamos National Laboratory, and the proposed PIP2-BD detector at Fermilab’s Facility for Dark Matter Discovery (F2D2). Using realistic neutrino fluxes and detector configurations corresponding to these facilities, we evaluate event rates and sensitivities to a range of observables. Within the SM, argon-based CEvNS detectors enable precision tests of electroweak parameters, including the weak mixing angle, at momentum transfers well below the electroweak scale. We also investigate the sensitivity of these experiments to neutrino electromagnetic properties, such as the magnetic moment and effective charge radius, as well as to possible non-standard neutrino interactions with quarks. Together, these studies highlight the potential of argon-based CEvNS experiments as a clean and versatile platform for precision exploration of non-standard neutrino properties.

Carey, Sam [Wayne State U.] (ORCID:000900029607531↗
DOE OSTI2024

Testing and modeling of CLT-to-CLT joints made using hardwood dowels and screws with varying spacing

With the increasing popularity of mass timber, utilizing large-diameter wood screws and hardwood dowels can offer an easy and cost-efficient way to join structural elements. However, the existing research on hardwood dowels is limited to single fastener joints under monotonic loading. No research is available on the group effects of joints with multiple hardwood dowels. Here, this paper presents the results of 96 monotonic and cyclic single-shear plane CLT-to-CLT joint tests with hardwood dowels (Red Oak and Yellow Birch, 25.4 mm diameter) and self-tapping wood screws (8 mm diameter) in Grand-Fir CLT, where screws were inserted both perpendicular and at a 45-degree angle to the grain for comparison to hardwood dowels. The study evaluated the mechanical properties of the joints by changing fastener spacing. Mechanical properties such as yielding and peak strength, displacement, elastic and yielding stiffness, and ductility were evaluated. The results showed equal strength and stiffness properties for hardwood dowels comparable with large-diameter screws. Additionally, the hardwood dowel joints demonstrated moderate ductility properties. The strength and stiffness of the tested joints were compared with analytical equations found in the literature by considering the group action factor of the fasteners. Lastly, a nonlinear force-displacement model was presented and compared to the experimental results and analytical calculations, which gives the possibility to predict and optimize the mechanical behavior of hardwood dowel and screw joints with varying numbers and spacing between fasteners.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A comprehensive review on valorization of chestnut processing wastes into bio‐based composites and bioplastics

Abstract This review examines the characterization and utilization of chestnut processing wastes (35%) in the production of bioplastics and biocomposites. In this review, a Web of Science search without any publishing year restriction on the biochemical compositions of all the components of Castanea sativa . The obtaining of bioplastics and biocomposites based on C. sativa was reviewed. First, it highlights the biochemical composition and antioxidant properties of chestnut fruit, shell, burrs, leaves, flowers, and wood focusing on the most important compounds, such as phenolic acids, flavonoids, carbohydrates, Klason lignin, cellulose, and glucan, which can enhance the properties of these materials. Then the review covers using several chestnut extracts and fillers in bioplastics production through solvent casting technique. The mechanical, structural, bioactive properties, and moisture content were optimized through the composition and production. The color, UV absorption, antioxidant, and antimicrobial activity were also discussed. Biocomposites reinforced with chestnut burs, shells, or wood flour increased the intended properties. The enhancements in tensile strength, elastic modulus, and the effects of a pre‐treatment were evaluated. Additionally, it discusses material recovery, recycling, and reuse, particularly how it affects the biodegradability of composites incorporating chestnut waste residues. Highlights Chestnut fruit, shells, and burrs are rich in starch, lignin, and cellulose. The waste of chestnut processing can be used in bioplastics and biocomposites. Chestnut‐based films and biocomposites exhibit promising mechanical properties. The antimicrobial activity, making films, and composites proper for food packaging. New techniques boost performance, offering alternatives to conventional plastics.

Silva, Simão B. [REQUIMTE/LAQV, ISEP, Polytechnic ↗
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↗