SEARCH · Search NASA

Results for “two phase heat transfer”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Paradigm for approaching the forbidden spontaneous phase transition in the one-dimensional Ising model at a fixed finite temperature

The Ising model describes collective behaviors such as phase transitions and critical phenomena in various physical, biological, economical, and social systems. It is well known that spontaneous phase transition at finite temperature does not exist in the Ising model with short-range interactions in one dimension. Yet, little is known about whether this forbidden phase transition can be approached arbitrarily closely—at fixed finite temperature. Here I use symmetry analysis of the transfer matrix to reveal the existence of spontaneous ultranarrow phase crossover (UNPC) at finite temperature in one class of one-dimensional Ising models on decorated two-leg ladders, in which the crossover temperature T 0 is determined solely by on-rung interactions and decorations, while the crossover width 2 δ T is independently, exponentially reduced ( δ T = 0 means a genuine phase transition) by on-leg interactions and decorations. These findings establish a simple ideal paradigm for realizing an infinite number of one-dimensional Ising systems with spontaneous UNPC at desirable T 0 , which would be characterized in routine laboratory measurements as a genuine first-order phase transition with large latent heat thanks to the ultranarrow δ T (say less than one nanokelvin), paving a way to push the limit in our understanding of phase transitions and the dynamical actions of frustration arbitrarily close to the forbidden regime. Published by the American Physical Society 2024

1-dimensional systems↗

DOE OSTI2020

COMBLE-ISLAS water isotopes in precipitation (COMBLEISLASISO)

During the measurement period, precipitation was sampled on daily to sub-daily basis with a sampling kit consiting of a clear plastic box (60x40x40 cm) mounted to a ground structure. Daily sampling lasted from 22 Feb to 24 Mar 2020. At a fixed time of day (10 UTC), samples were collected and the sampling box cleaned for the next sample. Additionally, during intense observations periods (IOPs) announced by ISLAS, ARMS scientists increased sampling frequency to sub-daily sampling. During the campaign period, 2 IOPs were conducted, with in total 12 samples taken. IOP1 lasted from 27 Feb 10 to 02 Mar 2020, and IOP2 lasted from 12 Mar to 14 Mar 2020. Solid precipitation in the box was then melted, transferred to an 8 ml glass vial, sealed and shipped to FARLAB (University of Bergen, Norway) for stable water isotope analysis. Measured stable isotope composition of precipitation samples will be related to conditions at other sampling locations, and at the evaporation site to inform about conservation of water isotope quantities. In the period from 22 Feb to 24 Mar, 23 daily samples were collected at 10 UTC for stable water isotope analysis. During the entire sampling period, a high variability of precipitation phase was encountered. Most samples were recovered in the form of graupel/snow (12 samples), as rain (7 samples), and the remaining samples fell as graupel, snow, rain or frozen rain. Precipitation amounts varied between <1 mm (26 Feb) and 140 mm (25 Feb) within the sampling box. Furthermore, higher-resolution sampling was conducted during 2 intense observation periods (IOPs). During IOP1, lasting from 27 Feb 10 UTC to 02 Mar 10 UTC, 7 samples were collected, all of them as solid phase precipitation. IOP2 lasted from 12 Mar 10 UTC to 14 Mar 09 UTC, and 5 samples (all as graupel/snow) were collected. Samples were processed according to FARLAB standard measurement procedures. In short, samples were transferred to 1.5 ml glass vials with rubber/PTFE septa (part #548-0907, VWR, USA). An autosampler (A0325, Picarro Inc) transferred ca. 2µl per injection into a high-precision vapourizer (A0211, Picarro Inc, USA) heated to 110°C. After blending with dry N2 (< 5 ppm H2O) the gas mixture was directed into the measurement cavity of a Cavity-Ring Down Spectrometer (L2140-i, Picarro Inc) for about 7 min with a typical water concentration of 20 000 ppm. Memory effects were reduced by two times measuring a vapour mixture at a mixing ratio of 50 000 ppm, obtained from 2 injections of 2 µl for 5 min at the beginning of each new sample vial. Thereafter, another 6 injections of 2 µl per sample were measured individually as described above, and averages of the last 5 injections were used for further processing. Three standards were measured at the beginning and end of each batch, including a drift standard DI2 (δD: -50.72±0.73 permil, δ18O: -7.63±0.10 permil), and for calibration the laboratory standards GLW (δD: -307.79±0.75 permil, δ18O: -40.02±0.07 permil) and EVAP2 (δD: 9.52±0.65 permil, δ18O: 1.81±0.13 permil). A detailed calibration report is included with the final uploaded data set.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

A field–particle correlation analysis of a perpendicular magnetized collisionless shock

Using the field–particle correlation technique, we examine the particle energization in a three-dimensional (one spatial dimension and two velocity dimensions; 1D-2V) continuum Vlasov–Maxwell simulation of a perpendicular magnetized collisionless shock. The combination of the field–particle correlation technique with the high-fidelity representation of the particle distribution function provided by a direct discretization of the Vlasov equation allows us to ascertain the details of the exchange of energy between the electromagnetic fields and the particles in phase space. We identify the velocity-space signatures of shock-drift acceleration of the ions and adiabatic heating of the electrons arising from the perpendicular collisionless shock by constructing a simplified model with the minimum ingredients necessary to produce the observed energization signatures in the self-consistent Vlasov–Maxwell simulation. We are thus able to completely characterize the energy transfer in the perpendicular collisionless shock considered here and provide predictions for the application of the field–particle correlation technique to spacecraft measurements of collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Plasma production and ion heating systems for the Material Plasma Exposure eXperiment

Plasma production and ion heating in the Material Plasma Exposure eXperiment (MPEX), whose design is nearing completion, is accomplished using continuous wave RF power with average power density up to 1.6 MW/m2 at the plasma interface. Plasma is produced using helicon waves coupled through a single helical antenna at 13.56 MHz with power supplied by three 100 kW fixed-frequency RF generators feeding a power combiner network followed by the matching network and launcher. The helicon source can utilize various gasses including hydrogen, deuterium, and helium, with a magnetic field strength in the source region up to 0.2 T, and maximum device |B| of 2.5 T. Power is coupled to ions via ion cyclotron heating at the fundamental resonance using a pair of phased helical antennas operating in the frequency range 4-9 MHz, that launch waves towards the resonance from the antenna region where ω > ωci. ICH power is supplied by a single 500 kW tunable RF transmitter through a 90° power splitter and matching/decoupling network. In the case of both the helicon and ICH systems the antennas are located external to the vacuum, with power transferred through novel water-cooled coaxial vacuum windows consisting of fused quartz outer cylinders and silicon nitride inner cylinders with forced convection water cooling between them. The antenna enclosures are pressurized with dry air to 3 bar absolute for voltage standoff.Several 3-D COMSOL models have been created to simulate the two systems. A model of the helicon region utilizing a cold plasma dielectric tensor with accurate magnetic field and realistic plasma density profiles has been used to calculate the plasma loading/complex antenna input impedance at the launcher feed for various ne and |B| values, for the purpose of estimating power handling. It also incorporates the geometry of all launcher structures relevant to this determination. A still more detailed model of the launcher together with a lossy dielectric plasma surrogate has been used to determine RF electric field values and power losses in the device components.Similar models have been produced to predict the performance and power handling of the ICH launcher. For this device the impedance matrix of the two-element antenna array is calculated using a warm plasma model, necessary to properly determine the wave propagation and absorption.

Goulding, Richard↗

DOE OSTI2021

Critical roles of pores and moisture in sustainable nanocellulose-based super-thermal insulators

In a recent report in Matter, Bergström and co-workers describe a novel thermal transport behavior that subverts this assumption of porous materials. In their work, the authors describe a nanocellulose-based foam that features ultrahigh porosity (> 99.6%) and aligned µm-scale pores (10–100 µm) whose radial thermal conductivity (i.e., perpendicular to the pore alignment) is close to that of free air in the dry state. Furthermore, the radial thermal conductivity of the cellulose-based foam can be reduced to ~14 mW/(m·K) when the relative humidity is ~35%. Super-thermal insulators, which demonstrate a thermal conductivity below that of stationary air (~ 25 mW/(m·K), 20 ºC, 1.0 atm), are needed to minimize heat loss in various applications (e.g., buildings, thermal energy storage tanks, cold chain packaging, etc.) to mitigate the energy crisis and reduce carbon emissions. Introducing pores into a material is a facile and effective way to achieve low thermal conductivity as pores can suppress thermal transport through solids by reducing the cross-sectional area and increasing the tortuosity of the heat transfer pathway. In porous structures there are two kinds of pores: open and closed. While increased porosity can reduce heat conduction through solids with an open porous structure, the improved gas conduction creates a competing effect that simultaneously elevates the heat transfer. Therefore, the thermal conductivity of a material with µm-scale open pores is usually larger than that of stationary air. Reducing the pore size to less than the mean free path of air (~ 70 nm, 20 ºC, 1.0 atm) can effectively reduce gas conduction, enabling the material to achieve a thermal conductivity below that of stationary air. However, high cost of the nanosized raw materials and time-consuming fabrication processes limit the large-scale applications of nanoporous thermal insulators. Meanwhile, closed pores can block heat transport only through the continuous gas phase. As a result, the thermal conductivity of a closed porous structure can theoretically be much smaller than that of stationary air if the solid conduction can also be suppressed by (1) lowering the solid content of the material, (2) reducing the thermal conductivity of the building blocks of the material, and/or (3) increasing interfacial thermal resistance between neighboring building blocks. Most processes used to generate porous structures (e.g., supercritical drying, freeze drying) involve a solvent that escapes the material. Therefore, it is difficult to create pores and isolate them simultaneously, preventing the fabrication of super-thermal insulators with closed pores.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Development & Experimental Validation of a Generalized Resistance-Capacitance Model for Numerical Simulation of Phase-Change Material Embedded Heat Exchangers

Latent heat thermal energy storage (LHTES) using phase change material (PCM) has attracted increased attention as a viable solution for overcoming the mismatch between energy supply and demand for renewable energy-based systems. PCM-embedded heat exchangers (PCM-HX) have the potential to significantly improve thermal performance due to high storage capacity and low temperature variation during the phase change process. Most models for simulating LHTES heat transfer use Computational Fluid Dynamics (CFD) simulations, which have high computational costs resulting from considering the complex and time-dependent physics relevant to PCM-HXs. In this paper, a Generalized Resistance Capacitance-based Model (GRCM) was developed to predict the thermal performance of arbitrary PCM-HXs in a computationally efficient manner without compromising modeling accuracy. The GRCM is exercised for three case studies: (i) verification for a single-slabbed finned PCM-HX, (ii) verification and validation for a copper foam/paraffin composite PCM-HX, and (iii) validation for a straight tube annular finned PCM-HX. The copper foam PCM-HX uses an electric heater at the top of HX, while the other two configurations utilize water as heat transfer fluid. For the single-slabbed finned PCM-HX melting case, the mean deviation in average PCM temperature predicted by the GRCM compared to the CFD model was between 0.56 – 0.73 K, with maximum temperature deviation of 2.68 K. For the HTF outlet temperature, the validation results showed that GRCM prediction matches very well with experimental data, with mean temperature deviation of 0.24 K during melting case, while for solidification case was 0.34 K. These results showcase the GRCM’s capability for accurately reproducing the thermal characteristics of PCM-HXs with considerably lower computational effort.

42 ENGINEERING↗

DOE OSTI2022

Numerical Modeling of Atmospheric Rime Ice Accretion on an Airfoil Using an Eulerian Approach

This paper presents an approach to numerically simulate the inherently unsteady rime ice accretion problem on a two-dimensional airfoil and elucidate the associated variations under different icing conditions. The airflow field and the water impingement on the airfoil are obtained based on an Eulerian two-phase model. A dynamic mesh strategy is employed to unsteadily account for the changes in the ice profile and its impact on the air and droplet flow by continuously reconstructing the computational grid at each time-step through smoothing and layering mechanisms. All main icing modules including the airflow field, droplet trajectory, icing thickness profile, and mesh management are fully coupled within the same computational framework without resorting to any external tools. Classical icing theory is employed to model the rime ice roughness, and it is assumed that the ice accretes in a direction normal to the airfoil surface. The governing Reynolds-averaged Navier–Stokes (RANS) conservation equations along with the energy and continuity equations are solved to produce the velocity and temperature fields. A convective film heat transfer coefficient is computed based on the surface heat flux and a recovery temperature which takes into account the dissipative heat release in the boundary layer in the vicinity of the airfoil surface. With the implemented strategy and calculating the convective heat transfer coefficient, the water film thickness is also calculated along with the ice shape. Furthermore, the model is validated by comparing the local collection efficiency distribution and ice shape with experimental data, and the results show that the implemented approach provides acceptable predictions of ice accretion profiles and rates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

Falling-film thermosyphons: Application to water harvesting from humid gas streams

Water harvesting from humid gas streams is central to fresh-water production, desalination, and particulate removal from combustion gas streams. Two-phase thermosyphons are well suited for these applications due to their very low thermal resistance, but have several operating limits at the high heat flows required for such applications. This work introduces the falling-film thermosyphon for large-scale condensation applications. Working fluid is pumped to the top of the evaporator to provide the evaporating liquid film on the inner tube wall, rather than by vapor condensation. The flooding, dry-out and pool-boiling limits are eliminated, resulting in significantly higher maximum heat fluxes for the same physical evaporator size. Furthermore, the condenser no longer needs to be located vertically above the evaporator, which allows for a standard steam condenser to be used. A model for the condensation process in humid air was developed that estimates the fluid and heat transfer and is confirmed experimentally. Further benefits include the use of a high-thermal-conductivity polymer material for the evaporator section to minimize corrosion, and the ability to impose a temperature boundary condition on the evaporator, which is made possible due to the elimination of a liquid pool.

42 ENGINEERING↗

DOE OSTI2020

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗

DOE OSTI2022

FY22 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY22. As an extension of the effort from FY21, this year’s work focuses primarily on the two-phase modeling of the NSTF using RELAP5-3D, particularly with the inclusion of the cavity model. The results from simulations were used to compare against experimental data for benchmarking purposes of the RELAP5 deck. Additionally, RELAP5 was used as a predictive tool to guide planned test operations and identify expected system behaviors. In the first part of this report, details are provided of the cavity omitted model where heat flux is applied directly as a boundary condition to the risers. The general trend predicted by the RELAP5 model matches that from the experimental data when a single-phase natural circulation flow is first established, followed by an oscillatory two-phase period and finally a stable two-phase flow. However, the onset of oscillations is predicted early by the model due to the smaller thermal mixing region in the tank. However, by expanding the simulated thermal mixing region in the tank, the onset of oscillations predicted by the model is able to match that from the experiment. These oscillations where studied in depth and are deduced to be flashing-induced instability. The model was then modified to simulate an accident scenario case where a representative heat load based on the full-scale Framatome’s 625 MW t SC-HTGR was applied directly to the riser channels. The simulated initial and boundary conditions were identical to those performed experimentally, facilitating direct comparisons between the predicted and experimental results. It was determined that the results showed some discrepancies remain, likely due to the overprediction of vapor generation rate by the computer model. In the second part of this report, the cavity model is re-introduced where it is observed that the RELAP5 prediction is now able to capture the major trends of the observed flow commonly observed during two-phase conditions. However, the onset of oscillations is once again predicted early by the model, possibly caused by the underprediction of heat loss from the heater and cavity. This is likely due to the omission of support structures in the cavity that can act as additional pathways for heat to escape to the environment. To overcome the underprediction of heat loss, part of the insulation surrounding the cavity side panels and the back of the heaters are removed to allow heat to escape directly to the environment, which then improves the RELAP5 prediction. Parametric studies are also performed to investigate the effects of heater power, tank inventory level, and tank gas space pressure on flow behaviors, also in direct comparison to conditions tested experimentally. User option-61 in the RELAP5-3D input deck, which changes the heat transfer coefficient correlations used for calculating the vapor generation, is also investigated where it is found that by enabling the option, the overall duration of oscillations is increased and matches that from the experiment better. The RELAP5 model is further benchmarked with a header inlet- throttling case where it is observed that the prediction from the model fails to capture some major features observed in the experiment. By using a modified loss coefficient curve for the valve, the accuracy of the prediction is improved where most of the major features observed in the experiment are predicted by the model. Lastly, the model is benchmarked with an inventory depletion scenario where it is observed that despite the modeling limitation of RELAP5, the prediction shows good agreement with the experimental data where major trends and features are captured by the model. Future work will see continued development of the current RELAP5-3D input deck of the NSTF to both improve the accuracy of the model’s predictive capability and continuing serving the experimental program. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2026

The high explosives & affected targets (HEAT) dataset

Artificial Intelligence (AI) surrogate models offer a computationally efficient alternative to full-physics simulations, yet no existing datasets are publicly available for training, testing, and validation of machine learning models of the dynamics of high-explosive driven shocks through multiple materials. Shock propagation through materials is a computationally challenging problem because simulations must include material-specific equations of state (EOS) along with descriptions of other physical processes such as plastic deformation, phase change, damage processes, fluid instabilities, and multi-material interactions. Shocks are typically initiated by high-velocity impacts or explosive loading. The latter case necessitates the addition of models of reactive materials to represent high-explosive (HE) detonation. Here, to address the lack of an expansive dataset for multi-material shock propagation in the AI/ML community, we present the High-Explosives and Affected Targets (HEAT) Dataset. HEAT is a physics-rich collection of two-dimensional, cylindrically symmetric, simulations generated using an Eulerian, multi-material, shock-propagation code developed at Los Alamos National Laboratory. The dataset includes two partitions: (1) the expanding shock-cylinder (CYL) simulations, Figs. 1, and (2) the Perturbed Layered Interface (PLI) simulations, Fig. 2. Entries in both partitions consist of time series of arrays of thermodynamic fields (pressure, density, and temperature), kinematic fields (position and velocity), and additional fields that depend on thermodynamic and/or kinematic fields (e.g., material stress). Materials in the CYL partition include solids (aluminium, copper, depleted uranium, stainless steel, tantalum, and a generic polymer), a liquid (water), gases (air, nitrogen), and a generic detonating material (high explosive, HE). The PLI partition spans a highly varying geometry but consists of fixed materials across entries: Copper, aluminium, stainless steel, generic polymer, and generic HE. HEAT captures critical phenomena such as momentum transfer, shock propagation, plastic deformation, and thermal effects, making HEAT a valuable benchmark for development of AI/ML emulation of multi-material shock propagation.

36 MATERIALS SCIENCE↗

DOE OSTI2024

High-Fidelity and High-Performance Computational Simulations for Rapid Design Optimization of Sulfur Thermal Energy Storage

Industrial process heating (IPH) accounts for approximately 70% of US manufacturing energy use and is primarily produced by fossil fuel combustion. Approximately 1500 TWht (approximately 60%) of IPH demand is in the temperature range of 100-300. Industrial applications in this temperature range include drying, hydrothermal processing, thermal enhanced oil recovery, food and beverage, bioethanol production, etc. Cost-effective thermal energy storage (TES) that increases the utilization of waste and renewable heat (solar, geothermal, etc.) could provide significant energy savings and reliable heat sources, decrease emissions, and increase US manufacturing competitiveness through reductions in fuel consumption. TES development has historically been dominated by technologies suitable for deployment with concentrating solar power (CSP). State-of-the-art thermal storage deployed commercially with power tower CSP plants uses a 60%/40% NaNO3/KNO3 molten salt and operates between temperatures of approximately 280 degrees Celsius and 570 degrees Celsius using a two-tank configuration. However, these nitrate salts are unsuitable for operation outside of this temperature range due to a high freezing point of approximately 220 degrees Celsius, and limits on high-temperature salt stability and corrosion resistance of containment alloys. Other materials being investigated for TES include those based on: (1) sensible energy storage (various molten salt compositions, inert solid particles, rocks or pebble beds, sulfur, water, concrete, graphite, etc.), (2) latent energy storage in materials that undergo solid-liquid phase change at relevant temperatures (organic materials for low-temperature applications, inorganic salts and/or metals for high-temperature applications), or (3) thermochemical energy storage (hydrides, hydroxides, carbonates, metal oxides, etc.). The application temperature and challenges pertaining to storage material and/or containment cost, energy density, long-term thermal and cyclic stability, and charge/discharge heat transfer effectiveness drive material selection for a given IPH or electricity generation application. Sulfur is a cheap commodity at $80/ton compared to $1100 - 1300/ton for conventional salts. When using a metric of storage cost per kWh, sulfur costs around 2-3 $/kWh. Previous sulfur TES development focused on high temperature (>600 degrees) concentrated solar power applications with sulfur encapsulated in pipes and flow of gaseous HTF (air) in the shell side. However, for lower-temperature IPH applications in the range of approximately 100-300 degrees Celsius Element 16 adopted a compact and scalable TES design with molten sulfur in the shell and HTF pipes submerged in the molten sulfur bath. The low-cost molten sulfur TES for dispatchable IPH has deployment potential for broad applications. The spatial and temporal evolution of the HTF and sulfur temperature is critical to the TES system performance, and thus detailed modeling can improve understanding of the performance and facilitate design improvements. Using high performance computing and computational fluid dynamics (CFD) a low-cost molten sulfur thermal energy storage (TES) system for industrial process heating (IPH) applications was developed. The unique challenges in CFD modeling of sulfur TES are the sharp property changes of sulfur relevant to the working temperatures. Above 159, liquid sulfur undergoes polymerization, and the viscosity of sulfur rapidly increases by several orders of magnitude between 159 degrees Celsius and 188 degrees Celsius, followed by a decrease in viscosity beyond 188 degrees Celsius due to thermal bound dissociation. In addition, various concentrations of H2S impurities can also modify sulfur viscosity. This numerical challenge is especially relevant to transient simulation of the sulfur TES charging and discharging processes as the extreme property variations limit the applicability of traditional heat transfer correlations. Transient CFD simulations including the temperature-dependent sulfur properties and geometric complexity of the TES design were used to predict the effect of natural convection during charging and discharging on the heat transfer process, sulfur temperature uniformity, charge/discharge rates, and performance of the storage devices. The CFD model was validated with experimental results for a full charge and discharge cycle. The work will show 3D and 2D simulation comparisons aimed to facilitate rapid design iterations and a machine learning based design optimization approach.

CFD↗

DOE OSTI2024

Improvement of Coal Power Plant Dry Cooling Technology through Application of Cold Thermal Energy Storage

The U.S. power infrastructure is currently heavily reliant on water cooling. The power plants in the U.S. account for approximately 40% of freshwater withdrawals, with 90% of it used in condenser cooling. The most used cooling technology in coal-fired power plants is a once-through condenser; however, this cooling method requires high water withdrawal rates and results in thermal pollution of the water source. Wet WCTs offer an alternative to once-through condensers due to much lower water withdrawals. However, these systems suffer from water consumption through evaporation making them undesirable options in areas subject to droughts and in arid areas. The direct dry cooled condensers (ACCs) and dry cooling towers (DCT) account for 1.8%, hybrid cooling (ACC + WCT) accounts for 0.5%, while other cooling technologies represent the rest (0.7%). The ACC/DCTs represent an attractive alternative for power plants; however, this technology has not been widely adopted in the US (less than 2% of power plants) due to its negative impact on plant performance. As a rule of thumb, dry cooling results in performance penalty equivalent to approximately 2%-point efficiency loss compared to wet cooling, although the actual magnitude varies with ambient dry bulb temperature, DBT which may vary considerably during the day. As DBT increases, the plant power output decreases, reaching a minimum at the hottest period of the day, which usually coincides with the highest electricity demand for air condition load. Therefore, performance of power plants using DCT/ACC cooling technology is the lowest during the summer mid-day when ambient temperature is the highest. For example, the decrease of the inlet air temperature to the DCT/ACC by 2 Deg C could generate up to 5% additional power at peak demand. It is, therefore, important to improve dry cooling technology to maintain the viability of coal-fired power plants in a carbon constrained future. The method for reducing the cooling air temperature and keeping it constant would mitigate this problem significantly. Objectives of this project were to develop, design, evaluate, and demonstrate a cost-effective system for improving performance of a DCT or ACC for thermal (coal-fired) power plant applications using a low-cost heat storage materials, such as pervious concrete (PC) and phase change material (PCM). Thus, the study focused on development of the system(s) that could be used to alleviate the difficulties in operating DCT/ACC during the summer months by storing cold energy during the nighttime in inexpensive materials PC and PCM and using it during the hottest period(s) of the day. Since very large quantities of cold energy need to be stored to make an impact on performance of a large power plant, it is essential that the storage materials and associated cold energy storage design(s) are inexpensive and the system is simple to build, maintain and operate. To achieve the project objectives, a comprehensive approach, including material development and characterization, component and system modeling, and laboratory- and prototype-scale experiments, was employed including modeling of the system components and of the entire system, development (engineering) of the materials for the heat storage modules of the Cold Thermal Energy Storage System (CTESS) and determination of their properties, design, manufacturing and setup of the laboratory- and prototype-scale test facility and testing, design, manufacturing and setup of the prototype-scale test facility. A modular design of CTESS was employed, where representative modules were designed as the integrated direct contact heat exchanger and thermal energy storage (TES) system. CTESS modules were manufactured and tested. Two prototype-scale designs of the CTESS modules were developed and tested. The use of CTESS increases plant generation increases since it lowers air temperature entering ACC/DCT and keeps it constant during the hottest time of the day. For a PCM-based CTESS, the ambient air temperature is lowered close to the PCM phase change temperature. The duration of the cooling effect depends on the latent heat and mass of PCM in the CTESS. For this project, commercial grade CaCl2 hexahydrate (CaCl2·6H2O or CC6) PCM with phase change temperature of 25 Deg. C was used due to its low cost. For practical reasons, the CTESS was designed to maintain the cooling effect for four hours. The low phase change temperature associated with the commercial grade PCM used in CTESS results in considerably higher improvement in net generation compared to the laboratory (pure) grade. The resistance to heat transfer results in lower net generation compared to the ideal case where resistance to heat transfer is zero. The results demonstrate that CTESS is effective in improving the performance of a dry cooling system. However, its effectiveness depends on the relationship between the ambient air conditions and PCM phase change temperature. As is the case with the heat rejection system, for the best performance, the PCM used in CTESS would need to be matched to the ambient air conditions. The results obtained in this report for selected geographical locations are valid for CC6 and demonstrate that the level of performance to be achieved by the technology will be location-dependent, as is the case with the air cooled condensers. The methodology for engineering of PC-PCM-based heat storage medium is applicable to other PCMs that may need to be used for other ambient air conditions and geographical locations.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2024

The Onset of Nucleate Boiling Prediction in a Mini Channel with the Eulerian Multiphase Flow

The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Thermodynamic investigation of the NaCl-KCl salt system from 25 to 950 °C

NaCl-KCl molten salt system has been proposed as a primary component of several promising heat transfer eutectics and fuel host for Generation IV molten salt reactors (MSR). In this work, several key thermodynamic parameters, including high temperature enthalpies, isobaric heat capacities (C p ) and molar enthalpies of mixing (ΔH mix ), of the NaCl -KCl system (NaCl, 75mol% NaCl – 25mol% KCl, 51mol% NaCl – 49mol% KCl, 25mol% NaCl – 75mol% KCl, and KCl) were measured by high temperature drop calorimetry (HTDC) utilizing laser sealed nickel and aluminum crucibles. Salts were checked for phase purity by X-ray diffraction (XRD), and water content by Karl Fischer coulometric titrimetry and thermogravimetric analysis (TGA). Thermo-mechanical analysis (TMA) utilizing custom-built boron nitride (BN) crucibles was employed to determine temperatures of phase transitions (solid $\rightarrow$ liquid) for the salts and plotted against the pseudobinary phase diagram for the NaCl-KCl system. Small deviations (3-10 %) were found among C p values obtained from HTDC, differential scanning calorimetry (DSC), and NIST, comparable to the average uncertainty based on two standard deviations of measured data that are ~6% for DSC and ~3% for HTDC. Additionally, the measured enthalpies and C p values of mixed NaCl-KCl salts (e.g., 51mol% NaCl – 49mol% KCl eutectic) were shown to behave as statistical mixtures of its pure endmembers, consistent with ΔH mix = ~ 0 kJ/mol at 800 °C, and a regular interaction parameter of Ω = 1.20 ± 0.76 kJ/mol. In conclusion, this study represents a step towards improved accuracy and precision in determinations of thermodynamic parameters of molten salts by using DSC, TMA, and HTDC, in completion of the Molten Salt Thermal Properties Database (MSTDB).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Multi-physics Topology OPtimization and Additive Manufacturing for High-temperature Heat Exchangers

This research significantly advances the understanding of high-temperature heat exchanger design through an integrated approach that combines topology optimization (TO), triply periodic minimal surface (TPMS) structures, additive manufacturing (AM) and thermohydraulic testing. Each of these components contributes uniquely to a unified, high-performance design, fabrication and testing workflow. Topology optimization serves as the foundation of the design methodology by providing a systematic way to determine the most effective material layout for separating hot and cold fluids while maximizing thermal performance. The researchers introduced a novel three-material optimization framework using two density fields to represent hot fluid, cold fluid, and solid domains. This approach enables automated discovery of optimal shapes and flow paths that cannot be intuitively designed, especially under constraints imposed by manufacturing technologies. Furthermore, constraints such as minimal wall thickness and overhang angles were embedded into the optimization process, ensuring that resulting designs are not only thermally efficient but also manufacturable using modern additive techniques. In parallel, the study delves into the use of Gyroid-based TPMS geometries for constructing the core of the heat exchanger. TPMS structures are known for their high surface area, excellent fluid mixing capabilities, and minimal pressure drop characteristics. The researchers applied a data-driven modeling framework using Heteroscedastic Sparse Gaussian Process Regression (HSGPR) combined with genetic algorithms. This allowed for the rapid evaluation and optimization of key geometric parameters such as frequency, iso-value, and phase shift. The result was a set of Gyroid structures tailored for high heat transfer and low flow resistance, demonstrating clear improvements over conventional straight-channel designs. After the designing process, additive manufacturing played a critical role by turning these highly complex, optimized geometries into physical components. Utilizing Laser Powder Bed Fusion (LPBF) with Haynes 282, the study demonstrated the feasibility of fabricating these heat exchangers at high precision. Post-processing methods, including dilation-erosion operations, were applied to ensure local features adhered to self-supporting constraints. The fabricated structures were then subjected to thermohydraulic testing under conditions representative of supercritical CO 2 Brayton cycles, validating the predicted performance and confirming the viability of the full design-to-fabrication pipeline. Finally, thermohydraulic testing across the above studies served as a crucial experimental validation of advanced heat exchanger. Under consistent high-temperature and high-pressure conditions using supercritical CO 2 , the testing demonstrated that both TO and Gyroid-based TPMS designs significantly outperformed conventional straight-channel HXs. The TO design achieved a 115% increase in UA and NTU and a 27.6% boost in gravimetric power density, while the data-driven optimized Gyroid design delivered a 166% increase in UA and NTU and improved effectiveness from 68.7% to 86.1%. These results validate the simulation models, confirm the manufacturability of complex geometries under AM constraints, and provide key insights into design-performance trade-offs, thereby advancing the development of high-efficiency, compact heat exchangers for extreme environments.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Chirality Transfer to the Magnetic Sublattice in the Hybrid Perovskite (R)-/(S)-3-Fluoropyrrolidinium Copper(II) Chloride

Incorporating chiral organic cations into organic–inorganic hybrid materials has been shown to enable the inorganic sublattice to display chiroptical properties. We report a new two-dimensional magnetic (S = 1/2) chiral metal halide perovskite, (R)- and (S)-(C 4 H 9 FN) 2 CuCl 4 (where (C 4 H 9 FN) + is 3-fluoropyrrolidinium), which consists of Cu–Cl inorganic layers separated by (C 4 H 9 FN) + organic cations. The presence of the chiral (C 4 H 9 FN) + organic cation induces the formation of chiral magnetic order, even though the inorganic sublattice itself is nearly structurally centrosymmetric. We also report the racemic variant, containing an equal amount of (R)- and (S)- cations, which shows no evidence of chiral magnetic order. When the magnetic susceptibility is measured perpendicular to the inorganic Cu–Cl layer propagation direction, an antiferromagnetic phase transition at Néel temperature T N = 2.23 K is observed in both the chiral and racemic materials, and the existence of the magnetic phase transition is supported by specific heat capacity measurements. Field-induced magnetic chirality is observed through the existence of a second-order magnetoelectric effect in the chiral variant, while no magnetoelectric signal is observed for the racemic material, indicating the absence of magnetic chirality. Our findings demonstrate that materials exhibiting chiral magnetic order can be created through the incorporation of a chiral cation into an organic–inorganic hybrid magnetic material, potentially allowing for the design of tailored materials that combine chiral magnetism with other desirable optical and electronic properties that come from structural chirality.

Cations↗

DOE OSTI2024

Novel PCM-based fin and tube heat exchanger system for building heating and cooling applications

In the United States, the building sector accounts for 40% of all energy use, and buildings are responsible for more than two-thirds of electricity consumption. Buildings remain the major driver of energy-related carbon emissions, and such emissions are projected to increase in the years ahead because of urbanization and population growth. To accomplish the low carbon energy goal in the building sector, phase change material (PCM)–based thermal energy storage (TES) is increasingly being adopted because it offers several advantages, such as reducing building peak load and energy consumption, enabling large scale deployment of renewables, and improving grid stability. The integration of TES with thermally anisotropic building envelopes (TABEs) is a promising solution because, TABE can redirect natural thermal energy from a building using hydronic loops to TES, and the stored energy can be used later for heating and cooling applications. In this study, we experimentally investigate the thermal performance of a novel fin-tube heat exchanger TES system designed for potential integration with a TABE. An experimental rig for a 5-gal PCM fin-tube heat exchanger system with water as the heat transfer medium is described which records temperatures at the heated and cooled boundaries of the system. Experimental observations provide insight into the role a TES system can play in offsetting a building’s heating and cooling demand and also offer a means to characterize the performance of PCMs for building applications. The scale and study presented in this work will help in the design of future thermal storage systems optimized for storage capacity while also accounting for overall system costs for building applications.

Tamraparni, Achutha↗