Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum 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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Collaborative Research: Vlasov-Maxwell Simulations to Resolve Electron Heating and Dissipation, in Quasi-Perpendicular Shocks (Final Report)

Collisionless shocks are a grand challenge problem in plasma physics and have been the subject of study for more than six decades. Shocks are ubiquitous phenomena in the universe and are responsible for transforming high energy and momentum flows into thermal energy and energetic particles. Understanding how shocks operate is of primary importance to understand the sun-earth coupling, protecting manned missions and spacecrafts from high energy particles, inertial confinement fusion, and radiation observed from astrophysical plasmas, such as supernova remnants and astrophysical jets. A major unanswered question on this frontier is, how does a collisionless plasma transform flow energy and momentum into electron thermal energy? Many potential mechanisms have been proposed to perform the conversion between flow and thermal energy in a collisionless plasma, but the answer has been elusive. To solve the mystery, non-relativistic perpendicular and quasi-perpendicular shocks are studied via a coordinated approach employing cutting-edge numerical simulations with the continuum kinetic code Gkeyll and spacecraft observations supplied by NASA's Magnetospheric Multiscale (MMS) mission. In both sets of data, we employ novel techniques to identify how and where energy is exchanged between electromagnetic fields and plasma particles. By combining cutting-edge spacecraft data and numerical simulations with these diagnostics, we have been able to identify some of the primary mechanisms responsible electron heating and dissipation in heliospheric shocks. In addition to electron dynamics, we have complete access to proton simulation and in situ data, enabling us to characterize proton heating and acceleration in quasi-perpendicular shocks. Many astrophysical quasi-perpendicular shocks have similar parameters to those in this study, including supernova remnants and shocks in galaxy clusters. Additionally, the simulation code, Gkeyll, and associated analysis scripts are openly available for community use. This award has supported a collaborative effort between Princeton University, the University of Arizona - Tucson, and the University of Maryland - College Park. The award funded three early career scientists, including the funding and training of one female postdoctoral researcher, who was involved in all aspects of the MMS data analysis and some portions of the simulation analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Computational investigation of water glasses using machine-learning potentials

The molecular origins of water’s anomalous properties have long been a subject of scientific inquiry. The liquid–liquid phase transition hypothesis, which posits the existence of distinct low-density and high-density liquid states separated by a first-order phase transition terminating at a critical point, has gained increasing experimental and computational support and offers a thermodynamically consistent framework for many of water’s anomalies. However, experimental challenges in avoiding crystallization near the postulated liquid–liquid critical point have focused attention to water’s canonical glassy states: low-density and high-density amorphous ice. Here, we use two Deep Potential machine-learning models, trained on the Strongly Constrained and Appropriately Normed density functional and the highly accurate Many-Body Polarizable potential, to conduct an investigation of water’s glassy phenomenology based on quantum mechanical calculations. Despite not being explicitly trained on amorphous ices, both models accurately capture the structure and transformation of the water glasses, including their interconversion along different thermodynamic paths. Isobaric quenching of liquid water at various pressures generates a continuum of intermediate amorphous ices and density fluctuations increase near the liquid–liquid critical pressure. The glass transition temperatures of the amorphous ices produced at different pressures exhibit two distinct branches, corresponding to low-density and high-density amorphous ice behaviors, consistent with experiment and the liquid–liquid transition hypothesis. Extrapolating transformation pressures from isothermal compressions to experimental compression rates brings our simulations into excellent agreement with data. Our findings demonstrate that machine-learning potentials trained on equilibrium phases can effectively model nonequilibrium glassy behavior and pave the way for studying long-timescale, out-of-equilibrium processes with quantum mechanical accuracy.

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

Exotic $\mathbb{Z}_N$ symmetries, duality, and fractons in 3+1-dimensional quantum field theory

Following our earlier analyses of nonstandard continuum quantum field theories, we study here gapped systems in 3+1 dimensions, which exhibit fractonic behavior. In particular, we present three dual field theory descriptions of the low-energy physics of the X-cube model. A key aspect of our constructions is the use of discontinuous fields in the continuum field theory. Spacetime is continuous, but the fields are not.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Exotic symmetries, duality, and fractons in 2+1-dimensional quantum field theory

We discuss nonstandard continuum quantum field theories in 2+1 dimensions. They exhibit exotic global symmetries, a subtle spectrum of charged excitations, and dualities similar to dualities of systems in 1+1 dimensions. These continuum models represent the low-energy limits of certain known lattice systems. One key aspect of these continuum field theories is the important role played by discontinuous field configurations. In two companion papers, we will present 3+1-dimensional versions of these systems. In particular, we will discuss continuum quantum field theories of some models of fractons.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Ion Size Effects on the Thermodynamic, Kinetic, and Mechanical Properties during Ion Exchange in Solid-State Electrolytes

Ion exchange offers a pathway to impose residual compressive stresses to mitigate the electro-chemo-mechanical cracking of solid-state electrolytes such as lithium lanthanum zirconium oxide. This study uses a coupled multiscale framework (integrating density functional theory (DFT), molecular dynamics (MD), and continuum modeling) to examine how exchange ion size influences stress, diffusion, fracture toughness, and electronic properties. Larger isovalent ions (Na + , Ag + , K + ) were exchanged with Li + , with DFT confirming their preference for octahedral 96h sites and a linear relationship between ion size and chemical free expansion coefficient. MD simulations reveal stress and concentration effects on exchange ion diffusivity at elevated temperatures, with Na + and Ag + maintaining favorable mobility while K + showing concentration-dependent clustering. Continuum modeling predicts the range of fracture strength improvements and the required ion exchange concentration profile. It was shown that a 5% surface exchange concentration can induce ∼0.6 GPa of surface compressive stress using Na + and ∼1.0 GPa of surface compressive stress using Ag + . On the other hand, larger ion exchange species may penalize Li + diffusivity by increasing the activation volume and activation energy. Interestingly, Na + has a negligible penalty on Li-ion diffusivity. The room temperature Li + ion diffusivity is reduced by ∼40% with Ag + ion exchange. Electronic band structure analysis shows no size-dependent change in the bandgap, though Ag + introduces localized defect states near the valence band maximum. This study highlights ion size as a key factor in optimizing LLZO properties, offering a framework to improve the solid-state battery performance.

Jagad, Harsh D. [Brown Univ., Providence, RI (Unit↗
DOE OSTI2024

Accelerating multiscale electronic stopping power predictions with time-dependent density functional theory and machine learning

Knowing the rate at which particle radiation releases energy in a material, the “stopping power,” is key to designing nuclear reactors, medical treatments, semiconductor and quantum materials, and many other technologies. While the nuclear contribution to stopping power, i.e., elastic scattering between atoms, is well understood in the literature, the route for gathering data on the electronic contribution has for decades remained costly and reliant on many simplifying assumptions, including that materials are isotropic. We establish a method that combines time-dependent density functional theory (TDDFT) and machine learning to reduce the time to assess new materials to hours on a supercomputer and provide valuable data on how atomic details influence electronic stopping. Our approach uses TDDFT to compute the electronic stopping from first principles in several directions and then machine learning to interpolate to other directions at a cost of 10 million times fewer core-hours. We demonstrate the combined approach in a study of proton irradiation in aluminum and employ it to predict how the depth of maximum energy deposition, the “Bragg Peak,” varies depending on the incident angle—a quantity otherwise inaccessible to modelers and far outside the scales of quantum mechanical simulations. The lack of any experimental information requirement makes our method applicable to most materials, and its speed makes it a prime candidate for enabling quantum-to-continuum models of radiation damage. The prospect of reusing valuable TDDFT data for training the model makes our approach appealing for applications in the age of materials data science.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Consistent linearization of micromorphic continuum theories

We outline the procedure of consistent linearization and apply it to the micromorphic, microstretch, and micropolar theories of continua. This yields tractable linear theories for nonlinearly elastic microstructured materials undergoing finite deformations. The results may be readily utilized in computational mechanics, stability and bifurcation analyses, and small-deformation problems in the context of these types of continua. Our results generalize those existing in current literature and facilitate their recovery upon incorporating appropriate kinematic and constitutive assumptions.

McAvoy, Ryan C.↗
DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗
DOE OSTI2023

A Three-Dimensional, Thermodynamically and Variationally Consistent, Fully Coupled, Electro-Chemo-Thermo-Mechanical Model of Solid-State Batteries

Here we present a theoretical framework for high fidelity modeling of three-dimensional electro-chemo-thermo-mechanical coupled systems such as solid-state batteries. The framework provides a useful means of utilizing molecular level information to predict the overall performance of the system at the continuum level. The central statement of the theoretical framework is a thermodynamically consistent inf-sup problem. The total potential of the coupled system consists of the Helmholtz free energy, the electrical kinetic potential, the chemical potential, the Fourier potential, the chemical reaction potential, and the external power expenditure. All governing equations derive from variational principles. We also show that the resulting initial-boundary-value problem (IBVP) is consistent with existing, widely adopted equations for modeling batteries. Moreover, the IBVP can be naturally discretized using the (variational) finite element method to simulate practical solid-state batteries with arbitrary geometries. Several examples illustrate multiphysics coupling between pairs of the physical field equations, culminating in a comprehensive three-dimensional model of charge-discharge in a solid state battery.

25 ENERGY STORAGE↗
DOE OSTI2024

Kinetics of HCP-BCC Phase Transition Boundary in Magnesium at High Pressure

Under high pressures, many crystalline metals undergo solid–solid phase transformations. In order to accurately model the behavior of materials under extreme loading conditions, it is essential to understand the kinetics of phase transition. Using molecular dynamics simulations, this work demonstrates the feasibility of characterizing the speeds of a moving phase boundary using atomistic simulations employing a suitable empirical potential for single-crystal magnesium. The model can provide temperature- and tensorial stress-dependent velocity of a moving phase boundary as a rate-limiting contribution to the kinetics of phase transformation in continuum codes. Results demonstrate that a nonlinear interaction exists between plasticity and phase transition, facilitating a jump in the velocity of a moving phase boundary, facilitated by activated plastic deformation mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Cristae formation is a mechanical buckling event controlled by the inner mitochondrial membrane lipidome

Cristae are high–curvature structures in the inner mitochondrial membrane (IMM) that are crucial for ATP production. While cristae–shaping proteins have been defined, analogous lipid–based mechanisms have yet to be elucidated. Here, we combine experimental lipidome dissection with multi–scale modeling to investigate how lipid interactions dictate IMM morphology and ATP generation. When modulating phospholipid (PL) saturation in engineered yeast strains, we observed a surprisingly abrupt breakpoint in IMM topology driven by a continuous loss of ATP synthase organization at cristae ridges. We found that cardiolipin (CL) specifically buffers the inner mitochondrial membrane against curvature loss, an effect that is independent of ATP synthase dimerization. To explain this interaction, we developed a continuum model for cristae tubule formation that integrates both lipid and protein–mediated curvatures. This model highlighted a snapthrough instability, which drives IMM collapse upon small changes in membrane properties. We also showed that cardiolipin is essential in low–oxygen conditions that promote PL saturation. These results demonstrate that the mechanical function of cardiolipin is dependent on the surrounding lipid and protein components of the IMM.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A novel peridynamics-based approach to predict pharmaceutical tablet robustness

The pharmaceutical drug product development process can be greatly accelerated through the use of modeling and simulation techniques to predict the manufacturability and performance of a given formulation. The anticipation and possible mitigation of tablet damage due to manufacturing stresses represents a specific area of interest in the pharmaceutical industry for predicting formulation and tableting performance. While the finite element method (FEM) has been extensively used for predicting the mechanical behavior of powder material in the compaction processes, a shortcoming of the approach is the inherent difficulty to predict discontinuities (e.g., damage or cracking) within a tablet as FEM is a continuum-based approach. In this work, we propose a novel method utilizing peridynamics (PD), a numerical method that can capture discontinuities such as tablet fracture, to predict the evolution of damage and breakage in pharmaceutical tablets. The approach links (1) the finite element method – to elucidate the behavior of powders during die compaction – with (2) the peridynamics modeling technique – to model the discontinuous nature of damage and predict tablet breakage during the critical stages of unloading and ejection from the compression die. This short communication presents a proof of concept including a workflow to calibrate the linked FEM-PD simulation models. Further, it demonstrates promising results from a preliminary experimental validation of the approach. Following further development, this approach could be used to guide the optimization of compression processes through targeted changes to formulation material properties, compression process conditions, and/or tooling geometries to deliver improved process efficiency and tablet robustness.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Scale Invariance of Hot Spot Formation in TATB High Explosives

Shock-induced detonation of insensitive high explosives based on 1,3,5-triamino-2,4,6-trinitrobenzene starts with formation of hot spots at microstructural defects but has eluded atomistic modeling treatment at micron length scales. To this end, we performed multimicron scale all-atom molecular dynamics (MD) simulations of hot spots that form during the collapse of cylindrical pores with diameters between 10 and 300 nm. Our MD simulations show that hot spots formed at pores larger than 20 nm exhibit temperature fields with scale-invariant features for sizes up to at least 300 nm. Through a continuum-based grain-scale modeling framework, we span and extend beyond the size scales currently accessible to MD and find that hot spot scale invariance is a general feature that arises when the mechanical strength is insensitive to strain rate. Finally, our results demonstrate the applicability of all-atom MD to simulate the complicated dynamical evolution of micron-sized systems and bolster confidence in insights from MD simulations of materials that exhibit strength with negligible rate dependence over the relevant intervals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stress-Dependent Chemo-Mechanical Performance of Amorphous Si Anodes for Li-Ion Batteries upon Lithiation

Alloying-type anodes are significantly governed by their chemo-mechanical performance during the electrochemical cycling. The reaction-induced huge volumetric change of these anodes may cause material degradation and failure under mechanical constraints. Here, we investigate the stress-dependent lithiation behavior of amorphous Si (a-Si) anodes using molecular dynamics simulations. It is indicated that a-Si anodes can sustain higher hydrostatic stress than biaxial/uniaxial ones without the occurrence of mechanical failure. Thermodynamic and electrochemical calculations demonstrate that although the lithiation procedure also affects the thermodynamic stability of a-Si anodes, it is mainly dominated by the external mean stresses. Compressive stress is confirmed to destabilize a-Si anodes and further trigger their capacity fading. Compared with our atomistic simulations, previous continuum models underestimate the open-cell potentials of a-Si anodes, due to their ignored large volumetric deformation at higher stresses and Li concentrations. Finally, this computational study provides the intensive atomic-level understanding of the stress-dependent lithiation behavior of a-Si anodes.

25 ENERGY STORAGE↗
DOE OSTI2025

Dislocation nano-hydrides in nickel: Nucleation, evolution and effects on dislocation behaviors

Nano-hydrides have been predicted to precipitate at the core of edge dislocations in the Ni-H system, a mechanism that may promote hydrogen embrittlement. However, nano-hydride nucleation, growth, and effects on dislocation behavior have seldom been explored. This work combines molecular dynamics grand canonical Monte Carlo (MD-GCMC) simulations and continuum modeling to uncover a wide range of phenomena linked to dislocation nano-hydrides. Simulations reveal that nano-hydrides can be stabilized at dislocation cores with all character angles, including screw segments, due to the hydrostatic stresses around the cores of the Shockley partials. Nano-hydride nucleation takes place in these regions, and growth is dictated by the character angle θ of the perfect dislocation. The equilibrium stacking fault width d eq varies dynamically to increase the local hydrostatic stress field and facilitate the formation of the nano-hydride, forming a constriction-like feature and leading to three distinct behaviors: d eq decreases for θ > 30°, d eq remains unchanged for θ = 30°, and d eq increases for θ < 30°. Remote hydrostatic and Escaig stresses are also shown to influence the nucleation stage, implying stress concentrations such as those ahead of crack tips may facilitate nano-hydride precipitation. Moreover, we identify a new hydrogen-induced 60° dislocation reaction that emits a Shockley partial on a conjugate plane, with potential implications for twin nucleation. Testable predictions from this study are then used to reinterpret previous results from the literature. In conclusion, these findings provide a comprehensive framework to assess nano-hydride formation and evolution at dislocations in nickel and other face-centered cubic metals, with important implications to hydrogen embrittlement.

Dislocations↗