Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Critical 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.

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

At least 1,081 records · Page 60

A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival

ABSTRACT The Shox2 homeodomain transcriptional regulator is known for its critical functions during mouse embryogenesis, enabling accurate development of limbs, craniofacial structures, neural populations and the cardiac conduction system. At the genomic level, the Shox2 gene is flanked by an extensive gene desert, a continuous non-coding genomic region spanning over 500 kilobases that contains a multitude of evolutionarily conserved elements with predicted cis -regulatory activities. However, the transcriptional enhancer potential of the vast majority of these elements in combination with the biological necessity of the gene desert have not yet been explored. Using transgenic reporter assays in mouse embryos to validate an extensive set of stringent epigenomic enhancer predictions, we identify several novel gene desert enhancers with distinct tissue-specific activities in Shox2 expressing tissues. 4C-seq chromatin conformation capture further uncovers a repertoire of gene desert enhancers with overlapping activities in the proximal limb, in a compartment essential for Shox2 -mediated stylopod formation. Leveraging CRISPR/Cas9 to delete the gene desert region contained in the Shox2 topologically associated domain (TAD), we demonstrate that this complex cis -regulatory platform is essential for embryonic survival and required for control of region-specific Shox2 expression in multiple developing tissues. While transcription of Shox2 in the embryonic limb is only moderately affected by gene desert loss, Shox2 expression in craniofacial and cardiac domains is nearly abolished. In particular, Shox2 transcripts in the sinus venosus (SV) encompassing the sinoatrial node (SAN) were depleted in embryos lacking the gene desert, likely accounting for the embryonic lethality due to Shox2 -dependency of the SAN pacemaker. Finally, we discover a 1.5kb SV enhancer within the deleted gene desert region, which may act as a genomic module controlling the development of the cardiac conduction system. In summary, our results identify a gene desert indispensable for pleiotropic patterning and highlight the importance of these extensive regulatory landscapes for embryonic development and viability.

Abassah-Oppong, Samuel↗

Perturbative instability of nonergodic phases in non-Abelian quantum chains

An important challenge in the field of many-body quantum dynamics is to identify nonergodic states of matter beyond many-body localization (MBL). Strongly disordered spin chains with non-Abelian symmetry and chains of non-Abelian anyons are natural candidates, as they are incompatible with standard MBL. In such chains, real space renormalization group methods predict a partially localized, nonergodic regime known as a quantum critical glass (a critical variant of MBL). Here, this regime features a treelike hierarchy of integrals of motion and symmetric eigenstates with entanglement entropy that scales as a logarithmically enhanced area law. We argue that such tentative nonergodic states are perturbatively unstable using an analytic computation of the scaling of off-diagonal matrix elements and accessible level spacing of local perturbations. Our results indicate that strongly disordered chains with non-Abelian symmetry display either spontaneous symmetry breaking or ergodic thermal behavior at long times. We identify the relevant length and timescales for thermalization: Even if such chains eventually thermalize, they can exhibit nonergodic dynamics up to parametrically long timescales with a nonanalytic dependence on disorder strength.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Analysis of hydrogen and lithium isotopes using Laser-Induced Breakdown Spectroscopy

The detection and quantification of hydrogen (1,2,3H) and lithium (6,7Li) isotopes are critical to several energy and defense application areas, including hydrogen storage, nuclear forensics, and safeguards/non-proliferation. In this context, laser-induced breakdown spectroscopy (LIBS) is very promising and it comes with certain advantages such as rapid detection, and standoff capability. Although LIBS provides experimental simplicity and is capable of detecting all elements in the periodic table in any phase (solid, liquid, or gas), there exist certain challenges for isotopic analyses of H and Li due to spectral broadening, the presence of closely spaced fine and hyperfine structures, and line distortion effects (e.g., self-absorption and self-reversal) present in laser-produced plasmas. This article reports recent developments in H and Li isotopic analysis using LIBS, and existing challenges.

(020.3260) Isotope shifts, (140.3440) Laser-induce↗

Towards exascale for wind energy simulations

We examine large-eddy-simulation modeling approaches and computational performance of two open-source computational fluid dynamics codes for the simulation of atmospheric boundary layer flows that are of direct relevance to wind energy production. The first code, NekRS, is a high-order, unstructured-grid, spectral element code. The second code, AMR-Wind, is a second-order, block-structured, finite-volume code with adaptive mesh refinement capabilities. The objective of this study is to co-develop these codes in order to improve model fidelity and performance for each. These features will be critical for running ABL-based applications such as wind farm analysis on advanced computing architectures. To this end, we investigate the performance of NekRS and AMR-Wind on the Oak Ridge Leadership Facility supercomputers Summit, using 4 to 800 nodes (24 to 4,800 NVIDIA V100 GPUs), and Crusher, the testbed for the Frontier exascale system, using 18 to 384 Graphics Compute Dies on AMD MI250X GPUs. We compare strong- and weak-scaling capabilities, linear solver performance, and time to solution. We also identify leading inhibitors to parallel scaling.

17 WIND ENERGY↗

Simulation Tools for Characterizing Stress Distribution in Laser Welded Dissimilar Joints

This project focuses on developing a thermo-metallurgical-mechanical modeling method to accurately predict the microstructural evolution and residual stress in laser welding between dissimilar metals, such as HSLA steel and high carbon equivalent (CE) gear steel. The method leverages a comprehensive material database to model the temperature and rate dependent phase transformations, along with their associated effects on material properties, such as thermal expansion and flow stress, throughout the welding process. A key innovation is the incorporation of phase transformation and phase-specific properties, which enhances the accuracy of residual stress predictions. The mixture material in the fusion zone due to the dissimilar metals will also be addressed in the numerical model. This is especially critical in scenarios involving phase transformations in the fusion zone and heat-affected zone (HAZ), where the phase changes can induce substantial residual stress variations. The material database has been generated using JMatPro. The modeling approach is implemented through a custom User Material (UMAT) subroutine, executed with the commercial finite element software Abaqus.

36 MATERIALS SCIENCE↗

Simulation Tools for Characterizing Stress Distribution in Laser Welded Dissimilar Joints

This project focuses on developing a thermo-metallurgical-mechanical modeling method to accurately predict the microstructural evolution and residual stress in laser welding between dissimilar metals, such as HSLA steel and high carbon equivalent (CE) gear steel. The method leverages a comprehensive material database to model the temperature and rate dependent phase transformations, along with their associated effects on material properties, such as thermal expansion and flow stress, throughout the welding process. A key innovation is the incorporation of phase transformation and phase-specific properties, which enhances the accuracy of residual stress predictions. The mixture material in the fusion zone due to the dissimilar metals will also be addressed in the numerical model. This is especially critical in scenarios involving phase transformations in the fusion zone and heat-affected zone (HAZ), where the phase changes can induce substantial residual stress variations. The material database has been generated using JMatPro. The modeling approach is implemented through a custom User Material (UMAT) subroutine, executed with the commercial finite element software Abaqus.

36 MATERIALS SCIENCE↗

A Comprehensive Review of Selected Major Categories of Lithium Isotope Separation Techniques

The ability to generate materials with an enriched isotopic abundance, that is one that differs from natural abundance, is critical in terms of technology. Isotopes are used in a wide range of applications, including radioactive tracers, nondestructive tests, radiotherapy in human medicine, nuclear fuels, and isotope‐substituted compounds for chemistry and biology research, as well as environmental geochemical signature tracers. A commonly discussed research topic is the isotopic separation of various elements—uranium, hydrogen, lithium, and iodine, among many others. Among these isotopes, lithium isotopes are especially important in various applications including strategic areas. These isotopes are light and hence their separation techniques are not very obvious. Herein, all the main categories of Li isotope separation are reviewed in depth with a focus on electrochemical technologies for stable lithium isotope separation.

Murali, Arun↗

Bayesian sequential optimal experimental design for nonlinear models using policy gradient reinforcement learning

We present a mathematical framework and computational methods for optimally designing a finite sequence of experiments. This sequential optimal experimental design (sOED) problem is formulated as a finite-horizon partially observable Markov decision process (POMDP) under a Bayesian setting and with information-theoretic utilities. The formulation is general and may accommodate continuous random variables, non-Gaussian posteriors, and nonlinear forward models. The sOED design policy incorporates elements of feedback and lookahead simultaneously, and we show it to generalize the commonly-used batch and greedy design strategies. We solve for the sOED policy using the policy gradient (PG) method from reinforcement learning, and provide a derivation for the PG expression in the sOED context. Adopting an actor-critic approach, the policy and value functions are parameterized using deep neural networks and improved via PG estimates produced from simulated episodes of designs and observations. The new PG-sOED algorithm is first validated on a linear-Gaussian benchmark, and then compared against other design baselines on a sensor movement problem for contaminant source inversion in a convection-diffusion field. As a result, we provide explanation for the policy behaviors using knowledge of the underlying physical process.

97 MATHEMATICS AND COMPUTING↗

Unveiling damage mechanisms of chromium-coated zirconium-based fuel claddings at LWR operating temperature by in-situ digital image correlation

Here, we investigate the coupled thermomechanical fracture mechanisms of coated nuclear fuel claddings at Light Water Reactor (LWR) operating temperatures. These coated claddings are a highly attractive, near-term solution, which addresses the demands for accident-tolerant fuel systems and provide greater oxidation resistance. However, the fracture mechanisms of these coatings, which may create channels for oxidation ingression, must be fully understood prior to implementation. Thus, high-temperature expanding plug experiments were conducted on coated cladding specimens at a temperature of 315 °C, consistent with the operating environment of LWRs. Furthermore, in-situ thermomechanical deformation was measured with stereo digital image correlation during heating and mechanical testing to separately resolve contributions of thermal and mechanical strain. Digital image correlation, supported by acoustic emissions (AE) detection, was also leveraged to track cracking activity during loading. Coating fracture was found to initiate at total hoop strains of 0.34%. The thermal deformation of the coated claddings was investigated via finite element simulations, revealing a bidirectional stress-state within the coating with axial and circumferential strains of 0.026 and 0.031%. This bidirectional stress-state was attributed with the generation of off-axis fracture pattern within the coating as identified via post-experiment scanning electron microscopy. Thus, this study unveiled critical, coupled thermomechanical mechanisms governing the coating fracture of coated claddings at LWR temperatures.

36 MATERIALS SCIENCE↗

Tailoring the Structure of Carbon Molecular Sieves Derived from an Aromatic Polyamide

Energy-efficient and selective separation technologies are required to reclaim aqueous waste streams for reuse. Carbon molecular sieves (CMS) are one material capable of separating organic solvent and gaseous mixtures due to the presence of permanent, rigid, and molecular sieving pores and could potentially be used in aqueous separations. Here, a novel CMS material is derived from an aromatic polyamide precursor with tailored fabrication techniques. Structural characteristics of this CMS have been probed with various analyses and are potentially ideal for selective separations as adsorbents. In this work, elemental analyses provide insight about the potential chemical structures of CMS materials that result from degradation mechanisms. Gravimetric sorption data is used to estimate sorption, diffusion, and permeability selectivity with water and N, N-dimethylformamide. While challenges to the development and implementation of such carbonaceous materials remain, highly selective materials are critical starting points for enabling ultra-challenging separation processes of the future.

36 MATERIALS SCIENCE↗

Challenges and Gaps in the Development of Pulsed Power for Fusion Applications: A Preroadmapping Perspective From Industry, Academia, and National Laboratory Experts

Fusion energy meets the twenty-first century World Grand Challenge of sustainable, ubiquitous, and safer energy sources. However, harnessing the promise of fusion energy has proven elusive. The competing approaches to fusion power plant design include inertial confinement fusion, National Ignition Facility (ICF-NIF, Z machine, etc.,) magnetic confinement fusion (MCF-Tokamak, stellarators, etc.), and other approaches that show promise in small- (flow stabilized Z pinches) or large-scale applications. These approaches are being accelerated with private and public funding and seek to demonstrate the feasibility of different approaches to fusion-based power plants. Yet, how can the necessary pulsed power technologies for these disruptive technology bases be accelerated with no clear “Dominant Design?” Roadmapping holds the promise to identify and develop common critical pulsed power components for laboratory, prototype, and commercial fusion, and can accelerate the commercialization of fusion reactor designs. A preroadmapping Workshop on Pulsed Power for Fusion was held at the IEEE International Pulsed Power Conference in San Antonio, TX, USA, in June 2023. The workshop had 177 attendees. Here, the common elements for many of the ICF technologies vying for dominant design were identified. The advancement of these technologies through roadmapping will enhance commercial expectations that require their rapid and innovative development in the next five years, as well as the next five to ten years. The key technologies identified that underpin and limit the advancement of fusion power include pulsed power technologies such as energy storage, high-voltage switching, additive manufacturing, and modular pulsed power circuit topologies. In conclusion, they are the focus of our effort in the following roadmap scenario, which will delineate potential paths to technology development.

Curry, Randy D. [I-Pulse Group, Albuquerque, NM (U↗

High-Power Electric Vehicle Charging Hub Integration Platform (eCHIP): Site Energy Management System Platform Development

A site energy management system is a critical component of future high-power EV charging hubs. It monitors and supervises all inter-hub operations, coordinates grid interactions, and implements energy management strategies and load-sharing controls at various levels and time scales. This report provides a comprehensive high-level overview and guidelines on the essential elements of such energy management systems. These elements primarily focus on the core functions, architecture, and monitoring aspects of site energy management. The core functions define the grid services and operational objectives that a charging hub is expected to provide and achieve. The architectural section investigates the coordination between hub assets, different control schemes, and their potential benefits and drawbacks. Lastly, the report discusses the needs, types, and functionalities of a monitoring and visualization system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Computational Modeling of Graphite Degradation due to Molten Salt Infiltration and Wear

Molten-salt reactors (MSRs) represent a promising next-generation reactor design, with graphite serving as a moderator and/or reflector in several designs. However, due to limited experimental data and operational experience, a technical understanding of the structural integrity of graphite in molten salt environments remains incomplete. This report presents a modeling-based evaluation of graphite degradation in MSR environments, focusing on the effects of salt infiltration in fuel salt-based designs and surface wear in pebble bed reactor designs. The objective of this study is to enhance understanding of the structural integrity challenges posed by these degradation mechanisms and to provide a framework for assessing graphite behavior in MSRs. The first part of the report investigates the phenomenon of molten salt infiltration into graphite. This infiltration occurs when molten salt permeates the interconnected pore structure of the graphite moderator, driven by factors such as pressure differentials and the physical properties of both the salt and graphite. The infiltration process is influenced by characteristics of the pore structure, viscosity of the molten salt, and the interfacial energies between the graphite, salt, and the atmosphere within the graphite pore. Utilizing a coupled multiphysics modeling approach with Grizzly software, the study evaluates the stress induced by internal heat sources due to infiltration, which can lead to structural concerns. This evaluation is crucial for understanding how infiltration affects the mechanical integrity of graphite components in MSRs. The study considers the Molten-Salt Reactor Experiment (MSRE) graphite stringer geometry due to the availability of relevant data. Through detailed finite element analysis, the study examines stress distributions at varying infiltration percentages, revealing that stress levels increase with higher amounts of infiltration. Rare-event simulations, using the parallel subset simulation (PSS) framework, further quantify the failure probabilities under input uncertainties, with a user-specified failure metric. The PSS framework also identifies critical input parameters that significantly affect the stress values, including infiltration amount, thermal conductivity, and power density. Additionally, considering realistic reactor scenarios, the analysis was performed to account for the combined effects of radiation and infiltration, and modeling strategies on how to analyze new reactor designs or new graphite grades are discussed. The second part of the report focuses on wear mechanisms in pebble bed-based MSRs. As graphite fuel pebbles interact with the graphite reflector block, wear can result in material loss and the formation of surface defects, which may act as stress concentrators. A similar multiphysics modeling framework is employed to assess the impact of wear on the structural integrity of graphite components. This study considers a generic fluoride-cooled high-temperature reactor (gFHR) design due to the availability of comprehensive data. Worst-case scenario dimensions of the reflector blocks were analyzed under thermal and radiation conditions. Subsequently, wear in the form of idealized pits and grooves is modeled on the inner surface of the graphite block, with the maximum stress from previous simulations. The simulations show that groove-type defects are more detrimental than pits, leading to higher stress concentrations. Considering worst-case simulation scenarios and experimental wear rates, it was determined that the formation of a surface defect critical enough to affect the stress may not be possible in a gFHR design. Overall, the findings of this research contribute to the development of robust modeling tools for predicting graphite behavior under various operational conditions in MSRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

History of the Chemical Stockpile Emergency Preparedness Program – Volume 2: Lessons Learned (Final Report 2026)

The CSEP Program involved multiple federal agencies, multiple states, and numerous local units of government. The involvement of diverse agencies at different levels of government led to conflict and difficulty in resolving issues. That, in turn, led to slow progress in achieving program goals. Program management was criticized repeatedly by the U.S. General Accounting Office (GAO). After other mechanisms for interagency coordination had failed to give the desired results, FEMA and the Army agreed to form site-specific and national-level IPTs. The IPTs worked well to break the logjam and accelerate implementation of the program. Key elements of the IPTs’ success included constant interagency communication, solving problems at the lowest level, and obtaining buy-in from higher authority. The IPT technique could be applied to other emergency management and homeland security programs that require interagency coordination.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Degradation science: Integrating modeling and experiments to predict localized corrosion processes (Annual Progress Report)

Additively manufactured (AM) eutectic high-entropy alloys (EHEAs), such as nano-lamellar AlCoCrFeNi 2.1 , have excellent strength, ductility, and wear resistance even at elevated temperatures, but their corrosion behavior in aggressive acids at different length scales remain poorly understood. This work investigates the corrosion behavior of laser powder bed–fused (L-PBF) AlCoCrFeNi 2.1 as a function of annealing temperatures, probing degradation mechanisms from nanoscopic to macroscopic length scales. The alloy is dual phase consisting of a ductile FCC L1 2 phase and a high-strength BCC B2 phase. Rapid solidification during L-PBF produces a far-from-equilibrium nano-lamellar structure with nearly homogeneous elemental distribution, which tends to evolve upon annealing toward Cr/Co/Fe-enriched FCC and Al/Ni-enriched B2. Three conditions were studied: as-printed, 600 °C/5 h, and 1000 °C/1 h, over which B2 lamellae coarsen, lamellar spacing increases, and elemental segregation becomes more prominent. Microstructure and chemistry were characterized by scanning electron microscopy (SEM) and energy-dispersive spectroscopy (EDS), while in-situ electrochemical atomic force microscopy (EC-AFM) was used to link early (<5 h) local dissolution to microstructure after exposure in sulfuric acid. EC-AFM highlights preferential dissolution of the BCC/B2 phase where the surrounding matrix is Cr-depleted and directly quantifies the dissolution rates within individual phases, tracks the transition from early nano-scale attack to partial repassivation, to correlate height differences with current and impedance responses. To monitor longer-term behavior (up to 96 h), ex-situ AFM, SEM, and confocal imaging were combined with conventional bulk electrochemical tests, bridging nanoscale observations to micro/meso-scale damage morphologies. At the meso-scale, the deepest dissolution channels align with the build-direction lamellae and melt-pool boundaries, indicating that printing directionality guides the propagation of these localized corrosion sites. Annealing modifies corrosion by restructuring BCC/FCC phase fractions, lamellar spacing, and Cr/Al segregation, thereby changing the cathode/anode ratio and passive film stability. The results clarify how as-printed and annealed nano-lamellar architectures differ in their susceptibility to selective dissolution; how elemental segregation competes with residual stresses along the build direction. With these insights, future work will use CALPHAD-guided alloy modification to stabilize higher Cr contents in the B2 phase while retaining a dominant FCC+B2/BCC microstructure, with the goal of designing mechanically robust, corrosion-resistant EHEAs for safety-critical applications to leverage the LLNL’s broader national and global security mission.

36 MATERIALS SCIENCE↗

Global 10 m Land Use Land Cover Datasets: A Comparison of Dynamic World, World Cover and Esri Land Cover

The European Space Agency’s Sentinel satellites have laid the foundation for global land use land cover (LULC) mapping with unprecedented detail at 10 m resolution. We present a cross-comparison and accuracy assessment of Google’s Dynamic World (DW), ESA’s World Cover (WC) and Esri’s Land Cover (Esri) products for the first time in order to inform the adoption and application of these maps going forward. For the year 2020, the three global LULC maps show strong spatial correspondence (i.e., near-equal area estimates) for water, built area, trees and crop LULC classes. However, relative to one another, WC is biased towards over-estimating grass cover, Esri towards shrub and scrub cover and DW towards snow and ice. Using global ground truth data with a minimum mapping unit of 250 m 2 , we found that Esri had the highest overall accuracy (75%) compared to DW (72%) and WC (65%). Across all global maps, water was the most accurately mapped class (92%), followed by built area (83%), tree cover (81%) and crops (78%), particularly in biomes characterized by temperate and boreal forests. The classes with the lowest accuracies, particularly in the tundra biome, included shrub and scrub (47%), grass (34%), bare ground (57%) and flooded vegetation (53%). When using European ground truth data from LUCAS (Land Use/Cover Area Frame Survey) with a minimum mapping unit of <100 m 2 , we found that WC had the highest accuracy (71%) compared to DW (66%) and Esri (63%), highlighting the ability of WC to resolve landscape elements with more detail compared to DW and Esri. Although not analyzed in our study, we discuss the relative advantages of DW due to its frequent and near real-time data delivery of both categorical predictions and class probability scores. We recommend that the use of global LULC products should involve critical evaluation of their suitability with respect to the application purpose, such as aggregate changes in ecosystem accounting versus site-specific change detection in monitoring, considering trade-offs between thematic resolution, global versus. local accuracy, class-specific biases and whether change analysis is necessary. We also emphasize the importance of not estimating areas from pixel-counting alone but adopting best practices in design-based inference and area estimation that quantify uncertainty for a given study area.

54 ENVIRONMENTAL SCIENCES↗

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↗