SEARCH · Engineering Papers

Results for “thermodynamic entropy”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Predicting phase behavior in high entropy and chemically complex alloys

The interest in high entropy alloys and other metallic compounds with four or more elements at near-equiatomic ratios has drawn attention to the ability to rapidly predict phase behavior of these complex materials, particularly where existing thermodynamic data are lacking. In this paper, we discuss aspects of this from the point of view of predicting without utilizing (or fitting) experimental data. Of particular interest are heuristic approaches that provide prediction of single-phase compositions, more rigorous approaches that tackle the thermodynamics from a more fundamental point of view, and simulation approaches that provide further insight into the behaviors. This paper covers cases of all three of these, in order to examine the strengths and weaknesses of each approach, and to indicate directions where these may be utilized and improved upon. Of particular interest is moving beyond “which composition may form a solid solution,” to recognizing the importance of underlying thermodynamic realities that affect the temperature- and composition-dependent transformations of these materials.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Localization dynamics in a centrally coupled system

In systems in which interactions couple a central degree of freedom and a bath, one would expect signatures of the bath's phase to be reflected in the dynamics of the central degree of freedom. This has been recently explored in connection with many-body localized baths coupled with a central qubit or a single-cavity mode - systems with growing experimental relevance in various platforms. Such models also have an interesting connection with Floquet many-body localization via quantizing the external drive, although this has been relatively unexplored. Here we adapt the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) method, a well-known tree tensor network algorithm, to numerically simulate the dynamics of a central degree of freedom, represented by a d-level system (qudit), coupled to a disordered interacting one-dimensional spin bath. ML-MCTDH allows us to reach ≈10 2 lattice sites, a far larger system size than what is feasible with exact diagonalization or kernel polynomial methods. From the intermediate time dynamics, we find a well-defined thermodynamic limit for the qudit dynamics upon appropriate rescaling of the system-bath coupling. The spin system shows similar scaling collapse in the Edward-Anderson spin-glass order parameter or entanglement entropy at relatively short times. At longer timescales, we see slow growth of the entanglement, which may arise from dephasing mechanisms in the localized system or long-range interactions mediated by the central degree of freedom. Similar signs of localization are shown to appear as well with unscaled system-bath coupling.

1-dimensional spin chains↗

NASA NTRS1976

Equilibrium properties of chemically reacting gases

The equilibrium energy, enthalpy, entropy, specific heat at constant volume and constant pressure, and the equation of state of the gas are all derived for chemically reacting gas mixtures in terms of the compressibility, the mol fractions, the thermodynamic properties of the pure gas components, and the change in zero point energy due to reaction. Results are illustrated for a simple diatomic dissociation reaction and nitrogen is used as an example. Next, a gas mixture resulting from combined diatomic dissociation and atomic ionization reactions is treated and, again, nitrogen is used as an example. A short discussion is given of the additional complexities involved when precise solutions for high-temperature air are desired, including effects caused by NO produced in shuffle reactions and by other trace species formed from CO2, H2O and Ar found in normal air.

Source record↗

DOE OSTI2023

Ethanol-fueled metal supported solid oxide fuel cells with a high entropy alloy internal reforming catalyst

High-performance metal supported solid oxide fuel cells (MS-SOFC) with an integrated high entropy alloy (HEA) internal reforming catalyst (IRC) are demonstrated for transportation applications using ethanol and methanol as fuels. Addition of the HEA IRC dramatically improves cell performance and stability when using ethanol/water blend fuel. Absence of carbon deposition predicted by thermodynamic calculations is confirmed by Raman spectroscopy analysis of posttest anodes. Optimal catalyst processing (deposition technique, loading, firing temperature) and cell operation conditions (flow rates, temperature, fuel compositions) are explored. Infiltrated HEA reforming catalyst provides a highly porous structure and low catalyst loading (6 mg cm –2 ). The designed structure and catalysts achieve small mass transport resistances in the fuel electrode (26.2 s m –1 ) and oxygen electrode (41.6 s m –1 ). The best ethanol concentration (60:40 v% ethanol: water) provides 0.83 W cm –1 at 700 °C, without carbon deposition. The ethanol-fueled MS-SOFC is operated for 500 h, including five thermal cycles. As a result, cell evolution is similar to that reported previously for hydrogen fuel; nickel aggregation and chromia deposition were the major observed changes, and carbon formation can be avoided even after long-term operation.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2022

Reversible Ligand Protonation in Noninnocent Constrained-Geometry-Like Group 4 Complexes

The proligands C 5 H 5 CMe 2 CHPhOx R (R = Me 2 , CHMe 2 , and CMe 3 ) react with M(NMe 2 ) 4 (M = Ti, Zr, and Hf) to give monodeprotonated (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 , doubly deprotonated (C 5 H 4 CMe 2 CPhOx R )M(NMe 2 ) 2 , or a mixture of both. The observed products depend on reaction conditions, the oxazoline substituent, and the metal center, with 4-t-butyl-oxazoline or Ti giving only the monodeprotonated (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 . Amine elimination from (Ox R CHPhCMe 2 C 5 H 4 )M(NMe 2 ) 3 to give (C 5 H 4 CMe 2 CPhOx R )M(NMe 2 ) 2 is reversible for 4,4-dimethyl- and 4-isopropyl-oxazoline-based ligands in Zr or Hf complexes. Temperature-dependent kinetic studies of the equilibration of (C 5 H 4 CMe 2 CPhOx Me2 )Zr(NMe 2 ) 2 , HNMe 2 and (Ox Me2 CHPhCMe 2 C 5 H 4 )Zr(NMe 2 ) 3 provide the experimental thermodynamic parameters ΔS° = -17.4 ± 2.6 cal·mol –1 K –1 and ΔH° = -6.8 ± 0.8 kcal·mol –1 . An Eyring plot of the rate constants, determined from the system as it approaches equilibrium, gives the activation entropy and activation enthalpy for the addition of (C 5 H 4 CMe 2 CPhOx Me2 )Zr(NMe 2 ) 2 and HNMe 2 as -36 ± 4 cal·mol –1 K –1 and 8.9 ± 1.3 kcal·mol –1 , respectively; the elimination of HNMe 2 from (Ox Me2 CHPhCMe 2 C 5 H 4 )Zr(NMe 2 ) 3 is characterized by ΔS ‡ = -19 ± 5 cal·mol –1 K –1 and ΔH ‡ = 15.7 ± 1.5 kcal·mol –1 . DFT computational models indicate a single-step, nonlinear transfer of the H between the benzylic position of the noninnocent, oxazoline-coordinated ligand and NMe 2 . Computations also confirm the negative activation entropy and the trends in the barriers support the experimental results. Together, these studies indicate the importance of steric effects from the oxazoline ligand, metal center, ancillary ligands, and leaving group on the shuttling of the proton between HNMe 2 and the noninnocent ligand. Finally, these effects suggest that coordination of oxazoline to the metal center is a key part of the benzylic deprotonation and noninnocent behavior of the cyclopentadienyl-oxazoline ligand.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

New High-Entropy Perovskite Oxides with Increased Reducibility and Stability for Thermochemical Hydrogen Generation

This project aims to design, synthesize, and test a transformative class of High-Entropy Perovskite Oxides (HEPOs) as redox oxides to enable thermochemical hydrogen generation with improved stability, kinetics, and efficiency. These developed HEPOs are expected to demonstrate improved kinetics with oxygen surface exchange coefficient ( k > 7.5×10 -4 cm/s) in Budget Period (BP) 1, retain its structural stability in a broad range of oxygen non-stoichiometry (Δδ > 0.15) at a low operating reduction temperature of T red < 1400°C in BP 2, and deliver a H 2 yield of over 400 µmol per gram of oxide and high stability with less than 20% degradation after at least 50 cycles in BP 3. This project is feasible due to the unique thermodynamic properties (simultaneously increased reducibility and phase stability) and kinetic characters (stability against particle coarsening and potentially enhanced oxygen transport and surface reaction kinetics) of such HEPOs, and it is enabled by a unique active learning computational design approach. Computational studies have been conducted to investigate the oxygen vacancy formation in complex perovskite systems. * Accurate prediction of V O .. concentration with disordered A-site cations in Fe-based high-entropy perovskite oxides * Combined MC/DFT computation elucidates the mechanism of Co preference on the redox due to the strain introduced by local distortion. In this project, we explored a large number (~150) of perovskite compositions, which are listed in Tables 2 – 4). * All perovskite specimens have been synthesized through a high-throughput high-energy ball milling process, followed by the conventional sintering process. * XRD, SEM/EDS and TGA were performed to confirm the crystal structure, phase homogeneity and oxygen non-stoichiometry for compositionally complex perovskite oxides (CCPOs). * 110 compositions show single-phase from XRD * Unusual aliovalent doping effects in medium-entropy perovskite compositions. * V-shape relation between Δδ vs. x (La 1-x Sr x )(Mn 1/3 Fe 1/3 Ti 1/3 )O 3 * The highest reported hydrogen production for the CCPOs made in this project ( T re = 1350 ºC 30 min, T Ox = 1100 ºC 30 min) * B-site mixing (La 0.8 Sr 0.2 )(Mn 0.2 Fe 0.2 Co 0.4 Al 0.2 )O 3 : 89.97 ± 2.73 mmol H2 /mol oxide (395 ± 10 μmol/g oxide ) (i) No phase transformation during reactions when Co molar ratio is less than 61% (ii) Balance between intrinsic kinetics (oxygen surface exchange) and thermodynamics (oxygen non-stoichiometry) (iii) Preference of Co identified by in-situ XPS * A-site mixing (La 1/6 Pr 1/6 Nd 1/6 Gd 1/6 Ba 1/6 Sr1/6)MnO 3 : 98.48 mmol H2 /mol oxid e (~415 μmol/g oxide )

08 HYDROGEN↗

DOE OSTI2020

The concept of spin ice graphs and a field theory for their charges

Originally detected in rare earth pyrochlores, spin ice physics is now being artificially extended to a variety of geometries that control collective behavior and exotic properties, making graph theory their proper framework. We relate spin ice notions, such as ice rule, ice manifold, Coulomb phases, charges, and monopoles, to graph-theoretical notions, such as balance, in/out-degrees, and Euler paths. We then propose a field-theoretical treatment in which topological charges and monopoles are the degrees of freedom, while the binary spins are subsumed in an entropic interaction among charges. We show that for a spin ice on a graph in a Gaussian approximation, the kernel of the entropic interaction is the inverse of the graph Laplacian, and we compute screening functions from the graph spectra as Green operators for the screened Poisson problem on a graph. We then apply the treatment to star graphs, tournaments, cycles, and regular spin ice in different dimensions. Our aim is twofold: to set spin ice physics in a proper graph setting, where only topological rather than geometrical notions hold, and to invite graph theorists to contribute their powerful tools to the field of spin ice.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Atomic-Level, Energy-Conversion Heat Transfer

Abstract Heat is stored in quanta of kinetic and potential energies in matter. The temperature represents the equilibrium and excited occupation (boson) of these energy conditions. Temporal and spatial temperature variations and heat transfer are associated with the kinetics of these equilibrium excitations. During energy-conversion (between electron and phonon systems), the occupancies deviate from equilibria, while holding atomic-scale, inelastic spectral energy transfer kinetics. Heat transfer physics reaches nonequilibrium energy excitations and kinetics among the principal carriers, phonon, electron (and holes and ions), fluid particle, and photon. This allows atomic-level tailoring of energetic materials and energy-conversion processes and their efficiencies. For example, modern thermal-electric harvesters have transformed broad-spectrum, high-entropy heat into a narrow spectrum of low-entropy emissions to efficiently generate thermal electricity. Phonoelectricity, in contrast, intervenes before a low-entropy population of nonequilibrium optical phonons becomes a high-entropy heat. In particular, the suggested phonovoltaic cell generates phonoelectricity by employing the nonequilibrium, low-entropy, and elevated temperature optical-phonon produced population—for example, by relaxing electrons, excited by an electric field. A phonovoltaic material has an ultranarrow electronic bandgap, such that the hot optical-phonon population can relax by producing electron-hole pairs (and power) instead of multiple acoustic phonons (and entropy). Examples of these quanta and spectral heat transfer are reviewed, contemplating a prospect for education and research in this field.

Engineering↗

NASA NTRS1973

WASP: A flexible FORTRAN 4 computer code for calculating water and steam properties

A FORTRAN 4 subprogram, WASP, was developed to calculate the thermodynamic and transport properties of water and steam. The temperature range is from the triple point to 1750 K, and the pressure range is from 0.1 to 100 MN/m2 (1 to 1000 bars) for the thermodynamic properties and to 50 MN/m2 (500 bars) for thermal conductivity and to 80 MN/m2 (800 bars) for viscosity. WASP accepts any two of pressure, temperature, and density as input conditions. In addition, pressure and either entropy or enthalpy are also allowable input variables. This flexibility is especially useful in cycle analysis. The properties available in any combination as output include temperature, density, pressure, entropy, enthalpy, specific heats, sonic velocity, viscosity, thermal conductivity, surface tension, and the Laplace constant. The subroutine structure is modular so that the user can choose only those subroutines necessary to his calculations. Metastable calculations can also be made by using WASP.

Hendricks, R. C.↗

NASA NTRS2007

Microphysical Timescales in Clouds and their Application in Cloud-Resolving Modeling

Independent prognostic variables in cloud-resolving modeling are chosen on the basis of the analysis of microphysical timescales in clouds versus a time step for numerical integration. Two of them are the moist entropy and the total mixing ratio of airborne water with no contributions from precipitating particles. As a result, temperature can be diagnosed easily from those prognostic variables, and cloud microphysics be separated (or modularized) from moist thermodynamics. Numerical comparison experiments show that those prognostic variables can work well while a large time step (e.g., 10 s) is used for numerical integration.

Zeng, Xiping↗

DOE OSTI2021

Asymmetric temperature equilibration with heat flow from cold to hot in a quantum thermodynamic system

A model computational quantum thermodynamic network is constructed with two variable temperature baths coupled by a linker system, with an asymmetry in the coupling of the linker to the two baths. It is found in computational simulations that the baths come to “thermal equilibrium” at different bath energies and temperatures. In a sense, heat is observed to flow from cold to hot. Additionally, a description is given in which a recently defined quantum entropy S univ Q for a pure state “universe” continues to increase after passing through the classical equilibrium point of equal temperatures, reaching a maximum at the asymmetric equilibrium. Thus, a second law account Δ S univ Q ≥ 0 holds for the asymmetric quantum process. In contrast, a von Neumann entropy description fails to uphold the entropy law, with a maximum near when the two temperatures are equal, then a decrease Δ S v N < 0 on the way to the asymmetric equilibrium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Role of water model on ion dissociation at ambient conditions

In this work, we study ion pair dissociation in water at ambient conditions using a combination of classical and ab initio approaches. The goal of this study is to disentangle the sources of discrepancy observed in computed potentials of mean force. In particular, we aim to understand why some models favor the stability of solvent-separated ion pairs vs contact ion pairs. We found that some observed differences can be explained by non-converged simulation parameters. However, we also unveil that for some models, small changes in the solution density can have significant effects on modifying the equilibrium balance between the two configurations. We conclude that the thermodynamic stability of contact and solvent-separated ion pairs is very sensitive to the dielectric properties of the underlying simulation model. In general, classical models are very robust in providing a similar estimation of the contact ion pair stability, while this is much more variable in density functional theory-based models. The barrier to transition from the solvent-separated to contact ion pair is fundamentally dependent on the balance between electrostatic potential energy and entropy. This reflects the importance of water intra- and inter-molecular polarizability in obtaining an accurate description of the screened ion–ion interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Solvent Reorganization and Surface–Species Interactions Control Transfer Hydrogenation Rates in Sn-Beta Zeolites

Reorganization of confined alcohol networks reshapes adsorption and activation thermodynamics in Lewis acid zeolites, yet the roles of adsorbate sterics and non-H-bonding cosolvents remain unclear. Here, we investigate transfer hydrogenation of methyl- (mCH) and tert-butylcyclohexanone (tbCH) with 2-butanol over hydrophobic Sn-Beta. First-order rates per defect-open Sn site are 2–5× higher for mCH than tbCH, but converge in the zero-order regime (333–393 K), indicating that solvent reorganization entropically drives adsorption, while bond activation is insensitive to ketone structure. Reactions in toluene or acetonitrile (with 1 M 2-butanol) show that intraporous solvent organization governs transition-state stabilization: on hydrophobic Sn-Beta, apparent enthalpies and entropies rise with decreasing solvent polarity, whereas hydrophilic zeolites display nearly constant barriers across solvents. Liquid-phase IR spectroscopy corroborates these trends, showing that hydrophilic pores enforce 2-butanol-like solvation environments even in toluene. Together, these results reveal how zeolite pores regulate the intraporous solvent structure, whose reorganization upon adsorption reshapes the thermodynamic landscape of Lewis acid catalysis

Adsorption↗

DOE OSTI2024

Thermal quasiparticle theory

The widely used thermal Hartree–Fock (HF) theory is generalized to include the effect of electron correlation while maintaining its quasi-independent-particle framework. An electron-correlated internal energy (or grand potential) is postulated in consultation with the second-order finite-temperature many-body perturbation theory (MBPT), which then dictates the corresponding thermal orbital (quasiparticle) energies in such a way that all fundamental thermodynamic relations are obeyed. The associated density matrix is of a one-electron type, whose diagonal elements take the form of the Fermi–Dirac distribution functions, when the grand potential is minimized. The formulas for the entropy and chemical potential are unchanged from those of Fermi–Dirac or thermal HF theory. The theory thus stipulates a finite-temperature extension of the second-order Dyson self-energy of one-particle many-body Green’s function theory and can be viewed as a second-order, diagonal, frequency-independent, thermal inverse Dyson equation. At low temperatures, the theory approaches finite-temperature MBPT of the same order, but it may outperform the latter at intermediate temperatures by including additional electron-correlation effects through orbital energies. Here, a physical meaning of these thermal orbital energies is proposed (encompassing that of thermal HF orbital energies, which has been elusive) as a finite-temperature version of Janak’s theorem.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1974

Tables and charts of equilibrium normal shock and shock tube solutions for pure CO2 with velocities to 16 km/second

Equilibrium thermodynamic and flow properties are presented in tabulated and graphical form for moving, standing, and reflected normal shock waves in pure CO2, representative of Mars and Venus atmospheres. Properties include pressure, temperature, density, enthalpy, speed of sound, entropy, molecular weight ratio, isentropic exponent, velocity and species mole fractions. Incident (moving) shock velocities are varied from 1 to 16 km/sec for a range of initial pressure of 5 Newtons per square meter to 500 kilo Newtons per square meter. The present results are applicable to shock tube flows, and to free-flight conditions for a blunt body at high velocities. Working charts illustrating idealized shock-tube performance with CO2 test gas and heated helium and hydrogen driver gases are also presented.

Miller, C. G., III↗

NASA NTRS2016

Isentropic Analysis of a Simulated Hurricane

Hurricanes, like many other atmospheric flows, are associated with turbulent motions over a wide range of scales. Here the authors adapt a new technique based on the isentropic analysis of convective motions to study the thermodynamic structure of the overturning circulation in hurricane simulations. This approach separates the vertical mass transport in terms of the equivalent potential temperature of air parcels. In doing so, one separates the rising air parcels at high entropy from the subsiding air at low entropy. This technique filters out oscillatory motions associated with gravity waves and separates convective overturning from the secondary circulation. This approach is applied here to study the flow of an idealized hurricane simulation with the Weather Research and Forecasting (WRF) Model. The isentropic circulation for a hurricane exhibits similar characteristics to that of moist convection, with a maximum mass transport near the surface associated with a shallow convection and entrainment. There are also important differences. For instance, ascent in the eyewall can be readily identified in the isentropic analysis as an upward mass flux of air with unusually high equivalent potential temperature. The isentropic circulation is further compared here to the Eulerian secondary circulation of the simulated hurricane to show that the mass transport in the isentropic circulation is much larger than the one in secondary circulation. This difference can be directly attributed to the mass transport by convection in the outer rainband and confirms that, even for a strongly organized flow like a hurricane, most of the atmospheric overturning is tied to the smaller scales.

convection↗

DOE OSTI2023

Thermodynamic properties and superconductivity of natural carrollite (CuCo 2 S 4 )

The thermodynamic properties of natural carrollite with a formula Cu 0.92 Co 2.07 S 4 are studied using low temperature heat capacity and high temperature oxidative solution calorimetry. The standard thermodynamic functions at 298.15 K are: C p,m °, S m °, H m °, and Φ m ° are 158.48 J∙K −1 ∙mol −1 , 176.33 J∙K −1 ∙mol −1 , 28.40 kJ∙mol −1 , and 81.07 kJ∙mol −1 with an estimated error of 1%. The enthalpies (Δ f H), entropies (Δ f S), and Gibbs energies of formation (Δ f G) from elements at 298.15 K are: −344.46 ± 12.87 kJ∙mol −1 , −45.22 ± 1.81 J∙K −1 ∙mol −1 , and −331.14 ± 1.89 J∙K −1 ∙mol −1 . The stability of carrollite relative to the elements is demonstrated with a negative Gibbs energy (ΔG r °) of formation for several temperatures between 0 and 300 K. Furthermore, a superconductivity transition, which has been observed previously for synthetic carrollite, is confirmed by both the low temperature heat capacity and magnetization measurements.

Chemistry↗

DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗