Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermal 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 217 records · Page 12

Evolution of a Geological Model for Co-Producing Electricity at the Blackburn Oil Field, Nevada: Preprint

Oil fields characterized by significant proportions of hot water in the reservoir offer a unique opportunity for geothermal energy co-production and utilization. Development of an accurate geological model is essential in understanding the subsurface. However, this is a challenge at the preliminary stages of resource assessment when data on the subsurface parameters and subterranean fault network is scarce. Our focus in this paper is to explore a methodology for developing an accurate conceptual geological model of the complexly faulted reservoir at the Blackburn Field, Nevada, to create a numerical simulation that will eventually be used to assess the potential of co-producing electricity. There are four steps to the evolution of this model: (1) collection of public data and research on the geological setting of the system; (2) digitization of seismic time-interpretation and cross-sections to create the most representative model of the subsurface; and (3) creation of a natural-state numerical model to simulate the thermal gradient through the subsurface. Utilizing publicly available data of the geological system and Leapfrog conceptual modeling software, we refined a preliminary model to be more representative of the subsurface geology of an active geothermal exploration stage project. The conceptual geological model developed has been imported into VOLSUNG reservoir modeling software to create a natural-state model, further characterizing the Blackburn Field reservoir's thermal gradient. In the future, new seismic data will be utilized to further evolve the conceptual geological model for numerical simulations.

Basin and Range province↗

Extrapolation of thermal conductivity in non-equilibrium molecular dynamics simulations to bulk scale

Predictions of the bulk scale thermal conductivity of solids using non-equilibrium molecular dynamics (MD) simulations have relied on the linear extrapolation of the thermal resistivity versus the reciprocal of the system length in the simulations. Several studies have reported deviation of the extrapolation from linearity near the micro-scale, raising a concern of its applicability to large systems. To investigate this issue, present work conducted extensive MD simulations of silicon with two different potentials (EDIP and Tersoff-II) for unprecedented length scales up to 10.3 μm and simulation times up to 530 ns. For large systems ≥0.35 μm in size the non-linearity of the extrapolation of the reciprocal of the thermal conductivity is mostly due to ignoring the dependence of the thermal conductivity on temperature. To account for such dependence, the present analysis fixes the temperature range for determining the gradient for calculating the thermal conductivity values. However, short systems ≤0.23 μm in size show significant non-linearity in the calculated thermal conductivity values using a temperature window of 500 ± 10 K from the simulations results with the EDIP potential. Since these system sizes are shorter than the mean phonon free path in EDIP (~0.22 μm), the nonlinearity may be attributed to phonon transport. For the MD simulations with the Tersoff-II potential there is no significant non-linearity in the calculated thermal conductivity values for systems ranging in size from 0.05 to 5.4 μm.

42 ENGINEERING↗

Thermal insulation versus capacitance: A simulation experiment comparing effects of shade and hyporheic exchange on daily and seasonal stream temperature cycles

In streams where water temperatures stress native biota, management of riparian shade or hyporheic exchange are both considered viable management strategies for reducing the peaks of daily and seasonal stream channel temperature cycles. Although shade and hyporheic exchange may have similar effects on stream temperatures, their mechanisms differ. Improved understanding of the heat-exchange mechanisms influenced by shade and hyporheic exchange will aid in the appropriate application of either stream temperature management strategy. To illustrate a conceptual model highlighting shade as ‘thermal insulation’ and hyporheic exchange imparting ‘thermal capacitance’ to a stream reach, we conducted an in-silico simulation modelling experiment increasing shade or hyporheic exchange parameters on an idealized, hypothetical stream. We assessed the potential effects of increasing shade or hyporheic exchange on a stream reach using an established process-based heat-energy budget model of stream-atmosphere heat exchange and incorporated an advection-driven hyporheic heat exchange routine. The model tracked heat transport through the hyporheic zone and exchange with the stream channel, while including the effects of hyporheic water age distribution on upwelling hyporheic temperatures. Results showed that shade and hyporheic exchange similarly damped diurnal temperature cycles and differentially altered seasonal cycles of our theoretical stream. In winter, hyporheic exchange warmed simulated channel temperatures whereas shade had little effect. In summer, both shade and hyporheic exchange cooled channel temperatures, though the effects of shade were more pronounced. Our simple-to-grasp analogies of ‘thermal insulation’ for shade effects and ‘thermal capacitance’ for hyporheic exchange effects on stream temperature encourage more accurate conceptualization of complex, dynamic heat exchange processes among the atmosphere, stream channel, and alluvial aquifer.

54 ENVIRONMENTAL SCIENCES↗

Incorporation of Thermal Hydraulic Models for Thermal Power Dispatch into a PWR Power Plant Simulator

This report describes the development, modeling, and results of a generic pressurized water reactor power plant simulator that incorporates coupled electrical and thermal power dispatch to an industrial process located approximately one kilometer from the nuclear power plant. The simulator is a commercial PWR simulator that has been modified to include thermal power dispatch as described in past milestone reports [ , ]. The commercial PWR simulator is a generic simulator available from GSE SYSTEMS® (Sykesville, MD, USA) that is built using RELAP5-HDTM Real-Time Solution and in-house software developed by GSE Systems. This generic PWR (GPWR) simulator performs real-time simulation of the complete power plant from the reactor neutronics to the electricity generation and distribution. All primary, secondary, and auxiliary systems are modeled including all control logic in order to provide the most accurate representation of actual nuclear power plant (NPP) operation, and the simulator results have been rigorously verified by an actual NPP operating at approximately 1 GWe. This report is a continuation of worked performed in previous years, and supplemental information from previous reports is included in the appendix for reference.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CFD/FEA Co-Simulation Framework for Analysis of the Thermal Barrier Coating Design and Its Impact on the HD Diesel Engine Performance

Thermal barrier coatings (TBCs) have been investigated both experimentally and through simulation for mixing controlled combustion (MCC) concepts as a method for reducing heat transfer losses and increasing cycle efficiency, but it is still a very active research area. Early studies were inconclusive, with different groups discovering obstacles to realizing the theoretical potential. Nuanced papers have shown that coating material properties, thickness, microstructure, and surface morphology/roughness all can impact the efficacy of the thermal barrier coating and must be accounted for. Adding to the complexities, a strong spatial and temporal heat flux inhomogeneity exists for mixing controlled combustion (diesel) imposed onto the surfaces from the impinging flame jets. In support of the United States Department of Energy SuperTruck II program goal to achieve 55% brake thermal efficiency on a heavy-duty diesel engines, this study sought to develop a deeper insight into the inhomogeneous heat flux from mixing controlled combustion on thermal barrier coatings and to infer concrete guidance for designing coatings. To that end, a co-simulation approach was developed that couples high-fidelity computational fluid dynamics (CFD) modeling of in-cylinder processes and combustion, and finite element analysis (FEA) modeling of the thermal barrier-coated and metal engine components to resolve spatial and temporal thermal boundary conditions. The models interface at the surface of the combustion chamber; FEA modeling predicts the spatially resolved surface temperature profile, while CFD develops insights into the effect of the thermal barrier coating on the combustion process and the boundary conditions on the gas side. The paper demonstrates the capability of the framework to estimate cycle impacts of the temperature swing at the surface, as well as identify critical locations on the piston/thermal barrier coating that exhibit the highest charge temperature and highest heat fluxes. In addition, the FEA results include predictions of thermal stresses, thus enabling insight into factors affecting coating durability. An example of the capability of the framework is provided to illustrate its use for investigating novel coatings and provide deeper insights to guide future coating design.

42 ENGINEERING↗

A time-parallel method for scalable heat transfer simulations of additive manufacturing

Here, a major challenge in simulating the thermal behavior in additive manufacturing processes is the disparate length and time scales between transport phenomena occurring in the melt pool and the component. A common simulation approach relies on spatial decomposition for parallel computing, but due to the nature of heat transfer in AM, where most of the computational expenditure is localized near the melt pool, the computational speedup from spatial parallelization saturates quickly. Therefore, additional parallelism by means of time-domain decomposition is needed to fully take advantage of high-performance computing (HPC) resources. This work introduces a time-parallel method to improve the computational scalability of additive manufacturing simulations on HPC systems, while maintaining high temporal resolution of heat transfer near the melt pool. The method, inspired by the nonlinear paraexp formalism, performs an iterative superposition of nonlinear solutions to the initial value problem, integrating the heat equation across overlapping time-parallel intervals. For a single layer of the NIST AMB2018–01 L7 benchmark problem, the method achieves a 38.51x speedup in wall-clock time with a maximum error in the global temperature solution of 0.99%. This reduces the total solution time from 196.72 min to 5.11 min on 128 nodes of the ORNL Frontier supercomputer. The tradeoff between accuracy and total wall-clock time is investigated and recommendations for time-parallel deployment for AM problems are made.

Additive manufacturing↗

Nonnuclear Experimental Capabilities to Support Design, Development, and Demonstration of Microreactors

This work provides a summary of selected experimental capabilities being developed to support nonnuclear testing and demonstration of technology in support of microreactors under the U.S. Department of Energy’s (DOE’s) Microreactor Program. Major capabilities include the Single Primary Heat Extraction and Removal Emulator (SPHERE) and the Microreactor Agile Non-nuclear Experimental Test Bed (MAGNET). The SPHERE facility allows for controlled testing of the steady-state and transient heat rejection capabilities of a single heat pipe using electrical heaters that simulate nuclear heating. The facility is capable of monitoring axial temperature profiles along the heat pipe and surrounding test articles during startup, steady-state operation, and transients. Instrumentation includes noncontact infrared thermal imaging, surface thermocouples, spatially distributed fiber optic temperature and strain sensors, electrical power meters, and a water-cooled, gas-gap calorimeter for quantifying heat rejection from the heat pipe. The facility can be operated under both vacuum and inert-gas conditions. The MAGNET facility is a large-scale, 250-kW electrically heated microreactor test bed to enable nonnuclear experimental evaluation of a variety of microreactor concepts. It can be supplied to electrically heat a scaled section of a microreactor and further test the capabilities of heat rejection systems. The initial MAGNET experiments will support technology maturation and reduce uncertainty and risk associated with the design, operation, and deployment of monolithic heat pipe–based reactors. However, this test bed can broadly be applied to multiple microreactor concepts to evaluate a wide range of thermal-hydraulic and structural phenomena such as interface coupling with power conversion units and other collocated systems. MAGNET can evaluate integral thermomechanical effects during electrical heating of an array of heat pipes in a larger test article. Examples of initial testing will include thermal stresses in the monolith and the impact of debonding of a heat pipe from the core block and how that failure could impact surrounding heat pipes, i.e., understanding the potential for cascading failure. This work also discusses some modeling capabilities that can support experiment design, analysis, and interpretation, including the heat pipe code Sockeye and a comparison of thermal-structural simulations performed using ABAQUS and STAR-CCM+.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Comment on “Brownian motion of droplets induced by thermal noise”

We simulate phase separated fluids using the Cahn-Hillard fluctuating hydrodynamic (CH-FHD) model and measure the statistical properties of capillary waves generated by thermal fluctuations. Our measurements are in good agreement with stochastic lubrication theory and molecular dynamics simulations but differ significantly from recent CH-FHD results by Zhang et al. [Phys. Rev. E 109, 024208 (2024)2470-004510.1103/PhysRevE.109.024208]. Specifically, we find that capillary wave statistics at thermodynamic equilibrium are independent of transport properties, namely viscosity and species diffusion.

Bell, J B↗

Reservoir Characterization and Techno-Economic Analysis of Enhanced Geothermal and Closed-Loop Systems in the Wattenberg Field of the Denver-Julesburg Basin, Colorado

This report provides details of reservoir characterization and techno-economic analysis for enhanced geothermal systems and closed-loop geothermal systems in the Denver-Julesburg Basin Wattenberg area. The analysis contributes to the Geothermal Limitless Approach to Drilling Efficiencies (GLADE) project, a multi-institutional research initiative focused on advancing geothermal energy development by reducing drilling costs and improving the rate of penetration. The GLADE project brings together national laboratories, academic institutions, and industry partners, including the U.S. Department of Energy's Geothermal Technologies Office, Occidental Petroleum, the National Laboratory of the Rockies (NLR), Los Alamos National Laboratory, Colorado School of Mines, Louisiana State University, Texas A&M University, and several drilling technology companies. As part of this effort, NLR developed and applied slender-body theory modeling tools and the GEOPHIRES simulator to support techno-economic analysis. The slender-body theory wellbore simulator evaluated thermal performance in U-loop, Eavor-type multilateral, and enhanced geothermal system well configurations. Results indicated that an optimized single U-loop design delivers an outlet temperature of 168.8 degrees Celsius,stabilizing at a final outlet temperature of 135.3 degrees Celsius. Heat production also improved significantly, with an average thermal output of 16.9 MWth and a final thermal output of 11.0 MWth. These results provided critical baseline data for selecting optimized systems for techno-economic analysis. Building on these outputs, NLR applied the GEOPHIRES techno-economic simulator to estimate the levelized cost of electricity for both the Eavor multilateral and enhanced geothermal systems. The sensitivity analysis showed that drilling costs exert a particularly strong impact on the levelized cost of electricity of multilateral systems. Therefore, advances in drilling technology and improved drilling rates could yield significant gains in the economic performance of these systems. In this context, the GLADE project directly supports the U.S. Department of Energy's mission to reduce the cost of geothermal energy and accelerate the deployment of next-generation geothermal technologies.

15 GEOTHERMAL ENERGY↗

Reduced-Order Modeling Method for Phase-Change Thermal Energy Storage Heat Exchangers

Thermal energy storage can facilitate the effective utilization of renewable energy. To speed up the design process of thermal energy storage devices, it is critical to develop fast and accurate modeling methods for phase change material embedded heat exchangers (PCM HXs). This study developed and compared two approximation-assisted reduced-order PCM HX models for the simulation of thermal storage components and systems, which were verified against a validated finite-volume model. They are a pure black-box model, and a grey-box model based on the Number of Transfer Units (or effectiveness-NTU) approach. We used the reduced-order models to predict the performance of standalone PCM HX for 1000 cases with different fluid inlet conditions and HX designs. We also integrated the models with a vapor compression system to predict the compressor energy consumption and total charging time under various conditions. The results show that overall, the black-box model gave more accurate results than the grey-box model. On average, the mean absolute deviation in the PCM HX fluid outlet temperature was 0.05 K and 0.1 K for the black-box model and grey-box model, respectively. The grey-box model gave larger temperature deviation toward the end of the phase change process, due to the simplified two-node PCM representation. On the system level, the mean absolute deviation in compressor energy consumption was 0.2% and 0.3% for the black-box model and grey-box model, respectively. Their mean absolute deviation in total charging time was 1.1% and 2.6%, respectively. In terms of computation efficiency, the system simulation speedup ratio gained by using the reduced-order models was 11 to 57. System simulation time was decreased from an average of 1465 s to an average of 59 s. This shows that the proposed reduced-order modeling methods can be used to predict system performance of PCM HXs with less than 3% accuracy penalty, and 25 times less computational time than finite-volume models, enabling faster design and evaluation of PCM thermal storage devices.

25 ENERGY STORAGE↗

COUPLED GRIFFIN AND PRONGHORN SIMULATION OF THE MOLTEN SALT FAST REACTOR (MSFR) FOR THE VIRTUAL TEST BED

The Virtual Test Bed (VTB) repository hosts a wide range of challenge problems to showcase modeling and simulation capabilities to support advanced reactor demonstration. An overview of a coupled multiphysics model for Molten Salt Reactors (MSR) is presented here. The analysis leverages MOOSE-based tools (notably Griffin and Pronghorn) for neutronic and thermal hydraulic simulation. Neutron precursor drifting capability is showcased, along with some initial turbulence models. Both steady-state and transient coupled multiphysics results for an MSR concept are discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Evaluation of Additively Manufactured Monolithic SiC and SiC Ceramic Matrix Composites for CSP- Fabrication and Testing of Receiver Design Feature Specimens in a Simulating Lab Test via Laser Heating

Critical to the creation of a multi-component SiC receiver design concept capable of meeting the SETO CST cost goal objectives of <$150/kWth is the optimal utilization of material for solar absorption and heat transfer to transport fluid. This is to maximize the solar-thermal efficiency in the techno-economic analysis of derived use case. As a first step in preparation of sub-component design for on-sun test evaluations we use different proposed SiC absorber element in a customized laser heat flux test to assess an effective heat transfer coefficient characteristic of the design. The test results are analyzed to extract the effective heat transfer coefficient, thus enable optimization of the geometry of receiver components for heat transfer, efficiency attributes and the thermo-structural management of the components. High intensity CO2 laser is used in combination with air flow through the test specimen to assess the effective heat transfer coefficient under heat flux conditions corresponding to high solar concentrations in the range of 1000-2000 suns.

14 SOLAR ENERGY↗

Influence of extreme fracture flow channels on the thermal performance of open-loop geothermal systems at commercial scale

Abstract Adequate stewardship of geothermal resources requires accurate forecasting of long-term thermal performance. In enhanced geothermal systems and other fracture-dominated reservoirs, predictive models commonly assume constant-aperture fractures, although spatial variations in aperture can greatly affect reservoir permeability, fluid flow distribution, and heat transport. Whereas previous authors have investigated the effects of theoretical random aperture distributions on thermal performance, here we further explore the influence of permeability heterogeneity considering field-constrained aperture distributions from a meso-scale field site in northern New York, USA. Using numerical models of coupled fluid flow and heat transport, we conduct thermal–hydraulic simulations for a hypothetical reservoir consisting of a relatively impervious porous matrix and a single, horizontal fracture. Our results indicate that in highly channelized fields, most well design configurations and operating conditions result in extreme rates of thermal drawdown (e.g., 50% drop in production well temperatures in under 2 years). However, some other scenarios that account for the risks of short-circuiting can potentially enhance heat extraction when mass flow rate is not excessively high, and the direction of geothermal extraction is not aligned with the most permeable features in the reservoir. Through a parametric approach, we illustrate that well separation distance and relative positioning play a major role in the long-term performance of highly channelized fields, and both can be used to help mitigate premature thermal breakthrough.

Rangel-Jurado, Nicolás↗

Effects of Galactic Irradiation on Thermal and Electronic Transport in Tungsten

The impact of irradiation on the thermal and electronic properties of materials is a persistent puzzle, particularly defect formation at the atomic and nanoscales. This work examines the nanoscale effects of low-energy irradiation on tungsten (W), focusing on defect-induced modifications to thermal and electronic transport. Using the site-projected thermal conductivity (SPTC) method [A. Gautam et al. PSS-RRL, 2400306, 2024], bulk and twin-grain boundary W with vacancy defects based on the Norgett–Robinson–Torrens displacements per atom model is analyzed. SPTC provides a detailed prediction of postcascade spatial thermal conductivity distribution. Electronic conductivity activity is estimated using the “N 2 method” [K. Nepal et al. Carbon, 119711, 2025] to explore the consequences of vacancies and grain boundaries, highlighting the defect-dependent nature of charge transport behavior. These findings offer high-resolution insights into irradiation-driven transport phenomena, with implications for space-exposed materials and nanoscale thermal/electronic management.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

RaDIATE Collaboration Thermal shock studies

As next-generation accelerator target facilities (High Energy Physics, Spallation Sources, ... ) become increasingly more powerful and intense, high power target systems face key technical challenges, such as radiation damage and thermal shock. Those ultimately degrade the performance and lifetime of targets and have been identified as the leading cross-cutting challenges of high-power target facilities. In order to operate reliable beam-intercepting devices in the framework of energy and intensity increase for next generation accelerators, the RaDIATE Collaboration (Radiation Damage In Accelerator Target Environment), established in 2012 and managed by Fermilab, brings together existing expertise in nuclear material and accelerator targets from 20 international institutions, including CERN, to execute a coordinated strategy for high power targetry R&D. In this context, several thermal shock studies were performed at CERN's HiRadMat (High-Radiation to Materials) facility, that took a key step towards improving our knowledge on target damage tolerance. The first experiment HRMT-24, supported by EURCARD2 and completed in 2015, successfully validated the Johnson-Cook strength model developed at SwRI on beryllium S200FH (used for beam window material), providing a better confidence in simulating the thermal shock response of current and future S200FH beryllium components. HRMT-43, supported by ARIES and completed in 2018, tested various materials (Be, C, SiC, Si, Ti and ceramic nanofiber). It was a first and unique test which included pre-irradiated specimens from high energy proton beam irradiation to identify thermal shock response differences between non-irradiated and previously irradiated materials. Real-time measurement of dynamic thermomechanical response of graphite helped to benchmark numerical simulations.

43 PARTICLE ACCELERATORS↗

Integrated Steady-State System Package for Nuclear Thermal Propulsion Analysis Using Multi-Dimensional Thermal Hydraulics and Dimensionless Turbopump Treatment

Nuclear thermal propulsion is an evolving technology that can be utilized for long-distance space travel. This technology yields the advantage of a high thrust and specific impulse, but requires an examination of the potential design adjustments necessary to enhance its feasibility. The development of nuclear thermal propulsion requires a comprehensive understanding of the system-level behavior during transient and steady-state operation. This paper extends our previous research by including the proper handling of turbomachinery with multi-channel thermal hydraulic simulations only for steady-state solutions. The system-level approach presented here enables the treatment of the turbopump components through non-dimensional analysis that eliminates the assumption of constant efficiencies. All the other components within the system (e.g., reflector and core) can be discretized to multiple channels and layers, in which the full thermal hydraulic solution is established. The approach chosen here enables the realistic modeling of the propellant flow within the expander cycle by capturing the pressure losses, mass flow rate splits, and enthalpy gain for various operational conditions. The verification of the package is completed through point comparisons of previous investigations into similar system designs. Furthermore, sensitivity studies are used to benchmark the capabilities of the package and investigate solution variations due to the perturbation of operational conditions and regimes. The sensitivity studies performed here are important to capture variation in flow characteristics (e.g., temperature, pressure, mass flow rates) for different design objectives such as the thrust and specific impulse. This work demonstrates that system-level simulations lacking multi-channel capability and proper turbomachinery treatment may yield higher uncertainties in understanding the engine’s response and characteristics to changing various requirements. This is extremely important when screening the design space of such propulsion systems and when transient simulations are required.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Design and Study of Inductively Coupled Plasma Chamber Components Using the SupRISE Test Device at DIII-D

The DIII-D National Fusion Facility aims to increase the auxiliary heating power for the tokamak by upgrading the Neutral Beam Injection (NBI) system. In collaboration with North Carolina State University, the conventional ‘arc and filament’ NBI ion sources will be converted to inductively coupled plasma (ICP) sources which utilize radio frequency (RF) coupling to maximize reliability for high power operation. In support of this initiative, a full scale test device (SupRISE, Superior Radiofrequency Ion Source Experiment) is currently under construction at the DIII-D Facility. In preparation for the construction of a full scale prototype that can be installed on the DIII-D NBI system, experiments on SupRISE will be conducted to determine the optimal RF frequency for high power coupling, Faraday shield slit configuration, and ICP chamber length. SupRISE is comprised of an approximately 30 x 70 cm quartz dielectric vessel with an internal Faraday shield enclosed in a secondary vacuum chamber to ensure structural stability of the dielectric. Actively cooled front and back plates are designed to reduce the thermal stresses on the plasma facing components and mate with the existing accelerator used by the NBI system at DIII-D. 50 kW of RF power will be coupled to the plasma through the quartz over a variable frequency range of 4-8 MHz to sustain a plasma density of ∼ 10 18 m −3 for a 10 s ON, 210 s OFF duty cycle. Various modelling efforts have been employed to simulate the thermal and stress profiles over the primary components of the SupRISE device as well as the inductance behavior of the RF antenna. These simulation results and the final design for SupRISE will be presented. An additional reduced-scale predecessor ICP source (called RISE) has been used to develop a predictive match model that will be applied to frequency optimization studies on SupRISE. Furthermore, the outcomes of this research and complementary efforts at North Carolina State University are essential for the incorporation of ICP NBI positive ion sources at the DIII-D facility.

DIII-D↗

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↗