Search NASA⌕ Search

SEARCH · Search NASA

Results for “continuum damage mechanics”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Determinants of the mechanical properties of bones

The mechanical properties of bones are governed by the same principles as those of man-made load-bearing structures, but the organism is able to adapt its bone structure to changes in skeletal loading. In this overview of the determinants of the strength and stiffness of bone, a continuum approach has been taken, in which the behavior of a macroscopic structure depends on its shape and size, and on the mechanical properties of the material within. The latter are assumed to depend on the composition (porosity and mineralization) and organization (trabecular or cortical bone architecture, collagen fiber orientation, fatigue damage) of the bone. The effects of each of these factors are reviewed. Also, the possible means of non-invasively estimating the strength or other mechanical properties of a bone are reviewed, including quantitative computed tomography, photon absorptiometry, and ultrasonic measurements. The best estimates of strength have been obtained with photon absorptiometry and computed tomography, which at best are capable of accounting for 90% of the strength variability in a simple in vitro test, but results from different laboratories have been highly variable.

Non-NASA Center↗
NASA NTRS2008

Multiscale Modeling of Damage Processes in fcc Aluminum: From Atoms to Grains

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, current analysis is limited to small domains and increasing the size of the MD domain quickly presents intractable computational demands. A preferred approach to surmount this computational limitation has been to combine continuum mechanics-based modeling procedures, such as the finite element method (FEM), with MD analyses thereby reducing the region of atomic scale refinement. Such multiscale modeling strategies can be divided into two broad classifications: concurrent multiscale methods that directly incorporate an atomistic domain within a continuum domain and sequential multiscale methods that extract an averaged response from the atomistic simulation for later use as a constitutive model in a continuum analysis.

Glaessgen, E. H.↗
DOE OSTI2020

Functional Photoresists for Energy Applications. Final Report

Monolithic ultralow-density porous bulk materials have recently attracted much interest due to many emerging applications in the areas of catalysis, energy storage and conversion, and thermal insulation. They are also important components of high energy density (HED) and inertial confinement fusion (ICF) targets. However, despite tremendous progress that has been made in the synthesis of porous materials, deterministic and independent control over microscopic architecture, density and composition remain key issues, and their integration in high precision devices requires cost and time-intensive mechanical machining that not only reduces reproducibility by generating debris but also limits the complexity of the 3D shapes that can be realized. In this project, we overcame these limitations by developing a universal templating capability that provides deterministic and independent control over density, composition, architecture, and macroscopic sample shape. This was achieved by developing the technology to 1) 3D print ultrahigh resolution, ultra-high precision polymeric micro-lattice templates, 2) coat these templates with the desired materials, and 3) removing the template (Fig. 1a). Atomic layer deposition (ALD) provides the atomic scale coating thickness accuracy required for precisely controlling density. While this templating approach had been demonstrated in prior work, limitations in suitable photoresists, 3D print technologies, print design, and template removal techniques did not allow the fabrication of millimeter-sized high-precision parts with sub-micron resolution. To enable this technology, we developed 1) two-photon polymerization (TPP) print designs that enable the fabrication of millimeter-sized, mechanically robust polymeric templates with sub-micron resolution and 2) a continuum level TPP printing simulation capability for additional print design guidance; 3) atomistic models to study photoresist polymerization kinetics and network topography, 4) refractive index matched polymeric and preceramic TPP photoresists, and 5) functional TPP photoresists including porous voxel structures and self-immolative polymer photoresist chemistries; and 6) damage free template removal techniques that enable the fabrication of defect-free high-precision low-density foam components. We also developed a templating approach for pure carbon microlattices with a unique tube-in-tube ligament morphology. As a test platform, we pursued the fabrication of foam liners that promise to further increase the neutron yield in indirect drive ICF experiments by improving implosion symmetry control and coupling between the laser and the deuterium-tritium fuel. This application requires fabrication and integration of a ultra-high precision, millimeter-sized, thin-walled (200-400 micrometer thick), low-density (10-30 mg/cc), high atomic number (high Z) cylindrical foam tube into the gold hohlraum of an indirect drive ICF target (Fig. 1b). While our hohlraum liner test case will mainly find application in HED and ICF experiments, the underlying science will also directly apply to previously developed nanoparticle and additive manufacturing technologies and will advance those techniques as well.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A better understanding of the mechanics of borehole breakout utilizing a finite strain gradient-enhanced micropolar continuum model

Borehole breakout denotes the failure in rock mass subjected to drilling, caused by stress concentrations exceeding the material strength. Depending on the material properties, the preexisting in situ stress state, and the borehole dimensions, different types of borehole breakout, such as spiral-shaped breakout or v-shaped breakout, are distinguished in the literature. In the present work, we address the influence of the material properties on the predicted borehole breakout mode in a comprehensive finite element study. To this end, herein we employ a gradient-enhanced micropolar damage-plasticity model based on the Mohr–Coulomb strength criterion, formulated in the finite strain regime, which is calibrated based on experimental results from plane strain compression tests on Red Wildmoor sandstone. In the numerical study, the influence of the in situ stress state, the material friction angle, plastic dilation, post peak residual strength, and the inherent material length scale parameters are investigated. Thereby, we demonstrate that depending on the material parameters, substantially different failure modes, characterized by strongly localized shear bands or diffuse failure zones, are predicted. It is shown that in particular the brittleness of the material in the post peak regime has a major influence on the predicted breakout type. Moreover, a statistical validation of the results is obtained by considering different random field distributions of the initial material strength.

58 GEOSCIENCES↗
DOE OSTI2024

Cohesive phase-field chemo-mechanical simulations of inter- and trans- granular fractures in polycrystalline NMC cathodes via image-based 3D reconstruction

The optimal design and durable utilization of lithium-ion batteries necessitates an objective modeling approach to understand fracture and failure mechanisms. This paper presents a comprehensive chemo-mechanical modeling study focused on elucidating fracture-induced damage and degradation phenomena in the polycrystalline Li $\mathcal{x}$ Ni 0.5 Mn 0.3 Co 0.2 O 2 (NMC532) cathode. An innovative approach that utilizes image-based reconstructed 3D geometry as finite element (FE) mesh input is employed to enhance the precision in capturing the convoluted architecture and morphological features. For accurately representing the intricate crack configurations within the polycrystalline system, we adopted the cohesive phase-field fracture (CPF) model. Through the integration of advanced image-based geometry reconstruction technique and the promising CPF modeling approach, lithium (de)intercalation induced crack evolution (e.g., nucleation, propagation, branching and diverse modes including inter-/trans-(intra-) granular patterns) and the resulting chemical degradation can be precisely captured, which is also compared and validated with numerical predictions using a continuum damage model. In particular, this model predicts fracture induced degradation under varying fracture properties of grain boundaries and charging rates; the conclusion that NMC particles comprised of larger grains are predicted to have less degradation than those with smaller grains can also be drawn. This comprehensive analysis provides valuable insights into the fracture and degradation within polycrystalline NMC cathodes.

25 ENERGY STORAGE↗
NASA NTRS2000

The Significance of Small Cracks in Fatigue Design Concepts as Related to Rotorcraft Metallic Dynamic Components

In this paper the significance of the "small" crack effect as defined in fracture mechanics will be discussed as it relates to life managing rotorcraft dynamic components using the conventional safe-life, the flaw tolerant safe-life, and the damage tolerance design philosophies. These topics will be introduced starting with an explanation of the small-crack theory, then showing how small-crack theory has been used to predict the total fatigue life of fatigue laboratory test coupons with and without flaws, and concluding with how small cracks can affect the crack-growth damage tolerance design philosophy. As stated in this paper the "small" crack effect is defined in fracture mechanics where it has been observed that cracks on the order of 300 microns or less in length will propagate at higher growth rates than long cracks and also will grow at AK values below the long crack AK threshold. The small-crack effect is illustrated herein as resulting from a lack of crack closure and is explained based on continuum mechanics principles using crack-closure concepts in fracture mechanics.

Everett, R. A., Jr.↗
NASA NTRS2019

Progressive Damage and Failure Analysis of Bonded Composite Joints at High Energy Dynamic Impacts

Both wing and fuselage structures utilize bonded composite joints for structural efficiency in modern commercial and military aircraft. To ensure compliance with certification requirements mechanical fasteners are typically used as a failsafe mechanism for appropriate strength in the event of complete stiffener disbond. However, the use of fasteners decreases the structural efficiency of the structure by adding weight. This establishes the requirement to better exploit the efficiency of bonded structures and fully understand the failure behavior of adhesively bonded composite structures, particularly when subjected to elevated loading rates due to high energy dynamic impacts (HEDI). For this reason, the NASA Advanced Composite Consortium (ACC) HEDI team developed an experimentation and numerical modeling program for high rate loading of composite joints. In the present work, the response of adhesively bonded composite joints subjected to elevated loading rates is studied numerically and validated against experimental results. Due to dynamic considerations of experiments, the idea of wedge insert was extended to use with Split Hopkinson Pressure Bar (SHPB) testing techniques. Mode-I and Mode-II test configurations were simulated to evaluate the capability of two continuum damage material (CDM) models in LS-DYNA, namely MAT162 and MAT261. Three different levels of fidelity were considered to investigate the level of detail required to numerically predict the failure behavior and the results from high fidelity analysis are presented.

MAT261↗
DOE OSTI2025

Shear Band Formation in Thin-Film Multilayer Columns Under Compressive Loading: A Mechanistic Study

Micro-pillar compression is a popular experimental technique used for characterizing the mechanical behavior of nano- and micro-laminates. The compressive stress–strain response of the column-shaped thin-film composite can be measured, and the deformation and damage features can be revealed by post-test cross-section microscopy. The development of plastic instability in the form of localized strain concentration (shear bands), leading to eventual failure, is frequently observed. In the present study, a computational approach is used to illustrate the commonality of shear band formation from a continuum standpoint. Systematic finite element analyses are conducted, showing that the strain field tends to become localized once plastic yielding commences. Distinct shear offsets of the layered structure can be revealed from the numerical model, which is similar to those observed in experiments. The actual appearance of shear bands depends on the materials’ constitutive behavior and precise geometries. Post-yield strain hardening reduces the propensity of shear band formation, while strain softening enhances it. Imperfections such as the undulated layer geometry, as well as the frictional characteristics between the specimen and test apparatus, can also influence the shear band morphology and overall stress–strain response.

finite element modeling↗
NASA NTRSDate not supplied

Application of OpenFOAM to Plume Impingement in Space Environments

After 30 years of continuous human presence in low-earth orbit, NASA is returning to the moon and eventually will go to Mars. Travelling beyond low earth orbit requires NASA to learn how humans can live in Deep Space environments – beyond the protection of Earth’s magnetosphere and at distances from Earth that prevent a quick return in case of trouble. To this end, NASA is constructing the Lunar Gateway, an ISS-like space station to be put in orbit around the moon to act as a home base for Lunar exploration for NASA astronauts. The Gateway Lunar outpost will be built incrementally, via modules which will arrive at separate times and dock to the existing structure. The incremental addition of Gateway modules, and the docking of visiting vehicles, is achieved via a sequence of firings from the approaching body’s onboard reaction control system (RCS) thrusters to achieve the required approach trajectory. The typical hypergolic chemical RCS thrusters work by firing hot gases to produce adverse thrust and the needed change in velocity to safely finish the docking process. The exhaust gas from the RCS thrusters form plumes that expand into the vacuum of space and can impinge onto the outer surfaces of the Lunar Gateway, causing unwanted forces and moments, heat loads, sediment deposition, and in extreme cases, even surface erosion - all mechanisms that can damage the Lunar Gateway and must be minimized. Both permanent and visiting modules will have this RCS thruster exhaust impingement problem. This research aims to establish existing OpenFOAM solvers as a methodology for improving simulation techniques of rocket exhaust plume impingement in space environments. The flow structure of a plume in a space environment is complex; a plume that originates from a hypergolic chemical RCS thruster and expands into a vacuum will experience several regimes of rarefication. This range includes the continuum flow in the rocket nozzle through the fully rarefied free molecular flow further from the nozzle. The flow physics is different at these two extremes, and as such, the simulation approach for plumes is generally divided into a traditional computational fluid dynamics (CFD) simulation in and near the nozzle which is coupled to a subsequent direct simulation Monte Carlo (DSMC) simulation. At this time, the scope of this research is developing, verifying, and validating a method using existing solvers in the OpenFOAM framework for performing coupled CFD/DSMC calculations to determine the extent of plume impingement loading on generic space structures. This presentation will detail code-to-code comparisons between the hyStrath dsmcFoam+ solver, developed using OpenFOAM and available as open-source, and NASA’s in-house DSMC Analysis Code (DAC). Comparisons to several open-source publication findings using DAC [3,4] are presented, and advantages of using an OpenFOAM based solver are also discussed. The presentation concludes with a discussion of future work, and a plan for coupling the dsmcFoam+ solver with CFD simulations of chemical rocket engines for unified coupled plume simulation.

DSMC↗
NASA NTRS2018

Plasticity and Damage Modeling of Stress Asymmetry and Dynamic Behavior of AFS Additive Manufactured Aluminum Alloy 2219

The Solid State Additive Manufacturing (AM) process referred as MELD that fabricated the samples in this study, provides a new path for repairing, coating, joining and additive manufacturing metals and metal matrix composites. This research will be the first application of a physics-based microstructure dependent internal state variable (ISV) plasticity and damage material model to capture the mechanical response of an AM Aluminum Alloy (AA) 2219 via the MELD process. In this research, a microstructure-based internal state variable (ISV) plasticity-damage model was used to capture the mechanical behavior of AFS 2219 aluminum alloy. Aeroprobe Corporation, creator and patent holder for the MELD process, fabricated the material by pushing a solid filler rod of AA2219-T861 material through a hollow rotating tool onto an AA2219 T851 plate substrate. As feedstock, solid or powder precursor metals are pushed through a nonconsumable rotating cylindrical tool. Herein, added layers are deposited and metallurgically bonded to substrate material or previously deposited layers by the heat generated from the rotating tool through plastic deformation of the filler material. Once a layer has been added, the tool height increases, and starts the deposition of the next layer. This process results in beneficial properties such as grain refinement, homogenization and reduced porosity (fully dense). This process will experience temperatures similar to those in the weld nugget zone (WNZ) in friction stir welding (FSW), ranging from 0.6-0.9 Tm, with Tm being the melting point of the material. MELD is highly scalable with AA deposition rates reaching over 1000 cm3/hr, which allows for MELD being used for repairs, coatings, and building components. A motivating factor driving the research for physics-based history dependent material modeling of MELD components is the ability to accurately capture the stress-state and strain rate dependence in the material caused by variations in material microstructure from the MELD processing of new or repaired components. The ISV model incorporates microstructural content and is consistent with continuum level kinematics, kinetics, and thermodynamics. These features allow the ISV model to capture large deformations at the structural scale using the kinematic and isotropic hardening, while microscale damage is obtained from the microstructural features. The benefits of the ISV model arise from the inclusion of structure-property relationships identified from microstructural characterization and experimentation. The Bauschinger effect (BE) is an important concept, vital in the accurate prediction of cyclic stress-strain response of ductile materials such as metals. The ISV model has been successfully used to capture the behavior and damage, and the BE of different aluminum alloys and steels. The ISV model uses kinematic and isotropic hardening to help capture deformations of the material at the macro scale. To understand this hardening relationship, calculating the kinematic and isotropic hardening relationship in the material is warranted for a high-fidelity model. Electron Backscattered Diffraction (EBSD) was used to characterize the as-fabricated microstructure, where a fully-dense equiaxed grain morphology with average grain size of 2.5 μm was observed. Microhardness mapping of the as-built structures, monotonic tension and compression experiments at both quasi-static (0.001/s) strain rates, tension-followed-by-compression and compression-followed-by-tension experiments were performed to obtain the set of plasticity and damage constants necessary to capture strain rate and stress state behavior of this additive material. To calibrate the plasticity-damage model, a single set of constants were determined to capture the different stress states the MELD AA2219. One set of the constants was determined from experimental true stress-strain curves for the tension and compression data. Additionally, microstructural information and data from the open literature were used as the other model constants. This research is a first of its kind for AFS AA2219, includes correlating the ISV model to the monotonic experimental results that capture the isotropic and kinematic plasticity mechanical response.

Rivera, O. G.↗
NASA NTRS2006

Multiscale Analysis of Delamination of Carbon Fiber-Epoxy Laminates with Carbon Nanotubes

A multi-scale analysis is presented to parametrically describe the Mode I delamination of a carbon fiber/epoxy laminate. In the midplane of the laminate, carbon nanotubes are included for the purposes of selectively enhancing the fracture toughness of the laminate. To analyze carbon fiber epoxy carbon nanotube laminate, the multi-scale methodology presented here links a series of parameterizations taken at various length scales ranging from the atomistic through the micromechanical to the structural level. At the atomistic scale molecular dynamics simulations are performed in conjunction with an equivalent continuum approach to develop constitutive properties for representative volume elements of the molecular structure of components of the laminate. The molecular-level constitutive results are then used in the Mori-Tanaka micromechanics to develop bulk properties for the epoxy-carbon nanotube matrix system. In order to demonstrate a possible application of this multi-scale methodology, a double cantilever beam specimen is modeled. An existing analysis is employed which uses discrete springs to model the fiber bridging affect during delamination propagation. In the absence of empirical data or a damage mechanics model describing the effect of CNTs on fracture toughness, several tractions laws are postulated, linking CNT volume fraction to fiber bridging in a DCB specimen. Results from this demonstration are presented in terms of DCB specimen load-displacement responses.

Riddick, Jaret C.↗
NASA NTRS1985

Electrical characterization of plasma-grown oxides on gallium arsenide

Plasma-grown GaAs oxides and their interfaces have been characterized by measuring the electrical properties of metal-oxide-semiconductor capacitors and of Schottky junctions. The current transport mechanism in the oxide at high electrical field was found to be Frankel-Poole emission, with an electron trap center at 0.47 eV below the conduction band of the oxide. The interface-state density, evaluated from capacitance and conductance measurements, exhibits a U-shaped interface-state continuum extending over the entire band gap. Two discrete deep states with high concentration are superimposed on this continuum at 0.40 and 0.70 eV below the conduction band. The results obtained from measurements on Schottky junctions have excluded the possibility that these two deep states originate from plasma damage. Possible origins of these states are discussed in this paper.

Hshieh, F. I.↗
NASA NTRS1990

Elevated temperature crack growth in aluminum alloys: Tensile deformation of 2618 and FVS0812 aluminum alloys

Understanding the damage tolerance of aluminum alloys at elevated temperatures is essential for safe applications of advanced materials. The objective of this project is to investigate the time dependent subcritical cracking behavior of powder metallurgy FVS0812 and ingot metallurgy 2618 aluminum alloys at elevated temperatures. The fracture mechanics approach was applied. Sidegrooved compact tension specimens were tested at 175, 250, and 316 C under constant load. Subcritical crack growth occurred in each alloy at applied stress intensity levels (K) of between about 14 and 25 MPa/m, well below K (sub IC). Measured load, crack opening displacement and displacement rate, and crack length and growth rate (da/dt) were analyzed with several continuum fracture parameters including, the C-integral, C (sub t), and K. Elevated temperature growth rate data suggest that K is a controlling parameter during time dependent cracking. For FVS0812, da/dt is highest at 175 C when rates are expressed as a function of K. While crack growth rate is not controlled by C (sub t) at 175 C, da/dt appears to better correlate with C (sub t) at higher temperatures. Creep brittle cracking at intermediate temperatures, and perhaps related to strain aging, is augmented by time dependent transient creep plasticity at higher temperatures. The C (sub t) analysis is, however, complicated by the necessity to measure small differences in the elastic crack growth and creep contributions to the crack opening displacement rate. A microstructural study indicates that 2618 and FVS0812 are likely to be creep brittle materials, consistent with the results obtained from the fracture mechanics study. Time dependent crack growth of 2618 at 175 C is characterized by mixed transgranular and intergranular fracture. Delamination along the ribbon powder particle boundaries occurs in FVS0812 at all temperatures. The fracture mode of FVS0812 changes with temperature. At 175 C, it is characterized as dimpled rupture, and at 316 C as mixed matrix superplastic rupture and matrix-dipersoid debonding.

Leng, Yang↗
NASA NTRS1991

Deformation mechanisms in negative Poisson's ratio materials - Structural aspects

Poisson's ratio in materials is governed by the following aspects of the microstructure: the presence of rotational degrees of freedom, non-affine deformation kinematics, or anisotropic structure. Several structural models are examined. The non-affine kinematics are seen to be essential for the production of negative Poisson's ratios for isotropic materials containing central force linkages of positive stiffness. Non-central forces combined with pre-load can also give rise to a negative Poisson's ratio in isotropic materials. A chiral microstructure with non-central force interaction or non-affine deformation can also exhibit a negative Poisson's ratio. Toughness and damage resistance in these materials may be affected by the Poisson's ratio itself, as well as by generalized continuum aspects associated with the microstructure.

Lakes, R.↗
DOE OSTI2022

Modeling Chemo-Mechanics with Electrolyte Infiltration to Quantify Degradation of Cathode Particles

One of the main goals in modeling lithium-ion batteries is to improve/predict longevity and resilience of new chemistries. To that end, this talk investigates the formation of stress-induced fracture within polycrystalline cathode particles and the impact on capacity loss. Physically based cathode aging dynamics is simulated in a 3D, continuum-level chemo-mechanical model. The model captures anisotropic Li diffusion within a single polycrystalline particle comprised of hundreds to thousands of randomly oriented grains. A recent addition to this model includes electrolyte infiltration, which occurs when the electrolyte seeps into surface cracks within the particle. The model predicts that particle fracture is primarily due to non-ideal grain interactions with slight dependence on high-rate charge demands. Essentially, when neighboring grains are misaligned, they expand a different rates relative to one another leading to high stresses and ultimately the formation of intraparticle cracks. The model predicts that small particles with large grains develop significantly less damage than larger particles with small grains. Finally, the model predicts most of the chemo-mechanical damage accumulates in the first charge after formation. This chemo-mechanical "damage saturation" effect indicates that initial particle fracture occurs within the first few cycles, while long-term cathode degradation is not solely chemo-mechanically induced. The principle contribution of this research is the use of an anisotropic chemo-mechanical model to test how particle geometry affect capacity fade, which predicts that particle size has a stronger effect on capacity fade than grain size and ultimately that small particles with large grains have the least capacity fade.

cathode↗
NASA NTRS1989

A study of material damping in large space structures

A constitutive model was developed for predicting damping as a function of damage in continuous fiber reinforced laminated composites. The damage model is a continuum formulation, and uses internal state variables to quantify damage and its subsequent effect on material response. The model is sensitive to the stacking sequence of the laminate. Given appropriate baseline data from unidirectional material, and damping as a function of damage in one crossply laminate, damage can be predicted as a function of damage in other crossply laminates. Agreement between theory and experiment was quite good. A micromechanics model was also developed for examining the influence of damage on damping. This model explicitly includes crack surfaces. The model provides reasonable predictions of bending stiffness as a function of damage. Damping predictions are not in agreement with the experiment. This is thought to be a result of dissipation mechanisms such as friction, which are not presently included in the analysis.

Highsmith, A. L.↗
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↗