Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Temperature-Dependent”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Numerical Simulation of Electron Magnetohydrodynamics with Landau-quantized Electrons in Magnetar Crusts

Abstract In magnetar crusts, magnetic fields are sufficiently strong to confine electrons into a small to moderate number of quantized Landau levels. This can have a dramatic effect on the crust's thermodynamic properties, generating field-dependent de Haas–van Alphen oscillations. We previously argued that the large-amplitude oscillations of the magnetic susceptibility could enhance the ohmic dissipation of the magnetic field by continuously generating small-scale, rapidly dissipating field features. This could be important to magnetar field evolution and contribute to their observed higher temperatures. To study this, we performed quasi-3D numerical simulations of electron MHD in a representative volume of neutron star crust matter, for the first time including the magnetization and magnetic susceptibility resulting from Landau quantization. We find that the potential enhancement in the ohmic dissipation rate due to this effect can be a factor ∼3 for temperatures of the order of 10 8 K, and ∼4.5 for temperatures of the order of 5 × 10 7 K, depending on the magnetic field configuration. The nonlinear Hall term is crucial to this amplification: without it, the magnetic field decay is only enhanced by a factor ≲2 even at 5 × 10 7 K. These effects generate a high wavenumber plateau in the magnetic energy spectrum associated with the small-scale de Haas–van Alphen oscillations. Our results suggest that this mechanism could help explain the magnetar heating problem, though due to the effect's temperature-dependence, full magneto-thermal evolution simulations in a realistic stellar model are needed to judge whether it is viable explanation.

Rau, Peter B. (ORCID:0000000152209277)↗
DOE OSTI2021

Astromers in the Radioactive Decay of r-process Nuclei

Certain nuclear isomers are well known to affect nucleosynthesis with important observable consequences (e.g., 26 Al and 180 Ta). We study the impact of nuclear isomers in the context of rapid neutron capture process (r-process) nucleosynthesis. Here we demonstrate that nuclear isomers are dynamically populated in the r process and that some are populated far from thermal equilibrium; this makes them astrophysical isomers, or "astromers." We compute thermally mediated transition rates between long-lived isomers and the corresponding ground states in neutron-rich nuclei. We calculate the temperature-dependent β-decay feeding factors, which represent the fraction of material going to each of the isomer and ground state daughter species from the β-decay parent species. We simulate nucleosynthesis following the decay of a solar-like r-process composition and include as separate species nuclear excited states with measured terrestrial half-lives greater than 100 μs. We introduce a new metric to identify those astromers most likely to be influential and summarize them in a table. Notable entries include many second peak nuclei (e.g., the Te isotopic chain) and previously overlooked isomers in stable nuclei (e.g., 119 Sn, 131 Xe, and 195 Pt). Finally, we comment on the capacity of isomer production to alter radioactive heating in an r-process environment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

AlGaN High Electron Mobility Transistor for High-Temperature Logic

Here we report on AlGaN high electron mobility transistor (HEMT)-based logic development, using combined enhancement- and depletion-mode transistors to fabricate inverters with operation from room temperature up to 500°C. Our development approach included: (a) characterizing temperature-dependent carrier transport for different AlGaN HEMT heterostructures, (b) developing a suitable gate metal scheme for use in high temperatures, and (c) over-temperature testing of discrete devices and inverters. Hall mobility data (from 30°C to 500°C) revealed the reference GaN-channel HEMT experienced a 6.9x reduction in mobility, whereas the AlGaN channel HEMTs experienced about a 3.1x reduction. Furthermore, a greater aluminum contrast between the barrier and channel enabled higher carrier densities in the two-dimensional electron gas for all temperatures. The combination of reduced variation in mobility with temperature and high sheet carrier concentration showed that an Al-rich AlGaN-channel HEMT with a high barrier-to-channel aluminum contrast is the best option for an extreme temperature HEMT design. Three gate metal stacks were selected for low resistivity, high melting point, low thermal expansion coefficient, and high expected barrier height. The impact of thermal cycling was examined through electrical characterization of samples measured before and after rapid thermal anneal. The 200-nm tungsten gate metallization was the top performer with minimal reduction in drain current, a slightly positive threshold voltage shift, and about an order of magnitude advantage over the other gates in on-to-off current ratio. After incorporating the tungsten gate metal stack in device fabrication, characterization of transistors and inverters from room temperature up to 500°C was performed. The enhancement-mode (e-mode) devices’ resistance started increasing at about 200°C, resulting in drain current degradation. This phenomenon was not observed in depletion-mode (d-mode) devices but highlights a challenge for inverters in an e-mode driver and d-mode load configuration.

42 ENGINEERING↗
DOE OSTI2025

Axion Perturbations: A General Analytical Treatment

Cosmological data provides us two key constraints on dark matter (DM): it must have a particular abundance, and it must have an adiabatic spectrum of density perturbations in the early universe. Many different cosmological scenarios have been proposed that establish the abundance of axion DM in qualitatively different ways. In this paper we emphasize that, despite this variety of backgrounds, the perturbations in axion DM can be understood from universal principles. How does a feebly interacting axion field acquire perturbations proportional to those of photons? How do the isocurvature power spectrum and non-Gaussianity depend on the background evolution of the universe? We answer these questions for a completely general choice of cosmological background and temperature-dependent axion potential. We show that the most general solution to the axion field equation on super-horizon scales is entirely determined by the family of background solutions for different initial field values . This holds for both the component in the field perturbation solution contributing to the DM isocurvature perturbation (enhanced at late times by the sensitivity of the DM abundance to the initial condition, , which can be large for initial conditions near the hilltop), and the other component that contributes to the DM curvature perturbation. In particular, we explain that an unperturbed axion field in the early universe evolving into one with nontrivial adiabatic perturbations is guaranteed by Weinberg's theorem on adiabatic modes. These results have been derived before with various assumptions, such as a radiation dominated background or a quadratic potential. Our aim is to give a clear, simple derivation that is manifestly independent of those assumptions, and thus can be applied to any cosmological axion scenario.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗
DOE OSTI2022

Source apportionment of VOCs, IVOCs and SVOCs by positive matrix factorization in suburban Livermore, California

Abstract. Gas- and particle-phase molecular markers provide highly specific information about the sources and atmospheric processes that contribute to air pollution. In urban areas, major sources of pollution are changing as regulation selectively mitigates some pollution sources and climate change impacts the surrounding environment. In this study, a comprehensive thermal desorption aerosol gas chromatograph (cTAG) was used to measure volatile, intermediate-volatility and semivolatile molecular markers every other hour over a 10 d period from 11 to 21 April 2018 in suburban Livermore, California. Source apportionment via positive matrix factorization (PMF) was performed to identify major sources of pollution. The PMF analysis identified 13 components, including emissions from gasoline, consumer products, biomass burning, secondary oxidation, aged regional transport and several factors associated with single compounds or specific events with unique compositions. The gasoline factor had a distinct morning peak in concentration but lacked a corresponding evening peak, suggesting commute-related traffic emissions are dominated by cold starts in residential areas. More monoterpene and monoterpenoid mass was assigned to consumer product emissions than biogenic sources, underscoring the increasing importance of volatile chemical products to urban emissions. Daytime isoprene concentrations were controlled by biogenic sunlight- and temperature-dependent processes, mediated by strong midday mixing, but gasoline was found to be the dominant and likely only source of isoprene at night. Biomass burning markers indicated residential wood burning activity remained an important pollution source even in the springtime. This study demonstrates that specific high-time-resolution molecular marker measurements across a wide range of volatility enable more comprehensive pollution source profiles than a narrower volatility range would allow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A chemical kinetics theory for interpreting the non-monotonic temperature dependence of enzymatic reactions

Abstract. One notable observation of enzymatic chemical reactions is that, for a given abundance of enzymes and substrates, temperature increases cause reaction rates to first increase consistent with the Arrhenius relationship, then plateau, and finally fall off quickly to zero at high temperatures. While many mathematical functions have been used to describe this pattern, we here propose a chemical kinetics theory which successfully replicates this observation and provides insights into the processes responsible for these dynamics. The chemical kinetics theory combines the law of mass action, von Smoluchowski's diffusion-limited chemical reaction theory, and Eyring's transition state theory. This new theory reveals that the thermally reversible enzyme denaturation ensured by the ceaseless thermal motion of molecules and ions in an enzyme solution explains the plateau and subsequent decrease in chemical reaction rates with increasing temperature. The temperature-dependent affinity parameter (K) that relates enzymes and substrates through their binding also affects the shape of the emergent temperature response. We demonstrate that with an increase in substrate availability, K shifts the optimal temperature, where reaction rates plateau, towards higher values. Further, we show that the chemical kinetics theory accurately represents 12 sets of published enzyme assay data and includes the popular mechanistic model by Ratkowsky et al. (2005) as a special case. Given its good performance and solid theoretical underpinning, we believe this new theory will facilitate the construction of more mechanistic-based environmental biogeochemical models.

Tang, Jinyun↗
DOE OSTI2022

A New Search for Low-mass Dark Matter and an Examination and Reduction of the Uncertainty due to the Photoelectric Absorption Cross Section using a Cryogenic Silicon Detector with Single-charge Sensitivity

Many astrophysical observations point to the abundant existence of dark matter (DM) in the universe composed of particles beyond the Standard Model. However despite numerous experiments using different techniques, DM particles have yet to be directly observed. The Super Cryogenic Dark Matter Search (SuperCDMS) SNOLAB experiment will employ silicon and germanium crystal detectors operated at temperatures as low as 15~mK to probe DM interactions using phonon and ionization signals. Recently, R\&D; facilities have developed gram-sized high-voltage eV-scale (HVeV) silicon detectors that achieve single-electron-hole-pair resolution. By probing effective DM-electron interactions, these devices can be used for searches of low-mass DM candidates.This dissertation presents a DM search experiment known as HVeV Run~2 that employs a second-generation HVeV detector operated in an above-ground laboratory at Northwestern University (IL, USA). Energy spectra are obtained from a blind analysis with 0.39 and 1.2\,g-days of exposure with the detector biased at 60 and 100\,V, respectively. The 0.93~gram detector achieves a 3\,eV phonon energy resolution, corresponding to a world-leading charge resolution of 3\,\% of a single electron-hole pair for a detector bias of 100\,V. With charge carrier trapping and impact ionization effects incorporated into the DM signal models, the resulting exclusion limits are reported for inelastic DM-electron scattering for DM masses from 0.5--$10^{4}$\,MeV$/c^{2}$; in the mass range from 1.2--50\,eV$/c^{2}$ the limits for dark photon and axion-like particle absorption are reported.Several DM search experiments, including HVeV Run~2, are sensitive to low-mass DM candidates that rely on the temperature-dependent photoelectric absorption cross section of silicon. However discrepancies in the underlying literature data result in dominating systematic uncertainties on the DM exclusion limits. In order to reduce these systematic uncertainties, this dissertation presents a novel method of making a direct, low-temperature measurement of the photoelectric absorption cross section of silicon at energies near the band gap (1.2--2.8\,eV).

Wilson, Matthew James↗
DOE OSTI2022

Development of a Laser Ultrasonics-based Approach for Rapid Screening of High Entropy Alloys

This project utilized a laser ultrasonic technique to systematically study temperature-induced evolution of material properties in a set of interrelated binary alloys and a high entropy alloy (HEA) fabricated using arc-melting and spark plasma sintering processes. This technique involved the use of a nanosecond duration, high-intensity pulsed laser to thermo-elastically generate ultrasonic waves that propagate in the bulk of the metal alloy. Sub-nanometer-scale displacements associated with the propagating bulk ultrasonic waves were detected along the epicenter on the opposite surface of the sample using a 1 GHz bandwidth photorefractive interferometer. Phase transformations and microstructural changes were inferred from the temperature-dependent trends in the bulk acoustic velocities and features in the ultrasonic epicentral waveforms measured in the binary alloys. These inferences were then correlated with electron/optical microscopy observations and predictions using calculations of phase diagrams (CALPHAD). The results showed that the laser-generated ultrasonic pulses were strongly influenced by changes in material microstructure and could accurately track thermally driven phase transformations and detect the presence of microscale heterogeneities (grain boundaries, dendritic structures, etc.) in the set of binary alloy samples. The measurement approach was then applied to a quinary HEA sample for estimating phase transition temperature and determining microstructural heterogeneity. The rapid, non-contact and non-destructive ultrasonic testing approach demonstrated here is amenable to high throughput combinatorial investigations that can be applied to graded composition HEAs produced using advanced manufacturing methods. When paired with atomistic simulations and CALPHAD modeling, this approach can overcome the bottlenecks faced by current material characterization methods in efficiently screening the vast discovery space of HEAs that spans over a hundred million unique quinary alloy compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Impact of Irradiation-induced Lattice Defects and Dopants on Thermal Transport in Oxide Nuclear Fuels

Advanced nuclear reactor concepts aim to use fuels that must withstand unprecedented temperature and radiation extremes. In these fuels, thermal energy transport under irradiation is directly related to fuel longevity, reactor safety, and is arguably one of the most important performance metrics. Here we provide a comprehensive, first-principles-informed treatment of phonon mediated thermal transport in a defect-bearing actinide oxide with direct comparison to experimental measurements. Pristine and proton irradiated thorium dioxide was chosen as a model system to treat the complexity of thermal transport in the presence of lattice defects. A thermal transport model is implemented using the linearized Boltzmann transport equation (LBTE) with input from first principles calculations and defect evolution models. The output of the LBTE is compared directly to mesoscopic measurements of thermal conductivity on length scales commensurate with defect accumulation. Parametric measurements of conductivity with irradiation dose and temperature reveal a monotonic decrease in conductivity with irradiation dose, in the range of 0.001 dpa to 0.1 dpa. Besides radiation-induced defects, doping the thorium dioxide crystal with small quantities of uranium atoms also results in a drastic reduction in thermal conductivity. Temperature-dependent measurements of thermal conductivity in uranium-doped thorium dioxide single crystals suggests a resonant scattering mechanism is responsible for the observed reduction. This comprehensive, atomistic- to meso-scale treatment provides the necessary basis to investigate thermal transport under irradiation in more complex systems that exhibit strong electron correlation. This is a poster presentation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Temperature Dependence of Band Gap Renormalization in High-T Sensor Materials via First-Principles and Experimental Corroboration

Understanding the temperature dependence of functional properties of high-T gas sensing materials is vital for their applications in combustion environments. The electron-phonon coupling that derives the electronic structure change with temperatures is a key property of interest as it affects other sensing responses. Herein, we assess the temperature dependence of band gap renormalization in metal oxides and perovskites by employing Allen-Heine-Cardona theory with first-principles simulations and corroborate with experimental observation. The calculated temperature-dependent band gap changes of these materials studied are in good agreement with in-house experimental data, proving that the theory can adequately predict renormalization on the band gap in the system of interest. The predicted and measured band gap variations are characterized using an analytical model, which can provide useful insights on the simulated zero-temperature band gaps. Based on the available data, a set of 53 metal oxides and perovskites were identified as potential high-T gas sensors. A machine learning model has been developed to predict the band-gap change by capturing the overall trend of the empirical parameters with respect to a reduced feature obtained by transforming the set of available physical features.

Park, Jongwoo↗
DOE OSTI2023

The Collaborative Effects of Intrinsic and Extrinsic Impurities in Low RRR SRF Cavities

The superconducting radio-frequency (SRF) community has shown that introducing certain impurities into high-purity niobium can improve quality factors and accelerating gradients. We question why some impurities improve RF performance while others hinder it. The purpose of this study is to characterize the impurity profile of niobium with a low residual resistance ratio (RRR) and correlate these impurities with the RF performance of low RRR cavities so that the mechanism of impurity-based improvements can be better understood and improved upon. The combination of RF testing and material analysis reveals a microscopic picture of why low RRR cavities experience low temperature-dependent BCS resistance behavior more prominently than their high RRR counterparts. We performed surface treatments, low temperature baking and nitrogen-doping, on low RRR cavities to evaluate how the intentional addition of oxygen and nitrogen to the RF layer further improves performance through changes in the mean free path and impurity profile. The results of this study have the potential to unlock a new understanding on SRF materials and enable the next generation of SRF surface treatments.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Exascale Multiphysics Nuclear Reactor Simulations for Advanced Designs

ENRICO is a coupled application developed under the U.S. Department of Energy's Exascale Computing Project (ECP) targeting the modeling of advanced nuclear reactors. It couples radiation transport with heat and fluid simulation, including the high-fidelity, highresolution Monte-Carlo code Shift and the Computational fluid dynamics code NekRS. NekRS is a highly-performant open-source code for simulation of incompressible and low-Mach fluid flow, heat transfer, and combustion with a particular focus on turbulent flows in complex domains. It is based on rapidly convergent high-order spectral element discretizations that feature minimal numerical dissipation and dispersion. State-of-the-art multilevel preconditioners, efficient high-order time-splitting methods, and runtime-adaptive communication strategies are built on a fast OCCA-based kernel library, libParanumal, to provide scalability and portability across the spectrum of current and future high-performance computing platforms. On Frontier, Nek5000/RS has recently achieved an unprecedented milestone in breaching over 1 billion spectral elements and 350 billion degrees of freedom. Shift has demonstrated the capability to transport upwards of 1 billion particles per second in full core nuclear reactor simulations featuring complete temperature-dependent, continuous-energy physics on Frontier. Shift achieved a weak-scaling efficiency of 97.8% on 8192 nodes of Frontier and calculated 6 reactions in 214,896 fuel pin regions below 1% statistical error yielding first-of-a-kind resolution for a Monte Carlo transport application.

Hamilton, Steven P.↗
DOE OSTI2023

Temperature-based reactive flow model for TATB plastic bonded explosives

A new reactive flow model is presented for TATB-based plastic bonded explosives, applicable to shock initiation and detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and detonation products. The model incorporates temperature-dependent decomposition reaction rates. Arrhenius model parameters were derived from subscale quantum-based molecular dynamics simulations of shock-induced TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the late-time slow release of energy inherent to detonation in TATB-based materials. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments on PBX 9502. Here, the model is found to perform equally well in predicting the size effect curve of ambient, cold, and hot PBX 9502 rate sticks.

Chemistry - Chemical explosives↗
DOE OSTI2024

Embedding Neural Thermal Scattering (NeTS) Modules in SERPENT for Higher Fidelity Advanced Reactor Analysis

When a neutron born in fission thermalizes to the order of $k$ $B$ $T$, it’s de-Broglie wavelength and energy approach the order of inter-atomic spacing and elementary lattice oscillations, respectively. $S$($a,β,t$) or the scattering law, uuantify these temperature-dependent crystallographic contributions to total cross section (or reaction rate). In a Monte Carlo analysis, cumulative distribution functions (CDFs) of $S$($a,β,t$) are loaded to memory from “A Compact ENDF” (ACE) files for stochastically selecting thermal scattered neutron trajectories. In this work, novel neural thermal scattering (NeTS) modules for $S$($a,β,t$) CDFs are designed, trained, serialized and embedded within SERPENT using Python’s limited C-API for on-the-fly deployment of crystalline graphite $S$($a,β,t$) sampling. Torchscript tracing and Numba just-in-time (JIT) compilation streamline neural inference on NVIDIA GPUs with CUDA libraries. Demonstrations of bare sphere thermalization of fast and thermal sources show excellent agreement between embedded NeTS in SERPENT and MCNP. With an explicit model of the reactor, NeTS can predict on-the-fly changes in TREAT neutron spectra as a function of local temperature, which can serve to improve transient and accident predictions in a multiphysics analysis framework. This framework can be further extended to account on-the-fly for changes in local graphitic microstructure to scattering cross sections, and outlines a novel coupling of modern machine learning with state-of-the-art reactor physics methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen

Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.

first-principles calculations↗
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 OSTI2025

Developing a Nuclear Quality Assurance Compliant Design Methodology for Neutronic Analysis of Xe-100 Design

The primary objective of this work is to develop a design methodology compliant with nuclear quality assurance standards for the Xe-100 neutronic design verification studies. To achieve this, a Monte Carlo model of the Xe-100 reactor was constructed using the exclusion principle, transformation technique, and universe-based level specification following Idaho National Laboratory (INL) NQA level-1 compliant standards and an NQA-1 compliant version of MCNP6. The model encompasses the entire reactor core structures, including the upper plenum, core region, and lower plenum sections, along with all sub-components. The active core section was represented using the spectral regions, each comprising a particular fuel composition and temperature averaged over the considered zone, calculated by X-energy using Very Superior Old Programs (VSOP). Additionally, a component-wise temperature map was implemented into the model, not only for the core region but also for the structural components. Temperature-dependent cross-section libraries, along with thermal scattering law libraries, generated using INL NQA-1 compliant version of NJOY21, were utilized for each isotope in the burnt fuel and the structural materials. Furthermore, the volume of each modeled component was estimated using a stochastic approach with the ray tracing method in MCNP and criticality calculations were performed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗