Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat Conduction”

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 541 records · Page 30

Reservoir Thermal Energy Storage Benchmarking (Rev. 3)

A benchmarking analysis of RTES research funded by GTO through the Beyond Batteries projects was conducted against the ESGC to see where they fit within the identified ESGC Use Cases. The projects were found to advance knowledge in multiple ESGC use cases, either directly or in some cases, indirectly as enabling technologies. This analysis is helpful to understand where RTES and associated research fits into the larger discussion around energy storage technologies. Also, a retrospective analysis of the Beyond Batteries projects was conducted to evaluate what the projects learned and how the results can be applied to advance the value of RTES. Major results of each of the studies are summarized in Table 2. Additionally, a comparative metrics analysis for RTES was completed to understand where RTES lies within the energy storage industry. Metrics for evaluation of RTES and its comparison to other storage technologies were selected and ranges of their values compiled. The selected metrics – LCOE (levelized cost of energy), capital costs, roundtrip efficiency, energy storage capacity, and storage time – were chosen based on data availability and have a particularly strong influence on the potential deployment of a storage technology. Charts which compare the metrics are presented in section 4.3 and show ranges for each of the 10 selected technologies. However, due to a lack of domestic operational facilities, values for RTES and for portions of the remaining technologies are based on theoretical modeling and studies of best-case scenarios. LCOE estimates for RTES fall within the lower reaches of Figure 15, but nevertheless amount to 2 – 5 times the ESGC Roadmap goal for LCOE, for example in the Facilitating and Evolving Grid Use Case. Capital costs for RTES sit on the higher end (Figure 16) but are expected to decrease as new projects are developed and the technology is refined. The theoretical roundtrip efficiency reported for RTES varies from mid to high percentages (Figure 17) with efficiencies upwards of 93% in modeled scenarios in the Portland Basin (Bershaw et al.,2020). RTES is also expected to have the largest energy storage capacities and longest storage times, likely matched only by lower efficiency hydrogen storage. To better assess the role that RTES could play in energy storage we examined it’s potential in the U.S. The potential depends on many factors. Recently, many researchers have started looking at deep sedimentary basins, depleted oil and gas fields, and basalt formations as potential targets for RTES development. The United States Geological Survey (USGS) has analyzed various cities and shown substantial RTES potential in the cooling sector (Pepin et al., 2021). By modeling RTES in low-quality groundwater (e.g., brackish), it is shown to be favorable across the U.S. with particular suitability in the Illinois Basin, Coastal Plains, and Basin and Range regions. Seasonal RTES operations have also been modeled in the Portland Basin by those at the USGS and Portland State University to simulate an RTES system supplying heating loads needed for the Oregon Health and Science University. Simulations suggest that high conductive heat loss in the initial years exists but tends to decrease with increasing time and development of the resource due to self-insulating nature of the basalts (Burns et al., 2020). Other national laboratory efforts are taking a close look at many of the technical issues involved with RTES (McLing et al., 2019, McLing et al., 2022). These include difficulties in understanding geochemical, hydrogeological, mechanical, and microbiological changes at such elevated temperatures and operational scenarios. Major gaps in research are identified and suggested for future work. With this increased focus to understand how to make RTES successful in the U.S., this technology could be a potential solution to many of the nation’s energy storage problems. For the energy independence of this country, the DOE should prioritize de-risking this technology by making future investments in pilot-scale demonstrations to attract potential investors.

15 GEOTHERMAL ENERGY↗

Heat transfer and transport property contrast effects on the compressible Rayleigh-Taylor instability

In extreme environments, the Rayleigh-Taylor instability (RTI) may occur under large variations in density and temperature and with fluid transport properties strongly dependent on temperature. Direct numerical simulations of the 3D fully compressible RTI are conducted, examining the idealized configuration of a hotter, less dense fluid pushing against a colder, denser fluid. Various temperature ratios and transport property configurations are explored to examine how heat conduction, large variations in transport properties, and sudden changes in transport properties can affect the evolution of the mixing layer. Nonuniform fluid expansion and contraction induced by heat transfer can significantly affect local density differences and overall instability growth, causing profile asymmetries about the initial interface for flow and mixing statistics. The departures from classical self-similar development of the instability along with misalignment between regions of mixing and regions of most intense turbulent activity caused by both heat transfer and transport property contrasts are examined. After sudden changes in fluid transport properties, which may occur as a result of rapid heating (e.g., in inertial confinement fusion), the flow quickly responds and begins to relax towards quasi-self-similar late-time evolution. For many dynamical quantities such as vorticity and dissipation, this late-time evolution resembles that of the configuration that already started with the final transport property magnitudes, suggesting that these quantities depend only on the transport properties and not on past flow history, provided that the density field distributions for the flows remain similar. On the other hand, the mixing evolution after the transport property change is unique, implying that both property magnitudes and previous history are impactful on the mixing. In conclusion, these simulations demonstrate how various temperature-related effects are extremely important to consider in compressible RTI flows with large temperature variations.

42 ENGINEERING↗

Interaction of magnetic fields with spinons in a fractionalized state

The 4d-electron trimer lattice Ba₄Nb₁₋ₓRu₃₊ₓO₁₂ exhibits either a quantum spin liquid (QSL) or a heavy-fermion strange metal (HFSM) phase, depending on Nb content. In the QSL state, itinerant spinons act as effective heat carriers, enhancing thermal conductivity. Strikingly, applying a magnetic field up to 14 T causes an abrupt, up-to-5000% increase in heat capacity below 150 mK, disrupting the linear temperature dependence typical of both phases. Meanwhile, AC susceptibility and electrical resistivity remain nearly unchanged, while thermal conductivity drops by up to 40% below 4 K. These results suggest spinons, despite being charge-neutral, are highly sensitive to magnetic fields at low temperatures. We propose that the magnetic field could induce Anderson localization of spinons, creating emergent non-magnetic two-level systems responsible for the Schottky-like anomaly in heat capacity. These findings point to a previously unexplored regime of spinon dynamics, potentially governed by field-induced localization and distinct from conventional magnetic or transport signatures.

Condensed-matter physics↗

Structural and thermal properties and the origin of the ultralow thermal conductivity in the defect stannite CuIn 2 Se 4

Phase-pure CuIn 2 Se 4 , a ternary metal chalcogenide that forms in a disordered stannite crystal structure, was synthesized to investigate the structure thermal property relationships as well as reveal the origin of the ultralow thermal conductivity this material possesses over a large temperature range. Modeling of the temperature-dependent heat capacity and thermal conductivity revealed distinctive thermal properties and large lattice anharmonicity. Electron localization function calculations highlight the asymmetric bonding inherent to CuIn 2 Se 4 , which together with lattice anharmonicity directly impacts the thermal properties. Our findings reveal the specific atomic arrangement and bonding governing the thermal properties of this ternary metal chalcogenide. Our findings underscore the specific atomic arrangement and bonding governing the thermal properties in this ternary chalcogenide. This study advances the fundamental understanding of stannites, and our findings can be applied to these and other multinary metal chalcogenides of interest for applications where low thermal conductivity is desirable.

36 MATERIALS SCIENCE↗

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗

Thermal conductivity of salvaged fused cast alumina used in the glass industry

Fused cast alumina (FCA) has been and continues to be used as a refractory material in energy intensive industries such as glass melting and chemical processing. In-service degradation due to high temperature exposure in harsh environment affects the designed furnace thermal profiles and energy consumption. Phase transformation may occur at the refractory hot face during glass melting altering the properties. Three FCA blocks recovered from industrial furnaces were investigated in this study. The as-received FCA consists primarily of a mixture of alpha (α) and beta (β) alumina that has a thermal conductivity value of 5–6 W/mK. The Hot Disk method was used to obtain thermal conductivity directly on the refractory blocks. At the hot face, a transformation from β to α alumina occurred and was confirmed by an X-ray diffraction study. Thermal conductivity measurements as a function of position also showed a clear transition from β to α alumina at both ends of a complete block with no voids. Thermal conductivity of the α alumina tripled compared to β alumina. This study provides important information of heat transfer and thermal conductivity evolution to refractory manufacturers and users.

36 MATERIALS SCIENCE↗

Experimental analysis of conductive drying of paper and board

The primary heat transfer source in the conventional multi-cylinder drying of paper and board is conduction. Conduction is facilitated by high-tension contact with steam-heated cylinders, while convective drying, the main mass transfer source, operates as heated air flows over the paper web in the pockets. This conduction process occurs at elevated temperatures and contact pressures to ensure effective contact between the wet paper web and the heated cylinders. The contact pressures and temperatures of the hot surface significantly influence the conductive drying characteristics. This research paper presents an experimental study that involves designing a simple lab-scale setup with in-situ and continuous sensing capabilities for various commercially available grades of paper and board. Embedded thermocouples measure the temperatures of the heated platen and the sheet, allowing the collection of flux data to determine heat transfer characteristics. Here, the goal of this study is to ascertain the instantaneous contact coefficients for different basis weights as a function of moisture content. The acquired data will provide valuable insights and information towards process development, design and simulation of the paper drying process.

42 ENGINEERING↗

Thermal conductivity and its relation to atomic structure for symmetrical tilt grain boundaries in silicon

We perform a systematic study of thermal resistance and conductance of tilt grain boundaries (GBs) in Si using classical molecular dynamics. The GBs studied are naturally divided into three groups according to the structural units forming the GB core. In this work we find that, within each group, the GB thermal conductivity strongly correlates with the excess GB energy. All three groups predict nearly the same GB conductivity extrapolated to the high-energy limit. This limiting value is close to the thermal conductivity of amorphous Si, suggesting similar heat transport mechanisms. While the lattice thermal conductivity decreases with temperature, the GB conductivity slightly increases. However, at high temperatures it turns over and starts decreasing if the GB structure undergoes a premelting transformation. Analysis of vibrational spectra of GBs resolved along different directions sheds light on the mechanisms of their thermal resistance. The existence of alternating tensile and compressive atomic environments in the GB core gives rise to localized vibrational modes, frequency gaps creating acoustic mismatch with lattice phonons, and anharmonic vibrations of loosely bound atoms residing in open atomic environments.

36 MATERIALS SCIENCE↗

Ductless Heat Pump Energy Monitoring in a Cold Climate

In 2018, the City and Borough of Juneau, Alaska set the goal of reaching 80% renewable energy for heating and transportation by 2045. An important factor to help reach this goal is finding energy efficiency home heating technologies such as ductless air source heat pumps (DHPs), which rely on the local hydropower electric grid to provide space heating for an entire home or a main living area. To help residents reduce energy costs and move toward the community renewable energy goal, a team of local, state, and national partners joined to run Thermalize Juneau 2021, an energy campaign that provided education, a simplified installation process, and a bulk purchase discount. Overall, the campaign enrolled 165 participants, resulting in 75 DHP installations and over 30 energy efficiency improvements. To evaluate the reliability and energy use of the DHPs installed during the Thermalize Juneau campaign, researchers installed energy monitors on heat pumps of 14 participants and conducted interviews with homeowners. Data from two heating seasons (2021-2022 and 2022-2023) has allowed researchers to examine the accuracy of pre-campaign energy savings predictions, establish energy savings for a small group of households, and provide insights on the energy use and occupant satisfaction of DHPs which will help to inform future heat pump deployment efforts in the area.

air source heat pump↗

Ductless Heat Pump Energy Monitoring in a Cold Climate

In 2018, the City and Borough of Juneau, Alaska, set the goal of reaching 80% renewable energy for heating and transportation by 2045. An important factor to help reach this goal is energy-efficient home heating technologies such as ductless air- source heat pumps (DHPs), which rely on the local hydropower electric grid to provide space heating for an entire home or a main living area. To help residents reduce energy costs and move toward the community renewable energy goal, a team of local, state, and national partners joined together to run Thermalize Juneau 2021, an energy campaign that provided education, a simplified installation process, and a bulk purchase discount. Overall, the campaign enrolled 165 participants, resulting in 75 DHP installations and over 30 other energy efficiency improvements. To evaluate the reliability and energy use of the DHPs installed during the Thermalize Juneau campaign, researchers installed energy monitors on heat pumps of 14 participants and conducted interviews with homeowners. Data from two heating seasons (2021-2022 and 2022-2023) allowed researchers to examine the accuracy of pre-campaign energy savings predictions, establish energy savings for a small group of households, and provide insights on the energy use and occupant satisfaction of DHPs, which will help inform future heat pump deployment efforts in the area.

air source heat pump↗

Processing Compressed Expanded Natural Graphite for Phase Change Material Composites

Phase change materials (PCMs) are used in various thermal energy storage applications but are limited by their low thermal conductivity. One method to increase conductivity involves impregnating organic PCMs into highly porous conductive matrix materials. Of these materials, compressed expanded natural graphite (CENG) matrices have received the most attention. Despite this attention, the effect that CENG processing has on PCM saturation and overall matrix thermal conductivity has not been fully investigated. Therefore, the effect of the heat treatment process used to expand intercalated graphite flakes is evaluated here. Higher heat treatment temperatures yielded higher saturation rates and overall saturation at similar matrix porosities. For example, increasing temperature from 300 °C to 700 °C resulted in approximately 60%-70% increase in pore saturation after 100 minutes of soaking. The exposure time to heat treatment had less of an effect on PCM saturation. The exposure time had negligible effect above 30 min and above 500 °C heating temperatures. However, because the expanded graphite was found to oxidize around 700 °C, the use of longer exposure time in manufacturing applications can be beneficial if a shortened impregnation time is needed. Heat treatment conditions did not impact thermal conductivity. The composite latent heat of fusion was also reduced approximately proportionally to the PCM mass fraction. A local maximum in axial thermal conductivity was observed at around 83% porosity, which is similar to previous studies. The observed conductivity at this maximum was a factor of 81 times greater than the conductivity of the PCM.

36 MATERIALS SCIENCE↗

Heat Transfer in Nuclear Waste Glasses: Measurements and Modeling of Thermal Radiation Properties

We measured and modeled near infrared extinction of nuclear waste glasses from 300 °C to 1150 °C to enable predictive radiation heat transfer and thermal conductivity estimates. A composition and redox informed model resolved contributions from key chromophores (Fe+2-O-Fe+3, V+4, free and bonded ?Si-OH groups) and, when present, spinel particles that can cause strong scattering. The model reproduced measured absorption from room temperature up to 1150 °C, with minor discrepancies near 1 µm (likely due to possible trace impurities) and 2.5 µm (linked to uncertainty in hydroxy groups). Spectra showed silicate melts were semitransparent mainly in the 0.5–4.0 µm window, responsible for radiation thermal conductivity that generally increases with increasing temperature. We quantified the dependence of effective thermal conductivity on dissolved water and provided distributions across >100 LAW/HLW/DFHLW melts at 1150 °C, supporting improved melter heat transfer modeling.

Ferkl, Pavel↗

Coupled Monte Carlo Transport and Conjugate Heat Transfer for Wire-Wrapped Bundles Within the MOOSE Framework

Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics (CFD) to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. We first provide a brief introduction to Cardinal's software design, data mapping, and coupling strategy to highlight our approach to overcoming common challenges in high-fidelity multiphysics simulations. Here we then present two Cardinal simulations for hexagonal pin bundles. The first is a validation of Cardinal's conjugate heat transfer coupling of NekRS's Reynolds-Averaged Navier Stokes model with MOOSE's heat conduction physics for a bare seven-pin Freon-12 bundle flow experiment. Predictions for pin surface temperatures under three different heating modes agree reasonably well with experimental data and similar CFD modeling from the literature. The second simulation is a multiphysics coupling of OpenMC, NekRS, and BISON for a reduced-scale, seven-pin wire-wrapped version of an Advanced Burner Reactor bundle. Wire wraps are approximated using a momentum source model, and coupled predictions are provided for velocity, temperature, and power distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Thermal-hydraulic Performance of Emerging Low GWP Refrigerant Mixture under Flow Boiling in Brazed Plate Heat Exchangers

The hydrofluorocarbon (HFC) refrigerants used in the current refrigeration systems are facing a phase-down due to their higher greenhouse effect resulting in global warming, and thus HVAC&R industry has undergone a transition to low Global Warming Potential (GWP) refrigerants. Refrigerant mixtures are attractive alternatives since their composition can be tailored to comply with environmental regulations while preserving favorable thermophysical properties. However, the new low-GWP zeotropic mixture refrigerants have two or more components with different saturation temperatures at the same pressure level, known as temperature glide, which can cause the degradation of the overall heat transfer performance. The brazed plate heat exchangers (BPHX) provide excellent heat transfer performance due to a compact design and are used in several air-conditioning and refrigeration applications. In this study, flow boiling heat transfer and the associated pressure drop of the refrigerant mixture in a vertical BPHX were experimentally investigated. The single-phase water-to-water experiments were conducted in the tested heat exchanger with a counter-flow configuration. The flow boiling experiments charged with R-134a and R-454C were then performed in a pumped refrigerant loop to evaluate its thermal-hydraulic performance. Furthermore, parametric studies of various heat fluxes, mass fluxes, vapor qualities, and saturation temperatures were also conducted.

Yang, Cheng-Min↗

The value of integrating a geothermal district heating system into a microgrid

As electrical grids increasingly rely on variable renewable energy, maintaining reliability and cost efficiency becomes more complex. To address these challenges, this study analyzed the integration of geothermal district heating as a grid-responsive thermal resource within a microgrid in Tuttle, Oklahoma. Building energy modeling using EnergyPlus estimated annual district heating demand at 2.9 GWh, with a peak load of 2.8 MW th . Techno-economic analyses were conducted to meet the heating demand under three geothermal scenarios, varying by production depth, flow rate, and thermal output, each supplemented by natural gas peaking boilers. In parallel, equivalent electrical load profiles were developed using typical coefficients of performance (COPs) for air-source heat pumps and electric boilers to establish an electrified baseline scenario. A complete end-use electrical load profile was also developed for the microgrid using Cambium dataset. The modeling results demonstrated reliable and economic operation of the geothermal systems over 30 years, with COPs ranging from 2.6 to 8.9 and the lowest levelized heating cost at $\$$54.6/MWh. Geothermal integration reduced electricity consumption by up to 94.7 % compared to the non-geothermal base case, yielding annual energy savings of up to $\$$803 k. Avoided grid costs ranged from $\$$65 k–$\$$147 k per year, with individual events avoiding up to $\$$4,863 per hour. Grid-responsive operation further reduced wholesale energy costs by 53–56 %. These findings demonstrate geothermal heating, traditionally treated as a non-grid-responsive thermal resource, can be reconfigured to support dynamic grid services, offering a scalable pathway to enhance reliability and reduce costs in renewable-rich microgrids and district heating networks.

15 GEOTHERMAL ENERGY↗

Risk Importance Ranking of Fire Data Parameters to Enhance Fire PRA Model Realism

Fire is historically and analytically a significant contributor to nuclear power plant risk. The level of fire risk and the methods, tools and data for modeling this risk is highly debated by experts. One area of debate is the input data used in fire modeling and how to deal with this data’s high uncertainty. This report outlines initial work performed for determining the key parameters causing this uncertainty and how it propagates into nuclear power plant models. This research paves the way for the development of methods to reduce fire data uncertainty used in modeling. The Nuclear Regulatory Commission has mandated that nuclear power plants perform fire risk modeling. However, there are several issues with the current risk modeling implementation that affect the results. Approved modeling methods can be overly conservative and often do not match plant experience. Also, the data used in the modeling can have high uncertainties and is influenced by expert judgement. To evaluate input data uncertainty, researchers performed an initial review of several fire experiments done at Sandia National Laboratories. Uncertainties for fire data can come from many sources, such as experiment design constraints, environmental conditions, or other plant-specific aspects. There are many different significant and insignificant parameters driving the uncertainty. Additionally, the uncertainty of the different input data used in the fire modeling could have a significant or insignificant effect on the entire plant risk. A four-step methodology was developed to perform Integrated Probabilistic Risk Assessment Importance Ranking. A demonstration case using these steps was set up and three of the four steps were completed in fiscal year (FY) 2019 and the fourth step done FY 2020. These steps are: 1. The qualitative analysis of potential sources was conducted with the following items identified for the demonstration. • Maximum heat release rate • Time to maximum heat release rate • Duration of max heat release rate • Time to decay • Thermal conductivity of concrete • Specific heat of concrete • Density of concrete • Cable jacket thickness 2. A quantitative characterization of dominant sources of uncertainty was performed. A list of distributions and determined values of the dominant sources is shown in Appendix A. 3. A quantitative screening of the potential sources of uncertainty using Morris Elementary Effects Analysis was completed. An experimental model using the physics-based fire modeling tool Fire Dynamics Simulator was developed and coupled with the Risk Analysis Virtual Environment. The Morris analysis identified at least two parameters that can be eliminated as significant contributors (specific heat of concrete and cable jacket thickness). 4. Global importance measure (Global IM) analysis to generate a comprehensive ranking based on their influence on the plant risk. In this research, a moment-independent Global IM is used since it can address (a) uncertainty in the input parameters of the fire model, (b) uncertainty in the risk outputs, and (c) non-linearity and interactions among input parameters in the fire model, more accurately than the correlation-based and variance-based global methods. The observations from the research showed that, depending on the initial and boundary conditions of the fire scenarios, fire-induced damage could have a very small probability and could be dominated by the tail of the uncertainty distribution; hence, the accuracy of the correlation-based and variance-based methods is questionable. The moment-independent Global IM analysis in this research provides a better understanding of how experimental uncertainty data affects industry’s plant models and where improvements in that data will have the largest benefit for improving fire modeling accuracy in causing core damage. Among the five unscreened parameters obtained from the Morris EE analysis, the Global IM analysis results for the case study indicated that max heat release rate and fire location are the most important parameters. The report also outlines benefits of using a unified computational platform that integrates the underlying simulations (e.g., a fire progression model), quantitative screening (using the Morris EE method), and the Global IM analysis. A unified platform can (i) facilitate the ranking of input parameters considering multiple key fire scenarios simultaneously, rather than considering one scenario at a time, (ii) contribute to more explicit and accurate treatment of dependencies at multiple levels of Fire PRA, (iii) facilitate the sampling-based uncertainty quantification for Fire PRA, and (iv) help generating both “industry-wide” and "plant-specific" ranking of uncertainty sources in Fire PRA. Future research should be done to include additional parameters such as detection/suppression or cable fire spread. Adding a PRA software such as SAPHIRE to the RAVEN platform would help with plant model integration and improve treatment of fire-induced dependency. The I-PRA risk importance ranking methodology offered in this report can provide valuable information for efficiently (a) enhancing the realism of Fire PRA for existing plants and (b) supporting the development of Dynamic Fire PRA for advanced reactors and new plants.

97 MATHEMATICS AND COMPUTING↗