SEARCH · Search NASA

Results for “heating elements”

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 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Tailoring Mechanical Properties of a Multi-Principle-Element Alloy by a Multi-Length-Scale Approach

A multiscale strengthening approach was employed to tailor the mechanical properties of a NiCoCr-based multi-principal element alloy (MPEA). Initially, the addition of Mo was utilized to strengthen the alloys by achieving grain size refinement, severe lattice distortion, and a reduction in stacking fault energy (SFE). Subsequently, the Suzuki segregation phenomenon was leveraged by carefully controlling stacking fault (SF) density via pre-straining and a subsequent 500°C heat treatment. Consequently, the Mo-added MPE alloys that were pre-strained and subsequently heat-treated exhibited a true yield stress of 2.3 GPa, with a corresponding fracture elongation of 12%. The Mo promotes complex SF interactions and the SFs formed during prior straining serve as sites for Suzuki segregation during subsequent heat treatment, varying the local SFE with Cr segregation, creating a roughened landscape and hindering dislocation dynamics. This multiscale strengthening approach presents a promising route for achieving desired mechanical properties in MPEAs.

mechanical properties↗

DOE OSTI2025

Development of a thermal creep model for aluminum alloy 6061 cladding in U-10Mo monolithic fuel plates

Plate-type fuel elements consisting of a high-density, low-enriched uranium (LEU) U–10Mo-based fuel foil encapsulated in an aluminum alloy (AA) cladding are fabricated using the hot isostatic pressing (HIP) technique. During the HIP process, the fuel plate system is heated to 560 °C, then cooled to room temperature. This heat cycle significantly affects the mechanical properties of the aluminum cladding, and experimental investigations have shown that, post-HIP bonding, the mechanical properties of the aluminum cladding transition from those of AA 6061-T6 to something closer to the O temper. More specifically, the ultimate strength of the cladding decreases while its ductility increases, making it challenging to capture the changes in mechanical behavior and material properties. Understanding the residual stresses generated during the HIP process is critical for assessing the fuel plate’s integrity under various temperature, pressure, and irradiation. To simulate the HIP bonding process, the elastic, plastic, and thermal properties of the cladding are assumed to be similar to those of AA 6061-O temper. However, the primary challenge lies in the lack of available data for the creep model of the AA 6061 cladding during this transient process of HIP. The present study focuses on developing a computational model that predicts the creep behavior of the aluminum cladding in the fuel plates during the HIP process, as cladding creep significantly influences the residual stresses generated in U-10Mo fuel plates during HIP fabrication. Furthermore, as HIP bonding occurs at high temperatures that are nearing the melting point of aluminum, the present work considered a temperature-dependent Arrhenius-type creep model. In particular, a hyperbolic sine creep model is employed to estimate the creep properties of the as-fabricated aluminum cladding. In conclusion, the residual stresses predicted in the U-10Mo fuel when using the newly calibrated creep model closely align with the experimental measurements, validating the model’s accuracy.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Physics-preserving enriched Galerkin method for a fully-coupled thermo-poroelasticity model

This paper proposes a new numerical method for a fully-coupled, quasi-static thermo-poroelasticity model in a unified enriched Galerkin (EG) method framework. In our method, the mechanics sub-problem is solved using a locking-free EG method, and the flow and heat sub-problems are solved using a locally-conservative EG method. The proposed method offers mass and energy conservation properties with much lower costs than other methods with the same properties, including discontinuous Galerkin methods and mixed finite element methods. The well-posedness and optimal a priori error estimates are carefully derived. Here, several numerical tests confirm the theoretical optimal convergence rates and the mass and energy conservation properties of the new method.

15 GEOTHERMAL ENERGY↗

DOE OSTI2026

Thermal cycling-driven microstructural changes of eutectic Al–Si phase change materials in SS304 containers revealed by multi-modal imaging

Aluminum-based Al–Si alloys are widely used as phase change materials (PCMs) in thermal energy storage (TES) systems owing to their high volumetric latent heat and superior thermal conductivity. However, their long-term reliability is limited by degradation processes that remain insufficiently understood. In this work, we employ a multimodal, correlative characterization framework to systematically resolve the degradation behavior of eutectic Al–Si PCMs in contact with SS304 containers under repeated thermal cycling. By integrating high-resolution electron microscopy, three-dimensional X-ray fluorescence (3D XRF) imaging, and differential scanning calorimetry (DSC), we directly link spatially resolved compositional and microstructural evolution to changes in thermophysical properties. The correlative analysis reveals that elemental leaching of Fe, Cr, and Ni from the stainless-steel container into the PCM drives the formation of intermetallic compounds (IMCs) both at the interface and within the bulk PCM, leading to pronounced compositional heterogeneity. These interfacial reactions and diffusion-induced transformations progressively destabilize the Al–Si eutectic, reducing the effective phase-transforming fraction. Consistent with these observations, DSC measurements show a decrease in melting temperature and latent heat of fusion with thermal cycling. These results underscore the critical influence of interfacial reactions and material compatibility on the stability, durability, and overall performance of Al–Si-based TES systems.

25 ENERGY STORAGE↗

DOE OSTI2025

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗

DOE OSTI2025

Predicting Quantum Criticality in Single-Crystalline Ce 2 Ru 3 Ge 5

Strongly correlated f-electron systems are known to host exotic quantum states, such as quantum criticality, complex order parameters, and unconventional superconductivity. However, the appearance of these exotic states is difficult to predict, making the study of quantum critical behavior challenging, especially in ferromagnetic materials. Herein, we report a structure–property map for Ce 2 M 3 X 5 (M = transition metal; X = main group element) that aids in the targeted design of materials likely to exhibit quantum criticality. Guided by this map, we report on the synthesis of single-crystalline Ce 2 Ru 3 Ge 5 and provide, for the first time, magnetic susceptibility, heat capacity, resistivity, and magnetoresistance measurements on single crystals. Here, we observe a weak ferromagnetic-like response at 7.5 K, which is contrasted with the bulk ferromagnetic ordering that appears in polycrystalline samples. Non-Fermi liquid behavior is seen in the temperature dependent electrical resistivity and heat capacity of the single crystals, suggesting proximity to a ferromagnetic quantum critical point without chemical or physical pressure. Given the contrast with previous reports of polycrystals, these results lead us to propose that single crystalline Ce 2 Ru 3 Ge 5 is intrinsically tuned into the vicinity of a ferromagnetic quantum critical point.

36 MATERIALS SCIENCE↗

DOE OSTI2024

LDRD poster - Identifying the speciation of salt-based actinides in the presence of contaminants

Molten salt reactors (MSRs) are among the most promising Generation IV nuclear reactors. The salts used in these reactors exhibit a high affinity for moisture and oxygen, and heating does not eliminate the oxides that are formed. Additionally, due to the high operating temperatures, reactor parts experience pronounced corrosion, and the resulting corrosion products integrate into the fuel salt system. It is hypothesized that the speciation of f-elements in molten salts would change in the presence of oxides and corrosion products. These changes in local structure can alter the rheological properties, which are critical to the operational and safety parameters of the reactors. Therefore, an in-depth understanding of how corrosion products influence lanthanide/actinide speciation is crucial. To test this hypothesis, lanthanides (as surrogates for actinides) and actinide chlorides were mixed with oxides (such as nickel oxide and lithium oxide) and/or corrosion products (e.g., nickel chloride) and reacted in molten salts. The resulting phases were analyzed using techniques such as powder X-ray diffraction, Raman spectroscopy, UV-Visible spectroscopy, and nuclear magnetic resonance spectroscopy. In the presence of oxides, due to the high oxophilicity of f-elements, these elements capture oxygen from the oxides, forming lanthanide/actinide oxide or oxychloride phases. Specifically, with cerium chloride, the introduction of lithium oxide to the formation of various final products depending on the oxide concentration. At lower oxide concentrations, cerium oxychloride (CeOCl) in the +3 oxidation state was observed, while an increase in the oxide concentration led to the formation of cerium dioxide (CeO2), where cerium is in the +4 oxidation state. In contrast, with nickel oxide, CeOCl formation was consistently observed, regardless of the oxide concentration. This project successfully determined that the speciation of f-elements change in the presence of oxides and different oxides in molten salts result in different final products. When the influence of nickel chloride was examined, changes in f-element speciation were not observed. However, it reacted with lithium chloride to generate a lithium-nickel-chloride (Li6NiCl8) phase, which was also observed when nickel oxide was reacted. Computational models were developed to simulate the phases anticipated under experimental conditions.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

DOE OSTI2024

The Evolution of Galaxies and Clusters at High Spatial Resolution with Advanced X-ray Imaging Satellite (AXIS)

Stellar and black hole feedback heat and disperse surrounding cold gas clouds, launching gas flows off circumnuclear and galactic disks, producing a dynamic interstellar medium. On large scales bordering the cosmic web, feedback drives enriched gas out of galaxies and groups, seeding the intergalactic medium with heavy elements. In this way, feedback shapes galaxy evolution by shutting down star formation and ultimately curtailing the growth of structure after the peak at redshift 2–3. To understand the complex interplay between gravity and feedback, we must resolve both the key physics within galaxies and map the impact of these processes over large scales, out into the cosmic web. The Advanced X-ray Imaging Satellite (AXIS) is a proposed X-ray probe mission for the 2030s with arcsecond spatial resolution, large effective area, and low background. AXIS will untangle the interactions of winds, radiation, jets, and supernovae with the surrounding interstellar medium across the wide range of mass scales and large volumes driving galaxy evolution and trace the establishment of feedback back to the main event at cosmic noon. This white paper is part of a series commissioned for the AXIS Probe mission concept; additional AXIS white papers can be found at the AXIS website.

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

Creep in multi-principal element materials –– A review

The ongoing push towards enhanced energy efficiency and reduced emissions has necessitated the creation of materials with superior performance, especially under extreme conditions. Modern industries, such as aerospace, energy production, and nuclear power, rely heavily on materials that can withstand elevated temperatures without compromising structural integrity. At these heightened temperatures, materials, even when subjected to mechanical stresses well below their yield strength, may experience slow deformation leading to eventual rupture — a phenomenon known as creep. With the expansive design space that comes with the high entropy concept and their reported excellent high temperature strength, multi-principal element materials (MPEMs) have attracted interest in the scientific community for high-temperature applications. Here, this review offers a comprehensive examination of existing studies on creep in MPEMs, which includes multi-principal element−alloys, −bulk metallic glasses, −ceramics, and −superalloys, comparing published findings on MPEMs with pure elements, traditional alloys, bulk metallic glasses, and superalloys. The sub-topics covered include a comparison among different creep-testing methods, creep mechanisms, creep exponents, creep strain rates, activation volume, and creep-activation energy. Modeling efforts for predicting creep behavior of MPEMs are also reviewed. Methods for improving creep resistance by performing heat treatments and/or modifying microstructures are discussed. Overall, the current state of MPEMs has not yet surpassed the creep performance of commercial alloys. Finally, directions for future efforts are suggested, such as experimenting in various controlled environments, expanding the number of compositions tested, exploring advanced manufacturing techniques, and using machine-learning to predict creep properties based on compositions and microstructures.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Synthesis of Si-Fe Chondrule-like Dust Analogues in RF Discharge Plasmas

Chondrules are tiny particles that occur in stony meteorites and are considered as the building blocks of early asteroids and planets. It is believed that they were formed by the fast heating of the dust in the solar nebula. To date, there is no lab-scale experimental study of the formation of chondrules from the initial gas phase precursors following fast heating and crystallisation. The motivation of this work is a pre-trial study of the formation of chnodrule-like particles. The formation of meteorites in the space environment is associated with the aggregation of small particles or molecular clouds under the influence of shock waves or high-energy gas discharges in the solar nebula. In this work, the properties of product formation at the nanoscale-level were investigated using different feedstock materials which are the dominant elements in the meteorite. The structural and morphological properties of the synthesised Si-Fe nanomaterials were analysed by scanning/transmission electron microscopy (SEM/TEM), and chemical composition was analysed by X-ray energy-dispersive spectroscopy (EDS). The identification of crystalline phases was carried out by X-ray diffraction (XRD), whereas the presence of an Fe-Si system in the synthesised particles was demonstrated by Mössbauer spectroscopy. The obtained materials were exposed to the relatively high-energy pulsed plasma beam on the substrate with the aim to emulate the possible fast heating and melting of the formed nanoparticles. The formation steps of growing synthetic (engineered) chondro-like particles and nanostructures in laboratory conditions is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Macro-micro multiscale modeling to assist the design of HPDC Al castings microstructure and alloys for EV super-large body structures (Phase 1)

Implementation of High Pressure Die Casting (HPDC) Aluminum (Al) body structures for high volume electrified vehicles (EV) to improve electric efficiency remains a key strategy within many original equipment manufacturers (OEM)s. In addition to high strength for safety requirements, superior Self Piercing Riveting (SPR) performance is demanded for HPDC Al alloys to be compatible with high volume SPR joining. In this work, it is proposed to extend and validate an existing Contractor finite element multiscale macro-micro modeling approach to quantify the influence of the microstructure of HPDC alloys on the fracture strain/displacement under 3-point bend and clinch testing. The success of this work will allow to replace solution treatment stage with low energy consumption heat treatment (HT) processes, or to design new non heat treatable (NHT) HPDC Al alloys to eliminate HT requirements. Ultimately, this project will facilitate the application of HPDC Al alloys for super-large vehicle structures to significantly reduce vehicle weight and thus improving energy efficiency. The purpose of this project is to extend and validate an existing f

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2024

PIP-II Linac Cryogenic Distribution System Design Challenges

The PIP-II linac Cryogenic Distribution System (CDS) is characterized by extremely small heat inflows and robust mechanical design. It consists of a Distribution Valve Box (DVB), Intermediate Transfer Line, Tunnel Transfer Line, comprising 25 Bayonet Cans, and ends with a Turnaround Can. Multiple helium streams, each characterized by distinct helium parameters, flow through each of these elements. The CDS geometry allows maintaining an acceptable pressure drop for each helium stream, considering the planned flows and helium parameters in different operation modes. This is particularly crucial for the return line of helium vapors, which return from cryomodules to the cold compressors and thus have very restrictive pressure drop requirements. On both sides of the DVB there are fixed supports for process pipes. One of the DVB design challenges was to route the process pipes in such a way that their shape provided sufficient compensation for thermal shrinkage. This ensures th at the forces resulting from thermal shrinkage acting on the cryogenic valves remain at a level acceptable to the manufacturer. The required thermal budget of the CDS was achieved by thermo-mechanical optimization of its components, like process pipes fixed supports in Bayonet Cans.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

MURR LEU structural and thermal hydraulics analyses: Part I – Preliminary irradiation thermo-mechanical behavior

The University of Missouri Research Reactor (MURR) is expected to be converted from highly enriched uranium (HEU, ≥ 20 wt% U-235) U-Al x dispersion fuel to low-enriched uranium (LEU, < 20 wt% U-235) with U-10Mo monolithic fuel. This work introduces high-fidelity irradiation thermo–mechanical (T-M) analysis of the MURR LEU focusing on changes in coolant channel gap thickness. Three-dimensional (3D) finite element (FE) models were developed to simulate the irradiation T-M behavior of the MURR LEU element with all 23 curved fuel plates, the two side plates, and the combs. It was shown that channel gap thickness changes were influenced not only by plate thickness variations due to fuel swelling and creep but also by the radial displacement of consecutive MURR LEU plates. Modeling the fuel element assembly captured side plate displacements, which were shown to reduce radial fuel plate displacements towards the convex side. The maximum local radial displacement in the element was predicted at the end of life (EOL) as 23.7 mil (602.0 µm) on the lateral centerline of plate 23 towards the convex side. The maximum stripe-averaged reduction in channel gap thickness, particularly relevant for thermal hydraulics (TH) safety analysis, was calculated as 15.9 mil (403.9 µm) in single-side heated channel 24 (the outermost channel). These results account for the thermal resistance from the oxide build-up on cladding surfaces which was shown to be up to 0.82 mil (20.8 µm) thick. It was demonstrated that accounting for oxide layer thermal resistance led to a 10 °C higher peak fuel temperature and a 4.4 mil (111.8 µm) greater maximum local radial displacement. The impact of the calculated channel gap thickness changes on the MURR LEU TH safety analysis is evaluated in Part II.

Irradiation thermo-mechanical analysis↗

DOE OSTI2025

Analysis Background & Noise in Stretched Wire Alignment Technique Measurements

The Stretched-Wire Alignment Technique (SWAT) is one method of magnet alignment for linear induction accelerators. The applications of SWAT have been implemented for aligning solenoid magnets on the Scorpius linear induction accelerator which will be sited at the Nevada National Security Site and the Flash X-Ray (FXR) linear induction accelerator at Lawrence Livermore National Laboratory’s Contained Firing Facility. This article describes both systematic (repeatable) and random sources of background and noise as well as practical ways to eliminate or reduce them to acceptable levels. Systematic sources include reflections from wire ends, rapid sag due to ohmic heating of the wire, magnetic materials, and shot rate. Random sources include air currents, vibration of nearby equipment, mechanical stability of test equipment, and the instruments used to measure the wire motion. Mitigations include curve fitting and adaptive noise signal cancellation, and mechanical damping. Finite Element Analysis (FEA) was used to identify and resolve a repeatable wire vibration frequency interfering with the signal resolution. Two stretched wire alignment technique set ups from Sandia National Labs and Lawrence Livermore National Lab have shown background noise sources and ways of mitigating them by either analysis methods or change of mechanical configuration. Conclusions that were drawn included the severe sensitivity of the deflection to even small external interferences of the SWAT wire such that it requires attention to detail in mechanical set up and analysis.

Linear Inductive Accelerator↗

DOE OSTI2024

High-Entropy Alloy R&D for Accelerator Beam Window Applications

High-Entropy Alloys are a class of novel material that can offer improved resistance to beam-induced radiation damage and thermal shock. Development of these new alloys to serve as beam windows in multi-megawatt accelerator target applications is ongoing at Fermilab. Currently we are investigating AlCoCrMnTiV alloy systems of 4-6 components for service as beam windows; these compositions are predicted by CALPHAD simulation to have a low density and single-phase BCC crystal structure. The microstructures of these systems are being studied by electron microscopy techniques such as energy dispersive X-ray spectroscopy (EDS) to determine elemental homogeneity and composition, electron backscatter diffraction (EBSD) to quantify grain structure and orientation, and transmission electron microscopy (TEM) to observe defect structures and precipitate formation; nanoindentation is used to probe microstructural mechanical properties. Initial bulk property characterization utilizes differential scanning calorimetry to quantify specific heat capacity (cp), dilatometry to determine the coefficient of thermal expansion (CTE), and time-domain thermoreflectance (TDTR) to measure thermal conductivity (K). A miniature tensile testing apparatus is also being developed to test tensile properties. Post-irradiation examination is currently ongoing for several of these compositions that have been irradiated by low-energy ions to damage levels and at a temperature relevant to beam window applications at future next-generation accelerator facilities. This talk will briefly describe the alloy design and synthesis before going into more depth covering microstructural pre-characterization, and post-irradiation examination results from low-energy ion irradiated specimens. This will be followed by the plans for alloy down selection and future prototypic proton irradiations.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

2025 Advances in NekRS: Supporting improved performance for nuclear applications

This report presents several 2025 advancements in NekRS, a high-fidelity spectral element CFD code developed at Argonne National Laboratory to support the NEAMS thermal-hydraulics program. The forthcoming v25 release consolidates several of these advances, adding new features for portability across heterogeneous GPU architectures, real-time in situ visualization, improved turbulence modeling, and conjugate heat transfer coupling. Over the past year, NekRS has demonstrated strong scalability and performance on DOE’s leading exascale platforms, including Aurora and Frontier, confirming its readiness for some of the largest and most complex simulations attempted to date. These achievements provide a powerful new platform for high-fidelity data generation, which in turn supports the development and validation of advanced closure models critical for reactor safety and design. Significant algorithmic innovations have also been introduced. A new global runtime h-refinement capability simplifies workflows by reducing mesh preparation burdens and enabling coarse-to-fine restarts. Building on this, a novel multigrid strategy was implemented to accelerate pressure and transport solves at scale, addressing long-standing bottlenecks in exascale CFD. Together, these developments improve both the efficiency and accessibility of high-fidelity simulations for reactor-relevant problems. Collectively, these enhancements represent a major step forward in simulation technology, positioning NekRS as a cornerstone of NEAMS efforts to enable accurate, efficient, and scalable high-fidelity analysis of advanced nuclear systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

A Recipe for a Great Meal: A Benchtop Route from Elemental Se to Superior Thermoelectric β-Ag 2 Se

The low-temperature modification of β-Ag 2 Se has proven to be useful as a near-room-temperature thermoelectric material. Over the past years, research has been devoted to interstitial, vacancy, and substitutional doping into the parent β-Ag 2 Se structure, aiming at tuning the material’s charge and heat transport properties to enhance thermoelectric performance. The transformation of β-Ag 2 Se into α-Ag 2 Se at ~134 °C and the low solubility of dopants are the main obstacles for the doping approach. Herein, we report a facile, safe, scalable, and cost-effective benchtop approach to successfully produce metal-doped β-Ag 2 Se. The doped materials display a remarkable enhancement of thermoelectric performance with a record-high peak zT of 1.30 at 120 °C and an average zT of ~1.15 in the 25–120 °C range for 0.2 at. % Zn-doped Ag 2 Se. The enhancement in zT is attributed to point defects created by Zn doping into Ag vacancies/interstitials, which enhances the scattering of phonons and tunes the charge carrier properties, leading to the significant suppression of thermal conductivity. In conclusion, the simplicity of the synthetic method developed herein and the high performance of the final products provide an avenue to produce high-quality Ag 2 Se-based thermoelectric materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗