Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The Interaction between the Nocturnal Amazonian Low-Level Jet and Convection in CESM

A nocturnal Amazonian low-level jet (ALLJ) was recently diagnosed using reanalysis data. This work assesses the ability of CESM1.2.2 to reproduce the jet and explores the mechanisms by which the ALLJ influences convection in the Amazon. The coupled CESM simulates the nocturnal ALLJ realistically, while CAM5 does not. A low-level cold air temperature bias in the eastern Amazon exists in CAM5, and thus the ALLJ is weaker than observed. However, a cold SST bias over the equatorial North Atlantic in the coupled model offsets the cold air temperature bias, producing a realistic ALLJ. Climate models significantly underestimate March–May (MAM) precipitation over the eastern Amazon. Here, we ran two sensitivity experiments using the coupled CESM by adding bottom-heavy diabatic heating at noon and midnight for 2.5 h along the coastal Amazon during MAM to mimic the occurrence of shallow precipitating convection. When heating is added during the early afternoon, coastal convection deepens and the ALLJ transports moisture inland from the ocean, preconditioning the environment for deep convective development during the ensuing hours. The increased convection over the eastern Amazon also moderately alleviates the equatorial Atlantic westerly wind bias, leading to deepening of the east Atlantic thermocline in the following months and partially improving the simulated June–August (JJA) Atlantic cold tongue in the coupled model. When heating is added at night, coastal convection does not strengthen as much and the ALLJ transports less moisture. Improvements in the simulated Atlantic winds and SST are negligible. Therefore, diurnal circulations matter to the organization of convection and rain across the Amazon, with impacts over the equatorial Atlantic.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1999

Inlet-Compressor Analysis using Coupled CFD Codes

Propulsion performance and operability are key factors in the development of a successful aircraft. For high-speed supersonic aircraft, mixed-compression inlets offer high performance but are susceptible to an instability referred to as unstart. An unstart occurs when a disturbance originating in the atmosphere or the engine causes the shock system to be expelled from the inlet. This event can have adverse effects on control of the aircraft, which is unacceptable for a passenger plane such as the high speed civil transport (HSCT). The ability to predict the transient response of such inlets to flow perturbations is, therefore, important to the proper design of the inlet and the control measures used to prevent unstart. Computational fluid dynamics (CFD) is having an increasing role in the analysis of individual propulsion components. Isolated inlet studies are relatively easy to perform, but a major uncertainty is the boundary condition used at the inlet exit to represent the engine - the so-called compressor face boundary condition. A one-dimensional (1-D) Euler inlet simulation showed that the predicted inlet unstart tolerance to free-stream pressure perturbations can vary by as much as a factor of about six, depending on the boundary condition used. Obviously a thorough understanding of dynamic interactions between inlets and compressors/fans is required to provide the proper boundary condition. To aid in this understanding and to help evaluate possible boundary conditions, an inlet-engine experiment was conducted at the University of Cincinnati. The interaction of acoustic pulses, generated in the inlet, with the engine were investigated. Because of the availability of experimental data for validation, it was decided to simulate the experiment using CFD. The philosophy here is that the inlet-engine system is best simulated by coupling (existing) specialized CFD component-codes. The objectives of this work were to aid in a better understanding of inlet-compressor interaction physics and the formulation of a more realistic compressor-face boundary condition for time-accurate CFD simulations of inlets. Previous simulations have used 1-D Euler engine simulations in conjunction with 1-D Euler and axisymmetric Euler inlet simulations. This effort is a first step toward CFD simulation of an entire engine by coupling multidimensional component codes.

Cole, Gary↗
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↗
NASA NTRS2015

Downscaling GISS ModelE Boreal Summer Climate over Africa

The study examines the perceived added value of downscaling atmosphere-ocean global climate model simulations over Africa and adjacent oceans by a nested regional climate model. NASA/Goddard Institute for Space Studies (GISS) coupled ModelE simulations for June- September 1998-2002 are used to form lateral boundary conditions for synchronous simulations by the GISS RM3 regional climate model. The ModelE computational grid spacing is 2deg latitude by 2.5deg longitude and the RM3 grid spacing is 0.44deg. ModelE precipitation climatology for June-September 1998-2002 is shown to be a good proxy for 30-year means so results based on the 5-year sample are presumed to be generally representative. Comparison with observational evidence shows several discrepancies in ModelE configuration of the boreal summer inter-tropical convergence zone (ITCZ). One glaring shortcoming is that ModelE simulations do not advance the West African rain band northward during the summer to represent monsoon precipitation onset over the Sahel. Results for 1998-2002 show that onset simulation is an important added value produced by downscaling with RM3. ModelE Eastern South Atlantic Ocean computed sea-surface temperatures (SST) are some 4 K warmer than reanalysis, contributing to large positive biases in overlying surface air temperatures (Tsfc). ModelE Tsfc are also too warm over most of Africa. RM3 downscaling somewhat mitigates the magnitude of Tsfc biases over the African continent, it eliminates the ModelE double ITCZ over the Atlantic and it produces more realistic orographic precipitation maxima. Parallel ModelE and RM3 simulations with observed SST forcing (in place of the predicted ocean) lower Tsfc errors but have mixed impacts on circulation and precipitation biases. Downscaling improvements of the meridional movement of the rain band over West Africa and the configuration of orographic precipitation maxima are realized irrespective of the SST biases.

Regional climate models↗
NASA NTRS2015

Development of a Grid-Independent Geos-Chem Chemical Transport Model (v9-02) as an Atmospheric Chemistry Module for Earth System Models

The GEOS-Chem global chemical transport model (CTM), used by a large atmospheric chemistry research community, has been re-engineered to also serve as an atmospheric chemistry module for Earth system models (ESMs). This was done using an Earth System Modeling Framework (ESMF) interface that operates independently of the GEOSChem scientific code, permitting the exact same GEOSChem code to be used as an ESM module or as a standalone CTM. In this manner, the continual stream of updates contributed by the CTM user community is automatically passed on to the ESM module, which remains state of science and referenced to the latest version of the standard GEOS-Chem CTM. A major step in this re-engineering was to make GEOS-Chem grid independent, i.e., capable of using any geophysical grid specified at run time. GEOS-Chem data sockets were also created for communication between modules and with external ESM code. The grid-independent, ESMF-compatible GEOS-Chem is now the standard version of the GEOS-Chem CTM. It has been implemented as an atmospheric chemistry module into the NASA GEOS- 5 ESM. The coupled GEOS-5-GEOS-Chem system was tested for scalability and performance with a tropospheric oxidant-aerosol simulation (120 coupled species, 66 transported tracers) using 48-240 cores and message-passing interface (MPI) distributed-memory parallelization. Numerical experiments demonstrate that the GEOS-Chem chemistry module scales efficiently for the number of cores tested, with no degradation as the number of cores increases. Although inclusion of atmospheric chemistry in ESMs is computationally expensive, the excellent scalability of the chemistry module means that the relative cost goes down with increasing number of cores in a massively parallel environment.

CTM↗
DOE OSTI2021

Design principles of tandem cascade photoelectrochemical devices

Cascade photoelectrocatalysis (PEC) is a possible method to improve the selectivity of solar-driven CO 2 reduction (CO 2 R). This concept can be realized by coupling different CO 2 R catalysts to different subcells in a multijunction photovoltaic (PV) stack. Efficient implementation will require careful tuning of the photocurrents and design of the photovoltages provided by the subcells to the CO 2 R catalysts in such a way as to facilitate the target reaction. Here, we outline the design principles of the tandem PEC approach using two-step conversion of CO 2 to ethylene in aqueous electrolyte, via a CO intermediate, as a model system. To perform this reaction, the first coupled PV-catalyst component should provide 4 electrons to reduce 2 molecules of CO 2 to CO; the second component should provide 8 electrons to reduce 2 CO molecules to C 2 H 4 . Based on known CO 2 R catalysts, the overpotential required to produce CO can be less than that required to reduce it to ethylene, creating the opportunity for improved efficiency. Cascade PEC can be realized in a three-terminal tandem (3TT) configuration using III–V-semiconductor based subcells coupled to Au (produces CO intermediate) and Cu (converts CO to ethylene). The current to each catalyst can be controlled by the area of the subcell exposed to the electrolyte, and the photovoltage is determined by the materials selected and device configuration. Operating conditions are found by simulating the coupled system using the open-source circuit simulator SPICE (Simulation Program with Integrated Circuits Emphasis). We identify conditions under which a 3TT configuration can have a higher solar to chemical conversion efficiency compared to a two-terminal two-junction tandem (2T 2J) with the same absorbers and a Cu catalyst only. We also show that 3TT PEC devices can be less sensitive to variations in catalyst activity compared to 2T devices. Lastly, we discuss the applications of cascade PEC to CO 2 reduction, using different intermediates, and to other chemical networks.

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

The PENELOPE Physics Models and Transport Mechanics. Implementation into Geant4

A translation of the PENELOPE physics subroutines to C++, designed as an extension of the GEANT4 toolkit, is presented. The Fortran code system PENELOPE performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. PENELOPE implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of GEANT4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original PENELOPE programs. The combined code, named PENG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of GEANT4.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗
DOE OSTI2023

Convectively coupled Rossby–Gravity waves in a field campaign: How they are captured in reanalysis products

Abstract Convectively coupled equatorial waves are a significant source of atmospheric variability in the tropics. Current numerical models continue to struggle in simulating the coupled diabatic heating fields that are responsible for the development and maintenance of these waves. This study investigates how the diabatic fields associated with Mixed Rossby–Gravity waves (MRGs) are represented in four reanalysis products by using a unique observational dataset from the TRMM‐KWAJEX (Tropical Rainfall Measuring Mission‐Kwajalein Experiment) field campaign. These reanalyses include ERA5, Japanese 55‐year Reanalysis (JRA‐55), Climate Forecast System Reanalysis (CFSR), and Modern‐Era Retrospective Analysis for Research and Applications (MERRA). We found that all four reanalyses captured the MRG structures in winds and temperature, and to a lesser degree in the humidity field except in the boundary layer. However, only the ERA5 and MERRA reanalyses captured the gradual rise and succession of the diabatic heating from boundary layer turbulence, shallow convection, cumulus congestus, and deep convection within the waves. ERA5 is the only product that also captured the gradual rise of the subgrid‐scale vertical transport of moist static energy. All reanalysis products underestimated the diabatic heating from cumulus congestus. Results provide observational basis on what aspects of MRG can be trusted and what cannot in the reanalysis products.

Wang, Xiaocong↗
DOE OSTI2025

Software For Advanced Large-scale Analysis Of Magnetic Confinement For Numerical Design, Engineering & Research (salamander)

As magnetic confinement fusion energy gains traction internationally to enable abundant energy production, designing components for fusion systems is a pressing challenge. During the planned lifetime of a fusion device, components evolve in extreme environments and must withstand large, repeated thermal loads and bombardment by 14 MeV neutrons, plasma ions, and neutral particles (deuterium, tritium, and helium), corrosive conditions, etc. All these physical processes take place simultaneously, interact in intricate ways, and impose important constraints that can affect performance. Experimental data is rare and costly to obtain, making design particularly challenging. Predictive computational frameworks must be an integral part of an accelerated and cost-effective design process by modeling fusion system performance in simulated environments. To better understand component degradation and operational impacts on their performance, the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER) is designed as an open-source, fully integrated, multiphysics, multiscale, NQA-1 compliant framework facilitating 3D, high-fidelity fusion system modeling. To that end, SALAMANDER is a MOOSE-based framework, and therefore leverages MOOSE upstream libraries such as PETSc and libMesh to deliver sophisticated finite element, finite volume, and nonlinear solver technology for fusion energy simulations. SALAMANDER couples MOOSE physics module capabilities—such as thermal hydraulics, heat conduction, Navier-Stokes, and thermomechanics—with tritium transport via TMAP8, neutronics via Cardinal, and nascent particle-in-cell capabilities. Direct simulation Monte Carlo methods will be used to address neutral transport near the walls. By coupling all these physics in an integrated application, SALAMANDER will enable high-fidelity modeling of irradiation levels and plasma exposure conditions of plasma facing components and their impact on heat and tritium distributions, as well as the resulting mechanical constraints experienced by the plasma facing components and performance of blanket systems. Furthermore, SALAMANDER will be particularly suited for engineering studies thanks to the stochastic tool module readily available in MOOSE, allowing for extended uncertainty quantification and risk analysis studies. It is also able to use computer-aided design (CAD) meshes to model complex geometries, which is indispensable for fusion systems. SALAMANDER therefore supports design, safety, engineering, and research projects for magnetic confinement fusion systems

Simon, Pierre-Clement [Idaho National Laboratory (↗
DOE OSTI2021

Long-term single-column model intercomparison of diurnal cycle of precipitation over midlatitude and tropical land

General Circulation Models (GCMs) have for decades exhibited difficulties in modelling the diurnal cycle of precipitation (DCP). This issue can be related to inappropriate representation of the processes controlling sub-diurnal phenomena like convection. In this study, 11 single-column versions of GCMs are used to investigate the interactions between convection and environmental conditions, processes that control nocturnal convections, and the transition from shallow to deep convection on a diurnal time-scale. Long-term simulations are performed over two continental land sites: the Southern Great Plains (SGP) in the USA for 12 summer months from 2004 to 2015 and the Manacapuru site at the central Amazon (MAO) in Brazil for two full years from 2014 to 2015. The analysis is done on two regimes: afternoon convective regime and nocturnal precipitation regime. Most models produce afternoon precipitation too early, likely due to the missing transition of shallow-to-deep convection in these models. At SGP, the unified convection schemes better simulate the onset time of precipitation. At MAO, models produce the heating peak in a much lower level compared with observation, indicating too shallow convection in the models. For nocturnal precipitation, models that produce most of nocturnal precipitation all allow convection to be triggered above the boundary layer. This indicates the importance of model capability to detect elevated convection for simulating nocturnal precipitation. Sensitivity studies indicate that (a) nudging environmental variables towards observations has a minor impact on DCP, (b) unified treatment of shallow and deep convection and the capability to capture mid-level convection can help models better capture DCP, and (c) the interactions of the atmosphere with other components in the climate system (e.g. land) are also important for DCP simulations in coupled models. Finally, these results provide long-term statistical insights on which physical processes are essential in climate models to simulate DCP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Implementation of the Windowed Multipole Method in Shift

The windowed multipole (WMP) method has been implemented in the Shift Monte Carlo (MC) radiation transport code with support for both CPU and GPU execution. With this method, small WMP data libraries (~100 MB) can be used to accurately Doppler broaden cross sections to arbitrary temperatures “on the fly” during an MC simulation. This approach yields significant memory savings relative to traditional methods, making it ideal for high-fidelity analysis such as coupled multiphysics simulations. This document provides the exact forms of the WMP equations used by Shift, as well as a detailed description of the structure of WMP HDF5 data files provided by the Massachusetts Institute of Technology (MIT). The Shift implementation has been validated against the OpenMC radiation transport code, with excellent agreement demonstrated for 70 nuclides across an operative range of temperatures. CPU and GPU performance testing using a small module reactor (SMR) problem demonstrated that this method decreases the neutron tracking rate by a factor of ~2 on the Summitdev machine. A new set of WMP data being developed in-house will employ novel methods to improve tracking rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Complementary Study of Radiative Heat Transfer and Flow Physics from Moderate-scale Hydrocarbon Pool Fire Simulations

As part of the Advanced Simulation and Computing Verification and Validation (ASCVV) program, a 0.3-m diameter hydrocarbon pool fire with multiple fuels was modeled and simulated. In the study described in this report, systematic examination was performed on the radiation model used in a series of coupled Fuego/Nalu simulations. A calibration study was done with a medium-scale methanol pool fire and the effect of calibration traced throughout the radiation model. This analysis provided a more detailed understanding of the effect of radiation model parameters on each other and on other quantities in the simulations. Heptane simulation results were also examined using this approach and possible areas for further improvement of the models were identified. The effect of soot on radiative losses was examined by comparing heptane and methanol results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2019

Equilibrium Climate Sensitivity Estimated by Equilibrating Climate Models

The methods to quantify equilibrium climate sensitivity are still debated. We collect millennial length simulations of coupled climate models and show that the global mean equilibrium warming is higher than those obtained using extrapolation methods from shorter simulations. Specifically, 27 simulations with 15 climate models forced with a range of CO2 concentrations show a median 17% larger equilibrium warming than estimated from the first 150 years of the simulations. The spatial patterns of radiative feedbacks change continuously, in most regions reducing their tendency to stabilizing the climate. In the equatorial Pacific, however, feedbacks become more stabilizing with time. The global feedback evolution is initially dominated by the tropics, with eventual substantial contributions from the midlatitudes. Time dependent feedbacks underscore the need of a measure of climate sensitivity that accounts for the degree of equilibration, so that models, observations, and paleo proxies can be adequately compared and aggregated to estimate future warming.

Maria Rugenstein↗
DOE OSTI2023

The MOOSE fluid properties module

The fluid properties module in MOOSE~\cite{lindsay2022moose} serves a variety of fluid simulation applications based on MOOSE, including the MOOSE Navier Stokes module~\cite{moose_ns}, Pronghorn~\cite{pgh}, the MOOSE thermal hydraulics module, SAM~\cite{sam}, RELAP-7~\cite{relap7}, Sockeye~\cite{sockeye}, Pronghorn-subchannel~\cite{subchannel} and the MOOSE porous flow module~\cite{porous}. These applications are used to solve coarse mesh multi-dimensional thermal-hydraulics~\cite{pgh}, 1D systems analysis~\cite{sam,relap7} in nuclear reactor analysis, heat pipe modeling~\cite{sockeye} and porous flow simulations~\cite{porous} for underground gas storage and water seepage. The use of consistent fluid properties across fluid flow applications facilitates coupled flow simulations~\cite{anl_sam_pgh,osti_1889653}. Each application has historically driven the implementation of several fluid properties, which were later extended to be compatible with other applications. The unique diversity of applications of the module, due to its presence in MOOSE, has driven its expansion to new fluids, such as advanced nuclear reactor coolants and, more recently, arbitrary functions or tables-based property definitions, as detailed in section~\ref{content}, as well as numerous thermophysical properties and variable sets, as detailed in subsection~\ref{sec:prop}. The need for different discretizations of flow equations based on the compressibility of the fluid has motivated support for both primitive (pressure- and temperature-based) and conservative (internal energy- and specific volume-based) flow variables in the module. Thermodynamic relations are used to convert between these two formulations, as needed. The module serves a dual purpose of providing fluid properties for direct use in flow simulations and facilitating the implementation of user-specific fluid properties. Contributions of new properties for existing fluids or new fluids are strongly encouraged.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Studies of the influence of heat pipe model and coupling strategy in heat pipe microreactor simulations using Sockeye

This work demonstrates various modeling approaches for thermally coupling heat pipes to a graphite monolithic block in a microreactor context, accounting for the presence of a gap. We show that heat pipe model selection has a significant impact on the transient behavior of a three-dimensional coupled microreactor assembly problem by comparing results obtained using different heat pipe models from the heat pipe code Sockeye. Additionally, we show that the choice of coupling strategy has significant consequences for both accuracy and performance. Furthermore, we found that our heat flux transfer strategy exhibited greater accuracy than our temperature transfer strategy, even when comparing a loose coupling of the heat flux transfer strategy to a tight coupling of the temperature transfer strategy. Although this research uses only theoretical test problems, it provides important insights in modeling thermal fluids phenomena in heat pipe microreactors.

Direwolf↗
DOE OSTI2022

Multi-fidelity modeling to predict the rheological properties of a suspension of fibers using neural networks and Gaussian processes

Unveiling the rheological properties of fiber suspensions is of paramount interest to many industrial applications. There are multiple factors, such as fiber aspect ratio and volume fraction, that play a significant role in altering the rheological behavior of suspensions. Three-dimensional (3D) numerical simulations of coupled differential equations of the suspension of fibers are computationally expensive and time-consuming. Machine learning algorithms can be trained on the available data and make predictions for the cases where no numerical data are available. However, some widely used machine learning surrogates, such as neural networks, require a relatively large training dataset to produce accurate predictions. Multi-fidelity models, which combine high-fidelity data from numerical simulations and less expensive lower fidelity data from resources such as simplified constitutive equations, can pave the way for more accurate predictions. Here, we focus on neural networks and the Gaussian processes with two levels of fidelity, i.e., high and low fidelity networks, to predict the steady-state rheological properties, and compare them to the single-fidelity network. High-fidelity data are obtained from direct numerical simulations based on an immersed boundary method to couple the fluid and solid motion. The low-fidelity data are produced by using constitutive equations. Multiple neural networks and the Gaussian process structures are used for the hyperparameter tuning purpose. Results indicate that with the best choice of hyperparameters, both the multi-fidelity Gaussian processes and neural networks are capable of making predictions with a high level of accuracy with neural networks demonstrating marginally better performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Overlapping Domain Coupling of Multidimensional and System Codes in NEAMS - Pronghorn and SAM

This report describes the implementation and testing of domain-overlapping coupling of the Pronghorn and SAM codes. Both Pronghorn and SAM are codes developed by the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program for analysis of advanced nuclear reactors. Pronghorn focuses on analyzing multi-dimensional core flow conditions, while SAM focuses on analysis of the entire system including, among other components, piping, pumps, and heat exchangers. It is desirable for many advanced reactor thermal-hydraulics simulations to couple Pronghorn and SAM to obtain self-consistent solutions in the core and the system. The domain-overlapping coupling approach provides a robust and numerically efficient candidate for this coupling. We describe the algorithm for obtaining self-consistent solutions between a SAM network of flow channels and an overlapped multi-dimensional Pronghorn domain. The algorithm is tested for a variety of simple canonical test problems and a model of the Molten Salt Reactor Experiment (MSRE). Results demonstrate that the domain overlapping coupling approach provides consistent pressure, mass flow rate, enthalpy, and passive scalar distributions between Pronghorn and SAM and converges reliably and efficiently (at most 15 iterations, usually much less than 10) for the considered test problems.

97 MATHEMATICS AND COMPUTING↗