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

Coupling thermal energy storage with a thermally anisotropic building envelope for building demand-side management across various US climate conditions

Here, the thermally anisotropic building envelope (TABE) is a novel active building envelope that enhances energy efficiency and thermal comfort in buildings by transferring heat and cold between building envelopes and hydronic loops. When coupled with thermal energy storage (TES) units, the TABE + TES enables the storage of both heat and cold energy captured by the TABE roof or exterior walls. This stored energy can be later released by the TABE floor for indoor heating and cooling, benefiting both the grid and the end user. This paper evaluates the merits of TABE + TES for building demand-side management across various US climate conditions, focusing on peak load shaving, annual energy savings, and cost savings under time-of-use (TOU) electric rate schedules. Simulations were conducted by integrating time-of-day–informed, rule-based control strategies in MATLAB, TABE components and TES units in COMSOL Multiphysics, and whole-building energy analysis in EnergyPlus. A case study using the US Department of Energy’s prototype single-family detached house model in Birmingham, Alabama; Los Angeles, California; Oak Ridge, Tennessee; and Denver, Colorado, showed that the TABE + TES system achieved (1) 70 % peak load shaving in Los Angeles and Denver and 20 % in Birmingham and Oak Ridge; (2) significant peak electricity savings of 351–497 kWh, reducing peak energy consumption by 38 %–78 %; and (3) annual heating cost savings of 0.79 $\$$/m2–1.17 $\$$/m 2 and cooling cost savings of 0.60 $\$$/m 2 –1.17 $\$$/m 2 using a normal utility rate or low-TOU rate. The benefits of employing the TABE + TES system are even more significant under high TOU rates.

25 ENERGY STORAGE↗
DOE OSTI2020

Impact of Grain Size on Performance Degradation of TREAT LEU

We argue that radiation damage induced degradation of thermal conductivity does not set a lower limit on fuel grain sizes for the low enriched uranium fuel design of the Transient Reactor Test Facility (TREAT). Earlier work reports that smaller grains cause a larger degradation of thermal conductivity than larger grains constraining the smallest feasible size of fuel grains. This work assesses TREAT’s transient performance in the presence of radiation damage. The difference between the two studies is in treating damaged and fresh graphite as serial (this work) or parallel (previous) heat resistors. We use a multiphysics model of TREAT fuel grains to compute the reduction in transient capability measured by the total deposited energy as a function of irradiation dose. We find that radiation damage has a negligible effect on energy deposition.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Multi-Objective Optimization of a Bi-Metal High-Temperature Recuperator for Application in Concentrating Solar Power

Supercritical CO 2 closed Brayton cycles are a major candidate for future power cycle designs in concentrating solar power applications, with high-temperature recuperators playing an essential role in realizing their high thermal efficiency. Here, printed circuit heat exchangers (PCHEs) are often chosen for this role due to their thermal-hydraulic and mechanical performance at high temperatures and pressures, all while remaining compact. However, PCHEs can be costly because of the high-performance materials demanded in these applications, and the heat exchanger internal geometry is restricted by their manufacturing process. Additively manufactured heat exchangers can address both of these shortcomings. This work proposes a modular bi-metal high-temperature recuperator with integrated headers to be produced with additive manufacturing. Beginning with existing PCHE channel geometries, a 1D heat exchanger model is developed. Then, multi-objective optimization is used to maximize the heat transfer effectiveness of a lab-scale device while limiting its size. Two distinct channel geometries emerge from the optimization. Optimal designs achieve up to 88% effectiveness with negligible pressure drop. Deterioration of effectiveness due to axial conduction of heat in the heat exchanger walls is found to be a notable problem for lab-scale PCHEs, and the optimal designs obtained here minimize its detrimental effects. A sensitivity analysis reveals that the effectiveness of the recuperator is much less sensitive to variation in mass flowrate in off-design operation when axial conduction is significant, while increasing the length of the device easily increases effectiveness.

14 SOLAR ENERGY↗
DOE OSTI2020

Utilizing Temperature and Brine Inflow Measurements to Constrain Reservoir Parameters During a Salt Heater Test

Brine availability in salt has multiple implications for the safety and design of a nuclear waste storage facility. Brine availability includes both the distribution and transport of brine through a damaged zone around boreholes or drifts excavated into the salt. Coupled thermal, hydrological, mechanical, and chemical processes taking place within heated bedded salt are complex; as part of DECOVALEX 2023 Task E this study takes a parsimonious modeling approach utilizing analytical and numerical one-dimensional simulations to match field measurements of temperature and brine inflow around a heater. The one-dimensional modeling results presented arrive at best-fit thermal conductivity of intact salt, and the permeability and porosity of damaged salt of 5.74 W/m∙K, 10 –17 m 2 , and ≈ 0.02, respectively.

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

hydro-hifem v.1.0

SAND2021-3180 O The hydro-hifem 1.0 software is used to simulate transient physical processes in large geologic environments consisting of complex fractures and man-made infrastructures with affordable computational cost and mesh refinement. The Earth models with such complex model features can be designed by considering the reduction of insignificant dimensions of fractures/infrastructures and embracing their corresponding hierarchical material properties, and their computational meshes are generated by using Cubit. The finite element solutions are automatically distributed across nodes of parallel cluster for efficient computation of transient fluid flow or heat conduction in geologic media. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Beskardes, GungorDidem↗
DOE OSTI2024

Challenges and Opportunities for Basic Efficiency Measures in Low-Income Homes: A Southeast Alaska Case Study

Juneau, Alaska is the state's capital city and has a renewable energy goal to reach 80% renewable energy for the space heating and transportation sectors by 2045. In practical terms, this indicates a need to electrify both sectors, to take advantage of the inexpensive hydropower available from the local electric utility, Alaska Energy Light & Power. In an effort to complement and enable electrification, researchers examined the feasibility to deploy storm windows via a case study of installing storm windows in two local low-income homes. Newer models of storm windows provide an extra layer of insulation over existing windows, while preserving operability and views. They can also improve comfort and reduce noise. Researchers worked with the regional housing authority and a local builder to install the storm windows and replace inoperable windows in the two houses in 2021, in addition to conducting pre- and post-installation air leakage tests, energy monitoring, and occupant interviews. The team encountered several challenges, including a lack of egress windows, energy data from a wide variety of heating systems, extremely leaky homes, and installation issues such as windows that were not square. These results point to several barriers to widespread deployment of storm windows in the area, but also open the door to opportunities to design deployment programs that improve safety as well as efficiency.

cold climate↗
DOE OSTI2020

Research three-dimensional hierarchical graphite architectures for anodes for fast charging

With current lithium ion batteries optimized for performance under relatively low charge rate conditions, implementation of XFC has been hindered by drawbacks including Li plating, kinetic polarization, and heat dissipation. This project will utilize model-informed design of 3-D hierarchical electrodes to tune key XFCrelated variables like 1) bulk porosity/tortuosity 2) vertical pore diameter, spacing, and lattice 3) crystallographic orientation of graphite particles relative to exposed surfaces 4) interfacial chemistry of the graphite surfaces through "artificial sEr' formation using ALD 5) current collector surface roughness (aspect ratio, roughness factor, etc.). A key aspect of implementing novel electrodes is characterizing them in relevant settings. For this project, ultimately led out of University of Michigan by Neil Dasgupta, that includes both coin cell and 2+ Ah pouch cell testing, as well as comparison testing against baselines. Sandia National Labs will be conducting detailed cell characterization on iterative versions/improvements of the model-based hierarchical electrodes, as well as COTS cells for baseline comparisons. Key metrics include performance under fast charge conditions, as well as the absence or degree of lithium plating. Sandia will use their unique high precision cycling and rapid EIS capabilities to accurately characterize performance and any lithium plating during 6C charging and beyond, coupling electrochemical observations with cell teardown. Sandia will also design custom fixturing to cool cells during rapid charge, to decouple any kinetic effects brought about by cell heating and allow comparisons between different cells and charge rates. Using these techniques, Sandia will assess HOH electrodes from the University of Michigan, as well as aiding in iterative model and electrode design.

25 ENERGY STORAGE↗
DOE OSTI2021

Flow Assisted Evaporative Cooling for Electric Motor

This paper examines a novel concept of flow assisted latent heat driven two-phase evaporative cooling (EC) confined in-between slot liner and active-winding of electric motor. Wicking micro-structure enhanced PDMS liner axially sucks coolant in the form of thin film between the PDMS liner and active-winding and eventually enables thin film evaporation on the outer surface of the active-winding. Therefore, EC based thermal management eliminates contact resistance between the winding and slot-liner and, enhances the heat extraction from the winding without compromising the electro-magnetic performance. Two-way coupled electro-magnetic (EM) – computational fluid dynamics/heat transfer (CFD/HT) and EM - lumped parameter thermal network (LPTN) models have been developed to assess the electro-thermal performance of the EC under steady and transient conditions. Taking a case study of a 125 kW jacket cooled BMW i3 motor and dielectric coolant FC-84, EC is shown to be capable of handling a maximum steady state rms current density of 26 A/mm2 at a evaporative heat transfer coefficient of 5,000 W/m2.K, which is about 78.7% higher compared to the traditional jacket cooling (JC). In case of EC, a maximum steady state and transient rms current density of 30 A/mm2 (106.2% higher compared to the JC), and 40.8 A/mm2 have been realized by using high thermal conductivity epoxy (1.9 W/m.K) impregnation material. Thermal performance of the EC is also assessed and compared with JC over a dynamic drive cycle. Lastly, a motorette testing has been performed to demonstrate the applicability of the proposed EC method and to validate the developed modeling framework.

42 ENGINEERING↗
DOE OSTI2026

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗
DOE OSTI2026

Computationally inexpensive part-scale thermal history of additive friction-stir deposition

This study presents an analytical model for steady-state power generation and tool heat loss in additive friction-stir deposition (AFSD), developed to enable part-scale thermal simulation while remaining computationally inexpensive. The model predicts total generated power, yielding 3.7–4.7 kW across deposition temperature setpoints of 400–460 °C for the deposition of AA6061 with a Be-Cu tool. This corresponds to 90–95% of the reported spindle power. Tool heat loss is experimentally determined by calibrating a steady-state energy balance between the generated power, the substrate-deposition thermal gradient, and a temperature dependent tool heat loss term: q tool (T) = a + b (T - 400°C) with a = 2.7 x 10 6 Wm -2 and b = 9.5 x 10 3 Wm -2 K -1 . The calibration indicates that about 69% of the generated heat is conducted into the tool for this configuration, which is much higher than previously reported. The calibrated heat-source is implemented in finite element software (Adamantine) to simulate the transient thermal history of a 100 cm 3 representative build in 8 min on a standard desktop (at 0.635 mm build-height resolution). For the first three layers, the substrate temperatures between simulation and experiment are within 10% mean absolute percentage error. Sensitivity analysis indicates that uncertainties in average deposition temperature and deformation localization (stir-zone geometry, depth, and spatial dependance of strain-rate and flow stress) dominate model variance, motivating additional experimental verification.

Additive Friction-Stir Deposition↗
DOE OSTI2023

Carbon and Water Fluxes of the Boreal Evergreen Needleleaf Forest Biome Constrained by Assimilating Ecosystem Carbonyl Sulfide Flux Observations

Gross primary production (GPP) by boreal forests is highly sensitive to environmental changes. However, GPP simulated by land surface models (LSMs) remains highly uncertain due to the lack of direct photosynthesis observations at large scales. Carbonyl sulfide (COS) has emerged as a promising proxy to improve the representation of GPP in LSMs. Because COS is absorbed by vegetation following the same diffusion pathway as CO2 during photosynthesis and not emitted back to the atmosphere, incorporating a mechanistic representation of vegetation COS uptake in LSMs allows using COS observations to refine GPP representation. Here, we perform ecosystem COS flux and GPP data assimilations to constrain the COS- and GPP-related parameters in the ORCHIDEE LSM for boreal evergreen needleleaf forests (BorENF). Assimilating ecosystem COS fluxes at Hyytiälä forest increases the simulated net ecosystem COS uptake by 14%. This increase largely results from changes in the internal conductance to COS, highlighting the need to improve the representation of COS internal diffusion and consumption. Moreover, joint assimilation of ecosystem COS flux and GPP at Hyytiälä improves the simulated latent heat flux, contrary to the GPP-only data assimilation, which fails to do so. Finally, we scaled this assimilation framework up to the boreal region and find that the joint assimilation of COS at Hyytiälä and GPP fluxes at 10 BorENF sites increases the modeled vegetation COS uptake up to 18%, but not GPP. Therefore, this study encourages the use of COS flux observations to inform GPP and latent heat flux representations in LSMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Thermal performance of an absorption-assisted two-stage ejector air-to-water heat pump

In this work, a gas-fired air-to-water ejector heat pump (EHP) for domestic hot water production is proposed and analyzed thermodynamically. The COP of the EHP is improved by using binary fluids (HFE7500 and water), and the high backpressure in the primary ejector (first-stage) is actively controlled by either a second-stage ejector and/or a LiBr-H 2 O absorption subsystem. A thermodynamic model is established for analysis and performance evaluation. A case study shows the backpressure of the first-stage ejector is reduced by 53.5%, and the backpressure control using the second-stage ejector increases the system COP by 21.0% compared with that using LiBr-H 2 O absorption. Parametric studies are conducted to investigate the roles of the first-stage ejector, the second-stage ejector, and the LiBr-H 2 O absorption subsystem. The first-stage ejector plays a dominant role in the EHP system performance. The EHP’s heating COP can reach 2.0 with an entrainment ratio of 0.12 in the first-stage ejector. EHP performance can be improved with a high entrainment ratio and a low-temperature primary working fluid in the first-stage and second-stage ejectors, a low condensation temperature in the first-stage ejector, and a high generating temperature and a low concentration of the LiBr-H 2 O solution.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Mapping Thermoelectric Transport in a Multicomponent Alloy Space

Abstract Interest in high entropy alloy thermoelectric materials is predicated on achieving ultralow lattice thermal conductivity κ L through large compositional disorder. However, here it is shown that for a given mechanism, such as mass contrast phonon scattering, κ L will be minimized along the binary alloy with highest mass contrast, such that adding an intermediate mass atom to increase atomic disorder can increase thermal conductivity. Only when each component adds an independent scattering mechanism (such as adding strain fluctuation to an existing mass fluctuation) is there a benefit. In addition, both charge carriers and heat‐carrying phonons are known to experience scattering due to alloying effects, leading to a trade‐off in thermoelectric performance. Analytic transport models are applied, based on perturbation and effective medium theories, to predict how alloy scattering will affect the thermal and electronic transport across the full compositional range of several pseudo‐ternary and pseudo‐quaternary alloy systems. To do so, a multicomponent extension is demonstrated to both thermal and electronic binary alloy scattering models based on the virtual crystal approximation. Finally, it is shown that common functional forms used in computational thermodynamics can be applied to this problem to further generalize the scattering behavior that is modeled.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Comprehensive Review of Latent Space Dynamics Identification Algorithms for Intrusive and Non-Intrusive Reduced-Order-Modeling

Numerical solvers of partial differential equations (PDEs) have been widely employed for simulating physical systems. However, the computational cost remains a major bottleneck in various scientific and engineering applications, which has motivated the development of reduced-order models (ROMs). Recently, machine-learning-based ROMs have gained significant popularity and are promising for addressing some limitations of traditional ROM methods, especially for advection dominated systems. In this chapter, we focus on a particular framework known as Latent Space Dynamics Identification (LaSDI), which transforms the high-fidelity data, governed by a PDE, to simpler and low-dimensional latent-space data, governed by ordinary differential equations (ODEs). These ODEs can be learned and subsequently interpolated to make ROM predictions. Each building block of LaSDI can be easily modulated depending on the application, which makes the LaSDI framework highly flexible. In particular, we present strategies to enforce the laws of thermodynamics into LaSDI models (tLaSDI), enhance robustness in the presence of noise through the weak form (WLaSDI), select high-fidelity training data efficiently through active learning (gLaSDI, GPLaSDI), and quantify the ROM prediction uncertainty through Gaussian processes (GPLaSDI). We demonstrate the performance of different LaSDI approaches on Burgers equation, a non-linear heat conduction problem, and a plasma physics problem, showing that LaSDI algorithms can achieve relative errors of less than a few percent and up to thousands of times speed-ups.

Computational Engineering, Finance, and Science (c↗
DOE OSTI2023

Thermal Analysis of the SIRIUS3 Nuclear Propulsion Fuel Calibration Experiment

NASA is considering Nuclear Thermal Propulsion (NTP) for long range extraterrestrial missions; these engines eject hot hydrogen gas heated by a nuclear reactor for rocket thrust. To economize the use of hydrogen to the greatest extent possible, the NTP engines will be expected to, in a very short time (i.e., on the order of a minute or less), go from warm$(\sim 300 K)$ zero power conditions to full operational power, with a coolant outlet temperature on the order of 2700–3000 K \cite{en15176181}. These conditions will introduce significant thermomechanical stresses on the NTP fuel. The objective of the SIRIUS series of experiments is to examine the performance of candidate NTP fuel materials when subjected to temperature ramp rates that are prototypical of NTP system startup and operation. The SIRIUS experiments are a series of experiments that will be irradiated in the TREAT reactor and are subjected to power ramps and cycles that are prototypical for NTP operation. The experiments will be accomplished by executing a series of shaped transients on the SIRIUS specimens while collecting in situ specimen temperature data. These tests will determine whether operational startup ramps and peak temperatures will result in detrimental fuel performance phenomena (i.e. fuel deformation, fragmentation and cracking) To this end, INL has been evaluating a number of SIRIUS experiments, and these evaluations include thermal analysis of the SIRIUS experiments. Thermal analysis were performed for the calibration irradiation of the SIRIUS-3 experiment, and the focus of this memo is to document the results from the calibration irradiation thermal analysis. The thermal analysis reveals that radiation heat emission of the outer fuel elements and conduction to different metal components such as the molybdenum element tubes remove significant quantities of heat from the fuel element specimen, and future experiment designs need to consider these heat transfer mechanics. This paper begins with a brief experiment overview with a discussion of the fuel sample and experiment configuration. The experiment and model description section is followed by a set of results with a brief discussion, and finally the memo concludes with suggestion for future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Implementation and Evaluation of a Unified Turbulence Parameterization Throughout the Canopy and Roughness Sublayer in Noah‐MP Snow Simulations

Abstract The Noah‐MP land surface model (LSM) relies on the Monin‐Obukhov (M‐O) Similarity Theory (MOST) to calculate land‐atmosphere exchanges of water, energy, and momentum fluxes. However, MOST flux‐profile relationships neglect canopy‐induced turbulence in the roughness sublayer (RSL) and parameterize within‐canopy turbulence in an ad hoc manner. We implement a new physics scheme (M‐O‐RSL) into Noah‐MP that explicitly parameterizes turbulence in RSL. We compare Noah‐MP simulations employing the M‐O‐RSL scheme (M‐O‐RSL simulations) and the default M‐O scheme (M‐O simulations) against observations obtained from 647 Snow Telemetry (SNOTEL) stations and two AmeriFlux stations in the western United States. M‐O‐RSL simulations of snow water equivalent (SWE) outperform M‐O simulations over 64% and 69% of SNOTEL sites in terms of root‐mean‐square‐error (RMSE) and correlation, respectively. The largest improvements in skill for M‐O‐RSL occur over closed shrubland sites, and the largest degradations in skill occur over deciduous broadleaf forest sites. Differences between M‐O and M‐O‐RSL simulated snowpack are primarily attributable to differences in aerodynamic conductance for heat underneath the canopy top, which modulates sensible heat flux. Differences between M‐O and M‐O‐RSL within‐canopy and below‐canopy sensible heat fluxes affect the amount of heat transported into snowpack and hence change snowmelt when temperatures are close to or above the melting point. The surface energy budget analysis over two AmeriFlux stations shows that differences between M‐O and M‐O‐RSL simulations can be smaller than other model biases (e.g., surface albedo). We intend for the M‐O‐RSL physics scheme to improve performance and uncertainty estimates in weather and hydrological applications that rely on Noah‐MP.

Meteorology & Atmospheric Sciences↗
DOE OSTI2026

Large Eddy Simulations of the Atmospheric Boundary Layer Over Satellite‐Sensed Sea Ice Maps

Surface heterogeneity in the marginal ice zone (MIZ) causes multiscale secondary atmospheric circulations that are challenging to model or observe. The absence or inadequate representation of these circulations in ocean‐atmosphere exchange schemes in climate models is partially responsible for the underestimation of Arctic sea ice loss. Observationally, such circulations obfuscate the interpretation of polar atmospheric chemistry measurements, among others. To address this open challenge, large‐eddy simulations are conducted over real‐world satellite‐sensed sea ice maps with an ice‐sea temperature contrast, as well as over idealizations of these maps that alter the ice pattern but conserve its fraction, showing that the ice fraction of a surface is not sufficient to predict the surface heat flux. In a second suite of simulations, three other heterogeneity metrics (representing the surface fragmentation, patch edge tortuosity, and patch size variability) are introduced to complement the ice fraction in describing the surface. Simulations varying these parameters suggest that they also significantly impact surface‐air interactions. A roughness contrast is then added to the surface temperature contrast, indicating that the contribution of roughness changes to the resulting atmospheric circulations is less pronounced than that of thermal heterogeneity. Based on these results, we illustrate, using a multi‐linear regression on these map features, that MIZ surface parameterizations in Earth Systems Models can be improved if they account for these various characteristics of the sea‐ice patterns.

54 ENVIRONMENTAL SCIENCES↗