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 253 records · Page 14

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↗

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↗

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↗

Model-Based Investigation of Multi-Fault Interactions and Performance Degradation in Residential Heat Pump Systems

Faults in heat pump systems can significantly degrade performance, reduce efficiency, and accelerate component wear, leading to higher operating costs and maintenance demands. While numerous studies have investigated the impact of individual faults, the interactions between multiple concurrent faults remain insufficiently understood, despite their common occurrence in real-world operation. This study conducts a comprehensive simulation analysis of multiple simultaneous faults in a vapor compression heat pump using a validated heat pump design model (HPDM) tool. Detailed component-level modeling methods are implemented to examine performance sensitivity under combinations of refrigerant flow and heat exchanger faults. The results reveal complex fault interactions that can mask or amplify system deviations, challenging conventional diagnostic approaches. Findings from this work provide meaningful insights for the development of more robust fault detection and diagnosis algorithms, supporting improved reliability and energy efficiency in next-generation heat pump technologies.

Hu, Yifeng [ORNL] (ORCID:0000000242875185)↗

Residual Stress Modeling with Phase Transformation for Wire Arc Additive Manufacturing of B91 Steel

Wire arc additive manufacturing (WAAM) is an energy-efficient manufacturing technique used for near-net-shape production of functional industrial components. However, heat accumulation during deposition and the associated mechanical and metallurgical changes result in complex residual stress profiles across the cross section of the fabricated components. These residual stresses are detrimental to the service life of the components. In this study, sequentially coupled thermomechanical analysis of WAAM B91 steel is conducted to quantify the residual stress variation across the component. The thermomechanical analysis includes a transient heat transfer model and a static stress model that incorporates the transformation-induced plasticity due to martensitic phase transformation. Here, the experimentally calibrated heat transfer model mirrors the temperature variation of the system during the deposition. The results from the stress model are validated via x-ray diffraction measurements, and the numerical results are in good agreement with the experimental data.

20 FOSSIL-FUELED POWER PLANTS↗

Thermal coupling-decoupling mechanism of heat transfer across van der Waals interfaces in n-eicosane

Progress toward improving thermal transport properties of n-eicosane, a promising phase change material for thermal energy storage near human body temperature, is hindered by lack of a general theoretical framework to model thermal interfacial conductance (TIC) of the van der Waals (vdW) force bonded molecular interfaces, which have little compositional/structural mismatch. Combining molecular dynamics simulations at temperatures up to the melting point and a time-domain analysis of interfacial heat currents in individual molecules, we unveil that (1) the heat flux across a molecular vdW interface has an alternating coupling-decoupling pattern at the single molecule level due to thermal fluctuations, and (2) the net heat transfer during the thermally coupled periods is orders of magnitude larger than that during the thermally decoupled periods. Furthermore, we demonstrate that a vdW interface's TIC correlates almost linearly with two atomistic-scale coupling-decoupling parameters: αc and σc, that represent the duration and strength of the thermal coupling, respectively. This thermal coupling-decoupling mechanism proposes a new paradigm that explicitly accounts for the effects of structural flexibility and dynamical fluctuation on heat transfer across molecular interfaces, and will likely provide new insights into modeling thermal transport properties in a wide-ranging soft materials, including self-assembled monolayers, paraffins or lipids.

42 ENGINEERING↗

Heat transfer model for moving packed-bed particle-to-sCO 2 heat exchangers integrated with metal foams

Particle-to-supercritical carbon dioxide (sCO 2 ) heat exchangers (HXs) play a vital role in coupling heat transfer fluid (HTF) from high-temperature thermal receivers to power cycle working fluids (WF). Heat transfer enhancement is essential for adopting particle-based moving packed-bed heat exchangers (MPBHXs) in next-generation thermal energy storage (TES) systems, as MPBHXs usually exhibit low particle bed-to-wall heat transfer coefficients. High-porosity metal foams have shown their effectiveness in heat transfer enhancement. Here, this work presents a continuum heat transfer model to demonstrate the heat transfer enhancement in MPBHXs when the particle bed channel is filled with high porosity metal foams compared to open channel MPBHXs. The presence of metal foams increases the effective thermal conductivity in the particle channel and enhances the interstitial heat transfer coefficient between the moving particle bed and the stationary metal foams. The present model considers coupled two-dimensional (2D) heat transfer in the particle channel with metal foams and 1D heat transfer in the sCO 2 channel and the dividing wall. The temperature profiles for the particle bed, metal foam, dividing wall, and sCO 2 stream, as well as the local heat flux profiles between those, are studied in detail for various foam porosities, from which the particle bed-to-wall heat transfer coefficient, overall heat transfer coefficient, and total heat transfer rate for the MPBHX are determined. The effects of major MPBHX design and operating parameters on the particle bed-to-wall heat transfer coefficient and overall heat exchange capacity are thoroughly examined; thus, the MPBHX performance improvement with metal foams is quantified for each case. The present heat transfer model can provide valuable insights into the metal-foam MPBHX design and optimization, scale-up, and operating parameters selection.

42 ENGINEERING↗

Building Model Calibration: Validation of GridLAB-D Thermal Dynamics Modeling

This report investigates how well GridLAB-D’s house model characterizes the thermal dynamics of buildings given the overpredicted diurnal electric load swing observed in the Distribution System Operation with Transactive (DSO+T) study. This study seeks to validate GridLAB-D’s house model by calibrating it to data from the well-instrumented Pacific Northwest National Laboratory Lab Homes in Richland, WA. The datasets chosen included multiple pre-cooling and pre-heating testing periods where the indoor air temperature was allowed to float over a multi-hour length of time to mimic diurnal behavior. The multi-season calibrations were evaluated by comparing the heating/cooling electric power, indoor air temperature, and the rise/decay time during temperature float periods with Lab Homes data. The default GridLAB-D assumptions for latent load fraction, air heat capacity, mass heat capacity, window-to-wall-ratio, overall envelope conductance (assumed as NORMAL thermal integrity level), and solar heat gain coefficient were compared with the calibrated model to confirm when the default assumptions were adequate and the impact of calibration on the accuracy of modeling the thermal dynamics of homes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Continuous stabilization of polyacrylonitrile (PAN) - carbon nanotube (CNT) fibers by Joule heating

We report stabilization is one of the important processes used to convert polyacrylonitrile (PAN) fibers into carbon fiber and is usually driven by energy-intensive and time-consuming convective heating. A continuous and energy-efficient alternative is desirable. Here, a continuous Joule heating process has been used for the first time to stabilize PAN fibers containing carbon nanotubes (CNT). Single-component PAN/CNT fibers and bi-component fibers with core-sheath geometry were fabricated through the gel spinning method. Joule heating was conducted continuously on a 200 filament bundle for single-component fibers and on bundles with up to 6000 filaments for the bi-component fibers. A model of heat transfer and thermodynamics occurring during Joule heating is developed and utilized to understand the reductions in stabilization time and energy consumed when compared to convective heating. The minimum energy needed for stabilization through Joule heating is estimated to be less than 1% of the energy used through convective heating.

42 ENGINEERING↗

Site specific porosity-thermal performance correlations in neutron irradiated U-10Zr fuel

In this work, we present a site-specific three-dimensional analysis of porosity evolution in neutron-irradiated U–10Zr metallic fuel using high-resolution synchrotron X-ray tomography. Focused ion beam cubic lift-outs from four radial positions—fuel center, middle, edge, and fuel-cladding interaction (FCCI) zones—were imaged, reconstructed, and segmented to quantify pore volume, density, morphology, and pore connectivity. Porosity increased modestly from 5.5% at the center to 10.5% at the edge, yet pore number density increased by over two orders of magnitude in the fuel portion near the FCCI interface (from 4.8 x 10 4 to 6.0 x 10 6 pores/mm 3 ). Morphological classification revealed a progression from small spherical pores at the center to equiaxed and tortuous networks at the periphery, with enhanced orientation along the radial direction. Connectivity and permeability analysis reveal that the FCCI region maintains dense, highly interconnected pores despite a moderate volume fraction, enabling rapid fission gas transport and lanthanide migration. Effective thermal conductivity models incorporating sodium logging confirm that these site-specific pore features critically influenced heat transport during reactor irradiation. These findings demonstrate that pore topology—not just porosity fraction—controls thermal performance and thermal pathways in U–10Zr fuels, impacting fuel thermal performance in a reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Computational and Experimental Study of Film-Cooling Effectiveness With and Without Downstream Vortex Generators

Downstream vortex generators that involve a pair of rectangular plates arranged in an open V-shape placed just downstream of each film-cooling hole were shown to create flow and vortical structures that entrain lifted film-cooling flow back to the surface and increase its lateral spreading on the surface (GT2020–14317). In this study, computations and measurements were performed to examine the flow mechanisms this vortex generator induces to improve film-cooling effectiveness of a flat plate with the cooling flow emanating from one row of inclined holes. Parameters studied include blowing ratio (BR = 0.75 and 1.0) and temperature ratio (TR = 1.07 and 1.9). Here, the computational study is based on steady Reynolds-averaged Navier–Stokes (RANS) closed by the shear-stress transport (SST) turbulence model with and without conjugate analysis. The experimental study was conducted by using a conjugate heat transfer test rig with a plenum, where cooling flow is introduced. Measurements made include velocity and temperature profiles upstream and downstream of the film-cooling holes as well as the temperature at several locations on the hot and cold sides of the film-cooled flat plate. The computational study was validated by comparing computed results with those from measurements at BR = 0.75 and 1.0 and TR = 1.9. Computational and experimental results are presented to show the effects of BR and TR on the flow structures and how those structures improve the effectiveness of film cooling with and without the downstream vortex generators and with and without conjugate heat transfer.

42 ENGINEERING↗

Oakridge Geothermal Resources Assessment

Through the Communities Local Energy Action Program (LEAP), NREL provided technical assistance to a coalition of stakeholders from Oakridge, Oregon. The coalition included community-based non-profit organizations, the city government, the local utility, and others. Technical assistance provides analysis and information to support Oakridge stakeholders with their goals to increase energy reliability and resilience in the community, while promoting economic development. The Geothermal Resource Assessment provides technical analysis related to the geothermal resource in the Oakridge area, possible technology applications, and information related to further geothermal exploration. Technical analysis includes a review of a previous geothermal study conducted for Oakridge, thermal modeling to estimate reservoir temperatures, and technoeconomic simulation of a combined heat and power application in the Oakridge area.

15 GEOTHERMAL ENERGY↗

A model to separate conduction and radiation in high temperature laser flash measurements for semi-transparent materials

Laser flash analysis (LFA) is a widely used technique in the scientific and industrial communities for measuring the thermal diffusivity and conductivity of materials. However, traditional fitting methods based solely on heat conduction do not account for radiation heat transfer in transparent and semi-transparent materials, which becomes significant at high temperatures. In this study, we have developed a coupled radiation-conduction model that considers the optical properties of materials to accurately extract thermal diffusivity from LFA measurements. We conducted LFA measurements on fused silica glass samples of varying thicknesses and within the temperature range of 25 to 1000 °C and fitted the data with the model. Since the coupled model requires optical properties as input parameters, which are often unavailable for new materials, we have also proposed a modified transparent model. Compared to the coupled model, results from the modified transparent model generally overestimate thermal conductivity, but the errors are acceptable in silica glass with small thickness. Here, we also investigated the applicability of the modified transparent model with materials of different optical properties.

14 SOLAR ENERGY↗

CFD modeling of crystallization fouling with CO 2 desorption incorporated for a falling-film evaporator in thermal desalination

Seawater flowing around horizontal tubes in falling-film evaporators is a common configuration for thermal desalination. Heat is transferred from in-tube steam condensation to the shell-side seawater by conduction through the tube wall and scaling layer, and conduction and convection to the evaporating liquid film. CO 2 is simultaneously released from seawater and mixed with the produced vapor. Fouling coupled with CO 2 desorption has deleterious impacts on both evaporation and condensation. The evaluation of spatiotemporal dependent crystallization fouling and CO 2 desorption is vital to selecting optimal operating conditions. Here, in this work, a CO 2 desorption model is integrated into CFD modeling for predictions of coupled heat and mass transfer, to understand and predict scale formation and CO 2 desorption with local and transient profiles of temperature, carbonate species concentrations, pH, and total alkalinity. The porosity of the scale layer is experimentally determined as 67.7% and invoked in the calculation of effective thermal conductivity. The simulation results reveal that high steam temperatures increase seawater total alkalinity and accelerate scale formation. The scale thickness on the bottom tube reaches 0.30 μm and 5.24 μm for steam temperatures of 60 °C and 80 °C, respectively. High salinity leads to a large CO 2 desorption rates. The CO 2 desorption rate increases 43% when the seawater salinity increases from 35 to 55 g/kg. The effects of operating conditions on carbonate speciation and pH have been compared and analyzed. This model can serve as a comprehensive tool for the design of thermal desalination systems and optimal operation.

42 ENGINEERING↗

Coupled Experimental and High-Temperature Discrete-Element Method Modeling Studies of Aluminosilicate Particle Handling in Concentrated Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrated solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrated solar power environments. Coupled experimental and numerical studies of aluminosilicate particles in rotary kilns and dense particle curtains were performed for bulk temperatures up to 1073 K. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Temperature, spatial, and velocity profile data were extracted from experimental runs using embedded K-type thermocouple probes and particle image velocimetry techniques. Experimental and numerical studies were compared using spatial temperature profiles, velocity fields, and shape profiles of the bulk, granular flows. Numerical models were developed using commercially available discrete element method modeling software, Aspherix®. Existing Aspherix® functionality was expanded by introducing coupled radiative exchange modeling tools. The laboratory-scale rotary kiln was developed to investigate the steady-state heat and mass transfer performance of aluminosilicate particles based on particle type, bulk handling temperature, and wall roughness. The rotational speed of the rotary kiln was varied to control the relative impact of collisional and frictional effects upon the granular flow behavior. Heat and mass transfer performance was categorized based on the Froude number and the observed flow regimes of slipping, rolling, cascading, and centrifuging. Coupled discrete element method modeling studies were used to evaluate the effects of temperature-dependent, particulate mechanical properties upon bulk flow behavior and upon the relative effects of radiative, advective, and/or conductive heat transfer. A high-temperature (< 1073 K) falling particle curtain was similarly fabricated to investigate the heat and mass transfer performance of aluminosilicate particles in particle handling situations dominated by inter-particle collisions. The impact of particle type, flow preheat temperatures (< 1073K), and bulk mass flow rates were investigated upon the particle curtain shape, temperature, and velocity profiles. Coupled discrete element method modeling studies were performed to evaluate the varying impact of temperature-dependent, particulate mechanical properties on the bulk flow behavior and the temperature profile of the particle curtain.

14 SOLAR ENERGY↗

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system↗

Bentonite buffer under high temperature: laboratory experiments and coupled process modeling

Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.

Zheng, Liange↗