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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Non-dimensional performance and safety parameters for heat pipes

The use of heat pipes in safety-critical systems such as nuclear microreactors dictates the development of generalized, practical, scalable performance and safety parameters. Traditional dimensional metrics, while informative, lack the universality required for comparative analysis across varying designs and operating regimes. Here, this work introduces a comprehensive set of non-dimensional parameters to characterize heat pipe performance and safety, including capillary performance, effective thermal conductivity, response time, exergetic efficiency, allowable temperature gradients, allowable rate of temperature change, priming coefficients, and factor of safety. A reference heat pipe design representative of microreactor applications was analyzed via the developed parameters using both traditional analytical models and Sockeye simulations under transient and steady-state conditions. Sodium, potassium, and water were evaluated as working fluids to demonstrate the applicability of the framework across a broad temperature range. The proposed non-dimensional parameters effectively captured key thermal-hydraulic behaviors and safety concerns, as was demonstrated via Sockeye simulations. This framework supports the development of design optimization strategies, operational protocols, and safety assurance practices for advanced reactor systems and other high-reliability applications.

42 - ENGINEERING↗
DOE OSTI2022

Application of the Redlich-Kister expansion for estimating the density of molten fluoride psuedo-ternary salt systems of nuclear industry interest

The development of next-generation molten salt reactors relies on accurate knowledge of the thermophysical properties of the candidate coolant and fueled molten salts. These thermophysical properties include density, viscosity, thermal conductivity, and heat capacity. Because of difficulties in measuring thermophysical properties of molten salts, there are many gaps in the current state of thermophysical property knowledge of these salts, particularly those that contain actinides or beryllium. Therefore, leveraging modeling techniques to estimate unknown molten salt thermophysical properties and guide future experimental measurements has high value for the nuclear industry. Here, the densities of molten fluoride pseudo-ternary salt systems, which are of interest to the nuclear industry, were estimated using Redlich-Kister expansion and Muggianu interpolation techniques. The pseudo-ternary systems considered for estimation in this study were NaF-LiF-ZrF 4 , LiF-BeF 2 -ZrF 4 , LiF-BeF 2 -ThF 4 , NaF-LiF-BeF 2 , NaF-KF-BeF 2 , NaF-ZrF 4 -UF 4 , and NaF-BeF 2 -UF 4 . This Redlich-Kister estimation approach accounts for nonideal mixing behavior based on pseudo-binary subsystem interaction parameters determined from experimentally measured pseudo-binary system density data sets. The Redlich-Kister estimation was compared with the method of additive molar volumes, which assumes ideal mixing. Additionally, the Redlich-Kister approach was used to determine previously unknown binary and ternary interaction parameters based on experimentally measured density data sets for select pseudo-ternary salt systems. The results of this study show improvement in density estimation using the Redlich-Kister approach for all systems considered compared with estimation by additive molar volumes. Furthermore, this analysis allowed for the estimation of nonideal density behavior in experimentally unstudied ZrF 4 -UF 4 and BeF 2 -UF 4 , as well as the quantification of ternary interaction in NaF-LiF-ZrF 4 , NaF-BeF 2 -UF 4 , and NaF-ZrF 4 -UF 4 .

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

BaCu 4/3 Si 2/3 P 2 and BaCu 2–( x + y ) Zn x Si y P 2 : Expanding the Semiconducting Landscape in the ThCr 2 Si 2 -Type Family

ThCr 2 Si 2 -type layered materials are a large family of compounds with applications ranging from thermoelectricity to magnetism, with the vast majority of the members exhibiting metallic behavior. Here, in this study, we synthesized a new group of materials with Cu-Si and Cu-Zn-Si square nets with the general formula BaCu 1.33 Si 0.67 P 2 and BaCu 2–(x+y) Zn x Si y P 2 (0 ≤ x ≤ 0.9; 0.3 ≤ y ≤ 0.7). Several synthesized compounds are charge-balanced semiconductors, which are rare in the ThCr 2 Si 2 family. All the reported compounds crystallize in the ThCr 2 Si 2 -type tetragonal I4/mmm space group, with Cu/Zn/Si jointly occupying the same 4d crystallographic site. In the Zn-free composition, BaCu 1.33 Si 0.67 P 2 , Ba, and P each occupy a single crystallographic site. The introduction of Zn results in the expansion of the unit cell and splitting the Ba atomic sites along the [001] direction. Such structural displacement of the Ba atoms was confirmed by the heat capacity measurements. Band structure and density-of-states calculations on ordered hypothetical structural models reveal either a small bandgap (∼0.2 eV) or semimetallic band structures. The compounds reported here exhibit high Seebeck coefficients and ultralow thermal conductivity, making them promising candidates for the development of thermoelectric materials.

crystal structure↗
DOE OSTI2025

Iso-cost-performance of thermal energy storage

As thermal energy storage (TES) systems gain increasing recognition as next-generation energy storage solutions, evaluating their techno-economic performance is crucial. This perspective analyzes the cost-performance of latent heat-based TES systems and introduces the concept of iso-cost-performance—maintaining a constant cost per unit energy ($\$$/kWh) despite material degradation. Using an empirical degradation model, we show that after 1200 thermal cycles, the thermal conductivity and volumetric energy density of a paraffin-based phase change material (PCM) decreased by 36.0% and 26.1%, respectively, resulting in a 31.2% reduction in the figure of merit (FOM) and, consequently, the system cost-performance. However, by introducing thermally conductive additives to enhance effective thermal conductivity, the TES system can recover its initial FOM, achieving iso-cost-performance operation. This framework quantitatively demonstrates how degradation-mitigation strategies—such as improving thermal conductivity—can offset material degradations and maintain long-term cost-effectiveness. Beyond PCM-based TES, the proposed FOM-based approach provides a generalized pathway for cost-performance optimization across various TES technologies.

25 ENERGY STORAGE↗
DOE OSTI2024

Comparison of Results between the Legacy and Refined RELAP5-3D Models of the High Temperature Test Facility in Exercises 1 and 2 of the HTTF Benchmark

Work conducted in FY23 identified that RELAP5-3D was capable of reproducing trends in HTTF data during experiment PG-27 but was incapable of reproducing measured values. The primary cause of this discrepancy between RELAP5-3D results and experimental data was hypothesized to be a distortion in power density that was introduced by the radial nodalization of the model. We further hypothesized that a new model would provide better results when compared to the experiments PG-27 and PG-29. Work this FY developed a new model that is better capable of capturing local heat generation rates and contains a representation of each 1/6 azimuthal sector of the core. We used this model to develop a new set of solutions to Exercises 1 and 2 of Problems 2 and 3 in the benchmark. In this report, we present the first comprehensive comparison of the results between the two models. We see that in Exercise 1A, which is common between problems 2 and 3, the results are similar, though the results from the new model show greater detail than those from the legacy model. In Problem 2 Exercise 1B and Problem 3 Exercise 1B, we see that heat removal is slower in the new model than the legacy model. Problem 3 Exercise 1C shows temperatures that are lower in most places in the new model than the legacy model, but the area with active heat generation has higher block temperatures in the new model than the legacy model. Problem 3 Exercise 1D further shows that long-term heat removal is lower in the new model. Problem 2 Exercise 1C demonstrated that the new model observes higher temperatures in the core regions than the legacy model, justifying the need to preserve the power density in HTTF. The validation of PG-27 and PG-29 also demonstrated the improved temperature agreement in the core regions, particularly with a calibrated model that implements an effective thermal conductivity for the core material. Overall, PG-27 models show reasonable to excellent agreement for steady-state temperatures and minimal to reasonable agreement for transients. PG-29 models showed minimal to insufficient agreement with the data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Modeling of Transport Processes in Liquid-Metal Fusion Blankets: Past, Present, and Future

The successful development of robust breeding blanket systems will strongly rely on computational tools for predicting the complex behavior of the electrically conducting liquid-metal (LM) breeder flowing in the complex-shaped blanket ducts in the presence of a strong plasma-confining magnetic field, volumetric heating, and tritium generation. Associated transport processes involve magnetohydrodynamic (MHD) flows, heat transfer, corrosion, and tritium transport. This paper is an overview of past and present efforts in the development, application, and verification and validation (V&V) of such computational tools. As a result of the ongoing campaign on V&V of computer codes for LM blankets, the international fusion community has identified several candidates that promise to become real blanket design and analysis tools in the near future. Among them are HIMAG, MHD-UCAS, COMSOL Multiphysics, ANSYS FLUENT, ANSYS CFX, and OpenFOAM. The progress, over the last decade, in the application of such codes in blanket studies is tremendous. This is illustrated with two examples for a dual-coolant lead-lithium (DCLL) blanket: (1) integrated computer modeling for the recently designed DCLL blanket in the United States and (2) application of the code MHD-UCAS to the analysis of PbLi flows and heat transfer in a generic DCLL blanket prototype at high Hartmann (Ha ~ 10 4 ) and Grashof numbers (Gr ~ 10 12 ). Here, this paper also presents an approach to the development of a new integrated computational tool called the virtual dual-coolant lead-lithium (VDCLL) blanket, which elaborates the existing U.S. MHD code HIMAG.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Structural Damage and Recrystallization Response of Garnet Crystals to Intense Electronic Excitation

Here, the susceptibility to irradiation-induced damage, the fundamental mechanism of the relaxation kinetics, and the corresponding recrystallization effect related to the cation radius ratio are comprehensively investigated for Y 3 Al 5 O 12 and Gd 3 Ga 5 O 12 garnet crystals under 645 MeV Xe 35+ irradiation with different fluences of 5 × 10 11 –3 × 10 12 ions cm –2 . Regarding different lattice distortion and swelling levels, the observed microstructure transformations to disordered and amorphous phases, and corresponding hillock dimensions, consistently confirm that Gd 3 Ga 5 O 12 has a higher susceptibility to radiation damage than Y 3 Al 5 O 12 . Combined with iTS model calculations, although Y 3 Al 5 O 12 has higher atomic temperature and energy deposition than Gd 3 Ga 5 O 12 under the same ion velocity and electronic energy loss, the relatively high thermal conductivity and specific heat coefficient of Y 3 Al 5 O 12 crystals enhance the conduction and dissipation of deposition energy, and Gd 3 Ga 5 O 12 , with a higher cation-radius-ratio (r A /r B ), is more easily damaged to amorphous phase due to the less favorable kinetics of ordering and recovery of a melted track region to the crystalline phase. Additionally, the significant bandgap modification in spectral ranges of 5.88–6.75 eV for Y 3 Al 5 O 12 and 4.83–5.41 eV for Gd 3 Ga 5 O 12 , and the enhancement of defect-assisted-related luminescence are achieved, providing a basis to design novel optoelectronic devices in microelectronics fields.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling of crystallization fouling on a horizontal-tube falling-film evaporator for thermal desalination

Scale formation within horizontal-tube falling-film evaporators is an important issue for thermal desalination, because of the deterioration of heat transfer performance and added maintenance costs. Predicting the fouling process is crucial for desalination plant design, operation, and maintenance. In this paper, a model is developed to predict the crystallization fouling process during seawater falling-film flow over a horizontal tube bundle. The scaling is estimated based on deposition theory, with couple heat transfer including the in-tube steam condensation, conduction through the tube wall and the scale layer, and falling-film evaporation. The spatial and temporal variations of temperature, heat transfer coefficient, and scale thickness are predicted, and the effects of variations in process parameters of steam, feed seawater, and tube properties on the scale thickness, evaporation rate, and heat transfer coefficient for falling-film evaporation are analyzed. Comparisons to existing experimental data show the scaling layer thickness falls in the general range of reported measurements, with a slight overestimation due to neglecting the delay of scaling onset under real conditions. The scale layer thickness increases dramatically as the in-tube steam pressure is increased or the seawater flow rate is decreased. In general, the scaling layer becomes thicker on the lower tubes in the tube bundle, due to the heat and mass transfer of the falling film. The effect of tube material on scaling appears mainly dependent on thermal conductivity; thus, using polymer rather than a stainless steel tube decreases scaling and evaporation rate 96% and 88%, respectively. Finally, this work has potential to guide thermal desalination plant design, material selection, and operating parameter optimization.

42 ENGINEERING↗
DOE OSTI2022

Magnetic Control of Soft Chiral Phonons in PbTe

We report PbTe crystals have a soft transverse optical phonon mode in the terahertz frequency range, which is known to efficiently decay into heat-carrying acoustic phonons, resulting in anomalously low thermal conductivity. Here, we studied this phonon via polarization-dependent terahertz spectroscopy. We observed softening of this mode with decreasing temperature, indicative of incipient ferroelectricity, which we explain through a model including strong anharmonicity with a quartic displacement term. In magnetic fields up to 25 T, the phonon mode splits into two modes with opposite handedness, exhibiting circular dichroism. Their frequencies display Zeeman splitting together with an overall diamagnetic shift with increasing magnetic field. Using a group-theoretical approach, we demonstrate that these observations are the result of magnetic field-induced morphic changes in the crystal symmetries through the Lorentz force exerted on the lattice ions. Thus, our Letter reveals a novel process of controlling phonon properties in a soft ionic lattice by a strong magnetic field.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Empower Wall: Active insulation system leveraging additive manufacturing and model predictive control

Buildings are one of the largest energy consumers worldwide, using large amounts of energy during their construction and for climate control during operation. Active insulation systems (AIS) have been shown to reduce the energy needed for climate control in buildings by dynamically regulating the heat transferred between a building’s interior and exterior. Infrastructure-scale additive manufacturing (AM) has the potential to reduce the resources needed for building construction. Combining these two technologies into a single building envelope would create a path towards more sustainable buildings. A test was conducted for the Federal Energy Management Program (FEMP) Energy Exchange training and trade show, in August 2021, to investigate a new building envelope design, termed the Empower Wall, that utilized an AIS and was constructed using AM. Model predictive control was implemented to manage operation of the Empower Wall in concert with the existing HVAC system. Finally, the prototype system demonstrated that the Empower Wall lowered total energy consumption and reduced the cost of energy used.

36 MATERIALS SCIENCE↗
DOE OSTI2020

GaAs growth rates of 528 μm/h using dynamic-hydride vapor phase epitaxy with a nitrogen carrier gas

Herein we demonstrate record GaAs growth rates approaching 530 μm/h using nitrogen carrier gas and 400 μm/h using hydrogen carrier gas in a dynamic-hydride vapor phase epitaxy reactor. We measured root mean square surface roughness below 1 nm using a 1 μm × 1 μm atomic force microscopy scan for GaAs growth rates up to 483 μm/h using a nitrogen carrier gas and 400 μm/h using a hydrogen carrier gas. We performed computational fluid dynamics modeling to study the effect of the carrier gas choice on the thermal profile within the reactor and how that influences the degree of AsH 3 decomposition. The modeling suggests that the lower thermal conductivity of the nitrogen carrier gas minimizes the amount of AsH 3 that thermally decomposes before reaching the wafer surface and the heavier atomic mass decreases the likelihood that AsH 3 will reach the heated reactor walls, leading to a growth rate enhancement relative to the hydrogen carrier case in a hydride-enhanced growth regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Transducerless time domain reflectance measurement of semiconductor thermal properties

The dominant techniques for laser-based thermal measurement, time-domain thermoreflectance (TDTR) and frequency-domain thermoreflectance, generally require that samples be coated with a metal transducer layer. The properties of the transducer and its interface to the sample are then important to the interpretation of results, and heat conduction in the transducer layer can limit measurements of ultra-thin samples or in-plane conductivity. Here, we demonstrate the feasibility of measuring semiconductor samples using standard TDTR instrumentation, without the use of a metal transducer. A mathematical model accounting for non-zero laser penetration depth as well as both thermal and charge carrier effects is presented, and the ability of the model to describe and fit to experimental data for silicon and germanium samples is demonstrated. Though at this stage measurement uncertainties are larger than for traditional TDTR, we anticipate that the use of this technique will expand the range of samples that can be measured with existing TDTR setups.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Urban Parameters Maricopa County 100m Grid Spacing

132 urban parameters for Maricopa County are provided in WRF-readable binary form at 100m grid spacing with index. These parameters allow regional weather modelers to account for the contribution of built structures in modeling the urban heat island. This data is also provided in CSV format. Data are generated using 2023 Model America Version 1 data (https://www.osti.gov/biblio/1774134) and will be used to represent the urban morphology in Maricopa County to the WRF model for experiments conducted for the Southwest Urban Corridor Integrated Field Laboratory (SW-IFL). Three files associated with these parameters are included in this repository: index, binary file 00001-01549.00001-01796, Maricopa_Parameters_03072024.csv . To use the binary file with the Weather Research and Forecasting (WRF) model, the binary file and the index file must be placed in their own directory in WRF_GEOG and accessed in the same way the National Urban Database and Access Portal Tool data (NUDAPT44) would be accessed. The name of the binary file must be in the format: 00001-02715.00001-01803. If, upon download, the binary filename has additional or different characters, please change them to match this format.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

CRADA Number NFE-18-07325 with TCPoly, Inc. (CRADA Final Report)

TCPoly is commercializing high performance composite materials for the 3D printing industry and the company focus is manufacturing high thermal conductivity plastics for 3D printing of heat exchangers, electronic cooling devices, and mold tooling. Through the Innovation Crossroads Program, I was able to collaborate with ORNL researchers to expand the company’s technology offering and revamp the company’s business model and focus.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of Decay Heat Sensitivity Analysis Capability in SCALE/ORIGEN

In this work, a decay heat sensitivity analysis capability was developed and implemented in the ORIGEN code of the SCALE nuclear modeling and simulation suite. This capability introduces improved numerical integration schemes, which overcome the challenges associated with accurately modeling the behavior of adjoint nuclide amounts during coarse time steps for both nuclide amount and decay heat sensitivity calculations. This capability significantly improves the accuracy of calculations without compromising computational efficiency compared to the existing method. Extensive verification was conducted for various benchmark problems, including a 238 Pu decay and an irradiation problem involving 135 Xe, evaluated with both coarse and fine time grids. The results show excellent agreement with reference direct perturbation solutions, reaffirming the computational accuracy of the newly proposed numerical integration methods. Furthermore, sensitivity analyses were performed for fission product inventories ( 147 Sm, 150 Sm, 155 Gd) in pressurized water reactor UO 2 and MOX fuel assemblies. These analyses demonstrated that the ORIGEN sensitivity analysis capability can capture detailed sensitivity coefficients and underlying physics in real applications. Additionally, a decay heat sensitivity analysis for high-assay low-enriched uranium fuel, including various initial 235 U enrichment and burnup points, highlights the extended capabilities of SCALE/ORIGEN in comprehensively assessing the factors influencing total decay heat. These advancements in ORIGEN offer valuable insights for reactor analysis, fuel design, and safety assessments, especially in the context of advanced nuclear fuel development and design changes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Electronic and thermal properties of the phase-change memory material, Ge 2 Sb 2 Te 5 , and results from spatially resolved transport calculations

Here, we report new insights into the electronic, structural, and transport (heat and charge) properties of the phase-change memory material amorphous Ge 2 Sb 2 Te 5 . Using realistic structural models of Konstantinou et al., (2019), we analyze the topology, electronic states, and lattice dynamics with density functional methods, including hybrid-functional calculations and machine-learned interatomic potentials. The Kohn–Sham orbitals near the Fermi level display a strong electron–phonon coupling, and exhibit large energy fluctuations at room temperature. The conduction tail states exhibit larger phonon-induced fluctuations than the valence tail states. To resolve transport at the atomic scale, we employ space-projected electronic conductivity and site-projected thermal conductivity methods. Local analysis of heat transport highlights the role of filamentary networks dominated by Te, with Sb and Ge making progressively smaller contributions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A Review of Modeling Approaches for Predicting Frost Growth and Defrosting on Tube-Fin Heat Exchangers: Preprint

Frost formation and growth on the evaporator surface is a common process that deteriorates the air-refrigerant heat transfer and restricts airflow. This degrades the performance of the vapor compression system by increasing temperature lift and air-side pressure drop. To accurately predict these effects during coil frosting, as well as the energy use and duration of the defrost process, there is a need to estimate the heat and mass transfer, momentum transport, and solid-liquid and liquid-vapor phase change. Therefore, in the past few decades, continuous effort has been made to model frosting and defrosting processes using approaches ranging from empirical correlations to computational fluid dynamic models. To provide a clearer overview for researchers, engineers, and manufacturers in this field, this paper provides a comprehensive literature review for frosting and defrosting models. The paper begins with theoretical background of frost formation and defrost processes, and then reviews the common modeling approaches in literature and their underlying assumptions when trying to account for various physical phenomenon. Based on the literature review, the most critical modeling effort for frost formation is the determination of frost densification rate and frost growth rate. Various methods to predict these two parameters are reviewed. Empirical correlations commonly used for frost density and thermal conductivity are presented and compared. For the defrost process, various multi-stage models have been proposed with different assumptions. Some assume the presence of air gap between the tube wall and the frost, while others consider the melted frost flow due to gravity. We also review physics-based and empirical approaches to integrate defrost models into heat pump models. We conclude by identifying research gaps and providing recommendations.

defrost↗
DOE OSTI2025

Comparative performance assessment of air-source and ground-source heat pumps using CO₂ and R-410A with water well integration: A simulation study

This study investigates the performance and feasibility of heat pump systems for residential space heating in cold climates, with a particular focus on ground source heat pumps (GSHPs) with water wells. Four configurations are modeled and compared, a CO 2 air-source heat pump (ASHP), an R-410A ASHP, a CO 2 GSHP with water well integration, and an R-410A ground-source heat pump with water well integration. System simulations are conducted at both equipment and whole-building levels, followed by a nationwide analysis across ten representative cities using EnergyPlus. Results show that water-well-coupled GSHPs maintain approximately 87 % of their efficiency and 85.5 % of heating capacity as ambient temperature drops from 47 °F/8.3 °C to −15 °F/-26.1 °C, whereas ASHPs retain only 55 % efficiency and 44.5 % capacity. R-410A systems achieve higher efficiency, while CO 2 systems provide significant environmental advantages, including reduced risks of groundwater contamination from refrigerant leakage. On average, the CO 2 GSHPs deliver 35–45 % heating energy savings relative to ASHPs and demonstrate strong potential as a replacement for gas furnaces in cold climates. This work represents a systematic comparative assessment of CO 2 and R-410A air-source and ground-source heat pumps with water well integration. In conclusion, the findings highlight the technical viability, environmental benefits and deployment potential of CO 2 GSHPs, offering a pathway toward efficient and sustainable residential heating solutions in diverse U.S. climates.

CO 2 GSHP↗