Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal system modeling”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Bentonite buffer under high temperature: laboratory experiments and coupled process modeling

Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.

Zheng, Liange↗
DOE OSTI2025

CFD modeling of natural circulation in LiCl-KCl molten salt closed loop

Characterizing flow within a molten salt closed-loop system is crucial for assessing system requirements, evaluating performance, and identifying potential flaws. Direct flow measurement using instrumentation is challenging due to extreme environmental conditions and the limitations associated with measuring molten salt flow under natural convection. Here, this study aims to provide comprehensive insights into the thermal-hydraulic behavior of a closed loop, with a particular focus on temperature distribution and velocity prediction. The Computational Fluid Dynamics (CFD) model demonstrated the capability to effectively simulate and predict both temperature distributions and flow velocities within the molten salt loop. The CFD model's predictive capability was validated by its ability to replicate temperature measurements under varying boundary conditions. The analysis revealed that the CFD model tends to underpredict temperatures in the cold leg and overpredict them in the hot leg, highlighting the need for continuous model refinement and acknowledging the limitations of using a steady-state approach. Furthermore, the potential of using external temperature measurements to estimate internal molten salt temperatures and predict flow velocity was explored, revealing that this approach could introduce up to a 5.5% error in flow velocity calculations. Line probes mapping temperature distributions across the tube's cross-section and molten salt provided valuable insights into temperature gradients, emphasizing the need for a thermal conductivity equation for molten salt with lower uncertainty to achieve more accurate temperature predictions of the system.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

First principles free energy model with dynamic magnetism for δ -plutonium

We present an ab initio free energy model derived from a fully relativistic density functional theory (DFT) electronic structure with dynamic magnetism for δ -plutonium (face-centered cubic, fcc). The DFT model is extended with orbital-orbital interaction in a parameter free orbital polarization (OP) mechanism consistent with previous modeling of plutonium. Gibbs free energy is built from components associated with the temperature dependence of the electronic structure and the corresponding electronic entropy, lattice vibrations within an anharmonic lattice dynamics model, and dynamical fluctuations of the magnetization density, i.e. magnetic fluctuations. The fluctuation model consists of transverse and longitudinal modes driven by temperature induced excitations of the DFT + OP electronic structure. The ab initio model thus incorporates fluctuating states beyond the electronic ground state. Thanks to the dynamic magnetism, the theory predicts excellent thermodynamic properties and a Gibbs free energy in accord with CALPHAD and semi-empirical modeling developed from the thermodynamic observables. The magnetic fluctuations further explain anomalous behaviors of the thermal expansion in plutonium. Specifically, a thermal expansion for the δ -plutonium system turning from positive to negative at temperatures above room temperature, a tendency for gallium to reduce and remove the negative thermal expansion depending on composition, and a positive thermal expansion for the high temperature ϵ phase.

dynamic↗
DOE OSTI2024

A nonlinear journey from structural phase transitions to quantum annealing

Motivated by an exact mapping between equilibrium properties of a one-dimensional chain of quantum Ising spins in a transverse field (the transverse field Ising (TFI) model) and a two-dimensional classical array of particles in double-well potentials (the “$\phi$ 4 model”) with weak inter-chain coupling, we explore connections between the driven variants of the two systems. Here, we argue that coupling between the fundamental topological solitary waves in the form of kinks between neighboring chains in the classical $\phi$ 4 system is the analog of the competing effect of the transverse field on spin flips in the quantum TFI model. As an example application, we mimic simplified measurement protocols in a closed quantum model system by studying the classical $\phi$ 4 model subjected to periodic perturbations. This reveals memory/loss of memory and coherence/decoherence regimes, whose quantum analogs are essential in annealing phenomena. In particular, we examine regimes where the topological excitations control the thermal equilibration following perturbations. This paves the way for further explorations of the analogy between lower-dimensional linear quantum and higher-dimensional classical nonlinear systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Open-cycle thermochemical energy storage for building space heating: Practical system configurations and effective energy density

Salt-hydrate thermochemical materials (TCM) are promising candidates for energy storage systems for building space heating due to their high theoretical energy density and the need for low regeneration temperature. However, water vapor is required to drive the hydration process of the TCM reactor, which poses a challenge during winter when water vapor is typically scarce. Using indoor air directly lowers the building's humidity to an unconformable level in practice, while the cold outdoor air contains limited moisture. Here we consider different integration strategies for open-cycle TCM reactors in buildings and develop a model to simulate their thermal performance across diverse buildings and climates, specifically for building space heating. The potential energy densities and the levelized cost of storage of the TCM reactor are evaluated in practical scenarios to demonstrate the load-shifting potential of TCM systems for heating applications. We use a strontium chloride (SrCl 2 )-based composite as the baseline and explore the impact of various reactor and material changes to the energy density and levelized cost of storage.

25 ENERGY STORAGE↗
DOE OSTI2021

Implementation of an Experimental Setup Based on the Three-Omega Method for Thermal Conductivity Measurements of Molten Salts

Molten salt reactors (MSRs) are one of the most promising next-generation nuclear reactor designs which use molten salts as the fuel and coolant. To design, model, predict, license, and operate MSRs, the behavior of the fuel and coolant must be understood as a function of temperature and irradiation, i.e., burnup. Of particular importance and complexity is the reliable and reproducible measurement of thermal conductivity of molten salts at high temperatures (> 500 °C). Thermal conductivity provides information on the rate at which heat passes through a specified material and is dependent on chemical composition. Currently, in the national laboratory complex, several methods for measuring thermal conductivity on molten salt systems are under development, however, to date, consistent data have not been produced. Typically, thermal conductivity measurements of high-temperature melts have been conducted by using the Transient Hot Wire (THW) method where thermal conductivity is determined by the change in temperature over time at a certain distance from a heating wire. Another time-dependent technique is the laser flash method, commonly used for solid samples. Laser flash requires prior measurements of density and heat capacity to calculate the thermal conductivity of the sample, introducing additional measurement errors. Steady state and time-dependent measurements such as THW and laser flash have errors from convection and ambient temperature fluctuations (heat radiation). In order to produce high fidelity thermal conductivity data, a measurement technique needs to be developed that will reduce the effect of convection at high temperatures by measuring in a very fast time window. In addition, to measure on actinide and irradiated salts the sample size should be small. Therefore, this project aims to develop a frequency-domain technique based on the three-omega (3?) method, previously proven to work on solids, nanofluids and low temperature solar salts (25 to 250 °C). If the 3?-device can be designed with materials that can withstand high temperatures (>500 °C) and can be configured to reduce errors due to convection, then it would prove to be a valuable technique for ensuring thermal conductivity on molten salt systems needed for MSR development. For this proposal, the 3? instrument will be designed, assembled, and tested on non-radioactive salts before being used on actinide bearing salts. In addition, a thermal model will be developed to determine the thermal conductivity of the sample from the measured thermal oscillation of the sensor. Because 3? is a direct method that requires a small volume of sample and works in a very fast time window, it will aid in suppressing the interference due to convection which is ideal when working with high temperature liquid samples.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Energy Model to Evaluate Thermal Energy Storage Integrated with Air Source Heat Pumps: Preprint

Full decarbonization in buildings requires the replacement of combustion appliances with electric ones, and air source heat pumps (ASHP) are a candidate alternative. However, technical limitations, such as the efficiency decrease when operating in cold weather, limit their adoption in the global heating market. Among several options to improve ASHP efficiency operating in colder climates, thermal energy storage (TES) has been considered, as it may provide heating when it is cold and shift ASHP operation to times when the weather is warmer. It may also take advantage of time of use electricity rates and support defrosting when necessary. The evaluation of ASHP-TES systems, however, is still limited because traditional metrics do not capture their full economic and environmental benefits. In this work, a python framework is presented to model the ASHP with and without the presence of TES. Metrics are proposed to analyze the system performance in terms of costs, equivalent CO2 emission, and efficiency metrics to evaluate and compare alternative systems. Model validation against experimental data obtained for a commercial heat pump is provided, as well as an application example using Denver, Colorado, to highlight the model capabilities.

cold climate↗
DOE OSTI2020

Modeling study of deep direct use geothermal on the West Virginia university campus-morgantown, WV

To reduce the geothermal exploration risk, a feasibility study is performed for a deep direct-use (DDU) system proposed at the West Virginia University (WVU) Morgantown campus. This study applies numerical simulations to investigate reservoir impedance and thermal production. Because of the great depth of the geothermal reservoir, few data are available to characterize reservoir features and properties. Consequently, the study focuses on the following three aspects: 1. model choice for predicting reservoir impedance and thermal breakthrough: after investigating three potential models (one single permeability model and two dual permeability models) for flow through fractured rock, it is decided only the single permeability model is needed; 2. well placement (horizontal vs. vertical) options: horizontal well placement seems to be more robust to heterogeneity and the impedance is more acceptable; 3. Prediction uncertainty: the most influential parameters are identified using a First-Order-Second-Moment uncertainty propagation analysis, and the uncertain range of the model predictions is estimated by performing a Monte Carlo simulation. Heterogeneity has a large impact on the prediction, therefore, heterogeneity is included in the predictive model and uncertainty analysis. The numerical model results and uncertainty analysis will be used for further economic studies.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Pumped Thermal Energy Storage Using Low-Cost Particles and a Fluid Bed Heat Exchanger for Maximum Power Efficiency

Robust, efficient, cost-effective long-duration electricity storage (LDES) solutions can enhance grid resiliency, support existing transmission and distribution infrastructure, and enable a greater share of low-cost, variable alternative energy sources to penetrate the market. To meet this need, the project team at the National Renewable Energy Laboratory developed a transformative LDES system based on pumped thermal energy storage (TES) using low-cost particles and a fluid bed heat exchanger for maximum power efficiency (PUMP). The PUMP system is composed of high-temperature, low-cost particle TES coupled with an advanced pressurized fluid bed heat exchanger (PFB HX) that supports a high-efficiency pumped thermal energy storage (PTES) system integrated with concentrating solar thermal power (CSP). The PUMP project developed and de-risked a PFB HX and particle CSP system intended to be integrated with reversible turbomachinery and a modeling tool to assess PTES cost and performance.

14 SOLAR ENERGY↗
DOE OSTI2021

Comprehensive framework for data-driven model form discovery of the closure laws in thermal-hydraulics codes

The two-phase two-fluid model is a basis of many thermal-hydraulics codes used in design, licensing, and safety considerations of nuclear power plants. Thermal-hydraulics codes rely on the closure laws to close the system of conservation equations and describe the interactions between phases. These laws, derived from years of experimental investigations, are semi-empirical correlations that lack generality and have a limited range of applicability. Increase of computational power, availability of new experiments, and development of high-fidelity simulations has increased the number of validation data. The discrepancies between the code predictions and the validation data are a great source of knowledge. Missing physics that are not included in the model but are important for the considered phenomena can be discovered by propagating the information from the experimental results through the model. Furthermore, physics-discovered data-driven model form (P3DM) methodology integrates available integral effect tests and separate effects tests to determine the necessary corrections to the model form of the closure laws. In contrast to existing calibration techniques, the methodology modifies the functional form of the closure laws. Based on the functional form of the correction, the missing physics that were not included in the original model can be discovered. The methodology provides the alternative to the machine learning approach, in which the model is discovered in the form of the intractable black-box relation. In this work, the methodology was applied to the CTF subchannel code to improve the prediction of the two-phase flow phenomena.

42 ENGINEERING↗
DOE OSTI2023

A high‐order discontinuous Galerkin approach for physics‐based thermospheric modeling

Abstract The accurate prediction of aerodynamic drag on satellites orbiting in the upper atmosphere is critical to the operational success of modern space technologies, such as satellite‐based communication or navigation systems, which have become increasingly popular in the last few years due to the deployment of constellations of satellites in low‐Earth orbit. As a result, physics‐based models of the ionosphere and thermosphere have emerged as a necessary tool for the prediction of atmospheric outputs under highly variable space weather conditions. This paper proposes a high‐fidelity approach for physics‐based space weather modeling based on the solution of the Navier–Stokes equations using a high‐order discontinuous Galerkin method, combined with a matrix‐free strategy suitable for high‐performance computing on GPU architectures. The approach consists of a thermospheric model that describes a chemically frozen neutral atmosphere in nonhydrostatic equilibrium driven by the external excitation of the Sun. A novel set of variables is considered to treat the low densities present in the upper atmosphere and to accommodate the wide range of scales present in the problem. At the same time, and unlike most existing approaches, radial and angular directions are treated in a nonsegregated approach. The study presents a set of numerical examples that demonstrate the accuracy of the approximation and validate the current approach against observational data along a satellite orbit, including estimates of established empirical and physics‐based models of the ionosphere‐thermosphere system. Finally, a one‐dimensional radial derivation of the physics‐based model is presented and utilized for conducting a parametric study of the main thermal quantities under various solar conditions.

Engineering↗
DOE OSTI2023

The impact of non-ideal phase change properties on phase change thermal energy storage device performance

Phase change materials have been known to improve the performance of energy storage devices by shifting or reducing thermal/electrical loads. While an ideal phase change material is one that undergoes a sharp, reversible phase transition, real phase change materials do not exhibit this behavior and often have one or more non-idealities – glide, hysteresis, supercooling – associated with them. Experimental and modeling techniques to characterize these non-ideal properties are reasonably well understood, however, their impact on the performance of a thermal energy storage system is not fully understood. Herein, we analyze the performance of a heat exchanger with a phase change material as a function of the different non-idealities, for a heating and a cooling application. We focus on the impact of these non-ideal behaviors on different modes of operation of the thermal energy storage device, which will serve as a useful guide to researchers on the relative importance of the non-ideal properties, along with the necessary accuracy required during an experimental characterization.

25 ENERGY STORAGE↗
DOE OSTI2022

Vapor deposition rate modifies anisotropic glassy structure of an anthracene-based organic semiconductor

We control the anisotropic molecular packing of vapor-deposited glasses of ABH113, a deuterated anthracene derivative with promise for future organic light emitting diode materials, by changing the deposition rate and substrate temperature at which they are prepared. We find that at substrate temperatures from 0.65 T g to 0.92 T g , the deposition rate significantly modifies the orientational order in the vapor-deposited glasses as characterized by x-ray scattering and birefringence. Both measures of anisotropic order can be described by a single deposition rate–substrate temperature superposition (RTS). This supports the applicability of the surface equilibration mechanism and generalizes the RTS principle from previous model systems with liquid crystalline order to non-mesogenic organic semiconductors. We find that vapordeposited glasses of ABH113 have significantly enhanced density and thermal stability compared to their counterparts prepared by liquidcooling. For organic semiconductors, the results of this study provide an efficient guide for using the deposition rate to prepare stable glasses with controlled molecular packing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Energy Use in Quantum Data Centers: Scaling the Impact of Computer Architecture, Qubit Performance, Size, and Thermal Parameters

As quantum computers increase in size, the total energy used by a quantum data center, including the cooling, will become a greater concern. The cooling requirements of quantum computers, which operate at temperatures near absolute zero, are determined by computing system parameters, including the number and type of physical qubits, the packaging efficiency of the system, and the split between circuits operating at cryogenic temperatures and those operating at room temperature. When combined with thermal system parameters such as cooling efficiency and cryostat heat transfer, the total energy use can be determined using a first-principles energy model. These models show that cooling of quantum computers differs in two fundamental ways from conventional data centers: (1) the energy required for cooling is much greater than the energy required for computation, and (2) the cooling loads are sensitive to the computational architecture. The temperature requirements for different qubit types can change energy requirements by orders of magnitude. Power use and computational power, as quantified by quantum volume, are analytically correlated. Approaches are identified for minimizing energy use in integrated quantum systems relative to computational power. Furthermore, designing a sustainable quantum computer will require both efficient cooling and system design that minimizes cooling requirements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Performance Assessment of a Dual-purpose HP-TES for a Typical Year Comparing Different Climate Zones

Complete electrification of heating and cooling systems in the United States would overload the existing electrical grid. To this end, recent research efforts have focused on energy storage as one way to reduce ever-increasing peak electricity demands. One such technology is heat pump integrated thermal-energy storage (HP-TES) systems, which can enable load-shifting and reduce demand during peak hours. In this work, a Modelica model of a 4-Ton, dual-mode (heating and cooling) HP-TES system using a room-temperature phase-change material (PCM) integrated through a secondary hydronic loop was developed and simulated for a typical meteorological year in Riverside, CA (ASHRAE Climate Zone 2B), and Chicago, IL (ASHRAE Climate Zone 5A), for a DOE small office prototype building. The COP improvement and demand reduction were analyzed during the time-of-use peak hours along with the required recharge cycle energy demand. A typical summer in Riverside, CA resulted in model-predicted peak COP increases of 50 – 100% compared to the baseline heat pump, especially at higher outdoor temperatures, coupled with maximum peak energy reductions of 22%. During shoulder and winter season simulations, heating energy savings were higher as outdoor temperatures dropped below the PCM phase change temperature. For Chicago, IL, backup heating was necessary to meet the required system capacity for very cold ambient conditions, resulting in winter peak energy savings between 16 – 60% when accounting for backup heating requirements. The recharge energy was found to be mostly below the baseline HP energy during the peak. With improved component selections and controls, the recharge energy may be reduced. The outcomes from the simulations can provide first-order estimates for potential peak energy savings, challenges, and required controls in any location before conducting any field testing.

25 ENERGY STORAGE↗
DOE OSTI2022

Optimal operation and sizing of pumped thermal energy storage for net benefits maximization

Abstract Current trends in the modern grid are leading to the development and deployment of energy storage to help integrate increasing variable renewable energy sources into the grid. This paper studies a pumped thermal energy storage (PTES) system for multiple grid services including energy arbitrage, frequency regulation, spinning and non‐spinning reserve, and resource adequacy. Optimal dispatch methods are proposed for individual services as well as value stacking from multiple services to maximize the economic benefits. Assessment results demonstrate the superiority of value stacking. Specifically, the study shows the maximum revenue from an individual grid service with a 30‐MWh PTES system was $522,520, while the value stacking could increase the benefits to $678,477. In addition, sensitivity analyses were conducted to explore the cost‐effectiveness of a PTES system with different combinations of power transfer limits and energy capacity. It was found that the power transfer limit had a greater impact than the energy capacity on the benefits. The proposed method could help determine the optimal duration of a future PTES system.

25 ENERGY STORAGE↗
DOE OSTI2026

Integral Nuclear Data and Benchmarking Needs for Fusion Energy Systems

Fusion energy systems are currently being designed and optimized using radiation transport codes. To deal with the unique environment inside a fusion-based system, many of these designs incorporate novel materials able to withstand the high radiation fields, ensure adequate cooling and thermal protection, and produce tritium. Validation plays a vital role in building trust in the predictive power of these models and computational methods. Validation of a code consists of modeling documented real-world experiments and comparing the code-predicted response to the measured response. Adequate validation requires measured responses from real-world experiments, also known as integral data, that mimic the system being designed, including materials, impinging radiation, and temperature, among other variables. The most trusted integral data are experimental responses that have been through a rigorous benchmarking process that develops a recommended computational model and evaluates all experimental uncertainties. Finally, there are a few research groups around the world that have been producing integral data for fusion applications, but a substantial investment is needed to address the unique validation needs of the fusion community.

Fusion↗
DOE OSTI2026

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗