Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale simulations”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

GrainPaint: A multi-scale diffusion-based generative model for microstructure reconstruction of large-scale objects

Simulation-based approaches to microstructure generation can suffer from a variety of limitations, such as high memory usage, long computational times, and difficulties in generating complex geometries. Generative machine learning models present a way around these issues, but they have previously been limited by the fixed size of their generation area. Here, we present a new microstructure generation methodology leveraging advances in inpainting using denoising diffusion models to overcome this generation area limitation. We show that microstructures generated with the presented methodology are statistically similar to grain structures generated with a kinetic Monte Carlo simulator, SPPARKS.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Bridging Scales to Model Reactive Diffusive Transport in Porous Media

Published by ECS. Two novel scale-bridging algorithms to model reaction-diffusion transport in porous media are presented. The algorithms are based on direct numerical simulations and couple the information of a micro-scale model, which accounts for the large field of view provided by micro X-ray computed tomography (X-ray CT), and a nano-scale model, which locally resolves transport in the fine structure extracted from nano X-ray CT. The micro-scale model is discretized in the through-plane direction into a 1D grid, where effective properties and internal boundaries are determined based on the results from the nano-scale model. The validated algorithms are used to examine transport of oxygen in precious group metal-free electrodes considering both zero- and first-order kinetics. Unlike conventional methods, the results show that the effective diffusivity is not a passive property but increases in regions where the reaction-rate coefficient is large. The proposed algorithms account for the multiscale coupling of reaction-diffusion transport and material microstructure, thus improving the predictions compared to conventional methods.

42 ENGINEERING↗
DOE OSTI2026

Revealing EDL-driven reduction mechanisms in binary, ternary, and quaternary fluorinated electrolytes via an integrated MD–DFT–ML framework

Accurately predicting solid electrolyte interphase (SEI) formation requires explicitly resolving the electric double layer (EDL) structure, which deviates significantly from that of the bulk electrolyte. Although an established molecular dynamics (MD) and Density Functional Theory (DFT) framework can model SEI formation by evaluating reduction reactions of local clusters in the EDL, it suffers from a combinatorial computational bottleneck. To overcome this limitation, we introduce a machine-learning-accelerated simulation workflow (MD–DFT–ML), integrating a gradient-boosted regression model trained on EDL composition data to efficiently predict reduction potentials. We apply this framework to seven fluorinated electrolytes comprising fluorinated anions, a fluorinated ester solvent, two types of diluent (ion-solvating ester vs. non-solvating ether), and an FEC additive. The analysis shows that the EDL selectively accumulates cation-binding species; consequently, the non–cation-binding ether diluent rarely enters the EDL and makes minimal contributions to SEI formation. DFT calculations on statistically representative EDL clusters provide reduction potentials and fluorine-release pathways, while the ML model, which substantially reduces the DFT workload, predicts cluster reduction energies with a mean absolute error of 0.1 eV. The combined MD–DFT–ML approach also quantifies contributions from different sources to LiF formation in the SEI. This methodology establishes a generalizable route for multiscale modeling electrolyte and interphase design for next-generation electrochemical energy-storage systems.

DFT-MD-ML workflow↗
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 OSTI2023

Design, fabrication, and testing of a water-cooled tungsten-foam high-heat flux plasma-facing module

Here, we present a comprehensive study detailing the design, fabrication and testing a new modular “finger” component, featuring a hexagonal tungsten foam armor tile bonded to the top surface of the water-cooled finger unit. The mechanical, thermal, and particle transport properties of the tungsten foam armor were determined for a wide range of porosity and cell size, and then combined in an equivalent solid material with “effective” properties suitable for large-scale simulations. A series of 3-D multiphysics simulations for coupled thermal and mechanical behavior of the “finger” module are used to optimize the design before fabrication. The simulation results show that, for the expected experimental conditions, a safe range of heat flux and foam porosity is possible for both water-cooling and helium-cooling. Based on these multiphysics simulations, an optimum finger cooling module has been fabricated at Ultramet, Inc. A series of experiments were designed to study the thermomechanical performance of the “finger” unit under steady-state and transient thermal loading at surface heat flux up to 16 MW/m 2 in the High Energy Flux Test facility (HEFTY) at UCLA. Testing conditions included low-power, steady-state plasma exposure at 6 MW/m 2 , followed by severe high heat flux cycles at 16 MW/m 2 . Results of these tests are explained through comparisons with multiphysics simulations. To eliminate the effects of thermomechanical damage carried out on a module on its thermo-mechanical response on a subsequent test, three different yet identically designed and fabricated finger modules were tested. The results show that the current design is capable of cooling the walls of fusion devices at surface heat flux up to 6 MW/m 2 without failure, and that the failure heat flux limit is below 16 MW/m 2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Group-Additivity–Embedded Multiscale Modeling for Electric Field-Enhanced Nanocatalysis

Elucidating structure-performance relationships remains a central challenge in field-enhanced catalysis, where nanoparticles exhibit nonuniform surface sites with site-dependent responses to electric fields. Low-coordination sites (edges, corners, and tips) are particularly electric field-sensitive (EF), leading to nonuniform charge distribution, adsorption energies, and catalytic activity. Here, using ammonia decomposition on a ruthenium cluster as a model system, we develop a transferable multiscale framework integrating density functional theory, group additivity (GA), Brønsted-Evans-Polanyi scaling, and microkinetic modeling to predict EF-dependent activity across nonuniform cluster sites. Across sites and fields, the nitrogen adsorption energy (E N ) emerges as the governing descriptor, yielding robust volcano relationships whose optimum shifts systematically with field: negative fields strengthen N binding via electron accumulation, while positive fields weaken N binding via charge depletion, moving the optimal E N toward weaker binding. Microkinetic analysis shows that N≡N bond formation remains the key kinetic bottleneck over most conditions; positive fields lower the effective barrier and, critically, increase the fraction of near-optimal active sites, leading to a net enhancement in overall activity relative to zero-field and negative-field cases. By capturing EF- and site-dependent energetics with high accuracy and low computational cost, this GA-embedded multi-scale simulation workflow provides a physically interpretable route to predict and design field-enhanced nanocatalysis.

ammonia decomposition↗
DOE OSTI2021

Damage Modeling of Solid Oxide Fuel Cells Accounting for Redox Effects

This work applies a multiscale mechanistic damage model developed for brittle ceramics and implemented in commercial finite element (FE) packages via user subroutines to study progressive damage in solid oxide fuel cells (SOFC) subjected to thermomechanical loading under normal operating and shutdown conditions including redox effects. The damage model captures the micromechanics of stiffness reduction due to material porosity change and microcracking and integrates the as-obtained stiffness reduction law into a continuum damage mechanics (CDM) formulation for the evolution of microcracks up to fracture. The volumetric “swelling” that occurs during redox is treated in constitutive modeling similarly to thermal expansion, but swelling strains are irreversible. Furthermore, this damage model was first validated through predictions of strength and stress-strain response for the SOFC electrode materials. Next, it has been applied to predict the potential for degradation in a generic planar SOFC stack with large active area cells. Multicell stack models were simulated in both co-flow and counter-flow configurations. In addition, a constant temperature redox cycle was also simulated to capture overall cell electrode damage due to volumetric swelling of the nickel (Ni)-based anode in the anode-supported cells.

25 ENERGY STORAGE↗
DOE OSTI2022

Coupled Monte Carlo and thermal-hydraulics modeling of a prismatic gas reactor fuel assembly using Cardinal

Cardinal is a MOOSE application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics to the MOOSE framework, closing the neutronics and thermal-fluid gaps in conducting tightly-coupled, high-resolution multiscale and multiphysics analyses. By leveraging MOOSE's interfaces for wrapping external codes, Cardinal overcomes many challenges encountered in earlier multiphysics coupling works, such as file-based I/O or overly-restrictive geometry mapping requirements. In this work, we leverage a subset of the multiphysics interfaces in Cardinal to perform coupling of OpenMC neutron transport, MOOSE heat conduction, and THM thermal-fluids for steady-state modeling of a prismatic gas reactor fuel assembly.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Digital Twins for Materials

Digital twins are emerging as powerful tools for supporting innovation as well as optimizing the in-service performance of a broad range of complex physical machines, devices, and components. A digital twin is generally designed to provide accurate in-silico representation of the form (i.e., appearance) and the functional response of a specified (unique) physical twin. This paper offers a new perspective on how the emerging concept of digital twins could be applied to accelerate materials innovation efforts. Specifically, it is argued that the material itself can be considered as a highly complex multiscale physical system whose form (i.e., details of the material structure over a hierarchy of material length) and function (i.e., response to external stimuli typically characterized through suitably defined material properties) can be captured suitably in a digital twin. Accordingly, the digital twin can represent the evolution of structure, process, and performance of the material over time, with regard to both process history and in-service environment. This paper establishes the foundational concepts and frameworks needed to formulate and continuously update both the form and function of the digital twin of a selected material physical twin. The form of the proposed material digital twin can be captured effectively using the broadly applicable framework of n-point spatial correlations, while its function at the different length scales can be captured using homogenization and localization process-structure-property surrogate models calibrated to collections of available experimental and physics-based simulation data.

36 MATERIALS SCIENCE↗
DOE OSTI2024

3.0 - MOOSE: Enabling massively parallel multiphysics simulations

The development of MOOSE has kept accelerating since the last release, with over 2,100 pull requests merged over the last 30 months that involved nearly fifty contributors across close to a dozen institutions internationally. The growth in MOOSE's capabilities and downstream applications is reflected in the growth of the community. User support provided on the GitHub discussions forum has steadily increased to nearly 50 daily interactions. New simulation projects, notably to model advanced nuclear reactor and fusion devices, are driving a significant expansion of the capabilities. This paper reports on these developments, with several major released features, new physics modules, and key improvements to the user experience and simulation workflow.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Bridging Scales in Bioenergy and Catalysis: A Review of Mesoscale Modeling Applications, Methods, and Future Directions

Between the molecular and reactor scales, which are familiar to the chemical engineering community, lies an intermediate regime, here termed the “mesoscale,” where transport phenomena and reaction kinetics compete on similar time scales. Bioenergy and catalytic processes offer particularly important examples of mesoscale phenomena owing to their multiphase nature and the complex, highly variable porosity characteristic of biomass and many structured catalysts. In this review, we overview applications and methods central to mesoscale modeling as they apply to reaction engineering of biomass conversion and catalytic processing. A brief historical perspective is offered to put recent advances in context. Applications of mesoscale modeling are described, and several specific examples from biomass pyrolysis and catalytic upgrading of bioderived intermediates are highlighted. Methods including reduced order modeling, finite element and finite volume approaches, geometry construction and import, and visualization of simulation results are described; in each category, recent advances, current limitations, and areas for future development are presented. Owing to improved access to high-performance computational resources, advances in algorithm development, and sustained interest in reaction engineering to sustainably meet societal needs, we conclude that a significant upsurge in mesoscale modeling capabilities is on the horizon that will accelerate design, deployment, and optimization of new bioenergy and catalytic technologies.

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

Perspective—Lithium Metal Nucleation and Growth on Conductive Substrates: A Multi-Scale Understanding from Atomistic, Nano-, Meso-, to Micro-Scales

Reversible lithium (Li) deposition and stripping on conductive substrates like copper is vital for long-life high-energy rechargeable Li batteries. The reversibility is believed hindered by solid electrolyte interphase (SEI) formation and Li dendritic growth. Via in situ, operando cryo-microscopy and molecular dynamics simulation, we discovered amorphous Li before a disorder-order phase transition to bcc phase. The kinetics plays a significant role in Li nucleation and morphology evolution. This perspective on Li nucleation and growth from atomistic to nano- (<20 nm), meso- (20–100 nm), and micro-scales (>100 nm) provides a practical guidance on regulating dense Li deposits reversibly for long life performance.

25 ENERGY STORAGE↗
DOE OSTI2025

AK112: Full Waveform Inversion Tomography of Alaska Improves Waveform Fits While Imaging Crustal, Mantle, and Slab Structure

We report a full waveform inversion tomography model of Alaska and the surrounding regions, inferring radially anisotropic shear and isotropic compressional wavespeeds by fitting complete waveforms from 120 regional earthquakes. Our multiscale approach inverted time–frequency phase misfits (maximum period of 100 s), starting with a minimum period of 40 s and ending at 20 s in 7 stages and 112 total iterations. The model (AK112) was evaluated by computing the misfits for 36 independent validation events. We find that misfit reductions were large and equal (∼55%) for both the inversion and validation data sets, providing confidence in the model. AK112 also provides much better waveform fits compared to other reported models for the region, including an isotropic version of itself, highlighting the importance of anisotropy. The model resolves known crustal, upper mantle, and slab structure to depths of 100 km with new detail: sedimentary basins in the Alaskan Shelf, Cook Inlet, and Colville basins, among others; discontinuous lithospheric structure across major terrane boundaries; and subducting slab geometry and back‐arc volcanic sources. In addition to tectonic interpretations, the model enables full waveform simulations for long‐period earthquake ground motions and source characterization (e.g., moment tensor and finite‐fault inversion).

Rodgers, Arthur [Lawrence Livermore National Labor↗
DOE OSTI2024

A spatially-resolved model of neutron-irradiated tungsten coupling stochastic cluster dynamics and finite deformation plasticity

Structural materials used in nuclear reactors face severe degradation in mechanical properties, such as hardening and embrittlement. At the microscopic scale, this occurs due to creation and accumulation of irradiation-induced defects and their interaction with system dislocations. Although techniques exist which can model evolution of irradiation defects, for instance kinetic transport theory-based models, their interaction with mechanical deformation of the bulk material has not been investigated extensively. In this work, we demonstrate a novel spatially-resolved multiscale coupling between microscopic irradiation defect evolution, modeled using Stochastic Cluster Dynamics (SCD) and macroscopic mechanical deformation modeled using a finite-deformation plasticity model. SCD is used to determine the statistically averaged defect cluster spacing, dependent on operating conditions such as irradiation dose and temperature. This acts as an initial condition that governs the critical resolved shear stress of dislocation glide in the macroscopic plasticity model. This framework is used to predict mechanical behavior in post-mortem test of irradiated Tungsten samples, which has found its importance as structural material used in nuclear reactors. The results obtained using the coupled approach are in good agreement with experimental data of uniaxial tension tests. The model is able to capture the effect of temperature and irradiation dose on the material hardening. Two methods are proposed to estimate hardness – using Tabor's Law relating uniaxial yield stress to hardness and from flat-punch simulations. The results are in reasonable agreement with hardness data from micro-indentation experiments of irradiated Tungsten samples. Finally, the model is also able to reveal microstructural details such as spatial variation in defect density and local stress.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Sensitivity Study of Multiscale and Phenomenological Elasto-Viscoplastic Grade 91 Material Models for Component-Scale Response

Many advanced nuclear reactor concepts currently being developed are targeting higher operating temperatures relative to the current fleet of light water nuclear reactors, for efficiency gains and other operational considerations. The design of high temperature structural components with reliable long-term operational performance will depend on material models that accurately capture the inelastic deformation mechanisms active in these environments. In this work, we perform a detailed parameter sensitivity analysis of two unified elasto-viscoplastic Grade 91 material models capable of capturing long term high temperature creep deformation. The first model is a phenomelogical material model from the Nuclear Engineering Material Library (NEML) developed at Argonne National Lab. The NEML model parameters and their uncertainty were fit to a range of Grade 91 experimental data using Bayesian Markov Chain Monte Carlo analysis. The second model is a LAROMance data-driven surrogate material model developed at Los Alamos National Lab. The LAROMance model is fit to a large database of responses produced by a mechanistic crystal plasticity based polycrystal model. Parameters for the LAROMance surrogate material model reflect the pedigree of the Grade 91 microstructure. Both material models have been integrated into the Grizzly code, based on the open-source MOOSE multiphysics simulation framework, to simulate both the progression of aging mechanisms and the effects of that aging on nuclear power plant structures. Grizzly is used analyze a three-dimensional Grade 91 piping system to compare the long-term inelastic response predicted by these two fundamentally different models and assess the sensitivity of the material model input parameters on this quantity of interest.

42 ENGINEERING↗
DOE OSTI2023

A Hybrid Domain Overlapping Method for Coupling System Thermal Hydraulics and CFD Codes

A hybrid multiscale coupling methodology based on a domain overlapping approach has been developed for coupling System Thermal Hydraulics (STH) and Computational Fluid Dynamics (CFD) codes. The method has been implemented between the modern STH code SAM and the CFD code NekRS, using the coupling tool Cardinal. The coupling aims to extend the STH code's applicability to scenarios where local momentum and energy transfers are important yet difficult for STH codes to capture, such as three-dimensional mixing. Two coupling strategies are implemented and compared: a hybrid domain overlapping method and the conventional domain decomposition method. Here, the strategies are applied to two closed-loop applications, and the present method shows superior stability behavior when compared to the domain decomposition method. Then, the present coupling method is validated against experimental data from a double T-junction experiment. The present STH/CFD coupling shows improved agreement with experimental data when compared to STH standalone simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Improving the particle dry deposition scheme in the CMAQ photochemical modeling system

Dry deposition of atmospheric aerosols in large-scale models is a critical, but highly uncertain, sink process with a strong dependence on particle size, meteorological conditions, and land surface properties. This study investigates the particle dry deposition scheme implemented in the standard Community Multiscale Air Quality (CMAQ) model v5.2.1, characterizes its underlying parameterized components with comparison to a similar scheme in a contemporary regional-scale model, and proposes two updated schemes that are then evaluated with available ambient particle deposition velocity (V d ) measurements. Both updated schemes reduce the surprisingly strong dependence of deposition velocity on the aerosol mode width, with one scheme further introducing a dependence on vegetation coverage that is broadly consistent with variability in observations between vegetated and non-vegetated surfaces. Compared to the base scheme, the updated scheme with vegetation dependence increases V d for submicron particles and decreases it for larger particles by an average of 37% and –66%, respectively. This scheme performs statistically better than the base scheme, reducing fractional biases by 56%–97% for vegetated land-use types and has roughly equivalent performance over water. The base and updated schemes are tested with three annual CMAQ (v5.2.1) simulations for the year 2011; predicted ambient aerosol concentrations are evaluated with routine monitoring network observations and predicted dry deposition fluxes are evaluated with data from the Clean Air Status and Trends Network (CASTNET). The updated scheme with vegetation dependence reduces negative fractional biases for PM 10 by 41% and positive fractional biases for PM 2.5 organic carbon by 15%. This scheme has been incorporated into the most recent publicly accessible versions of CMAQ (v5.3 and beyond) to replace the scheme used in previous versions of CMAQ (v4.5 through v5.2.1).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A High‐Resolution Tropical Mesoscale Convective System Reanalysis (TMeCSR)

Abstract Modern global reanalysis products have greatly accelerated meteorological research in synoptic‐to‐planetary‐scale phenomena. However, their use in studying tropical mesoscale convective systems (MCSs) and their regional‐to‐global impact has mostly been limited to supplying initial and boundary conditions for MCS‐resolving simulations and providing information about the large‐scale environments of MCSs. These limitations are due to difficulties in resolving tropical MCS dynamics in the relatively low‐resolution global models and that tropical MCSs often occur over poorly observed regions. In this work, a Tropical MCS‐resolving Reanalysis product (TMeCSR) was created over a region with frequent tropical MCSs. This region spans the tropical Indian Ocean, tropical continental Asia, Maritime Continent, and Western Pacific. TMeCSR is produced by assimilating all‐sky infrared radiances from geostationary satellites and other conventional observations into an MCS‐resolving regional model using the Ensemble Kalman Filter. The resulting observation‐constrained high‐resolution (9‐km grid spacing) data set is available hourly during the boreal summer (June‐August) of 2017, during which widespread severe flooding occurred. Comparisons of TMeCSR and European Center for Medium Range Weather Forecast Reanalysis version 5 (ERA5) against independent satellite retrievals indicate that TMeCSR's cloud and multiscale rain fields are better than those of ERA5. Furthermore, TMeCSR better captured the diurnal variability of rainfall and the statistical characteristics of MCSs. Forecasts initialized from TMeCSR also have more accurate rain and clouds than those initialized from ERA5. The TMeCSR and ERA5 forecasts have similar performances with respect to sounding and surface observations. These results indicate that TMeCSR is a promising MCS‐resolving data set for tropical MCS studies.

Chan, Man‐Yau↗