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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Arbitrage and Capacity Firming in Coordination with Day-Ahead Bidding of a Hybrid PV Plant

A hybrid PV plant (HPP) combines a photovoltaic (PV) plant with a battery energy storage system (BESS), which is considered a promising step towards the future of renewable power plants by the U.S. Department of Energy. When the renewable penetration reaches a significant level, a hybrid PV plant can bid in as a controllable thermal plant in the future electricity market. In this study, a bidding and BESS scheduling model is proposed for the HPP. The robust optimization (RO) technique has been utilized to identify the worst-case scenario of uncertainties during the bidding process. To address the overly conservative issue of the single-stage RO, we have decoupled the BESS schedule for arbitrage and PV capacity firming by a two-stage RO formulation. By comparing the output of single-stage RO and two-stage RO, the two-stage RO bids and schedules in a more aggressive manner, which increases the income of HPP. Also, the penalty of under-generation is considered in our model so that the day-ahead bidding decision and arbitrage schedules can be adjusted based on the potential UNDER-GENERATION penalty. Because the proposed model is non-convex and contains multi-stages, the Column-and-Constraint Generation (C&CG) algorithm is applied to the model as the solution. The proposed model has shown better economic performance compared to a state-of-art single-stage bidding method in case studies.

BESS scheduling↗
DOE OSTI2023

Techno-Economic Analysis and Market Potential of Geological Thermal Energy Storage (GeoTES) Charged With Solar Thermal and Heat Pumps into Depleted Oil/Gas Reservoirs and Shallow Reservoirs: A Technology Overview: Preprint

Depleted oil/gas reservoirs represent a waste of underground resources and investments of drilling, and also a potential risk to the earth's environment. Geologic thermal energy storage (GeoTES) is proposed as a solution to convert depleted oil/gas reservoirs into long-term seasonal energy storage. GeoTES can be hybridized with other techniques for viable commercial deployment, such as 1) concentrating solar power (CSP) collectors and 2) heat pumps with excess renewable energy. Here, a technology overview is given on the two GeoTES technologies, which includes system overview, techno-economic models, and case study. Both GeoTES technologies show great potential in the economics of individual project deployment and applicability across the US.

concentrating solar power↗
DOE OSTI2021

Modeling heat transport processes in enhanced geothermal systems: Validation study from EGS Collab Experiment 1

Heat recovery from enhanced geothermal systems (EGS) is a complex process involving heat transport in both fracture networks and rock formations. A comprehensive understanding of and the ability to model the underlying heat transport mechanisms is important for the success of EGS but remains challenging in practice due to the generally insufficient characterization of EGS reservoirs. In the present study, we analyze an extensively monitored intermediate-scale EGS field experiment performed in a well-characterized testbed. The high-resolution, high-quality measurements from the field experiment enable the development of a high-fidelity model incorporating a well-constrained fracture network. Based on the field experiment, we investigate the complex heat transport processes in an EGS-relevant environment and validate the capability of a numerical approach in simulating these inherently coupled heat transport processes. A series of numerical simulations were performed to study the effects of different heat transport mechanisms, including thermal convection with fracture flow, thermal conduction in rock formations, and the Joule-Thomson effect. The agreement of thermal responses between field measurements and simulation results indicates that our numerical approach can appropriately model the heat transport processes pertaining to heat recovery from EGS reservoirs.

Wu, Hui↗
DOE OSTI2025

Operando X-ray imaging reveals size-dependent evolution of cobalt oxide thermochemical material during thermal redox cycles

Multivalent metal oxides are promising thermochemical materials (TCMs) for energy storage and conversion owing to their high energy density, air compatibility, and high-temperature stability. Co 3 O 4 serves as a model system for examining particle-size- and structure-dependent redox behavior. While particle size and porosity are known to affect performance, their interplay and the kinetics of pore formation during cycling remain unclear. Here we show the chemical and 3D morphological evolution of Co 3 O 4 micro- and nanoparticles during redox cycles at 800–900 °C using thermal analysis, in-situ synchrotron transmission X-ray microscopy (TXM), and scanning electron microscopy. Thermal analysis shows that nanoparticles re-oxidize more rapidly than microparticles at 800 °C. In-situ nanotomography and chemical imaging reveals that nanoparticles undergo redox conversion without forming internal pores, whereas microparticles develop isolated porosity during reduction. These pores persist through re-oxidation, correlating to a lower conversion rate in subsequent cycles. Our results demonstrate distinct degradation kinetics in Co 3 O 4 micro- and nanoparticles, underscoring the critical role of particle size and porosity in redox performance and informing strategies to enhance the long-term efficiency of metal oxide TCMs.

25 ENERGY STORAGE↗
DOE OSTI2024

DFT Calculation of Thermal Expansion within Debye-Grüneisen Framework Made Simple and Accurate

The density functional theory (DFT) based Debye-Grüneisen approach is a computational method used to calculate the thermodynamic properties of materials by relating the volume dependence of the Debye frequency to the Grüneisen parameter. In this study, we propose a new model to describe the volume dependence of the Grüneisen parameter. This model offers a simple yet accurate approach to predict the thermal expansion of materials within the framework of the Debye-Grüneisen theory. Specifically, for pure elements or species used as a reference system, we calibrated the parameters to reproduce the experimental coefficient of thermal expansion, heat capacity, and entropy. For compounds with arbitrary compositions, the parameters are determined by taking the arithmetic or geometric average of those obtained from the pure elements or species based on their composition. We then compared the model prediction with experiments. This allowed for a comprehensive characterization of the thermal expansion behavior of various materials.

Wang, Yi↗
DOE OSTI2026

Organic Rankine Cycle Integration and Optimization for High Efficiency CHP Genset Systems (Final Technical Report)

This project successfully advanced the integration of Organic Rankine Cycle (ORC) technology with reciprocating engine–based combined heat and power (CHP) systems to improve electrical efficiency, total CHP efficiency, and grid-responsive operation. Over three budget periods, the work progressed from high-temperature ORC component development and thermodynamic model validation to next-generation system design, working fluid transition, and techno-economic analysis. Key technical accomplishments include development and validation of a thermodynamic model capable of accurately predicting ORC performance across an expanded temperature and pressure envelope; successful identification and validation of low-global-warming-potential (GWP) working fluids—most notably R1233zd(E)—as viable replacements for R245fa; and demonstration of scalable ORC architectures suitable for integration with 1–20 MW class reciprocating engines. These advances enable flexible CHP configurations that can increase electrical output while maintaining high overall utilization of available thermal energy. The project also produced a clean-sheet design for a next-generation ORC system targeting substantially higher power output per unit, supported by detailed component selection, heat exchanger evaluation, and system-level modeling. Techno-economic analyses indicate that ORC-enabled flexible CHP systems can meet or exceed Department of Energy (DOE) efficiency targets while providing value to both facility operators and the electric grid.. Late-stage testing of the largest next-generation ORC prototype identified limitations related to pump net positive suction head (NPSH) requirements and condenser flooding under certain operating conditions. Although these issues constrained full validation of that configuration within the project timeframe, they provided clear and actionable design guidance for future system refinements. Importantly, validated modeling, smaller-scale testing, and working fluid evaluations confirmed the technical viability of the overall approach. In aggregate, this project met its core objectives by establishing validated design tools, de-risking key ORC technologies for CHP applications, and defining a credible pathway toward commercialization of flexible, high-efficiency CHP systems. The results form a strong foundation for continued development and deployment beyond the conclusion of the DOE-funded effort.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

REDOTHERM: a thermodynamic modeling framework for redox-based thermochemical processes

Two-step thermochemical redox cycles are being developed as a potential pathway for the production of hydrogen and syngas. While there are many possible reactor and system configurations, moving oxide systems are considered promising in terms of the redox thermodynamics, due to the potential implementation of a countercurrent system that can achieve higher performance compared to other configurations. There is a lack of a robust thermodynamic modeling framework in the field, with multiple models incorporating incorrect thermodynamic assumptions that violate the second law of thermodynamics. We present in this work REDOTHERM, an open-source system model for moving oxides that incorporates the correct thermodynamic limits, as well as various options for the system auxiliary units including product separation, heat recovery, and oxygen removal. The model is agnostic to the energy source, and could be used for solar thermal or other configurations. We highlight the uses of this model, presenting some of the tradeoffs and challenges in redox-active material selection and how they affect the entire thermochemical hydrogen production process. This model could be easily adapted and used for material exploration, system/reactor design, and technoeconomic analysis.

08 HYDROGEN↗
DOE OSTI2024

Thermal Fluid Modeling Approaches in SAM for High Temperature Gas-Cooled Reactor Applications

This paper provides an overview of the unique modeling and simulation challenges and needs in the system and safety analysis of two common types of High Temperature Gas-cooled Reactor (HTGR) designs. The challenges are associated with the complex core geometry configurations and the change of the dominant heat removal mechanisms between normal operating conditions and decay heat removal transients. This requires that simulation tools utilize models with reasonable length scales and capture both the localized heat transfer and the core-scale heat transfer at the same time. To meet such analysis needs, different methodologies were developed. In conclusion, this includes adopting existing methods and further developing/improving them, as well as proposing innovative methods that leverage advanced computational frameworks such as the Multiphysics Object Oriented Simulation Environment (MOOSE).

HTGR↗
DOE OSTI2024

Techno-Economic Simulation Results Using dGeo for EGS-Based District Heating in the Northeastern United States

This dataset presents the results of techno-economic simulations performed using the Distributed Geothermal Market Demand Model (dGeo) to evaluate the feasibility of Enhanced Geothermal Systems (EGS)-based district heating in the Northeastern United States. Developed by the National Renewable Energy Laboratory (NREL), dGeo is a geospatially resolved, bottom-up modeling framework designed to explore the deployment potential of geothermal distributed energy resources. The dataset, created as part of the Cornell EGS Ground-Truthing Project, provides census tract-level data that includes inputs and outputs such as thermal demand, road length, energy prices, geothermal system sizing, annual energy contributions from geothermal and natural gas peaking boilers, system capital costs (CAPEX), operation and maintenance costs (OPEX), and the levelized cost of heat (LCOH). Key simulation parameters include geothermal gradients, measured well depths, production temperatures, and district heating piping lengths based on S1400 neighborhood road lengths. The simulations assume a target bottom hole temperature of 80C and the development of new district heating networks in each census tract.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Topology optimization of an airfoil fin microchannel heat exchanger using artificial intelligence

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, is analyzed and optimized using a design of experiments with artificial intelligence 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 printed circuit heat exchangers (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 is developed that employs computational fluid dynamics for a set of design points identified using Latin hypercube sampling. Computational fluid dynamics is used to analyze a simplified two-channel configuration where five design parameters are varied – inlet angle, fin scale, extent of staggering, transverse and longitudinal pitches. Two methods (a 5D polynomial and a regression neural network) are compared for generating surrogate models and the resulting response surface approximation is input to a genetic algorithm that is used to identify a set of optimal parameters. The optimal geometries are found across six channel Reynolds numbers ranging from 1000 to 5000, since inlet conditions affect flow through the heat exchanger. Additionally, a set of optimal designs that maximizes heat transfer and minimizes pressure drop is identified, and a thermal stress analysis is performed on the optimal design. Correlations for the Nusselt number and Darcy friction factor are developed that can be useful for thermal hydraulic analyses using system codes. Thermal stresses are analyzed 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↗
DOE OSTI2026

Interface and Thermophysical Properties of R 32 Refrigerant

Driven by the urgent demand for efficient cooling in microelectronics and advanced thermal management systems, difluoromethane (R32/CH 2 F 2 ) has emerged as a promising candidate owing to its favorable thermophysical properties, including high heat transfer efficiency and low viscosity. While bulk properties such as density, viscosity, and thermal conductivity have been widely studied, interfacial properties, including surface tension and interfacial thickness, remain comparatively underexplored, despite their importance in phase-transition dynamics. Here, we perform molecular dynamics (MD) simulations from 180 to 300 K using an optimized transferable force field for fluoropropenes with enhanced electrostatics to assess both bulk and interfacial behavior of R32. Simulations reproduced density within ±2.1%, viscosity within 3.05%, and thermal conductivity within 7.41% of NIST reference data. Heat capacities (C p and C v ) were predicted within 5%. For interfacial properties, surface tension trends were reproduced within 13.58% deviation, and the vapor–liquid coexistence curve closely matched reference data, yielding a critical temperature of 345.7 K (1.6% deviation) and a critical density of 0.397 g/cm 3 (6.4% deviation). Importantly, the vapor–liquid interface exhibited pronounced temperature-dependent broadening across the 180–290 K range. This behavior correlates with increasing molecular kinetic energy, reduction in intermolecular cohesive interactions, and a progressive loss of preferential dipole alignment, which collectively enhance thermal fluctuation amplitudes at elevated temperatures. These validated results provide predictive molecular-level insights, particularly for interfacial properties that remain less characterized. By reducing property prediction errors in key parameters such as critical temperature, this work provides reliable inputs for heat-exchanger and system models. Such correlations can support optimized component sizing, improved performance, and reduced refrigerant charge. Beyond R32, the methodology offers a transferable framework for blended and next-generation low-GWP refrigerants, contributing to sustainable thermal management aligned with the 2027 EU F-Gas regulation and 2030 Kigali Amendment.

Fluids↗
DOE OSTI2023

Scaling Outlooks on Integrated Energy Systems

The report provides research summaries of two developed theories: data projection to enable multi-scale facility digital interactions and introducing the Evaluation Model Development and Assessment Process (EMDAP) approach to Integrated Energy Systems (IES) design. For the data projection, a two-system scaling case including interactions between thermal storage and hydrogen production is provided. The corresponding facilities are the thermal energy distribution system (TEDS) facility and high temperature steam electrolysis (HTSE) facility at the Dynamic Energy Transport and Integration Laboratory (DETAIL) in Idaho National Laboratory (INL). For the EMDAP adaption to system design, a modified process of the EMDAP was developed and applied to indirect two-tank heat storage system design. System phenomena were identified, ranked, and scaled to discover helpful key relations that would determine the ideal component geometry.

08 HYDROGEN↗
DOE OSTI2024

A study on the impact of using a subchannel resolution for modeling of large break loss of coolant accidents

The nuclear industry is investigating the feasibility of transitioning from 18- to 24-month fuel cycles because of the positive impact it would have on the operational costs for the current fleet of light-water reactors. A challenge to making this change is the increased risk of fuel fragmentation, relocation, and dispersal (FFRD) due to the known potential for ceramic fuel to pulverize into fine particles at the higher discharge burnups. Previous work has been performed by the Nuclear Energy Advanced Modeling and Simulation program to assess FFRD risk in high-burnup cores using the BISON fuel performance code and a coarse mesh thermal hydraulics (T/H) solution for a loss-of-coolant accident (LOCA) using the TRACE system T/H code. Because of the importance of the T/H solution for FFRD assessment, this study seeks to investigate the impact of using higher-fidelity subchannel techniques for modeling of the LOCA transient. CTF was used to model a subregion of a high-burnup core that was depleted by the Virtual Environment for Reactor Applications (VERA) multiphysics core simulator. Both coarse-mesh and pin-resolved models were created in CTF, and a consistent coarse-mesh TRACE model was also developed to allow for benchmarking the code results. Further, a large-break loss-of-coolant accident (LBLOCA) reflood transient was simulated using these three models, and results were compared. Results showed some consistent differences between the CTF and TRACE coarse models, including a higher peak cladding temperature (PCT) prediction in CTF and later quenching in CTF; however, the transient clad temperature behavior was similar, and these differences are likely due to post-critical heat flux heat transfer modeling differences and minimum film boiling temperature model differences. The pin-resolved results indicate that the PCT in the lumped model is often under-predicted by as much as 70 °C and that PCT occurs at a different location than the high-power pin in the assembly. The lumped model predicts a difference of 10 °C or less between the average and hot pins in the assembly, whereas the pin-resolved model predicts a range of over 100 °C. These results indicate that higher-fidelity T/H results may have an impact on predicted core behavior during LOCA, which may be important to consider when assessing FFRD risk.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Exact representations of many-body interactions with restricted-Boltzmann-machine neural networks

Restricted Boltzmann machines (RBMs) are simple statistical models defined on a bipartite graph which have been successfully used in studying more complicated many-body systems, both classical and quantum. In this work, we exploit the representation power of RBMs to provide an exact decomposition of many-body contact interactions into one-body operators coupled to discrete auxiliary fields. This construction generalizes the well known Hirsch's transform used for the Hubbard model to more complicated theories such as pionless effective field theory in nuclear physics, which we analyze in detail. Finally, we also discuss possible applications of our mapping for quantum annealing applications and conclude with some implications for RBM parameter optimization through machine learning.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Adaptive behavior and different thermal experiences of real people: A Bayesian neural network approach to thermal preference prediction and classification

Various observed and unquantifiable factors affect the thermal comfort of occupants in indoor environments and can lead to high uncertainty in the prediction and classification of their thermal preferences. The behavioral adaptation of occupants, by operating window systems for example, changes their thermal experience and expectations and therefore contributes to even higher prediction uncertainty. In this study, we applied a Bayesian neural network (BNN) algorithm to build a predictive model for occupant thermal preference using the ASHRAE Global Thermal Comfort Database II. The Bayesian method allows us to synthesize prior knowledge and available measurements into a unified modeling framework. It also offers a way to express and quantify uncertainty. Here we have performed a systematic study to test the efficiency and robustness of different BNN model configurations. In this study, the results show that the BNN model outperforms conventional thermal comfort models such as Predicted Mean Vote (PMV) and adaptive comfort model. The BNN model tends to produce more confident “prefer cooler” predictions with high possibility and low uncertainty. In contrast, the BNN model produces less certain predictions for “prefer no change” and “prefer warmer” across all occupants. Our findings suggest that linking occupants’ subjective evaluation measures and window opening/closing behavior to thermal comfort modeling effectively improves predictive performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

TNSL Overview

Thermal neutron scattering law (TNSL) data describe low-energy neutrons scattering off of bound materials, and can have a significant impact on modeling any system with slow neutrons, including nuclear reactors. Previous work to introduce TNSL data to neutron transport codes at LLNL focused on COG and TART [1], with the limitation that these codes require highly specialized data processing and formatting. We have recently increased efforts to process TNSL data with the central LLNL nuclear data processing code FUDGE, to be stored in the generalized nuclear database structure (GNDS) for use in any general transport code with the ability to read GNDS data. The first step in this effort is to verify that the TNSL processing with FUDGE yields results comparable to results obtained using the LANL nuclear data processing code NJOY. The next step is to verify the transport of thermal neutrons in Mercury (a Monte Carlo code) and Ardra (a deterministic code) against one another, as well as against the LANL Monte Carlo neutron transport code MCNP. This verification step has not been completed, due to a number of discrepancies between results obtained using differently processed data. There is ongoing effort to understand differences between FUDGE and NJOY. Finally, we map out our current capability to validate TNSL data against benchmark systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Design of a prototypical natural circulation water-based reactor cavity cooling system (RCCS) for a pebble-bed generic FHR

A prototypical natural circulation water-based reactor cavity cooling system (RCCS) for the UC Berkeley Mark-1 pebble-bed gFHR is designed based on one-dimensional thermal hydraulics modeling and optimization implemented in an in-house MATLAB code. The model employs a lumped core consisting of fuel pebble and graphite pebble regions, while the graphite reflector, vessel, and RCCS are represented as separate but energy-coupled regions. The model is derived based on steady state energy balance equations accounting for conductive heat transfer from the lumped core to the vessel and subsequent radiative heat transfer to the RCCS and convective heat transfer in the water. Further, mass flow rate is calculated based on momentum and thermal energy balance in the RCCS. Starting with ANL’s water NSTF as baseline, effects of RCCS design parameters such as size of the plate between pipes, pipe diameter, source-to-sink distance, surface emissivity, and pitch from the core are examined in an integral effects framework which accounts for subsequent variation of the mass flow rate, number of pipes, convection coefficient, mean region and surface temperatures, radiative view factors, and water outlet temperature. The importance of including a prototypical reactor design in RCCS design calculations is emphasized as the calculations show that temperature drops from the core to the RCCS fluid to transfer the heat tend to be the limitation not the heat removal capacity of the fluid. A closed-loop design is obtained with capability to safely remove up to 0.72% of the nominal reactor power at maximum estimated peak conditions and 0.40% at shutdown. The design is based on physics calculations and does not account for economic optimization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Artificial Magnetic Tripod Ice

We study the collective behavior of interacting arrays of nanomagnetic tripods. These objects have six discrete moment states, in contrast to the usual two states of an Ising-like moment. Our experimental data demonstrate that triangular lattice arrays form a “tripod ice” that exhibits charge ordering among the effective vertex magnetic charges, in direct analogy to artificial kagome spin ice. The results indicate that the interacting tripods have effective moments that act as emergent local variables, with strong connections to the well-studied Potts and clock models. In addition, the tripod moments display a tendency toward a nearest neighbor alignment in our thermalized samples that separates this system from kagome spin ice. In conclusion, our results open a path toward the study of the collective behavior of nonbinary moments that is unavailable in other physical systems.

36 MATERIALS SCIENCE↗