Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale materials modeling”

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 325 records · Page 18

Thermal fluctuations in the Landau-Lifshitz-Bloch model

A formulation for thermal noise in the stochastic form of the Landau-Lifshitz-Bloch equation used for modeling the magnetization dynamics at elevated temperatures is presented here. The diffusion coefficients for thermal fluctuations are obtained via the Fokker-Plank equation using the mean-field approximation of the field in defining the free energy. The presented model leads to a mean magnetization consistent with the equilibrium magnetization for small and large particles. The distribution of the magnetization magnitude is of the Maxwell-Boltzmann type. The presented model was tested by studying the equilibrium magnetization in macrospin particles at high temperatures. The model is appealing for multiscale modeling, such as modeling heat assisted magnetic recording systems and all-optical magnetization reversal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multiscale, mechanistic modeling of irradiation-enhanced silver diffusion in TRISO particles

Tristructural isotropic (TRISO) particles are under consideration for use in several proposed advanced nuclear reactor concepts. The silicon carbide (SiC) layer in TRISO acts as a barrier to prevent the release of the fission products. However, despite remarkable retention, silver (Ag) release has been observed from intact particles, which requires investigation since the Ag isotope ( 110m Ag) has a long half-life. Previous work focused on developing a multiscale, mechanistic model for Ag diffusion accounting for temperature and microstructure effect and has been successfully validated. In this work, we expand the previous model to account for irradiation-enhanced Ag diffusivity in SiC and improve its accuracy over a wider grain size and temperature ranges relevant for advanced reactor conditions. A temperature, grain size, and flux dependent diffusivity is therefore derived using the mesoscale code MARMOT and implemented in the fuel performance code BISON. The irradiation-enhanced Ag diffusivity in SiC is compared against experimental data and validated using BISON against Ag release measurements from the Advanced Gas Reactor Fuel Development and Qualification Program (AGR-1 and AGR-2). Herein, we quantify the impact of SiC grain size, irradiation, and temperature on Ag release. In agreement with previous studies, we find accounting for SiC grain size improves agreement between BISON predictions and experimental observations for most cases. In conclusion, we also find that accounting for irradiation improves agreement for cases where Ag release was underestimated, but the impact was less significant than accounting for microstructure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Model for Humidity-Mediated Diffusion on Aluminum Surfaces and Its Role in Accelerating Atmospheric Aluminum Corrosion

Bare aluminum metal surfaces are highly reactive, which leads to the spontaneous formation of a protective oxide surface layer. Because many subsequent corrosive processes are mediated by water, the structure and dynamics of water at the oxide interface are anticipated to influence corrosion kinetics. Here, using molecular dynamics simulations with a reactive force field, we model the behavior of aqueous aluminum metal ions in water adsorbed onto aluminum oxide surfaces across a range of ion concentrations and water film thicknesses corresponding to increasing relative humidity. We find that the structure and diffusivity of both the water and the metal ions depend strongly on the humidity of the environment and the relative height within the adsorbed water film. Aqueous aluminum ion diffusion rates in water films corresponding to a typical indoor relative humidity of 30% are found to be more than 2 orders of magnitude slower than self-diffusion of water in the bulk limit. Connections between metal ion diffusivity and corrosion reaction kinetics are assessed parametrically with a reductionist model based on a 1D continuum reaction–diffusion equation. Our results highlight the importance of incorporating the properties specific to interfacial water in predictive models of aluminum corrosion.

36 MATERIALS SCIENCE↗

Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling

Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.

58 GEOSCIENCES↗

Application of NEAMS Multiphysics Framework for Species Tracking in Molten Salt Reactors

This report from Idaho National Laboratory (INL) summarizes the key modeling and simulation activities conducted under the Department of Energy (DOE) Molten Salt Reactor (MSR) Campaign during the Fiscal Year 2023 (FY23). The focus of the work was to leverage state-of-the-art modeling capabilities from the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) codes to enable novel multiphysics and multiscale modeling and simulation of MSRs. Through collaboration with NEAMS code developers, advanced multiphysics analysis capabilities for MSR systems were demonstrated by coupling depletion, thermal-hydraulics, and thermochemistry into an innovative framework for chemical species transport in MSRs. As a result, the framework can track nuclides throughout their lifetimes in the core, from production (depletion) to advection throughout the salt volume (thermal-hydraulics) and off-gassing or precipitation outside of the salt (thermochemistry). This work supports the near-term deployment of MSRs by integrating the synergistic efforts between the DOE’s MSR Campaign and NEAMS program. The resulting framework will help better connect system design modelers with experimentalists to better understand and predict complex physical behaviors in MSRs. Researchers and MSR developers alike can now leverage these new modeling and simulation capabilities to perform novel analyses with applications including: • MSR dynamics during normal operational transients and accident scenarios • Off-gas system design and performance for fuel cycle and depletion analysis • Corrosion and active chemistry control for reactor component health and lifetime determination • Source term, decay heat and activity determination in accident scenarios • Special nuclear material accountancy and chemical forensic analysis for safeguards • Digital twin development of experiments and experimental reactor demonstrations • Measurement requirements for instrumentation and control design • Uncertainty and sensitivity analysis of missing data to inform future experimental data collection.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Substitution-Mediated Calcination of Nickel-Based Cathodes: Decoupling Lithiation and Crystallization

Nickel-based layered cathodes such as LiNiO 2 offer high energy density for lithium-ion batteries, yet improvements in cycling performance and safety are required for practical use–often achieved through manganese and cobalt substitution as in LiNi 0.80 Mn 0.10 Co 0.10 O 2 (NMC811). However, how such substitution impacts calcination, the key process that governs lithiation, structural ordering, crystallization, and ultimately the resulting material properties, remains unclear. Here, in this study, we investigate substitution-mediated calcination dynamics in NMC811 compared to LiNiO 2 using multiscale-correlated in situ spectroscopy and atomistic-to-mesoscale modeling. While both systems progress through the same sequence of intermediates toward the thermodynamically favored layered phase, NMC811 exhibits an earlier onset of layering, concurrent with hydroxide decomposition followed by sluggish crystallization. Modeling reveals that Mn and Co lower the energy barrier for lithium incorporation and ordering but increase the penalty for interlayer gliding, thereby slowing crystal growth at elevated temperatures. This substitution-mediated decoupling of lithiation and crystallization explains the fine-grained microstructure observed in NMC811 versus coarsened particles in LiNiO 2 and establishes a mechanistic framework for predictive microstructure engineering of Ni-based cathodes.

Calcination↗

On the Theory and Numerical Simulation of Cohesive Crack Propagation with Application to Fiber-Reinforced Composites

The phenomenon of crack propagation is among the predominant modes of failure in many natural and engineering structures, often leading to severe loss of structural integrity and catastrophic failure. Thus, the ability to understand and a priori simulate the evolution of this failure mode has been one of the cornerstones of applied mechanics and structural engineering and is broadly referred to as "fracture mechanics." The work reported herein focuses on extending this understanding, in the context of through-thickness crack propagation in cohesive materials, through the development of a continuum-level multiscale numerical framework, which represents cracks as displacement discontinuities across a surface of zero measure. This report presents the relevant theory, mathematical framework, numerical modeling, and experimental investigations of through-thickness crack propagation in fiber-reinforced composites using the Variational Multiscale Cohesive Method (VMCM) developed by the authors.

Rudraraju, Siva Shankar↗

Microstructural heterogeneity of the buffer layer of TRISO nuclear fuel particles

Tristructural isotropic (TRISO) nuclear fuel particles contain a layered spherical shell designed to retain fission products; however, failure occurs in rare cases—commonly initiated in the porous pyrocarbon buffer layer. Achieving a comprehensive understanding of the buffer-initiated failure mechanisms requires detailed characterization of the buffer porosity and its heterogeneous distribution across multiple length scales. Here we performed FIB-SEM tomography across the buffer layer thickness (~100 µm) to produce 3D reconstructions of the buffer microstructure with 50 nm spatial resolution. We found an average overall porosity of ~14%, which does not solely account for the low density of the buffer (50% of the theoretical density). Additionally, the local porosity and its fluctuation increase from the kernel interface towards the inner pyrocarbon (IPyC) layer, which we attribute to the chemical vapor deposition process conditions during the TRISO particle fabrication. Detailed characterization of the porous microstructure—including analysis of the pore size, distribution, shape, and orientation—provides insight into the process-structure-property-performance relations of TRISO nuclear fuel particles and will inform multiscale models designed to predict the failure of TRISO particles under irradiation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Influences of shock imprinting on mix in a 3D-printed porous media

Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.

2 Photon Polymerization↗

Surface engineering to enhance heat generation and joint strength in dissimilar materials AZ31 and DP590 ultrasonic welding

A multiscale simulation approach was developed and employed to optimize the sheet surface conditions for higher interfacial temperature and joint strength in ultrasonic welding of magnesium alloy AZ31 and dual-phase steel DP590. First, a mesoscale model was used to study the relationship between friction coefficient and surface roughness, which can be modified by various engineering methods. Then a macroscopic process model was employed to study the effects of surface roughness on heat generation, indicating that a temperature increase can be achieved with rougher surfaces on two sides of both DP590 and AZ31 sheets. Samples prepared by sanding and filing, as well as grinding, were first characterized for surface roughness and then welded under ultrasonic vibration. An infrared camera was used to measure temperatures in situ for model validation. Overall, lap shear test results for the welded joint showed that the joint strength can be improved by 10~25% using filing and round grinding methods as a result of the enhanced heat generation and mechanical interlocking on the interface.

42 ENGINEERING↗

Influence of microstructure on replacement and porosity generation during experimental dolomitization of limestones

Replacement reactions commonly alter the multiscale pore structures of rocks during fluid-rock interactions. Analysis of these processes in various model fluid-rock systems during controlled laboratory experiments provides insights into the origins of microstructures found in natural materials. This study focused on understanding the effects of initial starting material permeability and resultant differences in transport pathways on porosity and mineralogical changes during limestone dolomitization.A series of replacement experiments (32–317 days in duration) have been conducted in which 1.59 cm (5/8 in.) diameter cores of two different limestones were reacted with saturated MgCl2 solutions at 200 °C. The Texas Cream (Austin Chalk) is a high-porosity, high-permeability limestone, whereas both the porosity and permeability of the Carthage Marble (Burlington Limestone) are relatively low. Altered limestones were imaged using scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM-EDX), Time-of-Flight Secondary Ion Mass Spectrometry (ToF–SIMS) and electron microprobe analysis (EMPA). A representative grain boundary of the low-porosity limestone was targeted for a focused ion beam (FIB) lift-out and characterized using transmission electron microscopy (TEM). These results were coupled with analyses of radial changes in the porosity distribution of the core derived from X-ray and neutron small- and ultra-small angle scattering ((U)SANS/(U)SAXS).The high-porosity/permeability limestone showed a four times faster bulk replacement rate than the lower-porosity/permeability material, and a different mechanism of porosity development. For the low-porosity limestone, a two-stage replacement occurred, with the reacted region of the core consisting of an inner rim in which the limestone was replaced by two calcite-dolomite solid solutions, and an outer rim in which the dolomite was replaced by magnesite. Elongated pores formed along grain boundaries at the initial limestone/dolomite reaction interface, and additional nanometer-scale porosity was formed at the secondary magnesite replacement rim. Grain boundaries were identified as preferential pathways for transport leading to dolomitization and a grain boundary diffusion rate was calculated based on microstructural characterization. In contrast, replacement in the high-porosity limestone was accompanied by porosity generation through replacement of individual grains by dolomite throughout the sample and, in longer runs, magnesite in outer parts of the core. These observations emphasize that both the mechanisms of the replacement reaction and the microstructure and chemistry of the replaced product are contingent on the initial structure of the starting material.

36 MATERIALS SCIENCE↗

Anomalous phase separation in a correlated electron system: Machine-learning–enabled large-scale kinetic Monte Carlo simulations

Significance Phase separation is crucial to the functionalities of many correlated electron materials with notable examples including colossal magnetoresistance in manganites and high- T c superconductivity in cuprates. However, the nonequilibrium phase-separation dynamics in such systems are poorly understood theoretically, partly because the required multiscale modeling is computationally very demanding. With the aid of machine-learning methods, we have achieved large-scale dynamical simulations in a representative correlated electron system. We observe an unusual relaxation process that is beyond the framework of classical phase-ordering theories. We also uncover a correlation-induced freezing behavior, which could be a generic feature of phase separation in correlated electron systems.

Zhang, Sheng↗

Coupling of Spark Plasma Sintering with Advanced Modeling to Enable Process Scale-Up: LDRD Final Project Poster

Electric-field assisted sintering (EFAS) manufacturing techniques reduce energy usage requirements, fabricate near-net-shape parts, and produce bulk nanostructured materials through rapid heating and cooling rates achievable through Joule heating. Process parameters of pulsed DC voltage and constant pressure are applied through stainless steel rams and graphite tooling to the ceramic yttrium oxide powder to be sintered. Our high-fidelity modeling code application advances manufacturing for extreme environments by enabling process-informed design and improving microstructure consistency through multiphysics multiscale process-structure-property-performance correlation simulations.

36 MATERIALS SCIENCE↗

Multiscale Modeling of Damage Processes in Aluminum Alloys: Grain-Scale Mechanisms

This paper has two goals related to the development of a physically-grounded methodology for modeling the initial stages of fatigue crack growth in an aluminum alloy. The aluminum alloy, AA 7075-T651, is susceptible to fatigue cracking that nucleates from cracked second phase iron-bearing particles. Thus, the first goal of the paper is to validate an existing framework for the prediction of the conditions under which the particles crack. The observed statistics of particle cracking (defined as incubation for this alloy) must be accurately predicted to simulate the stochastic nature of microstructurally small fatigue crack (MSFC) formation. Also, only by simulating incubation of damage in a statistically accurate manner can subsequent stages of crack growth be accurately predicted. To maintain fidelity and computational efficiency, a filtering procedure was developed to eliminate particles that were unlikely to crack. The particle filter considers the distributions of particle sizes and shapes, grain texture, and the configuration of the surrounding grains. This filter helps substantially reduce the number of particles that need to be included in the microstructural models and forms the basis of the future work on the subsequent stages of MSFC, crack nucleation and microstructurally small crack propagation. A physics-based approach to simulating fracture should ultimately begin at nanometer length scale, in which atomistic simulation is used to predict the fundamental damage mechanisms of MSFC. These mechanisms include dislocation formation and interaction, interstitial void formation, and atomic diffusion. However, atomistic simulations quickly become computationally intractable as the system size increases, especially when directly linking to the already large microstructural models. Therefore, the second goal of this paper is to propose a method that will incorporate atomistic simulation and small-scale experimental characterization into the existing multiscale framework. At the microscale, the nanoscale mechanics are represented within cohesive zones where appropriate, i.e. where the mechanics observed at the nanoscale can be represented as occurring on a plane such as at grain boundaries or slip planes at a crack front. Important advancements that are yet to be made include: 1. an increased fidelity in cohesive zone modeling; 2. a means to understand how atomistic simulation scales with time; 3. a new experimental methodology for generating empirical models for CZMs and emerging materials; and 4. a validation of simulations of the damage processes at the nano-micro scale. With ever-increasing computer power, the long-term ability to employ atomistic simulation for the prognosis of structural components will not be limited by computation power, but by our lack of knowledge in incorporating atomistic models into simulations of MSFC into a multiscale framework.

Hochhalter, J. D.↗

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗

Micromechanics-Based Structural Analysis (FEAMAC) and Multiscale Visualization within Abaqus/CAE Environment

A unified framework is presented that enables coupled multiscale analysis of composite structures and associated graphical pre- and postprocessing within the Abaqus/CAE environment. The recently developed, free, Finite Element Analysis--Micromechanics Analysis Code (FEAMAC) software couples NASA's Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC) with Abaqus/Standard and Abaqus/Explicit to perform micromechanics based FEA such that the nonlinear composite material response at each integration point is modeled at each increment by MAC/GMC. The Graphical User Interfaces (FEAMAC-Pre and FEAMAC-Post), developed through collaboration between SIMULIA Erie and the NASA Glenn Research Center, enable users to employ a new FEAMAC module within Abaqus/CAE that provides access to the composite microscale. FEA IAC-Pre is used to define and store constituent material properties, set-up and store composite repeating unit cells, and assign composite materials as sections with all data being stored within the CAE database. Likewise FEAMAC-Post enables multiscale field quantity visualization (contour plots, X-Y plots), with point and click access to the microscale i.e., fiber and matrix fields).

Arnold, Steven M.↗

Learning nonlinear operators in latent spaces for real-time predictions of complex dynamics in physical systems

Abstract Predicting complex dynamics in physical applications governed by partial differential equations in real-time is nearly impossible with traditional numerical simulations due to high computational cost. Neural operators offer a solution by approximating mappings between infinite-dimensional Banach spaces, yet their performance degrades with system size and complexity. We propose an approach for learning neural operators in latent spaces, facilitating real-time predictions for highly nonlinear and multiscale systems on high-dimensional domains. Our method utilizes the deep operator network architecture on a low-dimensional latent space to efficiently approximate underlying operators. Demonstrations on material fracture, fluid flow prediction, and climate modeling highlight superior prediction accuracy and computational efficiency compared to existing methods. Notably, our approach enables approximating large-scale atmospheric flows with millions of degrees, enhancing weather and climate forecasts. Here we show that the proposed approach enables real-time predictions that can facilitate decision-making for a wide range of applications in science and engineering.

97 MATHEMATICS AND COMPUTING↗

Multiscale modeling of solute diffusion in triblock copolymer membranes

We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model’s predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Furthermore, our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane.

36 MATERIALS SCIENCE↗