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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Systematically derived thermodynamic properties for alkane oxidation

Key combustion properties, such as ignition delay time, show strong sensitivity to the thermochemistry of the main species in the standard radical oxidation pathway at low temperatures (600-1000K). Significant uncertainties persist in current estimates of thermodynamic properties, particularly for the larger species of relevance to common practical fuels. In this article, we use advanced computational schemes to evaluate thermodynamic properties for the fuel, fuel radical, peroxy, hydroperoxy-alkyl radical, and hydroperoxide species for a set of 17 fuels containing up to 9 carbon atoms, and with various degrees of branching in the alkane backbone. The procedure, termed STAR-1D, combines conformer sampling to find the minimum geometry, B2PLYP-D3/cc-pVTZ harmonic frequency evaluations, and omega B97X-D/cc-pVTZ one-dimensional torsional mapping. It includes two physically-based scaling routines: a frequency-dependent scaling to B2PLYP-D3/cc-pVTZ anharmonic frequencies and a scaling of one-dimensional omega B97X-D/ccpVTZ torsional profiles to reproduce the product of the B2PLYP-D3/cc-pVTZ frequencies in the harmonic limit. Substantive comparisons with existing experimental databases, together with careful examinations of key theoretical assumptions, are used to explore the accuracy of the predictions. These computationally intensive explorations of 195 species were facilitated by automated thermochemistry software. Further, high accuracy 0 K heats of formation from a separate study are used with the STAR-1D computations to generate NASA polynomial representations. In addition to their intrinsic value, the present results also provide a reliable database for the optimization of group additivity or machine learning schemes for scaling to larger combustion systems. Towards this end a complementary extensive conformational analysis is carried out for the medium sized species and the thermodynamic properties of the lowest energy hydrogen-bonded and non-hydrogen bonded conformers are contrasted for larger hydroperoxy-alkyl radical species.

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

Ab initio studies on structural and thermodynamic properties of magnetic Fe

The present work systematically investigates the total energy, phonon spectra, and thermodynamic properties of different polymorphs of pure Fe, i.e., FCC, BCC, and HCP, with the ab initio approach, considering various magnetic configurations. In general, the calculated energy vs. volume curves and phonon spectra agree well with previous calculations and the experimental data. In addition, their thermodynamic properties are estimated by the quasiharmonic approximation (QHA). Specifically, a superposition approach based on the latest Zentropy theory was utilized to predict magnetic transition temperatures and thermodynamic properties of pure Fe. With the ensemble of the partition function considering the multiplicity of each magnetic microstate, the current work successfully reproduced the Curie/Néel temperature and the Schottky anomaly of heat capacity in FCC, BCC, and HCP Fe purely based on the ab initio input, which exhibits good agreement with the experimental data and CALPHAD modeling.

36 MATERIALS SCIENCE↗
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 OSTI2020

Thermodynamic analysis and optimization of supercritical carbon dioxide Brayton cycles for use with low-grade geothermal heat sources

The transcritical carbon dioxide Rankine cycle has emerged as an alternative for power generation in low-grade heat applications. However, the low heat rejection temperatures required for condensing the carbon dioxide is almost prohibitive for many locations. A simple solution is the adoption of a supercritical carbon dioxide Brayton cycle which does not need such low temperature heat sinks. This power cycle has been poorly studied previously in low-grade heat applications. Therefore, with the aim of determining the performance and feasibility of the supercritical carbon dioxide Brayton cycle from thermodynamic viewpoint, using a low-grade geothermal heat source, this paper presents a comparative thermodynamic analysis between four different supercritical carbon dioxide Brayton cycles. For this purpose, detailed models were developed and solved by coupling the thermodynamic model with the thermal-hydraulic model of a Printed Circuit Heat Exchanger. Subsequently, simulations of a base case, a parametric analysis, and an optimization using Genetic Algorithms, were carried out. Results show that there is a combination of operating conditions that maximizes the electrical output of the system. Besides, some values of the minimum pressure of the cycles can increase substantially the precooler cooling water pumping power. At the optimum design point, using a 20 kg/s geothermal brine stream at 150 °C as heat source with a minimum allowable reinjection temperature of 70 °C, the Intercooled Recuperated Brayton Cycle achieved the highest electric power output, energy and exergy efficiencies, obtaining values of 779.99 kW, 11.51%, and 52.49%, respectively. The other feasible alternatives were the Recuperated Brayton Cycle, the Simple Brayton Cycle and the Intercooled Brayton Cycle, ranked, in that order.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Thermodynamic properties and phase stability of the Ba-Bi system: A combined computational and experimental study

The thermodynamic properties and phase stability of the Ba-Bi system are investigated computationally and experimentally in the present work. The enthalpies of formation and the finite temperature thermodynamic properties of seven compounds are predicted by first-principles calculations based on density functional theory (DFT), indicating five compounds (BaBi 3 , Ba 11 Bi 10 , Ba 4 Bi 3 , Ba 5 Bi 3 , and Ba 2 Bi) to be stable. Phase relations at 773 K and 858 K with composition x Ba = 0.90 are established by isothermal annealing and powder X-ray diffraction (XRD) to clarify the previously observed phase transition at 796 K. The extremely low chemical activity of Ba in liquid for a wide range of temperatures and compositions indicates very strong short-range ordering in the liquid phase which is modeled in the present work by introducing the Ba 4 Bi 3 and BaBi 3 associates in the liquid phase. Both thermodynamic and phase equilibrium data are then used to evaluate the model parameters in Gibbs energy functions of the five stable compounds and three solution phases of liquid, bcc, and rhombohedral phases by the CALPHAD (CALculation of PHAse Diagram) technique. Here the present work shows that the Ba-Bi system consists of three eutectic reactions, two peritectic reactions, one peritectoid reaction, and two congruent reactions, as well as that the concentrations of associates are very high in the liquid phase with very low concentration of atomic Ba, which provides the fundamental understanding as to why Bi can be used to remove Ba ions from molten salt solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Thermodynamical effects and high resolution methods for compressible fluid flows

One of the fundamental differences of compressible fluid flows from incompressible fluid flows is the involvement of thermodynamics. This difference should be manifested in the design of numerical schemes. Unfortunately, the role of entropy, expressing irreversibility, is often neglected even though the entropy inequality, as a conceptual derivative, is verified for some first order schemes. In this paper, we refine the GRP solver to illustrate how the thermodynamical variation is integrated into the design of high resolution methods for compressible fluid flows and demonstrate numerically the importance of thermodynamic effects in the resolution of strong waves. As a by-product, we show that the GRP solver works for generic equations of state, and is independent of technical arguments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Practical thermodynamic quantities for aqueous vanadium- and iron-based flow batteries

A simple method for experimentally determining thermodynamic quantities for flow battery cell reactions is presented. Equilibrium cell potentials, temperature derivatives of cell potential (d E /d T ), Gibbs free energies, and entropies are reported here for all-vanadium, iron–vanadium, and iron–chromium flow cells with state-of-the-art solution compositions. Proof is given that formal potentials and formal temperature coefficients can be used with modified forms of the Nernst Equation to quantify the thermodynamics of flow cell reactions as a function of state-of-charge. Such empirical quantities can be used in thermo-electrochemical models of flow batteries at the cell or system level. In most cases, the thermodynamic quantities measured here are significantly different from standard values reported and used previously in the literature. The data reported here are also useful in the selection of operating temperatures for flow battery systems. Because higher temperatures correspond to lower equilibrium cell potentials for the battery chemistries studied here, it can be beneficial to charge a cell at higher temperature and discharge at lower temperature. As a result, proof-of-concept of improved voltage efficiency with the use of such non-isothermal cycling is given for the all-vanadium redox flow battery, and the effect is shown to be more pronounced at lower current densities.

25 ENERGY STORAGE↗
DOE OSTI2023

Kinetically versus thermodynamically controlled factors governing elementary pathways of GaP(111) surface oxidation

GaP and related III-V semiconductors have attracted interest as high-efficiency photoelectrochemical electrodes for solar water splitting. However, their efficacy is linked to the presence, identity, and integrity of native surface oxides, which are structurally and chemically complex and evolve during operation. Using ambient pressure X-ray photoelectron spectroscopy (APXPS) coupled with ab initio simulations, we track key chemical motifs expressed during evolution of GaP(111) surface oxides, their associated reaction kinetics, and their correlations with electronic properties. We identify two distinct thermal regimes corresponding to kinetically and thermodynamically controlled oxidation. Below 600 K, exposure to O2 generates kinetically facile Ga-O-Ga configurations, whereas higher temperatures cause activated oxygen to insert into Ga–P bonds as part of a thermodynamically driven transformation into a complex, heterogeneous 3D network of surface POx (1≤x ≤ 4) groups and Ga2O3 species, the latter of which eventually dominates upon depletion of surface phosphorus. Here our study highlights the critical competition between kinetic and thermodynamic factors during GaP oxidation, yielding insights for fabricating stable III-P-based photoelectrodes with precisely engineered surface properties.

08 HYDROGEN↗
DOE OSTI2017

Discrete Thermodynamics

Here, we consider the dependence of velocity probability distribution functions on the finite size of a thermodynamic system. We are motivated by applications to computational fluid dynamics, hence discrete thermodynamics. We then begin by describing a coarsening process that represents geometric renormalization. Then, based only on the requirements of conservation, we demonstrate that the pervasive assumption of local thermodynamic equilibrium is not form invariant. We develop a perturbative correction that restores form invariance to second-order in a small parameter associated with macroscopic gradients. Finally, we interpret the corrections in terms of unresolved kinetic energy and discuss the implications of our results both in theory and as applied to numerical simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Bayesian automated weighting of aggregated DFT, MD, and experimental data for candidate thermodynamic models of aluminum with uncertainty quantification

Atomic-scale modeling methods such as density functional theory (DFT) and molecular dynamics (MD) can predict the thermodynamic properties of materials at a lower cost than experimental measurements. However, their regular usage in thermodynamic model construction is hampered by the lack of quantitative agreement with experimental measurements and the lack of uncertainty estimates on the data. To make regular usage of this atomistic simulation data, it is important to assess whether the atomistic simulation datasets, by themselves or in combination with experimental measurements, result in the same physics-informed models best supported by experimental measurements alone. Here, models of aluminum thermodynamic properties are discussed using three data sources: atomistic calculations (DFT and MD), experiments, and a combination of atomistic calculations and experiments. The study shows that, after ensuring self-consistency in predicting key invariant points, both experimental measurements and atomistic calculations can significantly contribute to an optimal model.

36 MATERIALS SCIENCE↗
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 OSTI2017

Thermodynamic Routes to Novel Metastable Nitrogen-Rich Nitrides

Compared to oxides, the nitrides are relatively unexplored, making them a promising chemical space for novel materials discovery. Of particular interest are nitrogen-rich nitrides, which often possess useful semiconducting properties for electronic and optoelectronic applications. However, such nitrogen-rich compounds are generally metastable, and the lack of a guiding theory for their synthesis has limited their exploration. Here, we review the remarkable metastability of observed nitrides, and examine the thermodynamics of how reactive nitrogen precursors can stabilize metastable nitrogen-rich compositions during materials synthesis. We map these thermodynamic strategies onto a predictive computational search, training a data-mined ionic substitution algorithm specifically for nitride discovery, which we combine with grand-canonical DFT-SCAN phase stability calculations to compute stabilizing nitrogen chemical potentials. We identify several new nitrogen-rich binary nitrides for experimental investigation, notably the transition metal nitrides Mn3N4, Cr3N4, V3N4, and Nb3N5, the main group nitride SbN, and the pernitrides FeN2, CrN2, and Cu2N2. By formulating rational thermodynamic routes to metastable compounds, we expand the search space for functional technological materials beyond equilibrium phases and compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman–Rubin Convergence Diagnostics

The determination of the solvation free energy of ions and molecules holds profound importance across a spectrum of applications spanning chemistry, biology, energy storage, and the environment. Molecular dynamics simulations are powerful tools for computing this critical parameter. Nevertheless, the accurate and efficient calculation of the solvation free energy becomes a formidable endeavor when dealing with complex systems characterized by potent Coulombic interactions and sluggish ion dynamics and, consequently, slow transition across various metastable states. Here, in the present study, we expose limitations stemming from the conventional calculation of the statistical inefficiency g in the thermodynamic integration method, a factor that can hinder the determination of convergence of the solvation free energy and its associated uncertainty. Instead, we propose a robust scheme based on Gelman–Rubin convergence diagnostics. We leverage this improved estimation of uncertainties to introduce an innovative accelerated thermodynamic integration method based on the Gaussian Process regression. This methodology is applied to the calculation of the solvation free energy of trivalent rare-earth elements immersed in ionic liquids, a scenario in which the aforementioned challenges render standard approaches ineffective. The proposed method proves to be effective in computing solvation free energy in situations where traditional thermodynamic integration methods fall short.

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

Revealing the Phase Separation Behavior of Thermodynamically Immiscible Elements in a Nanoparticle

Phase-separation is commonly observed in multimetallic nanomaterials, yet it is not well understood how immiscible elements distribute in a thermodynamically stable nanoparticle. Herein, we studied the phase-separation of Au and Rh in nanoparticles using electron microscopy and tomography techniques. The nanoparticles were thermally annealed to form thermodynamically stable structures. HAADF-STEM and EDS characterizations reveal that Au and Rh segregate into two domains while their miscibility is increased. Using aberration-corrected HAADF-STEM and atomic electron tomography, we show that the increased solubility of Au in Rh is achieved by forming Au clusters and single atoms inside the Rh domains and on the Rh surface. Furthermore, based on the three-dimensional reconstruction of a AuRh nanoparticle, we can visualize the uneven interface that is embedded in the nanoparticle. Finally, the results advance our understanding on the nanoscale thermodynamic behavior of metal mixtures, which is crucial for the optimization of multimetallic nanostructures for many applications.

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

Thermodynamic Stability of Low- and High-Index Spinel LiMn 2 O 4 Surface Terminations

Density functional theory calculations are performed within the generalized gradient approximation (GGA+U) to determine stable terminations of both low- and high-index spinel LiMn2O4 (LMO) surfaces. A grand canonical thermodynamic approach is employed, permitting a direct comparison of offstoichiometric surfaces with previously reported stoichiometric surface terminations at various environmental conditions. Within this formalism, we have identified trends in the structure of the low-index surfaces as a function of the Li and O chemical potentials. The results suggest that, under a range of chemical potentials for which bulk LMO is stable, Li/O and Li-rich (111) surface terminations are favored, neither of which adopts an inverse spinel structure in the subsurface region. This thermodynamic analysis is extended to identify stable structures for certain high-index surfaces, including (311), (331), (511), and (531), which constitute simple models for steps or defects that may be present on real LMO particles. The low- and high-index results are combined to determine the relative stability of each surface facet under a range of environmental conditions. The relative surface energies are further employed to predict LMO particle shapes through a Wulff construction approach, which suggests that LMO particles will adopt either an octahedron or a truncated octahedron shape at conditions in which LMO is thermodynamically stable. These results are in agreement with the experimental observations of LMO particle shapes.

DFT↗
DOE OSTI2020

Decoupling Kinetics and Thermodynamics of Interfacial Catalysis at a Chemically Modified Black Silicon Semiconductor Photoelectrode

Understanding the interplay between the kinetics of interfacial catalytic reactions and the thermodynamics of an underlying semiconductor electrode is imperative for rational construction of efficient photoelectrocatalytic systems. Current understanding of the thermodynamic effects of molecular catalyst attachment to semiconductor electrodes is limited. We report the immobilization of a molecular cobalt bis(benzenedithiolate) proton reduction catalyst onto nanoporous black silicon (b-Si) electrodes through p–p interactions with a series of aromatic molecules covalently attached to the surface. Furthermore, intensity-modulated high-frequency resistivity and linear sweep voltammetry measurements are used to show that the kinetics of proton reduction are decoupled from the thermodynamic properties of the underlying b-Si photoelectrode.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗