Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Projecting Permafrost Thaw of Sub-Arctic Tundra With a Thermodynamic Model Calibrated to Site Measurements

Northern circumpolar permafrost thaw affects global carbon cycling, as large amounts of stored soil carbon becomes accessible to microbial breakdown under a warming climate. The magnitude of carbon release is linked to the extent of permafrost thaw, which is locally variable and controlled by soil thermodynamics. Soil thermodynamic properties, such as thermal diffusivity, govern the reactivity of the soil-atmosphere thermal gradient, and are controlled by soil composition and drainage. In order to project permafrost thaw for an Alaskan tundra experimental site, we used seven years of site data to calibrate a soil thermodynamic model using a data assimilation technique. The model reproduced seasonal and interannual temperature dynamics for shallow (5–40 cm) and deep soil layers (2–4 m), and simulations of seasonal thaw depth closely matched observed data. The model was then used to project permafrost thaw at the site to the year 2100 using climate forcing data for three future climate scenarios (RCP 4.5, 6.0, and 8.5). Minimal permafrost thawing occurred until mean annual air temperatures rose above the freezing point, after which we measured over a 1 m increase in thaw depth for every 1 °C rise in mean annual air temperature. Under no projected warming scenario was permafrost remaining in the upper 3 m of soil by 2100. We demonstrated an effective data assimilation method that optimizes parameterization of a soil thermodynamic model. The sensitivity of local permafrost to climate warming illustrates the vulnerability of sub-Arctic tundra ecosystems to significant and rapid soil thawing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Impact of a New Sea Ice Thermodynamic Formulation in the CESM2 Sea Ice Component

Abstract The sea ice component of the Community Earth System Model version 2 (CESM2) contains new “mushy‐layer” physics that simulates prognostic salinity in the sea ice, with consequent modifications to sea ice thermodynamics and the treatment of melt ponds. The changes to the sea ice model and their influence on coupled model simulations are described here. Two simulations were performed to assess the changes in the vertical thermodynamics formulation with prognostic salinity compared to a constant salinity profile. Inclusion of the mushy layer thermodynamics of Turner et al. (2013, https://doi.org/10.1002/jgrc.20171 ) in a fully coupled Earth system model produces thicker and more extensive sea ice in the Arctic, with relatively unchanged sea ice in the Antarctic compared to simulations using a constant salinity profile. While this is consistent with the findings of uncoupled ice‐ocean model studies, the role of the frazil and congelation growth is more important in fully coupled simulations. Melt pond drainage is also an important contribution to simulated ice thickness differences as also found in the uncoupled simulations of Turner and Hunke (2015; https://doi.org/10.1002/2014JC010358 ). However, it is an interaction of the ponds and the snow fraction that impacts the surface albedo and hence the top melt. The changes in the thermodynamics and resulting ice state modify the ice‐ocean‐atmosphere fluxes with impacts on the atmosphere and ocean states, particularly temperature.

Bailey, David A.↗
DOE OSTI2020

Thermodynamics, dynamics, and structure of supercritical water at extreme conditions

Molecular dynamics (MD) simulations to understand the thermodynamic, dynamic, and structural changes in supercritical water across the Frenkel line and the melting line have been performed. The two-phase thermodynamic model [J. Phys. Chem. B(2010), 114(24), 8191-8198] and the velocity autocorrelation functions are used to locate the Frenkel line and to calculate the thermodynamic and dynamic properties. The Frenkel lines obtained from the two-phase thermodynamic model and the velocity autocorrelation criterion do not agree with each other. Structural characteristics and the translational diffusion dynamics of water suggest that this inconsistency could arise from the two oscillatory modes in water, which are associated with the bending of hydrogen bonds and intermolecular collisions inside the first coordination shell. Finally, the overall results lead us to conclude that the universality of the Frenkel line as a dynamic crossover line from rigid to nonrigid fluids is preserved in water.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Thermodynamic guiding principles of high-capacity phase transformation materials for splitting H 2 O and CO 2 by thermochemical looping

Here, thermochemical looping splitting of water and carbon dioxide (CO 2 ) with greenhouse-gas-free (GHG-free) energy has the potential to help address the Gt-scale GHG emissions challenge. Reaction thermodynamics largely contributes to the main bottlenecks of cost reduction for thermochemical looping water/CO 2 splitting cycle. Here, we analyze thermodynamic driving forces in such cycles with two-phase ternary ferrites as model systems. We find that cation configurational entropy chiefly determines the change of partial molar entropy with oxygen stoichiometry. In addition, our phase diagram analysis accurately predicts the optimal Fe ratio for maximal water/CO 2 splitting capacity in thermal reduction and in chemical reduction based cycles, underlining the significance of phase boundary positions. With chemical reduction, >10% CO 2 conversion and high oxygen exchange capacity can both be achieved. Furthermore, our reduced Gibbs free energy model illustrates critical thermodynamic factors that influence the water/CO 2 splitting capacity. Our research reveals the thermodynamic driving forces underlying the unconventional high-capacity Fe-poor ferrites, further explained via phase diagrams of Fe–Co–O, Fe–Ni–O and Fe–Mg–O. Future materials improvements can be guided by our reduced Gibbs free energy model.

08 HYDROGEN↗
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 OSTI2014

An analytical coarse-graining method which preserves the free energy, structural correlations, and thermodynamic state of polymer melts from the atomistic to the mesoscale

Structural and thermodynamic consistency of coarse-graining models across multiple length scales is essential for the predictive role of multi-scale modeling and molecular dynamic simulations that use mesoscale descriptions. Our approach is a coarse-grained model based on integral equation theory, which can represent polymer chains at variable levels of chemical details. The model is analytical and depends on molecular and thermodynamic parameters of the system under study, as well as on the direct correlation function in the k → 0 limit, c{sub 0}. A numerical solution to the PRISM integral equations is used to determine c{sub 0}, by adjusting the value of the effective hard sphere diameter, d{sub HS}, to agree with the predicted equation of state. This single quantity parameterizes the coarse-grained potential, which is used to perform mesoscale simulations that are directly compared with atomistic-level simulations of the same system. We test our coarse-graining formalism by comparing structural correlations, isothermal compressibility, equation of state, Helmholtz and Gibbs free energies, and potential energy and entropy using both united atom and coarse-grained descriptions. We find quantitative agreement between the analytical formalism for the thermodynamic properties, and the results of Molecular Dynamics simulations, independent of the chosen level of representation. In the mesoscale description, the potential energy of the soft-particle interaction becomes a free energy in the coarse-grained coordinates which preserves the excess free energy from an ideal gas across all levels of description. The structural consistency between the united-atom and mesoscale descriptions means the relative entropy between descriptions has been minimized without any variational optimization parameters. The approach is general and applicable to any polymeric system in different thermodynamic conditions.

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

The role of material flexibility on the drying transition of water between hydrophobic objects: A thermodynamic analysis

Liquid water confined between hydrophobic objects of sufficient size becomes metastable with respect to its vapor at separations smaller than a critical drying distance. Macroscopic thermodynamic arguments predicting this distance have been restricted to the limit of perfectly rigid confining materials. However, no material is perfectly rigid and it is of interest to account for this fact in the thermodynamic analysis. We present a theory that combines the current macroscopic theory with the thermodynamics of elasticity to derive an expression for the critical drying distance for liquids confined between flexible materials. The resulting expression is the sum of the well-known drying distance for perfectly rigid confining materials and a new term that accounts for flexibility. Thermodynamic arguments show that this new term is necessarily positive, meaning that flexibility increases the critical drying distance. To study the expected magnitude and scaling behavior of the flexible term, we consider the specific case of water and present an example of drying between thin square elastic plates that are simply supported along two opposite edges and free at the remaining two. We find that the flexible term can be the same order of magnitude or greater than the rigid solution for materials of biological interest at ambient conditions. In addition, we find that when the rigid solution scales with the characteristic size of the immersed objects, the flexible term is independent of size and vice versa. Thus, the scaling behavior of the overall drying distance will depend on the relative weights of the rigid and flexible contributions.

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

Thermodynamic scaling of molecular dynamics in supercooled liquid state of pharmaceuticals: Itraconazole and ketoconazole

Pressure-Volume-Temperature (PVT) measurements and broadband dielectric spectroscopy were carried out to investigate molecular dynamics and to test the validity of thermodynamic scaling of two homologous compounds of pharmaceutical activity: itraconazole and ketoconazole in the wide range of thermodynamic conditions. The pressure coefficients of the glass transition temperature (dT{sub g}/dp) for itraconazole and ketoconazole were determined to be equal to 183 and 228 K/GPa, respectively. However, for itraconazole, the additional transition to the nematic phase was observed and characterized by the pressure coefficient dT{sub n}/dp = 258 K/GPa. From PVT and dielectric data, we obtained that the liquid-nematic phase transition is governed by the relaxation time since it occurred at constant τ {sub α} = 10{sup −5} s. Furthermore, we plotted the obtained relaxation times as a function of T{sup −1}v{sup −γ}, which has revealed that the validity of thermodynamic scaling with the γ exponent equals to 3.69 ± 0.04 and 3.64 ± 0.03 for itraconazole and ketoconazole, respectively. Further analysis of the scaling parameter in itraconazole revealed that it unexpectedly decreases with increasing relaxation time, which resulted in dramatic change of the shape of the thermodynamic scaling master curve. While in the case of ketoconazole, it remained the same within entire range of data (within experimental uncertainty). We suppose that in case of itraconazole, this peculiar behavior is related to the liquid crystals’ properties of itraconazole molecule.

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

Thermodynamic properties and transport coefficients of air thermal plasmas mixed with ablated vapors of Cu and polytetrafluoroethylene

Because the fault arc in aircraft electrical system often causes a fire, it is particularly important to analyze its energy and transfer for aircraft safety. The calculation of arc energy requires the basic parameters of the arc. This paper is mainly devoted to the calculations of equilibrium composition, thermodynamic properties (density, molar weight, enthalpy, and specific heat at constant pressure) and transport coefficients (thermal conductivity, electrical conductivity, and viscosity) of plasmas produced by a mixture of air, Cu, and polytetrafluoroethylene under the condition of local thermodynamic equilibrium. The equilibrium composition is determined by solving a system of equations around the number densities of each species. The thermodynamic properties are obtained according to the standard thermodynamic relationships. The transport coefficients are calculated using the Chapman-Enskog approximations. Results are presented in the temperature range from 3000 to 30 000 K for pressures of 0.08 and 0.1 MPa, respectively. The results are more accurate and are reliable reference data for theoretical analysis and computational simulation of the behavior of fault arc.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2016

Quantum work and the thermodynamic cost of quantum measurements

Quantum work is usually determined from two projective measurements of the energy at the beginning and at the end of a thermodynamic process. However, this paradigm cannot be considered thermodynamically consistent as it does not account for the thermodynamic cost of these measurements. To remedy this conceptual inconsistency we introduce a paradigm that relies only on the expected change of the average energy given the initial energy eigenbasis. In particular, we completely omit quantum measurements in the definition of quantum work, and hence quantum work is identified as a thermodynamic quantity of only the system. As main results we derive a modified quantum Jarzynski equality and a sharpened maximum work theorem in terms of the information free energy. Lastly, a comparison of our results with the standard approach allows one to quantify the informational cost of projective measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2018

Thermodynamic Properties of Strontium-Lead Alloys Determined by Electromotive Force Measurements

The thermodynamic properties of Sr-Pb alloys were determined by electromotive force (emf) measurements. A Sr(s)|CaF 2 -SrF 2 |Sr(in Pb) electrochemical cell was used to measure emf values at 773–1073 K for Sr-Pb alloys at mole fractions x Sr = 0.07–0.59. These emf measurements were used to determine thermodynamic properties of Sr-Pb alloys, including activity, partial molar entropy, and partial molar enthalpy. At 873 K, activity values of Sr in Pb were as low as a Sr = 1.72 × 10 –9 at mole fraction x Sr = 0.07, implying strong atomic interactions between Sr and Pb. Phase transition temperatures of Sr-Pb alloys, observed during emf measurements, were corroborated by thermal analysis (0.07 ≤ x Sr ≤ 0.34) by differential scanning calorimetry (DSC), and the phase constituents of Sr-Pb alloys (0.07 ≤ x Sr ≤ 0.75) were characterized using X-ray diffraction (XRD). Experimentally-determined thermodynamic properties were compared to the assessed thermodynamic properties of the Sr-Pb system, confirming the phase transition temperatures and highlighting discrepancies in solution properties (activity and excess Gibbs energy).

36 MATERIALS SCIENCE↗
DOE OSTI2019

Thermodynamic versus Dynamic Controls on Extreme Precipitation in a Warming Climate from the Community Earth System Model Large Ensemble

The moisture budget is evaluated as a function of the probability distribution of precipitation for the end of the twentieth century and projected end of the twenty-first century in the Community Earth System Model Large Ensemble. For a given precipitation percentile, a conditional moisture budget equation relates precipitation minus evaporation ( P − E) to vertical moisture transport, horizontal moisture advection, and moisture storage. At high percentiles, moisture advection and moisture storage cancel and evaporation is negligible, so that precipitation is approximately equal to vertical moisture transport, and likewise for projected changes. Therefore, projected changes to extreme precipitation are approximately equal to the sum of thermodynamic and dynamic tendencies, representing changes to the vertical profiles of moisture content and mass convergence, respectively. The thermodynamic tendency is uniform across percentiles and regions as an intensification of the hydrological cycle, but the dynamic tendency is more complex. For extreme events, per degree of warming, in the mid-to-high latitudes the dynamic tendency is small, so that precipitation approximately scales by the Clausius–Clapeyron 7% K −1 increase. In the subtropics, a drying tendency originating from dynamics offsets the thermodynamic wetting tendency, with the net effect on precipitation varying among regions. The effect of this dynamic drying decreases with increasing percentile. In the deep tropics, a positive dynamic tendency occurs with magnitude similar to or greater than the positive thermodynamic tendency, resulting in generally a 10%–15% K −1 precipitation increase, and with a >25% K −1 increase over the tropical east Pacific. This reinforcing dynamical tendency increases rapidly for high percentiles.

54 ENVIRONMENTAL SCIENCES↗