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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Two-zone stratified wetwell model development and implementation for RELAP-7

The Reactor Core Isolation Cooling (RCIC) system consists of a turbine, pump and wetwell and is used in many US Boiling Water Reactors (BWRs) as an important heat removal system. The wetwell is a major heat sink within containment. Experimental investigations have shown that thermal stratification may occur in the wetwell during RCIC operation. Current systems-level analysis codes lack stratified wetwell models. Here, this paper describes the development of a stratified wetwell model for the RELAP-7 code. This paper uses a simple two-zone model for the wetwell water that can capture thermal stratification. This model assumes that a buoyant plume transports heat and mass from the steam injection site to the upper water zone. The buoyant plume assumption requires that steam is injected at a low enough rate so it adds negligible momentum. Model validation against experimental data shows similar trends of temperature field development in the suppression pool.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Computational Evaluation of Thermal Response of Open-Cell Foam With Circular Pore (Computional Evaluation of Thermal Response of Open-Cell Foam with Circular Pore)

The evaluation of effective material properties in heterogeneous materials (e.g., composites or multicomponent structures) critically relevant to a wide spectrum of applications, including nuclear power, electronic packaging, flame retardants, hypersonics, and gas turbine power. The work described in this paper is centered around the numerical assessment of the thermal behavior of porous materials obtained from finite element thermal modeling and simulation. Here, two-dimensional, steady state analyses were performed on unit cells with centered, circular pores using a second order accurate Galerkin finite element method (FEM). The effective thermal conductivities of the porous systems were examined, encompassing a range of porosities from 4.9% to 60.1%. The geometries of the models were generated based on ordered circular pores for each modeled porosity level. The system response quantity (SRQ) under investigation was the dimensionless effective thermal conductivity across the unit cell. The dimensionless effective thermal conductivity was compared across all simulated cases, producing a trend between porosity and effective thermal conductivity. In the presented investigation, the method of manufactured solutions (MMS) was used to perform code verification, and the grid convergence index (GCI) was employed to estimate discretization uncertainty as solution verification. Code verification concluded an approximately second order accurate Galerkin FEM solver. It was found that the introduction of porosity to the unit cell material structure reduces effective thermal conductivity, as anticipated. Numerical results obtained in this study are compared to an analytical solution and to a sample of empirical data.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development and Experimental Optimization of High-Temperature Modeling Tools and Methods for Concentrated Solar Power Particle - Systems

A novel, open-source radiative modeling toolset was developed to extend the functionality of particle-based modeling software (e.g. discrete element method (DEM)) to environmental conditions relevant to concentrated solar power applications. This toolset was optimized for deployment on desktop workstations instead of high-performance computing systems, to render such tools more accessible to the research community. Both particle-based modeling and radiative exchange modeling are computationally expensive and often require specialized programming expertise, making these methods cumbersome to use. Recent developments in DEM software by DCS Computing have greatly reduced these challenges, providing a graphical-user-interface based platform and modeling optimization for desktop workstations, HPCs, and cloud computing. The University of Dayton leveraged the experience of DCS Computing in developing a user-friendly, open-source radiative heat transfer expansion for DEM modeling. The University of Dayton DEM+ radiative modeling toolset was developed using a combination of fundamental experimental measurements, modeling, and simplified flow experiments over a range of temperatures and flow conditions. The toolset provides researchers with access to multiple radiative models including an accelerated Monte-Carlo Ray Tracing (application agnostic, highly computationally expensive), an expanded database of distance-based approximations (application limited, computationally light), and a weighted blending of the two methods capable of achieving over 90% reduction in computation time with equivalent accuracy compared to Monte-Carlo Ray Tracing. Through a graphical user interface, users can customize the radiative models to match their desired accuracy and available computational resources, improving access to particle based modeling for the research community. Ceramic sintered bauxite proppants were used in modeling and experimentally as a baseline. Both the radiative heat transfer and flow properties for particulate systems were investigated at elevated temperatures up to 800 °C. The major accomplishments for this work include a verified, open-source radiative modeling toolset to be distributed amongst the research community and the fabrication of three small-scale test facilities to investigate particle behavior and tune DEM flow properties for operation up to 800 °C. The findings have been shared with the research community via conference modeling workshops, deployment of the tools in DCS Computing Aspherix®, and open-source access to the developed radiative modeling tool. The development of next-generation CSP facilities and thermal energy storage systems based on ceramic particles requires providing access to computationally efficient and accurate modeling tools. Particles will experience a wide range of environments (20-800 °C) and handling conditions (dilute curtains or dense packing), requiring specially designed and optimized equipment. Optimizing solid particle physics models and establishing best-practices for particle modeling in CSP environments will assist researchers with designing optimized equipment, accelerating the deployment of more economically-competitive CSP facilities.

14 SOLAR ENERGY↗
DOE OSTI2025

Integral Heat Transfer Test Method for Determining Thermal Resistance of Silicon Carbide Cladding and Liquid Tin System at Elevated Temperatures - Master's Thesis Presentation University of Idaho

Presentation for master's work at University of Idaho discussing Integral Heat Transfer Test Method for Determining Thermal Resistance of Silicon Carbide Cladding and Liquid Tin System at Elevated Temperatures. Presentation covers slight literature review results/previous work, discussions on system design and approach for thermal property measurement of SiC composite and Tin system, modeling efforts, and results.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Thermal coupled cluster theory for SU(2) systems

Coupled cluster (CC) has established itself as a powerful theory to study correlated quantum many-body systems. Finite-temperature generalizations of CC theory have attracted considerable interest and have been shown to work as nicely as the ground-state theory. However, most of these recent developments address only fermionic or bosonic systems. The distinct structure of the su(2) algebra requires the development of a similar thermal CC theory for spin degrees of freedom. In this paper, we provide a formulation of our thermofield-inspired thermal CC for SU(2) systems. Here, we apply the thermal CC to the Lipkin-Meshkov-Glick system as well as the one-dimensional transverse field Ising model as benchmark applications to highlight the accuracy of thermal CC in the study of finite-temperature phase diagrams in SU(2) systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

A predictive analytical model of electrical transport in multi-principal-element alloys

A predictive analytical model is presented for the electrical conductivity of multi-principal-element alloys (MPEAs), including those containing aluminum, transition metals, and refractory metals. Given that the lattice parameter of the Wigner-Seitz cell of an MPEA is similarly variable to a bulk metallic glass, it is postulated that electron scattering can be approximated by a series of two-level systems. Here, the resulting reduced-order model enabled an accurate determination of electrical resistivity and electron thermal conductivity based on the scattering of electrons in a two-level system across a Bloch-potential-based virtual crystal approximation. Model results are compared to experimental four-point probe electrical resistivity measurements between 300 K and 700 K for Al 0.3 CoCrCuFeNi, CoCrFeMnNi, (CoCrFeMnNi) 0.98 W 0.02 , (CoCrFeMnNi) 0.95 W 0.05 , and Nb 4 Ta 4 V 3 Ti, for model validation.

Analytical model↗
DOE OSTI2021

A fast reduced model for a shell-and-tube based latent heat thermal energy storage heat exchanger and its application for cost optimal design by nonlinear programming

Numerical simulation of latent heat thermal energy storage (LHTES) systems plays a fundamental role in studying the physical process and guiding the engineering design. Discretization of the PDEs describing the nonlinear solidification/melting process of phase change materials (PCMs) leads to a large-scale complex dynamics system, where the system behavior depends on a set of parameters. In a design setting, repeated model evaluations are required over the set of parameters results in significant computational burden. In this paper, an explicit analytic solution was built for the propagation of the solidification front in a cylindrical coordinate. The analytic solution approach is further employed to develop a low computational reduced model (RM) as a module for a shell-and-tube based LHTES heat exchanger. The levelized Cost of Energy (LCOE) is used as a design metric and the RM model is used to apply system-level constraints in the nonlinear programming formulation that facilitates efficient global optimal design of the PCM properties, flow conditions and tube geometries. The use of LCOE as the design metric prevents over design of the heat transfer rate and also establishes a fair ground for evaluation of different thermal storage technologies and their integrated applications with other systems. Optimal results showed that a higher effectiveness results in a higher LCOE; the velocity of the HTF and the length of the channel are highly correlated with each other; both larger PCM latent energy and conductivity result in lower LCOE.

25 ENERGY STORAGE↗
DOE OSTI2023

Variance Reduction within Implicit Monte Carlo Thermal Radiation Transport using the Local Importance Function Transform [Slides]

System of equations derived by Fleck and Cummings to model time-dependent thermal radiative transfer (TRT) problems: Solution of radiation specific intensity (I), non-linearly coupled to material internal energy (U); Absorption and re-emission approximated by effective scattering events. Useful for high-energy density physics simulations: i.e., astrophysics, inertial confinement fusion (ICF). Whereas deterministic methods are fully discretized and free of statistical noise, Monte Carlo methods allow for dynamic sampling of the phase space at the price of statistical noise; Monte Carlo methods are also characterized by a slow $\frac{1}{√Ν}$ convergence rate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

THERMAL MODELING OF HANFORD CESIUM AND STRONTIUM CANISTERS DURING SIMULATED LOADING

A computational fluid dynamics (CFD) model was built to simulate planned testing of heater assemblies within a canister and overpack for the Hanford Lead Canister (HLC) project. The HLC is a canister storage system that will contain heaters to simulate the decay heat of nuclear material and provide the canister storage system with environmental conditions equivalent to the operating conditions on a dry storage pad. The HLC will be equipped with long-term data collection and monitoring systems to provide an early warning of corrosion, pitting, cracking, or other signs of canister degradation that might threaten the integrity of the containment boundary over the potentially long term of dry storage. An important part of the HLC development is to make pretest numerical predictions for the behavior of the heated canister during the simulated radiolytic decay heat testing, which simulates the dry storage system during loading operations. The simulated radiolytic decay heat test is planned for mid-2024 in a configuration that includes the heater assembly, overpack, and canister, but with the lids removed to allow loading cesium and strontium capsules into the canister. One of the goals of the test is to evaluate the thermal behavior of the canister and overpack assembly in the ambient air of the test facility, which will provide data critical to validating the thermal models and understanding how the HLC will perform as a system once deployed. To best approximate real-world conditions, the CFD model includes the full air volume of the mock-up truck bay the heated canister test will be performed in, enabling detailed investigation of how the heated canister affects airflow around it. Rigorous pre-deployment testing of the complete HLC cask and canister system is intended to be completed before the HLC is deployed in the 2028 timeframe. This study presents the pre-test temperature predictions of the simulated radiolytic decay heat test. A description of the heater assembly, canister, and overpack system is presented. The model was developed with the commercial CFD software STAR-CCM+. An uncertainty analysis was run with the CFD model to determine the uncertainty in the temperature predictions and provide a range over which the predicted temperatures are expected to vary. The uncertainty analysis was preformed by coupling STAR-CCM+ with the software Dakota, which provides advanced parametric analyses, including quantification of margins and uncertainty with computational models. This work is expected to provide insight into SNF canister behavior.

Carpenter-Graffy, Dina E.↗
DOE OSTI2023

2023 Real-time Optimization Workflow Status Update

Nuclear integrated energy systems are composed of a diverse set of energy generation sources and exist in dynamic and competitive electricity markets. With the inclusion of thermal energy storage, nuclear power systems can store heat for future use through various processes, such as water desalination or hydrogen generation. This heat storage can be managed in such a way as to economically optimize its usage. By combining real-time price data from the day-ahead and real-time markets with predictive and intelligent models, the charging and discharging of the thermal energy storage may be determined and optimized. This research details an approach through models and systems for the real-time optimization (RTO) of capacity allocation. Virtual models of the energy system and its physics phenomenon and component interactions provide intelligence to verification and prediction of operations. An optimization framework can use data generated from both a set of physical assets as well as the virtual models to predict future performance and create optimization and control workflows. Data warehouse technologies can be used to combine data across models, optimization workflows, market price data, and sensor data to intuitively store various types of data and provide integrations to physical control systems as well as user visualizations. Put together, these components can create a system for the RTO of nuclear integrated energy systems. Various virtual models and bench-scale physical demonstrations have been successfully performed and verified using this system. Larger scale testbeds that include thermal energy storage systems have been identified as future opportunities. A gap analysis that details the steps necessary to reach a physical demonstration at this scale is provided, along with conclusions on the current effort and future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Community Geothermal: Borefield Design, Thermal Conductivity, and Subsurface Modeling Data - Chicago, IL

This dataset encompasses the development of a geothermal energy system for the West Woodlawn neighborhood in Chicago, Illinois. This project is part of a broader initiative to design and deploy geothermal heating and cooling systems at a community scale. The dataset includes thermal conductivity test results, calculations for borehole sizing based on the American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) method, as well as simulated thermal loads based on actual energy usage from individual buildings. Also provided here are files used for numerical modeling via COMSOL Multiphysics to simulate borefield design and subsurface thermal behavior. Two manuscripts are attached, which outline the broad objectives of the project and a description of the numerical modeling methodology and results.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Localization crossover and subdiffusive transport in a classical facilitated network model of a disordered interacting quantum spin chain

Here, we consider the random-field Heisenberg model, a paradigmatic model for many-body localization (MBL) and add a Markovian dephasing bath coupled to the Anderson orbitals of the model's noninteracting limit. We map this system to a classical facilitated hopping model that is computationally tractable for large system sizes and investigate its dynamics. The classical model exhibits a robust crossover between an ergodic (thermal) phase and a frozen (localized) phase. The frozen phase is destabilized by thermal subregions (bubbles), which thermalize surrounding sites by providing a fluctuating interaction energy and so enable off-resonance particle transport. Investigating steady-state transport, we observe that the interplay between thermal and frozen bubbles leads to a clear transition between diffusive and subdiffusive regimes. This phenomenology both describes the MBL system coupled to a bath and provides a classical analog for the MBL transition in the corresponding quantum model, in that the classical model displays long local memory times. It also highlights the importance of the details of the bath coupling in studies of MBL systems coupled to thermal environments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Leveraging Polymorphism in YbCuBi to Map Transport and Elastic Properties

AMX Zintl compounds with the hexagonal ZrBeSi structure have gained significant attention for their remarkable vacancy tolerance and low thermal conductivity. Their 2D honeycomb sublattice, composed of M–X covalent bonds, is believed to contribute to high anharmonicity and unusual thermal transport properties. In this study, we explore the temperature-dependent polymorphism of YbCuBi as a model system to investigate the relationship between the structure and elastic and thermal transport properties in AMX Zintls. YbCuBi undergoes a structural transition from the “flat” Cu–Bi layers in the ZrBeSi structure to corrugated layers in the LiGaGe structure below 410 K, resulting in a distortion of its centrosymmetric structure. To probe the effects of this crystallographic transition, we employ inelastic neutron scattering and temperature-dependent resonant ultrasound spectroscopy. These experimental findings, coupled with first-principles calculations and thermal conductivity measurements, allow us to elucidate a direct relationship between corrugation of the honeycomb lattice and the observed changes in elastic and thermal transport properties. These insights can be extended to other Zintl phases with similar structure types, providing a platform for the rational design of functional materials with tailored thermal properties.

Chemical structure↗
DOE OSTI2022

Investigation of the thermal conductivity of molten LiF-NaF-KF with experiments, theory, and equilibrium molecular dynamics

Molten salts are being proposed for numerous advanced energy applications, including advanced nuclear reactors, concentrating solar power plants, thermal energy storage, and fusion reactors. Accurate knowledge of the thermophysical properties of molten salts directly impact the performance of these energy systems and are essential for design and safety analyses. Thermal conductivity data for fluoride molten salts and mixtures are especially lacking. In this work, experimental measurements of thermal conductivity using the steady-state variable gap technique were performed on eutectic LiF-NaF-KF from 834 to 1195 K. The experiment accounts for radiative, convective, and conductive heat losses. In addition, theoretical and molecular dynamics models are used, from 750 K up to 1300 K, to estimate the thermal conductivity for comparison with the experimental results. The results of experiments show a weak negative deviation of thermal conductivity with temperature, unlike previous experimental results in the literature. The measured thermal conductivity magnitudes agree with the theoretical and molecular dynamics predictions, aside from the data above 1100 K, where heat losses and radiative errors are the most significant, having a 16% maximum deviation from theory. These experimental results provide new thermal conductivity data for the LiF-NaF-KF system and further validation of the predictive models. The theoretical model was used to map the composition and temperature dependent thermal conductivity of LiF-NaF-KF and the mapping’s deviation from a linear additivity estimation of thermal conductivity. Additionally, this mapping showed the highest deviations from linearity for KF-LiF rich mixtures and increasing deviation with temperature. Notably, the deviation from linearity near the LiF-NaF-KF eutectic composition was around 25%.

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

Real-time Simulation Model of a Vanadium Redox Flow Battery Energy Storage System with Power Electronics Integration

This paper presents a real-time (RT) model of a 20 kW multi-stack vanadium redox flow battery (VRFB) system including power electronics converters in a single real-time framework. Implemented on a Typhoon HIL 604 device, this model is designed to simulate the operation of the flow battery with system-level controls. To provide a realistic overview of a VRFB system, the model incorporates a closed-loop control system for managing the speeds of two centrifugal pumps and an electro-thermal model of the battery stacks, enabling realtime temperature estimation. Furthermore, a two-stage power electronics topology and associated control solution are presented, operating with distinct strategies for interfacing the battery model with the main grid and a local load. Results from simulations demonstrate that the VRFB model is well suited for RT environment, and system can manage power during both charging and discharging cycles, while a SCADA monitor is used to display critical variables for safe, efficient, and reliable operation of the VRFB systems.

Rezende Da Costa Reis Kimpara, Renata [ORNL] (ORCI↗
DOE OSTI2021

Cabin Thermal Management Analysis for SuperTruck II Next-Generation Hybrid Electric Truck Design

In this article, we present a multistage, coupled thermal management simulation approach, informed by physical testing where available, to aid design decisions for PACCAR's SuperTruck II hybrid truck cabin concept. Focus areas include cabin insulation, battery sizing, and sleeper curtain position, as well as heating, ventilating, and air-conditioning (HVAC) component and accessory configurations, to maintain or improve thermal comfort while saving energy. The authors analyzed weather data and determined the national vehicle miles traveled weighted temperature and solar conditions for long-haul trucks. Example weather day profiles were selected to approximate the 5th and 95th percentile weighted conditions. A daylong drive cycle was developed to impose appropriate external wind conditions during rest and driving periods. Using the National Renewable Energy Laboratory's vehicle HVAC modeling and simulation tool VTCab, HVAC load design trade-off studies for the new truck geometry concept were completed. Parameters analyzed included effects of paint color, insulation, glass transmissivity, and curtain location. Simulation results helped with early design material selections for efficient cabin climate control. A detailed three-dimensional computer-aided engineering (CAE), computational fluid dynamics (CFD), radiation, and human physiology co-simulation, referred to in this article as CAE Thermal-CFD, was used to evaluate thermal comfort and energy impacts of diffuser configurations and air supply settings in driving and hoteling modes. Analysis revealed that it is more difficult to heat the cabin in hoteling mode during the winter than to cool the space in the summer. This seasonal load profile drives the requirement of additional energy storage for heating comfort. To determine the battery capacity requirement, multiday HVAC operation drive cycle simulations were then completed, showing that a 15-kWh battery would be required for HVAC operation during hoteling. Results helped reduce cabin thermal loads, determine component sizing requirements, and improve occupant comfort to save fuel and contribute to the economic viability of the hybrid system.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Stability analysis of a one-dimensional multiphysics model of a molten salt fast reactor

Reactor designs with a fast neutron spectrum have recently moved into the research focus for possible Generation IV reactor concepts - a representative of which is the molten salt fast reactor (MSFR). In order to assess the safety and reliability of this reactor concept, it is necessary to obtain knowledge of the system behaviour. An important part of this research concerns the dynamical stability of the MSFR. Research regarding the MSFR stability so far has either been using linearized models for stability criteria or only been examining certain transients at fixed parameter values. This work delivers a comprehensive stability analysis for a non-linearized MSFR model, and a wide range of values for all relevant parameters. For this purpose, a one-dimensional MSFR model was set up, taking neutron kinetics and thermal hydraulics in the form of a system of differential equations into account. The stability of this model was investigated by means of the numerical tool MATCONT, which was used to monitor the evolution of a so-called fixed-point solution, here referring to the steady state at operating conditions of the system, while system parameters got varied. The MATCONT analyses showed no loss of stability for any of the considered parameter variations and no solutions that might co-exist in parallel with the steady-state fixed-point solution. Therefore these results indicate a stable fixed point to which all solution transients converge, making the considered MSFR model insusceptible to deviations from the equilibrium. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

A review on the sizing and selection of control valves for thermal hydraulics for reactor system applications

This study presents an overview regarding how to model/size, choose, and maintain control valves (CVs) for a reactor system thermal hydraulics experimentation—specifically simulating pipe break(s) and a loss of coolant accident (LOCA) analysis. In a thermal hydraulics test loop, CVs modulate fluid flow according to the control signal generated by the controller to control the process parameter. The selection of a CV in industrial applications involves many factors, including process diversity, safety, and reliability. Further, the key issues and challenges for performing proper mathematical modeling, installation, servicing, and calibration are also discussed for the benefit of operators and end-users. Additionally, this research includes case studies illustrating major plant-level accidents associated with the malfunctioning of valves. Therefore, the study will serve as a knowledge base or reference guide for young professionals and operators who wish to better understand process control applications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗