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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Depinning, melting, and sliding of generalized Wigner crystals in Moiré systems

We numerically examine the depinning, sliding, and melting of commensurate and incommensurate Wigner crystals on two-dimensional hexagonal periodic substrates near fillings of 1 / 3 , 1 / 2 , and 2 / 3 to model the dynamics of generalized Wigner crystals in moiré heterostructures. At low temperatures where thermal fluctuations are irrelevant, for commensurate fillings we find a strongly pinned state that depins into a sliding crystal, while at incommensurate fillings the depinning threshold is strongly reduced. For fillings below 1 / 2 , the depinning occurs in a two-step process. Above the first depinning threshold, there is an extended range of drives where the conduction occurs via the sliding of antikinks along the charge stripes, followed by a second threshold where all the charges begin to slide. At finite temperatures, the external driving reduces the effective melting temperature at commensurate fillings and enhances creep at incommensurate fillings. We show that depinning into different sliding states, such as moving fluids or moving crystals, produces nonlinear features in the transport curves. We also show that transport is asymmetric on either side of a commensurate filling due to the different dynamics for interstitials versus holes in the commensurate structure. A variety of sliding states should be accessible for generalized Wigner crystals at finite temperatures even for low drives. If experiments can realize sliding dynamics in moiré systems, it would open another class of commensurate and incommensurate phases for study. Additionally, the sliding dynamics and antikink flow could provide unique functionalities for charge transport in moiré heterostructures. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Density functional theory calculations of the mixing enthalpy of ternary uranium carbide compounds

The high melting point of uranium-zirconium carbides (U,Zr)C makes them an ideal fuel for nuclear thermal propulsion (NTP) reactors. Gaps remain in the current understanding of the U-Zr-C system due to the difficulty of conducting thermodynamic experiments at NTP operation conditions. Density functional theory calculations using the Hubbard U model (DFT+U) were performed using orbital matrix occupation (OMC) to obtain the mixing enthalpy for UC and ZrC for (U,Zr)C ternary compounds. Similarly, DFT+U calculations were also carried out for the (U,Nb)C and (U,Ta)C systems. In conclusion, the DFT results are envisioned to be used in thermodynamic assessments of the uranium carbide systems based on the CALPHAD approach to supplement the lack of experimental data for the mixing thermodynamics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Optimal electric-distribution-grid planning considering the demand-side flexibility of thermal building systems for a test case in Singapore

The planning of district-scale electric grids, i.e., distribution grids, has traditionally relied on finding the most cost-effective design such that they are able to supply the peak loads in a district. With the advent of electric demand side flexibility (DSF), there is the opportunity to reshape peak loads such that the investment cost of the electric grid decreases in exchange for a minor increase in the operation cost. This paper formulates an optimal planning approach for the electric grid at the district scale, which incorporates the DSF from thermal building systems, e.g., heating ventilation and air-conditioning (HVAC) systems. The problem is formulated as a mixed-integer linear program (MILP) and aims at minimizing the investment cost for the grid along with the operation cost of the flexible loads. This is subjected to the fixed electricity demand and thermal comfort constraints of building occupants. To this end, linear models for the thermal comfort in the buildings and the power flow in electric grid are considered. The approach is tested on a district planning test case based in Singapore, where the results show up to 30.9 % reductions in investment cost and up to 3.7 % reduction in total annualized cost. Urban planning authorities, developers and utility companies can all benefit from the presented approach to make optimized investment decisions. For building operators, the results point to the need of adopting their control systems for DSF.

Troitzsch, Sebastian↗
DOE OSTI2024

Demonstration of a multi-channel fluidized bed particle–supercritical carbon dioxide heat exchanger for concentrating solar applications

High-temperature thermal energy storage in oxide particles at temperatures above 600°C can couple concentrated solar energy with high-efficiency thermal power cycles to provide dispatchable solar-driven electricity. Challenges remain in developing cost-effective primary heat exchangers, which require expensive alloys, to extract the high-temperature thermal energy from the particles to power cycle fluids, such as supercritical CO 2 (sCO 2 ) in recuperated Brayton cycles. To explore one pathway for cost-effective, high-temperature particle heat exchangers, the current study demonstrates a shell-and-plate, particle–sCO 2 heat exchanger with narrow- channel fluidized beds coupled with micro-channel sCO 2 flows in the heat exchanger walls. This study evaluates the feasibility of multiple parallel, narrow-channel fluidized beds in shell-and-plate particle–sCO 2 HXs, to achieve high bed-wall heat fluxes at elevated temperatures. A reduced-order model simulates the narrow- channel, fluidized-bed particle–sCO 2 heat exchanger to design the fluidized bed geometry, in terms of depth, height, and number of channels,for a nominal 40-kWth heat exchanger at particle and sCO 2 inlet temperatures up to 600 °C and 400 °C respectively. The resulting shell-and-plate heat exchanger design operates with bubbling fluidization of the downward-flowing oxide particles to enhance bed-wall heat transfer. The heat exchanger core is fabricated with etched sCO 2 micro-channels in thin wall plates that are diffusion bonded to spacer frames to form the shell-and-plate structure with 12 parallel, fluidized bed channels, 10.4 mm deep. The heat exchanger is tested at the National Solar Thermal Test Facility at Sandia National Laboratories with CARBOBEAD HSP particles at design particle flow rates of 0.20 kg s –1 and inlet temperatures up to 525 °C. Results show that fluidization across multiple parallel channel beds can maintain uniform particle inventory with a common freeboard zone above the heat exchanger core. Bubbling fluidization improves particle–wall heat transfer coefficients but also increases axial dispersion of particle thermal energy, which lowers the log- mean temperature difference such that total heat transfer remains relatively constant to within ±10% over a broad range of fluidization gas velocities. The axial dispersion required particle and sCO 2 flow rates to be increased by 25% over model-designed conditions to achieve the targeted 40 kWth, which indicates the importance of incorporating axial dispersion into heat exchanger design models and of deploying bed structures to suppress it. Furthermore, this study demonstrates the feasibility and preferred fluidizing gas conditions for particle heat exchangers for releasing high-temperature thermal energy storage systems.

14 SOLAR ENERGY↗
DOE OSTI2026

Propane Dehydrogenation Catalyzed by Supported Group IV (Ti, Zr, Hf) Organometallics on Silicon Nitride

Mesoporous silicon nitride (Si 3 N 4 ) enables access to chemisorbed group IV organometallics catalysts active for propane dehydrogenation (PDH) compared to the organometallic analogues on mesoporous silica under the same reaction conditions. The series of Si 3 N 4 -supported materials are active catalysts, (Zr > Hf > Ti k f = 290, 232, and 162 mol mol Metal -1 h -1 at 450 °C with 2% C 3 H 8 in Ar, respectively) with selectivity above 95%, demonstrating additional examples of Ti and Hf systems for PDH. However, the underlying mechanism of the improved performance relative to oxide supported homologues is not well-understood. Characterization of thermally treated samples (DRIFTS, XAS and SSNMR) and computational modeling of this catalyst series was utilized to differentiate between potential amido - (C-H activation along the M-N bond) and imido - (C-H activation along the M=N bond) mechanisms. Due to remaining mechanistic ambiguity, a Ga analogue was synthesized and evaluated for PDH activity as an indirect probe to experimentally differentiate pathways. An inversion of the oxide/nitride performance trend is observed for the Ga congener which does not form a Ga=N bond, most consistent with different mechanisms dictating the performance of the group IV/Si 3 N 4 catalysts vs Ga/Si 3 N 4 .

heterolytic cleavage↗
DOE OSTI2020

Experimental and theoretical studies of the doubly substituted methyl-ethyl Criegee intermediate: Infrared action spectroscopy and unimolecular decay to OH radical products

In this work the infrared (IR) action spectrum of the doubly substituted methyl-ethyl Criegee intermediate (MECI) is observed in the CH stretch overtone region with detection of OH products. The MECI exhibits four conformers, all of which undergo unimolecular decay via a 1,4 H-atom transfer mechanism, followed by the rapid release of OH products. Conformers with different orientations of the carbonyl oxide group with respect to the methyl and ethyl substituents (i.e., anti and syn ) decay via distinct transition state barriers (16.1 kcal mol –1 and 15.4 kcal mol –1 , respectively). The observed IR action spectrum is in good agreement with the predicted anharmonic IR absorption spectrum, but exhibits significant congestion, which is attributed to couplings between spectroscopic bright states and nearby dark states. Energy-dependent OH appearance rates are measured upon IR excitation of the strongest features in the IR action spectrum and are found to be on the order of 10 6 –10 7 s –1 . The experimental rates are in good agreement with computed Rice–Ramsperger–Kassel–Marcus rates for the unimolecular decay of MECI at these energies, which incorporate quantum mechanical tunneling and sophisticated hindered rotor treatments, as well as high-level theoretical calculations of the TS barrier heights, rovibrational properties, and torsional barriers associated with the MECI conformers. Master equation modeling is used to predict thermal rates for the unimolecular decay of anti - and syn -MECI of 473 s –1 and 660 s –1 , respectively. Comparison with other previously studied Criegee intermediate systems provides insights into substituent effects on unimolecular decay under both energy-dependent and thermal conditions.

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

SAM: A Modern System Code for Advanced Non-LWR Safety Analysis

The System Analysis Module (SAM), developed at Argonne National Laboratory and by collaborators at other organizations, is for advanced non–light water reactor safety analysis. SAM aims to provide fast-running, modest-fidelity, whole-plant transient analysis capabilities that are essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. To facilitate code development, SAM utilizes the MOOSE object-oriented application framework, its underlying finite element library, and linear and nonlinear solvers to leverage modern advanced software environments and numerical methods. SAM aims to solve tightly coupled physical phenomena, including fission reaction, heat transfer, fluid dynamics, and thermal-mechanical responses in advanced reactor structures, systems, and components with high accuracy and efficiency. Finally, this paper gives an overview of the SAM code development, including goals and functional requirements, physical models, current capabilities, verification and validation, software quality assurance, and examples of simulations for advanced nuclear reactor applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Thermal Analysis of a Polypropylene Capacitor for Resonant Tuning Networks in WPT Applications

This study presents the thermal analysis of polypropylene capacitors operating at high frequencies and high currents used in the resonant tuning networks for electric vehicle (EV) wireless charging systems. The thermal analysis and performance results of the CELEM polypropylene capacitors are presented at different currents for a polyphase wireless power transfer (WPT) system while the study can be expanded to other capacitors by different manufacturers. Thermal equivalent circuit of the capacitor is derived and required cold plate size and design calculations are conducted by modeling the polypropene capacitor in MATLAB. The polypropene capacitor's finite element analysis (FEA) model is developed using COMSOL Multiphysics simulation software, and the thermal characteristics and analysis of the capacitor are presented in this paper.

Aktas, Ahmet↗
DOE OSTI2024

Compositional zoning of the Otowi Member of the Bandelier Tuff, Valles caldera, New Mexico, USA

The Otowi Member of the Bandelier Tuff erupted at ca. 1.60 Ma from the Valles caldera (New Mexico, USA). It consists of as much as 400 km 3 (dense rock equivalent) of strongly differentiated high-silica rhyolite and shows systematic upward variations in crystallinity, mineral chemistry, and trace element concentrations through its thickness, but the major element composition is almost constant and is near the low-pressure granite minimum. Incompatible trace elements in whole pumice fragments and glasses show well-correlated linear covariations. Upward zoning to lower abundances of incompatible trace elements is accompanied by development of overgrowths on quartz and alkali feldspar, although earlier-formed interiors of quartz and feldspar have near-constant compositions throughout the tuff, modified by cation diffusion in the case of feldspar. Melt inclusions in remnant quartz cores show diverse Pb isotope ratios, pointing to a wide range of distinct protoliths that contributed rhyolitic melt to the Otowi magma. Mineral thermometers suggest a modest temperature gradient through the melt body, perhaps of 40 °C, at the time of eruption. Chemical, textural, and mineralogical variations and volume-composition relations through the tuff are consistent with an origin for zoning by melting of a high-crystallinity cumulate layer beneath cognate supernatant liquid to produce denser, remobilized liquid of accumulative composition (i.e., the “modified mush model”). Melting may have occurred in several episodes. The latest of these episodes, probably thousands of years prior to eruption, introduced new rhyolitic liquid into the system and was associated with a thermal excursion, recorded in core compositions of pyroxene, during which much of the earlier crystal mass was dissolved. This left inherited cores and interiors of accumulated quartz and feldspar mantled with new growth having less-evolved compositions (higher Ti, Sr, and Ba). Changing solubility of zircon during cumulate melting produced a reversal of Zr concentrations. There is no clear petrologic evidence of a recharge eruption trigger; nonetheless, compositional zoning resulted mainly from repeated recharge-induced remobilization of quartz-feldspar cumulate. The Otowi system was built, evolved, and modified by several events over the course of a few hundred thousand years.

Wolff, J. A.↗
DOE OSTI2024

AutoCSM

A template system-of-systems modeling approach for automating the development, deployment, and integration of cooling system models (CSMs) for supercomputing facilities within the ExaDigiT framework. AutoCSM is a Python-based framework to assist in CSM developers in accelerating the creation and deployment of system-level thermal-hydraulic CSMs. The intention is for this tool specifically to help standardize digital twin workflows for ExaDigiT. However, this tool can be used independent of ExaDigiT (and even other systems besides CSMs).

Greenwood, MichaelScott [Oak Ridge National Labora↗
DOE OSTI2023

Physics-constrained fault diagnosis framework for monitoring a multi-component thermal hydraulic system

A method for diagnosing faults includes receiving a system description of a thermal hydraulic system, the system description indicating a plurality of components and sensors. The method also includes constructing, based on physical conservations laws and using the system description, a plurality of physics-based models for the plurality of components, each of the plurality of physics-based models including unknown parameters. The method further includes receiving historical measurements and calibrating the physics-based models by calculating the unknown parameters of each of the physics-based models using the historical measurements to produce calibrated models. The method also includes receiving sensor measurements of the sensors, and calculating residuals corresponding to differences between measurements predicted by the plurality of calibrated models and the sensor measurements. The method also includes determining, based on the calculated residuals, a fault of a component or a sensor, and generating an alert indicating the fault.

Nguyen, Tat Nghia↗
DOE OSTI2023

Design of a Geothermal Power Plant With Solar Thermal Topping Cycle: Preprint

Geothermal power plants are a reliable source of low-carbon power generation. However, modern electricity markets comprise relatively large proportions of variable renewable energy generation that may require power plants to dispatch energy flexibly. The power output, efficiency, and dispatch flexibility of a geothermal plant can be enhanced by integrating solar thermal energy into the system, as well as possibly compensating against ambient temperature variations. Concentrating Solar Thermal (CST) can generate temperatures much higher than conventional geothermal systems. Using a solar topping cycle is one way to efficiently convert high-temperature solar heat to electricity while also adding lower-temperature heat to the geothermal power cycle, thereby increasing its power output and possibly its efficiency. A hybrid power cycle design is proposed and is simulated using SimTech IPSEpro process modelling software. The design configuration depends upon the expected temperature of the geothermal resource and the quantity of solar heat added at the design point. These design considerations are described and expected performance is calculated. The solar heat addition varies throughout the day and year, therefore off-design models are necessary to assess the impact of solar availability (and ambient temperature) on the power plant performance. Off-design models are developed and combined with hourly weather data to facilitate an evaluation of annual system performance.

concentrating solar power↗
DOE OSTI2024

Hydrodynamic Modeling of Stratification and Mixing in a Shallow, Tropical Floodplain Lake

Abstract Floodplain lakes are widespread and ecologically important throughout tropical river systems, however data are rare that describe how temporal variations in hydrological, meteorological and optical conditions moderate stratification and mixing in these shallow lakes. Using time series measurements of meteorology and water‐column temperatures from 17 several day campaigns spanning two hydrological years in a representative Amazon floodplain lake, we calculated surface energy fluxes and thermal stratification, and applied and evaluated a 3‐dimensional hydrodynamic model. The model successfully simulated diel cycles in thermal structure characterized by buoyancy frequency, depth of the actively mixing layer, and other terms associated with the surface energy budget. Diurnal heating with strong stratification and nocturnal mixing were common; despite considerable heat loss at night, the strong stratification during the day meant that mixing only infrequently extended to the bottom at night. Simulations indicated that the diurnal thermocline up and downwelled creating lake‐wide differences in near‐surface temperatures and mixing depths. Infrequent full mixing creates conditions conducive to anoxia in these shallow lakes given their warm temperatures.

Environmental Sciences & Ecology↗
DOE OSTI2023

Automotive Silicon Carbide Power Module Cooling With A Novel Modular Manifold And Embedded Heat Sink

The next generation of integrated power electronics packages will implement wide-bandgap devices with ultrahigh device heat fluxes. Although jet impingement has received attention for power electronics thermal management, it is not used in commercial electric vehicles (EVs) because of the associated pressure drop and reliability concerns. Here, in this paper, we present a modular thermal management system designed for automotive power electronics. The system achieves superior thermal performance to benchmarked EVs, while adhering to reliability standards and with low pumping power. The system utilizes a low-cost and lightweight plastic manifold to generate jets over an optimized heat sink, which is embedded in the direct-bonded-copper (DBC) substrate. The embedded heat sink concept leverages additive manufacturing to add elliptical pin fins to the DBC substrate. The heat sink geometry is optimized for submerged jet impingement using a unit-cell model and an exhaustive search algorithm. The model predictions are validated using unit-cell experiments. A full-scale power module model is then used to compare the DBC-embedded heat sink against direct DBC cooling and baseplate-integrated heat sinks for single-sided (SS) and double-sided (DS) cooling concepts. Using the SS and DS DBC-embedded cooling concepts, the models predict a thermal resistance that represents a reduction of 75% and 85% compared to the 2015 BMW i3, respectively, for the same water-ethylene glycol inverter flow rate. We have shown that an inverter with a 100-kilo-Watt-per-liter power density is achievable with the proposed design.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Classical signatures of quenched and thermal disorder in the dynamics of correlated spin systems

Neutron scattering is frequently used to look for evidence of features indicative of quantum-entangled phases of matter such as continua from fractionalisation or quantised excitations. However, the non-specificity of these features and difficulty of both fully quantum treatments and semiclassical models of disorder, make the diagnosis of such states problematic. Here, we demonstrate the feasibility of semiclassical treatments of disordered systems for supercells of ~10 000 spins. By examining a number of classically disordered models we show the presence of quantised excitations, broad continua and anomalous damping originating from quenched disorder or large classical degeneracies.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Unifying femtosecond and picosecond single-pulse magnetic switching in Gd-Fe-Co

Many questions are still open regarding the physical mechanisms behind the magnetic switching in Gd-Fe-Co alloys by single optical pulses. Phenomenological models suggest a femtosecond scale exchange relaxation between sublattice magnetization as the driving mechanism for switching. The recent observation of thermally induced switching in Gd-Fe-Co by using both several picosecond optical laser pulse as well as electric current pulses has questioned this previous understanding. This has raised the question of whether or not the same switching mechanics are acting at the femtosecond and picosecond scales. In this work, we aim at filling this gap in the understanding of the switching mechanisms behind thermal single-pulse switching. To that end, we have studied experimentally thermal single-pulse switching in Gd-Fe-Co alloys, for a wide range of system parameters, such as composition, laser power, and pulse duration. We provide a quantitative description of the switching dynamics using atomistic spin dynamics methods with excellent agreement between the model and our experiments across a wide range of parameters and timescales, ranging from femtoseconds to picoseconds. Furthermore, we find distinct element-specific damping parameters as a key ingredient for switching with long picosecond pulses and argue that switching with pulse durations as long as 15 ps is possible due to a low damping constant of Gd. Our findings can be easily extended to speed up dynamics in other contexts where ferrimagnetic Gd-Fe-Co alloys have been already demonstrated to show fast and energy-efficient processes, e.g., domain-wall motion in a track and spin-orbit torque switching in spintronics devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

One-step production of renewable adipic acid esters from mucic acid over an Ir–ReO x /C catalyst with low Ir loading

The production of adipic acid, a large-volume platform chemical, from biomass represents a petroleum-free route to manufacturing nylon and other polymers more sustainably. In this study, a one-step conversion of bioderived mucic acid to adipates is reported over a heterogeneous, bifunctional Ir–ReO x /C catalyst via deoxydehydration (DODH) and catalytic transfer hydrogenation (CTH) using isopropanol as a green solvent and reductant. With very low Ir content (0.05 wt%), the catalyst provides a low-cost option for the tandem DODH–CTH process, while still enabling reuse and regeneration for at least five cycles. The reported catalytic system generates adipates in good yield (63%) without the need for additives (acids or halogens) or high pressure H 2 . Through model reactions and spectroscopic analyses, a bifunctional DODH–CTH mechanism was supported, with the Re VI/IV redox pair as the active DODH species and Ir0 saturating the dialkene intermediate via CTH. Thermal treatment in catalyst preparation and regeneration was optimized to limit metal leaching while maintaining good activity. Further, the present work outlines a one-step deoxygenation of an aldaric acid substrate from biomass using a bimetallic catalyst of low noble metal content to an important monomer for the polymer industry. This system facilitates use in future syntheses of sustainable chemicals from renewable oxygenates where the catalyst is reusable and affordable.

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

Preliminary Plan to Inform Testing of a Heat Exchanger Test Article

This report presents a preliminary plan to guide the qualification testing of advanced heat exchanger (HX) components for nuclear-to-industrial heat transfer applications. The objective is to establish a defensible, physics-based methodology that integrates computational modeling, targeted experimentation, and in-service inspection considerations to demonstrate component performance and reliability under representative reactor conditions. The analysis identifies Sodium-cooled Fast Reactor (SFR) and High-Temperature Gas-cooled Reactor (HTGR) systems as reference configurations in terms of temperature, pressure, and chemical environment. Within these operating envelopes, dominant degradation mechanisms— including creep–fatigue interaction, flow-induced vibration, corrosion, and diffusion-bond deterioration—were evaluated to define test requirements. A comprehensive computationalexperimental framework is proposed to support life prediction and qualification activities. The framework couples high-fidelity structural-mechanics, thermal-hydraulic, and fluid-structure interaction models with accelerated degradation testing to produce a traceable linkage between microstructural evolution, mechanical performance, and remaining useful life (RUL). The approach adheres to established Verification, Validation, and Uncertainty Quantification (VVUQ) standards (ASME V&V 10/20; NUREG-2152) and incorporates a digital-twin architecture for continuous model refinement through data assimilation. The plan further outlines testing methodologies, including pre-test analyses, test-loop design parameters, and sensor placement strategies that maximize information yield while maintaining mechanistic fidelity. Complementary sections describe in-service inspection (ISI), on-line monitoring (OLM), and structural-health-monitoring (SHM) techniques applicable to compact HX geometries typical of advanced reactors. Collectively, these activities establish the technical foundation for demonstrating 40-60-year equivalent service life of advanced heat exchangers in support of the U.S. Department of Energy’s Advanced Reactor and Integrated Energy Systems programs. The forthcoming phase will execute the defined pre-test analyses, initiate hardware fabrication, and implement the integrated testing campaign to validate the proposed qualification methodology.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗