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 1,135 records · Page 63

A drift kinetic model for the expander region of a magnetic mirror

Here, we present a drift kinetic model for the free expansion of a thermal plasma out of a magnetic nozzle. This problem relates to plasma space propulsion systems, natural environments such as the solar wind, and end losses from the expander region of mirror magnetically confined fusion concepts such as the gas dynamic trap. The model incorporates trapped and passing orbit types encountered in the mirror expander geometry and maps to an upstream thermal distribution. This boundary condition and quasineutrality require the generation of an ambipolar potential drop of ~5T e /e, forming a thermal barrier for the electrons. The model for the electron and ion velocity distributions and fluid moments is confirmed with data from a fully kinetic simulation. Finally, the model is extended to account for a population of fast sloshing ions arising from neutral beam heating within a magnetic mirror, again resulting in good agreement with a corresponding kinetic simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiscale Thermal-hydraulic analysis of the MARVEL micro-reactor using coupled MOOSE Subchannel (SCM) and SAM

MARVEL is a natural-convection-cooled sodium-potassium microreactor that is anticipated to generate 85 kilowatts of thermal energy. It will operate within Idaho National Laboratory’s Transient Reactor Test Facility and is being developed by the DOE Microreactor Program. MARVEL will be used to test microreactor applications, evaluate systems for remote monitoring, and develop autonomous control technologies. A thermal-hydraulic computational model of this facility is a valuable tool to study important transients and calculate the safety limits of the micro-reactor design. For this purpose, the authors have chosen to use a multiscale coupled simulation: SCM for modeling the reactor core and SAM for the reactor’s primary cooling system. SCM is MOOSE physics module for subchannel analysis, which was designed to model single-phase flows through liquid-metal cooled, wire-wrapped fuel pin sub-assemblies, ordered in a triangular lattice. The SCM code was modified to be able to model MARVEL’s unique geometry. SAM is a systems analysis module based on the MOOSE framework. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for various advanced reactor types. The coupling between the two SCM and SAM for MARVEL modeling is done implementing a domain over-lapping approach. The resulting coupled simulation can model transients such as reactor startup/shutdown and provide an intermediate fidelity picture of the temperature field and other variables, in the core. Results for the steady-state simulations are presented in the article as well as flow blockage transient.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Multiscale Thermal-hydraulic analysis of the MARVEL microreactor using a coupled SAM and SCM simulation. (PPTX)

This is .pptx document presenting a summary of the paper of the same name that has been submitted and accepted to NURETH-21: MARVEL is a natural-convection-cooled sodium-potassium microreactor that is anticipated to generate 85 kilowatts of thermal energy. It will operate within Idaho National Laboratory Transient Reactor Test Facility and is being developed by the DOE Microreactor Program. MARVEL will be used to test microreactor applications, evaluate systems for remote monitoring, and develop autonomous control technologies. A thermal-hydraulic computational model of this facility is a valuable tool to study important transients and calculate the safety limits of the micro-reactor design. For this purpose, the authors have chosen to use a multiscale coupled simulation: SCM for modeling the reactor core and SAM for the reactor's primary cooling system. SCM is MOOSE physics module for subchannel analysis, which was designed to model single-phase flows through liquid-metal cooled, wire-wrapped fuel pin sub-assemblies, ordered in a triangular lattice. The SCM code was modified to be able to model MARVEL?s unique geometry. SAM is a systems analysis module based on the MOOSE framework. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for various advanced reactor types. The coupling between the two SCM and SAM for MARVEL modeling is done implementing a domain over-lapping approach. The resulting coupled simulation can model transients such as reactor startup/shutdown and provide an intermediate fidelity picture of the temperature field and other variables, in the core. Results for the steady-state simulations are presented in the article as well as flow blockage transient.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Bearing tester data compilation, analysis and reporting and bearing math modeling, volume 1

Thermal and mechanical models of high speed angular contact ball bearings operating in LOX and LN2 were developed and verified with limited test data in an effort to further understand the parameters that determine or effect the SSME turbopump bearing operational characteristics and service life. The SHABERTH bearing analysis program which was adapted to evaluate shaft bearing systems in cryogenics is not capable of accommodating varying thermal properties and two phase flow. A bearing model with this capability was developed using the SINDA thermal analyzer. Iteration between the SHABERTH and the SINDA models enable the establishment of preliminary bounds for stable operation in LN2. These limits were established in terms of fluid flow, fluid inlet temperature, and axial load for a shaft speed of 30,000 RPM.

Marshall, D. D.↗

Experimental Effervescence and Freezing Point Depression Measurements of Nitrogen in Liquid Methane-Ethane Mixtures

NASA is designing an unmanned submarine to explore the depths of the hydrocarbon-rich seas on Saturn's moon Titan. Data from Cassini indicates that the Titan north polar environment sustains stable seas of variable concentrations of ethane, methane, and nitrogen, with a surface temperature near 93 K. The submarine must operate autonomously, study atmosphere/sea exchange, interact with the seabed, hover at the surface or any depth within the sea, and be capable of tolerating variable hydrocarbon compositions. Currently, the main thermal design concern is the effect of effervescence on submarine operation, which affects the ballast system, science instruments, and propellers. Twelve effervescence measurements on various liquid methane-ethane compositions with dissolved gaseous nitrogen are thus presented from 1.5 bar to 4.5 bar at temperatures from 92 K to 96 K to simulate the conditions of the seas. After conducting effervescence measurements, two freezing point depression measurements were conducted. The freezing liquid line was depressed more than 15 K below the triple point temperatures of pure ethane (90.4 K) and pure methane (90.7 K). Experimental effervescence measurements will be used to compare directly with effervescence modeling to determine if changes are required in the design of the thermal management system as well as the propellers.

Submarine↗

Solving Problems With SINDA/FLUINT

SINDA/FLUINT, the NASA standard software system for thermohydraulic analysis, provides computational simulation of interacting thermal and fluid effects in designs modeled as heat transfer and fluid flow networks. The product saves time and money by making the user's design process faster and easier, and allowing the user to gain a better understanding of complex systems. The code is completely extensible, allowing the user to choose the features, accuracy and approximation levels, and outputs. Users can also add their own customizations as needed to handle unique design tasks or to automate repetitive tasks. Applications for SINDA/FLUINT include the pharmaceutical, petrochemical, biomedical, electronics, and energy industries. The system has been used to simulate nuclear reactors, windshield wipers, and human windpipes. In the automotive industry, it simulates the transient liquid/vapor flows within air conditioning systems.

Source record↗

Refractory-based thermal energy storage for industrial process heat: one-dimensional modeling, control, and optimization

The variable and weather-dependent output of wind and solar power plants present a substantial challenge for planning and operating electricity-systems, particularly in the absence of cost-effective and dispatchable energy storage technologies. This study investigates a high-temperature, electrically heated, refractory-based thermal energy storage (RTES) system that stores electrical energy as sensible heat in dense ceramic bricks over the 950–1800 °C range. The stored heat can be discharged as a controlled hot-gas stream for industrial heating, fuel substitution in high-temperature processes, or electricity generation. The main novelty is a comprehensive modelling, control, mapping, and optimization framework that integrates one-dimensional transient gas–solid heat transfer, fan-assisted discharge, bypass-flow regulation, reheating logic, fan-power evaluation, insulation-loss assessment, and genetic-algorithm-based design optimization. The model uses feedback from outlet temperature and delivered power to regulate discharge, while a two-stage genetic algorithm optimizes brick-channel geometry, gas-flow operation, and multilayer insulation thicknesses. Storage capacities below 50 MWh and discharge powers of 5–30 MW are analyzed to evaluate hold time, thermal delivery, fan-power penalty, heat loss, state-of-charge evolution, and indicative capital cost. Results demonstrate that optimized and well-insulated refractory-based thermal energy storage units can provide stable, efficient, and repeatable heat delivery over multiple discharge cycles. The generated performance and cost maps support modular refractory thermal energy storage as a practical option for large-scale integration of wind and solar generation and for high-temperature industrial process heat.

25 ENERGY STORAGE↗

Particle-Based Fiber Models of Woven Materials for Earth Entry Thermal Protection

Applications requiring materials with layer-to-layer strength, from basketry to thermal protection systems, use interlaced, three-dimensional woven materials. NASA is developing and deploying woven material heat shields for missions, including Artemis I (3D-MAT for compression pads) and potentially Mars Sample Return - Earth Entry System. These materials are complex, hierarchical and must protect from extreme environments and phenomena, such as deformation, impact and high-enthalpy heating. Woven material performance depends on microstructure, damage and weave geometry. Therefore, fiber-specific models are needed to simulate fiber contacts within the weave hierarchical geometry (fiber to tow to yarn to weave) and the inherent directionality of fibers. Explicit fiber models can simulate how weave microstructure evolution affects thermal and mechanical properties. We parameterize a discrete element bonded particle models (DEM-BPM) of fibers to capture thermal and mechanical behavior within and between fibers, with bonded and contact forces, respectively. We study the proportion of heat transfer and stress via the contact network, fiber bonds and the weave geometry, for example, with respect to yarn warp-weft identity (whether it interlaces weave layers). Our results demonstrate the importance of explicit fiber modeling for connecting microstructure with thermal and mechanical properties.

woven material↗

Particle-Based Fiber Models of Woven Materials for Earth Entry Thermal Protection

Applications requiring materials with layer-to-layer strength, from basketry to thermal protection systems, use interlaced, three-dimensional woven materials. NASA is developing and deploying woven material heat shields for missions, including Artemis I (3D-MAT for compression pads) and potentially Mars Sample Return - Earth Entry System. These materials are complex, hierarchical and must protect from extreme environments and phenomena, such as deformation, impact and high-enthalpy heating. Woven material performance depends on microstructure, damage and weave geometry. Therefore, fiber-specific models are needed to simulate fiber contacts within the weave hierarchical geometry (fiber to tow to yarn to weave) and the inherent directionality of fibers. Explicit fiber models can simulate how weave microstructure evolution affects thermal and mechanical properties. We parameterize a discrete element bonded particle models (DEM-BPM) of fibers to capture thermal and mechanical behavior within and between fibers, with bonded and contact forces, respectively. We study the proportion of heat transfer and stress via the contact network, fiber bonds and the weave geometry, for example, with respect to yarn warp-weft identity (whether it interlaces weave layers). Our results demonstrate the importance of explicit fiber modeling for connecting microstructure with thermal and mechanical properties.

woven material↗

Computer Modeling of Thermal Convection in Melts to Explain Glass Formation in Low Gravity and on Earth

Experiments conducted up to this time on glass forming melts in the low gravity environment of space show that glasses prepared in low-g are more chemically homogeneous and more resistant to crystallization than the comparable glasses prepared at 1-g on Earth. This result is somewhat surprising and opposite to the accepted concept on glass formation for a melt. A hypothesis based on "shear thinning" of a melt, a decrease in viscosity with increasing shear stress, is proposed as an explanation for the observed low-gravity results. This paper describes detailed simulation procedures to test the role of thermal convection in introducing shear stress in glass forming melts, using a lithium disilcate melt as a model. The simulation system in its idealized version consists of a cylinder that is heated at one end and cooled at the other with gravity acting in a transverse direction to the thermal gradient. The side wall of the cylinder is assumed to be insulating. The governing equations of motion and energy are solved using variable properties for viscosity (Arrehenius and non-Arrehenius behaviors) and density (constant and temperature dependent). Other parametric variables in the calculations include gravity level and gravity vector orientation. The shear stress in the system are then computed as a function of gravity from the calculated values of maximum melt velocity, and its effect on melt viscosity (shear thinning) is predicted. Also included and discussed are the modeling efforts related to other potential convective processes in glass forming melts and their possible effects on melt viscosity.

Ray, Chandra S.↗

Electric Water Heaters for Transactive Systems: Model Evaluations and Performance Quantification

Electric water heaters (EWHs) are opportune appliances for implementing demand-side control. EWH models serve as a fundamental step toward accurately estimating EWH flexibility potential and designing proper control strategies. Existing studies have adapted numerous modeling approaches in evaluating the potential of EWHs for a variety of grid applications. This paper presents an analytical study that evaluates the performance of state-of-art EWH models in terms of accuracy and computational complexity for adaptation in evaluation studies for the transactive system. The work proposes a transactive control strategy that optimally utilizes the thermal inertia of EWHs for providing grid services. Here, the performance of the control strategy and the impact of modeling accuracy is evaluated for device-level and feeder-level use cases using the IEEE 123 node test distribution system appropriately populated with EWHs. The simulation results illustrate the effectiveness of the control strategy in reducing the feeder demand during peak period by 13% and also quantities the impact of using simplified modeling approaches for determining the potential of EWHs for providing grid services.

42 ENGINEERING↗

The application of active controls technology to a generic hypersonic aircraft configuration

Analytical methods are described for the prediction of aerothermoelastic stability of hypersonic aircraft including active control systems. Thermal loads due to aerodynamic heating were applied to the finite element model of the aircraft structure and the thermal effects on flutter were determined. An iterative static aeroelastic trim analysis procedure was developed including thermal effects. And active control technology was assessed for flutter suppression, ride quality improvement, and gust load alleviation to overcome any potential adverse aeroelastic stability or response problems due to aerodynamic heating. A generic hypersonic aircraft configuration was selected which incorporates wing flaps, ailerons, and all moveable fins to be used for active control purposes. The active control system would use onboard sensors in a feedback loop through the aircraft flight control computers to move the surfaces for improved structural dynamic response as the aircraft encounters atmospheric turbulence.

Gilbert, M. G.↗

Equilibration and “Thermalization” in the Adapted Caldeira–Leggett Model

I explore the processes of equilibration exhibited by the Adapted Caldeira–Leggett (ACL) model, a small unitary “toy model” developed for numerical studies of quantum decoherence between an SHO and an environment. I demonstrate how dephasing allows equilibration to occur in a wide variety of situations. While the finite model size and other “unphysical” aspects prevent the notions of temperature and thermalization from being generally applicable, certain primitive aspects of thermalization can be realized for particular parameter values. I link the observed behaviors to intrinsic properties of the global energy eigenstates, and argue that the phenomena I observe contain elements which might be key ingredients that lead to ergodic behavior in larger more realistic systems. The motivations for this work range from curiosity about phenomena observed in earlier calculations with the ACL model to much larger questions related to the nature of equilibrium, thermalization, and the emergence of physical laws.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Flow instabilities in helical-coil steam generators for small modular reactors: A review

Here, this study covers the research and discoveries in two-phase flow-boiling instabilities available in the literature—specifically for a helical-coil steam generator (HCSG), including experimental findings, theoretical research, computational models, and system code analyses—supporting research and development of representative small modular reactors (SMRs). Like other new and advanced reactor systems, water-cooled SMRs require experimental data from both integral and separate thermal-hydraulics test facilities for the verification and validation (V&V) of the computational models and computer codes in order to design and obtain regulatory approval. The complex dynamics of two-phase flow-boiling instabilities includes flow regimes physics phenomena, flow-channel geometries, heat-transfer behavior, and interactions among the solid–liquid-gas within the system boundary, all of which are pivotal for understanding the design and operational challenges of SMRs. This study focuses on identifying the relevant knowledge gaps on boiling instabilities—specifically for a HCSG—and provides insights about future research direction optimizing the transport of thermal energy, mass-flow rates, and boundary conditions that ensure the adequate heat-transfer performance, operational stability, and safety associated with SMR systems.

20 FOSSIL-FUELED POWER PLANTS↗

Multi-Scale Thermal Analysis for Design of SiC-Based Medium Voltage Motor Drive: Preprint

This paper describes a multi-scale thermal analysis approach for the design of an air-cooled 1.7-kV SiC MOSFET-based medium-voltage variable-speed motor drive. The National Renewable Energy Laboratory supported the project objectives by providing thermal management design support to the project partner The Ohio State University (OSU). The scope of the models and limited resources required efficient and flexible thermal models to be developed. Two modeling techniques are described that significantly reduced model run time and enabled more complex models to be run faster while retaining needed accuracy. The first technique uses the effectiveness-NTU method to extract convection boundary conditions from a computational fluid dynamics (CFD) model that can be applied to a fast-running finite element analysis model. The second is a porous media technique that enables system-level CFD simulations that incorporate effects from heat exchangers (e.g., pin fin heat sinks) that run in a fraction of the time required for fully resolved CFD simulations. The multi-scale approach to the thermal analysis enabled fast and accurate simulation for the converter design ranging from the die level up to the full system with 36 submodules. The modeling results were validated against experimental data from system tests performed by OSU.

DIRECT ENERGY CONVERSION↗

GridCoPilot for Thermal Events: An LLM-Based Platform for Power Grid Reliability Analysis

Large Language Models show promise for translating natural language into database queries, but deploying such systems in safety-critical domains requires high reliability. We present an application of GridCoPilot to thermal event analysis (heatwaves and coldwaves) that affect power grid reliability. Our approach uses a LangChain SQL Agent to translate natural language queries into auditable SQL statements, with deterministic visualization routines that parse the structured query results. We introduce structural framing as a design principle, we integrate a NERC-region-level event library with county-level meteorology and decompose the combined data into three relational tables (event metadata, county-level event details, and a county-to-NERC subregion mapping), using prompt-guided joins to direct the model toward correct multi-table queries. For two core analytical patterns (identifying worst events by region and by region-year), the system achieved 100% SQL accuracy across all 16 NERC subregions and both event types (64 queries total). These results validate the approach for target use cases, though performance on diverse natural language formulations requires further investigation. We discuss design trade-offs, failure modes including JSON output truncation, and pathways for extending this approach to other hazard domains.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Thermal Property Modeling and Assessment of the Physical Properties of FLiNaK

Here, the thermodynamic properties of the LiF–NaF–KF system and its subsystems were re-evaluated using the CALPHAD method, focusing on the deviations from ideality in the heat capacity [C p (T)] of mixtures, temperature-dependence of mixing enthalpies (Δ mix H), vapor pressures of mixtures, and enthalpies of fusion (ΔH fus ) using reported values and measurements. The results of this work confirmed a pseudoternary eutectic at 732.9 K and 46.5LiF–11.5NaF–42KF mol % (FLiNaK) with determined Cp(T) = 69.822 + 0.000679T + 2,137,152T –2 J mol –1 K –1 (732.9 < T < 1200 K) and ΔH fus = 18,051 J mol –1 in good agreement to the experimental values. The resulting accurate thermodynamic representation can be used to find relevant excess thermochemical properties in the whole pseudoternary space and provide insights into their composition and temperature dependence. A detailed evaluation of select thermophysical properties was also conducted to determine recommended density, thermal diffusivity, thermal conductivity, vapor pressures, and boiling point with calculated uncertainties over the temperature range of interest for FLiNaK–salt energy applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assessment of Molten Eutectic LiF-NaF-KF Density through Experimental Determination and Semiempirical Modeling

Molten salts have favorable material properties for use in high-temperature energy systems, including thermal energy storage systems, concentrating solar power plants, nuclear reactors, and various industrial manufacturing processes. Knowledge of chemical and thermophysical property data is essential for the design and optimization of these systems, yet data are often limited or uncertain for many candidate salts due to the difficulty of thermophysical property measurements at relevant temperatures (e.g., 500–900 °C). Here, the density of molten LiF-NaF-KF eutectic is reassessed through review of previous experimental data, new density measurements from 470 to 800 °C, and semiempirical modeling. Additionally, the density was measured using the displacement technique. Compositional and temperature-dependent density estimates were calculated with a multidimensional Redlich–Kister model. The results of the new experimental measurements agree within 2% of the modeled density of molten eutectic LiF-NaF-KF. The Redlich–Kister model’s prediction shows a near-ideal density behavior for the LiF-NaF-KF system and is promising for the estimation of off-eutectic LiF-NaK-KF densities. Finally, through review of the existing literature and comparison to new measurements, recommendations are made for the density of LiF-NaF-KF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗