Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermochemistry analysis”

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

CFD Analysis of Hypersonic Flowfields With Surface Thermochemistry and Ablation

In the past forty years much progress has been made in computational methods applied to the solution of problems in spacecraft hypervelocity flow and heat transfer. Although the basic thermochemical and physical modeling techniques have changed little in this time, several orders of magnitude increase in the speed of numerically solving the Navier-Stokes and associated energy equations have been achieved. The extent to which this computational power can be applied to the design of spacecraft heat shields is dependent on the proper coupling of the external flow equations to the boundary conditions and governing equations representing the thermal protection system in-depth conduction, pyrolysis and surface ablation phenomena. A discussion of the techniques used to do this in past problems as well as the current state-of-art is provided. Specific examples, including past missions such as Galileo, together with the more recent case studies of ESA/Rosetta Sample Comet Return, Mars Pathfinder and X-33 will be discussed. Modeling assumptions, design approach and computational methods and results are presented.

Henline, W. D.↗

Ablation, Thermal Response, and Chemistry Program for Analysis of Thermal Protection Systems

In previous work, the authors documented the Multicomponent Ablation Thermochemistry (MAT) and Fully Implicit Ablation and Thermal response (FIAT) programs. In this work, key features from MAT and FIAT were combined to create the new Fully Implicit Ablation, Thermal response, and Chemistry (FIATC) program. FIATC is fully compatible with FIAT (version 2.5) but has expanded capabilities to compute the multispecies surface chemistry and ablation rate as part of the surface energy balance. This new methodology eliminates B' tables, provides blown species fractions as a function of time, and enables calculations that would otherwise be impractical (e.g. 4+ dimensional tables) such as pyrolysis and ablation with kinetic rates or unequal diffusion coefficients. Equations and solution procedures are presented, then representative calculations of equilibrium and finite-rate ablation in flight and ground-test environments are discussed.

Milos, Frank S.↗

Theoretical analysis of oxygen diffusion at startup in an alkali metal heat pipe with gettered alloy walls

The diffusion of oxygen into, or out of, a gettered alloy exposed to oxygenated alkali liquid metal coolant, a situation arising in some high temperature heat transfer systems, was analyzed. The relation between the diffusion process and the thermochemistry of oxygen in the alloy and in the alkali metal was developed by making several simplifying assumptions. The treatment is therefore theoretical in nature. However, a practical example pertaining to the startup of a heat pipe with walls of T-111, a tantalum alloy, and lithium working fluid illustrates the use of the figures contained in the analysis.

Tower, L. K.↗

Updates and Modernization of the Chemical Equilibrium with Applications (CEA) Code

NASA’s Chemical Equilibrium with Applications (CEA) code is a foundational tool for propulsion system analysis. It provides equilibrium chemistry, rocket performance, shock, and detonation calculations used across NASA and the broader aerospace community. NASA Engineering and Safety Center (NESC) Activity TI-22-01730 modernized the legacy CEA2 Fortran code into CEA v3, a Fortran 2008, object-oriented software package with expanded interface support, updated thermochemical data, improved maintainability, and substantially improved workflow integration. The modernized code preserves backward compatibility with legacy CEA input workflows while enabling direct use from modern analysis environments, including Python, C, MATLAB, and automated design studies.

Combustion↗

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↗

Magnetohydrodynamics (MHD) Aerocapture System for Enabling Faster-Larger Planetary Science & Human Exploration Missions

Since our completing the NIAC Phase I NIAC on this Advanced Aerocapture System, NASA Langley Research Center has funded or supported a number of studies and code enhancements through its Center Innovation Fund (CIF) and NASA’s NSTGRO and Internship Programs to mature the analysis capabilities and quantify the merits of the MHD Aerocapture System technology. These efforts have resulted in a plug and play analysis capability for assessing MHD aerocapture system performance for arrival at many planetary bodies of interest. Our efforts have especially focused on the potential mass savings for improving the capacity for science observations at Neptune and Triton. A re-cent Forbes article published “‘Orbital mechanics is probably going to decide for us whether we go to Uranus or Neptune because we need to flyby Jupiter,’ said Kunio Sayanagi at Hampton University, Virginia, who also worked on the Neptune Odyssey proposal…. Exactly when a mission can be sent to Uranus, or Neptune, depends on the relative position of Jupiter, which can help give a spacecraft a gravitational slingshot. That drastically shortens the cruise phase.” [1] Since shortening the cruise phase is important for these science missions, any mass savings enabled by the MHD Aerocapture System could be reallocated to increasing Thermal Protection System mass to allow faster arrival speeds and/or for onboarding additional payloads such as science instruments, batteries, or propellant for conducting more science for longer durations in the desired orbits. The analysis steps and codes for conducting trades and sizing vehicles for aerocapture are as follows: Step 1 is to conduct aeroheating analysis using LAURA of the selected entry vehicle shape to identify locations on the forebody where ionization and flow velocity are sufficient for producing Lo-rentz forces. LAURA is a multiblock structured grid finite-volume CFD solver developed at the NASA Langley Research Center. [2] LAURA has been used for aerothermal analysis support of the entry, de-scent and landing (EDL) phase of interplanetary missions over the last three decades [3-7]. Step 2 is to port the LAURA results into CFDWARP to calcu-late electrical and thermal conductivities of ionized flow for sizing MHD patch system and calculating Lorentz forces needed for controls analysis. CFDWARP is a CFD code that uses advanced nu-merical methods that enable the simulation of the full coupling between the aerodynamics, the magne-tohydrodynamics, and the non-neutral plasma sheaths. CFDWARP has the unique capability to simulate efficiently the non-neutral sheaths (near the electrodes) in coupled form with the quasi-neutral bulk MHD flow [8-11]. Step 3 is to link re-sults from LAURA and CFDWARP into POST2 for calculating entry trajectories and comparing MHD control results with other aerodynamic control strategies. The Program to Optimize Simulated Tra-jectories II (POST2) is a generalized point mass, discrete parameter targeting and optimization pro-gram. POST2 provides the capability to target and optimize point mass trajectories for multiple pow-ered or un-powered vehicles near an arbitrary rotat-ing, oblate planet [12]. Step 4: TPS sizing was per-formed using the Fully Implicit Ablation and Ther-mal-response code (FIAT) tool which computes the transient one-dimensional thermal response and surface thermochemistry of a multilayer stackup of thermal protection, bonding, and structural materi-als subject to aeroheating on one surface [13]. The sizing and margining methodology used was based on the approach documented by Mahzari and Milos [14] for the dual-layered heatshield for extreme entry environment technology (DL-HEEET) TPS concept. TPS analysis utilizes trajectory information from POST2. Using this step-wise plug and play MHD Aerocapture performance assessment process, our analysis targets a Neptune aerocapture trajectory that will place the spacecraft in an observation orbit for Triton. [15]. Magnetohydrodynamic (MHD) control of a 4.5-meter diameter MSL-style capsule resulted in TPS mass savings of nearly 2000 kg when using an MHD system mass of under 200 kg. The flight path for a vehicle using the MHD control strategy has a much lower heat rate and heat load compared to the conventional aerodynamic aerocapture strategies known as bank angle con-trolled (BAC) and direct force controlled (DFC). Both BAC and DFC have heat rates significantly greater than 1500 W/cm2 typically used as an upper limit for PICA. Thus, DL-HEEET TPS concept was required for the BAC and DFC control strategies. However, considering the more benign environ-ments for the MHD case, additional TPS concepts with improved mass efficiency were also assessed. PICA was considered for the MHD controlled strat-egy since the maximum heat rate was well within the limits (<1500 W/cm2) of PICA. TPS sizing re-sulted in a significant mass reduction. The PICA layer for this sizing case was about 7.8 cm. As a point of reference, the Mars 2020 mission, which used this same PICA concept, had a PICA thickness of 3.18 cm [16]. The trajectories used for the TPS sizing originat-ed from the POST2 simulations. The current, I, to an electromagnet configuration can be manipulated to allow for active control of the vehicle. Manipula-tion of the current, I, changes the magnetic field, B, which affects the Lorentz force and therefore the MHD drag force on the vehicle. Our analysis in-cluded both open-loop and close-loop control. Closed-loop control will enable improved overall performance when taking into account mission level uncertainties, such as interplanetary delivery errors and atmospheric modeling uncertainties. The open-loop and closed-loop MHD control cases do not dip as deep into the atmosphere as the aerodynamic cases. Three types of aerodynamic-only approach-es are investigated: bank angle modulation (BAM), director force control (DFC), and Drag Modulated. BAM and DFC make use of vehicle aerodynamic angles to steer the vehicle. Thus, changing the aer-odynamic forces acting on the vehicle for control, aerodynamic drag modulated case requires a vary-ing drag area to modulate the drag force. The MHD drag modulated case modulates MHD generated drag force that adds to the aerodynamic drag. This higher atmospheric activation of drag forces by the MHD patch results in significantly less heat flux on the vehicle. The MHD technology will enable shorter cruise times and deceleration of larger payloads for increasing the capacity for science at the Ice Giants or for returning astronauts to Earth from cislunar space or from Mars. The purpose of this presentation is to provide more details about this work and to highlight plans for further research and development including a flight demonstration.

R. W. Moses↗

Study of thermal stability and degradation of fire resistant candidate polymers for aircraft interiors

The thermochemistry of bismaleimide resins and phenolphthalein polycarbonate was studied. Both materials are fire-resistant polymers and may be suitable for aircraft interiors. The chemical composition of the polymers has been determined by nuclear magnetic resonance and infrared spectroscopy and by elemental analysis. Thermal properties of these polymers have been characterized by thermogravimetric analyses. Qualitative evaluation of the volatile products formed in pyrolysis under oxidative and non-oxidative conditions has been made using infrared spectrometry. The residues after pyrolysis were analyzed by elemental analysis. The thermal stability of composite panel and thermoplastic materials for aircraft interiors was studied by thermogravimetric analyses.

Hsu, M. T. S.↗

Detailed course of the O + HO2 reaction

The elementary radical-radical reactions that are of such importance in laser, astrophysical, atmospheric, and combustion processes are increasingly open to direct rate measurement. Attention is presently given to the O + HO2 reaction, which has only one exothermic product channel: OH + O2. Isotopic product analysis was conducted spectroscopically using laser-induced fluorescence of OH. It is found that the simple atom-radical reaction takes place by way of the rapid formation and breakup of an unexpected radical intermediate, thereby raising questions as to the formation, thermochemistry, and decomposition of poorly known products in many other radical-radical reaction systems.

Sridharan, U. C.↗

Thermochemistry of Protective Coatings and Molten Silicate Debris

The durability of gas-turbine engine components can be significantly affected by the ingestion of siliceous particles, which can melt at high temperature and corrode protective coatings that are essential for long life requirements. The silicate debris consists mainly of CaO-MgO-Al2O3-SiO2 (CMAS) and is usually ingested by aircraft engines during and after take-off, sticking to their hot surfaces and resulting in the formation of calcium rare-earth silicate oxyapatites. The thermochemistry of coatings and their reaction products with molten silicate debris are crucial to understand in order to improve the durability of gas-turbine engines. Here we discuss results of high temperature drop solution calorimetry, drop-and-catch calorimetry (DnC) and differential thermal analysis (DTA) techniques for the thermodynamic properties of both thermal barrier coatings (TBCs) and environmental barrier coatings (EBCs) and their reaction with CMAS compositions. The enthalpies of solution of Y2Si2O7, Yb2Si2O7, 31YSZ, and 16RESZ based coatings and the oxyapatite are moderately positive. However, oxyapatite formation is only favorable over coating dissolution in terms of enthalpy for 7YSZ. The enthalpies of mixing between the coatings and the molten silicate are less exothermic for Yb2Si2O7 and CaYb4Si3O13 than for 7YSZ, indicating lower energetic stability of the latter against molten silicate corrosion. We also report for the first time the calorimetric measurements of the enthalpies of formation of rare-earth silicate based EBC coatings and oxyapatites (rare-earth, RE = Y, Yb, Gd, Dy, Er, Nd and Sm).

Costa, Gustavo↗

Supersonic flow of chemically reacting gas-particle mixtures. Volume 1: A theoretical analysis and development of the numerical solution

A numerical solution for chemically reacting supersonic gas-particle flows in rocket nozzles and exhaust plumes was described. The gas-particle flow solution is fully coupled in that the effects of particle drag and heat transfer between the gas and particle phases are treated. Gas and particles exchange momentum via the drag exerted on the gas by the particles. Energy is exchanged between the phases via heat transfer (convection and/or radiation). Thermochemistry calculations (chemical equilibrium, frozen or chemical kinetics) were shown to be uncoupled from the flow solution and, as such, can be solved separately. The solution to the set of governing equations is obtained by utilizing the method of characteristics. The equations cast in characteristic form are shown to be formally the same for ideal, frozen, chemical equilibrium and chemical non-equilibrium reacting gas mixtures. The particle distribution is represented in the numerical solution by a finite distribution of particle sizes.

Penny, M. M.↗

Space shuttle SRM plume expansion sensitivity analysis

The exhaust plumes of the space shuttle solid rocket motors can have a significant effect on the base pressure and base drag of the shuttle vehicle. A parametric analysis was conducted to assess the sensitivity of the initial plume expansion angle of analytical solid rocket motor flow fields to various analytical input parameters and operating conditions. The results of the analysis are presented and conclusions reached regarding the sensitivity of the initial plume expansion angle to each parameter investigated. Operating conditions parametrically varied were chamber pressure, nozzle inlet angle, nozzle throat radius of curvature ratio and propellant particle loading. Empirical particle parameters investigated were mean size, local drag coefficient and local heat transfer coefficient. Sensitivity of the initial plume expansion angle to gas thermochemistry model and local drag coefficient model assumptions were determined.

Smith, S. D.↗

Kinetics and Thermochemistry of Reversible Adduct Formation in the Reaction of Cl((sup 2)P(sub J)) with CS2

Reversible adduct formation in the reaction of Cl((sup 2)P(sub J)) with CS2 has been observed over the temperature range 193-258 K by use of time-resolved resonance fluorescence spectroscopy to follow the decay of pulsed-laser-generated Cl((sup 2)P(sub J)) into equilbrium with CS2Cl. Rate coefficients for CS2Cl formation and decomposition have been determined as a function of temperature and pressure; hence, the equilbrium constant has been determined as a function of temperature. A second-law analysis of the temperature dependence of Kp and heat capacity corrections calculated with use of an assumed CS2Cl structure yields the following thermodynamic parameters for the association reaction: Delta-H(sub 298) = -10.5 +/- 0.5 kcal/mol, Delta-H(sub 0) = -9.5 +/- 0.7 kcal/mol, Delta-S(sub 298) = -26.8 +/- 2.4 cal/mol.deg., and Delta-H(sub f,298)(CS2Cl) = 46.4 +/- 0.6 kcal/mol. The resonance fluorescence detection scheme has been adapted to allow detection of Cl((sup 2)P(sub J)) in the presence of large concentrations of O2, thus allowing the CS2Cl + Cl + O2 reaction to be investigated. We find that the rate coefficient for CS2Cl + O2 reaction via all channels that do not generate Cl((sup 2)P(sub J)) is less than 2.5 x 10(exp-16) cu cm/(molecule.s) at 293 K and 300-Torr total pressure and that the total rate coefficient is less than 2 x 10 (exp -15) cu cm/(molecule.s) at 230 K and 30-Torr total pressure. Evidence for reversible adduct formation in the reaction of Cl((sup 2)P(sub J)) with COS was sought but not observed, even at temperatures as low as 194 K.

Nicovich, J. M.↗

Thermochemistry and Kinetics of the Cl+O2 Association Reaction

Laser flash photolysis of Cl2/O2 mixtures has been employed in conjunction with Cl((sup 2)P(sub 3/2)) detection by time-resolved fluorescence spectroscopy to investigate equilibration kinetics for the reactions Cl + O2 + O is in equilibrium with ClOO + O2 at temperatures of 181-200 K and O2 pressures of 15-40 Torr. The third-order rate coefficient for the association reaction at 186.5 +/- 5.5 K is (8.9 +/- 2.9) x 10(exp -33) cm(exp 6) molecule(exp -2) s(exp -1) and the equilibrium constant (K(p)) at 185.4 K is 18.9 atm(exp -1) (factor of 1.7 uncertainty). A third law analysis of our data leads to a value for the Cl-OO bond dissociation energy of 4.76 +/- 0.49 kcal mol(exp -1).

Nicovich, J. M.↗

Validation of a Thermo-Ablative Model of Elastomeric Internal Insulation Materials

In thermo-ablative material modeling, as in many fields of analysis, the quality of the existing models significantly exceeds that of the experimental data required for their validation. In an effort to narrow this gap, a laboratory-scale internal insulation test bed was developed that exposes insulation samples to realistic solid rocket motor (SRM) internal environments while being instrumented to record real-time rates of both model inputs (i.e., chamber pressure, total surface heat flux, and radiative heat flux) as well as model outputs (i.e., material decomposition depths (MDDs) and in-depth material temperatures). In this work, the measured SRM internal environment parameters were used in conjunction with equilibrium thermochemistry codes as inputs to one-dimensional thermo-ablative models of the PBINBR and CFEPDM insulation samples used in the lab-scale test firings. The computed MDD histories were then compared with those deduced from real-time X-ray radiography of the insulation samples, and the calculated in-depth temperatures were compared with those measured by embedded thermocouples. The results of this exercise emphasize the challenges of modeling and testing elastomeric materials in SRM environments while illuminating the path forward to improved fidelity.

Martin, Heath T.↗

Hydrogen-silicon carbide interactions

A study of the thermochemistry and kinetics of hydrogen environmental attack of silicon carbide was conducted for temperatures in the range from 1100 C to 1400 C. Thermodynamic maps based on the parameters of pressure and oxygen/moisture content were constructed. With increasing moisture levels, four distinct regions of attack were identified. Each region is defined by the thermodynamically stable solid phases. The theoretically stable solid phases of Region 1 are silicon carbide and silicon. Experimental evidence is provided to support this thermodynamic prediction. Silicon carbide is the single stable solid phase in Region 2. Active attack of the silicon carbide in this region occurs by the formation of gases of SiO, CO, CH4, SiH4, and SiH. Analysis of the kinetics of reaction for Region 2 at 1300 C show the attack of the silicon carbide to be controlled by gas phase diffusion of H2O to the sample. Silicon carbide and silica are the stable phases common to Regions 3 and 4. These two regions are characterized by the passive oxidation of silicon carbide and formation of a protective silica layer.

Eckel, Andrew J.↗

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↗

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↗

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↗