Thermokinetic Modeling in an Aerospace Research and Production Environment
Thermokinetic modeling has proven its value for measuring thermal response and reaction rates of a wide variety of aerospace materials.
SEARCH · Engineering Papers
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.
Thermokinetic modeling has proven its value for measuring thermal response and reaction rates of a wide variety of aerospace materials.
Current air revitalization technology onboard the International Space Station (ISS) cannot provide complete closure of the oxygen and hydrogen loops. This makes re-supply necessary, which is possible for missions in low Earth orbit (LEO) like the ISS, but unviable for long term space missions outside LEO. In comparison, Bosch technology reduces carbon dioxide with hydrogen, traditionally over a steel wool catalyst, to create water and solid carbon. The Bosch product water can then be fed to the oxygen generation assembly to produce oxygen for crew members and hydrogen necessary to reduce more carbon dioxide. Bosch technology can achieve complete oxygen loop closure, but has many undesirable factors that result in a high energy, mass, and volume system. Finding a different catalyst with an equal reaction rate at lower temperatures with less catalyst mass and longer lifespan would make a Bosch flight system more feasible. Developmental testing of alternative catalysts for the Bosch has been performed using the Horizontal Bosch Test Stand. Nickel foam, nickel shavings, and cobalt shavings were tested at 500 C and compared to the original catalyst, steel wool. This paper presents data and analysis on the performance of each catalyst tested at comparable temperatures and recycle flow rates.
Oxidation of high-temperature aerospace materials is a universal issue for combustion-path components in turbine or rocket engines. In addition to the question of the consumption of material due to growth of protective scale at use temperatures, there is also the question of cyclic effects and spallation of scale on cooldown. The spallation results in the removal of part of the protective oxide in a discontinuous step and thereby opens the way for more rapid oxidation upon reheating. In experiments, cyclic oxidation behavior is most commonly characterized by measuring changes in weight during extended time intervals that include hundreds or thousands of heating and cooling cycles. Weight gains occurring during isothermal scale-growth processes have been well characterized as being parabolic or nearly parabolic functions of time because diffusion controls reaction rates. In contrast, the net weight change in cyclic oxidation is the sum of the effects of the growth and spallation of scale. Typically, the net weight gain in cyclic oxidation is determined only empirically (that is, by measurement), with no unique or straightforward mathematical connection to either the rate of growth or the amount of metal consumed. Thus, there is a need for mathematical modeling to infer spallation mechanisms. COSP is a computer program that models the growth and spallation processes of cyclic oxidation on the basis of a few elementary assumptions that were discussed in COSP: A Computer Model of Cyclic Oxidation, Oxidation of Metals, vol. 36, numbers 1 and 2, 1991, pages 81-112. Inputs to the model include the selection of an oxidation-growth law and a spalling geometry, plus oxide-phase, growth-rate, cycle-duration, and spall-constant parameters. (The spalling fraction is often shown to be a constant factor times the existing amount of scale.) The output of COSP includes the net change in weight, the amounts of retained and spalled oxide, the total amounts of oxygen and metal consumed, and the terminal rates of weight loss and metal consumption.
A computational study has shown that by injecting drops in jets of gas having square, elliptical, triangular, or other noncircular injection cross sections, it should be possible to increase (relative to comparable situations having circular cross section) the entrainment and dispersion of liquid drops. This finding has practical significance for a variety of applications in which it is desirable to increase dispersion of drops. For example, in chemical-process sprays, increased dispersion leads to increases in chemical- reaction rates; in diesel engines, increasing the dispersion of drops of sprayed fuel reduces the production of soot; and in household and paint sprays, increasing the dispersion of drops makes it possible to cover larger surfaces. It has been known for some years that single-phase fluid jets that enter flow fields through noncircular inlets entrain more fluid than do comparable jets entering through circular inlets. The computational study reported here was directed in part toward determining whether and how this superior mixing characteristic of noncircular single phase jets translates to a similar benefit in cases of two-phase jets (that is, sprays). The study involved direct numerical simulations of single- and two-phase free jets with circular, elliptical, rectangular, square, and triangular inlet cross sections. The two-phase jets consisted of gas laden with liquid drops randomly injected at the inlets. To address the more interesting case of evaporating drops, the carrier gas in the jets was specified to be initially unvitiated by the vapor of the liquid chemical species and the initial temperature of the drops was chosen to be smaller than that of the gas. The mathematical model used in the study was constructed from the conservation equations for the two-phase flow and included complete couplings of mass, momentum, and energy based on thermodynamically self-consistent specification of the enthalpy, internal energy, and latent heat of vaporization of the vapor.
We have been studying Fischer-Tropsch type (FTT) reactions as a source for organic materials both in the gas phase of the solar nebula and incorporated into primitive comets and asteroids for almost 10 years, and over this time our concept has evolved greatly from the standard "catalytic" model to a much more robust chemical scenario. Our simulations have been conducted at temperatures that are much higher than we like, primarily for practical reasons such as the timescale of individual reactions, and we are just starting a series of measurements to allow us to measure reaction rates at temperatures from 873K down to as low as 373K. We have preliminary data on the carbon (d13C = -50) & nitrogen (d15N = +9.5) isotopic fractionation at 873K, but not on materials produced at lower temperature. Isotope values are on the VPDB scale for carbon and vs. Air for nitrogen. We have also investigated the noble gas trapping efficiency of the FTT process by adding a small amount of a noble gas mix to our standard synthesis mix. The noble gas ratio is 49:49:1:1::Ne:Ar:Kr:Xe. Xe and Kr are trapped at 873K and are more efficiently trapped at 673K with no isotopic fractionation at either temperature. Ar trapping is detected at 673K, but not at 873K. Ne has not yet been observed in our samples. The solar nebula was an extremely complex system, mixing materials from the innermost regions out to well into the zones where comets formed and thus mixing highly processed nebular materials with grains and coatings formed before the nebula began to collapse. Laboratory studies may provide the means to separate such diverse components based on carbon or nitrogen isotopic fractionation or the quantities of noble gases trapped in grain coatings and their thermal release patterns, among other observables. The ultimate goal of laboratory synthesis of nebular analogs is to provide the means to identifY the conditions under which natural samples were formed and the signatures of subsequent metamorphic events.
The analysis of mutual ozone-temperature variations can provide useful information on their interdependencies relative to the photochemistry and dynamics governing their behavior. Previous studies have mostly been based on satellite measurements taken at a fixed local time in the stratosphere and lower mesosphere. For these data, it is shown that the zonal mean ozone amounts and temperatures in the lower stratosphere are mostly positively correlated, while they are mostly negatively correlated in the upper stratosphere and in the lower mesosphere. The negative correlation, due to the dependence of photochemical reaction rates on temperature, indicates that ozone photochemistry is more important than dynamics in determining the ozone amounts. In this study, we provide new results by extending the analysis to include diurnal variations over 24 hrs of local time, and to larger spatial regimes, to include the upper mesosphere and lower thermosphere (MLT). The results are based on measurements by the SABER instrument on the TIMED satellite. For mean variations (i.e., averages over local time and longitude) in the MLT, our results show that there is a sharp reversal in the correlation near 80 km altitude, above which the ozone mixing ratio and temperature are mostly positively correlated, while they are mostly negatively correlated below 80 km. This is consistent with the view that above -80 km, effects due to dynamics are more important compared to photochemistry. For diurnal variations, both the ozone and temperature show phase progressions in local time, as a function of altitude and latitude. For temperature, the phase progression is as expected, as they represent migrating tides. For day time ozone, we also find regular phase progression in local time over the whole altitude range of our analysis, 25 to 105 km, at least for low latitudes. This was not previously known, although phase progressions had been noted by us and by others at lower altitudes. For diurnal variations, we find that between about 40 and 65 km, the ozone amounts and temperatures are mostly negatively correlated or neutral, while below approx. 40 km they are mostly positively correlated or neutral. The correlations are less systematic and less robust than for correlations of the mean. At altitudes above approx.65 km, the correlations are more complex, and depend on the tidal temperature variations. For the diurnal case, consideration needs to be given to transport by thermal tides and to the efficacy of response times of ozone concentrations and temperature to each other.
In 1929 London 1 published a very approximate solution of the Schroedinger equation for a system of chemical interest: H3. To the extent that chemistry can be regarded as existing separately from physics, this was a landmark in the history of chemistry, comparable in importance to the landmark in the history of physics marked by the appearance of the Heitler-London equation for H2. The expression for H3, was, of necessity, even less accurate than that for H2, but chemists, like the habitual poor, were accustomed to this sort of misfortune. Together with the physicists they enjoyed the sensation of living in a renaissance. The physicists still could not calculate a great deal that was of interest to them, and the chemists could calculate less, but both could now dream. It would be too easy to say that their dreams were dreams of unlimited computer time. Their dreams were a lot more productive than that. Two years after London published his equation, H. Eyring and M. Polanyi obtained the first numerical energy surface for H3. They infused the London equation with a measure of empiricism to produce an energy surface which, whether or not it was correct in its details, provided a basis for further speculations of an important sort. The existence of a tangible energy surface in 1931 stimulated speculation along two different lines. The following year Pelzer and Wigner used this London-Eyring-Polanyi (LEP) energy surface for a thermodynamic treatment of the reaction rate in H + H2. This important development reached its full flowering a few years later. In these remarks I shall be concerned with another line of development. A second more-or-less distinct category of speculation that began with (and, indeed, in) the 1931 paper has to do with the dynamics of individual reactive encounters under the influence of specified interaction potentials.
Simulation codes often utilize finite-dimensional approximation resulting in numerical error. Some examples include, numerical methods utilizing grids and finite-dimensional basis functions, particle methods using a finite number of particles. These same simulation codes also often contain sources of uncertainty, for example, uncertain parameters and fields associated with the imposition of initial and boundary data,uncertain physical model parameters such as chemical reaction rates, mixture model parameters, material property parameters, etc.
A computational effort was undertaken to analyze the details of fluid flow in Lean-Direct Injection (LDI) combustors for next-generation LDI design. The National Combustor Code (NCC) was used to perform reacting flow computations on single-element LDI injector configurations. The feasibility of using a reduced chemical-kinetics approach, which optimizes the reaction rates and species to model the emissions characteristics typical of lean-burning gas-turbine combustors, was assessed. The assessments were performed with Reynolds- Averaged Navier-Stokes (RANS) and Time-Filtered Navier Stokes (TFNS) time-integration, with a Lagrangian spray model with the NCC code. The NCC predictions for EINOx and combustor exit temperature were compared with experimental data for two different single-element LDI injector configurations, with 60deg and 45deg axially swept swirler vanes. The effects of turbulence-chemistry interaction on the predicted flow in a typical LDI combustor were studied with detailed comparisons of NCC TFNS with experimental data.
This work describes the direct simulation Monte Carlo (DSMC) investigation of Saturn entry probe scenarios and the influence of non-equilibrium phenomena on Saturn entry conditions. The DSMC simulations coincide with rarefied hypersonic shock tube experiments of a hydrogen-helium mixture performed in the Electric Arc Shock Tube (EAST) at NASA Ames Research Center. The DSMC simulations are post-processed through the NEQAIR line-by-line radiation code to compare directly to the experimental results. Improved collision cross-sections, inelastic collision parameters, and reaction rates are determined for a high temperature DSMC simulation of a 7-species H2-He mixture and an electronic excitation model is implemented in the DSMC code. Simulation results for 27.8 and 27.4 kms shock waves are obtained at 0.2 and 0.1 Torr respectively and compared to measured spectra in the VUV, UV, visible, and IR ranges. These results confirm the persistence of non-equilibrium for several centimeters behind the shock and the diffusion of atomic hydrogen upstream of the shock wave. Although the magnitude of the radiance did not match experiments and an ionization inductance period was not observed in the simulations, the discrepancies indicated where improvements are needed in the DSMC and NEQAIR models.
This work describes the direct simulation Monte Carlo (DSMC) investigation of Saturn entry probe scenarios and the influence of non-equilibrium phenomena on Saturn entry conditions. The DSMC simulations coincide with rarefied hypersonic shock tube experiments of a hydrogen-helium mixture performed in the Electric Arc Shock Tube (EAST) at the NASA Ames Research Center. The DSMC simulations are post-processed through the NEQAIR line-by-line radiation code to compare directly to the experimental results. Improved collision cross-sections, inelastic collision parameters, and reaction rates are determined for a high temperature DSMC simulation of a 7-species H2-He mixture and an electronic excitation model is implemented in the DSMC code. Simulation results for 27.8 and 27.4 km/s shock waves are obtained at 0.2 and 0.1 Torr, respectively, and compared to measured spectra in the VUV, UV, visible, and IR ranges. These results confirm the persistence of non-equilibrium for several centimeters behind the shock and the diffusion of atomic hydrogen upstream of the shock wave. Although the magnitude of the radiance did not match experiments and an ionization inductance period was not observed in the simulations, the discrepancies indicated where improvements are needed in the DSMC and NEQAIR models.
With advances in computational power and availability of distributed computers, the use of even the most complex of turbulent chemical interaction models in combustors and coupled analysis of combustors and turbines is now possible and more and more affordable for realistic geometries. Recent more stringent emission standards have enticed the development of more fuel-efficient and low-emission combustion system for aircraft gas turbine applications. It is known that the NOx emissions tend to increase dramatically with increasing flame temperature. It is well known that the major difficulty, when modeling the turbulence-chemistry interaction, lies in the high non-linearity of the reaction rate expressed in terms of the temperature and species mass fractions. The transport filtered density function (FDF) model and the linear eddy model (LEM), which both use local instantaneous values of the temperature and mass fractions, have been shown to often provide more accurate results of turbulent combustion. In the present, the time-filtered Navier-Stokes (TFNS) approach capable of capturing unsteady flow structures important for turbulent mixing in the combustion chamber and two different subgrid models, LEM-like and EUPDF-like, capable of emulating the major processes occurring in the turbulence-chemistry interaction will be used to perform reacting flow simulations of a selected test case. The selected test case from the Volvo Validation Rig was documented by Sjunnesson.
Ozone is an important air pollutant at the surface, and the third most important anthropogenic greenhouse gas in the troposphere. Since 1980, anthropogenic emissions of ozone precursors methane, non-methane volatile organic compounds, carbon monoxide and nitrogen oxides (NOx) have shifted from developed to developing regions. Emissions have thereby been redistributed equatorwards, where they are expected to have a stronger effect on the tropospheric ozone burden due to greater convection, reaction rates and NOx sensitivity. Here we use a global chemical transport model to simulate changes in tropospheric ozone concentrations from 1980 to 2010, and to separate the influences of changes in the spatial distribution of global anthropogenic emissions of short-lived pollutants, the magnitude of these emissions, and the global atmospheric methane concentration. We estimate that the increase in ozone burden due to the spatial distribution change slightly exceeds the combined influences of the increased emission magnitude and global methane. Emission increases in Southeast, East and South Asia may be most important for the ozone change, supported by an analysis of statistically significant increases in observed ozone above these regions. The spatial distribution of emissions dominates global tropospheric ozone, suggesting that the future ozone burden will be determined mainly by emissions from low latitudes.
We report measurements of bromine monoxide (BrO) and use an observationally constrained chemical box model to infer total gas-phase inorganic bromine (Br(sub y)) over the tropical western Pacific Ocean (tWPO) during the CONTRAST field campaign (January-February 2014). The observed BrO and inferred Bry profiles peak in the marine boundary layer (MBL), suggesting the need for a bromine source from sea-salt aerosol (SSA), in addition to organic bromine (CBry ). Both profiles are found to be C-shaped with local maxima in the upper free troposphere (FT). The median tropospheric BrO vertical column density (VCD) was measured as 1.6 x 10(exp 13) molec cm(exp -2), compared to model predictions of 0.9 x 10(exp 13) molec cm(exp -2) in GEOS-Chem (CBr(sub y) but no SSA source), 0.4 x 10(exp 13) molec cm(exp -2) in CAM-Chem (CBr(sub y) and SSA), and 2.1 x 10(exp 13) molec cm(exp -2) in GEOS-Chem (CBry and SSA). Neither global model fully captures the Cshape of the Br(sun y) profile. A local Br(sub y) maximum of 3.6 ppt (2.9-4.4 ppt; 95% confidence interval, CI) is inferred between 9.5 and 13.5 km in air masses influenced by recent convective outflow. Unlike BrO, which increases from the convective tropical tropopause layer (TTL) to the aged TTL, gas-phase Br(sub y) decreases from the convective TTL to the aged TTL. Analysis of gas-phase Br(sub y) against multiple tracers (CFC-11, H2O/O3 ratio, and potential temperature) reveals a Br(sub y) minimum of 2.7 ppt (2.3-3.1 ppt; 95% CI) in the aged TTL, which agrees closely with a stratospheric injection of 2.6 +/- 0.6 ppt of inorganic Br(sub y) (estimated from CFC-11 correlations), and is remarkably insensitive to assumptions about heterogeneous chemistry. Bry increases to 6.3 ppt (5.6-7.0 ppt; 95% CI) in the stratospheric "middleworld" and 6.9 ppt (6.5-7.3 ppt; 95% CI) in the stratospheric "overworld". The local Br(sub y) minimum in the aged TTL is qualitatively (but not quantitatively) captured by CAM-Chem, and suggests a more complex partitioning of gas-phase and aerosol Br(sub y) species than previously recognized. Our data provide corroborating evidence that inorganic bromine sources (e.g., SSA-derived gas-phase Br(sub y) ) are needed to explain the gas-phase Br(sub y) budget in the upper free troposphere and TTL. They are also consistent with observations of significant bromide in Upper Troposphere-Lower Stratosphere aerosols. The total Br(sub y) budget in the TTL is currently not closed, because of the lack of concurrent quantitative measurements of gas-phase Br(sub y) species (i.e., BrO, HOBr, HBr, etc.) and aerosol bromide. Such simultaneous measurements are needed to (1) quantify SSA-derived Br(sub y) in the upper FT, (2) test Br(sub y) partitioning, and possibly explain the gas-phase Br(sub y) minimum in the aged TTL, (3) constrain heterogeneous reaction rates of bromine, and (4) account for all of the sources of Br(sub y) to the lower stratosphere.
Glycine and methylamine are meteoritic water-soluble organic compounds that provide insights into the processes that occurred before, during, and after the formation of the Solar System. Both glycine and methylamine and many of their potential synthetic precursors have been studied in astrophysical environments via observations, laboratory experiments, and modeling. Despite these studies, the synthetic mechanisms for their formation leading to their occurrence in meteorites remain poorly understood. Typical 13C-isotopic values ((delta)13C) of meteoritic glycine and methylamine are 13C-enriched relative to their terrestrial counterparts; thus, analyses of their stable carbon isotopic compositions (13C/12C) may be used not only to assess terrestrial contamination in meteorites but also to provide information about their synthetic routes inside the parent body. Here, we examine potential synthetic routes of glycine and methylamine from a common set of precursors present in carbonaceous chondrite meteorites, using data from laboratory analyses of the well-studied CM2 Murchison meteorite. Several synthetic mechanisms for the origins of glycine and methylamine found in carbonaceous chondrites may be possible, and the prevalence of these mechanisms will largely depend on (a) the molecular abundance of the precursor molecules and (b) the levels of processing (aqueous and thermal) that occurred inside the parent body. In this work, we also aim to contextualize the current knowledge about gas-phase reactions and irradiated ice grain chemistry for the synthesis of these species through parent body processes. Our evaluation of various mechanisms for the origins of meteoritic glycine and methylamine from simple species shows what work is still needed to evaluate both the abundances and isotopic compositions of simpler precursor molecules from carbonaceous chondrites as well as the effects of parent body processes on those abundances and isotopic compositions. The analyses presented here combined with the indicated measurements will aid a better interpretation of quantitative analysis of reaction rates, molecular stability, and distribution of organic products from laboratory simulations of interstellar ices, astronomical observations, and theoretical modeling.
An accurate estimate of global hydroxyl radical (OH) abundance is important for projections of air quality, climate, and stratospheric ozone recovery. As the atmospheric mixing ratios of methyl chloroform (CH3CCl3) (MCF), the commonly used OH reference gas, approaches zero, it is important to find alternative approaches to infer atmospheric OH abundance and variability. The lack of global bottom-up emission inventories is the primary obstacle in choosing a MCF alternative. We illustrate that global emissions of long-lived trace gases can be inferred from their observed mixing ratio differences between the Northern Hemisphere (NH) and Southern Hemisphere (SH), given realistic estimates of their NH-SH exchange time, the emission partitioning between the two hemispheres, and the NH versus SH OH abundance ratio. Using the observed long-term trend and emissions derived from the measured hemispheric gradient, the combination of HFC-32 (CH2F2), HFC-134a (CH2FCF3, HFC-152a (CH3CHF2), and HCFC-22 (CHClF2), instead of a single gas, will be useful as a MCF alternative to infer global and hemispheric OH abundance and trace gas lifetimes. The primary assumption on which this multispecies approach relies is that the OH lifetimes can be estimated by scaling the thermal reaction rates of a reference gas at 272 K on global and hemispheric scales. Thus, the derived hemispheric and global OH estimates are forced to reconcile the observed trends and gradient for all four compounds simultaneously. However, currently, observations of these gases from the surface networks do not provide more accurate OH abundance estimate than that from MCF.
Of the C3Hx hydrocarbons, propane (C3H8) and propyne (methylacetylene, CH3C2H) were first detected in Titan’s atmosphere during the Voyager 1 flyby in 1980. Propene (propylene, C3H6) was first detected in 2013 with data from the Composite InfraRed Spectrometer (CIRS) instrument on Cassini. We present the first measured abundance profiles of propene on Titan from radiative transfer modeling, and compare our measurements to predictions derived from several photochemical models. Near the equator, propene is observed to have a peak abundance of 10 ppbv at a pressure of 0.2 mbar. Several photochemical models predict the amount at this pressure to be in the range 0.3–1 ppbv and also show a local minimum near 0.2 mbar which we do not see in our measurements. We also see that propene follows a different latitudinal trend than the other C3 molecules. While propane and propyne concentrate near the winter pole, transported via a global convective cell, propene is most abundant above the equator. We retrieve vertical abundances profiles between 125 km and 375 km for these gases for latitude averages between 60°S–20°S, 20°S–20°N, and 20°N–60°N over two time periods, 2004 through 2009 representing Titan’s atmosphere before the 2009 equinox, and 2012 through 2015 representing time after the equinox. Additionally, using newly corrected line data, we determined an updated upper limit for allene (propadiene, CH2CCH2, the isomer of propyne). We claim a 3-σ upper limit mixing ratio of 2.5 × 10(exp −9) within 30° of the equator. The measurements we present will further constrain photochemical models by refining reaction rates and the transport of these gases throughout Titan’s atmosphere.
Previous works on Membrane Aerobic Biological Reactors (MABR) CoMANDR 1.0, CoMANDR 2.0, and R-CoMANDAR have demonstrated their ability to stabilize various space based waste streams over operating periods of ~1 year. Biological pretreatment by MABR systems can stabilize space based waste streams. Biological stabilization includes reducing the pH, conversion of organic N to NOx- and oxidation of dissolved organic matter to CO2. These processes produce a more stable waste product (brine), facilitate the distillation processes, and enable evaporative or membrane based systems. An alternative to aerobic operation would be to include anoxic operation to promote denitrification and production of N2 gas. This results in a reduced O2 demand and increases ammonia oxidation efficiency. Denitrifictaion can be accomplished in either a single reactor (Simultaneous Nitrification Denitrification) or in a two-stage system with separate aerobic and anoxic reactors. We evaluated the performance of both architectures in pilot scale systems (1-2 crew/d). Each system was continuously operated for over 2 years during which they processed a variety of habitation waste streams (ISS, Transit, and EPB) in both a continuous and on production feed mode. Here we report the results of the two stage system. Results indicate that the two stage system can successfully remove organic carbon, lower pH and convert organic N to N2 gas. Organic carbon and organic N oxidation reaction rates for the two stage system are similar to past studies for single stage aerobic systems. The two stage system is more complex and requires an additional pump. While no maintenance was required on the system during the nearly two year period of operation, the packed bed did produce N2 gas for many operational test points. The performance and comparison of operational conditions are detailed below