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 811 records · Page 45

Follow-On Status Report on FY2022 Model Development within the Integrated Energy Systems Hybrid Repository

This publication details newly created energy storage and reactor models developed within the HYBRID modeling repository as part of the Department of Energy Office of Nuclear Energy (DOE-NE) Integrated Energy Systems (IES) program, led by Idaho National Laboratory (INL). Model development to-date includes creation of dynamic systems-level models of a pebble bed high temperature gas reactor (HTGR), compressed air energy storage (CAES), and Modelica standard library based two-tank sensible heat storage (SHS) in the IES-based HYBRID repository. Models are developed using the latest publicly available data and incorporate the possibility of control strategy inclusion for use with the existing IES modeling, analysis, and optimization toolset. Simulations showcase the abilities of each technology to flexibly operate in ways consistent with IES operation expectations. When these models are available, they can be utilized within different integrated energy park concepts to understand optimal system operation, control, and dispatching. Moreover, given the generic nature of the models, industrial partner technologies can be quickly added to the repository using the existing models as a basis. Additional dynamic models for thermal energy storage concepts can be developed and added to the HYBRID repository as needed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Dust Formation and Destruction in Core Collapse Supernovae Explosions [Dissertation]

Dust is a crucial component of the Interstellar Medium (ISM) and an important driver in many astrophysical processes. Dust formed from evolved stars increases the metallicity of the ISM, enriching it with heavier elements. Dust makes up the source material of collapsing clouds in the ISM leading to stellar and galactic formation. The stabilization of these systems are aided by dust through absorption and scattering of excess energy. As stellar systems evolve and undergo supernova and mass loss events, enriched dust is produced, further enhancing the metallicity of the ISM. However, despite the importance of dust, the processes that produce dust are not well understood. Core-Collapse Supernovae are one source of dust. Dying massive stars explode, sending their constituent material into the Interstellar Medium, which as the ejecta expands and cools, material condenses and nucleates forming dust grains. I investigate the formation, survivability, and properties of dust. In Chapter 2, dust grain nucleation in CCSNe is discussed. A revised formulation of Kinetic Nucleation Theory (KNT) is used to track the nucleation and growth of dust grains in the ejecta. The affects of the progenitor system on the amount and type of dust is explored, showing a strong dependence on explosion energy. I discuss the code and methods used to model dust nucleation and destruction in Chapter 3. The sputtering and erosion of dust grains in CCSNe shocks are explored in Chapter 4. The erosion and survival of this nucleated dust is modeled through thermal and non-thermal sputtering by the surrounding gas and shocks. The optical properties and scattering of light by surviving dust grains is modeled by Mie scattering in Chapter 5. A discussion of observations of dust in variable star systems and supernovae is outlined in Chapter 6.

79 ASTRONOMY AND ASTROPHYSICS↗

A Hybrid Biophysical‐Machine Learning Framework for Diurnal Surface Energy Flux Estimation Using Proximal Sensing

Thermal infrared-based remote sensing of surface energy fluxes has traditionally relied on high spatial resolution satellite data with revisit frequencies on the order of weeks. In this study, we evaluate a biophysics-based analytical surface energy balance model for predicting latent energy ( LE ) and sensible heat ( H ) fluxes using proximal sensing observations. The Surface Temperature Initiated Closure (STIC1.2) model has been extensively validated across a wide range of spatial and temporal scales using various satellite-derived thermal infrared data sets. Here we extend this validation by applying STIC at sub-hourly temporal resolution over multiple growing seasons for four distinct agricultural systems. We further develop and evaluate novel STIC variants that incorporate machine learning (ML) techniques to eliminate the need for surface energy balance observations, specifically net radiation and soil heat flux, thereby enhancing model applicability in data-sparse settings. The integration of a ML component to estimate surface available energy is shown to have strong predictive performance for both LE (R 2 = 0.81–0.94) and H (R 2 = 0.46–0.72) across all agricultural systems examined here, demonstrating the potential of hybrid biophysical-machine learning approaches for surface energy balance modeling with minimal data requirements. This study concludes with a novel application of explainable machine learning (exML) to diagnose sources of model error. This exML framework attributes residual prediction errors to both model input variables and environmental drivers not explicitly included in the simulation experiments. This approach provides a new pathway for improving model design and integrating previously overlooked yet influential variables into future model iterations.

evapotranspiration↗

Design and Optimization of a Gas-Cooled, Airfoil Fin Microchannel Heat Exchanger

High-performance microchannel heat exchangers are needed to supply heat for power conversion for nuclear microreactors. An airfoil fin microchannel design, constructed of Alloy 617 with helium as the working fluid, was analyzed and optimized using a design of experiments with artificial intelligence and machine learning techniques. The use of airfoil fins offers the potential to reduce pressure drop across the heat exchanger, as compared to other types of channel configurations. A framework for topology optimization of airfoil fin PCHEs has been developed that can be readily extended to different fin sizes and shapes, as well as different inlet and operating conditions, materials of construction, and working fluids. An optimization procedure was developed that employs computational fluid dynamics for a set of design points identified using Latin hypercube sampling. STAR-CCM+ was used to analyze a simplified two-channel configuration where five parameters were varied – inlet angle, fin scale, extent of staggering, transverse and longitudinal pitches. Two methods were compared for generating surrogate models – a 5D polynomial and a regression neural network. A response surface approximation was created from the surrogate models and input to a genetic algorithm. The genetic algorithm identified a set of optimal points on the Pareto front. The optimal geometry was found across six channel Reynolds numbers ranging from 1000 to 5000 to analyze how varying inlet conditions affects the optimal design. A set of optimal designs that maximizes heat transfer and minimizes pressure drop was identified, and a thermal stress analysis was performed on the optimal design. This work has developed a digital framework for the expedient topology design and evaluation of PCHE designs for gas-cooled microreactor applications. Correlations for the Nusselt number and Darcy friction factor were developed that can be useful for thermal hydraulic analyses using system codes. A thermal stress analysis was conducted and a brief discussion of the status of code cases of PCHEs for nuclear applications is given. Testing and thermomechanical modeling is needed to facilitate future code compliance of PCHEs for high pressure and high temperature applications.

42 ENGINEERING↗

Steady-state off-design modeling of the supercritical carbon dioxide recompression cycle for concentrating solar power applications with two-tank sensible-heat storage

Concentrating solar power researchers are evaluating the potential of the supercritical carbon dioxide recompression cycle to improve the thermal efficiency and decrease the capital costs of next-generation systems. This analysis investigates the steady-state off-design performance of a recompression cycle integrated with two-tank sensible-heat thermal energy storage as the ambient temperature and heat-transfer fluid (HTF) inlet conditions to the cycle change. This paper presents off-design component models and then cycle convergence and control models to maximize net cycle power output while constraining the high-side pressure and air-cooler fan power to their respective design values and fixing the cycle HTF outlet temperature to its design value. Results show that inventory control and air-cooler fan power are important control parameters that can be optimized to maximize off-design power output as the ambient temperature and HTF mass flow rate diverge from design. The high-side pressure and fan power constraints cause the cycle net power to degrade when the ambient temperature is warmer than the design value, and this study calculates a maximum mass flow rate for hot days above which the cycle cannot achieve the design HTF outlet temperature. Finally, the analysis shows that optimizing compressor shaft speeds can improve cold-day performance by around 1.5 percentage points and part-load performance by up to 2.5 percentage points versus a baseline case with no active compressor control.

14 SOLAR ENERGY↗

The MOOSE Thermal Hydraulics Module

The Multiphysics Object-Oriented Simulation Environment (MOOSE) is an open-source object-oriented finite element framework written in C++ (Lindsay et al., 2022). The Thermal Hydraulics Module (THM) is an optional MOOSE physics module that provides capabilities for studying thermal hydraulic systems. Its core capability lies in assembling a network of coupled components, for instance, pipes, junctions, and valves. THM provides several new systems to MOOSE to enable and facilitate thermal hydraulic simulations, most notably the Components system, which provides a higher-level syntax to MOOSE’s lower-level objects. This system is extensible by the user, but the current library primarily includes components based on a one-dimensional, single-phase, variable-area, compressible flow model, as well as heat conduction.

42 ENGINEERING↗

Deep reinforcement learning control for co-optimizing energy consumption, thermal comfort, and indoor air quality in an office building

With the recent demand for decarbonization and energy efficiency, advanced HVAC control using Deep Reinforcement Learning (DRL) becomes a promising solution. Due to its flexible structures, DRL has been successful in energy reduction for many HVAC systems. However, only a few researches applied DRL agents to manage the entire central HVAC system and control multiple components in both the water loop and the air loop, owing to its complex system structures. Moreover, those researches have not extended their applications by incorporating the indoor air quality, especially both CO2 and PM2.5concentrations, on top of energy saving and thermal comfort, as achieving those objectives simultaneously can cause multiple control conflicts. What's more, DRL agents are usually trained on the simulation environment before deployment, so another challenge is to develop an accurate but relatively simple simulator. Therefore, we propose a DRL algorithm for a central HVAC system to co-optimize energy consumption, thermal comfort, indoor CO2 level, and indoor PM2.5 level in an office building. To train the controller, we also developed a hybrid simulator that decoupled the complex system into multiple simulation models, which are calibrated separately using laboratory test data. The hybrid simulator combined the dynamics of the HVAC system, the building envelope, as well as moisture, CO2, and particulate matter transfer. Three control algorithms (rule-based, MPC, and DRL) are developed, and their performances are evaluated on the hybrid simulator environment with a realistic scenario (i.e., with stochastic noises). The test results showed that, the DRL controller can save 21.4 % of energy compared to a rule-based controller, and has improved thermal comfort, reduced indoor CO2 concentration. The MPC controller showed an 18.6 % energy saving compared to the DRL controller, mainly due to savings from comfort and indoor air quality boundary violations caused by unmeasured disturbances, and it also highlights computational challenges in real-time control due to non-linear optimization. Finally, we provide the practical considerations for designing and implementing the DRL and MPC controllers based on their respective pros and cons.

Guo, Fangzhou↗

A conservative, interface-resolved, compressible framework for the modeling and simulation of liquid/gas phase change

This paper presents a method for simulating evaporation in a compressible, interface-resolved framework appropriate for modeling problems of engineering interest. In order to achieve robustness and broad applicability, the method has been designed to discretely enforce consistent mass and thermal energy transport at the phase interface, to globally conserve mass, momentum, and energy, and to be capable of modeling compressible and incompressible systems. Verification is performed via the Sod-shock test, one-dimensional heat conduction, evaporation from a planar interface, and evaporation of three-dimensional droplets. Convergence with increasing mesh resolution is demonstrated in all tested configurations, and conservation is maintained near machine precision for a translating droplet. Conservation and accurate phase change rates are preserved at the low numerical resolutions commonly encountered in engineering calculations. Following verification, the method is validated by comparison to an empirical correlation for evaporating droplets in high temperature crossflow, and the presentation concludes with the simulation of an iso-octane spray at conditions representative of gasoline direct injection. In conclusion, successful verification, validation, and demonstrated practical utility suggest the method to be an accurate, efficient, and robust approach for the study of phase change in engineering systems.

97 MATHEMATICS AND COMPUTING↗

Hot-tail electrons’ impact on assimilation and injection penetration of D2 shattered pellet injections

The assimilation and the penetration of injected materials play significant roles in the mitigation efficiency of shattered pellet injection (SPI) as a disruption mitigation system (DMS), hence an accurate description of the fragment ablation rate is desirable. Current mainstream 3D MHD codes modelling SPIs mostly assume instantaneous thermalization between the previously hot ambient electrons and the newly released cold electrons. This could result in underestimation of the ablation rate if the hot electron thermalization time is comparable or even longer than the fragment flying time across flux tubes, putting the validity of these 3D SPI simulations into question. To resolve this doubt, we, in this study, investigate the thermalization dynamics and the overall hot-electron impact. The finite-time collisional thermalization of hot-tail electrons in a rapidly cooling plasma, as well as the so-called ‘self-limiting’ effect are considered. The former effect tends to deplete the colder population within a hot-tail species, while the latter is found to preferentially deplete the higher energy population. The combined result is found to cause an almost self-similar decay of the hot electron distribution function, while its shape does not deviate much from that of Maxwellian distribution and the mean energy does not change much during the thermalization process. Based on this observation, axisymmetric JOREK D2 SPI simulations were carried out with additional hot-tail contribution to evaluate their overall impact onto the injection assimilation and penetration. It is found that the hot-tail effect indeed causes enhanced assimilation and shallower penetration, although the overall effect depends on the exact injection configuration, with the slow injection showing negligible hot-tail effect while the fast single non-shattered pellet case shows drastic hot-tail ablation enhancement. For ITER-like SPI parameters, there is no significant deviation in the total assimilation, but some deviation in the injection penetration is observed for the fast injection velocity cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Confinement of many-body Bethe strings

Based on the Bethe ansatz approach and inelastic neutron scattering experiments, we reveal the evolution of confinement of many-body Bethe strings in ordered regions of the quasi-one-dimensional antiferromagnet YbAlO 3 . In the antiferromagnetic phase, the spin dynamics is dominated by confined length-1 Bethe strings, whose dominancy in the high-energy branch of the excitation spectrum yields to confined length-2 Bethe strings when the material is tuned to the spin-density-wave phase. In the thermal-induced disordered region, the confinement effect disappears, and the system restores the conventional quantum integrable physics of the one-dimensional Heisenberg model. Finally, our results establish a unified picture based on a Bethe string for the spin dynamics in different magnetic phases of YbAlO 3 , and thus provide profound insight into many-body quantum magnetism.

1-dimensional spin chains↗

Reinforcement Learning Control for Buildings Co-Optimizing Energy, Comfort, and Indoor Air Quality: An Annual Assessment

Efficient control of Heating, Ventilation, and Air Conditioning (HVAC) systems is crucial for optimizing energy use and maintaining indoor comfort in buildings. Traditional control methods, such as PID control, cannot handle energy use trade-offs among multiple components in the building energy system at a supervisory level. Reinforcement learning (RL) presents a promising solution, offering adaptive and data-driven control strategies that optimize performance over time. However, RL also faces several challenges, including the conflicts encountered in co-optimizing energy savings, occupant comfort, and indoor air quality, and the requirement for extensive interactions with the environment in training. We proposed a flexible simulation platform that integrates a hybrid model for RL training and designed an RL agent to control the entire central HVAC system, focusing on co-optimizing energy consumption, thermal comfort, and indoor air quality ($\text{CO}_{2}$ and PM2.5 concentrations). Finally, we evaluated the RL agent's performance over an annual cycle. Our findings indicate that the RL agent can effectively manage the HVAC system with 14.7 % energy savings annually and balance multiple objectives, which demonstrates significant potential for improving HVAC system control and sustainability in buildings.

Guo, Fangzhou↗

WHOLESCALE: Coordinates of wells at San Emidio, Nevada

This dataset includes position coordinates and elevation information for wells at the WHOLESCALE San Emidio project location. Well positions in the attached file are characterized by UTM coordinates (Easting, Northing) in meters, and WHOLESCALE coordinates (Easting, Northing) relative in meters from a chosen reference point. The elevation of the top and bottom of open intervals within each well (representing locations of perforated or open-hole sections) are measured in meters (positive) above WGS84 geoid (mean sea level). The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States.

15 GEOTHERMAL ENERGY↗

Modeling Approach for the Aluminum-clad Dry Storage Pilot using HFIR Fuel

To confirm that the dry storage of aluminum-clad research reactor spent nuclear fuel (ASNF) will remain within the safety envelope after applied drying schemes and that the resulting evolution of the gas space composition, temperature, and pressure conditions are understood, a dry storage pilot project is being established. The pilot will incorporate an instrumented lid for discrete interval or for on-demand gas composition and temperature monitoring of two DOE Standard Canisters (DSCs) loaded with three High Flux Isotope Reactor (HFIR) inner cores per DSC. Each DSC would be subjected to a separate alternative candidate drying scheme. Canisters will undergo 1 to 5 years of monitoring, including internal temperature and gas sampling to track pressure and composition changes. This report outlines the approach for modeling the ASNF-in-canister behavior in terms of evolving gas space conditions for the ASNF dry storage pilot using HFIR fuel. The ASNF has an adherent surface oxyhydroxide layer comprised of boehmite/bayerite that generates hydrogen when subjected to irradiation. Three-dimensional multi-physics computational fluid dynamics simulations will be executed to compute the thermal field within the DSC and provide inputs to a chemical model employed to compute pressure buildup as hydrogen is generated in the system. Implemented in Cantera, the chemical model solves gas phase and aluminum oxyhydroxide surface-mediated radiolysis reactions. Gas phase reactions are sourced from Wittman and Hanson (2015), whereas surface-mediated reactions are incorporated by fitting experimental data using an optimization algorithm (Abboud, 2023). Water radiolysis reactions from Wren and Ball (2001) are adopted with modifications as described in Abboud (2023c). Understanding the effect of the hydrogen buildup over time is important for long-term storage safety considerations. Modeling results will include the canister pressure, temperature, and composition evolution from the initial helium backfill with the addition of radiolytically-evolved chemical species (e.g., hydrogen and oxygen). The specific HFIR cores for the pilot program have not yet been selected, and the overall design is still in development. The CFD-chemical model used for this work will be based on prior models with necessary updates to allow for improved accuracy and efficiency. The experimental data obtained from the HFIR demonstration will be used to improve and validate the computational models to predict the ASNF-in-canister behavior.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Simulated wind speed and initial conditions over the WFIP2 region: Cold-front (D01)

The purpose of this work is to assess the sensitivity of the forecast for turbine height wind speed to initial condition (IC) uncertainties over the Columbia River Gorge and Columbia River Basin for two typical weather phenomena: a local thermal gradient induced by a marine air intrusion and passage of a cold front. The Weather Research and Forecasting (WRF) model data assimilation system (WRFDA) was used to generate ensemble ICs from the North American Regional Analysis (NARR) for the WRF model initialization. The simulated turbine-height wind speeds were categorized into four types using the self-organizing map (SOM) technique. This work advances understanding of IC uncertainties impacts on wind speed forecasts and locates the high-impact regions.

17 WIND ENERGY↗

Simulated wind speed and initial conditions over the WFIP2 region: Cold-front (D02)

The purpose of this work is to assess the sensitivity of the forecast for turbine height wind speed to initial condition (IC) uncertainties over the Columbia River Gorge and Columbia River Basin for two typical weather phenomena: a local thermal gradient induced by a marine air intrusion and passage of a cold front. The Weather Research and Forecasting (WRF) model data assimilation system (WRFDA) was used to generate ensemble ICs from the North American Regional Analysis (NARR) for the WRF model initialization. The simulated turbine-height wind speeds were categorized into four types using the self-organizing map (SOM) technique. This work advances understanding of IC uncertainties impacts on wind speed forecasts and locates the high-impact regions.

17 WIND ENERGY↗

Simulated wind speed and initial conditions over the WFIP2 region: Sea-breeze (D01)

The purpose of this work is to assess the sensitivity of the forecast for turbine height wind speed to initial condition (IC) uncertainties over the Columbia River Gorge and Columbia River Basin for two typical weather phenomena: a local thermal gradient induced by a marine air intrusion and passage of a cold front. The Weather Research and Forecasting (WRF) model data assimilation system (WRFDA) was used to generate ensemble ICs from the North American Regional Analysis (NARR) for the WRF model initialization. The simulated turbine-height wind speeds were categorized into four types using the self-organizing map (SOM) technique. This work advances understanding of IC uncertainties impacts on wind speed forecasts and locates the high-impact regions.

17 WIND ENERGY↗

Simulated wind speed and initial conditions over the WFIP2 region: Sea-breeze (D02)

The purpose of this work is to assess the sensitivity of the forecast for turbine height wind speed to initial condition (IC) uncertainties over the Columbia River Gorge and Columbia River Basin for two typical weather phenomena: a local thermal gradient induced by a marine air intrusion and passage of a cold front. The Weather Research and Forecasting (WRF) model data assimilation system (WRFDA) was used to generate ensemble ICs from the North American Regional Analysis (NARR) for the WRF model initialization. The simulated turbine-height wind speeds were categorized into four types using the self-organizing map (SOM) technique. This work advances understanding of IC uncertainties impacts on wind speed forecasts and locates the high-impact regions.

17 WIND ENERGY↗

High quality superconducting Nb co-planar resonators on sapphire substrate

We present measurements and simulations of superconducting Nb co-planar waveguide resonators on sapphire substrate down to millikelvin temperature range with different readout powers. In the high temperature regime, we demonstrate that the Nb film residual surface resistance is comparable to that observed in the ultra-high quality, bulk Nb 3D superconducting radio frequency cavities while the resonator quality is dominated by the BCS thermally excited quasiparticles. At low temperature both the resonator quality factor and frequency can be well explained using the two-level system models. Through the energy participation ratio simulations, we find that the two-level system loss tangent is $\sim 10^{-2}$, which agrees quite well with similar studies performed on the Nb 3D cavities.

Zhu, S.↗