Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat conduction modeling”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Weak-form latent space dynamics identification

Recent work in data-driven modeling has demonstrated that a weak formulation of model equations enhances the noise robustness of a wide range of computational methods. In this paper, we demonstrate the power of the weak form to enhance the LaSDI (Latent Space Dynamics Identification) algorithm, a recently developed data-driven reduced order modeling technique. We introduce a weak form-based version WLaSDI (Weak-form Latent Space Dynamics Identification). WLaSDI first compresses data, then projects onto the test functions and learns the local latent space models. Notably, WLaSDI demonstrates significantly enhanced robustness to noise. With WLaSDI, the local latent space is obtained using weak-form equation learning techniques. Compared to the standard sparse identification of nonlinear dynamics (SINDy) used in LaSDI, the variance reduction of the weak form guarantees a robust and precise latent space recovery, hence allowing for a fast, robust, and accurate simulation. We demonstrate the efficacy of WLaSDI vs. LaSDI on several common benchmark examples including viscid and inviscid Burgers', radial advection, and heat conduction. For instance, in the case of 1D inviscid Burgers' simulations with the addition of up to 100% Gaussian white noise, the relative error remains consistently below 6% for WLaSDI, while it can exceed 10,000% for LaSDI. Similarly, for radial advection simulations, the relative errors stay below 15% for WLaSDI, in stark contrast to the potential errors of up to 10,000% with LaSDI. Moreover, speedups of several orders of magnitude can be obtained with WLaSDI. For example applying WLaSDI to 1D Burgers' yields a 140X speedup compared to the corresponding full order model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Analysis of HolosGen Sub-Scale Simulator with Plant Dynamics Code

The Subcritical Power Module Sub-scale Simulator (SPM-SS) has been designed and constructed by HolosGen LLC under the ARPA-E MEITNER program to simulate the thermal-hydraulic and heat transfer behavior of the full-scale Holos-Quad Subcritical Power Modules (SPMs). Four coupled SPMs, each rated at 5.5MW, form the Holos-Quad gas-cooled microreactor design. The SPM-SS represents a substantially scaled-down system with a power rating less than 40 kW, equipped with an electrically heated fuel cartridge heat exchanger, an electrically heated compressor heat exchanger, and a valve actuated turbine heat exchanger, in addition to a recuperator and a cooler heat exchanger. The fuel cartridge represents a portion of the full-scale SPM core, the compressor heat exchanger mimics the temperature changes resulting from the compressor’s turbomachinery inefficiencies, the turbine heat exchanger mimics the expansion process normally occurring through the turbine, while the recuperator and cooler heat exchangers complete the subscale simulator loop. The heaters equipping the fuel cartridge and the compressor heat exchangers are electronically controlled to simulate normal and off-normal SPM operating conditions. The full-scale Holos-Quad SPM design eliminates the traditional balance of plant and executes thermal-to-electric energy conversion by means of an intercooled Brayton cycle with decoupled compressor-turbine turbomachinery. The Holos-Quad full-scale design is equipped with a multi-stage axial Low- and High-Pressure compressor, and a multistage axial turbine. The SPM-SS is designed for testing and validation of selected components which are instead coupled by a traditional balance of plant. The SPMSS is not equipped with turbo-machinery (compressor and turbine) as the development of these components were excluded from the scope of work under the ARPA-E MEITNER funding program. The SPM-SS balance of plant enables modifications, replacement and testing of individual components with different working fluids and is designed to include the turbo-machinery components that will be developed in future research . The SPM-SS can be operated with different gases, variable mass-flow-rates, pressures, and temperatures to obtain test data for selected components, whose performance can be scaled to validate the computer model of the full-scale SPM at various conditions (e.g., start-up, transients conditions). The SPM-SS can operate at the maximum Holos-Quad design pressure of 7 MPa, and a maximum temperature limited to 650 °C by the electrical heaters. Several SPM-SS tests have been conducted and analyzed with the Plant Dynamics Code (PDC) developed at the Argonne National Laboratory (ANL). These tests aimed at validating the PDC modeled predictions of the full-scale Holos-Quad design with data from selected SPM-SS components. In order to address SPM-SS specific characteristics, such as components heat losses and absence of turbomachinery components, some modifications to the PDC have been implemented to factor the design differences from the full-scale Holos-Quad SPM to the SPM-SS. As the PDC offers capabilities to analyze systems with different working fluids, air, nitrogen, and helium were utilized as the SPM-SS working fluids. Air was utilized to fine-tune the SPM-SS Systems Structures and Components (SSCs), nitrogen was utilized to pressure test the SPM-SS loop at the SPMs design maximum pressure of 7MPa. Helium was utilized as the working fluid circulating through the SPM-SS SSCs for specific tests to validate the PDC predictions of the fuel cartridge heat exchanger. SPM-SS tests data were also analyzed with both the steady-state and transient analysis capabilities offered by the PDC. This report describes the PDC analysis of the SPM-SS tests data, including the necessary code modifications and comparison of the code results with the experimental data. Based on the results, a discussion is presented on how the analysis supports design and transient calculations of the full-scale Holos-Quad microreactor. Also based on the results of this work, recommendations are made for future optimizations of the SPM-SS components and PDC model development needs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Drift kinetic electrostatic simulations of the edge localized mode heat pulse

In the present work, electrostatic drift kinetic simulations of parallel plasma transport within the tokamak scrape-off layer (SOL) are conducted using the COGENT code. The SOL configuration is represented in one-dimensional slab geometry, incorporating a heat source localized in the midplane. The heat source parameters correspond to those characterizing edge-localized modes observed in the Joint European Torus (JET) tokamak. The numerical model includes kinetic treatment of both ions and electrons, a simplified model for the gyrokinetic Poisson equation that allows one to step over short time scales associated with fast electrostatic shear Alfvèn waves, and the logical sheath boundary condition (LSBC) that enforces global system quasineutrality. A third-order accurate LSBC is derived to be consistent with the third-order accurate upwind advection scheme utilized in the code, and it was shown to noticeably impact the simulation results, especially parallel heat flux at the target plate. The findings of this study are in agreement with results from preceding fluid and kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Noble Gas Transport in the MSRE

This study explores the relationship between the physicochemical parameter known as Henry’s gas constant and gas transport across a two-layer film interface. The investigation utilized the Gibbs free energy, incorporating surface and volume terms to elucidate trends in enthalpy and entropy. Notably, our findings align with experimental data and offer predictive insights into the Henry’s gas constants for helium and krypton, which hold significance for future experiments and theoretical developments. Furthermore, this study enhances the Gibbs free energy theory pertaining to the liquid–gas interface. It underscores the substantial contribution of noble gases in this region to volumetric energy as temperature increases. Additionally, we employed Monte Carlo simulations to analyze the effective thermal neutron multiplication factor, denoted as k eff . Our analysis reveals a linear correlation between graphite density and uniform density as a function of temperature. For the 1D Molten Salt Reactor Experiment (MSRE) system, we employed the Mole code to conduct heat and mass transfer calculations. These computations enable us to ascertain the distribution of fuel temperature based on coefficients and thermal properties. We also studied delayed neutron precursors during fuel cycling, taking into account the drift of cycling fuel through Mole–Griffin coupling. Our model represents k eff and β eff across various volume flow rates and salinity compositions. Finally, this study leveraged xenon-135 for continuous on-line monitoring of fuel salts to investigate the impact of steady-state xenon-135 on the MSRE and to better understand its distribution. These efforts build upon previous research related to removal processes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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 OSTI2023

United States cool surfaces deployment plan

Solar-reflective building envelope surfaces, such as cool roofs and walls, can be especially helpful in disadvantaged communities that often have poorly insulated older homes, aging or absent air conditioning units, steep utility bills, polluted air, and high vulnerability and exposure to extreme heat. With support from the U.S. Department of Energy, our project seeks to dramatically increase the climate-appropriate deployment of cool surfaces across the United States with an emphasis on their application to disadvantaged communities. First, we sought to identify cool-surface deployment barriers, opportunities, and models by (a) reviewing the history of cool-surface deployment activities, (b) interviewing cool-surface stakeholders, (c) researching successful energy-efficiency/green building deployment models; and (d) interviewing the actors who have implemented these models. Second, we conducted a workshop to engage stakeholders in development of a deployment plan. Third, we asked several U.S. federal agencies (a) how cool surfaces and cool surface stakeholders could support their missions and (b) how agency activities could support cool-surface deployment. Fourth, we identified a set of transformative ideas that form the core of the deployment plan. Transformative ideas include but are not limited to initiatives to (a) launch an educational campaign to make the general public and building professionals aware of how cool roofs and walls exclude unwanted solar heat; (b) create a “Cool Roof Prize” stimulating the development of affordable, high-performance cool asphalt roofing shingles; (c) conduct high-profile, large-scale demonstration programs that bring cool surfaces to disadvantaged heat-vulnerable communities; and (d) support local, regional, and state climate action (heat mitigation) plans with cool surfaces.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

3D TRISO particle-explicit compact meshing

The TRI-structural ISOtropic (TRISO) layered fuel particle is a robust nuclear fuel form offering enhanced safety and performance for advanced reactor concepts, including high-temperature gas-cooled reactors and other Generation IV designs. These poppy-seed-sized particles are embedded in a graphite matrix to form fuel elements that must withstand elevated temperatures and high burn-up levels. The heterogeneous nature of these fuel elements — comprising thousands of randomly distributed TRISO particles — produces complex stress fields and thermal gradients that one- and two-dimensional models cannot accurately capture. While three-dimensional modeling has improved predictions of dimensional changes, internal pressure buildup, and fission product transport under irradiation, current approaches rely on homogenized material properties that are known to have considerable divergence from experimental observations. This work presents a methodology for optimized random packing of TRISO fuel compacts and full three-dimensional mesh generation within the BISON fuel performance code, with each particle coating layer individually discretized. The resulting mesh was demonstrated through heat conduction simulations under representative in-reactor operating conditions, showing strong agreement with expected behavior. This capability enables detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating layer failures — all of which directly govern fuel performance and safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Initial study on cross section generation requirements for a PBR equilibrium core

A Serpent model of the HTR-PM equilibrium core was developed for use in cross section prepa-ration studies in order to guide methods development for the Griffin reactor multiphysics applica-tion. The model includes detailed isotopics for 10 distinct pebble burnup groups in 126 core zoneswith unique fuel and moderator temperatures obtained from a coupled neutronics-thermal-fluidsequilibrium core calculation using Griffin-Pronghorn. A sensitivity study of the fuel and mod-erator temperatures for various core regions was performed with the MOOSE stochastic tools.The results show that the uncertainties are, not unexpectedly, dominated by the value of the fluidtemperature and that the power level, heat transfer coefficient and effective conduction to neigh-boring pebbles and fluid constitute, at best, second order effects. The temperature uncertaintyrange varies from 28 K to 57 K at the core entry and exit planes, respectively, but these val-ues are probably higher. We still have to quantify the significance of these uncertainties in thepreparation of cross sections, which will be postponed for future work. In addition, we verifythat the single effective pebble approach works well for the preparation of region averaged crosssections in the infinite domain approximation. Nevertheless, there are significant discrepanciesin the cross sections when compared to the multi-pebble model. This could affect the predictionof peak values and in the depletion calculation. We conclude that is highly desirable for futurestudies with Griffin to be able to handle both the ?effective? pebble approximation and the multi-pebble approach for the various pebble burnup groups. This enables Griffin users the flexibilityto perform higher-fidelity studies. Finally, we initiate the preparation of cross sections for variouscore regions from the full core Serpent reference model. We quantify the differences in 26 groupcross sections from infinite domain models. These reference cross sections will serve to validatethe double heterogeneity, self-shielding, and spectrum-correction methods in Griffin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Cost optimization of high recovery single stage gap membrane distillation

Estimating the cost of membrane distillation (MD) processes at scale is confounded by the vast array of possible process configurations, module configurations, and operating conditions and the extreme sensitivity of MD costs to these parameters. This work develops a comprehensive MD cost optimization framework for minimizing the levelized cost of water (LCOW) as a function of feed salinity and temperature. Our single stage gap MD model fully captures solution properties, heat transfer, and mass transfer relationships in the modules to determine the cost optimal design as a function of multiple decision variables, including process configuration (feed location and connectivity between the hot and cold-side), the gap type (air, permeate, conductive), equipment sizes (membrane, heater, chiller, and heat exchanger units), and operating conditions (flow rates, hot and cold side temperatures, and heating and cooling duties). We estimate that the minimum LCOW of conductive gap MD ranges between $\$$10 and $\$$16 per m 3 for feed salinities of 25 g/L to 200 g/L and water recoveries of 30%–75%. Further, we update several heuristic assumptions for cost optimal system operation, including: increased thermal conductivity of the gap decreases system costs; cost optimal MD operates at the highest allowable temperature; and the LCOW is a strong function of brine salinity. Finally, we find that improving membrane properties and lowering heating and cooling costs will provide only modest reductions in costs (<25%).

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

A non-fueled nuclear-heated rod for in-pile transient boiling studies

Separate-effects boiling experiments have recently been conducted in the Transient Reactor Test Facility at Idaho National Laboratory to investigate transient heating and irradiation effects on cladding-to-coolant heat transfer. Specifically, transient critical heat flux (CHF) remains an important area of research, and better understanding of this phenomenon has potential for improving predictive models related to operational and safety limits. Consequently, this knowledge is expected to improve efficiency of light-water reactor operations. A novel borated nuclear-heated rodlet (BNHR) was designed to enable observation of transient cladding-to-coolant heat transfer phenomena. The final BNHR design takes a surrogate approach, wherein nuclear heating is induced by 10 B(n, α) reactions rather than derived from fissions in a fueled specimen. The structure of the BNHR consists of a hollowed out borated (B nat ~ 2.05 wt %) stainless steel tube with an hourglass-shaped outer surface, capped at both ends with non-borated stainless steel. This geometry allows for inner-rodlet instrumentation and generation of the highest nuclear heating rates near the center of the rodlet to ensure onset of boiling near instrumentation for real-time observation. A novel approach to measuring the nuclear energy deposition rate in the BNHR separate and apart from the influence of the coolant, termed the n-a thermometer, is also detailed in this paper. This device has demonstrated excellent repeatability, and measurements indicate predictive modeling results for energy deposition in the BNHR rod agree within a 10% margin of the experiment measurements. In conclusion, these results give confidence that the BNHR design has successfully met experiment objectives.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Variable Laser Ignition Pathways in Al/Pt Reactive Multilayers across 10 Decades of Pulse Duration

Pulsed laser irradiation is used to investigate the local initiation of rapid, self-propagating formation reactions in Al/Pt multilayers. The single pulse direct laser ignition of these 1.6 μm thick freestanding foils was characterized over 10 decades of pulse duration (10 ms to 150 fs). Finite element, reactive heat transport modeling of the near-threshold conditions has identified three distinct ignition pathways. For milli- to microsecond pulses, ignition occurs following sufficient absorption of laser energy to enable diffusion of Al and Pt between layers such that the heat released from the corresponding exothermic reaction overcomes conductive losses outside the laser-irradiated zone. When pulse duration is decreased into the nanosecond regime, heat is concentrated near the surface such that the Al locally melts, and a portion of the top-most bilayers react initially. The favorable kinetics and additional heat enable ignition. Further reducing pulse duration to hundreds of femtoseconds leads to a third ignition pathway. In this work, while much of the energy from these pulses is lost to ablation, the remaining heat beneath the crater can be sufficiently concentrated to drive a transverse self-propagating reaction, wherein the heat released from mixing at each interface occurs under kinetic conditions capable of igniting the subsequent layer.

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

Uncertainty Quantification of Monolithic Tungsten Target Block Fatigue

The current STS spallation target design consists of 21 water cooled, tantalum clad, monolithic tungsten blocks. The target wheel assembly rotates such that each target block receives a direct proton pulse every 1.4 seconds (21 blocks / 15 Hz pulses). The quasi-steady heating induced by the repetitive beam pulses creates a significant mean stress distribution in the target. On top of this, the rapid energy deposition from a proton pulse induces dynamic stress waves that are of the same order of magnitude as the mean stresses. Over the design goal of 5,000 operational hours per year for a decade, each target must withstand over 130 million repetitions of these dynamic stress events. The fatigue life prediction of the target blocks is a complex sequence of simulations that depend upon a set of input parameters. For the tungsten block, the end goal is the fatigue factor-of-safety (FOS). This study documents a new structural analysis routine using software available to U.S. national laboratories (e.g. CUBIT, Sierra Multiphysics, Dakota) that enables more rapid target fatigue analysis. The increased throughput is amenable to parameter sensitivity and uncertainty quantification (UQ) studies. A set of input parameters (ultimate strength, fatigue strength, density, elastic modulus, Poisson’s ratio, thermal conductivity, coefficient of thermal expansion, specific heat, flow rate) with expected ranges are documented based on literature data, recent experimental testing, and engineering judgement. UQ bounds the resulting fatigue FOS at [0.5, 3.4] using epistemic interval analysis. These bounds are expected to be wider due to model effects (discretization) and beam effects (focus, offset) that were not included. Parameter sensitivity underscores the importance of ultimate strength and fatigue strength of irradiated tungsten on the fatigue calculations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Advanced analysis of laser-driven pulsed magnetic diffusion based on quantum molecular dynamics simulation

Magnetic diffusion plays an important role in inertial confinement fusion with strong magnetic fields. In this paper, we improve a previous analysis of the generation and diffusion of the magnetic field [Morita et al., Phys. Plasmas 25, 094505 (2018)]. For the generation process, we calculate the temporal evolution of the coil current using a self-consistent circuit model. The results show that the peak of the calculated magnetic field is delayed by 1.2 ns compared with that of the incident laser pulse. For the diffusion process, we evaluate the electrical conductivity of warm dense gold over a wide temperature range (300 K–100 eV) by combining the Kubo–Greenwood formula based on a quantum molecular dynamics simulation with the modified Spitzer model. Our simulation shows that the maximum magnetic field (530 T) that penetrates the cone is delayed by 2.5 ns compared with the laser peak. This result is consistent with experiments [Sakata et al., Nat. Commun. 9, 3937 (2018)] that showed that applying a strong magnetic field improved the heating efficiency of fusion fuel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Wide-Bandgap Semiconductor Amplifiers for Fusion Plasma Heating and Control

This paper discusses power electronics developed under the ARPA-E GAMOW program to support nuclear fusion power production. The goal of this project was to develop and assess the potential for wide-bandgap (WBG) semiconductor devices in power electronics to enable high-efficiency and high-voltage solid-state systems for fusion plasma generation, heating, and control. The power electronics use an architecture in which multiple high-power boards can be combined to produce megawatt-level power, where using multiple boards provides high reliability. Two main areas of power electronics boards are developed in this project for fusion plasma heating and control applications: (1) pulse generation and control and (2) radiofrequency generation. The first area is for boards capable of driving high-voltage millisecond pulses at high duty cycles. The envisioned application of these pulses is in plasma control of magnetohydrodynamic instabilities, plasma position, and edge-localized modes. Pulse-width modulation allows for the implementation of a wide variety of linear and nonlinear control systems. The boards developed for this project could actuate control coils based on digital input signals and can be parallelized to provide megawatts of output power. The design of the pulse generator is a low-side load switch. A load switch was designed and constructed that utilized 2-kV-rated field-effect transistor (FET)-based cascodes developed by Qorvo under this project to perform initial testing of these cascodes. The second area is being implemented using class E amplifiers with WBG devices and a reactance steering network to handle inductive or capacitive plasma loads. Applications include ion cyclotron resonance heating (ICRH) and high-harmonic fast-wave (HHFW) heating. A class E reactance steering network is demonstrated in modeling and experiment with a resistive-inductive load that models an inductively-coupled plasma. Power combining of boards with class E reactance steering networks is also simulated and demonstrated experimentally, to enable scaling up to high power. Modeling of high-power-density cooling and remaining useful life is conducted to enable reliable, effectively cooled high-power electronics for fusion applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Influence of dissolved salts and impurities on seawater on heat transfer and fluid flow through fuel channels and debris bed

This project aims to better understand this mechanism and its effect on heat transfer. Specifically, experimental research was conducted to investigate the long term coolability issues in scaled test setups. Advanced diagnostic and measurement techniques were investigated and incorporated into an annular test facility and a debris bed test facility. The diagnostic techniques investigated included: radioactively tagged salt solution, Frisch collar radiation detector array, neutron and x-ray imaging systems, and a Rayleigh backscatter-based fiber-optic temperature system. Models were checked against the in-situ experimental data and recommendations were made for modeling purposes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Simulation of Channel Flow with Square Ribs for Blanket First-Wall Cooling: Geometry-Specific Tuning of k-ω Model Using Adjoint Method

Cooling of the plasma-facing first wall is challenging in the design of blanket components because of the high heat flux (on the order of 𝑀𝑊/𝑚2) from the plasma, especially when a low thermal mass medium like helium is chosen as the coolant. Therefore, heat transfer enhancement in which the convective heat transfer rate is augmented by the addition of turbulence-promoting structures becomes a key initiative for providing sufficient cooling capability with helium. Previously, computational fluid dynamics simulations had been performed on pipe flows with different transverse and longitudinal ribbed geometries at Oak Ridge National Laboratory to compare the enhancement performance among different ribbed geometries. Rib shape morphing had been conducted to obtain an optimized rib profile. In the work presented here, the adjoint method is adopted in the ANSYS Fluent solver for turbulence model augmentation, and the Generalized k-ω (GEKO) turbulence model is employed because of its ability of tuning the turbulence model. The Nusselt number and pressure drop obtained from the channel flow with bottom ribbed wall experiments are used as the targets. Sensitivity analysis provides information as guidance to improve the turbulence model accuracy. The augmented GEKO model is tuned for the studied ribbed channel geometry and flow conditions, providing improved predictive accuracy within this context. Extension to other configurations offers potential but may require additional tuning and validation.

Xu, Tracy [ORNL] (ORCID:0009000193700887)↗
DOE OSTI2022

Initial study on cross-section generation requirements for a PBR equilibrium core

A Serpent model of the equilibrium core HTR-PM small modular nuclear reactor in China, was developed for use in cross-section preparation studies in order to guide methods development for the Griffin reactor multiphysics application. The model includes detailed isotopics for 10 distinct pebble burnup groups in 126 core zones with unique fuel and moderator temperatures obtained from a coupled neutronics-thermal-fluids equilibrium core calculation using Griffin-Pronghorn. A sensitivity study of the fuel and moderator temperatures for various core regions was performed with the MOOSE stochastic tools. The results show that the uncertainties are, not unexpectedly, dominated by the value of the fluid temperature and that the power level, heat transfer coefficient and effective conduction to neighboring pebbles and fluid constitute, at best, second order effects. The temperature uncertainty range varies from 28 K to 57 K between the core entry and exit planes, respectively, but these values are probably higher. We still have to quantify the significance of these uncertainties in the preparation of cross-sections in future work. In addition, we verified that the effective pebble approximation used in the PEBBED and V.S.O.P. computer codes works well for the preparation of region averaged cross-sections. Nevertheless, there are some discrepancies in the cross-sections when compared to the multi-pebble model, which could affect the prediction of peak values and the depletion calculation. We conclude that is highly desirable for future studies with Griffin to be able to handle both the 'effective' pebble approximation and the multi-pebble approach for various pebble burnup groups. This enables Griffin users with the flexibility to perform higher-fidelity studies. Finally, we initiated the preparation of cross-sections for various core regions from the full core Serpent reference model. We quantified the differences in 26 group cross-sections from infinite domain models versus the full core approach. These reference cross-sections will serve to verify the double heterogeneity, self-shielding, and spectrum-correction methods in Griffin. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Demonstration of RELAP5-3D for transient analysis of a dual coolant lead lithium fusion blanket concept

Our research aims to build a multiphysics framework for transient analysis of the Dual- Coolant Lead-Lithium (DCLL) blanket design of the proposed Fusion Nuclear Science Facility (FNSF) for all Onticipated Operational Occurrences (AOOs) using Reduced Order Models (ROMs). This paper presents an assessment of using RELAP5-3D for transient thermal-hydraulic blanket analysis of the DCLL design. Preliminary validation studies for the RELAP5-3D properties were conducted using a representative vertical flow loop model and a comparison study using experimental data from the thermal convection corrosion loops at ORNL. These studies inspected the thermal hydraulic response of Lead-Lithium Eutectic (PbLi) systems within the RELAP5-3D code; mainly pressure drop and heat transfer, as compared to analytical calculations and experimental data. The Magnetohydrodynamic (MHD) pressure drop effect for liquid metal under the effects of a uniform magnetic field was implemented into our models using an equivalent forms loss coefficient based on the phasic momentum equations. We verified our implementation of MHD pressure drop within RELAP5-3D using a systematic approach consisting of two studies. Models were developed to perform comparison studies based on experimental data from the Argonne Liquid metal EXperiment (ALEX) facility and the Magnetohydrodynamic PbLi Experiment (MaPLE) facility. Using the developed validation basis, we built a simplified DCLL channel model within RELAP5-3D, inspected MHD pressure drop and performed transient temperature analysis. The DCLL channel model used radial heating data from the International Thermonuclear Experimental Reactor (ITER) which was normalized to the FNSF design. Here, a representative startup transient was performed using a relative power curve from the initial ramp up of a long-term power pulse of the Tore Supra reactor. This study demonstrated that RELAP5-3D ROMs are capable of accurate transient analysis of the DCLL blanket. We have since utilized MCNP to develop heating profiles based on the FNSF design and are working toward implementing them into our DCLL model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗