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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Cylindricity Sensitivity Thermal Model of the AGR-5/6/7 Experiment in the Advanced Test Reactor

The AGR-5/6/7 experiment is currently being irradiated in the Advanced Test Reactor (ATR) at the Idaho National Laboratory and is approximately 70% complete. Several fuel and material irradiation experiments have been planned for the U.S. Department of Energy Advanced Gas Reactor Fuel Development and Qualification Program, which supports the development and qualification of tristructural isotropic (TRISO) coated particle fuel for use in high-temperature gas-cooled reactors. The goals of these experiments are to provide irradiation performance data to support fuel process development, qualify fuel for normal operating conditions, support development of fuel performance models and codes, and provide irradiated fuel and materials for post-irradiation examination and safety testing. Originally planned and named as separate fuel experiments, but subsequently combined into a single test train, AGR-5/6/7 is testing low-enriched uranium oxycarbide TRISO fuel. The AGR-5/6/7 test train has 5 capsules with thermocouples and independent gas control mixtures. Unique to this paper is a sensitivity study concerning the cylindricity of the graphite holders containing the fuel compacts and their eccentricity in relation to the stainless-steel capsule walls. Each capsule has small nubs on the outside used for centering the graphite holder inside the stainless-steel capsule with a small gas gap used to control temperature. Due to machining tolerances of these nubs and vibration wearing the nubs down when the experiment is running in the reactor, the possibility exists that the holder may move around radially. Each capsule is equipped with several thermocouples placed at various radii and depths within each graphite holder. This paper will show the sensitivity of offsetting the graphite holder for various radii in 45-degree increments around the circle with the objective of minimizing the difference between the measured thermocouples and the modeled thermocouple temperatures. Separate gas mixtures of helium/neon are introduced into this gas gap between the holder and capsule wall and changed as necessary to maintain the desired thermocouple temperatures to keep the fuel compacts at constant temperature as the nuclear reactor conditions change. The goal of the sensitivity study is to find a radius and an angle to offset the holder from perfectly centered for each of the five capsules separately. The complex thermal model includes fission heating, gamma heating, radiation heat transfer, and heat transfer via conduction and radiation across the control gaps. Subroutines linked to the thermal model offer an easy method to offset the graphite holder from the capsule walls without remeshing the entire model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Demonstration of Pronghorn’s Subchannel Code Modeling of Liquid-Metal Reactors and Validation in Normal Operation Conditions and Blockage Scenarios

Pronghorn-SC is a subchannel code within the Multiphysics Object-Oriented Simulation Environment (MOOSE). Initially designed to simulate flows in water-cooled, square lattice, subchannel assemblies, Pronghorn-SC has been expanded to simulate liquid-metal-cooled flows in triangular lattices, hexagonal subchannel assemblies. For this purpose, the algorithm of Pronghorn-SC was adapted to solve the subchannel equations as they are applicable to a hexagonal wire-wrapped sodium-cooled fast reactor. Cheng–Todreas models for pressure drop and cross-flow models were adopted and a coolant heat conduction term was added. To solve these equations, an improved implicit algorithm was developed robust enough to deal with the numerical issues, associated with low flow and recirculation phenomena. To confirm the prediction capability of Pronghorn-SC, calculations and comparisons with available experimental data of 19- and 37-pin assemblies were performed, as well as other subchannel codes. Finally, a flow blockage modeling feature was added. This capability was validated for both water-cooled square sub-assemblies and sodium-cooled hexagonal sub-assemblies, using experimental data of partially and fully blocked cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

High-Fidelity Multiphysics Modeling of a Heat Pipe Microreactor Using BlueCrab

Researchers who are actively developing nuclear microreactors are planning to employ innovative designs and features using traditional commercial modeling tools that may be inadequate for their design and licensing activities. The codes developed under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program provide flexibility in terms of geometry modeling and multiphysics coupling and are particularly well suited for modeling novel microreactor concepts. To test the maturity of these codes, this paper introduces a conceptual heat pipe microreactor (HP-MR) designed to gather various technologies of interest to microreactor developers such as control drums, heat pipes, and hydride moderators. Here, the objective of this effort is to demonstrate NEAMS tools capability to perform high-fidelity multiphysics simulations, using coupled neutronics (via the Griffin code), heat conduction (via the BISON code), heat pipe modeling (via the Sockeye code), and hydrogen redistribution in hydride metal moderator (via the SWIFT code). Codes are coupled in-memory through the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, which permits flexible multiphysics data transfer schemes. The analysis confirmed two key aspects of the HP-MR concept: (1) its ability to follow the power load requested from the heat pipe and (2) its ability to avoid heat pipe cascading failure unless designed with high power close to operating failure limits of its heat pipes. The developed computational model was distributed publicly on the Virtual Test Bed for training purposes to accelerate adoption by industry and to provide a high-fidelity multiphysics solution for benchmarking against other tools. Additional multiphysics analyses including other transients and coupled physics were identified as necessary future work, together with a focus on validating multiphysics behavior against experiments.

Microreactor↗
DOE OSTI2025

Moltensaltpropnet

MoltenSaltPropnet is a physics-informed machine learning framework that aims to predict the thermophysical properties of molten fluoride and chloride salt mixtures, which are crucial for the design and safety of Generation IV molten salt reactors. The code processes data from the Molten-Salt Thermal Properties Database (MSTDB-TP) and the Janz compendium, converting critically evaluated correlations into fast, differentiable surrogate models for density, viscosity, thermal conductivity, and heat capacity across 448 distinct salt systems. The implementation consists of several key components: 1. Data Curation: The code parses and cleans the raw data, normalizing elemental mole fractions and extracting relevant regression coefficients for various thermophysical properties. 2. Feature Engineering: It generates fixed-length numerical descriptors that encapsulate the composition and temperature, incorporating polynomial interaction terms and dimensionality-reduction techniques to optimize model performance. 3. Coefficient Learning: Four different machine learning architectures are employed: a deep residual network (ResNet), a Kolmogorov–Arnold network (KAN), a sparsity-inducing neural network (SNN), and classical regression models. Each model learns to predict coefficients that define the temperature-dependent correlations for the thermophysical properties. 4. Property Reconstruction: The predicted coefficients are used to compute temperature-dependent property values, ensuring positivity and monotonic trends through a composite loss function that enforces physical constraints. 5. User Interface: An open-source web application enables users to filter the database, train task-specific models, and visualize the results, allowing for rapid exploration of candidate salt mixtures. MoltenSaltPropnet bridges the gap between limited experimental data and high-fidelity reactor simulations, providing a powerful tool for researchers in the field of molten salt reactors and advanced nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗
DOE OSTI2023

Effect of Lateral Conduction in a Rectangular Fuel Plate and Optimum Stripe Widths for the OPAL Reactor

The impetus for this work is the Proliferation Resistance Optimization (PRO-X) program, which was created by the National Nuclear Security Administration (NNSA) of the US Department of Energy (DOE) to provide a framework for developing reactor designs that minimize the production of special nuclear materials while maximizing performance for peaceful uses. The Argentine company INVAP and Argonne National Laboratory collaborated to share knowledge and expertise regarding the thermal analysis of research reactors. Working together, they designed a series of test problems based on the Open Pool Australian Light Water (OPAL) reactor, whose fuel is in the form of flat plates cooled by water flowing through narrow rectangular channels. ANL used their PLTEMP/ANL code. INVAP used their TERMIC code (and its new multi-plate version, TERMIC-MP). Then, INVAP and ANL compared the two sets of results. 1-D models were used during the collaboration since for research reactor analysis, largely 1-D thermal-hydraulics models are commonly used because they are simple to apply and technically defendable. The collaboration between the two organizations aims for improvements in thermal-hydraulic models. Improvements to methods and models can reduce modeling uncertainties and excessive conservatism, allowing greater reactor performance without reducing predicted safety margins. Therefore, these improvements potentially expand the design space, making additional designs feasible, which, in turn, may achieve better fuel utilization and proliferation resistance in the designs. This report provides justification for the lateral node size (stripe width) chosen for the 1-D computer models during the collaboration. The lateral heat conduction in OPAL fuel plates is evaluated using the two-dimensional conduction capability of PLTEMP/ANL to determine a representative hot-stripe width that can be conservatively used for one-dimensional steady-state thermal-hydraulic analysis. Basically, a representative hot-stripe width needs to be selected so that the actual heat flux averaged over this stripe width is greater than or equal to the maximum heat flux obtained by the two-dimensional analysis. Detailed results in this report show that a hot-stripe width of 8.125 mm can be conservatively used for the one-dimensional thermal-hydraulics analysis of all plates. Therefore, eight equal stripes (8.125 mm each) are recommended to be modeled in the MCNP and PLTEMP/ANL full core models of the OPAL reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Potassium carbonate decomposition modeling within rotating detonation engines for direct power extraction applications

We report that recently, there has been a significant interest in detonation-based combustion systems, such as rotating detonation engines (RDE), due to potential performance advantages in propulsion and energy applications. This includes improvements in thermodynamic efficiency, inherently high gas velocities, and the possibility of achieving an increase in total pressure (i.e. pressure gain) through the combustion process. Due to approximating a constant volume combustion process, these devices also tend to exhibit extremely high local gas temperatures relative to a comparable constant pressure combustor. These advantages overlap well with desirable performance characteristics of direct power extraction (DPE) technologies, such as a magnetohydrodynamic (MHD) generator. Typically, in DPE systems hot combustion products are seeded with an easily ionizable material such as potassium carbonate (K 2 CO 3 ) in order to boost the electrical conductivity. However, due to the short gas residence times within an RDE, it was unclear whether forming an electrically conductive combustion plasma would be feasible for integration with a downstream MHD generator. A model is presented which describes the heating, decomposition, and ionization of solid K 2 CO 3 particles and aqueous solutions of K 2 CO 3 in water, for a given initial particle diameter. This model was combined with available computational fluid dynamics (CFD) data for an oxygen-methane RDE in a one-way coupled Eulerian-Lagrangian framework to predict particle trajectories and the corresponding heating, decomposition, and ionization histories. Electrical conductivities were computed using a previously developed model, and a method was proposed to determine an equivalent average electrical conductivity. Results show that particle sizes below ~30 µm are able to fully decompose before reaching the exit of the RDE. While the one-way coupled nature of the simulations precluded rigorous evaluation of the effects of seed material on detonability, a substantial temperature reduction is expected at the detonation wave front due to heating and decomposition. A preliminary comparison is presented between the RDE and an equivalent constant pressure adiabatic combustor, showing a potential performance advantage for the RDE.

42 ENGINEERING↗
DOE OSTI2020

Composite Structured Surfaces for Durable Dropwise Condensation

Dropwise condensation of steam on hydrophobic substrates has a 10X higher heat transfer coefficient compared to filmwise condensation. To promote dropwise condensation, low surface energy hydrophobic coatings (polymers) are typically utilized. The low intrinsic thermal conductivity (k < 1 W/(m · K)) of polymers, coupled with high heat transfer coefficient of dropwise condensation (100 kW/(m 2 · K)), necessitates that the coating be thin (< 1µm) in order to avoid reducing the overall heat exchanger conductance. However, thin polymeric films easily degrade. The two opposing requirements result in the need for optimization between the durability (thick coating) and the heat transfer (thin coating). To enable high thermal conductivity in thicker coatings, we develop metal-polymer structured surfaces. By using porous structures as inter-connected heat-conducting backbones that are filled with hydrophobic materials, we enable tuning of the coating effective thermal conductivity and surface energy. Three metal structures were studied; micro/nanowires, inverse opals, and sintered spheres. Heat transfer performance was calculated using three-dimensional finite element method simulations with two distinct boundary conditions; convection at the walls and isothermal walls. Interestingly, the overall conductance shows up to 40% difference depending on the boundary condition used in calculating the composite coating effective thermal conductivity. We use our model to predict the heat transfer performance as a function of metal fraction by volume and by surface area for condensation. By coupling our thermal simulations with a previously verified analytical model for predicting wetting behavior on heterogeneous surfaces, we propose a regime map to predict dropwise-to-filmwise transition. Furthermore, our work not only forms a starting point for the development of durable dropwise condensing surfaces, it identifies important considerations needed for computing effective thermal conductivity of composites.

42 ENGINEERING↗
DOE OSTI2022

Impact of electron transport models on capillary discharge plasmas

Magnetohydrodynamics (MHD) can be used to model capillary discharge waveguides in laser-wakefield accelerators. However, the predictive capability of MHD can suffer due to poor microscopic closure models. Here, we study the impact of electron heating and thermal conduction on the capillary waveguide performance as part of an effort to understand and quantify uncertainties in modeling and designing next-generation plasma accelerators. To do so, we perform two-dimensional high-resolution MHD simulations using an argon-filled capillary discharge waveguide with three different electron transport coefficients models. The models tested include (i) Davies et al., (ii) Spitzer, and (iii) Epperlein–Haines (EH). We found that the EH model overestimates the electron temperature inside the channel by over 20% while predicting a lower azimuthal magnetic field. Furthermore, the Spitzer model, often used in MHD simulations for plasma-based accelerators, predicts a significantly higher electron temperature than the other models suggest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

DECOVALEX-2023: Task E Final Report

This is the Task E final report for DECOVALEX-2023. Task E is focused on understanding thermal, two-phase hydrological, and mechanical (TH2M) processes, especially related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and essentially non-porous in the far field (away from excavations). Investigation of the short-term (days to years) near-field (centimeters to tens of meters) behavior of salt is important for radioactive waste disposal because this early period strongly controls the amount of brine in a salt repository. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment. The main test case used in Task E is the ongoing Brine Availability Test in Salt (BATS) heater test located underground at the Waste Isolation Pilot Plant (WIPP) near Carlsbad, New Mexico, USA. The Task was divided into a series of Steps. Step 0 was an introduction to processes in salt, that included matching historical unheated brine inflow data from boreholes at WIPP and matching temperature observations during BATS heater test 1a. Step 1 included validation of models against a thermo-poroelastic analytical solution relevant to heated boreholes in salt, and two-phase flow around an excavation in salt. Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the BATS 1a heater test. There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. The key learning points from Task E have been: • Heat conduction through salt typically requires non-linear thermal conductivity (as a function of temperature), but most models do a good job matching observations, given appropriate adjustments to the applied power and some thermocouple locations. • Thermal pressurization requires coupled thermal-hydrological-mechanical (THM) responses that consider the thermal expansion of the fluid and solid phases. • Initialization of two-phase flow models around a borehole or excavation in salt are more realistically represented as “wetting up”, rather than “drying down” (i.e., the initial state after excavation is mostly dry, rather than mostly wet). • The BATS 1a heater test includes a significant release of brine after the end of heating, which requires a large increase in permeability to recreate. Task E has been a great learning experience for all the teams involved, and feedback from the modeling teams has led to changes in the design of follow-on BATS experiments, which are now ongoing underground at WIPP. There was a balance throughout the task between freedom to model phenomena how each team saw fit, and prescriptiveness in problem design to bring the modeling teams closer together to allow attribution of smaller differences between models to different modeling choices. The modeling approaches seem to go through two phases: an early phase of discovery or testing, and a later phase of refinement and improvement. In future modeling efforts, different field data could be used (e.g., BATS 2) and more time should be included in the processes for teams to make multiple model refinement or even significant changes to their conceptual model or setup, based on lessons learned from the modeling exercise.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Thermo-Hydro-Mechanical Modeling of Brine Migration in a Heated Borehole Test in Bedded Salt

Abstract This research paper focuses on the thermo-hydro-mechanical (THM) modeling of brine migration in a heated borehole test conducted as part of the ongoing Brine Availability Test in Salt (BATS) at the Waste Isolation Pilot Plant (WIPP) in New Mexico. It is a component of the international collaboration project DECOVALEX-2023 (DEvelopment of COupled models and their VALidation against EXperiments), which aims to understand the THM processes governing brine flow in heated rock salt repositories through collaborative analysis by multiple research teams. Using the TOUGH–FLAC simulator, THM simulations were performed and compared with data from the BATS phase 1a. This experiment involved two identical horizontal-borehole arrays, one heated and one serving as a control, both equipped with sensor arrays. Analysis of measurements revealed water flow rate surges during heater power transitions, with the highest jump observed during cooling. Acoustic emission activity exhibited distinct patterns in response to heater power changes, suggesting that damage to rock salt is particularly pronounced during the cooling phase. The THM simulations successfully captured these phenomena, highlighting the significance of thermal effects, brine migration, and mechanical behavior in predicting brine availability in heated and damaged rock salt. Our modeling also revealed the critical interplay between heating and cooling-induced damage and its influence on flow properties, particularly affecting brine inflow estimation. Notably, we found that cooling-induced brine inflow spikes result from increased permeability due to tensile dilatancy. These findings have important implications for the development of robust containment strategies and enhance our understanding of the complex processes involved in repository performance.

58 GEOSCIENCES↗
DOE OSTI2022

Coupling subsurface and above-surface models for optimizing the design of borefields and district heating and cooling systems in the presence of varying water-table depth. In: Proceedings, 46th Workshop on Geothermal Reservoir Engineering

Dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes, without flow of groundwater. On the other hand, detailed simulation tools for subsurface heat and mass transfer exist, but these fall short in simulating above-surface energy systems. To support the design and operation of such systems, we have developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the ground. For the first, we use the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, we use the TOUGH simulator. TOUGH can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. In previous work, we described the coupling of these software packages, including how time-dependent boundary conditions for the borehole walls are synchronized for use in Modelica and TOUGH. We verified that the coupled Modelica/TOUGH code produced consistent results with the original Modelica code for an idealized problem in which heat transfer was purely by conduction in a uniform geologic medium. Here, we examine less idealized problems for which TOUGH’s advanced capabilities for modeling fluid flow are required. The first problem has a shallow vadose zone, and the second problem has a thicker vadose zone with a water-table depth that varies in time, which requires a fine vertical grid discretization for the TOUGH model.

Doughtry, Christine↗
DOE OSTI2020

FAST-1.0 Software Release Document

The purpose of this document is to detail the release of FAST-1.0 and the documents that accompany the development and testing of FAST-1.0. Design requirements were for FAST-1.0 developed in the NRC Statement of Work (SOW) for NRC Agreement Number NRC-HQ-25-14-D-001 and are transcribed to this document as part of the Software Quality Assurance Plan for the FAST Code System, PNNL-28767. FAST 1.0 represents the merger of the FRAPCON and FRAPTRAN codes and as such FAST-1.0 performs steady state and transient fuel performance calculations described below. FAST-1.0 calculates the steady state response of light-water reactor fuel rods during long-term burnup. The code calculates temperature, pressure, and deformation of a fuel rod as functions of time-dependent fuel rod power and coolant boundary conditions. The phenomena modeled by the code include: 1) heat conduction through fuel and cladding to the coolant; 2) cladding elastic and plastic deformation; 3) fuel-cladding mechanical interaction; 4) fission gas release from the fuel and rod internal pressure; and 5) cladding oxidation. FAST is used to perform independent calculations for regulatory evaluations of fuel performance under normal operation and anticipated operational occurrences (AOOs).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

CHARACTERIZATION ON ANISOTROPIC THERMAL CONDUCTIVITY FOR BIG AREA ADDITIVE MANUFACTURING WITH POLYMERS

Additive manufacturing with polymers has been used mainly for prototyping. A recent development of Big Area Additive Manufacturing (BAAM) at Oak Ridge National Laboratory has opened its applications in the mold and die industry. A numerical simulation and prediction for a mold heating performance requires accurate anisotropic thermal properties of the printed material, which are challenging to obtain, and often requires the use of multiple techniques. The transient plane source (TPS) technique has been widely used due to its ability to measure the thermal properties of an extensive range of materials (solids, liquids, powder). Despite the capability to characterize thermal conductivity k of isotropic and anisotropic materials, the measurements of latter materials are limited to the cases, where the samples have the same thermal conductivity k along x- and y-axis that form the radial plane. In this work, the method for a characterization of k in all three dimensions is developed, and the application of TPS is extended to the determination of thermal properties along the x-, y-, and z-axis individually. The materials are represented by additively manufactured polymers including polylactic acid (PLA) and styrene maleic anhydride (SMA). The developed method consists of (1) a determination of the heat capacity of the polymers by means of TPS in combination with the developed in this work data analysis procedure, (2) a machining three types of cylindrical samples from the same material, with the height corresponding either to x-, y-, or z-direction of printing, and (3) a determination of axial thermal conductivity employing anisotropic model and using previously determined heat capacity

Trofimov, Artem↗
DOE OSTI2020

CFD-DEM and PR-DNS studies of low-temperature densely packed beds

Over the past few decades, granular media is gaining attention as a viable option for heat transfer fluids (HTFs). Several research efforts are studying the use of particle-based heat transfer fluids in a wide variety of applications. With this motivation, the current work focusses on analyzing the different heat transfer mechanisms in low-temperature mono-sized densely packed granular media. To study the heat transfer behavior of granular media at different scales, the current work employs a two-way coupled computational strategy. The motion of particles is solved using the Discrete Element Method (DEM) and the interstitial air is solved using a Finite-Volume (CFD) approach. The Open-Source library CFDEM Coupling® is used in the current study to join the Finite Volume PISO solver of OpenFOAM® and the DEM solver of LIGGGHTS®. Typically, particle-particle contact conduction and particle-air convection are the most popular closure models. But recent research identified a different heat transfer phenomenon in packed beds that cannot be identified by conduction or convection models. While closure models were developed to implement this on a CFD-DEM framework, they did not capture the effect of intra-particulate thermal gradients on this phenomenon. Hence the current work also employs Particle-Resolved Direct Numerical Simulations (PR-DNS) to gain valuable insights allowing for the modification of existing models. A new closure model is then proposed here and is implemented in the CFD-DEM framework. This model provides key insights into the different heat transfer mechanism of packed beds.

42 ENGINEERING↗
DOE OSTI2020

A dynamic thermal algal growth model for pilot-scale open-channel raceways

Here this study describes the makeup of and results from a simulation consisting of a lumped thermal model integrated with a dynamic algae growth model to simulate the microalgae productivity of an open-channel raceway cultivation system. The thermal model considers the dynamic effects of weather, light absorption, convective heat transfer, radiation heat transfer, conductive heat transfer, thermal capacitance, and water control strategies. The dynamic algae growth model solves a set of ordinary differential equations consisting of growth functions dependent on incident radiation, temperature, nutrient availability, basal metabolism, and losses due to dark- and photo-respiration. Relative errors in the predicted ash-free dry weight of Nannochloropsis oceanica are 12.5%, -6.1%, and 4.4% for three separate replicates of cultivation cycles (~4.2 m 2 ) performed at AzCATI during the ATP 3 Unified Field Studies in Fall, Spring, and Summer. This research demonstrates that thermal modeling is an essential contributor to the validation of microalgae growth models.

09 BIOMASS FUELS↗
DOE OSTI2020

Completion Effects on Diagnosing Multistage Facture Treatments with Distributed Temperature Sensing

Summary Distributed temperature sensing (DTS) is a valuable tool to diagnose multistage hydraulic fracture treatments. When a stage interval is shut in, the clusters that take more fluid during pumping warm up more slowly. Therefore, the fluid volume injected into each cluster can be quantitatively interpreted by numerical inversion of the warm-back temperature behavior. This general concept assumes that the different warm-back behavior is controlled by only the injected fluid volume; however, recent observations of DTS data indicate that completion configurations significantly influence the warm-back behavior. This paper investigates the completion effects on the DTS interpretation. In ideal conditions, when a stage is fractured, the upstream stage intervals should show an almost uniform temperature that is close to the injected fluid temperature. This is due to the high fluid velocity of injected fluid in the wellbore, and the upstream intervals have not been perforated (noncommunicating intervals), so the only heat transfer is heat conduction between the wellbore fluid and the surrounding reservoir. But the field DTS data show considerably irregular variations in temperature along the upstream stage intervals. These variations are caused by the completion effects. The nonuniform temperature profile is caused by different heat transfer behavior induced by completion hardware along the production casing string, such as joints, clamps, and blast protectors, and by the sensing cable location in the cement, as well as the cement quality. Because the varying heat transfer behavior impacts the warm-back behavior as well as the temperature profile, the completion effects need to be considered in DTS interpretation. A method of DTS interpretation considering the completion effects to diagnose multistage fracture treatments was developed. Because the heat transfer between a wellbore and a reservoir depends on the overall heat transfer coefficient describing heat conduction through the completion in a forward model, this parameter needs to be tuned all along the wellbore. To calibrate the completion effect, the temperature inversion is conducted using the temperature measured at a stage interval that is upstream of a stage interval currently being treated. Because the interpreted stage interval is not perforated at that time, the thermal behavior at the noncommunicating interval is governed by only the heat conduction through the completion environment. Once the effective values of the overall heat transfer coefficient are estimated along the interpreted stage interval, they can be assumed to be constant physical parameters. Then, the fluid volume distribution is interpreted by using the effective overall heat transfer coefficient profile along each interval. This study provides a field application of the developed interpretation method. The new interpretation method provides more accurate diagnosis of fracture treatments by DTS interpretation.

Engineering↗
DOE OSTI2026

Statistical inference of anomalous thermal transport with uncertainty quantification for interpretive 2D SOL models

The critical task of inferring anomalous cross-field transport coefficients is addressed in simulations of boundary plasmas with fluid models. A workflow for parameter inference in the UEDGE fluid code is developed using Bayesian optimization with parallelized sampling and integrated uncertainty quantification. In this workflow, transport coefficients are inferred by maximizing their posterior probability distribution, which is generally multidimensional and non-Gaussian. Uncertainty quantification is integrated throughout the optimization within the Bayesian framework that combines diagnostic uncertainties and model limitations. As a concrete example, we infer the anomalous electron thermal diffusivity $\chi_\perp$ from an interpretive 2D model describing electron heat transport in the conduction-limited region with radiative power loss. The workflow is first benchmarked against synthetic data and then tested on H-, L-, and I-mode discharges to match their midplane temperature and divertor heat flux profiles. We demonstrate that the workflow efficiently infers diffusivity and its associated uncertainty, generating 2D profiles that match 1D measurements. Future efforts will focus on incorporating more complicated fluid models and analyzing transport coefficients inferred from a large database of experimental results.

Bayesian optimization↗
DOE OSTI2022

Representative-cell-based Modeling of HTTF in SAM

Although the current nuclear power market is primarily occupied by the light water reactors (LWRs), the TRISO-fueled helium-cooled graphite-moderated high temperature gas reactors (HTGRs) are drawing growing attention as the nuclear power industry marches towards more advanced systems. Compared to that of conventional LWRs, the thermal hydraulics of HTGR cores show extra complexity from their multi-scale heat transfer mechanisms with varying importance during different operation modes or transient stages. The modeling of HTGR cores for the system analysis purpose therefore faces the challenge of reaching reasonable fidelity and accuracy while maintaining sufficient simplicity. Previous development and validation efforts have demonstrated that a “2-D ring model” with reasonable performance can be implemented in the System Analysis Module (SAM) for prismatic HTGR cores. In the current project, an alternative modeling methodology based on “representative cells” is proposed for the typical prismatic core of HTGRs. Different from the previous ring model approach, the proposed methodology first separately models and then combines the small- and large-scale thermal conductions, by connecting representative local heat transfer units (cells) with effective core- wise thermal resistance. A modeling practice for an integral high-temperature test facility (HTTF) elaborates the modeling details. Steady-state validation of the resultant candidate model is performed against a higher-resolution benchmark from a 3D-1D coupled simulation, which shows satisfactory prediction with reasonably captured global parameters and well-represented temperature fields. A postulated pressurized conduction cooldown (PCC) is also simulated and analyzed, demonstrating the model’s capability of transient prediction with physically captured phenomena resolved in both small and large temporal and spatial scales. In general, the current work achieves a preliminary success in proposing a methodology using representative cells to model prismatic HTGR cores in SAM. Future efforts are envisioned with more validation activities and with potential modeling improvements to eventually achieve the high confidence on a high-fidelity robust modeling methodology with reasonable accuracy.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗