Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal dynamics”

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

Evaporation of Water Nanodroplets on Heated Surfaces: Does Nano Matter?

While experiments and continuum models have provided a relatively good understanding of the evaporation of macroscopic water droplets, elucidating how sessile nanodroplets evaporate is an open question critical for advancing nanotechnological applications where nanodroplets can play an essential role. In this study, using molecular dynamics simulations, we find that evaporating nanodroplets, in contrast to their macroscopic counterparts, are not always in thermal equilibrium with the substrate and that the vapor concentration on the nanodroplet surface does not reach a steady state. As a result, the evaporative behavior of nanodroplets is significantly different. Regardless of hydrophobicity, nanodroplets do not follow conventional evaporation modes but instead exhibit dynamic wetting behavior characterized by huge, non-equilibrium, isovolumetric fluctuations in the contact angle and contact radius. For hydrophilic nanodroplets, the evaporation rate, controlled by the vapor concentration, decays exponentially over time. Hydrophobic nanodroplets follow stretched exponential kinetics arising from the slower thermalization with the substrate. The evaporative half-lifetime of the nanodroplets is directly related to the thermalization time scale and therefore increases monotonically with the hydrophobicity of the substrate. Finally, the evaporative flux profile along the nanodroplet surface is highly nonuniform but does not diverge at the contact line as the macroscopic continuum models predict.

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

High-temperature ion-thermal behavior from average-atom calculations

Atom-in-jellium calculations of the Einstein frequency were used to calculate the mean displacement of an ion over a wide range of compression and temperature. Expressed as a fraction of the Wigner-Seitz radius, the displacement is a measure of the asymptotic freedom of the ion at high temperature, and thus of the change in heat capacity from six to three quadratic degrees of freedom per atom. A functional form for free energy was proposed based on the Maxwell-Boltzmann distribution as a correction to the Debye free energy, with a single free parameter representing the effective density of potential modes to be saturated. Finally, this parameter was investigated using molecular dynamics simulations, and found to be ~0.2 per atom. In this way, the ion-thermal contribution can be calculated for a wide-range equation of state (EOS) without requiring a large number of molecular dynamics simulations. Example calculations were performed for carbon, including the sensitivity of key EOS loci to ionic freedom.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

An Efficient 1-D Thermal Stratification Model for Pool-Type Sodium-Cooled Fast Reactors

Investigating thermal stratification in the upper plenum of a sodium fast reactor (SFR) is presently a technology gap in SFR safety analysis. Understanding thermal stratification will promote safe operation of the SFR before its commercial deployment. Stratified layers of liquid sodium with a large vertical temperature gradient could be established in the upper plenum of an SFR during a down-power or a loss-of-flow transient. These stratified layers are unstable and could result in uncertainties for the core safety of an SFR. In order to predict the occurrence of the thermal stratification efficiently, we developed a one-dimensional (1-D) transport model to estimate the temperature profile of the ambient fluid in the upper plenum. This model demands much less computational effort than computational fluid dynamics (CFD) codes and provides calculations with higher fidelity than historical system-level codes. Two flow conditions were considered separately in the current study depending on if in-vessel components are presented in the upper plenum. For the condition where in-vessel components, specifically the upper internal structure, are presented, we assumed that the impinging sodium was evenly dispersed in the ambient fluid within the distance between the bottom of the in-vessel component and the jet inlet surface. For the condition where no in-vessel components are presented, we assumed that the impinging sodium was evenly dispersed in the ambient fluid within the jet length, which was determined through data-driven trainings. The newly developed 1-D model showed similar performance with the CFD model in both cases. However, due to the assumption of flat profiles of the impinging jet axial dispersion rate, nonnegligible discrepancies between the 1-D prediction and the measured data were observed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Validation of Pronghorn for Natural-Circulation Molten Salt Loops

This paper presents the development and validation of a high-fidelity thermal-hydraulic model of a molten salt natural circulation flow loop, designed for integration within a digital twin framework. The study evaluates the performance of Idaho National Laboratory’s Pronghorn against experimental data from Texas A&M University Molten Salt Flow Loop (MSFL) four Hitec-salt test benchmark data. Natural circulation of high-Prandtl-number fluids exhibits complex, counter-intuitive flow patterns that make pointwise thermocouple readings unreliable. Experimental work at TAMU’s MSFL provides benchmark data, including flow visualization at a test-section and centerline steady-state temperature measurements along the loop. Validation includes four single-phase natural circulation test cases with Hitec salt. Key metrics include flow profile agreement and steady-state temperature accuracy. Pronghorn results for two-dimensional single-phase agree qualitatively with the experimental flow profile. This paper illustrates the importance of Computational Fluid Dynamics (CFD) in elucidating the behavior of high-Prandtl-number thermal-hydraulics, along with how misleading centerline temperature measurements can be. Pronghorn reproduces the axial and radial stratification that makes single thermocouple readings unreliable. Future research will focus on reduced-order modeling techniques to enable rapid simulation suitable for real-time digital twin applications. The validated cases provide a basis for developing reduced-order surrogates aimed at real-time digital-twin applications.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

LC-Opt: Benchmarking Reinforcement Learning and Agentic AI for End-to-End Liquid Cooling Optimization in Data Centers

Liquid cooling is critical for thermal management in high-density data centers with the rising AI workloads. However, machine learning-based controllers are essential to unlock greater energy efficiency and reliability, promoting sustainability. We present LC-Opt, a Sustainable Liquid Cooling (LC) benchmark environment, for reinforcement learning (RL) control strategies in energy-efficient liquid cooling of high-performance computing (HPC) systems. Built on the baseline of a high-fidelity digital twin of Oak Ridge National Lab's Frontier Supercomputer cooling system, LC-Opt provides detailed Modelica-based end-to-end models spanning site-level cooling towers to data center cabinets and server blade groups. RL agents optimize critical thermal controls like liquid supply temperature, flow rate, and granular valve actuation at the IT cabinet level, as well as cooling tower (CT) setpoints through a Gymnasium interface, with dynamic changes in workloads. This environment creates a multi-objective real-time optimization challenge balancing local thermal regulation and global energy efficiency, and also supports additional components like a heat recovery unit (HRU). We benchmark centralized and decentralized multi-agent RL approaches, demonstrate policy distillation into decision and regression trees for interpretable control, and explore LLM-based methods that explain control actions in natural language through an agentic mesh architecture designed to foster user trust and simplify system management. LC-Opt democratizes access to detailed, customizable liquid cooling models, enabling the ML community, operators, and vendors to develop sustainable data center liquid cooling control solutions.

Naug, Avisek [Hewlett Packard Enterprise]↗
DOE OSTI2023

Towards Quantum Computing Phase Diagrams of Gauge Theories with Thermal Pure Quantum States

The phase diagram of strong interactions in nature at finite temperature and chemical potential remains largely theoretically unexplored due to inadequacy of Monte-Carlo–based computational techniques in overcoming a sign problem. Quantum computing offers a sign-problem-free approach, but evaluating thermal expectation values is generally resource intensive on quantum computers. To facilitate thermodynamic studies of gauge theories, we propose a generalization of the thermal-pure-quantum-state formulation of statistical mechanics applied to constrained gauge-theory dynamics, and numerically demonstrate that the phase diagram of a simple low-dimensional gauge theory is robustly determined using this approach, including mapping a chiral phase transition in the model at finite temperature and chemical potential. Quantum algorithms, resource requirements, and algorithmic and hardware error analysis are further discussed to motivate future implementations. Thermal pure quantum states, therefore, may present a suitable candidate for efficient thermal simulations of gauge theories in the era of quantum computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Scanning ultrafast electron microscopy reveals photovoltage dynamics at a deeply buried $p^- \mathrm{Si/SiO_2}$ interface

The understanding and control of charge carrier interactions with defects at buried insulator/semiconductor interfaces is essential for achieving optimum performance in modern electronics. In this work, we report on the use of scanning ultrafast electron microscopy (SUEM) to remotely probe the dynamics of excited carriers at a Si surface buried below a thick thermal oxide. Our measurements illustrate a previously unidentified SUEM contrast mechanism, whereby optical modulation of the space-charge field in the semiconductor modulates the electric field in the thick oxide, thus affecting its secondary electron yield. By analyzing the SUEM contrast as a function of time and laser fluence we demonstrate the diffusion mediated capture of excited carriers by interfacial traps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Temporal and spatial resolution of distal protein motions that activate hydrogen tunneling in soybean lipoxygenase

The enzyme soybean lipoxygenase (SLO) provides a prototype for deep tunneling mechanisms in hydrogen transfer catalysis. This work combines room temperature X-ray studies with extended hydrogen–deuterium exchange experiments to define a catalytically-linked, radiating cone of aliphatic side chains that connects an active site iron center of SLO to the protein–solvent interface. Employing eight variants of SLO that have been appended with a fluorescent probe at the identified surface loop, nanosecond fluorescence Stokes shifts have been measured. We report a remarkable identity of the energies of activation ( E a ) for the Stokes shifts decay rates and the millisecond C–H bond cleavage step that is restricted to side chain mutants within an identified thermal network. These findings implicate a direct coupling of distal protein motions surrounding the exposed fluorescent probe to active site motions controlling catalysis. While the role of dynamics in enzyme function has been predominantly attributed to a distributed protein conformational landscape, the presented data implicate a thermally initiated, cooperative protein reorganization that occurs on a timescale faster than nanosecond and represents the enthalpic barrier to the reaction of SLO.

Science & Technology - Other Topics↗
DOE OSTI2024

Thermally frustrated phase transition at high pressure in B 2-ordered FeV

X-ray diffraction measurements of equiatomic B2-ordered FeV were performed in a diamond-anvil cell at room temperature at several pressure points up to 80 GPa that showed the cubic phase to be stable with no indication of structural phase transitions. Density functional theory at 0 K predicts Fermi surface nesting, an electronic topological transition, and a phonon dynamical instability within the experimentally investigated pressure range. Nevertheless, the instability is absent in phonon dispersion curves extracted from ab initio molecular dynamics simulations below the critical volume at temperatures as low as 10 K, indicating that thermal atomic displacements can frustrate the phase transition by renormalizing the phonon dispersion curves. Ferrimagnetism is critical for the stability of the cubic phase at low temperature, but thermal atomic displacements are enough to support the structure at and above the Néel temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

In-Situ Sensing and Dynamics Predictions for Electrothermally-Actuated Soft Robot Limbs

Untethered soft robots that locomote using electrothermally-responsive materials like shape memory alloy (SMA) face challenging design constraints for sensing actuator states. At the same time, modeling of actuator behaviors faces steep challenges, even with available sensor data, due to complex electrical-thermal-mechanical interactions and hysteresis. This article proposes a framework for in-situ sensing and dynamics modeling of actuator states, particularly temperature of SMA wires, which is used to predict robot motions. A planar soft limb is developed, actuated by a pair of SMA coils, that includes compact and robust sensors for temperature and angular deflection. Data from these sensors are used to train a neural network-based on the long short-term memory (LSTM) architecture to model both unidirectional (single SMA) and bidirectional (both SMAs) motion. Predictions from the model demonstrate that data from the temperature sensor, combined with control inputs, allow for dynamics predictions over extraordinarily long open-loop timescales (10 min) with little drift. Prediction errors are on the order of the soft deflection sensor’s accuracy. This architecture allows for compact designs of electrothermally-actuated soft robots that include sensing sufficient for motion predictions, helping to bring these robots into practical application.

Sabelhaus, Andrew P.↗
DOE OSTI2022

FLASSH 1.0: Thermal Scattering Law Evaluation and Cross Section Generation

The Full Law Analysis Scattering System Hub ( FLASSH ) is an advanced code which evaluates the thermal scattering law (i.e. TSL, S(α,β)) for thermal scattering cross sections and resonance Doppler broadening. The ability to accurately capture these two key cross section features is dependent on accurate, high fidelity TSL evaluations. FLASSH 1.0 provides advanced physics capabilities resulting in an improved, generalized TSL to most accurately represent the lattice dynamics within any material. This improved TSL will allow for consistent analysis in both the thermal and epithermal energy ranges. The features for TSL analysis are packaged within the FLASSH GUI for easy user interface along with data output in many file formats including ENDF File 7 and ACE files.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Harnessing Ocean Thermal Gradients Using Thermoelectric Based Submersibles for Ocean Power Applications

The urgent need for energy solutions in marine environments has accelerated the development of innovative technologies capable of leveraging natural resources for power generation. This study introduces a buoyancy-driven submersible system designed to harness ocean thermal gradients using thermoelectric generators (TEGs) and phase change materials (PCMs). The technology aims to provide autonomous power to offshore aquaculture farms, unmanned underwater vehicles (UUVs), offshore platform illumination, and ocean sensors, significantly reducing dependence on fossil fuels. Ocean thermal gradients, especially prevalent in mid-latitude regions, exhibit temperature differences between surface and deep waters ranging from 7 degrees Celsius to 30 degrees Celsius depending on seasonal variations. The proposed submersible technology utilizes TEGs to convert thermal energy from these gradients into electrical power, generating between 0.2 and 0.5 watts, while PCMs are employed to store and regulate this energy, ensuring a stable and continuous power supply. The buoyancy-driven mechanism of the submersible enhances its capability to navigate through varying depths, optimizing its exposure to different thermal gradients and maximizing energy harvesting. The performance of this submersible system is analyzed through detailed thermodynamic assessments and computational fluid dynamics (CFD) modeling focused on heat transfer. These analyzes consider real-world ocean temperature profiles and seek to refine the interaction between TEGs and PCMs to optimize energy extraction. The evaluation encompasses several key performance metrics, including power output and energy efficiency. Results confirm the potential of this innovative technology to provide a continuous and reliable power source for marine applications. By demonstrating the feasibility of using ocean thermal gradients for energy generation, this study contributes to the broader efforts of innovation in energy technologies for harsh, remote marine environments. The implementation of such promises is significant advancements in the autonomy of marine operations. The ongoing research will further investigate scalability ensuring that the system can be effectively adapted to various marine settings and operational demands.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Experimental observation of localized interfacial phonon modes

Interfaces impede heat flow in micro/nanostructured systems. Conventional theories for interfacial thermal transport were derived based on bulk phonon properties of the materials making up the interface without explicitly considering the atomistic interfacial details, which are found critical to correctly describing thermal boundary conductance. Recent theoretical studies predicted the existence of localized phonon modes at the interface which can play an important role in understanding interfacial thermal transport. However, experimental validation is still lacking. Through a combination of Raman spectroscopy and high-energy-resolution electron energy-loss spectroscopy in a scanning transmission electron microscope, we report the experimental observation of localized interfacial phonon modes at ~12 THz at a high-quality epitaxial Si-Ge interface. These modes are further confirmed using molecular dynamics simulations with a high-fidelity neural network interatomic potential, which also yield thermal boundary conductance agreeing well with that measured in time-domain thermoreflectance experiments. Simulations find that the interfacial phonon modes have an obvious contribution to the total thermal boundary conductance. Our findings significantly contribute to the understanding of interfacial thermal transport physics and have impact on engineering thermal boundary conductance at interfaces in applications such as electronics thermal management and thermoelectric energy conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗
DOE OSTI2020

Dimensional crossover in spin Hall oscillators

Auto-oscillations of magnetization driven by direct spin current have been previously observed in multiple quasi-zero-dimensional (0D) ferromagnetic systems such as nanomagnets and nanocontacts. Recently, it was shown that pure spin Hall current can excite coherent auto-oscillatory dynamics in quasi-one-dimensional (1D) ferromagnetic nanowires but not in quasi-two-dimensional (2D) ferromagnetic films. Here we study the 1D to 2D dimensional crossover of current-driven magnetization dynamics in wire-based Pt/Ni 80 Fe 20 bilayer spin Hall oscillators via varying the wire width. We find that increasing the wire width results in an increase of the number of excited auto-oscillatory modes accompanied by a decrease of the amplitude and coherence of each mode. We also observe a crossover from a hard to a soft onset of the auto-oscillations with increasing wire width. The amplitude of auto-oscillations rapidly decreases with increasing temperature suggesting that interactions of the phase-coherent auto-oscillatory modes with incoherent thermal magnons play an important role in suppression of the auto-oscillatory dynamics. Our measurements in this work set the upper limit on the dimensions of an individual spin Hall oscillator and elucidate the mechanisms leading to suppression of coherent auto-oscillations with increasing oscillator size.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Local Sn Dipolar-Character Displacements behind the Low Thermal Conductivity in SnSe Thermoelectric

The local atomic structure of SnSe was characterized across its orthorhmbic-to-orthorhombic structural phase transition using x-ray pair distribution function analysis. Substantial Sn displacements with a dipolar character persist in the high symmetry high temperature phase, albeit with a symmetry different from that of the ordered displacements below the transition. The analysis implies that the transition is neither order-disorder nor displacive, but rather a complex crossover. Robust ferro-coupled SnSe intra-layer distortions suggest a ferroelectric-like instability as the driving force. Furthermore, these local symmetry-lowering Sn displacements are likely integral to the ultra-low lattice thermal conductivity mechanism in SnSe.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Particle Acceleration by Magnetic Reconnection in Geospace

Particles are accelerated to very high, non-thermal energies during explosive energy-release phenomena in space, solar, and astrophysical plasma environments. While it has been established that magnetic reconnection plays an important role in the dynamics of Earth’s magnetosphere, it remains unclear how magnetic reconnection can further explain particle acceleration to non-thermal energies. Here we review recent progress in our understanding of particle acceleration by magnetic reconnection in Earth’s magnetosphere. With improved resolutions, recent spacecraft missions have enabled detailed studies of particle acceleration at various structures such as the diffusion region, separatrix, jets, magnetic islands (flux ropes), and dipolarization front. With the guiding-center approximation of particle motion, many studies have discussed the relative importance of the parallel electric field as well as the Fermi and betatron effects. However, in order to fully understand the particle acceleration mechanism and further compare with particle acceleration in solar and astrophysical plasma environments, there is a need for further investigation of, for example, energy partition and the precise role of turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A methodology for domain overlapping coupling of thermal-hydraulic systems

Multi-scale coupling has increasingly drawn attention as a promising approach for modeling thermal systems. Traditional system codes provide fast and robust predictions at the plant scale, while high-fidelity computational fluid dynamics (CFD)-based tools resolve localized flow and heat transfer phenomena with greater accuracy. By combining these complementary strengths, co-simulations enable multi-scale analysis that would otherwise be computationally prohibitive for a standalone CFD code. Here, this work introduces a robust and problem-agnostic domain overlapping (DO) coupling between the system thermal-hydraulic (STH) code System Analysis Module (SAM) and the coarse-mesh CFD code Pronghorn. Both applications belong to the Comprehensive Reactor Analysis Bundle (BlueCRAB) code suite, a code suite in active development at the Idaho National Laboratory (INL), tailored for multi-physics analysis of advanced reactors. Unlike previous approaches, BlueCRAB supports an agnostic interface between codes based on different fidelity, while its coupling formulation can address arbitrary flow geometries with multiple inlets and outlets in coupled components. The implemented method leads to consistent pressure drops, enthalpies, and scalar concentrations between coupled SAM and Pronghorn simulations. The methodology is demonstrated through two verification tests, which ensure the numerical consistency and conservation across the codes, and through one validation test against experimental data. The proposed problems explore different physical aspects inherent to thermal systems, with particular attention given to nuclear reactor analysis. These include buoyancy-driven flows, complex flow patterns, and setups with multiple inlets and outlets, representing challenges in advanced reactor applications.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗