Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamics”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Insights into controlling bacterial cellulose nanofiber film properties through balancing thermodynamic interactions and colloidal dynamics

In recent years, nanocellulose has emerged as a sustainable and environmentally friendly alternative to traditional petroleum-derived structural polymers. Sourced either from plants, algae, or bacteria, nanocellulose can be processed into colloid, gel, film and fiber forms. However, the required fundamental understanding of process parameters that govern the morphology and structure–property relationships of nanocellulose systems, from colloidal suspensions to bulk materials, has not been developed and generalized for all forms of cellulose. This further hinders the more widespread adoption of this biopolymer in applications. Our study investigates the dispersion of cellulose nanofibers (CNFs) produced by a bacterial–yeast co-culture, in solvents, highlighting the role of thermodynamic interactions in influencing their colloidal behavior. By adjusting Hansen solubility parameters, we controlled the thermodynamic relationship between CNFs and solvents across various concentrations, studying the dilute to semi-dilute regimes. Rheological measurements revealed that the threshold at which a concentration-based regime transition occurs is distinctly solvent-dependent. Complementing rheological analysis with small angle X-ray scattering and zeta potential measurements, our findings reveal that enhancing CNF–solvent interactions increases excluded volume in the dilute regime, emphasizing the importance of the balance between fiber–fiber and fiber–solvent interactions. Moreover, we investigated the transition from colloidal to solid state by creating films from dispersions with varying interaction parameters in semi-dilute regimes. Through mechanical testing and scanning electron microscopy imaging of the fracture surfaces, we highlight the significance of electrokinetic effects in such transitions, as dispersions with higher electrokinetic stabilization gave rise to stronger and tougher films despite having less favorable thermodynamic interaction parameters. Finally, our work provides insights into the thermodynamic and electrokinetic interplay that governs bacterial CNF dispersion, offering a foundation for future application and a deeper understanding of nanocellulose's colloidal and structure-property relationships.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Computing chemical potentials with machine-learning-accelerated simulations to accurately predict thermodynamic properties of molten salts

The successful design and deployment of next-generation nuclear technologies heavily rely on thermodynamic data for relevant molten salt systems. However, the lack of accurate force fields and efficient methods has limited the quality of thermodynamic predictions from atomistic simulations. Here we propose an efficient free energy framework for computing chemical potentials, which is the central free energy quantity behind many thermodynamic properties. We accelerate our simulations without sacrificing accuracy by using machine learning interatomic potentials trained on density functional theory (DFT) data. Using lithium chloride as our model system, we compute chemical potentials with DFT-accuracy for solid and liquid phases by transmuting ions into noninteracting particles. Notably, in the liquid phase, we demonstrate consistency whether we transmute one ion pair or the entire system into ideal gas particles. By locating the temperature where the chemical potential of solid and liquid phases cross, we predict a melting point of 880 ± 18 K for lithium chloride, which is remarkably close to the experimental value of 883 K. With this successful demonstration, we lay the foundation for high-throughput thermodynamic predictions of many properties that can be derived from the chemical potentials of the minority and majority components in molten salts.

Gibson, Luke D. [Oak Ridge National Laboratory (OR↗
DOE OSTI2024

Learning thermodynamically constrained equations of state with uncertainty

Numerical simulations of high energy-density experiments require equation of state (EOS) models that relate a material’s thermodynamic state variables—specifically pressure, volume/density, energy, and temperature. EOS models are typically constructed using a semi-empirical parametric methodology, which assumes a physics-informed functional form with many tunable parameters calibrated using experimental/simulation data. Since there are inherent uncertainties in the calibration data (parametric uncertainty) and the assumed functional EOS form (model uncertainty), it is essential to perform uncertainty quantification (UQ) to improve confidence in EOS predictions. Model uncertainty is challenging for UQ studies since it requires exploring the space of all possible physically consistent functional forms. Thus, it is often neglected in favor of parametric uncertainty, which is easier to quantify without violating thermodynamic laws. This work presents a data-driven machine learning approach to constructing EOS models that naturally captures model uncertainty while satisfying the necessary thermodynamic consistency and stability constraints. We propose a novel framework based on physics-informed Gaussian process regression (GPR) that automatically captures total uncertainty in the EOS and can be jointly trained on both simulation and experimental data sources. A GPR model for the shock Hugoniot is derived, and its uncertainties are quantified using the proposed framework. We apply the proposed model to learn the EOS for the diamond solid state of carbon using both density functional theory data and experimental shock Hugoniot data to train the model and show that the prediction uncertainty is reduced by considering thermodynamic constraints.

Sharma, Himanshu (ORCID:000900050235718X)↗
DOE OSTI2024

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
DOE OSTI2025

Thermodynamics of calcium binding to heparin: Implications of solvation and water structuring for polysaccharide biofunctions

Heparan sulfates are found in all animal tissues and have essential roles in living systems. This family of biomacromolecules modulates binding to calcium ions (Ca 2+ ) in low free energy reactions that influence biochemical processes from cell signaling and anticoagulant efficacy to biomineralization. Despite their ubiquity, the thermodynamic basis for how heparans and similarly functionalized biomolecules regulate Ca 2+ interactions is not yet established. Using heparosan (Control) and heparins with different positions of sulfate groups, we quantify how SO 3 − and COO − content and SO 3 − position modulate Ca 2+ binding by isothermal titration calorimetry. The free energy of all heparin-Ca2+ interactions (ΔG rxn ) is dominated by entropic contributions due to favorable water release from polar, hydrophilic groups. Heparin with both sulfate esters (O-SO 3 − ) and sulfamides (N-SO 3 − ) has the strongest binding to Ca 2+ compared to heparosan and to heparin with only O-SO 3 − groups (~3X). By linking Ca 2+ binding thermodynamics to measurements of the interfacial energy for calcite (CaCO 3 ) crystallization onto polysaccharides, we show molecule-specific differences in nucleation rate can be explained by differences in water structuring during Ca 2+ interactions. A large entropic term (-TΔS rxn ) upon Ca 2+ –polysaccharide binding correlates with high interfacial energy to CaCO 3 nucleation. Combining our measurements with literature values indicates many Ca 2+ –polysaccharide interactions have a shared thermodynamic signature. The resulting enthalpy–entropy compensation relationship suggests these interactions are generally dominated by water restructuring involving few conformational changes, distinct from Ca 2+ –protein binding. Our findings quantify the thermodynamic origins of heparin-specific interactions with Ca 2+ and demonstrate the contributions of solvation and functional group position during biomacromolecule-mediated ion regulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Ab initio anharmonic thermodynamic properties of cubic Ca Si O 3 perovskite

We present an ab initio study of the thermodynamic properties of cubic CaSiO 3 perovskite (CaPv) over the pressure and temperature range of the Earth's lower mantle. We compute the anharmonic phonon dispersions throughout the Brillouin zone by utilizing the phonon quasiparticle approach, which characterizes the intrinsic temperature dependence of phonon frequencies and, in principle, captures full anharmonicity. Such temperature-dependent phonon dispersions are used to calculate ab initio free energy in the thermodynamic limit (N → ∞) within the framework of the phonon gas model. Accurate free energy calculations enable us to investigate cubic CaPv's thermodynamic properties, e.g., thermal expansivity, Grüneisen parameter, bulk modulus, heat capacity and thermal equation of state, where anharmonic effects are demonstrated. Furthermore, the present methodology provides an important theoretical approach for exploring phase boundaries, thermodynamic, and thermoelastic properties of strongly anharmonic materials at high pressures and temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Increasing the Thermodynamic Driving Force of the Phosphofructokinase Reaction in Clostridium thermocellum

Glycolysis is an ancient, widespread, and highly conserved metabolic pathway that converts glucose into pyruvate. In the canonical pathway, the phosphofructokinase (PFK) reaction plays an important role in controlling flux through the pathway. Clostridium thermocellum has an atypical glycolysis and uses pyrophosphate (PP i ) instead of ATP as the phosphate donor for the PFK reaction. The reduced thermodynamic driving force of the PP i -PFK reaction shifts the entire pathway closer to thermodynamic equilibrium, which has been predicted to limit product titers. Here, in this study, we replace the PP i -PFK reaction with an ATP-PFK reaction. We demonstrate that the local changes are consistent with thermodynamic predictions: the ratio of fructose 1,6-bisphosphate to fructose-6-phosphate increases, and the reverse flux through the reaction (determined by 13 C labeling) decreases. The final titer and distribution of fermentation products, however, do not change, demonstrating that the thermodynamic constraints of the PP i -PFK reaction are not the sole factor limiting product titer.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Thermosteam: BioSTEAM's Premier Thermodynamic Engine

Thermosteam is a thermodynamic engine capable of solving mass and energy balances, estimating mixture properties, solving thermodynamic phase equilibria, and modeling stoichiometric reactions. All chemical data in Thermosteam is imported from the chemicals library, an open-source compilation of data and functions for the estimation of pure component chemical and mixture properties. Thermosteam’s fast and flexible platform has enabled the evaluation of conceptual and emerging biochemical production processes. The Biorefinery Simulation and Techno-Economic Analysis Modules (BioSTEAM) — capable of modeling reactors, distillation columns, heat exchangers, and other unit operations — has adopted Thermosteam as its premier thermodynamic engine. Published biorefinery designs modeled in BioSTEAM implement thermodynamic property packages created with Thermosteam, including a cornstover biorefinery for the production of cellulosic ethanol, a lipid-cane biorefinery for the co-production of ethanol and biodiesel, and a wheatstraw biorefinery for the production of cellulosic ethanol.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

M3SF-24LL010302062-NEA-TDB Management and International Collaborations in Sorption and Thermodynamic Modeling

This progress report (Level 3 Milestone Number M3SF-24LL010302062) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Crystalline International Collaborations Work Package Number SF-24LL01030206. The activity is focused on our long-term commitment to engaging our partners in international nuclear waste repository research. This includes participation in the Nuclear Energy Agency Thermochemical Database (NEA-TDB) Project and development of methodologies for integrating US and international thermodynamic databases for use in SFWST Generic Disposal System Assessment (GDSA) efforts. A continuing focus for FY24 efforts is to support the US participation in the NEA-TDB effort. The focus of FY24 activities was the development of an agreement for a Phase 7 activity that will start in Q1 of 2025. Mavrik Zavarin is now the US representative on both the Management Board and the Executive group to the NEA-TDB. He is also the POC for the Cements State of the Art Report that is undergoing peer review in FY24. In FY24, we used our position on the NEA-TDB MB and EG to facilitate the integration of NEA-TDB thermochemical data with LLNL’s SUPCRTNE thermodynamic database that supports the SFWST GDSA activities. This effort is coordinated with the Argillite work package SUPCRTNE database development efforts (Wolery, 2024). The goal is to provide a downloadable database that will be hosted on LLNL’s thermodynamics website which incorporates NEA-TDB data into the LLNL database where appropriate. We also began engagement with the EURAD2 program that was initiated in FY24 by our European collaborators at the Karlsruhe Institute of Technology (KIT), Germany. The primary focus of the engagement is with WP20: DITUSC Thermodynamic database evaluation program. A kickoff meeting for this activity is planned for early FY25. Finally, we have been selected to co-host (with Clemson University) the International Conference on Chemistry and Migration Behaviour of Actinides and Fission Products in the Geosphere in 2025 (Migration2025). The meeting will be held September 21-26, 2025, in New Orleans, Louisiana, and will focus on international efforts to understand the risks of radionuclide releases into the environment. This central focus of this conference is on international efforts to develop safe disposal options for nuclear wastes. As such, we are developing a theme focused on US underground nuclear waste repository science.

58 GEOSCIENCES↗
DOE OSTI2020

As a Matter of State: The Role of Thermodynamics in Magnetohydrodynamic Turbulence

Turbulence simulations play a key role in advancing the general understanding of the physical properties of turbulence and in interpreting astrophysical observations of turbulent plasmas. For the sake of simplicity, however, turbulence simulations are often conducted in the isothermal limit. Given that the majority of astrophysical systems are not governed by isothermal dynamics, we aim to quantify the impact of thermodynamics on the physics of turbulence, through varying adiabatic index, γ, combined with a range of optically thin cooling functions. Here, we present a suite of ideal magnetohydrodynamics simulations of thermally balanced stationary turbulence in the subsonic, super-Alfvénic, high ${\beta }_{{\rm{p}}}$ (ratio of thermal to magnetic pressure) regime, where turbulent dissipation is balanced by two idealized cooling functions (approximating linear cooling and free–free emission) and examine the impact of the equation of state by considering cases that correspond to isothermal, monatomic, and diatomic gases. We find a strong anticorrelation between thermal and magnetic pressure independent of thermodynamics, whereas the strong anticorrelation between density and magnetic field found in the isothermal case weakens with increasing γ. Similarly, the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing γ. This suggests that there exists a degeneracy in these relations with respect to thermodynamics and Mach number in this regime, which is dominated by slow magnetosonic modes. These results have implications for attempts to infer (e.g.,) Mach numbers from (e.g.,) Faraday rotation measurements, without additional information regarding the thermodynamics of the plasma. However, our results suggest that this degeneracy can be broken by utilizing higher-order moments of observable distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Thermodynamic properties of oxygen and nitrogen III

The final equation for nitrogen was determined. In the work on the equation of state for nitrogen, coefficients were determined by constraining the critical point to selected critical point parameters. Comparisons of this equation with all the P-density-T data were made, as well as comparisons to all other thermodynamic data reported in the literature. The extrapolation of the equation of state was studied for vapor to higher temperatures and lower temperatures, and for the liquid surface to the saturated liquid and the fusion lines. A new vapor pressure equation was also determined which was constrained to the same critical temperature, pressure, and slope (dP/dT) as the equation of state. Work on the equation of state for oxygen included studies for improving the equation at the critical point. Comparisons of velocity of sound data for oxygen were also made between values calculated with a preliminary equation of state and experimental data. Functions for the calculation of the derived thermodynamic properties using the equation of state are given, together with the derivative and integral functions for the calculation of the thermodynamic properties using the equations of state. Summary tables of the thermodynamic properties of nitrogen and oxygen are also included to serve as a check for those preparing computer programs using the equations of state.

Stewart, R. B.↗
NASA NTRS1973

Thermodynamic and transport properties of gaseous tetrafluoromethane in chemical equilibrium

Equations and in computer code are presented for the thermodynamic and transport properties of gaseous, undissociated tetrafluoromethane (CF4) in chemical equilibrium. The computer code calculates the thermodynamic and transport properties of CF4 when given any two of five thermodynamic variables (entropy, temperature, volume, pressure, and enthalpy). Equilibrium thermodynamic and transport property data are tabulated and pressure-enthalpy diagrams are presented.

Hunt, J. L.↗
NASA NTRS1979

Free energy of formation of Mo2C and the thermodynamic properties of carbon in solid molybdenum

As part of a study of the thermodynamical properties of interstitial elements in refractory metals, the free energy of formation of Mo2C is determined, and the thermodynamical properties of C in solution in solid Mo evaluated. The activity of C in the two-phase region Mo + Mo2C is obtained from the C content of iron rods equilibrated with metal + carbide powder mixtures. The free energy of formation of alpha-Mo2C is determined from the activity data. The thermodynamic properties of C in the terminal solid solution are calculated from available data on the solid solubility of C in Mo. Lattice distortion due to misfit of the C atoms in the interstitial sites appears to play a significant role in determining the thermodynamic properties of C in solid Mo.

Seigle, L. L.↗
NASA NTRS1981

Thermodynamic equilibrium-air correlations for flowfield applications

Equilibrium-air thermodynamic correlations have been developed for flowfield calculation procedures. A comparison between the postshock results computed by the correlation equations and detailed chemistry calculations is very good. The thermodynamic correlations are incorporated in an approximate inviscid flowfield code with a convective heating capability for the purpose of defining the thermodynamic environment through the shock layer. Comparisons of heating rates computed by the approximate code and a viscous-shock-layer method are good. In addition to presenting the thermodynamic correlations, the impact of several viscosity models on the convective heat transfer is demonstrated.

Zoby, E. V.↗
NASA NTRS1988

The development of flux-split algorithms for flows with non-equilibrium thermodynamics and chemical reactions

A finite-volume method for the numerical computation of flows with nonequilibrium thermodynamics and chemistry is presented. A thermodynamic model is described which simplifies the coupling between the chemistry and thermodynamics and also results in the retention of the homogeneity property of the Euler equations (including all the species continuity and vibrational energy conservation equations). Flux-splitting procedures are developed for the fully coupled equations involving fluid dynamics, chemical production and thermodynamic relaxation processes. New forms of flux-vector split and flux-difference split algorithms are embodied in a fully coupled, implicit, large-block structure, including all the species conservation and energy production equations. Several numerical examples are presented, including high-temperature shock tube and nozzle flows. The methodology is compared to other existing techniques, including spectral and central-differenced procedures, and favorable comparisons are shown regarding accuracy, shock-capturing and convergence rates.

Grossman, B.↗
NASA NTRS1992

Evolution of a terrestrial magma ocean: Thermodynamics, kinetics, rheology, convection, differentiation

The evolution of an initially totally molten magma ocean is constrained on the basis of analysis of various physical problems in the magma ocean. First of all an equilibrium thermodynamics of the magma ocean is developed in the melting temperature range. The equilibrium thermodynamical parameters are found as functions only of temperature and pressure and are used in the subsequent models of kinetics and convection. Kinematic processes determine the crystal size and also determine a non-equilibrium thermodynamics of the system. Rheology controls all dynamical regimes of the magma ocean. The thermal convection models for different rheological laws are developed for both the laminar convection and for turbulent convection in the case of equilibrium thermodynamics of the multiphase system. The evolution is estimated on the basis of all the above analysis.

Solomatov, V. S.↗
NASA NTRS1992

Thermodynamics, Solubility, and Diffusivity of Oxygen in Titanium and Ti-Al Alloys

Titanium aluminides and titanium aluminide-based composites are attractive candidate materials for high-temperature structural applications. As these materials may be exposed to oxidizing environments durine their use at elevated temperatures, it is essential that they possess a good oxidation resistance. Previous studies have shown that the oxidation resistance of Al-rich alloys in the Ti-Al system is superior to that of the Ti-rich alloys. The scales formed on the surface of the Al-rich and Ti-rich alloys have been reported to be predominantly Al2O3 and TiO2, respectively. Since the relative stabilities of the oxides of Al and Ti at various temperatures and oxygen pressures can be assessed from their thermodynamic data, it is possible, With the help of thermodynamic calculations, to determine the compositions of the alloys which would form scales of Al2O3, TiO(x) or a ternary oxide such as TiAl2O5 during oxidation at a given temperature. The thermodynamic calculations require reliable activity data for the Ti-Al system. These data have not been determined for the entire composition and temperature range of interest. Using the data available in the literature, recently performed thermodynamic calculations and concluded that the stable oxide changed from TiO to Al2O3 in the existence region of the tial phase. In the case of titanium aluminide-based composites, another major concern is the mutual chemical compatibility of the matrix material with the reinforcement phase. Fibers of SiC, TiB2 and Al2O3 are currently being investigated for reinforcement of titanium aluminide matrices.

Mehrotra, Gopal M.↗
NASA NTRS2004

An Interpolation Method for Obtaining Thermodynamic Properties Near Saturated Liquid and Saturated Vapor Lines

The two most common approaches used to formulate thermodynamic properties of pure substances are fundamental (or characteristic) equations of state (Helmholtz and Gibbs functions) and a piecemeal approach that is described in Adebiyi and Russell (1992). This paper neither presents a different method to formulate thermodynamic properties of pure substances nor validates the aforementioned approaches. Rather its purpose is to present a method to generate property tables from existing property packages and a method to facilitate the accurate interpretation of fluid thermodynamic property data from those tables. There are two parts to this paper. The first part of the paper shows how efficient and usable property tables were generated, with the minimum number of data points, using an aerospace industry standard property package. The second part describes an innovative interpolation technique that has been developed to properly obtain thermodynamic properties near the saturated liquid and saturated vapor lines.

Nguyen, Huy H.↗