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 469 records · Page 26

Thermodynamic properties of asymptotically Reissner–Nordström black holes

Motivated by possible relation between Born–Infeld type nonlinear electrodynamics and an effective low-energy action of open string theory, asymptotically Reissner–Nordström black holes whose electric field is described by a nonlinear electrodynamics (NLED) are studied. We take into account a four dimensional topological static black hole ansatz and solve the field equations, exactly, in terms of the NLED as a matter field. The main goal of this paper is investigation of thermodynamic properties of the obtained black holes. Moreover, we calculate the heat capacity and find that the nonlinearity affects the minimum size of stable black holes. We also use Legendre-invariant metric proposed by Quevedo to obtain scalar curvature divergences. We find that the singularities of the Ricci scalar in Geometrothermodynamics (GTD) method take place at the Davies points. -- Highlights: •We examine the thermodynamical properties of black holes in Einstein gravity with nonlinear electrodynamics. •We investigate thermodynamic stability and discuss about the size of stable black holes. •We obtain analytical solutions of higher dimensional theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Thermodynamical properties of graphene in noncommutative phase–space

We investigated the thermodynamic properties of graphene in a noncommutative phase–space in the presence of a constant magnetic field. In particular, we determined the behaviour of the main thermodynamical functions: the Helmholtz free energy, the mean energy, the entropy and the specific heat. The high temperature limit is worked out and the thermodynamic quantities, such as mean energy and specific heat, exhibit the same features as the commutative case. Possible connections with the results already established in the literature are discussed briefly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Thermodynamic modeling of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems

Molten salt electrolysis can recover the actinides from spent nuclear fuels, and it involves a eutectic LiCl-KCl in molten form as an electrolyte. During reprocessing, the concentration of fission products such as La, Nd, Pr, etc., increases and affects the recovery efficiency of the electrolyte. In this work, thermodynamic modeling of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems was carried out using the CALPHAD (Calculation of Phase Diagrams) approach for the first time. The thermodynamic functions for the pure salts were taken from the SGTE (Scientific Group Thermodata Europe) Substances (SSUB) database. The experimental thermochemical and phase equilibria data available in the literature were used as input for the assessment of KCl-PrCl 3 and KCl-LiCl-PrCl 3 systems. The model parameters for the KCl-LiCl system were adjusted to include the new Gibbs energy descriptions for the pure salts. In addition, the sublattice model for the liquid phase in the LiCl-PrCl 3 system was modified and reassessed to ensure the model compatibility for higher-order extrapolation. There is a good agreement between the experimental and calculated thermochemical and phase diagram data for all the optimized constituent binaries and the ternary system. Furthermore, this work will be beneficial for determining the solubility limit of PrCl 3 in molten LiCl-KCl electrolytes and their thermodynamic properties for improving the efficiency of the pyrochemical process.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermodynamic modeling with uncertainty quantification using the modified quasichemical model in quadruplet approximation: Implementation into PyCalphad and ESPEI

The modified quasichemical model in the quadruplet approximation (MQMQA) considers the first- and the second-nearest-neighbor coordination and interactions, particularly useful in describing short-range ordering (SRO) in complex liquids such as molten salts, slag in metal processing, and electrolytic solutions. Here, the present work implements the MQMQA into the Python based open-source software PyCalphad for thermodynamic calculations. This endeavor facilitates the development of MQMQA-based thermodynamic database with uncertainty quantification (UQ) and propagation (UP) using the open-source software ESPEI. A new database structure based on Extensible Markup Language (XML) is proposed for ESPEI evaluation of MQMQA model parameters. Using the KF-NiF 2 , KCl-NaCl-MgCl 2 , and CaCl 2 -CaF 2 -LiCl-LiF salt systems as examples, we demonstrate the successful implementation of MQMQA in PyCalphad through thermodynamic calculations of Gibbs energy, equilibrium quadruplet fractions, and phase diagram, as well as database development with UQ and UP using ESPEI. Furthermore, as an application of the present implementation, both the LiF–TbF 3 and LiF-HoF 3 systems have been modeled by MQMQA for the first time, which are in good agreement with experiments. The present implementation hence offers an open-source capability for performing CALPHAD modeling for complex liquids with SRO using MQMQA plus a new XML database structure.

36 MATERIALS SCIENCE↗

Thermodynamics and kinetics of Pb intercalation under graphene on SiC(0001)

SiC-supported graphene intercalated by a two-dimensional Pb monolayer can provide an appealing platform for spintronic applications. Such a monolayer structure is thermodynamically ultrastable, as observed in recent experiments. However, important fundamentals such as the structure of intercalated phases, locations of intercalated atoms, thermodynamic preference for intercalation, and intercalation pathways for this system have not yet been understood conclusively. Here, in this work, extensive density functional theory calculations are performed to assess Pb intercalation thermodynamics and kinetics under graphene on SiC(0001). We find that intercalation of isolated Pb atoms is strongly disfavored over adsorption on top of graphene. However, intercalation of complete Pb layers in the gallery between SiC and graphene buffer layer is strongly favored over supported Pb monolayers and moderately favored over formation of supported large three-dimensional Pb islands. We also find that initiation of intercalation either by individual Pb atoms directly penetrating graphene or by hopping under a static graphene step edge is energetically prohibitive at experimental temperatures. Consequently, more complex intercalation pathways are operative and further analyzed. We demonstrated that once an intercalated Pb monolayer forms around a graphene step edge, a facile Pb mass transport by Pb vacancy-mediated diffusion can be triggered for continued growth of the intercalated monolayer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Global optimization algorithms to compute thermodynamic equilibria in large complex systems with performance considerations

Several global optimization methods are reviewed that attempt to ensure that the integral Gibbs energy of a closed isothermal isobaric system is a global minimum to satisfy the necessary and sufficient conditions for thermodynamic equilibrium. In particular, the integral Gibbs energy function of a multicomponent system containing non-ideal phases may be highly non-linear and non-convex, which makes finding a global minimum a challenge. Consequently, a poor numerical approach may lead one to the false belief of equilibrium. Furthermore, confirming that one reaches a global minimum and that this is achieved with satisfactory computational performance becomes increasingly more challenging in systems containing many chemical elements and a correspondingly large number of species and phases. Several numerical methods that have been used for this specific purpose are reviewed with a benchmark study of three of the more promising methods using five case studies of varying complexity. A modification of the conventional Branch and Bound method is presented that is well suited to a wide array of thermodynamic applications, including complex phases with many constituents and sublattices, and ionic phases that must adhere to charge neutrality constraints. Also, a novel method is presented that efficiently solves the system of linear equations that exploits the unique structure of the Hessian matrix, which reduces the calculation from a O(N 3 ) operation to a O(N) operation. As a result, this combined approach demonstrates efficiency, reliability and capabilities that are favorable for integration of thermodynamic computations into multi-physics codes with inherent performance considerations.

97 MATHEMATICS AND COMPUTING↗

Thermodynamics of concentrated solid solution alloys

This study reviews the three main approaches for predicting the formation of concentrated solid solution alloys (CSSA) and for modeling their thermodynamic properties, in particular, utilizing the methodologies of empirical thermo-physical parameters, CALPHAD method, and first-principles calculations combined with hybrid Monte Carlo/Molecular Dynamics (MC/MD) simulations. In order to speed up CSSA development, a variety of empirical parameters based on Hume-Rothery rules have been developed. Herein, these parameters have been systematically and critically evaluated for their efficiency in predicting solid solution formation. The phase stability of representative CSSA systems is then illustrated from the perspectives of phase diagrams and nucleation driving force plots of the σ phase using CALPHAD method. The temperature-dependent total entropies of the FCC, BCC, HCP, and σ phases in equimolar compositions of various systems are presented next, followed by the thermodynamic properties of mixing of the BCC phase in Al-containing and Ti-containing refractory metal systems. First-principles calculations on model FCC, BCC and HCP CSSA reveal the presence of both positive and negative vibrational entropies of mixing, while the calculated electronic entropies of mixing are negligible. Temperature dependent configurational entropy is determined from the atomic structures obtained from MC/MD simulations. Current status and challenges in using these methodologies as they pertain to thermodynamic property analysis and CSSA design are discussed.

36 MATERIALS SCIENCE↗

Thermodynamic Modeling of Aqueous LiCl, LiBr, LiI, and LiNO 3 Solutions

Thermodynamic models are essential to facilitate the advancement of process design, optimization, and operation of electrolyte systems. In this work, a comprehensive thermodynamic framework based on the Electrolyte Nonrandom Two-Liquid model is developed to calculate phase equilibria behavior and salt solubility of aqueous LiCl, LiBr, LiI, and LiNO 3 solutions. The model describes the non-ideality of the electrolyte solutions by using two binary interaction parameters for each electrolyte-molecule pair in the system. Each binary interaction parameter is further expressed with up to three temperature coefficients which are regressed from experimental data. To take into account the hydration of lithium ion, two separate chemistries for the dissociation of lithium salts are investigated. In the first case, the lithium ion is considered as a bare ion, Li + , while in the second case hydration of the lithium ion from Li + to Li(H 2 O) + is considered. Here, the calculated thermodynamic properties compare adequately with the experimental data for both sets of chemistries for concentrations up to saturation and temperatures from 273.15 K up to 623.15 K. Moderate to significant improvements are observed with the incorporation of the hydration chemistry for aqueous LiCl, LiBr, and LiI solutions when compared to the non-hydrated lithium ion model results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamic modeling of adsorption at the liquid-solid interface

Adsorptive separation techniques are significantly energy efficient in comparison to conventional thermal separation techniques such as distillation. Despite extensive research and development activities undertaken for mixed-gas adsorption, the use of adsorption techniques for the separation of multicomponent liquid mixtures is still limited. A major barrier is the lack of accurate adsorption thermodynamic models, which form the scientific foundation of process simulation of such systems, making the translation to industrial scale challenging. In this work, we have rigorously computed the surface excess of adsorption for six binary liquid mixtures on silica gel at 303 K using the frameworks of the generalized Langmuir isotherm (gL) and the adsorbed solution theory (AST). The six binary liquid mixtures were formed by the pair-wise combinations of four components: benzene, 1,2-dichloroethane, cyclohexane, and n-heptane. We have based our calculations by considering simultaneous equilibria of three phases: saturated vapor phase, bulk liquid phase, and adsorbed phase. Further, the composition of the corresponding saturated vapor phase was estimated by the Nonrandom Two-Liquid activity coefficient model and experimental vapor-liquid equilibria data. The activity coefficients of the adsorbed phase, the central issue of multicomponent adsorption thermodynamics, were calculated using the adsorption Nonrandom Two-Liquid activity coefficient model. Devoid of simplifying assumptions, gL and AST provide rigorous thermodynamic frameworks for adsorption equilibria of multicomponent liquid mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Heat transfer and pseudo phase transition for low-Reynolds, mixed-convection channel flow in the supercritical thermodynamic regime

Fluids at supercritical thermodynamic conditions are inherently complex due to the large variations in thermodynamic and transport properties. Multiple recent numerical and experimental investigations indicate ongoing interest for these fluids, especially supercritical CO and supercritical water, for a variety of applications. The potential heat transfer benefits in this regime are directly influenced by the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. Limited studies have been conducted for low-Reynolds/intermediate-Rayleigh numbers mixed-convection channel flow of supercritical water, which is the focus of this study. To investigate the thermally driven hydrodynamic instabilities in this regime, we use a high-order fully-implicit numerical method. As expected, in this channel configuration and flow regime, a competition arises between the forced and the natural convection. Buoyancy, driven by the thermal gradient at the wall, pushes the fluid upwards in the channel, while inertial forces drive the fluid from inlet to outlet. We observe significant differences in the heat transfer coefficient between the top and bottom walls of the channel flow, due to natural convection, and unstable plumes. We further quantify the buoyancy forces along the channel flow through evaluation of the Richardson number. Simulations show that with a higher negative Richardson number, there are more unstable structures along the channel, resulting in enhanced mixing. Furthermore, we demonstrate the existence of gas-like, liquid-like, and the possibility of pseudo-two-phase regimes in the supercritical fluid. The pseudo-two-phase regime occurs when the range of fluid temperatures is crossing the pseudo-critical line, also known as the Widom line. The density of the fluid varies significantly near the Widom line, and further influences the natural convection and heat transfer in the flow.

42 ENGINEERING↗

Thermodynamics and kinetics of interconversion between platinum nanoparticles and cations in zeolites

In metal-containing zeolites, sintering and redispersion processes exercise control over the identity of metal site structures, but the thermodynamic and kinetic factors that influence these molecular-level processes are not completely understood. Furthermore, we assess the ability of first principles informed free energy models (for supported and unsupported nanoparticles) and kinetic models integrating Ostwald ripening with atom trapping to describe the interconversion between Pt cations and nanoparticles encapsulated in chabazite (CHA) zeolites. Density functional theory-derived thermodynamic phase diagrams show that the interconversion between cations, favored in oxidizing environments, and particles, favored in reducing environments, is fully reversible within a wide range of their respective conditions (temperatures and pressures) for CHA and several other zeolite topologies. Kinetic Monte Carlo simulations of Pt redispersion are consistent with experimentally observed redispersion kinetics of encapsulated Pt nanoparticles in CHA zeolites, and model results suggest the zeolite host imparts additional stability for Pt nanoparticles. We envision our thermodynamic and kinetic models for Pt-CHA are also capable of describing nanoparticle and cation interconversions for other zeolite frameworks under similar conditions.

36 MATERIALS SCIENCE↗

The low-temperature heat capacity and thermodynamic properties of greigite (Fe 3 S 4 )

Heat capacity measurements provide important insights into the energetic, thermodynamic, and magnetic properties of materials. Herein we report the heat capacity of greigite (Fe 3 S 4 ) from 1.8 to 300 K. Greigite is a magnetic spinel mineral and through a ferromagnetic magnon term, C fsw = B fsw T 3/2 , ferrimagnetic ordering is observed in the low-temperature heat capacity. Using a set of theoretical fits of the experimental data, we calculate the thermodynamic functions, including the standard entropy ($Δ_{0}^{T}$S m °). Greigite is important in iron sulfide formation and reaction pathways in environmental, ore-forming, and technological settings and previous work has measured enthalpies $ΔH_{r}°$ of formation and decomposition to neighboring phases. In this work, the stability of greigite relative to the elements is demonstrated with a negative Gibbs energy ($ΔG_{r}°$) of formation and the stability relative to decomposition products of pyrrhotite (FeS 1.092 ) and pyrite (FeS 2 ) is demonstrated with a positive Gibbs energy ($ΔG_{r}°$) of decomposition. Values of the standard thermodynamic functions C p,m °, $Δ_{0}^{T}$S m °, $Δ_{0}^{T}$H m °, and Φ m ° are tabulated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New perspectives on supercritical methane adsorption in shales and associated thermodynamics

Understanding methane adsorption behavior in shales is fundamental for optimizing shale gas development as the adsorbed methane is a large portion of the subsurface shale gas resource. However, the adsorption mechanism of supercritical methane in shales and associated thermodynamics are poorly understood because the equation of state of the adsorbed methane is unmeasurable. This work analyzed adsorption equilibria (up to 32 MPa and 393.15 K) using a rigorous framework that can account for non-ideal gas properties and accurately extrapolate absolute adsorption uptakes from measured adsorption isotherms. The framework also allows a straightforward calculation of thermodynamic potentials relevant to adsorption such as enthalpy and entropy. Modelling results show that methane adsorption isotherms in shale under different pressures and temperatures are represented by a two dimensional adsorption isotherm surface. The density of the adsorbed methane in shales depends on temperature and pressure, which is always lower than the liquid methane density but higher than the corresponding gaseous methane density. The temperature-dependent and pressure-dependent characteristics of adsorbed methane density leads to the corresponding temperature-dependent and pressure-dependent measured/absolute adsorption isotherms. The maximum adsorption uptake of shales is independent of temperature and pressure. The isosteric enthalpy/entropy of adsorption and enthalpy/entropy of adsorbed methane are found to be temperature- and surface coverage-dependent. These new findings not only clarify some historical misunderstandings of methane adsorption in shales for engineering application, but also provide a novel framework for interpreting methane adsorption behavior in shales and for determining the associated thermodynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamic assessments and inter-relationships between systems involving Al, Am, Ga, Pu, and U

We present a newly developed self-consistent CALPHAD thermodynamic database involving Al, Am, Ga, Pu, and U. A first optimization of the slightly characterized Am-Al and completely unknown Am-Ga phase diagrams is proposed. To this end, phase diagram features as crystal structures, stoichiometric compounds, solubility limits, and melting temperatures have been studied along the U-Al → Pu-Al → Am-Al, and U-Ga → Pu-Ga → Am-Ga series, and the thermodynamic assessments involving Al and Ga alloying are compared. In addition, two distinct optimizations of the Pu-Al phase diagram are proposed to account for the low temperature and Pu-rich region controversy. We included the previously assessed thermodynamics of the other binary systems (Am-Pu, Am-U, Pu-U, and Al-Ga) in the database and is briefly described in the present work. In conclusion, predictions on phase stability of ternary and quaternary systems of interest are reported to check the consistency of the database.

36 MATERIALS SCIENCE↗

Oxygen potentials and phase equilibria in the system Ca–Co–O and thermodynamic properties of Ca{sub 3}Co{sub 2}O{sub 6} and Ca{sub 3}Co{sub 4}O{sub 9.163}

Oxygen potentials established by the equilibrium between three condensed phases, CaO{sub ss}+CoO{sub ss}+Ca{sub 3}Co{sub 2}O{sub 6} and CoO{sub ss}+Ca{sub 3}Co{sub 2}O{sub 6}+Ca{sub 3}Co{sub 3.93+α}O{sub 9.36−δ}, are measured as a function of temperature using solid-state electrochemical cells incorporating yttria-stabilized zirconia as the electrolyte and pure oxygen as the reference electrode. Cation non-stoichiometry and oxygen non-stoichiometry in Ca{sub 3}Co{sub 3.93+α}O{sub 9.36−δ} are determined using different techniques under defined conditions. Decomposition temperatures and thermodynamic properties of Ca{sub 3}Co{sub 2}O{sub 6} and Ca{sub 3}Co{sub 4}O{sub 9.163} are calculated from the results. The standard entropy and enthalpy of formation of Ca{sub 3}Co{sub 2}O{sub 6} at 298.15 K are evaluated. Using thermodynamic data from this study and auxiliary information from the literature, phase diagram for the ternary system Ca–Co–O is computed. Isothermal sections at representative temperatures are displayed to demonstrate the evolution of phase relations with temperature. - Graphical abstract: Isothermal section of the phase diagram of the system Ca–Co–O at 1250 K. - Highlights: • Improved definition of cation and oxygen nonstoichiometry of Ca{sub 3}Co{sub 3.93+α}O{sub 9.36−δ}. • Measurement of Δμ{sub O{sub 2}} associated with two 3-phase fields as a function of temperature. • Use of solid-state electrochemical cells for accurate measurement of Δμ{sub O{sub 2}}. • Decomposition temperatures and thermodynamic properties for ternary oxides. • Characterization of ternary phase diagram of the system Ca–Co–O.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Density functional thermodynamic description of spin, phonon and displacement degrees of freedom in antiferromagnetic-to-paramagnetic phase transition in YNiO 3

This work herein demonstrates a direct density functional description of the finite-temperature thermodynamic properties of solids exhibiting phase transitions through positional and spin symmetry breaking degrees of freedom. A classic example addressed here is the rare-earth (R) nickelates RNiO 3 where the ground state is characterized by crystallographic and magnetic (e.g., antiferromagnetic) long-range order (LRO), whereas the higher temperature paramagnetic phase manifests a range of local spin and positional symmetry breaking motifs with short-range order (SRO). Unlike time-dependent simulations of spin and positional degrees of freedom, in the present work, phases are described via a superposition of static configurations constructed by populating a periodic base lattice supercell allowing for the formation of energy lowing distribution of positional and spin local motifs. The thermal populations of the configurations in such a superposition phase are obtained from the energy-minimized Density Functional Theory (DFT)-calculated partition functions at different temperatures. This approach offers flexible inclusion of different physical contributions to the free energy, such as elastic, electronic and phonon free energies, all obtained from the same underlying DFT total energy calculations of periodic structures. The thermodynamic and magnetic properties of both LRO and SRO crystallographic and spin phases, including antiferromagnetic (AFM) to paramagnetic (PM) Néel phase transition in YNiO3 are studied. Including spin and phonon contributions, we find a DFT-calculated Néel temperature to be 144 K in satisfactory agreement with the experimental value of 145 K; whereas omitting the phonon contribution, one obtains a Néel temperature of 81 K. We present phonon contributions to the DFT-calculated temperature-dependent SRO, heat capacities, and the polymorphous distribution of nonzero local magnetic moments in the PM phase. This approach thus extends to finite temperatures the symmetry-broken DFT description of both the AFM and PM phases, demonstrating that a thermodynamic superposition approach based on symmetry broken configurations evaluated by a mean-field like DFT is sufficient to obtain a consistent description of the thermal physics of the AFM, PM phases and their interconversion in 3d oxides illustrated by YNiO 3 .

36 MATERIALS SCIENCE↗

A pseudo-thermodynamic description of dispersion for nanocomposites

Dispersion in polymer nanocomposites is determined by the kinetics of mixing and chemical affinity. Compounds like reinforcing filler/elastomer blends display some similarity to colloidal solutions in that the filler particles are close to randomly dispersed through processing. It is attractive to apply a pseudo-thermodynamic approach taking advantage of this analogy between the kinetics of mixing for polymer compounds and thermally driven dispersion for colloids. In order to demonstrate this pseudo-thermodynamic approach, two polybutadienes and one polyisoprene were milled with three carbon blacks and two silicas. These samples were examined using small-angle x-ray scattering as a function of filler concentration to determine a pseudo-second order virial coefficient, A 2 , which is used as an indicator for compatibility of the filler and polymer. It is found that A 2 follows the expected behavior with lower values for smaller primary particles indicating that smaller particles are less compatible and more difficult to mix. The measured values of A 2 can be used to specify repulsive interaction potentials for coarse grain DPD simulations of filler/elastomer systems. In addition, new methods to quantify the filler percolation threshold and filler mesh size as a function of filler concentration are obtained. Moreover, the results represent a new approach to understanding and predicting compatibility in polymer nanocomposites based on a pseudo-thermodynamic approach.

36 MATERIALS SCIENCE↗

Alternating crystalline lamellar structures from thermodynamically miscible poly(ε-caprolactone) H/D blends

Thermodynamic interaction and non-isothermal crystallization behavior in a series of protiated (H-) and deuterated (D-) poly(ε-caprolactone) (PCL) blends have been systematically investigated. The blends were thermodynamically miscible in the melt. The Flory–Huggins interaction parameter (χ) between H- and D-PCL segments was estimated. The hydroxyl groups in the PCL chain-ends were found to contribute significantly to the negative χ values. Combined characterization of small-angle X-ray scattering (SAXS) and small-angle neutron scattering (SANS) revealed that upon slow cooling, unique alternating H-rich and D-rich PCL lamellar structures are formed regardless of the blend ratio, which is attributed to the T c difference between H- and D-PCLs. While upon rapid cooling, mixed crystals of H- and D-PCLs are predominantly formed. These results provide insightful information on the melt thermodynamics as well as detailed chain arrangements in lamellar crystals for semi-crystalline H/D polyester blends.

36 MATERIALS SCIENCE↗