Search NASA⌕ Search

SEARCH · Search NASA

Results for “total enthalpy”

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 19 records

Assessment of Enthalpy-Based Conditional Moment Closure Models in Predicting Ignition of Lean and Stoichiometric PRF-Air Mixtures with Temperature Inhomogeneity

Accurate prediction of the ignition delay in turbulent flows is important in developing advanced combustion technology for internal combustion engines. Enthalpy-based conditional moment closure (CMC) models are assessed with 2-D direct numerical simulation (DNS) data for ignition of homogeneous PRF-air mixtures with temperature inhomogeneity in homogeneous isotropic turbulence. The sensible-enthalpy-based CMC model is compared with the total-enthalpy-based formulation, with the focus being placed on assessing the capability of the CMC models in predicting ignition processes under different engine conditions. A priori analysis shows that the total-enthalpy-based CMC gives good predictions of the heat release rates in the lean mixtures when the temperature inhomogeneity level is low, but leads to significant overprediction under the stoichiometric condition regardless of the thermal stratification levels. On the other hand, the sensible enthalpy formulation performs well for all conditions considered in this study due to the suppressed conditional temperature fluctuations. A mapping method to model conditional scalar dissipation rates and probability density function of sensible enthalpy is also assessed.

Wang, Wei↗

Heat of Mixing During Fast Charge/Discharge of a Li-Ion Cell: A Study on NMC523 Cathode

Predicting temperature rise accurately during fast charge/discharge of a Li-ion cell is essential to avoid thermal runaway and extend battery life. While modeling the temperature rise, it is necessary to model the heat generation correctly. Heat of mixing, one of the four main sources of heat generation in a Li-ion battery, has often been considered insignificant and therefore excluded from modeling. When included, it is modeled using the expression from a Taylor expansion approximation. In this work, we have shown the conditions when including the heat of mixing becomes important and quantified the error associated with using the Taylor expansion, especially under high charge/discharge conditions. Consequently, we carry out the calculation of the rate of enthalpy change and the heat generation rate from the most fundamental equation for the rate of the total enthalpy change without simplifying assumptions or approximations. The heat generation rate calculated doing so naturally includes irreversible, reversible and mixing heat. We then exclusively separate out the heat of mixing by subtracting the rate of enthalpy change by reaction from the rate of the total enthalpy change. Results show that the contribution of heat of mixing in the total heat generated increases with the charge/discharge rate and is as large as 23% for a 6 C discharge. This result suggests that while modeling heat generation for fast charge/discharge, it is necessary to include the heat of mixing and avoid calculating it using the Taylor expansion approximation.

25 ENERGY STORAGE↗

Sensible-Enthalpy-Based Conditional Moment Closure Model for Homogeneous Charge Compression Ignition with Temperature Inhomogeneity

Abstract This paper presents a sensible-enthalpy-based conditional moment closure (CMC) model for homogeneous charge compression ignition (HCCI) with temperature inhomogeneity. The focus of study is placed on the effects of a conditioning variable on the prediction of heat release rates and developing the mixing models for sensible enthalpy, a reacting scalar. The sensible-enthalpy-based CMC model is validated a priori with 2-D direct numerical simulation (DNS) data for HCCI, and compared with the total-enthalpy-based formulation. Sensible enthalpy is shown to be strongly correlated with temperature, which makes sensible enthalpy well suited for a conditioning variable in CMC for HCCI. Results show that sensible-enthalpy-based CMC performs well for the prediction of heat release rates for all temperature stratification levels investigated in this study, while total-enthalpy-based CMC leads to the overprediction of heat release rates when the temperature stratification level is high. While conditional fluctuations of species mass fractions are substantial, sensible-enthalpy-based CMC is found to accurately estimate heat release rates during ignition, as well as reaction rates of key species, due to substantially suppressed conditional temperature fluctuations. The mapping closure models for the probability density function (PDF) and conditional dissipation rates of sensible enthalpy are proposed. The proposed models are found to capture the characteristics of reactive scalar mixing.

42 ENGINEERING↗

Investigating the Role of Energy Density in Thermal Runaway of Lithium-Ion Batteries with Accelerating Rate Calorimetry

This work uses accelerating rate calorimetry to evaluate the impact of cell chemistry, state of charge, cell capacity, and ultimately cell energy density on the total energy release and peak heating rates observed during thermal runaway of Li-ion batteries. While the traditional focus has been using calorimetry to compare different chemistries in cells of similar sizes, this work seeks to better understand how applicable small cell data is to understand the thermal runaway behavior of large cells as well as determine if thermal runaway behaviors can be more generally tied to aspects of lithium-ion cells such as total stored energy and specific energy. We have found a strong linear correlation between the total enthalpy of the thermal runaway process and the stored energy of the cell, apparently independent of cell size and state of charge. We have also shown that peak heating rates and peak temperatures reached during thermal runaway events are more closely tied to specific energy, increasing exponentially in the case of peak heating rates.

25 ENERGY STORAGE↗

An Automated Thermochemistry Protocol Based on Explicitly Correlated Coupled-Cluster Theory: The Methyl and Ethyl Peroxy Families

An automated computational thermochemistry protocol based on explicitly correlated coupled-cluster theory was determined to produce highly accurate enthalpies of formation and atomization energies for small- to medium-sized molecular species (3–12 atoms). Each possible source of error was carefully examined, and the sizes of contributions to the total atomization enthalpies were used to generate uncertainty estimates. The protocol was first used to generate total atomization enthalpies for a family of four molecular species exhibiting a variety of charges, multiplicities, and electronic ground states. The new protocol was demonstrated to be in good agreement with the Active Thermochemical Tables database for the four species: the methyl peroxy radical, methoxyoxoniumylidene (methyl peroxy cation), methyl peroxy anion, and methyl hydroperoxide. Updating the Active Thermochemical Tables to include those results yielded significantly improved accuracy for the formation enthalpies of those species. The derived protocol was then used to predict formation enthalpies for the larger ethyl peroxy family of species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mass, enthalpy, and chemical‐derived emission flows in mineral processing

Abstract The production of materials from mineral resources is a significant contributor to anthropogenic CO 2 emissions. This contribution is driven primarily by chemical CO 2 emissions from the conversion of mineral resources and emissions tied to energy demands for material processing. In this work, we synthesize the thermodynamically required enthalpy and chemically derived emissions of mineral processing and consumption in the United States. We quantify mass, enthalpy, and emissions flows for minerals described by the US Geological Survey, with 882 mass flows and 155 chemical reactions analyzed. In total, 503 PJ of enthalpy is thermodynamically required for 398 Mt of chemically converted material consumption in the United States, resulting in 129 Mt of chemically derived CO 2 emissions. Additionally, 249 PJ of fuel resources such as coke are stoichiometrically required for the chemical conversion of minerals. These enthalpy requirements and CO 2 emissions are primarily from high‐mass consumption materials such as cement, carbon steel, fertilizer, and aluminum. Cumulatively, the dataset synthesized in this work provides a complete view of the chemical requirements of mineral processing and can aid in guiding decarbonization or sustainable growth in critical minerals sectors, including construction materials and materials for energy storage or generation.

Kane, Seth↗

Determination of the regeneration energy of direct air capture solvents/sorbents using calorimetric methods

Regeneration energy and recyclability are important parameters that dictate the viability of a CO 2 direct air capture (DAC) process. In this study, differential scanning calorimetry (DSC) is used in conjunction with thermogravimetric analysis (TGA) and Fourier Transform Infrared Spectroscopy (FTIR) to evaluate the regeneration energy required for CO 2 -loaded 3 M aqueous potassium sarcosinate (K-SAR) solvent and crystalline methyl-glyoxal-bisiminoguanidine (MGBIG) sorbent. Based on calorimetric measurements, for aqueous K-SAR, the sensible heat amounts to 1.5 GJ/tCO 2 and the enthalpy of desorption is estimated at 3.68 GJ/tCO 2 . The total energy of regeneration will also include the enthalpy of solvent vaporization, which will depend on the thermodynamic conditions and the level of regeneration. For MGBIG, the total regeneration energy required was measured at approximately 7.0 GJ/tCO 2 . Further, FTIR measurements were used to observe CO 2 and H 2 O release during the thermal regeneration process. For K-SAR, CO 2 release is seen to take place at temperatures greater than 80 °C and, for MGBIG, the CO 2 release takes place at temperatures greater than 100 °C. The moderate temperatures required for K-SAR and MGBIG regenerations compared to other DAC solvents/sorbents suggest that low-cost sources of heat, such as geothermal heat, could potentially be used for regeneration to drive down the cost of direct air capture.

42 ENGINEERING↗

Data-driven approach to parameterize SCAN + U for an accurate description of 3 d transition metal oxide thermochemistry

Semilocal density-functional theory (DFT) methods exhibit significant errors for the phase diagrams of transition-metal oxides that are caused by an incorrect description of molecular oxygen and the large self-interaction error in materials with strongly localized electronic orbitals. Empirical and semiempirical corrections based on the DFT + U method can reduce these errors, but the parameterization and validation of the correction terms remains an on-going challenge. Here we develop a systematic methodology to determine the parameters and to statistically assess the results by considering interlinked thermochemical data across a set of transition metal compounds. We consider three interconnected levels of correction terms: (1) a constant oxygen binding correction, (2) Hubbard-U correction, and (3) DFT/DFT + U compatibility correction. The parameterization is expressed as a unified optimization problem. We demonstrate this approach for 3d transition metal oxides, considering a target set of binary and ternary oxides. With a total of 37 measured formation enthalpies taken from the literature, the dataset is augmented by the reaction energies of 1710 unique reactions that were derived from the formation energies by systematic enumeration. To ensure a balanced dataset across the available data, the reactions were grouped by their similarity using clustering and suitably weighted. The parameterization is validated using leave-one-out cross validation (CV), a standard technique for the validation of statistical models. We apply the methodology to the strongly constrained and appropriately normed (SCAN) density functional. Based on the CV score, the error of binary (ternary) oxide formation energies is reduced by 40% (75%) to 0.10 (0.03) eV/atom. A simplified correction scheme that does not involve SCAN/SCAN + U compatibility terms still achieves an error reduction of 30% (25%). The method and tools demonstrated here can be applied to other classes of materials or to parameterize the corrections to optimize DFT + U performance for other target physical properties.

36 MATERIALS SCIENCE↗

Partial molar properties from molecular simulation using multiple linear regression

Partial molar volumes, energies, and enthalpies can be computed from NpT-Gibbs ensemble simulations through a post-processing procedure that leverages fluctuations in composition, total volume, and total energy of a simulation box. By recording the instantaneous box volumes V and instantaneous number of molecules N i of each of n species for M frames, a large M ×n matrix N is constructed, as well as the M × 1 vector V. The 1 × n vector of partial molar volumes V¯ may then be solved using N • V¯=V. A similar construction permits calculation of partial molar energies using M instantaneous measurements of the total energy of the simulation box, and N • U¯=U. Partial molar enthalpies may be derived from U¯, V¯, and pressure p. These properties may be used to calculate enthalpy and entropy of transfer (absorption, extraction, and adsorption) for species in complex mixtures. The method is demonstrated on three systems in the NpT-Gibbs ensemble: a highly compressible natural gas condensate of methane, n-butane, and n-decane, the liquid-phase adsorption of 1,5-pentanediol and ethanol onto the MFI zeolite, and a relatively incompressible mixture of ethanol, n-dodecane, and water at liquid-liquid equilibrium. Property predictions are compared to those from numerical differentiation of simulation data sequentially changing the composition and from equations of state. Here, the method can also be extended to reaction equilibria in a closed system and is applied to a reactive first-principles Monte Carlo simulation of compressed nitrogen/oxygen.

74 ATOMIC AND MOLECULAR PHYSICS↗

Thermodynamics of MgO Atomic Layer Deposition Surface Reactions

The realities of atomic layer deposition (ALD) surface reactions often deviate from the simple ligand exchange frequently used to illustrate the technique. A detailed understanding of these reactions is necessary to develop greater surface synthetic control including chemical selectivity for patterning or to target defects. Here, the thermodynamics of surface reactions relevant to MgO ALD were investigated using pyroelectric calorimetry to measure the time-resolved heat generation. These reactions show exothermic heat generation of 0.12 mJ/cm 2 and 0.15 mJ/cm 2 for alternating Mg(CpEt) 2 and H 2 O reactions, respectively. The total reaction heat closely matches the standard reaction enthalpy for bulk MgO, supplemented with first-principles molecular calculations. First-principles models further reveal that while simple ligand exchange is favorable during surface reactions, the required increase in Mg-coordination number from two in the precursor to six in bulk MgO requires additional coverage-dependent surface reactions, which depend on the availability and stability of proximal surface hydroxyls.

36 MATERIALS SCIENCE↗

Large-Scale Molecular Dynamics Simulation of the Dehydration of a Suspension of Smectite Clay Nanoparticles

Fine-grained sediments and sedimentary rocks play important roles in a variety of modern energy technologies from petroleum geology to geological carbon sequestration and radioactive waste management. However, despite their utility and ubiquity, many of their properties remain poorly understood. In particular, the ability to predict the permeability and mechanics of these media remains a persistent fundamental challenge in the geosciences. In the present work, we show how large-scale classical molecular dynamics (MD) simulations can help interpret the properties of fine-grained sedimentary material. All-atom simulations containing 30 discrete clay particles are utilized to understand the evolution of a clay nanoparticle suspension during its progressive dehydration. Microstructural (pore size distribution, tortuosity, anisotropy), thermodynamic (enthalpy and free energy of hydration, anion exclusion), mechanical (total suction), and transport properties (diffusion coefficient tensors of water and sodium) are calculated and compared to the experiment. Overall, our results provide new insight into the coupled chemistry, mechanics, and transport properties of disordered nanoparticle assemblages and shed light upon the important role of water films in controlling these properties.

02 PETROLEUM↗

Conductivity, structure, and thermodynamics of Y 2 Ti 2 O 7 –Y 3 NbO 7 solid solutions

The defect fluorite yttrium niobate Y 3 NbO 7 and pyrochlore yttrium titanate Y 2 Ti 2 O 7 solid solutions have been synthesized via a solid state synthesis route. The resulting stoichiometry of the oxides is Y 2+x Ti 2-2x Nb x O 7 , where x = 0 to x = 1. All of the samples were single-phase; however, for those with a predominant fluorite phase, a small amount of additional pyrochlore phase was detected. The volume of the solid solution unit cells linearly increases with increase in yttrium niobate content. The water uptake increases with (x) and the protonic defect concentration reaches almost 4.5 × 10 -3 mol mol -1 at 300 °C. The calculated enthalpy of formation from oxides suggests strong stability for all of the compositions, with the values of enthalpy ranging from -84.6 to -114.3 kJ mol -1 . The total conductivity does not have a visible dependence on Y 3 NbO 7 content. For each compound, the total conductivity is higher in wet air. Interestingly, for samples where x < 0.5, the ratio of conductivity in hydrogen to air increases with increasing temperature, while for x > 0.5, the trend is the opposite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transported PDF Modeling of Compressible Turbulent Reactive Flows by using the Eulerian Monte Carlo Fields Method

Although the transported probability density function (PDF) method has been developed for decades, its application has been mainly focused on the low-Mach number flow problems. This work extends the transported PDF method to compressible flow problems. The Eulerian Monte Carlo fields (EMCF) solution method is employed to solve the transported PDF equation for compressible flow problems. A pseudo stagnation enthalpy is introduced and its stochastic partial differential equation is derived to ensure total energy conservation numerically. A new mixing model called interaction by partial exchange with mean (IPEM) is introduced to expand the available choices of mixing models for the EMCF method. The consistency of the EMCF method is examined for solving the transported PDF equation. Numerical implementation details are discussed, such as the density coupling between the compressible flow solver and the EMCF solver, discretization schemes for the mixing terms and the stochastic terms. The implemented compressible flow solver coupled with the EMCF solver is verified and validated in a series of test cases with increasing level of complexity, ranging from a statistically one-dimensional turbulent mixing layer to a self-excited resonance model rocket combustor. It is observed that in general with the increase of compressibility, there is an increase in the sensitivity of the modeling results to the different models and algorithms. This makes it necessary to develop a thorough understanding of the model sensitivity in order to develop a robust and accurate simulation solver for highly compressible turbulent reactive flows. The thermo-acoustic instability inside the model rocket combustor case is captured reasonably, which demonstrates the overall capability of the developed compressible turbulent combustion solver based on the transported PDF method.

42 ENGINEERING↗

Characterization of the arc in crossflow using a two-temperature nonequilibrium plasma flow model

Diverse industrial applications such as circuit breakers and wire arc spraying involve the interaction between an electric arc and a stream of gas impinging perpendicular to it, a configuration commonly referred to as the arc in crossflow. The arc in crossflow is simulated using a three-dimensional time-dependent two-temperature (heavy-species and electrons) plasma flow model to better capture plasma-gas interactions and deviations from Local Thermodynamic Equilibrium (LTE). The coupled fluid-electromagnetic flow model is solved in a monolithic manner using Variational Multiscale Finite Element Method. Simulation findings are validated with experimental findings and contrasted against results obtained with a LTE model. Results from the two-temperature model corroborate experimental observations while providing quantification of the deviation between heavy-species and electron temperatures. The model is used to characterize the arc in crossflow as a function of the Reynolds and Enthalpy dimensionless numbers, which encapsulate the inter-dependence among the main parameters total current, inflow velocity, and inter-electrode spacing. Furthermore, the characterization revealed the behavior of arc shape, voltage drop, arc power, the degree of nonequilibrium, as well as the characteristic plasma front thickness, with varying controlling parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First principle study of structural, electronic, optical, and transport properties of ternary compounds NaGaX{sub 2} (X = S, Se, and Te) in tetragonal chalcopyrite phase

Eemploying the all electron full potential linearized augmented plane wave plus local orbital (FP-LAPW + lo) method based on density functional theory (DFT) structural, electronic, optical, and transport properties of compounds NaGaX{sub 2} (X = S, Se, and Te) in tetragonal phase have been investigated at ground state. The exchange-correlation effects have been included using the Wu Cohen (WC) generalized gradient approximation (GGA) and modified Becke–Johnson (mBJ) exchange potentials. Enthalpy of formation values ensure that crystal structures are stable in tetragonal chalcopyrite phase. Total and partial density of states, and band structures have been plotted to calculate electronic properties which exhibits direct band gap for all compounds. The electronic band gap of NaGaS{sub 2} (2.1 eV) is found to be higher than the electronic band gaps of NaGaSe{sub 2} (1.7 eV) and NaGaTe{sub 2} (1.2 eV) compounds. Using real and imaginary parts of the dielectric function, optical properties such as refractive index (n), extinction coefficient (k), reflectivity (R), absorption coefficient (α), and energy loss spectrum (L) were determined. Transport properties like Seebeck coefficient, thermal and electrical conductivities, figure of merit, and power factors have also been reported by combining the results from DFT and Boltzmann transport theory. The high absorption peaks in visible and ultraviolet regions of the energy spectrum and direct band gaps of these compounds ensure that three compounds are promising materials for photovoltaic applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Phase Transformation Contributions to Heat Capacity and Impact on Thermal Diffusivity, Thermal Conductivity, and Thermoelectric Performance

Abstract The accurate characterization of thermal conductivity κ, particularly at high temperature, is of paramount importance to many materials, thermoelectrics in particular. The ease and access of thermal diffusivity D measurements allows for the calculation of κ when the volumetric heat capacity, ρc p , of the material is known. However, in the relation κ = ρc p D , there is some confusion as to what value of c p should be used in materials undergoing phase transformations. Herein, it is demonstrated that the Dulong–Petit estimate of c p at high temperature is not appropriate for materials having phase transformations with kinetic timescales relevant to thermal transport. In these materials, there is an additional capacity to store heat in the material through the enthalpy of transformation Δ H . This can be described using a generalized model for the total heat capacity for a material where φ is an order parameter that describes how much latent heat responds “instantly” to temperature changes. Here, C pφ is the intrinsic heat capacity (e.g., approximately the Dulong–Petit heat capacity at high temperature). It is shown experimentally in Zn 4 Sb 3 that the decrease in D through the phase transition at 250 K is fully accounted for by the increase in c p , while κ changes smoothly through the phase transition. Consequently, reports of κ dropping near phase transitions in widely studied materials such as PbTe and SnSe have likely overlooked the effects of excess heat capacity and overestimated the thermoelectric efficiency, zT .

Agne, Matthias T.↗

Thermodynamics of Liquid Uranium from Atomistic and Ab Initio Modeling

We present thermodynamic properties for liquid uranium obtained from classical molecular dynamics (MD) simulations and the first-principles theory. The coexisting phases method incorporated within MD modeling defines the melting temperature of uranium in good agreement with the experiment. The calculated melting enthalpy is in agreement with the experimental range. Classical MD simulations show that ionic contribution to the total specific heat of uranium does not depend on temperature. The density of states at the Fermi level, which is a crucial parameter in the determination of the electronic contribution to the total specific heat of liquid uranium, is calculated by ab initio all electron density functional theory (DFT) formalism applied to the atomic configurations generated by classical MD. The calculated specific heat of liquid uranium is compared with the previously calculated specific heat of solid γ-uranium at high temperatures. The liquid uranium cannot be supercooled below T sc ≈ 800 K or approximately about 645 K below the calculated melting point, although, the self-diffusion coefficient approaches zero at T D ≈ 700 K. Uranium metal can be supercooled about 1.5 times more than it can be overheated. The features of the temperature hysteresis are discussed.

36 MATERIALS SCIENCE↗

Particle Image Velocimetry and Total Temperature Characterization of Arc-Heated Mach 4.5 Free Jet

Particle image velocimetry (PIV) measurements were conducted in the Arc-Heated Combustion Tunnel (ACT-II) on a high-enthalpy Mach 4.5 nozzle flow. The primary motivation for the current work is to establish the feasibility of PIV as a high-quality flow diagnostic technique in hypersonic arc-heated flows and produce a high-resolution hypersonic nozzle flow characterization, including nozzle free jet shear layer velocity profiles. Due to the uncertain nature of the gas composition in the flowfield of interest, PIV measurements of a supersonic wedge flow were used to couple the flow velocities to the oblique shock relations. This enabled a PIV-derived total temperature calculation that does not require a priori assumptions of the exact gas composition. A PIV expansion module was developed for this study, and proper seed sizing and material were selected for the present flow. An analysis of the seed particle response from the PIV measurements provided an initial validation of the seed performance while revealing potential issues with consistent seed sizing. These measurements and derived conditions represent a means of acquiring high-resolution velocimetry and temperature datasets that consider uncertainties in the experimental generation of hypersonic high-enthalpy flows.

Engineering↗