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 451 records · Page 25

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↗

Filament Wound Composite Analysis Using the NASA Multiscale Analysis Tool (NASMAT) and Finite Element Analysis

Fiber reinforced composite materials, owing to their tailorable thermomechanical and functional properties, allow one to produce a structure that is stronger, stiffer, and lighter than its metal counterpart while performing the same function, yielding a more efficient structure. This not only allows for the improvement of current technologies like aircraft structures, but also enables new technologies like gaseous hydrogen storage for mobility applications, which are otherwise impractical when manufactured using traditional metals due to weight and space restrictions or material embrittlement. However, the use of composites imposes greater design and manufacturing challenges on an engineer, since they are heterogenous, having a distinct structure across multiple length scale, behave generally anisotropically at the structural level and require complex manufacturing and processing methods. Capturing this complex behavior requires detailed numerical simulations, including the modeling of microstructural features like undulations, voids, and fiber alignment. In this paper, multiple repeating unit cells (RUCs), representing filament wound composites, are developed (via a script provided in the Appendix) and analyzed. The refinement of these RUCs is varied, and the analyses are performed using both the Abaqus finite element software and the NASA Multiscale Analysis Tool (NASMAT). A study is undertaken to compare the predicted effective elastic properties of the wound RUC to a laminate representation of the wound RUC, which neglects the undulations. Additionally, two different sets of periodic boundary conditions (PBCs) have been examined. One approximates the real boundary conditions using a standard approach and the other represents the PBCs exactly through the use of an offset. Lastly, a comparison of the local elastic stress fields is made among the models and approaches. Since wound structures are often approximated as laminated structures, it is important to understand the degree to which this assumption is valid, namely by first comparing the elastic constants and local elastic fields. This will provide, on the one hand, information concerning the bulk mechanical behavior and, on the other hand, insights concerning local load distributions and likely damage initiation sites.

Composites↗

Coupled fvGCM-GCE Modeling System, 3D Cloud-Resolving Model and Cloud Library

Recent GEWEX Cloud System Study (GCSS) model comparison projects have indicated that cloud-resolving models (CRMs) agree with observations better than traditional singlecolumn models in simulating various types of clouds and cloud systems from Merent geographic locations. Current and future NASA satellite programs can provide cloud, precipitation, aerosol and other data at very fine spatial and temporal scales. It requires a coupled global circulation model (GCM) and cloudscale model (termed a super-parameterization or multiscale modeling framework, MMF) to use these satellite data to improve the understanding of the physical processes that are responsible for the variation in global and regional climate and hydrological systems. The use of a GCM will enable global coverage, and the use of a CRM will allow for better and more sophisticated physical parameteridon NASA satellite and field campaign cloud related datasets can provide initial conditions as well as validation for both the MMF and CRMs. A seed fund is available at NASA Goddard to build a MMF based on the 2D Goddard cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM). A prototype MMF in being developed and production nms will be conducted at the beginning of 2005. In this talk, I will present: (1) A brief review on GCE model and its applications on precipitation processes, (2) The Goddard MMF and the major difference between two existing MMFs (CSU MMF and Goddard MMF), (3) A cloud library generated by Goddard MMF, and 3D GCE model, and (4) A brief discussion on the GCE model on developing a global cloud simulator.

Tao, Wei-Kuo↗

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↗

Predicting Tropical Cyclogenesis with a Global Mesoscale Model: Hierarchical Multiscale Interactions During the Formation of Tropical Cyclone Nargis(2008)

Very severe cyclonic storm Nargis devastated Burma (Myanmar) in May 2008, caused tremendous damage and numerous fatalities, and became one of the 10 deadliest tropical cyclones (TCs) of all time. To increase the warning time in order to save lives and reduce economic damage, it is important to extend the lead time in the prediction of TCs like Nargis. As recent advances in high-resolution global models and supercomputing technology have shown the potential for improving TC track and intensity forecasts, the ability of a global mesoscale model to predict TC genesis in the Indian Ocean is examined in this study with the aim of improving simulations of TC climate. High-resolution global simulations with real data show that the initial formation and intensity variations of TC Nargis can be realistically predicted up to 5 days in advance. Preliminary analysis suggests that improved representations of the following environmental conditions and their hierarchical multiscale interactions were the key to achieving this lead time: (1) a westerly wind burst and equatorial trough, (2) an enhanced monsoon circulation with a zero wind shear line, (3) good upper-level outflow with anti-cyclonic wind shear between 200 and 850 hPa, and (4) low-level moisture convergence.

Shen, B.-W.↗

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↗

Advances in Design Capabilities for Planetary Missions from the NASA Entry Systems Modeling and Instrumentation Portfolio

The Entry Systems Modeling project (ESM) is supported by both the NASA Space Technology and the Science Mission Directorates and focuses on developing simulation tools and validated models for characterizing the performance of entry systems tailored to planetary destinations across the Solar System. ESM is organized into six technical capability areas that together address all relevant factors related to spacecraft entry, as well as some aspects of descent: Thermal Protection System (TPS) Materials; Aerothermodynamics; Entry & Descent Vehicle Dynamics; Guidance, Navigation, and Control; Vehicle Systems Analysis; and Advanced Tools and Numerical Methods. Development within the capability areas is undertaken explicitly with a focus on transition and infusion to science missions, human exploration missions, and commercial space activities. The present talk details developments that specifically impact science missions, including simulation tool capabilities that aid in mission design and model development to understand entry system performance at a given destination. Examples of the successful infusion and transition of such project outcomes to science missions also are provided. Several simulation tool development efforts within ESM have resulted in new design capabilities for missions. One such outcome is improved toolsets for mission trajectory and concept of operations design. Specifically, an initiative to couple a leading tool for entry, ascent/descent, and orbital trajectory optimization (Program to Optimize Simulated Trajectories II or POST2) to those used within the Agency for interplanetary trajectory optimization (Copernicus and Monte) has made substantial progress, with the outcomes to date promising to allow efficient trajectory optimization across mission phases. Additionally, toolchains for the evaluation of vehicle performance during entry and descent have been developed that allow assessment of multi-dimensional aeroheating on detailed vehicle geometries, characterization of deployment and inflation of parachutes, and assessment of vehicle dynamic stability during descent. These capabilities are achieved by coupling diverse sets of physics together – material response, computational fluid dynamics, radiation, and vehicle dynamics – to suitably describe complex entry and descent phenomena. Several model development and validation efforts for specific destinations and entry regimes also are underway within the ESM project. For instance, new experimental capabilities to validate radiation models at low densities/high altitudes recently have been established with project support, specifically the Low-Density Shock Tube (LDST) at the NASA Ames Research Center Electric Arc Shock Tube (EAST) facility. The LDST is being leveraged to develop improved models of shock layer kinetics and radiation in Titan atmospheres, while future studies will be conducted in the LDST and the existing high velocity shock tube to provide validation data for radiation models of Venus, Ice Giants, and Mars atmospheres. Models describing the aerothermal and thermo-structural performance of Thermal Protection System (TPS) materials has been another focus, with multiscale modeling activities on-going for the two leading TPS materials applicable to a range of entry conditions and science missions: the Phenolic-Impregnated Carbon Ablator (PICA) and woven materials like 3D Mid-Density Carbon Phenolic (3MDCP). A continual effort is made to infuse and transition outcomes from ESM simulation tool and model development activities into relevant science missions. Significant progress has been made on this front, with missions such as Dragonfly, DAVINCI, and Mars Missions benefitting from project outcomes. The groundwork also is being laid to provide insights into forward looking missions to Gas/Ice Giants as well as for potential sample returns.

Justin Haskins↗

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 Modeling of Fracture Strength in Fibrous Thermal Protection System Materials

This work presents a multiscale modeling approach to predict the fracture strength of fibrous Thermal Protection System (TPS) materials. The model assumes that system failure is initiated at the joints between individual fibers. We investigated three distinct TPS compositions: amorphous silica, alumina and aluminosilicate fibers. Molecular dynamics (MD) simulations were employed to determine the fracture strength values of these fiber joints for both material systems. These fracture strength values were then integrated into simulations of 3D randomly populated fiber structures, where tensile load transfer occurs through the fiber joints. These microscale properties are upscaled through a renormalization approach [1] to predict macroscale tensile strength of 3D random fiber networks, accounting for joint-dominated failure and effective load-bearing area. The study concludes by demonstrating the resulting strength variation as a function of material composition, fiber density, and morphology. We also show validation of results by comparing them against explicit fiber finite element (FE) modeling [2] where fiber joint fracture is represented by cohesive elements.

Jaehyun Cho↗