Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Plastic flow”

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

Superior plastic flow stability of self-patterned carbide – amorphous ceramic nanostructures

Amorphous ceramics and carbides exhibit superb strength but poor plasticity. Here, we synthesized TiC-SiOC nanostructures with TiC-nanocarbides embedded in amorphous ceramic SiOC by co-sputtering followed by high-temperature annealing and/or irradiation. TiC-SiOC nanostructures exhibit high strength and good plastic flow stability even after heavy irradiation, 7 GPa at room temperature and 3.6 GPa at 700 ℃ with a uniform strain of about 10%∼18%. The uniform deformation is accommodated by the shearing of amorphous ceramic and the rotation of nanocarbides. Nanocarbides inhibit the propagation of shear banding in amorphous SiOC, and amorphous-crystal interfaces act as sinks to manage irradiation-induced defects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Rules for the crystallite size and dislocation density evolution in phases during α-ω transformation in Zr under high-pressure and severe plastic flow

The first in-situ X-ray diffraction (XRD) study of the evolution of the crystallite size and dislocation density in phases during plastic strain-induced phase transformation (PT) is performed utilizing α-ω PT in strongly pre-deformed commercially pure Zr as an important example. Rough diamond anvils (rough-DA) are introduced to intensify all occurring processes during heterogeneous compression of Zr in a diamond anvil cell (DAC). The main rule is found that during α-ω PT the crystallite size and dislocation density in ω-Zr depend solely on the volume fraction of ω-Zr and are independent of pressure, plastic strain tensor, its path, and initial nanostructure. Crystallite size in ω-Zr increases from 10 to 60 nm during the PT, while dislocation density reduces from 1.83×10 15 /m 2 to 0.6×10 15 /m 2 . Rough-DA produce a steady nanostructure in α-Zr before PT with smaller crystallite size and larger dislocation density than smooth-DA, leading to a reduction of the minimum pressure for α-ω PT to a record value 0.67 GPa, 9 times smaller than under hydrostatic loading and 5.1 times lower than the phase equilibrium pressure. In addition to strain, the kinetics of strain-induced PT unexpectedly depends on time. Also, strain-controlled part of kinetics is zero order, in contrast to the first-order kinetics with smooth-DA. The obtained results open a new window for understanding the mutual effects of nanostructure evolution and PT during severe plastic flow in various technological and natural processes. In conclusion, they may bring up economic strategies of producing nanocomposites and single-phase nanostructured materials with optimal properties.

Crystalline size↗
DOE OSTI2026

Active learning of a crystal plasticity flow rule from discrete dislocation dynamics simulations

Continuum-scale material deformation models, such as crystal plasticity (CP), can significantly enhance their predictive accuracy by incorporating input from lower-scale (i.e. mesoscale) models. The procedure to generate and extract the relevant information is however typically complex and ad hoc, involving decision and intervention by domain experts, leading to long development times. In this study, we develop a principled approach for calibration of continuum-scale models using lower scale information by representing a CP flow rule as a Gaussian process model. This representation allows for efficient parameter space exploration, guided by the uncertainty embedded in the model through a process known as Bayesian optimization (BO). We demonstrate a semi-autonomous BO loop which instantiates discrete dislocation dynamics simulations whose initial conditions are automatically chosen to optimize the uncertainty of a model CP flow rule. Our self-guided computational pipeline efficiently generated a dataset and corresponding model whose error, uncertainty, and physical feature sensitivities were validated with comparison to an independent dataset four times larger, demonstrating a valuable and efficient active learning implementation readily transferable to similar material systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Using principal component analysis to distinguish different dynamic phases in superconducting vortex matter

Vortices in type-II superconductors driven over random disorder are known to exhibit a remarkable variety of distinct nonequilibrium dynamical phases that arise owing to the competition between vortex-vortex interactions, the quenched disorder, and the drive. These include pinned states, elastic flows, plastic or disordered flows, and dynamically reordered moving crystal or moving smectic states. The plastic flow phases can be particularly difficult to characterize since the flows are strongly disordered. Here, we perform principal component analysis (PCA) on the positions and velocities of vortex matter moving over random disorder for different disorder strengths and drives. We find that PCA can distinguish the known dynamic phases as well as or better than previous measures based on transport signatures or topological defect densities. In addition, PCA recognizes distinct plastic flow regimes, a slowly changing channel flow and a moving amorphous fluid flow, that do not produce distinct signatures in the standard measurements. In conclusion, our results suggest that this position and velocity-based PCA approach could be used to characterize dynamic phases in a broader class of systems that exhibit depinning and nonequilibrium phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Experimental Validation of Exact Burst Pressure Solutions for Thick-Walled Cylindrical Pressure Vessels

Burst pressure is one of the critical strength parameters used in the design and operation of pressure vessels because it represents the maximum pressure that a vessel can withstand before failing. Historically, the Barlow formula was used as a design base for estimating burst pressure. However, it does not consider the plastic flow response for ductile steels and is applicable only to thin-walled cylinders (i.e., the diameter to thickness ratio D/t ≥ 20). A new multiaxial plastic yield theory was developed to consider the plastic flow response, and the associated theoretical (i.e., Zhu–Leis) solution of burst pressure was obtained and has gained extensive applications in the pipeline industry because it was validated by different full-scale burst test datasets for large-diameter, thin-walled pipelines in a variety of steel grades from Grade B to X120. The Zhu–Leis flow theory of plasticity was recently extended to thick-walled pressure vessels, and the associated exact flow solution of burst pressure was obtained and is applicable to both thin and thick-walled cylindrical shells. Many full-scale burst tests are available for thin-walled line pipes in the pipeline industry, but limited pressure burst tests exist for thick-walled vessels. To validate the newly developed exact solutions of burst pressure for thick-walled cylinders, this paper conducts a series of burst pressure tests on small-diameter, thick-walled pipes. In particular, six burst tests are carried out for three thick-walled pipes in Grade B carbon steel. These pipes have a nominal diameter of 2.375 inches (60.33 mm) and three nominal wall thicknesses of 0.154, 0.218, and 0.344 inches (3.91, 5.54, and 8.74 mm), leading to D/t = 15.4, 10.9, and 6.9, respectively. With the burst test data, comparisons show that the Zhu–Leis flow solution of burst pressure matches well the burst test data for thick-walled pipes. Thus, these burst tests validate the accuracy of the Zhu–Leis flow solution of burst pressure for thick-walled cylindrical vessels.

42 ENGINEERING↗
DOE OSTI2025

2024 Second Half Semi Annual Report: Modeling plasticity-mediated flow in metals with pressurized cavities

The objective is to better predict the bulk-scale mechanical behavior of porous metals that have over pressurized cavities (e.g., irradiated metals with helium bubbles) by quantifying the complex coupling among cavity aspects (e.g., size distribution, inhomogeneous overpressure values, spatial arrangement) and metal properties (e.g., rate-dependency, crystallographic lattice). This requires up-scaling local mechanical fields from the single crystal scale and will be accomplished using a homogenization approach that combines full-field numerical simulations, analytical formalisms, and physics-informed machine learning to produce symbolically-defined constitutive equations (e.g., gauge functions). These equations will satisfy the objective because they enable computationally efficient predictions that approach the accuracy of computationally expensive full-field numerical simulations, abide by theoretical requirements (e.g., conservation of energy, work conjugacy), and retain the transparency of analytical models.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling plasticity-mediated void growth at the single crystal scale: A physics-informed machine learning approach

Modeling the evolution of voids during plastic flow as well as their effects on plastic dissipation is critical for both component manufacturing and lifetime estimation purposes. To this end, we propose a rate-dependent constitutive model to homogenize the effects of semi-randomly distributed voids on single crystal plasticity whilst capturing void interaction and plastic anisotropy. Here, this present work focuses on the case of face centered cubic crystals to introduce an anisotropic gauge function applicable within the crystal plasticity formalism. The approach combines analytical methods to describe the micromechanics of the system in combination with symbolic regression to capture analytically intractable mechanisms from data. The hybrid framework uses a physics-informed genetic programming-based symbolic regression algorithm to solve a multiform optimization problem simultaneously producing a new gauge function and a new strain rate equation. This is also a multi-objective optimization problem with many competing objectives. A new search and selection step is introduced to the genetic algorithm that promotes convergence toward a global solution that better satisfies all the objectives. Overall, the symbolic equations produced leverage data-driven methods to achieve greater accuracy than comparable alternatives on an analytically intractable problem while maintaining model transparency.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microstructural, Oxidation, and Mechanical Behavior of NbTi-Based Refractory Alloys with 5 to 10 Pct Co, Cr, and Ni Additions

NbTi-based refractory alloys with additions of Co, Cr, and Ni represent an interesting medium-entropy alloy system with potential for protective oxide film formation, high strength, and ductility. This study investigates the microstructural evolution, oxidation behavior, and mechanical properties of NbTi-based alloys containing 5 to 10 at. pct Co, Cr, and Ni. CALPHAD predictions suggest that this composition range can be heat treated to obtain a predominantly body-centered cubic matrix phase. Mechanical properties, including microhardness, yield strength, maximum strength, and specific strength are evaluated through isothermal compression tests conducted between room temperature and 800 °C. The oxidation kinetics of these alloys are assessed through discontinuous oxidation tests. Parabolic oxidation kinetics were observed for NbTi–10Ni and NbTi–5Co, while linear oxidation kinetics were found for NbTi–10Cr and NbTi–10(CoCrNi). Microstructures and oxide layers are characterized using X-ray diffraction, electron backscatter diffraction, energy-dispersive X-ray spectroscopy, and scanning electron microscopy. All alloys exhibit significant mechanical softening between room temperature and 800 °C, with elastic-perfectly plastic flow observed at 800 °C. The addition of 10 pct Cr to NbTi resulted in two BCC phases up to 1050 °C, conflicting with CALPHAD predictions of a single-phase solid solution at this temperature, and resulting in higher flow stress at 800 °C. NbTi–10(CoCrNi) exhibited the lowest flow stress at 800 °C despite having more ‘cocktail effect’ potential and insufficient molar fractions of Co, Cr, or Ni to form a desirable protective oxide film.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Active microrheology and dynamic phases for pattern-forming systems with competing interactions

We consider the driven dynamics of a probe particle moving through an assembly of particles with competing long-range repulsive and short-range attractive interactions, which form crystal, stripe, labyrinth, and bubble states as the ratio of attraction to repulsion is varied. We show that the probe particle exhibits a depinning-like threshold from an elastic regime, where the probe particle is trapped by interactions with the other particles, to a plastic flow regime, where the probe particle can break bonds in the surrounding medium. For a fixed particle density, the depinning threshold and sliding velocity of the probe particle vary nonmonotonically as the attraction term is increased. A velocity minimum appears near the crystal to stripe crossover, and there is a significant increase in the depinning threshold in the bubble regime when the probe particle is strongly confined inside the bubbles. For fixed attractive interaction but increasing particle density, the behavior is also nonmonotonic, and there are jumps and drops in the velocity and depinning threshold corresponding to points at which the system transitions between different structures. There are also several distinct flow states that can be characterized by the amount of plastic deformation induced by the probe particle in the surrounding medium. Each flow state generates a different amount of effective drag on the probe particle, and there can be jumps in the velocity-force curve at transitions between the states. In conclusion, we also find that when oriented stripes are present, the probe particle can move along the stripe in an edge transport state that has a finite Hall angle.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A variational phase-field framework for thermal softening and dynamic ductile fracture

A variational phase field model for dynamic ductile fracture is presented. Here, the model is designed for elasto-viscoplastic materials subjected to rapid deformations in which the effects of heat generation and material softening are dominant. The variational framework allows for the consistent inclusion of plastic dissipation in the heat equation as well as thermal softening. It employs a coalescence function to degrade fracture energy during regimes of high plastic flow. A variationally consistent form of the Johnson–Cook model is developed for use with the framework. Results from various benchmark problems in dynamic ductile fracture are presented to demonstrate capabilities. In particular, the ability of the model to regularize shear band formation and subsequent damage evolution in two- and three-dimensional problems is demonstrated. Importantly, these phenomena are naturally captured through the underlying physics without the need for phenomenological criteria such as stability thresholds for the onset of shear band formation.

42 ENGINEERING↗
DOE OSTI2026

Towards the Ultimate Strength of Medium‐Entropy Alloys Through Pulsed Lasers

The tensile strength of metals at extreme strain rates is a key predictor of their performance in ballistic and structural impact applications. An important experimental method to reach these extreme strain rates is the use of high-amplitude, short-duration pulsed lasers. The Jupiter Laser Facility at the Lawrence Livermore National Laboratory enabled probing for the first time the mechanical response of several promising High Entropy Alloys at times on the order of nanoseconds (strain rates of ∼10 7 and ∼10 9 s −1 ). The measured strength is in the range of 6 to 10 GPa, ten times the quasistatic value. The mechanisms of plastic deformation and failure were identified and quantified through analysis and molecular dynamics simulation. The reflected wave amplitudes, obtained by VISAR, were used to determine the tensile (spall) stress. The high tensile strength obtained is due to two factors: the strain-rate dependence of plastic flow and the kinetics of void nucleation, growth, and coalescence. The experimental results are compared with an analytical prediction considering both grain-interior and grain-boundary void initiation. Molecular dynamics simulations, conducted at strain rates of 10 8 and 10 9 s −1 , rationalize the experimental results. In conclusion, they provide valuable information about the process of failure evolution, and reveal that grain boundary separation plays a pivotal role in spalling.

Materials science↗
DOE OSTI2024

Evaluating Pressure‐dependent Discharge Behavior of Foil Versus In situ Plated Lithium Metal Anodes in Solid‐State Batteries

Abstract Anode‐free manufacturing of solid‐state batteries (SSBs) shows promise to maximize energy density by eliminating excess lithium (Li) and simplifying battery production. However, high reversibility during discharge (stripping of Li) is necessary for long‐lifetime SSBs with a limited Li reservoir. Further, the plastic flow of Li changes depending on the Li thickness, leading to possible differences in discharge performance under stack pressure. This work investigates the pressure‐dependent discharge performance of anode‐free manufactured SSBs with in situ plated Li and compares the performance to that of conventional thick Li foil cells. Distinct stripping behavior is observed at low pressures (0–1 MPa), where Li diffusivity and initial interfacial contact may control accessible capacity, compared to high pressures (3–10 MPa) where mechanical deformation of Li likely governs stripping behavior. Analysis of impedance spectra collected during stripping shows that additional stack pressure delays the formation of deep, as opposed to lateral, voids in the Li anode. These results provide insights to guide the transition from thick Li foil anodes to anode‐free manufactured SSBs.

Haslam, Catherine G.↗
DOE OSTI2025

Chromium enhances the mechanical performance of 3d transition metal high entropy alloys

A comprehensive examination of the compositional effects on the deformation behavior of high entropy alloys (HEAs) was conducted through nanoindentation, indentation creep, and stress relaxation experiments at ambient temperature. The alloys investigated included NiCoFe, NiCoCr, NiCoCrFe, and NiCoCrFeMn, all of which underwent identical thermomechanical processing. From the experimental results, NiCoCr exhibited the highest maximum shear stress for dislocation nucleation (13 GPa) and nanoindentation hardness (3.8 GPa) among the alloys tested. It also showed the lowest stress sensitivity (0.02) and activation volume (3–4 b 3 ) under creep and stress relaxation conditions, indicative of its superior plastic flow properties. Density functional theory (DFT) calculations further revealed an uneven charge density and larger bonding-length variation in NiCoCr due to the presence of Cr, which increased lattice distortion, impeding dislocation movement. Furthermore, in this study, post-deformation microscopy using transmission electron microscopy (TEM) revealed a high density of stacking faults and twins in NiCoCr that effectively enhanced its creep strength. The results highlight the significance of specific elemental effects, particularly from Cr in this study, over configurational entropy effect (i.e., compositional complexity) in governing the deformation microstructure and mechanical properties of HEAs. Finally, these insights will be instrumental for the design and optimization of advanced alloys for load-bearing applications.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Efforts to stabilize composite localization elements

Localization finite elements seek to provide a robust framework for modeling ductile failure. They utilize the same constitutive model as the bulk material through the introduction of a length scale in a specialized deformation gradient that regularizes displacement discontinuity. Similar to many other elements, localization elements exhibit locking and associated pressure oscillations under incompressible plastic flow, which is a critical issue when attempting to model pressure-driven damage evolution. These issues can be drastically improved through what are essentially reduced integration techniques for the Jacobian and pressure, but there seem to be pressure-related instabilities that persist and are specific to localization elements. This memo summarizes recent efforts to mitigate and understand this problem, mostly for the 12-node composite wedge localization element in particular. At this point, it remains unclear whether the pressure fields within any localization element can be sufficiently stabilized in order to properly model failures that include softening or fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Temperature-based reactive flow model for triaminotrinitrobenzene (TATB) plastic bonded explosives

A new reactive flow model is presented for triaminotrinitrobenzene (TATB)-based plastic bonded explosives, applicable to shock initiation and steady detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and reaction products. The model incorporates temperature-dependent decomposition reaction rates. Particularly, Arrhenius model parameters were derived from quantum-based molecular dynamics simulations of TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the slow release of energy inherent to detonation in TATB-based materials. Model parameters were calibrated against gas gun shock initiation experiments and steady detonation rate stick tests. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments. The model is found to perform equally well in predicting the size-effect curve of ambient, cold, and hot rate sticks. The present work demonstrates the viability of incorporating results from subscale simulations into a continuum-scale reactive flow model.

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

Development of a Liquid Bridge Model for Particle Agglomeration and Defluidization in Plastic Pyrolysis

Molten plastic that forms during the pyrolysis of plastic or municipal solid waste feedstock can lead to particle agglomeration. A liquid bridge model is developed in MFiX-DEM and validated against a cold flow experiment with glass beads coated with silicone oil. The liquid bridge model is then extended to support an evolving liquid layer thickness for pyrolysis applications. The extended model is used to study the sensitivity of the pyrolysis reactor to solids holdup, flow conditions, and plastic properties.

Banerjee, Subhodeep↗
DOE OSTI2025

Deconvoluting thermomechanical effects in X-ray diffraction data using machine learning

X-ray diffraction is ideal for probing the sub-surface state during complex or rapid thermomechanical loading of crystalline materials. However, challenges arise as the size of diffraction volumes increases due to spatial broadening and because of the inability to deconvolute the effects of different lattice deformation mechanisms. Here, we present a novel approach that uses combinations of physics-based modeling and machine learning to deconvolve thermal and mechanical elastic strains for diffraction data analysis. The method builds on a previous effort to extract thermal strain distribution information from diffraction data. The new approach is applied to extract the evolution of the thermomechanical state during laser melting of an Inconel 625 wall specimen which produces significant residual stress upon cooling. A combination of heat transfer and fluid flow, elasto-plasticity and X-ray diffraction simulations is used to generate training data for machine-learning (Gaussian process regression, GPR) models that map diffracted intensity distributions to underlying thermomechanical strain fields. First-principles density functional theory is used to determine accurate temperature-dependent thermal expansion and elastic stiffness used for elasto-plasticity modeling. The trained GPR models are found to be capable of deconvoluting the effects of thermal and mechanical strains, in addition to providing information about underlying strain distributions, even from complex diffraction patterns with irregularly shaped peaks.

36 MATERIALS SCIENCE↗