Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiscale Modeling and Simulation”

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 343 records · Page 19

Predicting and Controlling Corrosion (Abbreviated Final Report)

Corrosion accumulates significant replacement costs in the transportation, utility, manufacturing, and infrastructure sectors. It also underpins several of LLNL’s core missions in stockpile stewardship, defense, and energy security. Whereas most corrosion and aging models are empirically parameterized to describe well-defined conditions late in the reaction, the factors that determine the early stages of corrosion—during which mitigation could be most impactful— are poorly understood. This project addressed the critical need for new approaches to predict the kinetics of corrosion initiation based on firm physical and chemical understanding. The activities encompassed degradation of relevant metals in both hydrogen-rich and environmental corrosion scenarios. The team integrated state-of-the-art multiscale simulation, in situ characterization, and data science within three technical thrusts: hydriding of Ti alloys; aqueous corrosion of Al and Ni-Cr alloys; and degradation of additively manufactured 316L stainless steel. In each case, novel capabilities were developed to identify and track the impacts of key atomistic, compositional, and microstructural features on the metal systems. For hydriding, protocols were developed to tightly integrate multiscale models, advanced multimodal characterization, and machine learning to determine how hydrogen interacts with native passivating surface oxides and nucleates new undesired phases, shedding new light on the critical role of grain boundaries, interfaces, and atomically disordered regions. For aqueous corrosion, the project demonstrated methods to predict dissolution rates of metal surfaces in corrosive solutions, to measure and understand microstructural and grain orientation effects on corrosion susceptibility, and to investigate competing growth and dissolution kinetics of surface oxides. For additively manufactured metals, analysis using state-of-the-art microscopy techniques revealed the role of specific heterogeneities invoked during laser processing, including cellular structure, dislocations, and precipitates, on corrosion susceptibility. In addition to new capabilities and understanding, the project provided an avenue for workforce development, as well as key partnerships with stakeholders in corrosion science.

08 HYDROGEN↗

Validation Framework for Post-Combustion Carbon Capture CFD Simulations

First-principles based computational fluid dynamics (CFD) simulations are proposed as a fundamental tool for investigating solvent-based CO2 absorption in packed columns, due to the ability to accurately represent the underlying non-linear multiscale dynamics. In this work, we employ such models to investigate hydrodynamics of columns with structured by assessing the key hydrodynamic metrics, such as pressure drop and liquid holdup. Our models are validated with experimental data from a specifically designed column for this line of work. The test cases of gas and liquid flowrates and operating conditions were selected through a comprehensive sequential design of experiments approach offered by the CCSI2 toolset.

Panagakos, Grigorios↗

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗

MOOSE ProbML: Parallelizable Probabilistic Machine Learning and Uncertainty Quantification Capabilities

The Multiphysics Object Oriented Simulation Environment (MOOSE) is a widely used open- source finite element software for performing multiphysics multiscale simulations in a massively parallel fashion. Recently, the computational team at Idaho National Laboratory (INL) has implemented Probabilistic Machine Learning (ProbML) capabilities in MOOSE—in a parallelized fashion—and enable active learning with large-scale computational models for tasks such as surrogate model development, scale bridging, forward/inverse uncertainty quantification (UQ), Bayesian optimization, etc. This presentation summarizes these developments in MOOSE along with demonstrations on several real applications relevant to nuclear energy. At the fundamental level, samplers like Monte Carlo/Latin Hypercube, variance reduction, parallelized Markov Chain Monte Carlo (MCMC) support uncertainty propagation in both forward and inverse settings. These samplers can be integrated with the Gaussian processes (GP) suite in MOOSE, which offer several variants like scalar GPs, multi-output GPs, and deep GPs, to enable active learning. These GPs can be tuned using gradient-based optimization methods like Adam and its variants or gradient-free methods like the elliptical slice sampler (a variant of MCMC adept under Gaussian settings) for more complex covariance kernels or likelihoods whose gradient computations can be cumbersome. A variety of batch acquisition functions permit parallelized evaluation of the computational model and support different learning objectives with high efficiency like Bayesian inference, global surrogate development, optimization, etc. Furthermore, libtorch integration supports training, evaluation, and re-training of neural networks and other complex machine learning models in active learning settings. The impacts of these developments are shown on several real applications: (1) nuclear fuel inverse UQ and model inadequacy assessment using the Kennedy O’Hagan framework; (2) uncertainty aware surrogate modeling for additive manufacturing to predict field quantities; (3) nuclear reactor rare events analysis; and (4) complex fluid flow prediction using a global surrogate with quantified prediction uncertainty. Finally, the outlook of MOOSE ProbML is discussed for both outer-loop and inner-loop computations in the broad view to accelerate fuels and materials qualification, address gaps in knowledge and data, and assess new reactor/fuel systems.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Coupled Monte Carlo and thermal-fluid modeling of high temperature gas reactors using Cardinal

Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computa-tional fluid dynamics to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. Here, we provide an introduction to Cardinal's software design, data mapping, and multi -physics coupling strategy to highlight our approach to overcoming common challenges in multiphysics simulation. We then describe an application of Cardinal to prismatic High Temperature Gas Reactors (HTGRs) with various combinations of NekRS, OpenMC, BISON, and THM. A high-resolution coupling of NekRS, OpenMC, and BISON provides a reference solution at the unit cell level and shows excellent agree-ment with a lower-resolution coupling of THM, OpenMC, and BISON. A full core coupling of THM, OpenMC, and BISON resolving the three-dimensional conjugate heat transfer and sub-pin power distri-bution then provides detailed predictions of HTGR temperatures and the fission distribution.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

LDRD-DR: Accelerating Nuclear Fuel Qualification through Integrated Multiscale and Multiphysics Models [Slides]

Nuclear fuel qualification implies ensuring that a new fuel can safely operate in a reactor under both normal and accident conditions. Historically, extensive irradiation campaigns probe fuel performance relationships and interconnections in an integral sense, which is very time-consuming and costly. The presentation focus is on introducing advanced modeling and simulation into the qualification process in order to reduce the qualification timeline by targeting fewer integral tests

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Novel fractal-textured solar absorber surfaces for concentrated solar power

The solar absorptance of a thermal receiver surface significantly affects the photothermal conversion efficiency of concentrated solar power (CSP) plants. The development of low-cost Gen3 CSP systems calls for increasing the solar absorptance of the thermal receivers at operating temperatures above 750 °C. This study presents an innovative approach to the fractal, multiscale texturing of absorber surfaces to significantly enhance solar energy absorption for Gen3 CSP applications. The absorber surface is described in terms of its fractal parameters that are uniquely determined from surface profile measurements. The interaction of solar radiation with the fractal surface is numerically simulated by solving the governing Maxwell’s equations for electromagnetic wave propagation to investigate the effect of texturing on the absorptance of the surfaces. It is shown that absorption of the solar spectrum increases with an increase in the fractal dimension and the multiscale asperity height of the surface texturing. Here, the theoretical model is demonstrated to be in close agreement with experimental measurements of spectral absorptance of electrodeposited copper (Cu), copper mono-oxide (CuO), and copper-manganese oxide (CuMnO) surfaces that are textured to produce a range of fractal parameters by tailoring the deposition parameters. Fractal surface texturing is shown to reduce reflectance by over one order of magnitude, yielding an absorptance of greater than 0.985 for CuMnO. For the first time, the study presents a significant effective means of increasing solar absorptance and a fundamental theoretical description of the underlying physics.

14 SOLAR ENERGY↗

From Cell to System: Accelerated hpc Simulations of BESS Aging under Frequency Regulation and Arbitrage use cases

Lithium-ion battery energy storage systems (BESS) packs have emerged as a leading solution for grid-scale energy storage, enhancing resiliency and balancing load fluctuations. Yet, experimental characterization of large-format LIB packs-particularly to assess performance and degradation over hundreds of cycles - demands substantial hardware investment and multi-year testing campaigns. In this work, we couple a hierarchical, physics-based modeling framework agnostic to electrode chemistries with high-performance computing to accelerate systems level evaluation by upto two orders of magnitude. Building on the open-source liionpack platform, we implement cell, module, and pack-scale electrochemical models enriched with mechanistic aging mechanisms and deploy them on an HPC cluster to simulate 150−200kWh systems over 500 - 1,000 cycles with in days. We subject these virtual B ESS to both constant-current cycling and realistic grid service profiles spanning frequency regulation, ramp-rate support, and energy arbitrage-and quantify the resulting degradation patterns. Our results reveal that localized cell aging can induce substantial nonuniformity at module and pack levels, with service-specific cycling protocols driving distinct aging modes. This rapid, multiscale modeling approach provides a powerful design-space exploration tool for optimizing electrical architecture, control strategies, and operational schedules to prolong pack lifetime and lower total cost of ownership.

Ayalasomayajula, Surya [ORNL] (ORCID:0009000860788↗

Water Cycle-Driven Infectious Diseases as Multiscale, Reliable, Continuously Updating Water Cycle Sensors

We propose a novel synthesis of observational and simulated data (climatological and biological) to enhance understanding of the real-world interplay between climate (here, the water cycle) and the epidemiology of water cycle-driven infectious disease. Aligning with Focal Area 2, predictive modeling using AI techniques, we will develop systems of hierarchical models to discover latent features of the water cycle. Our approach will leverage state-of-the-art in Artificial Intelligence (AI) to measure the degree to which climate change-driven shifts in water cycle can be predicted by supplementing sparse and irregular climate data with water cycle-driven infectious disease resources.

54 ENVIRONMENTAL SCIENCES↗

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗

Mesh-based super-resolution of fluid flows with multiscale graph neural networks

A graph neural network (GNN) approach is introduced in this work which enables mesh-based three-dimensional super-resolution of fluid flows. In this framework, the GNN is designed to operate not on the full mesh-based field at once, but on localized meshes of elements (or cells) directly. To facilitate mesh-based GNN representations in a manner similar to spectral (or finite) element discretizations, a baseline GNN layer (termed a message passing layer, which updates local node properties) is modified to account for synchronization of coincident graph nodes, rendering compatibility with commonly used element-based mesh connectivities. Furthermore, the architecture is multiscale in nature, and is comprised of a combination of coarse-scale and fine-scale message passing layer sequences (termed processors) separated by a graph unpooling layer. The coarse-scale processor embeds a query element (alongside a set number of neighboring coarse elements) into a single latent graph representation using coarse-scale synchronized message passing over the element neighborhood, and the fine-scale processor leverages additional message passing operations on this latent graph to correct for interpolation errors. Demonstration studies are performed using hexahedral mesh-based data from Taylor–Green Vortex and backward-facing step flow simulations at Reynolds numbers of 1600 and 3200. Through analysis of both global and local errors, the results ultimately show how the GNN is able to produce accurate super-resolved fields compared to targets in both coarse-scale and multiscale model configurations. Reconstruction errors for fixed architectures were found to increase in proportion to the Reynolds number. Geometry extrapolation studies on a separate cavity flow configuration show promising cross-mesh capabilities of the super-resolution strategy.

Backward-facing step↗

Multidimensional modeling of fuel-cladding friction in an LWR fuel rod

Solving finite element problems with friction often increases the level of difficulty to obtain properly converged solutions. This type of challenge becomes more salient when advanced, possibly multiscale, material models are employed to capture the thermomechanical behavior of fuel and cladding materials. Here the present work details our recent developments in a nuclear fuel performance finite element code for the systematic consideration of friction in nuclear reactor finite element simulations. We show the application of friction and its effects on the mechanics of light-water reactor rods accounting for various fuel constitutive modeling techniques, model dimensionalities, pellet assumed geometries, and power conditions. In particular, we focus on the fuel rod mechanical behavior as it relates to fuel constitutive models, sensitivity to the coefficient of friction, pellet states of stress, and rod elongation. We discuss the trade-offs between the various multidimensional modeling options and highlight the relevance of frictional effects in the prediction of the fuel rod deformation and interfacial stresses. To relate our modeling results with actual reactor operation, simulations including frictional effects are compared with fuel rod elongation experimental data and some challenges for carrying out a full validation of the axial mechanics are discussed.

42 ENGINEERING↗

Multiscale Framework Simulates Utility-Scale Wind Plant in Its Natural Environment. First-ever simulation could lead to a vital understanding of wind plant performance

Wind plant performance depends on local environmental conditions—and understanding those conditions is crucial to research on wind plant dynamics. Factors such as sloping terrain, surface type and cover, atmospheric heating and cooling, and weather events all play a role in determining the winds that turbines convert into power. However, most models used for wind plant design and operation do not sufficiently account for environmental effects. This can limit our understanding of wind plant dynamics and, thus, our ability to design, site, and operate wind plants for maximum power production and reliability. To meet this challenge, Arthur and other LLNL researchers partnered with scientists at the National Center for Atmospheric Research; the University of California, Berkeley; the University of Colorado, Boulder; and Texas Tech University to develop and demonstrate a novel, mesoscale-tomicroscale wind plant modeling framework.

17 WIND ENERGY↗

Full Waveform Inversion Tomography for Improved Seismic Waveform Simulations and Lithospheric Structure of Northwest China and Surrounding Regions

We report a new model of three-dimensional (3D) anisotropic seismic wavespeed structure of the lithospheric of northwestern China, western Mongolia and east }|ern Kazakhstan from full waveform inversion tomography using openly available broadband data. To investigate the efficacy of a conservative multiscale inversion methodology, we used a simple isotropic starting model.

58 GEOSCIENCES↗

NeuroSEM: A hybrid framework for simulating multiphysics problems by coupling PINNs and spectral elements

Multiphysics problems that are characterized by complex interactions among fluid dynamics, heat transfer, structural mechanics, and electromagnetics, are inherently challenging due to their coupled nature. While experimental data on certain state variables may be available, integrating these data with numerical solvers remains a significant challenge. Physics-informed neural networks (PINNs) have shown promising results in various engineering disciplines, particularly in handling noisy data and solving inverse problems in partial differential equations (PDEs). However, their effectiveness in forecasting nonlinear phenomena in multiphysics regimes, particularly involving turbulence, is yet to be fully established. Here, this study introduces NeuroSEM, a hybrid framework integrating PINNs with the highfidelity Spectral Element Method (SEM) solver, Nektar++. NeuroSEM leverages the strengths of both PINNs and SEM, providing robust solutions for multiphysics problems. PINNs are trained to assimilate data and model physical phenomena in specific subdomains, which are then integrated into the Nektar++ solver. We demonstrate the efficiency and accuracy of NeuroSEM for thermal convection in cavity flow and flow past a cylinder. The framework effectively handles data assimilation by addressing those subdomains and state variables where the data is available. We applied NeuroSEM to the Rayleigh-B´enard convection system, including cases with missing thermal boundary conditions and noisy datasets. Finally, we applied the proposed NeuroSEM framework to real particle image velocimetry (PIV) data to capture flow patterns characterized by horseshoe vortical structures. Our results indicate that NeuroSEM accurately models the physical phenomena and assimilates the data within the specified subdomains. The framework’s plug-and-play nature facilitates its extension to other multiphysics or multiscale problems. Furthermore, NeuroSEM is optimized for efficient execution on emerging integrated GPU-CPU architectures. This hybrid approach enhances the accuracy and efficiency of simulations, making it a powerful tool for tackling complex engineering challenges in various scientific domains.

42 ENGINEERING↗

Electro-reduction of organics on metal cathodes: A multiscale-modeling study of benzaldehyde on Au (111)

Here, we present a multiscale modeling study on the electro-reduction of benzaldehyde (BZY) on Au(111) under realistic electrochemical conditions. To model the electrochemical cell, we adopt a capacitor model in which complex solvents are confined between a Au cathode and a carbon anode. Classical molecular dynamics simulations reveal that electrode charge density and the presence of alcohol show strong effects on the density, adsorption geometry and dynamics of benzaldehyde on the Au electrode. Under charging conditions, the surface concentration of benzaldehyde on the Au electrode decreases, while the content of other species increases. Finally, we proposed a scheme that correlates the electric current running through the Au/solvent interface with the applied bias. This study provides a molecular level understanding of how solvent composition, in this case water/alcohol content, controls the activity of electrocatalytic hydrogenation

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiscale Control of Generic Second Order Traffic Models by Driver-Assist Vehicles

We study the derivation of generic high order macroscopic traffic models from a follow-the-leader particle description via a kinetic approach. First, we recover a third order traffic model as the hydrodynamic limit of an Enskog-type kinetic equation. Next, we introduce in the vehicle interactions a binary control modeling the automatic feedback provided by driver-assist vehicles and we upscale such a new particle description by means of another Enskog-based hydrodynamic limit. The resulting macroscopic model is now a generic second order model (GSOM), which contains in turn a control term inherited from the microscopic interactions. We show that such a control may be chosen so as to optimize global traffic trends, such as the vehicle flux or the road congestion, constrained by the GSOM dynamics. By means of numerical simulations, we investigate the effect of this control hierarchy in some specific case studies, which exemplify the multiscale path from the vehiclewise implementation of a driver-assist control to its optimal hydrodynamic design.

GSOM↗