Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermochemistry”

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 73 records · Page 4

Thermochemistry of Silicates

The thermodynamic properties of vapor and condensed phases of silicates are crucial in many fields of science. These quantities address fundamental questions on the formation, stability, transformation, and physical properties of silicate minerals and silicate coating compositions. Here the thermodynamic activities of silica and other species in solid solution have been measured by the analysis of the corresponding high temperature vapors using Knudsen Effusion Mass Spectrometry (KEMS). In first set of experiments KEMS has been used to examine the volatility sequence of species (Fe, SiO, Mg, O2 and O) present in the vapor phase during heating of fosterite-rich olivine (Fo93Fa7) up to 2400 C and to measure the Fe, SiO and Mg activities in its solid solution. The data of fosterite-rich olivine are essential for thermochemical equilibrium models to predict the atmospheric and surface composition of hot, rocky exoplanets (Lava Planets). In the second set of experiments the measured thermodynamic activities of the silica in Y2O3-SiO2 and Yb2O3-SiO2 systems are used to assess their reactivity and degradation recession as environmental barrier coatings (EBCs) in combustion environments (e.g. non-moveable parts of gas turbine engine).

Silicates↗

Thermochemistry of CaO-MgO-Al2O3-SiO2 (CMAS) and Advanced Thermal and Environmental Barrier Coating Systems

CaO-MgO-Al2O3-SiO2 (CMAS) oxides are constituents in a broad number of materials and minerals which have recently inferred to discussions in materials science, planetary science, geochemistry and cosmochemistry communities. In materials science, there is increasing interest in the degradation studies of thermal (TBC) and environmental (EBC) barrier coatings of gas turbines by molten CMAS. These coatings have been explored to be applied on silicon-based ceramics and composites which are lighter and more temperature capable hot-section materials of gas turbines than the current Ni-based superalloys. The degradation of the coatings occurs when CMAS minerals carried by the intake air into gas turbines, e.g. in aircraft engines, reacts at high temperatures (1000C) with the coating materials. This causes premature failure of the static and rotating components of the turbine engines. We discuss some preliminary results of the reactions between CMAS and Rare-Earth (RE Y, Yb and Gd) oxide stabilized ZrO2 systems, and stability of the resulting oxides and silicates.

CaO-MgO-Al2O3-SiO2 (CMAS)↗

Thermochemistry of CaO-MgO-Al2O3-SiO2 (CMAS) and Advanced Thermal and Environmental Barrier Coating Systems

CaO-MgO-Al2O3-SiO2 (CMAS) oxides are constituents in a broad number of materials and minerals which have recently inferred to discussions in materials science, planetary science, geochemistry and cosmochemistry communities. In materials science, there is increasing interest in the degradation studies of thermal (TBC) and environmental (EBC) barrier coatings of gas turbines by molten CMAS. CMAS minerals usually are carried by the intake air into gas turbines, e.g. in aircraft engines, and their deposits react at high temperatures (1000C) with the coating materials. This causes degradation and accelerated failure of the static and rotating components of the turbine engines. We discuss some preliminary results of the reactions between CMAS and Rare-Earth (RE Y, Yb, Dy, Gd, Nd and Sm) oxide stabilized ZrO2 or HfO2 systems, and the stability of the resulting oxides and silicates. Plasma sprayed hollow tube samples ( 2.2 mm and 26 mm height) were half filled with CMAS powder, wrapped and sealed with platinum foil, and heat treated at 1310 C for 5h. Samples were characterized by differential scanning calorimetry, X-ray diffraction and cross section electron microscopy analysis.

degradation↗

Thermochemistry of Calcium-Magnesium-Aluminum-Silicate (CMAS) and Components of Advanced Thermal and Environmental Barrier Coating Systems

There is increasing interest in the degradation mechanism studies of thermal and environmental barrier coatings (TEBCs) of gas turbines by molten CaO-MgO-Al(exp. 2)O(exp. 3)-SiO(exp. 2) CMAS). CMAS minerals are usually referred as silicon-containing sand dust and volcano ash materials that are carried by the intake air into gas turbines, e.g. in aircraft engines, and their deposits often react at high temperatures (greater than 1200 degrees C) with the engine turbine coating systems and components. The high temperature reactions causes degradation and accelerated failure of the static and rotating components of the turbine engines. We discuss some results of the reactions between the CMAS and Rare-Earth (RE = Y, Yb, Dy, Gd, Nd and Sm) - oxide stabilized ZrO(exp. 2) or HfO(exp. 2) systems, and the stability of the resulting oxides and silicates. Plasma sprayed hollow tube samples (outside diameter = 4.7 mm, wall thickness = 0.76 mm and = 26 mm height) were half filled with CMAS powder, wrapped and sealed with platinum foil, and heat treated at 1310 degrees C for 5h. Samples were characterized by differential scanning calorimetry (DSC), X-ray diffraction, and cross-section electron microscopy analysis and energy dispersive X-ray spectroscopy. It was found that CMAS penetrated the samples at the grain boundaries and dissolved the TEBC materials to form silicate phases containing the rare-earth elements. Furthermore, it was found that apatite crystalline phases were formed in the samples with total rare-earth content higher than 12 mol% in the reaction zone for the ZrO(exp. 2) system. In general, samples with the nominal compositions (30YSZ), HfO(exp. 2)-7Dy(exp. 2)O(exp. 2) and ZrO(exp. 2)-9.5Y(exp. 2)O(exp. 3)-2.25Gd(exp. 2)O(exp. 3)-2.25Yb(exp. 2)O(exp. 3) exhibited lower reactivity or more resistance to CMAS than the other coating compositions of this work.

degradation↗

Characterization of CO Thermochemistry in Incident Shockwaves

Incident shock waves in pure CO have been characterized in the Electric Arc Shock Tube facility at NASA Ames Research Center. Spectrally and spatially resolved emission spectra characterize radiative signatures from CO in the VUV (Vacuum Ultra-Violet) and mid-infrared and atomic carbon and C2 in the visible. CO absorption of a single vibrational line is also measured with a tunable diode laser. The experimental data analyzed here are at a pressure of 0.25 Torr in the driven section and span a shock velocity range from 3.4-9.5 km/s. The emission and absorption signals are analyzed to extract temperature relaxation behind the shock which is used to derive the rate of CO dissociation. The emission spectra are compared to results using different kinetic parameters for CO dissociation and C2 dissociation and exchange. Different rates from the literature are found to match the data from 3.4-6.6 km/s and 6.6-9.5 km/s. Areas for improvement in CO and C2 radiation modeling are suggested on the basis of the analysis.

Entry↗

Characterization of CO Thermochemistry in Incident Shockwaves

Incident shock waves in pure CO have been characterized in the Electric Arc ShockTube facility at NASA Ames Research Center. Spectrally and spatially resolved emissionspectra characterize radiative signatures from CO in the VUV and mid-infrared and atomiccarbon and C2 in the visible. CO absorption of a single vibrational line is also measuredwith a tunable diode laser. The experimental data analyzed here are at a pressure of 0.25Torr in the driven section and span a shock velocity range from 3.4-9.5 km/s. The emissionand absorption signals are analyzed to extract temperature relaxation behind the shockwhich is used to derive the rate of CO dissociation. The emission spectra are compared toresults using different kinetic parameters for CO dissociation and C2 dissociation andexchange. Different rates from the literature are found to match the data from 3.4-6.6 km/sand 6.6-9.5 km/s. Areas for improvement in CO and C2 radiation modeling are suggested onthe basis of the analysis.

CO↗

Thermochemistry of Strong Air Plasmas for Hypervelocity Earth Entry of Asteroids

The specific heat, enthalpy and free energy of nitrogen and oxygen are calculated for pressures between 1 and 100 bars and temperatures between 10,000 and 100,000 degrees Kelvin. The data are presented in a form suitable for input into computational fluid dynamics computer codes. This range of temperature and pressure are typical of the conditions expected in the bow shock layer of a meteor or small asteroid entering Earth's atmosphere. The extent of double and triple ionization of nitrogen and oxygen atoms under these conditions will be estimated.

Ionization Energy Depression↗

An Assessment of Turbulence Models for MSR Sample Retrieval Lander Aerothermal Design

Several cases from past entries in the T9 and T5 hypersonic test facilities are examined using various turbulence models and shock-layer thermochemistries with the objectives of verifying/validating turbulence models and determining appropriate multiplicative factors on predicted turbulent heating for use in aerothermal design of the heatshield of the Sample Retrieval Lander (SRL) of the Mars Sample Return (MSR) Flagship mission. The zero-equation turbulence models of Baldwin-Lomax and Cebeci-Smith are shown to replicate measurements to within experimental uncertainty for the T9 cases. For T5 cases, the zero-equation model of Baldwin-Lomax can also replicate measurements provided the shock-layer thermochemistry is modified. The modified thermochemistry resolves a dichotomy from past CFD that required different wall catalycity assumptions in the laminar and turbulent flow regimes. The present work shows that such assumptions are not necessary. Finally, an initial estimate for a multiplicative factor on turbulent heating is estimated for the SST turbulence model – the chosen model for SRL aerothermal design simulations.

Mars entry↗

“Best” Iterative Coupled-Cluster Triples Model? More Evidence for 3CC

To follow up on the unexpectedly good performance of several coupled-cluster models with approximate inclusion of 3-body clusters we performed a more complete assessment of the 3CC method for accurate computational thermochemistry in the standard HEAT framework. New spin-integrated implementation of the 3CC method applicable to closed- and open-shell systems utilizes a new automated toolchain for derivation, optimization, and evaluation of operator algebra in many-body electronic structure. We found that with a double-ζ basis set the 3CC correlation energies and their atomization energy contributions are almost always more accurate (with respect to the CCSDTQ reference) than the CCSDT model as well as the standard CCSD(T) model. The mean absolute errors in cc-pVDZ {3CC, CCSDT, and CCSD(T)} electronic (per valence electron) and atomization energies relative to the CCSDTQ reference for the HEAT data set, were {24, 70, 122} μE h /e and {0.46, 2.00, 2.58} kJ/mol, respectively. The mean absolute errors in the complete-basis-set limit {3CC, CCSDT, and CCSD(T)} atomization energies relative to the HEAT model reference, were {0.52, 2.00, and 1.07} kJ/mol, The significant and systematic reduction of the error by the 3CC method and its lower cost than CCSDT suggests it as a viable candidate for post- CCSD(T) thermochemistry applications, as well as the preferred alternative to CCSDT in general.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accelerating Thermochemical Equilibrium Calculations for Nuclear Reactor Applications

Thermochemical properties play a key role in modeling and simulation of several key phenomena in nuclear reactors. There has been an increasing interest in incorporating CALPHAD-based formulations in multiphysics simulations including for Molten Salt Reactors where knowledge of phase evolution of the salt and the chemical potentials of various elements are of utmost importance in source term analyses and redox control. However, the size of such simulations is often limited by the high computational cost of full thermodynamic equilibrium calculations. This work discusses the current efforts aimed at accelerating thermochemical equilibrium calculations for multiphysics simulations performed using the open-source finite element / finite volume code Multiphysics Object Oriented Simulation Environment (MOOSE) [1]. While several methods have been proposed for accelerating phase equilibrium calculations [2], most focus on relatively small systems and often rely on a- priori knowledge of the state-space of the system. Nuclear materials, however, are often multi-component systems owing to the evolution of composition under irradiation and an approach based on a-priori mapping of phase diagram is often not enough. This work is aimed at demonstrating an on-the-fly surrogate modeling framework that uses active learning to reduce the number of full equilibrium calculations that must be performed. By combining with efficient coupling approaches, the surrogate framework helps in reducing the computational cost of thermodynamic equilibrium informed multiphysics simulations of nuclear materials. The performance is benchmarked against full coupling with the thermochemistry library Thermochimica [3]. This work uses a machine learning based approach for constructing surrogate models to predict the stable phases in a multicomponent system. The surrogates were constructed using neural networks and Gaussian process classification. In this work, we compare the relative performance of the two methods. We also demonstrate the use of caching previous calculations by interpolating the values from nearest neighbors. References [1] Lindsay, A.D., et al. "2.0 – MOOSE: Enabling massively parallel multiphysics simulation", SoftwareX, 20 (2022): 101202. [2] Roos, W.A. and Zietsman J.H. "Accelerating complex chemical equilibrium calculations – A Review", Calphad, 77 (2022): 102380. [3] Piro, M.H.A., et al. "The thermochemistry library Thermochimica", Computational Materials Science, 67 (2013): 266-272.

36 MATERIALS SCIENCE↗

Association Kinetics for Perfluorinated n -Alkyl Radicals

Radical-radical reaction channels are important in the pyrolysis and oxidation chemistry of perfluoroalkyl substances (PFAS). In particular, unimolecular dissociation reactions within unbranched n-perfluoroalkyl chains, and their corresponding reverse barrierless association reactions, are expected to be significant contributors to the gas-phase thermal decomposition of families of species such as perfluorinated carboxylic acids and perfluorinated sulfonic acids. Unfortunately, experimental data for these reactions are scarce and uncertain. Furthermore, obtaining reliable theoretical predictions for such reactions is a laborious and computationally intensive task. Here, in this work, the chemical kinetics of the various association/decomposition reactions producing/decomposing the C 2 -C 4 series of unbranched n-perfluoroalkanes (C 2 F 6 , C 3 F 8 , and C 4 F 10 ) are examined using state-of-the-art ab initio transition-state-theory-based master-equation calculations. The variable-reaction-coordinate transition-state theory (VRC-TST) formalism is employed in computing the microcanonical and canonical rates for the association reactions. Reaction thermochemistry is obtained via composite quantum chemistry calculations and the laddering of error-canceling reaction schemes via a connectivity-based hierarchy approach employing ANL1/ANL0-style reference energies. Lennard-Jones collision model parameters for the considered systems were estimated by a direct dynamics approach, and collisional energy transfer parameters were obtained from analogies to systems of similar size and heavy-atom connectivity. A one-dimensional master equation approach was used to convert the microcanonical rate coefficients from the VRC-TST analysis into temperature- and pressure-dependent rate constants for the association reactions and the reverse dissociation reactions. The data are reported in standardized formats for usage in comprehensive chemical kinetic models for PFAS thermal destruction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.

Correlation energy↗

Catalytic production of δ-valerolactone (DVL) from biobased 2-hydroxytetrahydropyran (HTHP) – Combined experimental and modeling study

δ-Valerolactone (DVL) is a five-carbon (C5) cyclic ester that can undergo ring-opening polymerization to yield high-performance, biocompatible polyesters. But current market prices of C5 chemicals like DVL are very high due to poor availability of C5 feedstock in petroleum. Herein, we demonstrate a novel route to DVL synthesis via dehydrogenation of biomass-derived 2-hydroxytetrahydropyran (HTHP) over Cu/SiO 2 without the use of toxic reagents. Since HTHP exists in thermal equilibrium with 3,4-dihydropyran (DHP) via dehydration, and with 2,2’-oxybis(tetrahydropyran) and 5-(tetrahydropyran-2-yloxy)pentanal via acetalization, we have also determined the thermochemistry (ΔH rxn and ΔG rxn ) of each competing reaction using density functional theory (DFT) calculations at the M06-2X/cc-pVTZ level. Finally, by developing a kinetic model of all 8 reactions involved, we have achieved 84% selectivity to DVL at 150°C in a packed bed reactor for over 72 hours of time on stream.

42 ENGINEERING↗

An experimental, theoretical and kinetic modeling study of the N 2 O-H 2 system: Implications for N 2 O + H

The reaction of N 2 O with His the key step in consumption of nitrous oxide in thermal processes. The major product channel is N 2 + OH, while NH + NO constitute minor products. In addition, a pathway involving HNNO, initiated by N 2 O + H (+M)$\rightleftarrows$HNNO (+M) (R3, R4), has been inferred from experiment and theory by Burke and coworkers. At longer reaction times, the reaction may reach partial equilibration, and in addition to k 3 and k 4 the importance of this channel depends on the thermodynamic properties of HNNO and its consumption reactions, mainly HNNO + H. In the present work, we re-examined the thermochemistry of HNNO and calculated rate constants and branching fractions for the HNNO + H reaction. Experiments on the N 2 O-H 2 system were conducted in a high-pressure flow reactor at 100 atm as a function of temperature (600-925 K) and stoichiometry and explained in terms of an updated chemical kinetic model. The results support the importance of the HNNO pathway, which results in inhibition of N 2 O consumption and formation of NH 3 . In addition, selected literature results on the N 2 O-H 2 system are re-examined and the implications for the other product channels of N 2 O + H, in particular NH + NO, are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretically Informed Kinetics (ThInK): Establishing a modern C 0 -C 3 mechanism for combustion modeling

In contrast to the adage “Models are to be used, not believed”, combustion kinetics models have been intended to be predictive in nature. Theoretical chemical kinetics is now understood to provide a firm foundation for the reaction parameters, thereby facilitating predictive simulations of chemical reactivity, even in regimes that are poorly characterized by chemical kinetic and/or combustion experiments. In this study, we describe a theory-informed kinetics model (ThInK) for small molecule combustion chemistry (H 2 and C 1 – C 3 species) that is based on the prodigious use of theoretical predictions for reaction rate coefficients, thermochemistry, and transport parameters. The distinct features of this kinetics model, which was developed over the course of several decades, are illustrated through simulations of flame propagation and auto-ignition.

Energy transfer↗

Vanadium Substitution Dictates H Atom Uptake at Lindqvist-type Polyoxotungstates

Understanding how modification of molecular structures changes the thermochemistry of H atom uptake can provide design criteria for the formation of highly active catalysts for reductive transformations. Herein, we describe the effect of doping an atomically precise polyoxotungstate with vanadium on proton-coupled electron transfer (PCET) reactivity. The Lindqvist-type polyoxotungstate [W 6 O 19 ] 2– displays reversible redox chemistry, which was found to be unchanged in the presence of acid, indicating an inability to couple reduction with protonation. However, the incorporation of a single vanadium center into the structure significantly changes the reactivity, and the potential required for one-electron reduction of [VW 5 O 19 ] 3– was shown to vary with the strength of the acid added. Construction of a potential-pK a diagram allowed assessment of the thermodynamics of H atom uptake, indicating BDFE(O–H) ≈ 64 kcal/mol, while chemical synthesis of the reduced/protonated derivative (TBA) 3 [VW 5 O 19 H] was used to probe the position of protonation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗