Results for “continuum mechanics of solids” — Page 4

DOE OSTI2021

Sierra/SolidMechanics 5.0 User's Guide

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra/SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra/SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes. This document describes the functionality and input syntax for Sierra/SM.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Sierra/SolidMechanics 5.2 User's Guide

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra/SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra/SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes. This document describes the functionality and input syntax for Sierra/SM.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Sierra/SolidMechanics 5.4 User's Guide

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra/SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra/SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes. This document describes the functionality and input syntax for Sierra/SM.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Sierra/SolidMechanics 5.8 User's Manual

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra/SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra/SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes. This document describes the functionality and input syntax for Sierra/SM.

42 ENGINEERING↗

DOE OSTI2022

Sierra/SolidMechanics 5.10 User's Guide

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra/SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra/SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes. This document describes the functionality and input syntax for Sierra/SM.

42 ENGINEERING↗

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 OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called â€¢ BCC_opt â€¢ FCC_opt â€¢ HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string â€œA_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters â€œA_B_Câ€, three files are provided: â€¢ A_B_C_opt.xyz: The optimized geometries in xyz format â€¢ A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. â€¢ A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. â€¢ A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Sierra/SolidMechanics 4.56 User's Guide

Siera/SolidMechanics (Sierra / SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra / SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra / SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes . This document describes the functionality and input syntax for Sierra/SM.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Sierra/SolidMechanics 4.58 User's Guide

Sierra / SolidMechanics (Sierra / SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra / SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra / SM has a versatile library of continuum and structural elements, an d a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes . This document describes the functionality and input syntax for Sierra / SM.

36 MATERIALS SCIENCE↗

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 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 OSTI2021

Quantifying the hierarchy of structural and mechanical length scales in granular systems

Continuum modeling of granular media is made possible by the existence of a length scale at and above which grain-resolved properties can be meaningfully homogenized. Progress has been made in identifying such length scales relevant to local structural properties such as porosity. However, a systematic analysis of scales above which different mechanical properties can be homogenized has yet to emerge. Here, X-ray tomography and 3D X-ray diffraction data are examined to identify such length scales. The data was obtained in-situ in compressed granular materials with rigid and flexible confinement. The experimental data are supplemented with validated discrete element simulations which examine different system sizes and different boundary conditions. Overall, our study reveals a hierarchy in the length scales of granular solids, with lengths governing structural variables being the shortest, lengths of stress variables being intermediate, and lengths of energy dissipation being the longest. All structural and mechanical length scales obey a power law based on the theory of Geostatistics, implying that the length scales can be found by analyzing samples significantly smaller than the length scales themselves. The length scales are also found to be sensitive to boundary conditions, implying that they are extrinsic features of granular media.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

A Brief Overview of Reactive Transport Codes Used in CO 2 Applications

Reactive transport modelling consists of computational and numerical models that describe the coupled physical, chemical, mechanical, and biological processes interacting with each other over a broad range of spatial and temporal scales. This report reviews some of the reactive transport codes available in the literature for carbon dioxide (CO 2 ) applications and can assist the scientific community with the applicability of each code to solve a specific problem. Approaches implemented in these codes can be categorized in three groups: (i) continuum-scale, (ii) pore-scale, and (iii) hybrid-scale approaches. Theoretical foundations, numerical implementations, and application examples using the models are described. In the continuum approach, flow and transport are formulated in terms of a representative control volume of the medium and coupled with geochemical reactions. The medium is characterized by bulk parameters such as porosity, permeability, or reactive surface area. In the pore-scale approach, each point of space in the pore network is occupied by either a fluid or solid phase. The pore-scale approach requires an exact knowledge of the spatial and temporal phase distribution. In the hybrid-scale approach, media are described by multiple characteristic length-scales, with some regions using a pore-scale approach while others used a continuum approach.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2026

Ionic Interdiffusion at Cathode|Solid-Electrolyte Interface: A Machine Learning–Assisted Multiscale Investigation and Mitigation Strategies

Future lithium batteries are expected to use solid electrolytes to achieve higher energy density and fast charge capabilities. However, most solid electrolytes are thermodynamically unstable against layered oxide cathodes. In this study, the stability of LiCoO2 (LCO) cathode with Li10GeP2S12 (LGPS) solid electrolyte is investigated using ab initio molecular dynamics (AIMD) and machine learning molecular dynamics (MLMD). The propensity of ionic interdiffusion, formation of a passivating interphase layer, and corresponding decay in cell performance is addressed using a continuum model. Large-scale MLMD simulations confirm that the LCO|LGPS interface permits interdiffusion of cobalt (Co) and other ionic species, leading to the formation and growth of a resistive interphase and to dramatic capacity fade even in the first cycle. We examine the literature evidence that incorporating a thin layer of LiNb0.5Ta0.5O3 (LNTO) between LCO and LGPS prevents the interdiffusion of ions. Atomistic simulations suggest that substituting lithium (Li) in LNTO with Co is thermodynamically unfavorable, thereby inhibiting ionic interdiffusion. The stable Nb5+/Ta5+ states form a rigid metal-oxide framework, which consequently also prevents the substitution of niobium (Nb) or tantalum (Ta). However, continuum-level analysis suggests that the higher mechanical stiffness of LNTO can lead to interfacial delamination between the LCO and LNTO. This phenomenon reduces the effectiveness of the protective layer. This paper, therefore, highlights the need to develop novel interlayers that balance low ionic interdiffusion with low mechanical stiffness.

Ncube, Musawenkosi K.↗

DOE OSTI2024

MatCal Users Guide: Release 1.3.0

Any continuum mechanics model will require three components: (1) a discretized geometry of the boundary value problem being studied, (2) the partial differential equations to be solved, and (3) the initial conditions and boundary conditions for the problem. To describe material behavior in these computational models, material models contribute to (2) the underlying equations and, occasionally, to (3) the initial conditions for the simulation. These material models can exhibit a mathematical form that is empirically based, based on first principles, or developed from both empirical observations and known physics. In general, these models are meant to represent a class of materials with well understood behavior. As a result, material models have parameters that must be tuned or calibrated so that the model response matches characterization data available for the specific material it is intended to represent when used to simulate a specific system. For simple models, such as isotropic, linear elastic materials in solid mechanics, this calibration process can be a simple analytical calculation directly extracting the parameters from experimental measurements. For complex models that have many inputs and require many characterization datasets to adequately identify the material behavior, the model calibration process can require an inverse problem approach where an optimization is performed to tune the model parameters to the available data.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

A comparative study of solute trapping in Fe-(33–45 at%) Cu alloys manufactured by laser directed energy deposition

The high cooling rates in laser directed energy deposition (DED-LB) of alloys lead to substantial amounts of solute trapping as solute atoms cannot diffuse away from the solid/liquid interface before it advances. In some concentrated alloys, this results in supersaturated solid phases that form nanoscale hierarchical microstructures when the solute atoms precipitate out during reheating from subsequent laser passes. We choose the Iron-Copper (Fe-Cu) binary alloy as model system as it is chemically homogeneous in the liquid phase and has negligible solid solubility at room temperature. Two alloys with nominal compositions in atomic (at.) %, Fe67Cu33 and Fe55Cu45, were manufactured using DED-LB. Scanning transmission electron microscopy (STEM), energy dispersive spectroscopy (EDS) and wavelength dispersive spectroscopy (WDS) were used to characterize the nanostructures and heterogeneous chemical compositions. A non-equilibrium solute partitioning model was used to compute the supersaturated chemistries of the constituent phases and validated with experimentally measured compositions. The measured phase compositions of the two alloys were very similar, at roughly 12 at.% Cu and 4 at.% Fe in the α(bcc)-Fe and ε(fcc)-Cu phases respectively, despite having different processing parameters and mechanical behavior. This indicates that the total thermal history, that depends on both the processing parameters and laser scan pattern, plays a stronger role on the final microstructure evolution than just the initial quantity of trapped solute. In addition, we find that current non-equilibrium solute partitioning models applied on the continuum scale fall short of predicting accurate quantitative phase compositions in concentrated alloys, although the qualitative trends are captured correctly.

36 MATERIALS SCIENCE↗