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

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗
DOE OSTI2025

Detecting thermodynamic phase transition via explainable machine learning of photoemission spectroscopy

Identifying thermodynamic signatures of electronic phases, such as superconductivity, is challenging in low-dimensional materials due to strong fluctuations and low probing volume. Spectroscopic methods are often used to identify new bulk phases, but their main measurable quantity—electronic energy gaps—is no longer an effective order parameter in low-dimensional and fluctuating systems. Combining angle-resolved photoemission with a domain-adversarial neural network, we report a data-driven method to identify thermodynamic phase transitions solely based on single-particle spectra. We demonstrate 97.6% accuracy in cuprate superconductor Bi 2 Sr 2 CaCu 2 O 8+δ with strong superconducting fluctuations. This model notably compensates for the scarcity of experimental data by leveraging virtually inexhaustible simulated data. Further, its explainability reveals the crucial role of in-gap spectral weight in detecting phase fluctuations and thermodynamic transitions. Our work pinpoints the spectroscopic signatures of fluctuating orders and enables using spectroscopy for machine-learning-assisted material discovery for low-dimensional and strong coupling systems.

2D materials↗
DOE OSTI2024

Adsorption Thermodynamics for Process Simulation

Adsorption has rapidly evolved in recent decades and is an established separation technology extensively practiced in gas separation industries and others. However, rigorous thermodynamic modeling of multicomponent adsorption equilibrium remains elusive, and industrial practitioners rely heavily on expensive and time-consuming trial-and-error pilot studies to develop adsorption units. Here, this article highlights the need for rigorous adsorption thermodynamic models and the limitations and deficiencies of existing models such as the extended Langmuir isotherm, dual-process Langmuir isotherm, and adsorbed solution theory. It further presents a series of recent advances in the generalization of the classical Langmuir isotherm of single-component adsorption by deriving an activity coefficient model to account for the adsorbed phase adsorbate–adsorbent interactions, substituting adsorbed phase adsorbate and vacant site concentrations with activities, and extending to multicomponent competitive adsorption equilibrium, both monolayer and multilayer. Requiring a minimum set of physically meaningful model parameters, the generalized Langmuir isotherm for monolayer adsorption and the generalized Brunauer–Emmett–Teller isotherm for multilayer adsorption address various thermodynamic modeling challenges including adsorbent surface heterogeneity, isosteric enthalpies of adsorption, BET surface areas, adsorbed phase nonideality, adsorption azeotrope formation, and multilayer adsorption. Also discussed is the importance of quality adsorption data that cover sufficient temperature, pressure, and composition ranges for reliable determination of the model parameters to support adsorption process simulation, design, and optimization.

09 BIOMASS FUELS↗
DOE OSTI2024

Role of Metal–Organic Framework Topology on Thermodynamics of Polyoxometalate Encapsulation

Polyoxometalates (POMs) are discrete anionic clusters whose rich redox properties, strong Brønsted acidity, and high availability of active sites poise them as potent catalysts for oxidation reactions. Here, metal–organic frameworks (MOFs) have emerged as tunable, porous platforms to immobilize POMs, thus increasing their solution stability and catalytic activity. While POM@MOF composite materials have been widely used for a variety of applications, little is known about the thermodynamics of the encapsulation process. Here, we utilize an up-and-coming technique in the field of heterogeneous materials, isothermal titration calorimetry (ITC), to obtain full thermodynamic profiles (ΔH, ΔS, ΔG, K a ) of POM binding. Six different 8-connected hexanuclear Zr-MOFs were investigated to determine the impact of MOF topology (csq, scu, the) on POM encapsulation thermodynamics.

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

First-Principles Insights into the Thermodynamics of Variable-Temperature Ammonia Synthesis on Transition-Metal-Doped Cu (100) and (111)

Ammonia (NH 3 ) is one of the most produced chemicals worldwide. NH 3 synthesis predominantly utilizes the Haber–Bosch (HB) process, requiring high temperatures and pressures. Despite significant process advances, ample opportunity remains for improving the rate, selectivity, catalyst stability, and energy efficiency. Inspired by a recently developed programmable heating and quenching (PHQ) technique, we present in this paper a first-principles screening of candidate single-atom alloy catalysts generated from doping (111) and (100) surfaces of copper (Cu), an ineffective HB catalyst in its pure form. We predict the thermodynamics of two rate-limiting reactions, N 2 dissociative adsorption and the final hydrogenation step leading up to NH 3 release, at 400 and 900 K. Thermodynamically, the former reaction is favored at low temperatures, while the latter is favored at high temperatures. Vanadium-, chromium-, and molybdenum-doped Cu surfaces, due to intermediate M–N covalent bonding character, emerge as appealing candidate catalysts for PHQ NH 3 synthesis, as they balance the thermodynamics of the above-mentioned reaction steps at their respective optimal temperatures.

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

First-principles thermodynamics of Al 10 ⁢V: An analytical treatment of localized anharmonic modes

Many complex intermetallic structures possess cagelike environments that can host additional guest atoms. In Al 10 ⁢V, these atoms give rise to low-frequency, localized vibrations (Einstein modes) that dominate the thermodynamic response at low temperature. They become imaginary under volume expansion as temperature rises, invalidating the harmonic approximation. Here, we develop a framework to incorporate these strongly anharmonic vibrational modes into first-principles thermodynamic calculations. By explicitly modeling the cage potential and solving the associated Schrödinger equation numerically, we compute the full anharmonic free energy contribution and demonstrate its impact on the thermodynamic behavior of Al 10⁢ V. This allows us to examine structures with different cage fillings and construct the Al-V phase diagram in the relevant composition range. Our results reproduce key experimental signatures, including the anomalous rise in the thermal expansion coefficient and heat capacity at low temperatures, and reveal that the presence and the extent of cage filling by guest atoms is essential to stabilizing the Al 10 ⁢V phase at elevated temperatures.

anharmonic lattice dynamics↗
DOE OSTI2025

Metriplectic four-bracket algorithm for constructing thermodynamically consistent dynamical systems

A unified thermodynamic algorithm is presented for constructing thermodynamically consistent dynamical systems, i.e., systems that have Hamiltonian and dissipative parts that conserve energy while producing entropy. The algorithm is based on the metriplectic 4-bracket given in Morrison and Updike [Phys. Rev. E 109, 045202 (2024)]. A feature of the unified thermodynamic algorithm is the force-flux relation 𝐉 𝛼 =−𝐿 𝛼⁢𝛽 𝛁(𝛿⁢𝐻⁡/𝛿⁢𝜉 𝛽 ) for phenomenological coefficients 𝐿 𝛼⁢𝛽 , Hamiltonian 𝐻, and dynamical variables 𝜉 𝛽 . The algorithm is applied to the Navier-Stokes-Fourier, the Cahn-Hilliard-Navier-Stokes, and Brenner-Navier-Stokes-Fourier systems, and significant generalizations of these systems are obtained.

Multiphase flows↗
DOE OSTI2026

Thermodynamics of grain-boundary phases

The past decade has seen a significant increase in research efforts aimed at understanding the thermodynamics of low-dimensional phases existing in many materials systems, ranging from two-dimensional materials to core regions of extended defects in crystalline solids. We review the current status of theoretical, computational, and experimental research on the “defect phases,” focusing on grain boundaries (GBs) in elemental and multicomponent polycrystalline materials. After reviewing the generalized concept of a phase of any dimensionality, we discuss recent progress in atomistic computer simulations of GB phase transformations and phase coexistences, including the observation of one-dimensional defects separating GB phases (defects in defects). Computational predictions compare well with experimental observations of multiple GB phases and segregation-induced phase transformations. An intriguing open question of GB thermodynamics is whether the GB free energy can be driven to a zero value by increasing solute segregation. We review recent efforts to understand this ultimate thermodynamic stabilization of GB phases and the possible polycrystalline microstructures that may arise. An outlook for future research in the field is discussed.

Materials science↗
DOE OSTI2025

First-Principles Thermodynamic Assessments of Sr-Containing Secondary Phase Formation in La1-xSrxMnO3±δ Perovskites for Solid Oxide Cell Applications

Sr-secondary phase formation is a potentially significant degradation mode threatening solid-oxide cell (SOC) commercial viability. A first-principles thermodynamic study was performed for rhombohedral perovskite (La1-xSrx) MnO3±δ (LSM) to assess its stability against Sr secondary phase formation in SOC applications. In this work, the Sr secondary phase formation reaction free energies were determined by combining ab initio lattice dynamics calculations for the solid phases and an ab initio thermodynamics approach for the gas phases. Furthermore, this approach goes beyond previous thermodynamic modeling studies by integrating first-principles based point-defect equilibria into the analyses. The modeling results indicate an increased tendency to form SrO oxide from LSM upon decreasing the oxygen partial pressure. Additionally, enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with contaminant gas species as a function of temperature and gas pressure.

Defect and phase stability↗
DOE OSTI2025

Thermodynamic Profiling Through ASSIST Observations and TROPoe Retrievals

This report reviews the most relevant theoretical aspects of thermodynamic profiling techniques based on spectral observations from ASSIST-II infrared radiometers and TROPoe retrievals. The ASSIST+TROPoe system is a cutting-edge remote sensing technology deployed during the AWAKEN and WFIP3 field campaigns to estimate high-frequency profiles of temperature and humidity in the atmosphere. These profiles are highly valuable for characterizing atmospheric stratification, improving wind models, and understanding the impacts of wind plants on the climate. In this document, we discuss the operating principles of ASSIST, the physics of atmospheric infrared radiation, and the mathematical framework and capabilities of TROPoe. Sources of uncertainties in both the instrument and the retrieval method are also thoroughly addressed. This guide is designed to help users of ASSIST, TROPoe, and thermodynamic data in collecting, estimating, and applying thermodynamic profiles rigorously and scientifically.

17 WIND ENERGY↗
DOE OSTI2025

Aerosol Thermodynamics (Aerosol Liquid Water Concentration and Aerosol pH) at S2 during CoURAGE

These data report aerosol liquid water content (ALWC) and aerosol pH for the S2 (Mt. Airy) site during CoURAGE for April - June 2025. The ISORROPIA-II aerosol thermodynamic equilibrium model was used to compute ALWC and pH. The model was run in forward mode, with solids formation disabled (metastable mode), according to Pye et al. (Atmospheric Chemistry and Physics, 2020). Information about the model can be found at https://www.epfl.ch/labs/lapi/models-and-software/isorropia/, and in Fountoukis and Nenes (Atmospheric Chemistry and Physics, 2007). Inputs to the model were: (1) measured Temperature and relative humidity, both from the DoE ARM “metwxt” data products at S2, (2) aerosol chemical composition from the ARM ACSM, and (3) gas-phase ammonia concentrations. For ISORROPIA-II, ambient RH values above 0.995 were input to the model as 0.995; therefore, caution should be exercised interpreting ALWC when RH was above 0.995. Aerosol composition or NH3(g) concentrations below the LOD were input into the model as 0.5*LOD. Model runs were only conducted for data points at 30-min resolution with all three inputs (meteorology, ACSM, and NH3). If any of the three model inputs were missing, aerosol thermodynamic outputs are flagged as -999. Similarly, aerosol thermodynamic outputs are flagged as -999 if measured NH3(g) and aerosol ammonium (NH4+) were simultaneously below the respective LODs. The ALWC includes liquid water from inorganic aerosol components, calculated directly by ISORROPIA-II, as well as liquid water from organics, estimated from the ACSM organics assuming a kappa value of 0.12 according to Guo et al. (Atmospheric Chemistry and Physics, 2015).

Aerosol Liquid Water Content↗
NASA NTRS1973

Thermodynamic Properties of Nitrogen Including Liquid and Vapor Phases from 63K to 2000K with Pressures to 10,000 Bar

Tables of thermodynamic properties of nitrogen are presented for the liquid and vapor phases for temperatures from the freezing line to 2000K and pressures to 10,000 bar. The tables include values of density, internal energy, enthalpy, entropy, isochoric heat capacity, isobaric heat capacity velocity of sound, the isotherm derivative, and the isochor derivative. The thermodynamic property tables are based on an equation of state, P=P (p,T), which accurately represents liquid and gaseous nitrogen for the range of pressures and temperatures covered by the tables. Comparisons of property values calculated from the equation of state with measured values for P-p-T, heat capacity, enthalpy, latent heat, and velocity of sound are included to illustrate the agreement between the experimental data and the tables of properties presented here. The coefficients of the equation of state were determined by a weighted least squares fit to selected P-p-T data and, simultaneously, to isochoric heat capacity data determined by corresponding states analysis from oxygen data, and to data which define the phase equilibrium criteria for the saturated liquid and the saturated vapor. The vapor pressure equation, melting curve equation, and an equation to represent the ideal gas heat capacity are also presented. Estimates of the accuracy of the equation of state, the vapor pressure equation, and the ideal gas heat capacity equation are given. The equation of state, derivatives of the equation, and the integral functions for calculating derived thermodynamic properties are included.

Thermodynamic properties↗
NASA NTRS2020

A Machine Learning-Based Cloud Detection and Thermodynamic Phase Classification Algorithm using Passive Spectral Observations

We trained two Random Forest (RF) machine-learning models for cloud mask and cloud thermodynamic phase detection using spectral observations from VIIRS on Suomi NPP (SNPP). Observations from CALIOP were carefully selected to provide reference labels. The two RF models were trained for all-day and daytime-only conditions using a 4-year collocated VIIRS/CALIOP dataset from 2013 to 2016. Due to the orbit difference, the collocated CALIOP and SNPP VIIRS training samples cover a broad viewing zenith angle range, which is a great benefit to overall model performance. The all-day model uses 3 VIIRS infrared (IR) bands (8.6,11, and 12 μm) and the daytime model uses 5 Near-IR (NIR) and Shortwave-IR (SWIR) bands (0.86, 1.24, 1.38, 1.64 and 2.25 μm) together with the 3 IR bands to detect clear, liquid water, and ice cloud pixels. Up to 7 surface types, namely, ocean/water, forest, cropland, grassland, snow/ice, barren/desert, and shrubland, were considered separately to enhance performance for both models. Detection of cloudy pixels and thermodynamic phase with the two RF models were compared against collocated CALIOP products from 2017. It is shown that, with a conservative screening process that excludes the most challenging cloudy pixels for passive remote sensing, the two RF models have high accuracy rates in comparison with the CALIOP reference for both cloud detection and thermodynamic phase. Other existing SNPP VIIRS and Aqua MODIS cloud mask and phase products are also evaluated, with results showing that the two RF models and the MODIS MYD06 optical property phase product are the top 3 algorithms with respect to lidar observations during the daytime. During the nighttime, the RF all-day model works best for both cloud detection and phase, in particular for pixels over snow/ice surfaces. The present RF models can be extended to other similar passive instruments if training samples can be collected from CALIOP or other lidars. However, the quality of reference labels and potential sampling issues that may impact model performance would need further attention.

cloud detection↗
NASA NTRS2021

Machine Learning the COSMO Model for Predicting Thermodynamics of Electrolyte Mixtures

Bottom-up design of electrolyte mixtures for battery systems requires predicting macro thermodynamic properties from molecular constituents. For instance, molten salt electrolyte batteries require conditions far above room temperature to operate. Therefore, discovering mixtures with increasingly lower eutectic melting points is desirable. A model that can approximate chemical activity is a valuable tool to search through the vast compositional design space. Machine learning can predict properties of materials such as vibrational free energies, electronic energy gaps, and thermal conductivities. Moreover, they can learn physical models such as interatomic potentials. The COSMO-SAC model uses theory and empirical parameterization to predict liquid-vapor and liquid-solid properties using first-principles calculations. However, obtaining activity coefficients required for parameterizing the COSMO-SAC model is costly and limited to a select chemical space. In this work, we explored if machine learning methods could improve the COSMO-SAC model and bridge density functional theory calculations to liquid phase thermodynamic properties. Our data-driven approach uses existing databases for sigma-profiles of organic solvents and reconciles their methodological differences via ensemble averaging. First, an optimal machine learning model is constructed for each dataset. Our machine learning algorithms use the sigma-profile as an input feature to predict binary mixtures' activity coefficients using multi-output regression. Each dataset uses different choices of functionals, methods, and basis sets. Therefore, our ensemble model attempts to predict corrected activity coefficients given the combination of all the model outputs. The activity coefficients used for training are generated using the COSMO-SAC model. This approach enables the extraction of meaningful information from the existing datasets to improve the COSMO-SAC model for obtaining thermodynamic properties of electrolyte mixtures. With the liquid phase activities, we can identify electrolyte mixtures that meet desired phase equilibria conditions.

Thermodynamics↗
NASA NTRS2021

Evaluation and Applications of Multi-Instrument Boundary-Layer Thermodynamic Retrievals

Recent reports have highlighted the need for improved observations of the boundary layer. In this study, we explore the combination of ground-based active and passive remote sensors deployed for thermodynamic profiling to analyze various boundary-layer observation strategies. Optimal-estimation retrievals of thermodynamic profiles from Atmospheric Emitted Radiance Interferometer (AERI) observed spectral radiance are compared with and without the addition of active sensor observations from a May–June 2017 observation period at the Atmospheric Radiation Measurement–Southern Great Plains Site. In all, three separate thermodynamic retrievals are considered here: retrievals including AERI data only, retrievals including AERI data and Vaisala water vapour differential absorption lidar data, and retrievals including AERI data and Raman lidar data. First, the three retrievals are compared to each other and to reference radiosonde data over the full observation period to get a bulk understanding of their differences and characterize the impact of clouds on these retrieved profiles. These analyses show that the most significant differences are in the water vapour field, where the active sensors are better able to represent the moisture gradient in the entrainment zone near boundary layer top. We also explore how differences in retrievals may impact results of applied analyses including land–atmosphere coupling, convection indices, and severe storm environmental characterization. Overall, adding active sensors to the optimal-estimation retrieval showed some added information, particularly in the moisture field. Given the costs of such platforms, the value of that added information must be weighed for the application at hand.

boundary-layer observation↗
DOE OSTI2025

Surface Nanostructure Control and Thermodynamic Stability Analysis of Femtosecond Laser-Ablated CuCoMn 1.75 NiFe 0.25 Nanoparticles

Surface nanostructure control is the key to functionalizing nanomaterials. This paper presents a characterization with thermodynamic stability analysis of CuCoMn 1.75 NiFe 0.25 high-entropy alloy (HEA) nanoparticles synthesized by femtosecond laser ablation in ethanol and liquid nitrogen (LN2). Using multimodal electron microscopy and spectroscopy, we examine phase, particle size, defect structure, chemical distribution, and surface composition and relate them to HEA stability. Elemental distributions are uniform in both media, but LN2 produces smaller particles with a narrower size distribution and mainly single- or few-domain interiors, whereas ethanol yields larger particles built from 2–4 nm crystallites with domain aggregation. Edge defects appear in both but energy-dispersive X-ray spectroscopy (EDS) is broadly uniform with local fluctuations in ethanol. X-ray photoelectron spectroscopy (XPS), supported by an attenuation model, indicates an ∼1 nm oxide overlayer that suppresses Mn 2p intensity; correcting for it returns Mn toward the bulk value. UV–NIR and photoluminescent spectra independently support a thin oxide shell. Composition-based thermodynamic descriptors place LN2 closer to bulk mixing parameters, while ethanol raises ΔH_mix and lowers Ω. Cooling simulations are consistent (LN2 ∼ 0.1 μs quench, ethanol ∼1 μs). In conclusion, these results connect solvent-controlled kinetics and thermodynamics to crystalline state and surface chemistry, informing surface control of HEA nanoparticles.

Femtosecond Laser Ablation↗
DOE OSTI2026

Thermodynamic Control of Interface Directs MnO 2 Nucleation Chemistry for Dense and Conformal Electrodeposition

Manganese dioxide (MnO 2 ) is widely recognized as a promising material for high-energy-density energy storage systems due to its broad applicability and facile electrodeposition. However, achieving uniform, thin, and high-mass-loading MnO 2 coatings on high-surface-area electrodes remains a significant challenge. Conventional electrodeposition methods typically yield nonuniform, thick layers with poor conductivity and limited material utilization, restricting their practical use. Here, we uncover a thermodynamically engineered vanadyl/pervanadyl (VO 2+ /VO 2 + ) interface that fundamentally reshapes MnO 2 electrodeposition chemistry, enabling highly uniform and dense coatings. Here, combining in situ AFM measurement, Classical Nucleation Theory, and Johnson–Mehl–Avrami–Kolmogorov modeling, we show that this interface reduces early-stage detectable MnO 2 island size by 35-fold and shifts the MnO 2 growth from diffusion-limited to reaction-limited progressive nucleation. This thermodynamically controlled interface yields highly dense and conformal MnO 2 films with record-high mass loading of 241 mg cm –2 (1607 mg cm –3 ) on 3D-printed graphene aerogels, without compromising porosity or inducing thickness gradient. As a prototype demonstration, the resulting MnO 2 electrodes deliver record-setting volumetric performance in both capacitors (106 F cm –3 ) and Zn//MnO 2 pouch cells (162 mAh cm –3 ). Beyond energy storage, our findings demonstrate the significance of thermodynamic interface control in MnO 2 nucleation chemistry for achieving dense and uniform coatings on various substrates, with implications for electrocatalysis, semiconductor processing, and advanced materials manufacturing.

Batteries↗
DOE OSTI2025

Circuit complexity and functionality: A statistical thermodynamics perspective

Circuit complexity, defined as the minimum circuit size required for implementing a particular Boolean computation, is a foundational concept in computer science. Determining circuit complexity is believed to be a hard computational problem. Recently, in the context of black holes, circuit complexity has been promoted to a physical property, wherein the growth of complexity is reflected in the time evolution of the Einstein-Rosen bridge (“wormhole”) connecting the two sides of an anti-de Sitter “eternal” black hole. Here, we are motivated by an independent set of considerations and explore links between complexity and thermodynamics for functionally equivalent circuits, making the physics-inspired approach relevant to real computational problems, for which functionality is the key element of interest. In particular, our thermodynamic framework provides an alternative perspective on the obfuscation of programs of arbitrary length—an important problem in cryptography—as thermalization through recursive mixing of neighboring sections of a circuit, which can be viewed as the mixing of two containers with “gases of gates.” This recursive process equilibrates the average complexity and leads to the saturation of the circuit entropy, while preserving functionality of the overall circuit. The thermodynamic arguments hinge on ergodicity in the space of circuits which we conjecture is limited to disconnected ergodic sectors due to fragmentation. The notion of fragmentation has important implications for the problem of circuit obfuscation as it implies that there are circuits of same size and functionality that cannot be connected via a polynomial number of local moves. Furthermore, we argue that fragmentation is unavoidable unless the complexity classes NP and coNP coincide, a statement that implies the collapse of the polynomial hierarchy of computational complexity theory to its first level.

Science & Technology - Other Topics↗