SEARCH · Engineering Papers

Results for “THERMAL HYDRAULICS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Investigation of Gas Turbine Internal Cooling Using Supercritical CO 2 —Effect of Surface Roughness and Channel Aspect Ratio

In this paper, an experimental and numerical investigation of internal cooling channels with rib turbulators is presented with sCO 2 as the working fluid at process conditions (pressure-20.7 MPa and temperature up to 150 °C). The effect of channel aspect ratio up to 2:1 on thermal-hydraulic performance is explored in additively manufactured rectangular channels and square channels, both with and without 60 deg ribs on the top and bottom sides. The Wilson-plot method is employed to experimentally measure channel-averaged Nusselt number over a Reynolds number range up to 370,000. The friction factor is calculated from pressure drop and mass flow rate and additionally, the overall thermal performance factor (TPF) is reported. A companion computational fluid dynamics (CFD) simulation is performed for the rib turbulated cooling configurations reported in the experiments using the Reynolds average Navier–Stokes-based turbulence model. The objective of the numerical study is to gain insight into the local heat transfer augmentation in the ribbed channels as a result of varying the aspect ratio, channel configuration (square versus rectangular), operating conditions (Reynolds number) and the surface roughness, an inherent outcome of the additive manufacturing process. Surface roughness is simulated using sand grain roughness height (K S ) calculated from the experimental data, and a comparison is presented with the corresponding channel configuration with varying surface roughness heights starting from smooth surfaces (K S = 0). Experimental results indicate that the heat transfer augmentation is negligible in the rectangular channels with ribs on the long side compared to the square channel. However, it is enhanced by 60% in comparison to placing ribs on the shorter side. The TPF remains constant at around 1 for the entire range of Reynolds numbers consistent with prior work at the National Energy Technology Laboratory (NETL). Here, the simulation results highlight that increased surface roughness can have a favorable considerable influence on Nusselt number and overall thermal performance enhancement.

42 ENGINEERING↗

DOE OSTI2026

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗

NASA NTRS2007

High Fidelity, Fuel-Like Thermal Simulators for Non-Nuclear Testing: Analysis and Initial Test Results

Non-nuclear testing can be a valuable tool in the development of a space nuclear power system, providing system characterization data and allowing one to work through various fabrication, assembly and integration issues without the cost and time associated with a full ground nuclear test. In a non-nuclear test bed, electric heaters are used to simulate the heat from nuclear fuel. Testing with non-optimized heater elements allows one to assess thermal, heat transfer. and stress related attributes of a given system, but fails to demonstrate the dynamic response that would be present in an integrated, fueled reactor system. High fidelity thermal simulators that match both the static and the dynamic fuel pin performance that would be observed in an operating, fueled nuclear reactor can vastly increase the value of non-nuclear test results. With optimized simulators, the integration of thermal hydraulic hardware tests with simulated neutronic response provides a bridge between electrically heated testing and fueled nuclear testing. By implementing a neutronic response model to simulate the dynamic response that would be expected in a fueled reactor system, one can better understand system integration issues, characterize integrated system response times and response characteristics and assess potential design improvements at relatively small fiscal investment. Initial conceptual thermal simulator designs are determined by simple one-dimensional analysis at a single axial location and at steady state conditions; feasible concepts are then input into a detailed three-dimensional model for comparison to expected fuel pin performance. Static and dynamic fuel pin performance for a proposed reactor design is determined using SINDA/FLUINT thermal analysis software, and comparison is made between the expected nuclear performance and the performance of conceptual thermal simulator designs. Through a series of iterative analyses, a conceptual high fidelity design is developed: this is followed by engineering design, fabrication, and testing to validate the overall design process. Test results presented in this paper correspond to a "first cut" simulator design for a potential liquid metal (NaK) cooled reactor design that could be applied for Lunar surface power. Proposed refinements to this simulator design are also presented.

Bragg-Sitton, Shannon M.↗

DOE OSTI2026

Ultra-low thermal resistance and pressure drop copper and copper-tungsten diamond-shaped pin fin cold plates for liquid cooling of electronics

Modern and future data centers face increasing cooling challenges due to increasing chip thermal design power and die size, along with the need to reduce energy consumption used for cooling. High performance cooling solutions that maintain a low chip junction temperature are needed to ensure electronics reliability. This work develops an ultra-low thermal resistance and low pressure drop 75 mm × 75 mm cold plate, intended for next-generation electronics cooling. The cold plate features an array of diamond-shaped pin fins and integrated copper tungsten heat spreader, selected for its low coefficient of thermal expansion which reduces thermomechanical deformation and allows for closer integration of the cold plate with silicon dies. Starting with 300 candidate designs, three-dimensional computational fluid dynamics simulations predict the thermal-hydraulic performance of cold plate subsections. The highest performing geometries are evaluated with high fidelity simulations. Four cold plates are manufactured for experiments: three with diamond-shaped pin fins and one with straights fins for comparison purposes. The cold plates are fabricated from copper-tungsten (CuW), copper (Cu), or aluminum-silicon-magnesium alloy (AlSi10Mg). The diamond-shaped pin fins achieve a roughly 15 % lower thermal resistance compared to the conventional straight fin microchannel. The highest performing design achieves a chip-to-coolant (including thermal interface material) thermal resistance of 9.0 K/kW in CuW and 6.9 K/kW in Cu under a 1 kW heat load with an inlet-to-outlet pressure drop of 9.0 kPa and water as the working fluid. This work demonstrates ultra-low thermal resistance and pressure drop cold plates for large die, high heat load applications, and shows that CuW is an attractive cold plate material for improved reliability in next generation data center cooling.

Coefficient of thermal expansion↗

DOE OSTI2024

Performance evaluation of finned tube heat exchanger using curved wavy delta winglet vortex generators with circular perforations

Vortex generation is recognized as an effective passive approach to improve the heat transfer rate in fin and tube heat exchangers (FTHEs). The current study proposed innovative designs of curved wavy delta winglet vortex generators (CWDWVGs), both without and with circular perforations, to improve the heat transfer efficiency of FTHEs. There is potential to increase heat transfer performance further through various CWDWVG designs. Here, this study explores seven unique CWDWVG configurations, from 1-wave to 7-wave. A 3-D computational numerical model is utilized to evaluate the Thermo-hydraulic performance of FTHEs fitted with these different CWDWVG configurations across Reynolds numbers from 400 to 2000. This comparative analysis of the Thermo-hydraulic performance of FTHEs featuring four parallel circular tube layouts assesses configurations both with and without vortex generators (VGs) and various hole configurations. The evaluation of Thermo-hydraulic performance involves different parameters, including the London area goodness factor (LAGF), Colburn factor (j), friction factor (f), pressure drop (?P), and Nusselt number (Nu). Results demonstrate that the various CWDWVG configurations and the number of holes in them substantially affect the efficiency, as evaluated by the dimensionless Performance Evaluation Criteria (PEC). Notably, the 7-wave CWDWVGs surpassed other configurations, and integrating circular punched perforations further improved the thermal-hydraulic performance of FTHE. Specifically, the 7-wave CWDWVGs without holes demonstrated superior performance over other configurations, showing a significant increase in Nusselt number by 75.18% and 85.16% at Reynolds numbers of 2000 and 400, respectively, alongside an increase in pressure drop by 216.38% to 224.96%. Meanwhile, the 7-wave CWDWVGs with eight holes, in comparison to those without holes, exhibited a Nusselt number increase of 0.85%, a pressure drop decrease of 7.31%, and a reduction in the friction factor by 5.82%.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

Assessment of Thermal Stratification in Versatile Test Reactor Transients

The Versatile Test Reactor (VTR) is a fast-spectrum test reactor currently being developed in the United States under the direction of the US Department of Energy. The conceptual design of the 300 MWth pool-type sodium-cooled fast reactor (SFR) has been led by the US National Laboratories in collaboration with General Electric-Hitachi and Bechtel National Inc. Safety performance analysis for the VTR conceptual design is being performed with the systems thermal-hydraulics (SYSTH) module of the SAS4A/SASSYS-1 liquid-metal reactor safety analysis code system. Since the current model of the VTR employs a simple perfect mixing model for large plena like the hot pool, it is not able to predict temperature variations that may develop during the transient. Prior work simulating the response of SFRs to postulated events like the Protected Station Blackout (PSBO) has shown that the phenomenon of thermal stratification, where stable thermal layers accumulate in the hot pool, may delay the transition to natural circulation and thus impact the predicted transient progression. Thus, an effort has begun to model this transient by integrating a Computational Fluid Dynamics (CFD) model of the hot pool into the SAS4A/SASSYS-1 model of the Primary Heat Transport System during the simulation of the PSBO event. A three-dimensional Volume-Of-Fluid CFD model of the VTR hot pool has been developed for the co-simulation of SAS4A/SASSYS-1 with CFD. In this work, the standalone SAS4A/SASSYS-1 calculation and the standalone CFD calculation based on the SAS4A/SASSYS-1 calculation result were produced on the Idaho National Laboratory High Performance Computing cluster, SAWTOOTH. At this time, only the standalone CFD and SAS4A/SASYS-1 simulation results are provided.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2021

Massachusetts Institute of Technology Reactor LEU Fuel Element Flow Test Conceptual Design – Hydraulic Reactor Design Parameters

The Massachusetts Institute of Technology Reactor (MITR-II, also referred to as MITR) is one of six U.S. high performance research reactors (USHPRR), including one critical facility, that is actively collaborating with the U.S. National Nuclear Security Administration (NNSA) Material Management and Minimization (M 3 ) Reactor Conversion Program to convert to the use of low-enriched uranium (LEU, < 20 wt% 235 U) fuel. The MIT Nuclear Reactor Laboratory has been working with the USHPRR Reactor Conversion (RC) Pillar at Argonne National Laboratory to perform fuel element design and fuel cycle performance analyses, steady-state thermal hydraulics safety analyses, and accident safety analyses in preparation for the conversion of MITR and support a preliminary Safety Analysis Report (SAR) for conversion to LEU fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Application of NEAMS Codes to Capture MSR Phenomena

This report documents work completed in FY21 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program’s Molten Salt Reactor (MSR) Application Drivers activity at Argonne. The common focus was on identifying the modeling and simulation functional requirements for designing, licensing, and operating MSRs and applying those capabilities already developed in NEAMS tools to example problems of interest. The four main parts of this report each focuses on a specific area of simulation physics as it relates to MSR phenomena: fuel evolution, chemistry, computational fluid dynamics (CFD), and systems analysis. In terms of fuel evolution, which includes depletion, decay, on-line separations, and transmutation, the current state of computational capabilities for modeling this behavior in liquid-fueled molten salt reactor is discussed. Some of the functional requirements to accomplish the various applications of MSR fuel depletion modeling are highlighted, followed by a summary of recent approaches and code development activities. The chemistry functional requirements were discussed in the context of several applications of high importance for MSRs, such as corrosion, salt chemistry, and species transport. Each of these types of chemistry modeling have considerable impact on various aspects of reactor applications, including informing on reactor designs, improving operational efficiencies, analyzing safety and reactivity concerns, and estimating the mechanistic source term of the reactor. A brief overview is also provided on code development activities ongoing under NEAMS relevant to chemistry modeling of MSRs. In terms of CFD applications, the state-of-the-art spectral element code, Nek5000 was used to model the fluid dynamics within a full core of the Molten Salt Fast Reactor (MSFR) concept designed as part of the Euratom EVOL project. This concept was selected as the challenge problem because of its similar features to several U.S. industry concepts. The goal was to model some of the fast MSR design challenges, including potential large internal re-circulations, the need of accurate tracking of delayed neutron precursors (DNP), and the lack of relevant thermal-hydraulics models/correlations, etc. Therefore, a series of CFD models were created for the MSFR core cavity using a k – τ model two-equation model for the turbulence. These first full core results demonstrated that any potential recirculation zones could be properly identified with the current NEAMS CFD capabilities. The development of these models will also set the stage for future testing of Nek5000’s functionalities to model other MSR thermal-fluid phenomena. Lastly, the validation of SAM against experimental data from the Molten Salt Reactor Experiment, which started in FY19, continues with the inclusion of modeling the reactivity insertion tests. This involved using SAM and its point kinetics model for flowing fuel salt to recreate the time dependent power changes and response after positive reactivity insertions at the 1, 5, and 8 MWt power levels. Through these exercises, several code modifications were suggested to the SAM development team and accommodated to enable closer agreement with solid technical and physical justifications. These include adding a moderator reactivity feedback coefficient as an available input and modifying the solution approach for the point kinetics model. Additional SAM development suggestions for flowing fuel MSRs include adding the capability to allow the moderator power change proportionally with the reactor power and enabling specification of the power and DNP distributions separately.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Design of a neutral thermal scattering (NeTS) module for hydrogen in light water

The accurate representation of thermal scattering law (TSL) data is integral to the design and characterization of many modern nuclear systems, particularly those using light water as a moderator/coolant. As a material-dependent distribution over energy-momentum phase space, the TSL may exhibit a variety of unique and relevant conditional (e.g., temperature, pressure) and compositional (e.g., porosity, stoichiometry, radiation damage) dependencies. Currently, there are various approaches to incorporate temperature dependence, which is especially important in coupled neutronic-thermal hydraulic simulations. In each approach, there is an inherent tradeoff between memory consumption and accuracy. Some techniques require tens to hundreds of MBs or more, while others fail to reproduce the underlying data to within 10% error over the considered input domain, despite having a reduced storage burden. This work aims to address both sides of the tradeoff simultaneously by implementing a novel deep learning (DL) approach to TSL representation. The neural thermal scattering (NeTS) concept, which is amenable to an arbitrary number of dependencies, is demonstrated via the inclusion of temperature dependence into a highly compact, highly accurate functional form of the multi- variate TSL (i.e., S(α, β, T)) for hydrogen in light water. Resulting storage requirements are on the order of 100 kB, and median and maximum percent deviations are on the order of 0.1% and 1%, respectively. These measures represent a step improvement over previous techniques. Notably, the developed neural network and feature methodology build on those employed in prior work on beryllium oxide. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Design and Analysis of a 250 MW Plate-fuel Reactor for Nuclear Thermal Propulsion

Nuclear thermal propulsion (NTP) system is a viable option for deep space missions considering its high thrust and lightweight. Promising fuel options being considered for the high operating temperature of NTP include the graphite composite fuel and the tungsten cermet fuel. As there remains uncertainties in the fabrication and performance of these two fuel elements, alternative designs using fuel plates and more standard fuels that are based on proven and tested nuclear technologies are being considered. This paper provides a summary of thermal hydraulic and neutronic analyses of the straight fuel plate concept used for a 250 MW NTP core. The summary of technical findings from a series of scoping studies might be useful for the future NTP engine designs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Solid-to-Fluid Radiative Heat Transfer Modeling for System Analysis Module

System Analysis Module (SAM) is a system-level thermal hydraulics code being developed at Argonne National Laboratory for advanced nuclear reactor analysis. In addition to a wide range of interests from the advanced reactor design community, SAM has also been adopted by the United States Nuclear Regulatory Commission's suite of codes purposed for advanced reactor licensing. Nevertheless, the code is still under active development and new capabilities are being added to address various modeling and simulation challenges for advanced reactor analysis. One such phenomenon important to the thermal behavior of some advanced reactor concepts is radiative heat transfer (radHT). Conditions, such as high temperatures and long optical paths, increase the radiative contributions from solids and coolants alike. This paper discusses the development of radHT modeling in SAM and describes the new capabilities provided for thermal analysis. Depending on the geometry and temperatures of the system at hand, as well as the coolant in question, thermal transfer due to thermal radiation will vary dramatically. Therefore, the ability to model variable radiative systems was maintained as a priority during development of SAM radHT modeling. This newly developed simulation feature provides a flexible solid-to-fluid radiative heat transfer framework necessary for SAM to perform accurate analysis for advanced reactor designs. Test cases are also presented and shown to match analytical solutions, which demonstrate the radHT model's efficacy.

Radiative heat transfer↗

DOE OSTI2025

SAM Two-Phase Flow Model Development and Applications for Operational Transients in Advanced Reactors

As advanced nuclear technologies continue to develop, the need for the flexible operation and generation of these advanced reactors becomes necessary to maximize economic potential. As large-scale experiments are not always feasible, modeling and simulations of advanced reactors play a crucial role in design optimization and analysis. The SAM (System Analysis Module) code developed at Argonne National Laboratory is a state-of-the-art system-level thermal-hydraulic code aimed at simulating advanced reactor systems. Recent code developments have implemented two-phase flow modeling using the homogeneous equilibrium model, and a new steam generator component has been developed to utilize the two-phase flow implementation. In addition to verification tests, a load-following simulation was performed to model a realistic load-following transient in a proposed integrated system consisting of a conceptual advanced reactor known as the Advanced Burner Test Reactor (ABTR) and thermal energy storage (TES) tanks. The integrated system model uses two large TES tanks designed for sodium and a model helical coil steam generator to simulate the operational load-following transient. The flow rates of the feedwater and secondary loops are regulated to meet a prescribed steam generator load consistent with the electricity demand over a 24-h period. In conclusion, the results found the ABTR system was able to maintain stable reactor conditions and primary- and secondary-side characteristics over the course of the load-following transient.

Advanced Burner Test Reactor (ABTR)↗

DOE OSTI2024

Design and evaluation of a dilute flow particle-to-air heat exchanger for energy storage applications

The use of inert and redox-active particles for high-temperature energy storage requires the development of components that can efficiently transfer energy to high-pressure working fluids like supercritical carbon dioxide (sCO 2 ). Dilute flow reactors can enable high working fluid outlet temperatures and minimal parasitic losses compared to moving packed bed and fluidized bed reactors. This research uses both computational and experimental methods to explore the design trade-offs and practical challenges of a novel component for transferring energy from dilute flows of hot, reduced metal oxide (MO x ) particles to sCO 2 in tubes. A discretized thermal resistance network model, which accounts for particle hydrodynamics, multi-mode heat transfer, and reaction equilibrium, guides the design of a prototype device. This device is experimentally tested with a surrogate heat transfer fluids and inert particle temperatures up to 400°C and a heat duty exceeding 1 kW. The data are used to validate the thermal hydraulic sub-models, allowing for the simulation of reacting particle scenarios. Under nominal design conditions, the flow rate of reactive particles is predicted to be 30% lower than that of inert particles for the same energy recovered, with over 70% of the stored particle energy transferred to the sCO 2 . Furthermore, these findings can inform the design of more efficient energy recovery reactors for particle-based systems and can be integrated into system-level concentrated solar power models with thermal storage to optimize operating conditions.

14 SOLAR ENERGY↗

NASA NTRS1998

A Report from the Thermal Science Research Center (TSRC)

A vertical flow loop was designed and assembled to determine the local (circumferential and axial) and mean wall temperature distributions for single-phase and two-phase (subcooled and saturated) downward flow in both uniformly-heated and single-side heated vertical channels. Freon-11 was used as the working fluid in order to directly relate and compare the results with a previous experimental campaign which employed this same working fluid. For a given steady-state experiment, the following parameters were held constant: (1) exit pressure, (2) inlet temperature, and (3) mass velocity. For a given configuration of the 2.2 m long cylindrical channel test section, which had a 1.2 m long heated section, the applied heat rate was varied from zero through successive quasi-steady states to a level which corresponded to localized film boiling in the test section. The measurements showed that the boiling curve changes significantly at higher mass velocities with respect to both the circumferential and axial directions. The slope of the boiling curve changes in a non-monotonic fashion with respect to the circumferential directions. The slope of the boiling curve changes in a non-monotonic fashion with respect to the circumferential direction. The measurements point to the existence of a dry-out phenomenon occurring at multiple levels of the applied heat for the single-side heated channel. In comparing the heat transfer for horizontal channel flow with a vertically downward flow, the results show that significantly lower heat transfer occurs in the horizontal flow. However, this trend reverses as both the Reynolds number and the applied heat rate increase. Both the Liu-Winterton and Shah correlations were compared with the experimental data. The Shah correlation predicted the uniformly heated tube data better. When a thermal hydraulic diameter approach was used for the single-side heated case, the data at upstream locations for Z/L less than 0.5 was bounded above by the Liu-Winterton correlation and below by the Shah correlation. At Z/L = 0.5, the Shah correlation bounded the data; and for Z/L greater than 0.5, both correlations overpredicted the data with the Shah correlation being closest to the data. The present results indicate that additional correlational development is needed. In addressing some of the advanced space thermal management objectives concerning accommodating high heat fluxes in non-uniformly heated systems, a large battery of experiments 88 have been completed where local two-dimensional wall temperature variations were measured for both single-phase and two-phase flow in a single-side heated circular tube. As noted above, the results show significant axial and circumferential variations. Accurately accounting for such variations can result in optimized future advanced space, enhanced (high heat flux) thermal management systems.

Boyd, Ronald D.↗

NASA NTRS1992

An improved heat transfer configuration for a solid-core nuclear thermal rocket engine

Interrupted flow, impingement cooling, and axial power distribution are employed to enhance the heat-transfer configuration of a solid-core nuclear thermal rocket engine. Impingement cooling is introduced to increase the local heat-transfer coefficients between the reactor material and the coolants. Increased fuel loading is used at the inlet end of the reactor to enhance heat-transfer capability where the temperature differences are the greatest. A thermal-hydraulics computer program for an unfueled NERVA reactor core is employed to analyze the proposed configuration with attention given to uniform fuel loading, number of channels through the impingement wafers, fuel-element length, mass-flow rate, and wafer gap. The impingement wafer concept (IWC) is shown to have heat-transfer characteristics that are better than those of the NERVA-derived reactor at 2500 K. The IWC concept is argued to be an effective heat-transfer configuration for solid-core nuclear thermal rocket engines.

Clark, John S.↗

DOE OSTI2025

Thermomechanics coupling to Monte Carlo particle transport on unstructured mesh geometries using Cardinal

Geometry deformation due to thermal expansion influences neutron transport in many systems. Studying this phenomenon involves coupling models for neutronics, thermal hydraulics, and solid mechanics. To enable high fidelity modeling of these coupled physics, new capabilities were introduced in Cardinal, coupling OpenMC Monte Carlo particle transport models with MOOSE thermomechanical physics on unstructured moving-mesh geometries. In this work, we present a fully open-source capability leveraging on-the-fly mesh skinning to automatically regenerate OpenMC geometry, which allows multiphysics feedback from temperature, density, and geometry changes. The new capability is verified using an analytic benchmark slab problem, which couples S 2 neutron transport with thermal conduction, convective boundary conditions, Doppler-broadened cross sections, and nonlinear thermal expansion effects along the heated slab. Cardinal reproduces the analytic solutions for the neutron flux, heating, k eff , and temperature with demonstrated convergence in various error terms including mesh resolution and cross section temperature library spacing. For the nominal benchmark conditions and with a fine mesh, maximum relative errors for neutron flux, temperature, and heating are lower than 1%, while errors in integral quantities such as k eff and slab length are within 1 pcm and 48 µm, respectively. This work (i) presents a new numerical approach to thermomechanics coupling with OpenMC models, (ii) is the first (to our knowledge) to utilize a mechanical partial differential equation (PDE) solution to solve the (Griesheimer and Kooreman, 2022) analytic benchmark, and (iii) develops this verified capability within an open-source package.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Physics analysis and design of heavy water reflected thermal test reactor

Here, this work investigates the option of modifying the Advanced Test Reactor by replacing the current beryllium reflector with heavy water. Such a change may provide some potential benefits for not only increasing the thermal irradiation capabilities but also resolving other problems such as reflector integrity issues due to fast fluence damage, which is always a limiting factor in the lifetime of the current beryllium reflector. This paper presents the analysis and estimation of the ATR core physics parameters by replacing the current beryllium reflector with heavy water (D 2 O). The paper first describes the details of two selected conceptual designs, which are partially reflected with either beryllium or graphite, and how they are derived from the baseline beryllium reflector concept. Then, reactor physics performance parameters for the two new concepts are assessed by comparing with those of the baseline concept. The performance parameters considered in this paper include in-pile tube neutron and gamma fluxes and heating rates, maximum loop voiding reactivity, core power behavior with different power splits, predicted cycle length with a given fuel loading, and thermal hydraulic analysis with a higher lobe power split. It is important to note that this study focuses on the reactor physics aspects and does not delve into the engineering challenges associated with such a design modification.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Design of a prototypical natural circulation water-based reactor cavity cooling system (RCCS) for a pebble-bed generic FHR

A prototypical natural circulation water-based reactor cavity cooling system (RCCS) for the UC Berkeley Mark-1 pebble-bed gFHR is designed based on one-dimensional thermal hydraulics modeling and optimization implemented in an in-house MATLAB code. The model employs a lumped core consisting of fuel pebble and graphite pebble regions, while the graphite reflector, vessel, and RCCS are represented as separate but energy-coupled regions. The model is derived based on steady state energy balance equations accounting for conductive heat transfer from the lumped core to the vessel and subsequent radiative heat transfer to the RCCS and convective heat transfer in the water. Further, mass flow rate is calculated based on momentum and thermal energy balance in the RCCS. Starting with ANL’s water NSTF as baseline, effects of RCCS design parameters such as size of the plate between pipes, pipe diameter, source-to-sink distance, surface emissivity, and pitch from the core are examined in an integral effects framework which accounts for subsequent variation of the mass flow rate, number of pipes, convection coefficient, mean region and surface temperatures, radiative view factors, and water outlet temperature. The importance of including a prototypical reactor design in RCCS design calculations is emphasized as the calculations show that temperature drops from the core to the RCCS fluid to transfer the heat tend to be the limitation not the heat removal capacity of the fluid. A closed-loop design is obtained with capability to safely remove up to 0.72% of the nominal reactor power at maximum estimated peak conditions and 0.40% at shutdown. The design is based on physics calculations and does not account for economic optimization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗