Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interfaces rotating interfaces”

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 OSTI2026

Epitaxy of Beta-Ga2O3 on Highly-Offcut (Greater than 10 Degrees) Substrates

Beta-Ga2O3 has emerged as a leading candidate for next-generation power electronics, radio frequency (RF) switches, and extreme environment electronics due to a wide band gap (4.6 - 4.9 eV), high dopability (approximately 40 meV activation energy for an isolated silicon donor), and melt growth characteristics resulting in commercially available 4-inch substrates and commercial demonstrations of 6-inch substrates by multiple techniques. The (100) surface of Ga2O3 is highly desirable from a device and epitaxy standpoint - bulk growth of (100) material is more scalable than (010), the surface is nearly lattice-matched to p-type partner NiO, and Al2O3 incorporates at higher concentrations without phase separation. However, the epitaxial growth rate on (100) surfaces is less than 10% of other faces due to weak bonding and favorable desorption. Recent demonstrations have shown growth rate improvements from 0.4 nm/min to 1.5 nm/min by growing on (100) wafers that are offcut 6 degrees in the -c direction.1 These films show step-flow growth from (-201) step-edges and high electron mobility. Despite these exciting results, offcuts greater than 6 degrees have not been explored due to the waste associated with grinding and polishing large offcuts. In this talk we will discuss the molecular beam epitaxy (MBE) growth and properties of Beta-Ga2O3 grown on (100) substrates offcut in the -c direction up to 13.4 degrees. These large offcuts are enabled by edge-fed film-defined growth (EFG) where the offcut is grown into the surface by pulling the crystal through the EFG die with the seed crystal rotated by the desired offcut angle. We will demonstrate that 13.4 degrees offcut substrates still exhibit a terraced (100) surface, and that a >10x increase (4.8 nm/min) in growth rate is achieved. As previously reported on lower offcuts, we observe 100% reversal of substrate twin domains around the (001) direction at the substrate-epilayer interface. We will discuss electrical properties including record-low (by MBE) unintentional doping densities of < 5E15 cm-3.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular beam scattering from flat jets of liquid dodecane and water

Abstract Molecular beam experiments in which gas molecules are scattered from liquids provide detailed, microscopic perspectives on the gas–liquid interface. Extending these methods to volatile liquids while maintaining the ability to measure product energy and angular distributions presents a significant challenge. The incorporation of flat liquid jets into molecular beam scattering experiments in our laboratory has allowed us to demonstrate their utility in uncovering dynamics in this complex chemical environment. Here, we summarize recent work on the evaporation and scattering of Ne, CD 4 , ND 3 , and D 2 O from a dodecane flat liquid jet and present first results on the evaporation and scattering of Ar from a cold salty water jet. In the evaporation experiments, Maxwell–Boltzmann flux distributions with a cos θ angular distribution are observed. Scattering experiments reveal both impulsive scattering (IS) and trapping followed by thermal desorption (TD). Super‐specular scattering is observed for all four species scattered from dodecane and is attributed to anisotropic momentum transfer to the liquid surface. In the IS channel, rotational excitation of the polyatomic scatterers is a significant energy sink, and these species accommodate more readily on the dodecane surface compared to Ne. Our preliminary results on cold salty water jets suggest that Ar atoms undergo some vapor‐phase collisions when evaporating from the liquid surface. Initial scattering experiments characterize the mechanisms of Ar interacting with an aqueous jet, allowing for comparison to dodecane systems. Key Points Molecular beam scattering from flat liquid jets is a powerful technique to elucidate mechanistic detail at the gas–liquid interface. Previous dodecane scattering experiments have uncovered angularly‐resolved thermal desorption fractions and energy transfer at the interface for several small molecule scatterers. Preliminary results on scattering from cold salty water reveal mechanisms of interaction between argon and an aqueous jet.

Yang, Walt↗
DOE OSTI2025

A fully-integrated lattice Boltzmann method for fluid–structure interaction

Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Pseudotunnel Magnetoresistance in Twisted van der Waals Fe 3 GeTe 2 Homojunctions

Abstract Twistronics, a novel engineering approach involving the alignment of van der Waals (vdW) integrated two‐dimensional materials at specific angles, has recently attracted significant attention. Novel nontrivial phenomena have been demonstrated in twisted vdW junctions (the so‐called magic angle), such as unconventional superconductivity, topological phases, and magnetism. However, there have been only few reports on integrated vdW layers with large twist angles θ t , such as twisted interfacial Josephson junctions using high‐temperature superconductors. Herein, vdW homojunctions of the thin‐magnetic flakes, Fe 3 GeTe 2 (FGT), with large θ t ranging from 0° to 90°, without inserting any tunnel barriers are assembled. Nevertheless, these vdW homojunctions exhibit tunnel‐magnetoresistance (TMR) like behavior (pseudo‐TMR (PTMR) effect) with the ratios highly sensitive to the θ t values, revealing that the vdW gap at the junction interface between the twisted FGT layers behaves like a tunnel barrier and the θ t serves a control parameter for PTMR by drastically varying magnitudes of the lattice‐mismatch and the subsequent appearance of antiferromagnetic (AFM) spin alignment. First‐principles calculations considering vacuum gaps indicate strong dependence of TMR on the θ t driven by the sixfold screw rotational symmetry of bulk FGT. The present homojunctions hold promise as a platform for novel AFM spin‐dependent phenomena and spintronic applications.

Obata, Reiji [Faculty of Science and Engineering A↗
DOE OSTI2024

Synergistic Effects of the Electric Field Induced by Imidazolium Rotation and Hydrogen Bonding in Electrocatalysis of CO 2

The roles of the ionic liquid, 1-ethyl-3-methylimidazolium tetrafluoroborate ([EMIM][BF₄]), and water in controlling the mechanism, energetics, and electrocatalytic activity of CO₂ reduction to CO on silver in non-aqueous electrolytes were investigated. The first electron transfer occurs to CO₂ at reduced overpotentials when it is trapped between the planes of [EMIM]⁺ ring and electrode surface due to cation re-orientation as determined from voltammetry, in-situ surface-enhanced Raman spectroscopy, and density functional theory calculations. Within this interface, water up to 0.5 M does not induce significant faradaic activity, opposing the notion of it as a free proton source. Instead, water acts as a hydrogen bond donor and the proton is sourced from [EMIM]⁺. Furthermore, this study demonstrates that alcohols with varying acidity tune the hydrogen bonding network in the interfacial microenvironment to lower the energetics required for CO₂ reduction. Here, the hydrogen bonding suppresses the formation of inactive carboxylate species, thus preserving the catalytic activity of [EMIM]⁺. The ability to tune the hydrogen bonding network opens new avenues for advancing IL-mediated electrocatalytic reactions in non-aqueous electrolytes.

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

How Silica Surface Chemistry Modulates Interfacial Water: Insights from Machine Learning Molecular Dynamics

Controlling water structure and dynamics at silica interfaces are central to a wide range of technologies, including protective oxide layers for solar water splitting and nanoporous membranes. In this work, we develop a machine learning interatomic potential, trained via active learning, to achieve ab initio accuracy for water confined between hydroxylated silica surfaces over a range of silanol coverages and slit widths. We find that partially hydroxylated surfaces (50 and 75% OH) support stronger water−surface hydrogen bonding and more extended interfacial density profiles than fully hydroxylated (100% OH) surfaces, indicating that increasing OH coverage does not necessarily strengthen interfacial hydrogenbond networks. Translational diffusion decreases approximately linearly with slit width and OH coverage, whereas rotational dynamics respond nonlinearly. In particular, at the smallest slit width of 5 Å, 75% OH coverage produces an enhanced local tetrahedral ordered interfacial network that strongly suppresses reorientation, while 100% coverage yields a crowded, disordered interfacial layer that also hinders rotation. In contrast, the 50% OH coverage is sufficiently sparse that it does not markedly alter water structure or dynamics under confinement. These results show that coupled control of pore size and surface chemistry enables nonlinear tuning of interfacial water structure and transport, providing a design strategy for optimizing porous silica for either enhanced interfacial stability and controlled reactivity or rapid and selective transport.

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

A Guide to Current Limiting and Stability with Grid-Forming Inverters

The electric power grid is changing. For many decades, the synchronous generator - an electromechanical device invented in 1887 - has been the workhorse and backbone of power grids across the globe. Our entire infrastructure is built around it. But times change, and so does the power grid. Since the early 21st century, we have seen a gradual shift in modern power grids away from synchronous generators to ones dominated by power electronic inverter-based resources (IBRs). Sources such as photovoltaics, wind turbines, battery storage, fuel cells, and other technologies like high-voltage DC transmission interconnections all rely on an inverter to connect and interface with the grid. They are increasingly being installed on the grid to augment, or even replace, traditional energy sources. This change is a fundamental shift that brings tremendous technical challenges and questions: Can a power grid remain stable with many more inverter-based resources? How do we avoid more blackouts on the grid? How do we keep the grid secure and resilient during disturbances? After all, power electronic inverters are nothing like the big, rotating, iron-and-copper machines that the grid heavily relies on. Many of these questions can be answered by using grid-forming (GFM) inverters, yet many research challenges remain. This document explores GFM inverters and how they can help stabilize the future grid, especially during disturbances and contingencies. It summarizes a two-year research and development fellowship program at NREL. We point interested readers to more detailed works developed during the project along the way. Let's dive in

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

ACDC (Automated Campbell Diagram Code) [SWR-26-042]

This application provides a web-based graphical user interface to generating Campbell Diagrams and visualizing mode shapes for OpenFAST turbine models. Determining the aeroelastic stability and dynamic characteristics of wind turbines is a critical step in turbine design and analysis. Historically, extracting natural frequencies and mode shapes from OpenFAST—the industry-standard whole-turbine simulation code—has been a fragmented and tedious process. It required manual model configuration, command-line linearization execution, and complex post-processing via proprietary scripts to handle rotating-frame dynamics. To address these workflow bottlenecks, we present the Automated Campbell Diagram Code (ACDC), an open-source graphical software tool developed by the National Laboratory of the Rockies (NLR) under the DOE-funded Distributed Wind Aeroelastic Modeling (dWAM) project. ACDC streamlines the end-to-end linearization and stability analysis workflow into a single, intuitive cross-platform application. The software guides users through OpenFAST model configuration, definition of operating points, and the automated execution of steady-state trim and linearization simulations. Under the hood, ACDC automates the complex mathematical post-processing steps required for rotating systems, including Multi-Blade Coordinate (MBC) transformations, eigenanalysis, and advanced modal tracking utilizing the Modal Assurance Criterion (MAC) and spectral clustering. Finally, ACDC processes these results to automatically generate Campbell diagrams and features a robust 3D visualization engine to animate full-system mode shapes. By eliminating the reliance on external post-processing environments and manual data manipulation, ACDC significantly accelerates dynamic analysis and lowers the barrier to entry for wind energy researchers and engineers.

Summerville, Brent [National Laboratory of the Roc↗
DOE OSTI2025

Quantifying Twist Angles in Cuprate Heterostructures with Anisotropic Raman Signatures

Artificially engineered twisted van der Waals (vdW) heterostructures have unlocked new pathways for exploring emergent quantum phenomena and strongly correlated electronic states. Many of these phenomena are highly sensitive to the twist angle, which can be deliberately tuned to tailor the interlayer interactions. This makes the twist angle a critical tunable parameter, emphasizing the need for precise control and accurate characterization during device fabrication. In particular, twisted cuprate heterostructures based on Bi 2 Sr 2 CaCu 2 O 8 + x (BSCCO) have demonstrated angle-dependent superconducting properties, positioning the twist angle as a key tunable parameter. However, the twisted interface is highly unstable under ambient conditions and vulnerable to damage from conventional characterization tools such as electron microscopy or scanning probe techniques. In this work, a fully non-invasive, polarization-resolved Raman spectroscopy approach is introduced for determining twist angles in artificially stacked BSCCO heterostructures. By analyzing twist-dependent anisotropic vibrational Raman modes, particularly utilizing the out-of-plane A 1g vibrational mode of Bi/Sr at ≈116 cm −1 , clear optical fingerprints of the rotational misalignment between cuprate layers are identified. The high-resolution confocal Raman setup, equipped with polarization control and RayShield filtering down to 10 cm −1 , allows for reliable and reproducible measurements without compromising the material's structural integrity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Structural characterization of influenza group 1 chimeric hemagglutinins as broad vaccine immunogens

Chimeric hemagglutinins (cHA) appear to be promising for the design and development of universal influenza vaccines. Influenza A group 1 cHAs, cH5/1, cH8/1, and cH11/1, comprising an H1 stem attached to either an H5, H8, or H11 globular head, have been used sequentially as vaccine immunogens in human clinical trials and induced high levels of broadly protective antibodies. Using X-ray crystallography and negative-stain electron microscopy, we determined structures of cH5/1, cH8/1, and cH11/1 HAs in their apo (unliganded) and antibody Fab-bound states. Stem-reactive antibodies 3E1 and 31.b.09 recognize their cognate epitopes in cH5/1, cH8/1, and cH11/1 HAs. However, with cH5/1, the head domains are rotated by 35 to 45° around the threefold axis of the HA trimer compared to native HA with a more splayed-open conformation at the stem base. cH11/1 with 3E1 is structurally more native-like but resembles cH5/1 with 31.b.09, whereas cH8/1 with 31.b.09 exhibited a range of closed-to-open stem configurations with some separation of head and stem domains. Furthermore, all of these group 1 cHAs effectively bound a broad head trimer interface antibody and other broad stem antibodies. Thus, the cHAs exhibit structural plasticity without compromising the stem and head trimer interface epitopes for elicitation of influenza A group 1 cross-reactive antibodies.

Science & Technology - Other Topics↗
DOE OSTI2025

Unified ELM Suppression on KSTAR and DIII-D via Adaptive Feedback Control Strategies

This paper reports on the extension of our amplitude-based resonant magnetic perturbation (RMP) edge localized mode (ELM) controller to support phasing control (relative toroidal phases of RMP waveforms between rows of coils), multiple toroidal mode numbers, and new ‘jump’ and ‘probing’ strategies, all deployed on KSTAR and DIII-D. By treating the control algorithm as device-independent and adjusting only the real-time interfaces to sensors and power supplies, we have confirmed that the same finite state machine—based feedback logic can be ported between machines with minor modifications. In experiments using n = 2 RMPs on KSTAR and n = 3 on DIII-D, the controller successfully modulated RMP amplitudes in real time to sustain ELM suppression while minimizing confinement degradation. Phasing control broadened the suppression window, as it permitted the system to avoid locked-mode regions and safely access ELM-free conditions. A rotating RMP phasing scheme, integrated into the same framework, distributes divertor heat loads more uniformly, making it a promising strategy for protecting plasma-facing components during long discharges. New ‘jump’ and ‘probing’ techniques demonstrate the possibility for the controller to preempt imminent ELMs and refine the minimum required RMP amplitude without returning to ELMy conditions. Taken together, these upgrades enable extended ELM-free operation while mitigating confinement degradation, providing a practical framework for real-time ELM control in future high-performance tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Active-noise-induced dynamic clustering of passive colloidal particles

Active fluids generate spontaneous, often chaotic mesoscale flows. Harnessing these flows to drive soft materials embedded within an active fluid into structures with controlled length scales and lifetimes is a key challenge at the interface between the fields of active matter and nonequilibrium self-assembly. Here, we present a simple and efficient computational approach to model soft materials advected by active fluids, by simulating particles moving in a spatiotemporally correlated noise field. To illustrate our approach, we simulate the dynamical self-organization of repulsive colloids within such an active noise field. The colloids form structures whose sizes and dynamics can be tuned by the correlation time and length of the active fluid, and range from small rotating droplets to clusters with internal flows and system-spanning sizes that vastly exceed the active correlation length. Our results explain how the interplay between active fluid time and length scales and emergent driven assembly can be used to rationally design functional assemblies. More broadly, our approach can be used to efficiently simulate diverse active fluids and other systems with spatiotemporally correlated noise.

Brownian dynamics↗
DOE OSTI2025

A Geometric Volume of Fluid-Based Multiphase Flow Solver Extension to the Reacting Flow Solver, PeleLM

A new algorithm is presented to simulate multiphase flows with surface tension in a pathway for spray combustion simulation. The algorithm combines capabilities from two open-source packages, including the interface reconstruction library (IRL), a library of computational geometry routines to enable the volume of fluid (VOF) method, and PeleLM, a solver for the reacting Navier-Stokes equations. Additionally, surface tension is implemented using the continuum surface force (CSF) model with an improved height function technique in the volume fraction field. Spurious errors in volume fraction arising from our combined strategy are corrected through a topology-based method that improves both numerical stability and accuracy. Multiple validation simulations are conducted, including (i) translations and rotations of Zalesak's disk, (ii) a stationary circular droplet with surface tension, (iii) an oscillating elliptical droplet, and (iv) three-dimensional deformation of a spherical droplet. Results indicate that the combined scheme retains the favorable properties of each of the component algorithms.

42 ENGINEERING↗
DOE OSTI2025

Polarization rotation in a ferroelectric BaTiO 3 film through low-energy He-implantation

Domain engineering in ferroelectric thin films is crucial for next-generation microelectronic and photonic technologies. Here, a method is demonstrated to precisely control domain configurations in BaTiO 3 thin films through low-energy He ion implantation. The approach transforms a mixed ferroelectric domain state with significant in-plane polarization into a uniform out-of-plane tetragonal phase by selectively modifying the strain state in the film’s top region. This structural transition significantly improves domain homogeneity and reduces polarization imprint, leading to symmetric ferroelectric switching characteristics. The demonstrated ability to manipulate ferroelectric domains post-growth enables tailored functional properties without compromising the coherently strained bottom interface. The method’s compatibility with semiconductor processing and ability to selectively modify specific regions make it particularly promising for practical implementation in integrated devices. This work establishes a versatile approach for strain-mediated domain engineering that could be extended to a wide range of ferroelectric systems, providing new opportunities for memory, sensing, and photonic applications where precise control of polarization states is essential.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Polarization rotation in a ferroelectric BaTiO$_3$ film through low-energy He-implantation

Domain engineering in ferroelectric thin films is crucial for next-generation microelectronic and photonic technologies. Here, a method is demonstrated to precisely control domain configurations in BaTiO$_3$ thin films through low-energy He ion implantation. The approach transforms a mixed ferroelectric domain state with significant in-plane polarization into a uniform out-of-plane tetragonal phase by selectively modifying the strain state in the film's top region. This structural transition significantly improves domain homogeneity and reduces polarization imprint, leading to symmetric ferroelectric switching characteristics. The demonstrated ability to manipulate ferroelectric domains post-growth enables tailored functional properties without compromising the coherently strained bottom interface. The method's compatibility with semiconductor processing and ability to selectively modify specific regions make it particularly promising for practical implementation in integrated devices. This work establishes a versatile approach for strain-mediated domain engineering that could be extended to a wide range of ferroelectric systems, providing new opportunities for memory, sensing, and photonic applications where precise control of polarization states is essential.

Luo, Liang [Ames Laboratory (AMES), Ames, IA (Unit↗
DOE OSTI2026

First-principles DFT modeling of nitrobenzene adsorption on the Ag(111) surface at varying monolayer coverages

Here, the adsorption behavior of nitrobenzene on the Ag(111) surface as a function of coverage is investigated using density functional theory. Adsorption energies and optimized geometries are analyzed together with isolated intermolecular interaction calculations and electronic structure analysis, including Bader charge partitioning and projected density of states, to disentangle the roles of adsorbate–surface bonding and through-space adsorbate–adsorbate interactions. At low and intermediate coverages (θ = 1/9 and θ = 2/9), bidentate adsorption at top sites is favored due to strong adsorbate–surface interactions, with additional stabilization at θ = 2/9 arising from favorable intermolecular separations. At higher coverage (θ = 1/3), bidentate adsorption is destabilized by strong intermolecular repulsion, and monodentate adsorption becomes energetically preferred as rotational freedom allows more favorable intermolecular spacing. Charge density differences, Bader charge, and density of states analyses show that charge transfer from the Ag surface is localized primarily on the nitro group and increases with coverage and adsorption denticity, although this increase does not directly correlate with adsorption strength at high coverage due to competing intermolecular interactions. These results demonstrate that surface coverage can induce a transition in preferred adsorption denticity driven by intermolecular interactions, highlighting the importance of adsorbate packing in organic molecule adsorption at metal interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2024

LibERI—A portable and performant multi-GPU accelerated library for electron repulsion integrals via OpenMP offloading and standard language parallelism

A portable and performant graphics processing unit (GPU)-accelerated library for electron repulsion integral (ERI) evaluation, named LibERI, has been developed and implemented via directive-based (e.g., OpenMP and OpenACC) and standard language parallelism (e.g., Fortran DO CONCURRENT). Offloaded ERIs consist of integrals over low and high contraction s, p, and d functions using the rotated-axis and Rys quadrature methods. GPU codes are factorized based on previous developments with two layers of integral screening and quartet presorting. In this work, the density screening is moved to the GPU to enhance the computational efficacy for large molecular systems. Here, the L-shells in the Pople basis set are also separated into pure S and P shells to increase the ERI homogeneity and reduce atomic operations and the memory footprint. LibERI is compatible with any quantum chemistry drivers supporting the MolSSI Driver Interface. Benchmark calculations of LibERI interfaced with the GAMESS software package were carried out on various GPU architectures and molecular systems. The results show that the LibERI performance is comparable to other state-of-the-art GPU-accelerated codes (e.g., TeraChem and GMSHPC) and, in some cases, outperforms conventionally developed ERI CUDA kernels (e.g., QUICK) while fully maintaining portability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗