SEARCH · Search NASA
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.
A review of the energy level and heat of formation data for the species C2H, C3H, and C4H is given. The procedure for computing thermochemical data from partition function is also reviewed. Working expressions for approximate partition functions, free-energy function, enthalpy function, and heat capacity for a linear polyatomic species are presented.
Flame-spray-pyrolysis (FSP) is a robust and scalable process to synthesize particles at the commodity-scale. FSP has been used to produce the precursor powders which were converted to the layered structure $(R\bar{3m} \text { phase})$ by a postannealing step in making nickel-rich cathode materials (NCMs). Theoretically, the high flame temperature (normally >1500 K) in FSP can provide adequate energy for the phase conversion from rock-salt to layered structures and potentially enables one-step synthesis. However, the high flame temperature is a critical issue to cause lithium loss and structural degradation, preventing the formation of the layered phase. Here, guided by the gaseous nucleation theory, we implemented several FSP processes with different solution recipes. The layered phase concentration in the as-burned products can be increased with the solution enthalpies. By adding a rapid quench step to suppress the lithium loss and phase degradation, the layered phase can be further increased. This work contributes new ideas to innovating process regarding the process efficiency and throughput of manufacturing cathode materials at a large scale.
A generic method to estimate the relative feasibility of formation of high-entropy compounds in a single phase, directly from first principles, is developed. As a first step, the relative formation abilities of 56 multicomponent, AO, oxides were evaluated. These were constructed from five cation combinations chosen from A = {Ca, Co, Cu, Fe, Mg, Mn, Ni, Zn}. Candidates for multicomponent oxides are predicted from descriptors related to the enthalpy and configurational entropy obtained from the mixing enthalpies of two-component oxides. The utility of this approach is evaluated by comparing the predicted combinations with the experimentally realized entropy-stabilized oxide, (MgCoCuNiZn)O. In the second step, Monte Carlo simulations are utilized to investigate the phase composition and local ionic segregation as a function of temperature. This approach allows for the evaluation of potential secondary phases, thereby making realistic predictions of novel multicomponent compounds that can be synthesized.
A computer program is described which (1) calculates thermodynamic functions (heat capacity, enthalpy, entropy, and free energy) for several optional forms of the partition function, (2) fits these functions to empirical equations by means of a least-squares fit, and (3) calculates, as a function of temperture, heats of formation and equilibrium constants. The program provides several methods for calculating ideal gas properties. For monatomic gases, three methods are given which differ in the technique used for truncating the partition function. For diatomic and polyatomic molecules, five methods are given which differ in the corrections to the rigid-rotator harmonic-oscillator approximation. A method for estimating thermodynamic functions for some species is also given.
A phenomenological multiple-relaxation theory of the deactivation rate constant for the nu-2 (1 - 0) bending mode of water vapor is presented which incorporates the role not only of the excited monomer but also of the bound molecular complex, in particular the dimer. The deactivation takes place by means of three parallel processes: (1) collisional deexcitation of the excited monomer, (2) a two-step reaction involving association and spontaneous redissociation of an H2O collision complex, and (3) spontaneous dissociation of the stably bound H2O dimer. Oxygen, but not nitrogen or argon, serves as an effective chaperon for the formation of the activated complex. This observation explains the impurity dependence of the self-deactivation rate constant of water vapor. Analysis of an ultrasonic absorption peak based on the third process yields values for the standard entropy and enthalpy of dissociation of the stably bound H2O dimer.
A set of 58 group additivity values (GAV) for the calculation of the heat of formation is derived from an extensive and accurate database of 192 ab initio heats of formation. The ab initio values are from companion calculations at the CCSD(T)–F12/cc-pVTZ-F12//B2PLYPD3/cc-pVTZ level of theory employing second order connectivity based hierarchy reference reactions with ANL energies for the reference species. This database of 192 species consists of alkanes (RH), alkyl radicals ($\dot{\text {R}}$), alkyl hydroperoxides (RO 2 H), alkyl-peroxy (R$\dot{\text {O}}$ 2 ) radicals, and hydroperoxy-alkyl ($\dot{\text {Q}}$OOH) radicals of all of the isomers of C 2 –C 5 alkane fuels and a select set of C 6 –C 8 isomers. The GAV and ab initio based enthalpies for this dataset show excellent agreement, with an estimated 2σ uncertainty of 0.9 kcal mol –1 . Among the 58 GAV terms, 40 are refinements of previously reported terms, while 18 are newly developed. These new GAV terms are mainly non-next-nearest neighbor interactions (NNI) and β-corrections. The inclusion of these new groups significantly improves the accuracy of the GAV estimates of the heats of formation. Finally, the updated GAVs can be used with increased confidence to estimate the heats of formation of combustion relevant hydrocarbons and oxygenated hydrocarbons and their corresponding radicals, which are important in predicting low-temperature chemistry and are critical in the development of accurate chemical kinetic models.
CALPHAD uncertainty quantification (UQ) is the foundation of materials design with quantified confidence. We report a framework and software packages to enable CALPHAD UQ assessment and calculation using commercial CALPHAD software (Thermo-Calc). This Bayesian inference framework is coupled with a Markov chain Monte Carlo algorithm to establish uncertainty traces with a given thermodynamic database file (TDB) and corresponding experimental data points. This general framework is demonstrated with the Ni–Cr binary system. The algorithm is firstly validated on synthetic data with known ground truth. Then it is applied to real experimental data to generate posterior traces. We develop a file format named TDBX, which provides a single source of truth by combining the original TDB content and the traces for each assessed Gibbs energy parameter. CALPHAD UQ calculations are performed based on the TDBX file, from which uncertainties for phase boundaries, enthalpy curves, and solidification range are collected as examples of basic design parameters. This TDBX file with corresponding scripts are made open-source. Finally, the combination of CALPHAD UQ assessments and calculations connected by TDBX supports uncertainty-assisted modeling, enabling the integrated application of modern design with uncertainty methodologies to computational materials design.
Solar thermochemical hydrogen is one of the few potential routes towards direct fuel production from renewable energy sources, but the thermodynamic boundary conditions for efficient and economic energy conversion are challenging. Success or failure of a given oxide working material depends on the subtle balance between enthalpy and entropy contributions in the redox processes. Developing a mechanistic understanding of the behavior of materials on the basis of atomistic models and first-principles calculations is an important part of advancing the technology. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox-active defects start to interact with each other, thereby impeding the formation of additional defects. This problem is of more general importance to applications that rely on high levels of off-stoichiometry or doping, including, for example, batteries, thermoelectrics, and ceramic fuel cells. To account for such repulsive defect interactions, we introduce a statistical mechanics approach, defining an expression for the free energy of defect interaction based on limited sampling of defect configurations in density functional theory supercell calculations. The parameterization of this energy contribution as a function of defect concentration and temperature allows on-the-fly simulation of thermochemical equilibria. The approach consistently incorporates finite temperature effects by including the leading contributions to the temperature-dependent free energy for the case at hand, i.e., the ideal gas and configurational enthalpies and entropies. We demonstrate the capability and utility of the approach by simulating the water splitting redox processes for Sr 1 − x Ce x Mn O 3 − δ alloys. Published by the American Physical Society 2024
Arguably the most ubiquitous construction material in modern civilization, concrete is enabling the development of megacities around the globe together with increasing living standards in developing nations. However, the exact formation mechanisms of the strength-giving calcium-rich gels remain a topic of debate. Using density functional modeling, we simulate the fundamental solution-based building blocks of cement hydrates (calcium ions and silicate and aluminate monomers) and their propensity to form pair-wise complexes with bonding environments characteristic of those found in calcium-silicate-hydrate, calcium-alumino-silicate-hydrate and sodium-containing calcium-alumino-silicate-hydrate gels, as assessed from Gibbs free energies of chemical reactions. By accurately simulating the high pH pore solution chemistry in Portland cements and related systems, along with discrete solvation of the species, we hypothesize potential early age formation routes of the gels and discuss limitations and future work associated with this approach.
The entropy coefficient of a battery cell is the property that governs the amount of reversible heat that is generated during operation. In this work, we propose an extension of the Multi-Species, Multi-Reaction (MSMR) model to capture the entropy coefficient of a large format lithium-ion battery cell. We utilize the hybridized time-frequency domain analysis (HTFDA) method using a multi-functional calorimeter to probe the entropy coefficient of a large format pouch type lithium-ion battery with a NMC 811 cathode and a graphite anode. The measured entropy coefficient profile of the battery cell is deconvoluted into an entropy coefficient for each active material, which is then estimated using an extension of the MSMR model. Finally, we extend the entropy of a material to individual entropy for each gallery as treated by the model.
An effort is currently underway at NASA Lewis to develop two and three dimensional Navier-Stokes codes, called Proteus, for aerospace propulsion applications. Proteus solves the Reynolds-averaged, unsteady, compressible Navier-Stokes equations in strong conservation law form. Turbulence is modeled using a Baldwin-Lomax based algebraic eddy viscosity model. In addition, options are available to solve thin layer or Euler equations, and to eliminate the energy equation by assuming constant stagnation enthalpy. An extensive series of validation cases have been run, primarily using the two dimensional planar/axisymmetric version of the code. Several flows were computed that have exact solution such as: fully developed channel and pipe flow; Couette flow with and without pressure gradients; unsteady Couette flow formation; flow near a suddenly accelerated flat plate; flow between concentric rotating cylinders; and flow near a rotating disk. The two dimensional version of the Proteus code has been released, and the three dimensional code is scheduled for release in late 1991.
In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).
Multi-principal-element alloys like high-entropy alloys (HEAs) have potential applications in many engineering fields due to their unique mechanical/functional properties. While HEAs are generally considered random solid solutions, recent studies revealed that they are prone to short-range-ordering (SRO) due to the complex multi-pair-wise interactions among the constituent elements. Meanwhile, SROs' evolution can sometimes be deleterious, and it is necessary to have control over their evolution. Examining the AlCoCrFe-Zr model alloy, long-range ordering occurs following the expectation of enthalpic predictions. Advanced characterization techniques—transmission electron microscopy, high-energy synchrotron X-ray diffraction/pair distribution function, and atom probe tomography, reveal that SRO is suppressed in as-milled and GB-decorated NC-(AlCoCrFe)100-xZrx (x = 0–1.5 atomic %). Warren-Cowley coefficient calculations are further used to validate the suppression of SRO. Besides the low segregation enthalpies of Cr, Fe, and Zr, and the high-mixing enthalpy of Cr and Fe, the short diffusion path to GBs due to high-GB density in the NC-HEAs and the higher energy state of the GBs than the matrix promotes GB-segregation that further alters the matrix chemistry and consequently disfavors SRO formation within the matrix. Despite the GB-segregation of Cr, Fe, and Zr, the matrices and GBs remain in a random solid solution.
We present a new mesoscale model for ionic liquids based on a low Mach number fluctuating hydrodynamics formulation for multicomponent charged species. The low Mach number approach eliminates sound waves from the fully compressible equations leading to a computationally efficient incompressible formulation. The model uses a Gibbs free energy functional that includes enthalpy of mixing, interfacial energy, and electrostatic contributions. These lead to a new fourth-order term in the mass equations and a reversible stress in the momentum equations. We calibrate our model using parameters for [DMPI+][F6P-], an extensively-studied room temperature ionic liquid (RTIL), and numerically demonstrate the formation of mesoscopic structuring at equilibrium in two and three dimensions. In simulations with electrode boundaries the measured double layer capacitance decreases with voltage, in agreement with theoretical predictions and experimental measurements for RTILs. Finally, we present a shear electroosmosis example to demonstrate that the methodology can be used to model electrokinetic flows.
Exothermic reactions during the annealing of laboratory synthesized amorphous magnesium-bearing silicate particles used as grain analogs of cosmic dust were detected by differential scanning calorimetry (DSC) in air. With infrared spectroscopy and transmission electron microscopy, we show that cosmic dust could possibly undergo fusion to larger particles, with oxidation of magnesium silicide and crystallization of forsterite as exothermic reactions in the early solar system. The reactions begin at approximately 425, approximately 625, and approximately 1000 K, respectively, and the reaction energies (enthalpies) are at least 727, 4151, and 160.22 J per gram, respectively. During the crystallization of forsterite particles, the spectral evolution of the 10 micrometer feature from amorphous to crystalline was observed to begin at lower temperature than the crystallization temperature of 1003 K. During spectral evolution at lower temperature, nucleation and/or the formation of nanocrystallites of forsterite at the surface of the grain analogs was observed.
A group of multi-component oxides based on BaZrO 3 have been prepared using a solid-state reaction method and examined in terms of their water uptake and thermodynamics of formation. Depending on the type and amount of acceptor substitution, the synthesized compounds exhibit various proton defect concentrations, reaching up to 0.2 mol/mol for a compound containing 10 different elements in the B-sublattice, where 50% of them are acceptors. Further, for the most promising materials, van’t Hoff plots were created and the enthalpies and entropies of hydration were calculated. At higher temperatures, these parameters do not differ from the values for the reference yttrium doped barium zirconate. However, at lower temperatures they are more negative, indicating a more exothermic process of proton incorporation.
Equilibrium thermodynamic properties for pure ammonia were generated for a range of temperature from 500 to 50,000 K and pressure from 0.01 to 40 MN/sq m and are presented in tabulated and graphical form. Properties include pressure, temperature, density, enthalpy, speed of sound, entropy, molecular-weight ratio, specific heat at constant pressure, specific heat at constant volume, isentropic exponent, and species mole fractions. These properties were calculated by the method which is based on minimization of the Gibbs free energy. The data presented herein are for an 18-species ammonia model. Heats of formation and spectroscopic constants used as input data are presented. Comparison of several thermodynamic properties calculated with the present program and a second computer code is performed for a range of pressure and for temperatures up to 30,000 K.
Using two-dimensional particle-in-cell (PIC) together with magnetohydrodynamic (MHD) Q1 simulations of magnetotail dynamics, we investigate the evolution toward onset of reconnection and the subsequent energy transfer and conversion. In either case, reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. The boundary deformation leads to the formation of thin embedded current sheets, which are bifurcated in the near tail, converging to a single sheet farther out in the MHD simulations. The thin current sheets in the PIC simulation are carried by electrons and are associated with a strong perpendicular electrostatic field, which may provide a connection to parallel potentials and auroral arcs and an ionospheric signal even prior to the onset of reconnection. The PIC simulation very well satisfies integral entropy conservation (intrinsic to ideal MHD) during this phase, supporting ideal ballooning stability. Eventually, the current intensification leads to the onset of reconnection, the formation and ejection of a plasmoid, and a collapse of the inner tail. The earthward flow shows the characteristics of a dipolarization front: enhancement of Bz, associated with a thin vertical electron current sheet in the PIC simulation. Both MHD and PIC simulations show a dominance of energy conversion from incoming Poynting flux to outgoing enthalpy flux, resulting in heating of the inner tail. Localized Joule dissipation plays only a minor role.