SEARCH · Engineering Papers

Results for “dynamic loading”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47

All sources NASA NTRS DOE OSTI

DOE OSTI2025

First-Principles Molecular Dynamics Simulations of Ammonia Adsorption onto MFI Zeolite Nanosheets

MFI zeolite nanosheet membranes are promising candidates for ammonia separation from nitrogen and hydrogen gases, yet questions remain on the origin of their high selectivity. Silanols, Si–OH, are present in high concentration at the surface of zeolite nanosheets, and force-field-based simulations indicate that surface adsorption at the silanols contributes to selectivity. Acidic silanol groups can chemically react with ammonia via transfer of a proton to form ammonium ions, which may further contribute to the ability of zeolite nanosheet membranes to separate ammonia from other gases. In this work, we used first-principles molecular dynamics techniques to simulate ammonia’s behavior within stacked MFI zeolite nanosheets. We found that at 523 K and a loading corresponding to 35 bar, conditions desired for the ammonia separation, about 30% of ammonia reacts with surface silanols. Our work explored H-bonding and proton transfer within this system.

ammonia↗

DOE OSTI2025

Distribution System Blackstart and Restoration Using DERs and Dynamically Formed Microgrids

Extreme weather events have led to long-duration outages in the distribution system (DS), necessitating novel approaches to blackstart and restore the system. Existing blackstart solutions utilize blackstart units to establish multiple microgrids (MGs), sequentially energize non-blackstart units, and restore loads. However, these approaches often result in isolated MGs. In DERs-aided blackstart, the continuous operation of these MGs is limited by the finite energy capacity of commonly used blackstart units like battery energy storage (BES)-based gridforming inverters (GFMIs). To address this issue, this article proposes a holistic blackstart and restoration framework that incorporates synchronization between dynamic MGs and the entire DS with the transmission grid (TG). To support synchronization, we leveraged virtual synchronous generator-based control for GFMIs to estimate their frequency response to load pick-up events using only initial/final quasi-steady-state points. Subsequently, a synchronization switching condition is developed to model synchronizing switches, aligning them seamlessly with a linearized branch flow problem. Finally, we designed a bottomup blackstart and restoration framework that considers the switching structure of the DS, energizing/synchronizing switches, DERs with grid-following inverters, and BES-based GFMIs with frequency security constraints. In conclusion, the proposed framework is validated in IEEE-123-bus system, considering cases with two and four GFMIs under various TG recovery instants.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Understanding the dynamic evolution of atomically dispersed Cu catalyst for CO 2 electrochemical conversion using integrated XANES analysis and mechanistic studies

Direct electrochemical conversion of CO 2 to ethanol (CH 3 CH 2 OH) offers a promising strategy to lower CO 2 emission while storing energy from renewable electricity. Our recent study reported a carbon-supported atomically dispersed Cu catalyst that achieved the highest reported selectivity for CH 3 CH 2 OH formation (91%) at a relatively low potential (-0.6 V), however, the active site structure that is responsible for such high activity and selectivity has yet to be understood. Here, in this paper, we demonstrate a computational investigation combining X-ray absorption near edge structure (XANES) simulations and a mechanistic study via density functional theory (DFT) to understand the catalyst structures of this Cu catalyst during electrocatalysis and the corresponding reaction mechanisms of the key products. An integrated computational and experimental XANES analysis depicted the dynamic evolution of the catalytic site during electrocatalysis. The as-prepared, atomically dispersed Cu catalyst aggregates and forms metallic clusters/nanoparticles under electrochemical condition, which then break down to smaller oxidized clusters after electrocatalysis. The formed Cu clusters/nanoparticles showed distinct catalytic activity and selectivity as a function of particle size based on the mechanistic investigation using DFT, which is consistent with experimental observations for catalyst samples with different Cu loadings. This comprehensive study which combines experimental and computational XANES investigation, mechanistic study via DFT calculations, and experimental performance of the catalysts, provides unprecedented dynamic and mechanistic insights into the supported atomically dispersed metal catalysts for CO 2 reduction. Such strategy and details gained can further guide discovery of novel catalyst materials for CO 2 electrochemical reduction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Performance Demonstration of an Occupancy Sensor-enabled Integrated Solution for Commercial Buildings

Traditionally, a single-loop fixed-gain controller is applied to supply fan (SF) and cooling coil (CC) valve controls while a fixed-damper position control is applied to outdoor air (OA) damper control at air handling units (AHUs) in commercial buildings. With the increasing application of occupancy sensors, the information generated by occupancy sensors is applied to not only reduce the electricity loads from lighting and controllable plug loads, but also reset OA intake and minimum supply airflow setpoints. Meanwhile, these intermittent operation actions greatly elevate the dynamics of AHU systems, which may introduce unstable SF and CC valve operations and inaccurate OA flow control at AHUs and consequently degrade maximum energy efficiency gains. With virtual fan and valve flow meter technologies, two advanced controls, including cascade control and gain scheduling control, can be implemented on both the SF and CC valve, and an advanced control using a virtual OA flow meter can be implemented on OA damper integrated with occupancy sensors. The goal of this project is to demonstrate the savings, cost, and performance of an integrated solution that integrates the three advanced HVAC controls with occupancy sensors to allow accurate and stable AHU operations in real buildings. The project objectives are to: 1) develop and validate an advance SF control algorithm; 2) develop and validate an advanced CC valve control algorithm; 3) validate an algorithm to implement a virtual OA flow meter; and 4) demonstrate the savings, cost, and performance of the integrated solution in real buildings. The technical approaches are to: 1) select a test system at the University of Oklahoma; 2) develop and implement the algorithms of advanced SF and CC valve controls and validate the performance; 3) develop and implement the advanced OA control using a virtual OA flow meter and validate the performance; 4) demonstrate the savings, cost, and performance of the proposed integrated solution with and without three advanced HVAC controls; and 5) disseminate the project results through publications and presentations. For the SF control, both the gain scheduling and cascade controls can improve the fan energy performance by reducing the fan power during the transient period and the fan control performance at lower speeds by reducing fan speed variation. Moreover, the gain scheduling control provides a simple and low-cost solution and is recommended. The fan power savings can reach 30% in a transient period. For the CC valve control, the gain scheduling control can considerably reduce the supply air temperature oscillation range and frequency under both higher and lower load conditions and the control valve response is much more stable. As a result, the gain scheduling control is recommended. The projected pump energy consumption can be reduced by 68.5%. With the developed virtual OA flow meter, the OA can be accurately controlled at its setpoint, which is determined based on the actual number of occupants in the building provided by occupancy sensors. The RMSE of the proposed OA control is 15.9 L/s. The energy data shows that the fan power and CC cooling energy were significantly reduced. On the other hand, the energy savings majorly results from the occupancy sensors and the energy savings by the advanced HVAC controls is minimal because that the controllers in the test AHU were tuned with very slow response. An annual technical savings potential is estimated as 0.5 quads in the commercial sector.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Towards a predictive, physics-based friction model for the dynamics of jointed structures

Bolted connections are ubiquitous in mechanical designs and pose a significant challenge to understanding and predicting the vibration response of assembled structures. The present paper develops a physics-based rough contact model of the frictional interactions within a joint. This model sums over the probable interactions of asperities – defined as locally maximum surface features – to determine the contact forces. Here, the tangential contact forces vary smoothly between sticking and slipping and allow the model to better capture the qualitative trends of experimental amplitude dependent frequency and damping than previous studies. Furthermore, the novel model is generalized to allow for arbitrarily varying normal pressure including potential separation to better represent the interfacial dynamics. This includes developing a new, computationally tractable approximation to the analytical Mindlin partial slip solution for tangential loading of contacting spheres. Further, the results highlight the importance of accurately characterizing the as-built topology of the interface, the plastic behavior of the contacting asperities, the relevant length scale of asperities, and the eccentricity of asperities. A predictive friction coefficient based on plasticity provides a poor match to experiments, so fitting the friction coefficient is also considered. Numerical results are compared to experiments on the Brake-Reuß Beam to assess the predictive potential of the models. While blind predictions over-predict the slip limit, the current model presents a significant improvement in physics-based modeling and highlights areas for ongoing research.

42 ENGINEERING↗

DOE OSTI2020

Influence of Silane Coupling Agents on Filler Network Structure and Stress-Induced Particle Rearrangement in Elastomer Nanocomposites

Filled rubber materials are key in many technologies having broad impact on the economy and sustainability, the most obvious being tire technology. Adding filler dramatically improves the strength of rubber by reinforcement and tailoring the type of filler and the chemistry of the interface between filler and rubber matrix is important for optimizing performance metrics such as fuel efficiency. In a highly loaded, silica-filled, crosslinked model rubber closely mimicking commercial materials, both the filler network structure and the dynamics of the silica filler particles change when the silica surface is modified with silane coupling agents. Reduction in size scales characteristic of the structure are quantified using ultra-small-angle X-ray scattering (USAXS) measurements and the particle dynamics probed with X-ray photon correlation spectroscopy (XPCS). While the structure averaged over the scattering volume changes little with aging after step strain, the dynamics slow appreciably in a manner that varies with the treatment of the silica filler. The evolution of filler particle dynamics depends on the chemical functionality at the silica surface and observing these differences suggests a way of thinking about the origins of hysteresis in nanoparticle-reinforced rubbers. These microscopic filler dynamics are correlated with the macroscopic stress relaxation of the filled materials. The combination of static and dynamic X-ray scattering techniques with rheological measurements is a powerful approach for elucidating the microscopic mechanisms of rubber reinforcement.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Analytical Identification Method of Generalized Short‐Circuit Ratio Using Phasor Measurement Units

This paper introduces a novel analytical approach for the identification of the admittance matrix and the generalized short-circuit ratio (gSCR) in power systems integrated with renewable energy sources. The proposed method leverages voltage and current measurements from phasor measurement units (PMUs) to construct a least squares objective function, which is then solved using matrix calculus and partial derivatives. Unlike conventional optimization algorithms, this approach provides an analytical solution that substantially reduces data requirements, enabling the efficient and accurate identification of the gSCR with smaller datasets. Additionally, its fixed computational complexity allows for real-time updates as new data are collected, ensuring continuous refinement of the system of equations and enabling rapid, precise gSCR calculations. The method also exhibits strong robustness against measurement noise, making it well-suited for practical applications in dynamic power systems. The combination of reduced data requirements, real-time adaptability, noise robustness and fixed computational load establishes this method as a highly efficient and reliable tool for real-time power system stability analysis. Case studies on an EPRI 36-bus system demonstrate the method's effectiveness, highlighting its accuracy in closely matching true gSCR values, even under diverse disturbances and noisy conditions.

Han, Zelei [Hohai University, Nanjing (China)] (OR↗

DOE OSTI2024

Advanced co-simulation framework for assessing the interplay between occupant behaviors and demand flexibility in commercial buildings

With buildings contributing significantly to electricity usage, enabling demand flexibility becomes a challenge, especially when accounting for occupant comfort. This study introduces an innovative co-simulation framework integrating multiple models: heating, ventilation, and air conditioning (HVAC) system, building zone load, indoor airflow, supervisory control, and occupant comfort and behavior. Uniquely, this framework allows for a comprehensive and dynamic analysis of building systems and occupant interactions in demand response events. Using this framework, we conducted a case study using a typical small office building model. Specifically, we focused on three areas: (1) the impact of indoor airflow modeling on energy use, occupant comfort, and behaviors forecasting, (2) the impact of occupant behaviors on demand flexibility, and (3) occupant comfort and behaviors under demand response events. Key performance indicators such as energy use, flexibility factor, durations of occupant discomfort and occupant behaviors were analyzed. Our findings indicated variations in energy usage and occupant comfort within demand flexibility events, marked by uncertainty boundaries, with variability in demand shedding up to 57.9%. Here, we concluded that this framework is suitable for analyzing typical commercial buildings and their HVAC systems in terms of demand flexibility potential under the impact of occupant behaviors.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

Improved Control Strategy of Grid Interface for EV High-Power Dynamic Wireless Charging

This paper proposes a control strategy for the grid interface converter in high-power dynamic wireless charging system (DWCS) to address two issues on distribution network integration. Due to the unique pulsating load profile of DWCS, load transient response capability is critical for the grid interface to maintain the dc-bus voltage stable. Besides, the inherent unbalanced situation of distribution network would lead to 2nd-order oscillations on the dc-bus voltage, which would further affect the stable operation of the entire system. In this paper, the DWCS model is developed, and the relationship between the dc-bus voltage and the input/output power is analyzed. Based on the developed model, a control strategy based on direct power control is presented. Both simulation results and hardware-in-the-loop (HIL) results demonstrate that the proposed control strategy not only improves load transient response capability, but also eliminates the 2nd-order oscillations on the dc-bus voltage under imbalanced distribution network conditions.

Zeng, Rong↗

DOE OSTI2025

Integrating Energy-Efficient Computing with Computational Research to Accelerate Energy Technology

NREL's computational sciences center hosts the largest high performance computing (HPC) capabilities dedicated to energy research while functioning as a living laboratory for energy-efficient computing. NREL's HPC capabilities support the research needs of the Department of Energy's Office of Energy Efficiency and Renewable Energy (EERE). In ten years of operation, HPC use in EERE-sponsored research has grown by a factor of 30, including work in electricity generation, energy efficiency, transportation, and energy system modeling. This paper analyzes this research portfolio, providing examples of individual use cases. The paper documents NREL's history of operating one of the world's most energy-efficient data centers while examining pathways to reduce economic and environmental impact beyond reduction of Power Usage Efficiency (PUE). This paper concludes by examining the unique opportunities created for accelerating improvements in data center efficiency created by combining an HPC system dedicated to energy research and a research program in energy-efficient computing.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Semi-Annual Report for Horizontal Compact High Temperature Gas Reactor (HC-HTGR) Development during Performance Period April 2023 – September 2023

Horizontal Compact High Temperature Gas Reactor (HC-HTGR) is being designed by a multi-disciplinary team of nuclear, mechanical, and structural engineers under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed HC-HTGR in 3 years and support its commercialization as a safe, low-cost HTGR. Argonne National Laboratory (Argonne) is responsible for the design and analysis of the reactor cavity cooling system (RCCS) as a safety system for passive decay heat removal of the reactor concept. Additionally, Argonne is providing analysis of the primary coolant system to ensure temperatures within the core remain below safety margins during steady-state and potential accident scenarios. This fourth semi-annual report summarized the progress made at Argonne on the two tasks during the second half of FY23. As a part of the RCCS design task, recent efforts have been made to complete a conceptual design of the RCCS for the HC-HTGR, including the design update of the water panel and system configuration favorable in point of view of fabrication and system operation. Design calculations were conducted under various heat load conditions to validate the system design. Transient simulations using RELAP5-3D were conducted to investigate system dynamics under transients of interest and to evaluate the system performance in the design condition. The results demonstrated the overall system feasibility that the RCCS design maintains structures temperatures lower than maximum allowable temperature with sufficient system inventory without any active heat sink. In the primary system thermal hydraulics task, the preliminary analysis was performed for a long term pressurized conduction cooldown (PCC) transient. This analysis used a combination of a fully resolved and coarse homogenized mesh to predict the temperature distribution for the steady-state initial condition and the PCC transient. The steady-state initial condition was determined using a fully resolved full core model with 3D solid to 1D fluid coupling. The fully resolved mesh was also used to model the first 20 seconds of the PCC. The temperature difference between fuel pins and the graphite matrix becomes minimal and the dominant heat transfer shifts to a larger scale radially towards the RCCS. After 20 seconds, a homogenized coarse mesh is used, greatly reducing the computational costs of the model. These results demonstrate that the core is designed to passively remove enough decay heat in a protected loss of primary coolant flow to prevent an unsafe rise of core temperatures.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2026

Development of Wireless Power Technology for Use in Electrochromic Window Retrofit Applications

Electrochromic (EC) windows, often referred to as smart windows, have the ability to adjust their tint in response to electrical signals, thereby controlling the amount of sunlight and heat entering a building. This dynamic modulation can lead to significant energy savings by reducing the reliance on artificial lighting and minimizing the load on heating, ventilation, and air conditioning (HVAC) systems. However, the widespread adoption of electrochromic window technology, especially in existing buildings, faces challenges, particularly concerning the installation cost associated with the need for electrical wiring to power and control the EC glass (tinting). In this project, GuRu Wireless demonstrated that it’s mmWave-based wireless power energy transfer technology can provide the necessary power for operation of electrochromic windows, wirelessly, and address the wiring challenges of electrochromic windows in retrofit applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

A Hybrid Heavy-Duty Diesel Power System for Off-Road Applications - Concept Definition

A multi-year Power System R&D project was initiated with the objective of developing an off-road hybrid heavy-duty concept diesel engine with front end accessory drive-integrated energy storage. This off-road hybrid engine system is expected to deliver 15-20% reduction in fuel consumption over current Tier 4 Final-based diesel engines and consists of a downsized heavy-duty diesel engine containing advanced combustion technologies, capable of elevated peak cylinder pressures and thermal efficiencies, exhaust waste heat recovery via SuperTurbo™ turbocompounding, and hybrid energy recovery through both mechanical (high speed flywheel) and electrical systems. The first year of this project focused on the definition of the hybrid elements using extensive dynamic system simulation over transient work cycles, with hybrid supervisory controls development focusing on energy recovery and transient load assist, in Caterpillar’s DYNASTY™ software environment. Three key off-road applications were the focus of the hybrid concept definition with an aim of understanding the system’s modular capability for the diverse off-road heavy-duty market. Core engine performance 1D and 3D simulations isolated the efficiency contributions from the downsized engine, turbocompounding, and in-cylinder thermal barrier coatings. A fuel consumption improvement range of 14 to 24% was predicted, resulting in successful project progression to the design and experimental validation phase. Furthermore, an overview of the experimental engine and hybrid system status concludes the discussion along with the multi-year project’s next steps.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2026

Measurement-informed Dynamic Aggregation of Distribution Systems

This paper proposes a measurement-informed dynamic aggregation methodology in order to create equivalent representations of distribution systems that are compatible with large-scale transmission analysis. By optimizing an equivalent feeder parameters using time-series measurements of active power, reactive power, and voltage at the Point of Interconnection (POI), the approach yields simplified yet dynamically accurate equivalents. Implemented in PSCAD with models of photovoltaic–battery systems, three-phase motors, and static loads, the method employs hybrid differential evolution and bounded least-squares optimization laying the foundation for for real-time state estimation and optimized sensor placement in distribution networks.

Ahmed, Kazi Ishrak [University of Tennessee, Knoxv↗

DOE OSTI2025

Integration of a grey-box refrigerated case model in EnergyPlus via Python plugin

Commercial buildings, in particular grocery stores (due mainly to their large refrigeration load), provide opportunities for energy cost reductions. Grocery stores could offer substantial load flexibility to the power grid through participation in demand response programs because of their usage patterns and relatively high energy intensity. This load flexibility could come from modifying the control of heating, ventilation, and air conditioning (HVAC) systems, refrigeration systems, or both. Although estimation of the HVAC system’s load flexibility potential is relatively targeted in the literature, estimating load flexibility of refrigeration systems is nascent and has been a challenge, in part because of the lack of proper simulation tools that capture the dynamics in the refrigeration cases. The existing refrigerated case model within EnergyPlus, a whole building energy simulation program, assumes a constant case temperature throughout the simulation period and does not explicitly model the cycling of the compressor serving the refrigerated case. In addition, it does not encompass modeling of temperatures of the product inside the refrigerated case. This difference between modeled and actual operation can be a barrier to the development of demand control algorithm and accurate analysis of load flexibility potential. In this paper, we present a grey-box model for modeling refrigerated cases in grocery stores, which include medium temperature and low temperature. Four cases are modeled; two are low-temperature closed cases and two are medium-temperature cases with one closed and one open. Data from an experimental facility are used to train and test the models. Results demonstrate the efficacy of the grey-box models in predicting the temperatures. This model is integrated into EnergyPlus to capture the dynamic effects of case temperature on the environment and enhance the calculation of sensible and latent heat exchange with the environment (case credits). These enhancements can be leveraged more broadly to model advanced refrigeration controls such as defrost, develop and test unique algorithms that could affect refrigeration interactions with HVAC, and refine store design for any commercial building with refrigeration.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2022

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↗

DOE OSTI2024

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↗