Search NASA⌕ Search

SEARCH · Search NASA

Results for “THERMOPHYSICS”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Liquid–Vapor Phase Equilibrium in Molten Aluminum Chloride (AlCl 3 ) Enabled by Machine Learning Interatomic Potentials

Molten salts are promising candidates in numerous clean energy applications, where knowledge of thermophysical properties and vapor pressure across their operating temperature ranges is critical for safe operations. Due to challenges in evaluating these properties using experimental methods, fast and scalable molecular simulations are essential to complement the experimental data. In this study, we developed machine learning interatomic potentials (MLIP) to study the AlCl 3 molten salt across varied thermodynamic conditions (T = 473–613 K and P = 2.7–23.4 bar), which allowed us to predict temperature-surface tension correlations and liquid–vapor phase diagram from direct simulations of two-phase coexistence in this molten salt. Two MLIP architectures, a Kernel-based potential and neural network interatomic potential (NNIP), were considered to benchmark their performance for AlCl 3 molten salt using experimental structure and density values. The NNIP potential employed in two-phase equilibrium simulations yields the critical temperature and critical density of AlCl 3 that are within 10 K (∼3%) and 0.03 g/cm 3 (∼7%) of the reported experimental values. An accurate correlation between temperature and viscosities is obtained as well. In doing so, we report that the inclusion of low-density configurations in their training is critical to more accurately represent the AlCl 3 system across a wide phase-space. The MLIP trained using PBE-D3 functional in the ab initio molecular dynamics (AIMD) simulations (120 atoms) also showed close agreement with experimentally determined molten salt structure comprising Al 2 Cl 6 dimers, as validated using Raman spectra and neutron structure factor. Furthermore, the PBE-D3 as well as its trained MLIP showed better liquid density and temperature correlation for AlCl 3 system when compared to several other density functionals explored in this work. Overall, the demonstrated approach to predict temperature correlations for liquid and vapor densities in this study can be employed to screen nuclear reactors-relevant compositions, helping to mitigate safety concerns.

Ab initio molecular dynamics↗
DOE OSTI2026

Ionic Liquid-Enhanced Interfaces to Boost Reactive C O2 Capture

The addition of ionic liquids (ILs) to a mixture containing a molecular solvent and other ionic species can induce the heterogeneous redistribution of cations and anions at the gas–liquid interface. This nonuniform redistribution of cations and anions driven by the differences in the solvophilicity of ions can improve the thermophysical and interfacial properties of such mixtures, creating a local chemical environment that is conducive to some reactions. In this work, ILs are added to a mixture of potassium hydroxide (KOH) and ethylene glycol (EG), used as a reactive absorbent and electrolyte in the migration-assisted moisture-gradient (MAMG) process for CO 2 capture. Molecular dynamics (MD) simulations are employed to probe into the effects of complex ion–ion and ion–solvent interactions and to examine the chemical composition at the gas–liquid interface. A total of 12 systems are investigated using molecular simulations to identify trends in the performance of IL additives based on the choice of cation, anion, and IL concentration. The cation effects are studied using IL additives based on 1-ethyl-3-methylimidazolium ([EMIM] + ) and 1-butyl-3-methylimidazolium ([BMIM] + ), while the impact of anions is examined using additives based on dicyanamide [DCA] − , triflate [TfO] − , bistriflimide [NTf 2 ] − , and hexafluorophosphate [PF 6 ] − anions, respectively. The influence of the IL concentration is also evaluated at molar concentrations between 1% and 4%. The simulation results indicate that the use of IL additives can affect the physical CO 2 solubility, surface tension, and the localization of CO 2 around the [OH] − ions at the gas–liquid interface. It is also evident that the choice of cations, anions, and IL concentration determines the extent to which the IL additives impact the local physicochemical properties. Physical dissolution, diffusive transport, and interaction with [OH] − are critical intermediate steps toward reactive CO 2 capture using a liquid absorbent. Hence, the improvement in one or more of these properties, aided by IL additives, is expected to improve the overall CO 2 capture performance. Experiments reaffirmed the impact of IL additives on CO 2 capture performance and the sensitivity to the choice of the cation, anion, and concentration of the IL additive.

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

Surface Segregation of Liquid Metal Plasma-Facing Component Alloys: A ReaxFF Investigation

Engineering liquid metal alloys offers a transformative pathway for plasma-facing components by enabling chemically tailored surfaces that can simultaneously optimize plasma-material interactions, reduce divertor heat flux, and enhance core plasma confinement, thereby advancing the commercial viability of nuclear fusion power plants. This study, employing an atomistic simulation approach, provides direct evidence that incorporating nonmetal surface-active agents (such as O and H, or their combination) enables strong surface segregation. This capability makes tin−aluminum (Sn−Al) and tin−lithium (Sn−Li) alloys, with suitable compositions, good candidates for PFC applications. Specifically, the presence of low-Z solutes (Li, Al) leads to preferential surface enrichment, which imparts low-Z sputtering characteristics, while the Sn solvent maintains thermophysical stability. To systematically examine this behavior, we developed a ReaxFF force field spanning the full Sn/Al/Li/O/H chemistry, validated it against formation energies and elastic constants, and applied it in reactive molecular dynamics simulations at fusion-relevant temperatures. We also introduced an overlapbased segregation index that captures interfacial compositional separation directly from atomistic density distributions. Here, this metric reveals a clear hierarchy of segregation regimes and provides a unified view across all systems studied. Together, these findings establish a mechanistic link between nonmetal chemistry and interfacial structure, providing a predictive framework for designing self-adaptive, low-sputtering liquid metal alloys for fusion applications.

Alloys↗
DOE OSTI2026

Elucidating the Microscale Behavior and Phase Separation Kinetics of Thermally Responsive Ionic Liquid–Water Mixtures

Thermally responsive ionic liquids (ILs) exhibit liquid-liquid phase separation into a water-rich (WR) and ionic-liquid-rich (ILR) phase when heated above a lower critical solution temperature (LCST). This phase behavior has been leveraged for applications ranging from forward osmosis (FO) desalination, where the IL acts as a draw solute, to refrigeration and dehumidification cycles, where the IL acts as a liquid desiccant. While significant effort has been devoted to characterizing the thermodynamic and thermophysical properties of LCST ILs, their phase separation kinetics have not been investigated. In this work, we describe the macroscale phase separation kinetics (phase separation time) by gleaning insight into the microscale colloidal behavior of aqueous mixtures of four different materials, P 4444 TFA (tetrabutylphosphonium-2,4-trifluoroacetate), P 4444 DMBS (tetrabutylphosphonium-2,4-dimethyl-benzenesulfonate), N 4444 Sal (tetrabutylammonium salicylate), and P 4444 Sal (tetrabutylphosphonium salicylate) as a function of IL concentration at a separation temperature of 70 °C. We report the discontinuous microscale size distributions for each material and correlate their theoretical settling velocities to experimental phase separation times. The results indicate that a simple Stokes' law model can predict the phase separation time within reasonable accuracy. Overall, this work lays the foundation for understanding the micro- to macroscale phase separation behavior and kinetics of LCST ILs for various water-energy applications.

LCST↗
DOE OSTI2024

Pore‐Scale Modeling of Carbon Dioxide and Hydrogen Transport During Geologic Gas Storage

Abstract Geologic storage of CO 2 and H 2 are climate‐positive techniques for meeting the energy transition. While similar formations could be considered for both gases, the flow dynamics could differ due to differences in their thermophysical properties. We conduct a rigorous pore‐scale study of water/CO 2 and water/H 2 systems at relevant reservoir conditions in a Bentheimer rock sample using the lattice Boltzmann method to quantify the effects of capillary, viscous, inertial, and wetting forces during gas invasion. At similar conditions, H 2 invasion is weaker compared to CO 2 due to unfavorable viscosity ratios. Increasing flow rate, however, increases the breakthrough saturation for both gas systems in the range of capillary numbers studied. At isolated conditions of flow rate, viscosity ratio, and wettability, local inertial effects are found to be critical and show consistent increase in the invaded gas saturation. The effect of inertial forces persits for both gases across all field conditions tested.

58 GEOSCIENCES↗
DOE OSTI2026

Passive Microwave Radiometry and Active Radar Sounding as Complementary Tools for Geophysical Investigations of Icy Ocean Worlds

Juno Microwave Radiometer (MWR) observations of Europa and Ganymede offer critical insights into the icy shells of these moons ahead of NASA's Europa Clipper and ESA's JUpiter ICy moons Explorer (JUICE) missions. Both missions are equipped with active radar sounders designed to address key unknowns such as ice shell thickness, thermal state, and composition. In this study, we explore how passive microwave radiometry and active radar sounding can constrain ice shell properties, focusing on Europa. Using modeled microwave brightness temperature observations at 0.6 and 1.2 GHz alongside simulated radar attenuation rate observations, we show that each instrument can independently produce robust ice shell thickness constraints under idealized conditions. We then relax these assumptions, quantifying how uncertainties from non-ideal properties—including convective layers, freezing-point depression, and chloride-doped ice—affect thickness estimates. Finally, we demonstrate how combining observations from these complementary techniques breaks degeneracies between ice shell properties, enabling more robust constraints than either method alone. This approach will maximize the science return of Europa Clipper and JUICE, advancing our understanding of the thermophysical structure and habitability of icy ocean worlds.

58 GEOSCIENCES↗
DOE OSTI2025

Nuclear quantum effects in molecular liquids across chemical space

Abstract Nuclear quantum effects (NQEs) influence many physical and chemical phenomena, particularly those involving light atoms or occurring at low temperatures. However, their impact has been carefully quantified in few systems-like water-and is rarely considered more broadly. Here we use path-integral molecular dynamics to systematically investigate NQEs on thermophysical properties of 92 organic liquids at ambient conditions. Depending on chemical constitution, we find substantial impact across thermal expansivity, compressibility, dielectric constant, enthalpy of vaporization, and notably molar volume, which shows consistent, positive quantum-classical differences up to 5%; similar, less pronounced trends manifest as isotope effects from deuteration. Using data-driven analysis, we identify three features-molar mass, classical hydrogen density, and classical thermal expansivity-that accurately predict NQEs and facilitate understanding of how characteristics like branching and heteroatom content influence behavior. This work highlights the broad relevance of NQEs in molecular liquids, while also providing a conceptual and practical framework to anticipate their impact.

Science & Technology - Other Topics↗
DOE OSTI2025

Attention-based functional-group coarse-graining: a deep learning framework for molecular prediction and design

Machine learning (ML) offers considerable promise for the design of new molecules and materials. In real-world applications, the design problem is often domain-specific, and suffers from insufficient data, particularly labeled data, for ML training. In this study, we report a data-efficient, deep-learning framework for molecular discovery that integrates a coarse-grained functional-group representation with a self-attention mechanism to capture intricate chemical interactions. Our approach exploits group-contribution concepts to create a graph-based intermediate representation of molecules, serving as a low-dimensional embedding that substantially reduces the data demands typically required for training. Using a self-attention mechanism to learn the subtle but highly relevant chemical context of functional groups, the method proposed here consistently outperforms existing approaches for predictions of multiple thermophysical properties. In a case study focused on adhesive polymer monomers, we train on a limited dataset comprising only 6,000 unlabeled and 600 labeled monomers. The resulting chemistry prediction model achieves over 92% accuracy in forecasting properties directly from SMILES strings, exceeding the performance of current state-of-the-art techniques. Furthermore, the latent molecular embedding is invertible, enabling the design pipeline to automatically generate new monomers from the learned chemical subspace. We illustrate this functionality by targeting several properties, including high and low glass transition temperatures (Tg), and demonstrate that our model can identify new candidates with values that surpass those in the training set. The ease with which the proposed framework navigates both chemical diversity and data scarcity offers a promising route to accelerate and broaden the search for functional materials.

Han, Ming [Univ. of Chicago, IL (United States)]↗
DOE OSTI2025

Thermodynamics and transport in molten chloride salts and their mixtures

Relationship between thermophysical properties and phonon mean free path. Heat capacity, viscosity, and thermal conductivity in ionic liquids decrease as mean free path decreases and dynamics become less “solid-like” and collect motion diminishes.

Cockrell, C. [Nuclear Futures Institute, Bangor Un↗
DOE OSTI2024

Synthesis and properties of anhydrous rare-earth phosphates, monazite and xenotime: a review

The synthesis methods, crystal structures, and properties of anhydrous monazite and xenotime (REPO 4 ) crystalline materials are summarized within this review. For both monazite and xenotime, currently available Inorganic Crystal Structure Database data were used to study the effects of incorporating different RE cations on the unit cell parameters, cell volumes, densities, and bond lengths. Domains of monazite-type and xenotime-type structures and other AXO4 compounds (A = RE; X = P, As, V) are discussed with respect to cation sizes. Reported chemical and radiation durabilities are summarized. Different synthesis conditions and chemicals used for single crystals and polycrystalline powders, as well as first-principles calculations of the structures and thermophysical properties of these minerals are also provided.

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

Open-source generation of sigma profiles: impact of quantum chemistry and solvation treatment on machine learning performance

The combination of machine learning (ML) models with chemistry-related tasks requires the description of molecular structures in a machine-readable way. The nature of these so-called molecular descriptors has a direct and major impact on the performance of ML models and remains an open problem in the field. Structural descriptors like SMILES strings or molecular graphs lack size-independence and can be memory intensive. Machine-learned descriptors can be of low dimensionality and constant size but lack physical significance and human interpretability. Sigma profiles, which are unnormalized histograms of the surface charge distributions of solvated molecules, combine physical significance with low dimensionality and size-independence, making them a suitable candidate for a universal molecular descriptor. However, their widespread adoption in ML applications requires open access to sigma profile generation, which is currently not available. This work details the development of OpenSPGen – an open-source tool for generating sigma profiles. Also presented are studies on the effect of different settings on the efficacy of the generated sigma profiles at predicting thermophysical material properties when used as inputs to a Gaussian process as a simple surrogate ML model. We find that a higher level of theory does not translate to more accurate results. We also provide further recommendations for sigma profile calculation and use in ML models.

Salih, Fathya Y. M. [University of Notre Dame, IN ↗
DOE OSTI2024

Measuring the density, viscosity, and surface tension of molten titanates using electrostatic levitation in microgravity

Rare earth and barium titanates are useful as ferroelectric, dielectric, and optical materials. Measurements of their thermophysical properties in the liquid state can help guide melt processing technologies for their manufacture and advance understanding of fragile liquids' behavior and glass formation. Here, we report the density, thermal expansion, viscosity, and surface tension of molten BaTi 2 O 5 , BaTi 4 O 9 , and 83TiO 2 -17RE 2 O 3 (RE = La or Nd). Measurements were made using electrostatic levitation and droplet oscillation techniques in microgravity, which provide access to quiescent liquid droplets and deep supercooling of 510–815 K below the equilibrium melting points. Densities were measured over 900–2400 K. Viscosities were similar for all four compositions, increasing from ~10 mPa s near 2100 K to ~30 mPa s near 1750 K. Surface tensions were 450–490 dyn cm -1 for the rare earth titanates and 383–395 dyn cm -1 for the barium titanates; surface tensions of all compositions had small or negligible temperature dependence over 1700–2200 K. In conclusion, for solids recovered after melt quenching, x-ray microtomography revealed the fracture mechanics in crystalline products and minimal internal porosity in glass products, likely arising from entrapped gas bubbles. Internal microstructures were generally similar for products processed either in microgravity or in a terrestrial aerodynamic levitator.

42 ENGINEERING↗
DOE OSTI2024

Sound speed determination in copper shock compressed to 190 GPa

Sound speed measurements in shock compressed solids have long been valuable for the development of equations of state at extreme conditions, shock-induced phase transformations, and a comprehensive characterization of the thermophysical response of high-pressure standards. We present results from plate impact experiments to 190 GPa to determine the longitudinal sound speed in copper—an important high-pressure standard. Surprisingly, the sound speeds determined using the two most common experimental techniques—the front surface impact (FSI) approach and the release wave overtake (RWO) approach—diverge significantly for stresses greater than ~100 GPa. Further analyses, including numerical simulations, show that the FSI experiments provide the correct sound speeds and that fundamental assumptions underlying the RWO method are likely violated due to the complex release response of shock compressed copper. The sound speeds determined using the FSI approach provide for a more accurate high-pressure description of copper in dynamic compression experiments. The present findings are in contrast to the results for shock compressed silver [Wallace et al., Phys. Rev. B 104, 214106 (2021)], where both methods provided consistent sound speed results. Thus, the findings presented here demonstrate the need to experimentally verify the validity of the RWO method on a case-by-case basis. Finally, we note that even at the high stresses in the present work, the copper unloading response shows a time-dependent, quasielastic response often observed in metals at lower stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Density of molten oxides measured in an aero-acoustic levitator

Knowing the thermophysical properties of high-temperature melts can aid the design of melt processes and validate atomic structural models, such as those used in studying glass formation. Property measurements on such melts are challenging, however, due to container-related contamination and heterogeneous nucleation. Containerless processing techniques that employ levitation can be used to avoid these obstacles. In that context, we demonstrate here the application of silhouette imaging to measure the density of molten oxides in an aero-acoustic levitation instrument (AAL). The AAL combines gas jet levitation with actively controlled acoustic positioning to enable full optical access to samples ca. 2–4 mm in diameter, which are laser beam heated and melted. The cross sections of molten drops are imaged using a monochromatic light source and narrowband-filtered camera. Melt volume is calculated from the cross sections and used to find density at several temperatures ranging 1530–1920 K, including up to 350 K of supercooling. We report densities for CaAl 2 O 4 , Ca 12 Al 14 O 33 , CaSiO 3 , their Fe 2 O 3 -bearing analogs, and 83TiO 2 -17RE 2 O 3 (RE = La or Nd). These provide important benchmarks of the capabilities, measurement uncertainties, and future outlook for this technique.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Code Validation of SAM Using Forced and Natural Circulation Data from NACIE-UP Benchmark

Heavy liquid metals (HLMs) are promising candidates as coolants of Generation IV fast reactors due to their thermophysical properties. In the last decade, experimental work has been proposed as part of research and development efforts to develop such systems. In this context, researchers from the Brasimone Research Center have conducted many experiments using the Natural Circulation Experiment Upgrade (NACIE-UP) facility to study the thermofluid dynamic behavior of HLMs in rod bundle configurations with or without wire wrappers. This facility consists of a rectangular loop operated with lead-bismuth eutectic. Sensors across the loop monitor relevant parameters, i.e. temperatures, heat transfer, and flow conditions. Here, in the present work, we carefully select published data from NACIE-UP to validate the System Analysis Module (SAM), a modern system analysis code developed at Argonne National Laboratory. We developed one SAM model using specifications of the facility geometry and materials existing in relevant papers and reports. On top of that, these references provided the boundary conditions for simulating natural circulation and forced convection experiments in either steady or transient conditions. The SAM model simulates five test cases with diverse operating conditions. Ultimately, the code is proven to predict temperatures and mass flow rates that closely match the experiments. The discrepancies between numerical predictions and diverse transients are limited to a few degrees Celsius, showcasing that SAM is well suited for analyzing nuclear systems relying on HLM coolants.

advanced reactors↗
DOE OSTI2025

Toward machine learning interatomic potentials for modeling uranium mononitride

Uranium mononitride (UN) is a promising accident-tolerant fuel because of its high fissile density and high thermal conductivity. In this study, we developed the first machine learning interatomic potentials for reliable atomic-scale modeling of UN at finite temperatures. We constructed a training set using density functional theory (DFT) calculations that was enriched through an active learning procedure, and two neural network potentials were generated. Both potentials successfully reproduce key thermophysical properties of interest, such as temperature-dependent lattice parameter, specific heat capacity, and bulk modulus. We also evaluated the energy of stoichiometric defect reactions and defect migration barriers and found close agreement with DFT predictions, demonstrating that our potentials can be used for modeling defects in UN. Additional tests provide evidence that our potentials are reliable for simulating diffusion, noble gas impurities, and radiation damage.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Impact of dynamic Jahn-Teller effect on magnetic excitations, lattice vibration, and thermal conductivity in U 𝑥 ⁢T⁢h 1−𝑥 ⁡O 2 system

Vibrational and magnetic properties of single-crystal uranium-thorium dioxide (U 𝑥 ⁢T⁢h 1−𝑥⁡ O 2 ) with a full range of 0 < 𝑥 < 1 are investigated. Thorium dioxide is a diamagnet whose thermal properties are governed by lattice vibration. The addition of paramagnetic uranium ions leads to the emergence of magnetic effects that alter the thermophysical properties noticeably even at room temperature. The interaction of phonons with magnetic moments of uranium 5⁢𝑓 electrons mediated by magnetoelastic coupling results in an anomalous low-temperature thermal conductivity profile. Analysis of the magnetic susceptibility measurements indicates a departure from the Curie-Weiss relationship characteristic of noninteracting paramagnetic ions, previously associated with the dynamic Jahn-Teller (DJT) effect characterized by coupling between spin and the oxygen sublattice. The T 2⁢g Raman peak position follows a nonlinear trend as a function of uranium concentration and hints that these Raman active optical modes play a role in either DJT or in mediating quadrupole-quadrupole interactions. A first-principle-based thermal transport model is implemented to explain the low-temperature transport measurements, where the anomalous reduction is attributed to phonon-spin resonant scattering. The interplay between spins and phonons is also captured using high-resolution inelastic x-ray scattering (IXS) measurements of phonon linewidths. Our results provide insights into the phonon interactions with the magnetic excitations governing DJT effect and impacting the low-temperature thermal transport processes in this material system. Furthermore, these findings have implications for understanding low-temperature thermal transport and magnetic properties in advanced materials for information processing and energy applications.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Equation of state of boron carbide B 4 ⁢C

We present the results of recent experiments conducted on the Sandia Z machine and a new tabular equation of state for B 4 ⁢C. The equation of state was calibrated to a combination of density functional calculations reported here and fits to preexisting data. It was constructed partly to recover the effects of a shock-driven, polymorphic phase transition of unknown character beginning at particle velocities of just under 3 km/s (shock pressures of 95 GPa). Some of the Z experiments included sound speeds determined by the overtaking rarefaction method, from which we calculate the Grüneisen parameter and compare with previous experiments conducted at the OMEGA laser [Fratanduono et al ., Phys. Rev. B 94 , 184107 (2016)], our own first principles calculations, and another recent tabular equation of state [Zhang et al ., Phys. Rev. E 102 , 053203 (2020)]. We also compare our results with previous static compression, thermophysical, and melt studies, finding mixed consistency. We predict the onset and completion of shock melting at 225 and 265 GPa, respectively, and predict a melt curve that is largely flat to pressures of several hundred GPa.

36 MATERIALS SCIENCE↗