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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Analytical homogenization techniques applied to the Fickian diffusion: Effective diffusivity coefficient

For multiple applications in nuclear energy, the ability to accurately represent material behavior with a simplified model is important to facilitate practical engineering-scale simulations. In this work, we focus on the homogenized thermal response of a medium containing spherical inclusions, similar to a fuel form (compact or pebble) containing TRISO particles. An extensive survey on effective thermal conductivity modeling was performed in our previous study, considering a random distribution of mono-sized spherical inclusions in a continuous matrix. Using the analogy between heat conduction and the simplified Fickian diffusion (or fission product species conservation), we can use the same analytical homogenization methods to obtain ETC as for the effective diffusivity coefficient (EDC). We performed several numerical experiments at varying conditions to assess the validity of our hypothesis for EDC calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Techno-Economic Analysis for a Potential Geothermal District Heating System in Tuttle, Oklahoma: Preprint

Geothermal deep direct use (DDU) has potential across a wide swath of the United States but is underutilized due to challenging project economics associated with developing a deep geothermal resource for a large-scale and variable heat demand. The National Renewable Energy Laboratory (NREL) and University of Oklahoma evaluated the feasibility of a geothermal district heating (GDH) and cooling system in two schools and 250 houses by utilizing existing oil and gas (O&G) wells in Tuttle, Oklahoma. Heating and cooling demand in the two schools and a typical single-family house were modeled using EnergyPlus building energy simulation software. The modeling results indicated that annual heating demand in two schools and 250 houses is approximately 2.61 GWhth, and cooling demand in the two schools is approximately 2.65 GWhth. In this scope, the techno-economic analysis (TEA) was conducted using the GEOPHIRES tool combined with the TOUGH2 reservoir simulator. The reservoir performance, including geothermal heat production capacity, was modeled by the reservoir simulator TOUGH2. Then, levelized cost of heat (LCOH) was calculated using GEOPHIRES version 3.0, which includes new features such as hourly heat load optimization and peak performance evaluation. Geothermal reservoir temperature was estimated as 90.5 degrees C at a total depth of 3.3 km by the regional average temperature gradient of 22.8 degrees C/km and validated by cation geothermometer calculations. Four production scenarios with two different well configurations and two different heat load profiles have been developed for well flow rates ranging between 3.1 kg/s and 9.3 kg/s. The LCOH of the district heating and cooling system was calculated between $95 and $210/MWh ($28/MMBtu to $62/MMBtu) for four different production scenarios. Typical natural gas prices for residential customers in Oklahoma have ranged from 9 to 19 $/MMBtu over the past decade, which indicates a challenge for deployment of such a GDH system.

deep direct-use↗

DOE OSTI2023

Techno-Economic Analysis for a Potential Geothermal District Heating System in Tuttle, Oklahoma

deep direct-use↗

DOE OSTI2024

Evaluating Heat Extraction Performance of Closed-Loop Geothermal Systems with Thermally Conductive Enhancements in Conduction-Only Reservoirs: Preprint

We investigated the impact of disc-shaped and linear thermal enhancements on performance of closed-loop geothermal systems using steady-state and transient COMSOL numerical simulations. Thermal enhancements refer to thermally conductive material introduced in the rocks surrounding a wellbore to compensate for the relatively low rock thermal conductivity and increase heat extraction. Materials proposed include composite cements with thermally conductive additives such as metals or graphite to obtain overall cement thermal conductivities of 10 to 100x the rock thermal conductivity. Our simulation approach was modeling a subsegment (e.g., 100 m long) of a wellbore with idealized (e.g., perfect disc-shaped) thermal enhancements, and we report results as relative increase in heat extraction with respect to systems without thermal enhancements. We did not evaluate technical feasibility or cost of installing thermal enhancements. Simulations indicate that for 5-mm thick disc-shaped thermal enhancements with 5-m radius, repeated every 1 m along the wellbore and with thermal conductivity of 100x the rock thermal conductivity, the thermal output increases roughly 20% with respect to a closed-loop system without thermal enhancements. For fishbone structure thermal enhancements, we estimate a roughly 7% increase in thermal output for 5-m long radially outward pipes with half the radius of the main wellbore, repeated every 2.5 m, with thermal conductivity of 100x the rock thermal conductivity.

advanced geothermal systems↗

DOE OSTI2025

FY-25 Progress on Computational Modeling of the Water Based NSTF

This report summarizes the system level modeling using RELAP5-3D of the Natural Convection Shutdown Heat Removal Test Facility (NSTF) completed in FY25. This year’s work focuses on a new tank configuration where the inlet of the tank was lowered in elevation by 45”. The stability boundaries of the NSTF are thoroughly studied and stability maps are constructed based on the stability and the oscillation patterns of the system. Five distinct operational modes are identified, namely single-phase liquid, uniform double peak oscillations, uniform sinusoidal oscillations, stable two-phase flow, and non-uniform oscillations. Next, the riser inlet throttling case of experimental test Run-104 is simulated with the RELAP5 model where good agreement is obtained between the model and the experimental data. The simulation also highlights the effects of backflow of water from the tank to the upper region of the chimney. Additionally, the decay heat removal test of Run-99 is simulated with the RELAP5 model. Comparison is carried out between this run and a similar run with the mid-tank inlet of Run-74 performed in FY22. With the lower tank inlet, the RELAP5 model is able to predict the experimental data more accurately than the previous mid tank inlet configuration. The discrepancy in model prediction accuracy highlights the non-symmetrical spatial effects in the tank that would otherwise be more easily captured with higher fidelity models. Lastly, two exploratory studies are conducted to investigate the behaviors of the NSTF when 1) heating is provided to the downcomer and 2) a bypass channel is added between the horizontal chimney section to the downcomer.

42 ENGINEERING↗

DOE OSTI2026

Thermal-Hydraulic-Mechanical Modeling and Simulation of Sodium-Potassium–Cooled MARVEL Microreactor Core

The U.S. Department of Energy's Microreactor Program, with Idaho National Laboratory's development of a nuclear microreactor applications test bed named MARVEL, aims to support R&D for the deployment of small, transportable reactors across civilian, industrial, and defense sectors. The MARVEL microreactor, an 85-kWth thermal fission reactor, incorporates TRIGA nuclear fuel and a sodium-potassium eutectic as its primary coolant, designed for safety and efficiency, with natural circulation eliminating the risk of critical heat flux conditions. The reliance on natural circulation for primary cooling means the reactor avoids using fuel spacers to minimize core pressure drop, which could disrupt the primary coolant's natural flow. However, the reactor core’s tight P/D ratio of 1.05, in the absence of fuel spacers, could pose a risk of fuel rod contact and increased peak cladding temperatures. To ensure the reactor safety, this study conducted computational modeling and simulations to investigate the reactor's thermal-hydraulic-mechanical characteristics, including the reactor core heat transfer coefficients, the potential for rod-to-rod contact, and assessed its impact on peak cladding temperature and overall reactor safety. The computational analyses of the MARVEL microreactor core revealed that the thermal deformation of fuel rods under worst-case scenario may lead to the fuel rod contact, but the peak cladding temperatures will remain significantly lower than the safety criteria, ensuring the safety operation reactor without fuel spacers under normal operating conditions.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

DOE OSTI2024

Modeling the Impacts of Material Properties on Oscillatory Neuron Behavior

In this study, neuromorphic computing, which mimics the functions of biological brains, offers improvements in both latency and energy efficiency over typical von Neumann computing architectures. Spiking neural networks can be especially power-efficient because they encode information temporally and can use more sparse electrical inputs. Here, we study the design of volatile memristors (variable resistors with memory) for neuronal devices, with particular consideration toward the feasibility of all-on-chip oscillation using built-in capacitance. We use circuit simulations to model the behavior of oscillator neurons with a range of realistic material properties. We find that energy inputs increase with insulating-phase resistivity, thermal conductivity, and device aspect ratio. However, we also find that the minimum capacitance needed for oscillation decreases with increasing insulating-phase resistivity, which opposes the constraints for power efficiency. Based on published data on NbO 2 , VO 2 , and EuNiO 3 , we find that existing materials can be engineered for all-on-chip spiking using their parasitic capacitance.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Hydro-mechanical behavior of heated bentonite buffer for geologic disposal of high-level radioactive waste: A bench-scale X-ray computed tomography investigation

The behavior of heated bentonite buffer is critical for the security and long-term performance of a geological repository for high-level radioactive waste (HLW). While laboratory column experiments have been conducted to investigate compacted bentonite and coupled THMC (thermal-hydro-mechanical and chemical) processes for a moderate temperature range of up to 100 °C, data for a higher temperature range are limited. Understanding bentonite behavior and coupled THMC processes under higher temperatures (e.g., up to 200 °C) could allow for a more economic repository design and would expand the data and knowledge base for more reliable modeling. In this study, a bench-scale experiment was conducted in a compacted bentonite column experiencing both heating up to 200 °C in the center and hydration from a sand-clay boundary surrounding the column. During the experiment run for 1.5 years, frequent X-ray computed tomography (CT) scanning of bentonite provided insights into the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. After the experiment, a comprehensive post-dismantling characterization of bentonite samples was conducted. Results showed that the bentonite hydration was axi-symmetrical despite the initial heterogeneity due to packing, confirming the ability of bentonite to seal fast flow/transport paths. Further, compared to a non-heated control experiment, the heated column showed greater CT density variations along the radial distance, indicating that homogenization of bentonite might be more difficult if a temperature gradient is maintained in the repository. Precipitation of an anhydrite layer occurred in the inner hot zone, pointing to potential concerns about salt precipitation causing canister corrosion. Ultimately, the experiments provided a high-resolution window into the strongly dynamic and coupled behavior of bentonite exposed to heating, hydration and swelling, which will be valuable for improving modeling of coupled processes, especially for the early state of a HLW repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Computational Methodology to Simulate Pyrometallurgical Processes in a Secondary Lead Furnace

Pyrometallurgical recovery of nonferrous metals involves a combination of thermally intensive transformations during exothermic gas-phase reactions, endothermic decomposition of solid charge, and melting of simpler solids. In the recovery of secondary lead, simultaneous thermal effects in a reverberatory-style furnace cause a melt pool to accumulate at the bottom, with lighter solids (slag) floating above and gaseous products from decomposition of the charge diffusing through the gas–slag interface. Species from oxy-fuel combustion of natural gas, species profiles from smelting reactions, and the formation of a melt pool consisting primarily of lead are simulated via a time-averaged formulation. Predictions of outflow are compared with preset inflow profiles to ensure conservation of mass. Thermal profiles for solid, liquid, and gas phases are presented by species. A novel method is implemented to model the latent heat of fusion using a heterogeneous chemical reaction. The simulation is conducted in Simcenter STAR-CCM+ v. 16.02.009-R8.

Rao, Vivek↗

DOE OSTI2020

Determining Design Criteria for Feeding Biomass into a Fluidized Bed Using a Feed Screw

Many pathways have been discovered for the thermochemical conversion of biomass to products. These pathways often show very promising results on the benchtop scale and could be economical with the appropriate scale up. However, many of the fundamental challenges of scaling up thermochemical reactions are not well understood. One primary problem area that rarely receives much upfront attention is feeding large quantities of biomass into thermochemical reactors. As biomass transitions from ambient to reactor conditions, it undergoes chemical and physical changes that affect its flowability. Often these changes lead to agglomeration of feed particles within the transfer screw, leading to large clumps of feed entering and disrupting the fluidized bed, or more critically leading to a plug in the feed screw causing costly downtime and maintenance to remove. This transitional phase must be better understood for scale up of conversion technologies to be successful. The Feedstock-Conversion Interface Consortium (FCIC) is a collaboration of eight National Laboratories. FCIC is currently conducting research on screw feeding of biomass into thermochemical fluidized bed reactors. Heat transfer models have been constructed for the screw/reactor interface. Mechanical forces on biomass within screw feeder are being determined, as well as how those forces change by varying feed screw design. Feedstock physical and chemical properties are then being measured along these temperature and mechanical force gradients. This work will better inform feed screw design for scale up, increasing the success rate of future biorefineries. Discussion will include research approach, testing, and preliminary results.

09 BIOMASS FUELS↗

DOE OSTI2021

Thermodynamic Modeling of Heat Engines Including Heat Transfer and Compression–Expansion Irreversibilities

In this paper, a thermodynamic model based on an endoreversible engine approach is developed to analyze the performance of heat engines operating under different thermodynamic cycles. The model considers finite heat transfer rate, variable heat source and sink temperatures, and irreversibilities associated with the expansion and compression. Expressions for the maximum power and efficiency at maximum power output are obtained as a function of hot and cold reservoir temperatures, the equivalent isentropic efficiency of compression and expansion components, and the effective conductance ratio between heat exchangers. In all cases, the Curzon–Ahlborn efficiency is retrieved at constant reservoir temperatures and neglected compression–expansion irreversibilities. The proposed model allows assessing the effect of isentropic efficiencies and heat exchanger design and operation characteristics for different thermodynamic cycles.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

Additive manufacturing and testing of a ceramic heat exchanger for high-temperature and high-pressure applications for concentrating solar power

Heat exchangers with excellent corrosion and oxidation resistances are essential for the next-generation concentrating solar power (CSP) plants using a molten salt heat transfer fluid and a supercritical CO 2 Brayton power cycle that can operate at temperatures > 700 °C for higher efficiencies. Techniques were developed for additively manufacturing ceramic materials for applications involving high temperature, high pressure, and high corrosion resistance as needed for the CSP application. Based on a previous ceramic heat exchanger design with its cross-section geometry of the flow channels optimized for heat transfer and mechanical strength, lab-scale prototype heat exchangers, with integrated headers and incorporated flow channel dimensional compensations, were fabricated by using the binder jetting process. The printed prototype heat exchangers were successfully densified through the processes of liquid polymer infiltration and pyrolysis. In addition, thermophysical properties of the densified silicon carbide parts were measured to provide necessary design information. Experimental heat transfer testing of the lab-scale prototype heat exchanger was conducted, and the experimental data agreed reasonably well with the simulated results. Furthermore, this agreement provides validation of the simulation models for their applicability to the development of a full-scale ceramic heat exchanger.

14 SOLAR ENERGY↗

DOE OSTI2021

In-Situ Ion Irradiation to Add Irradiation Assisted Grain Growth to the Marmot Tool

In the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the MARMOT mesoscale fuel performance tool is used to inform the development of mechanistic materials models for the BISON fuel performance tool. The grain size of the fuel has a large impact on its performance, directly impacting heat conduction, fission gas release, creep, and fracture. Thus, atomistic and mesoscale MARMOT simulations have been used to investigate grain boundary migration and grain growth in UO 2 . However, the current grain growth models in MARMOT does not consider the effect of irradiation on causing grain growth. This research project experimentally studied and implemented irradiation effects on grain growth of UO 2 into MARMOT.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Effectiveness of street trees in reducing air temperature and outdoor heat exposure in Las Vegas

Urban greening and especially tree planting strategies are being widely planned and implemented to mitigate urban overheating and thermal stress in many urban areas. However, the effectiveness of these strategies depends on multiple factors, including urban morphology, environmental conditions, and tree characteristics. This study investigates the effects of street tree planting strategies on air temperature and outdoor heat exposure in Las Vegas by combining an urbanized mesoscale climate model coupled to a multi-layer street tree model (Weather Research and Forecasting-BEP-Tree; Δ x = 900 m) with a microscale pedestrian heat exposure model (Temperatures of Urban Facets for Pedestrian; Δ x = 1 m). A series of simulations are conducted for July and August 2022. Large city-wide increases of a drought tolerant tree species cool air temperature mainly during the nighttime (up to 1.5 °C), with daytime effects being limited due to leaves shedding sensible instead of latent heat as stomata close in response to high vapor pressure deficits. Increased evaporative cooling is achieved with a different tree species (double at night, and reaching 0.4 °C during the day), but water requirements increase threefold. Despite their relatively small effect on air temperature during the day, trees provide significant shade by intercepting solar radiation, reducing mean radiant temperature (up to 16 °C) and enhancing outdoor thermal comfort, a major benefit of street trees in hot arid climates. The nighttime cooling of trees and the daytime reduction of radiant loading show potential to reduce the heat-related health impacts. Our results highlight the need to evaluate the effects of street trees on a case-by-case basis.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

U-rich U-Mo Solidus, Liquidus, Enthalpy, and Thermal Diffusivity

Several material properties have been studied as a function of composition as part of the USHPRR effort to use U-10Mo for LEU monolithic fuel foils, with special attention to areas where simulation outcomes can be improved. The liquidus and solidus phase boundaries are important for engineering design of processing parameters and for predicting the degree of microsegregation within initial castings. The numerous phase diagram constructions available in the literature differ significantly in their solid/liquid boundary definitions, therefore, an experimental strategy was employed in order to resolve some of the discord. That strategy focuses not on proving any singular point, but instead to provide data with highly characterized uncertainty. The data presented here is the most consistent with the U-rich liquidus and solidus boundaries provided by H. Okamoto (2012), Berche et al. (2011), and Mardon et al. (1959). The recently suggested change in the peritectic reaction temperature (to 1302 K) suggested by Berchè and adopted by Okamoto is also supported, but with lower confidence. The enthalpy of fusion is a basic characteristic of the solidification reaction property for which is required for both research-oriented predictions of behavior and development of applied design strategies for casting engineering. It has been determined that an admixture model for the heat of fusion between U BCC and Mo BCC is a reasonable assumption within the composition region of interest, though limitations of the measurement accuracy impacts the ability to support any trends more detailed than this baseline behavior. A limited number of thermal diffusivity measurements were also performed into the gamma region to provide additional information for solidification modeling and simulations as well as a check on the component properties of thermal conductivity, density, and heat capacity.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Modeling ammonia-fueled co-flow dual-channel protonic-ceramic fuel cells

This paper reports the model development for a dual-channel protonic-ceramic fuel cell (PCFC) operating on ammonia fuel. The model considers the coupled interactions of several physical and chemical processes, including three-dimensional heat conduction within the bipolar plates and the membrane-electrode assembly (MEA), one-dimensional flow within the fuel and air channels, detailed heterogeneous catalytic reactions within the porous composite anode structure, Butler–Volmer representation of the charge-transfer chemistry, and Nernst–Planck transport of three charged defects (protons, oxygen vacancies, and small polarons) within the dense electrolyte membrane. The membrane-electrode assembly is composed of a Ni-BCZYYb (BaCe 0.7 Zr 0.1 Y 0.1 Yb 0.1 O 3-δ ) anode, a BCZYYb electrolyte membrane, and a BCFZY (BaCo 0.4 Fe 0.4 Zr 0.1 Y 0.1 O 3-δ ) cathode. Chemical and physical parameters for the MEA model are established using previously published button-cell data. One aspect of the study is to investigate the partial ammonia decomposition upstream of the fuel cell. Such fuel cracking increases the H 2 content of the fuel entering the PCFC, which may have benefits. However, endothermic ammonia pyrolysis within the composite anode structure assists with thermal control of the cell. The dual-channel model can be considered as the unit cell of a full fuel-cell stack.

25 ENERGY STORAGE↗

DOE OSTI2025

Improved Weld Residual Stress Modeling System in BlackBear

This report presents enhancements to the MOOSE-based BlackBear application aimed at improving its capability to simulate welding and other thermo-mechanical manufacturing processes. Two primary avenues of improvement are pursued. First, to enhance user accessibility, we introduce a centralized default block restriction mechanism that ensures coverage checks are performed within user-specified default blocks. This default setting is applied consistently to all block-describable objects, such as variables, kernels, and more. In addition, we develop a modular action for moving heat source simulations, which integrates path file parsing, subdomain modification, and heat source kernel enforcement into a single, streamlined configuration. Second, to improve solver robustness, we implement an alternative method for assigning initial conditions to the updated active domain during the simulation, thereby enhancing convergence behavior. To validate the framework, we design and conduct several benchmark simulations, including heat conduction with progressive material addition, linear elasticity with time-dependent material deposition, and viscoplasticity model with isotropic hardening under similar conditions. Finally, we demonstrate the effectiveness of the proposed framework through large-scale thermo-mechanical welding simulations in both two and three dimensions.

42 ENGINEERING↗