Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Chemistry”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Modeling of Alkane Oxidation Using Constituents and Species

It is currently not possible to perform simulations of turbulent reactive flows due in particular to complex chemistry, which may contain thousands of reactions and hundreds of species. This complex chemistry results in additional differential equations, making the numerical solution of the equation set computationally prohibitive. Reducing the chemical kinetics mathematical description is one of several important goals in turbulent reactive flow modeling. A chemical kinetics reduction model is proposed for alkane oxidation in air that is based on a parallel methodology to that used in turbulence modeling in the context of the Large Eddy Simulation. The objective of kinetic modeling is to predict the heat release and temperature evolution. This kinetic mechanism is valid over a pressure range from atmospheric to 60 bar, temperatures from 600 K to 2,500 K, and equivalence ratios from 0.125 to 8. This range encompasses diesel, HCCI, and gas-turbine engines, including cold ignition. A computationally efficient kinetic reduction has been proposed for alkanes that has been illustrated for n-heptane using the LLNL heptane mechanism. This model is consistent with turbulence modeling in that scales were first categorized into either those modeled or those computed as progress variables. Species were identified as being either light or heavy. The heavy species were decomposed into defined 13 constituents, and their total molar density was shown to evolve in a quasi-steady manner. The light species behave either in a quasi-steady or unsteady manner. The modeled scales are the total constituent molar density, Nc, and the molar density of the quasi-steady light species. The progress variables are the total constituent molar density rate evolution and the molar densities of the unsteady light species. The unsteady equations for the light species contain contributions of the type gain/loss rates from the heavy species that are modeled consistent with the developed mathematical forms for the total constituent molar density rate evolution; indeed, examination of these gain/loss rates shows that they also have a good quasi-steady behavior with a functional form resembling that of the constituent rate. This finding highlights the fact that the fitting technique provides a methodology that can be repeatedly used to obtain an accurate representation of full or skeletal kinetic models. Assuming success with the modified reduced model, the advantage of the modeling approach is clear. Because this model is based on the Nc rate rather than on that of individual heavy species, even if the number of species increases with increased carbon number in the alkane group, providing that the quasi-steady rate aspect persists, then extension of this model to higher alkanes should be conceptually straightforward, although it remains to be seen if the functional fits would remain valid or would require reconstruction.

Bellan, Jasette↗
NASA NTRS2014

Model Analysis of the Factors Regulating Trends and Variability of Methane, Carbon Monoxide and OH: 1. Model Validation

Methane (CH4) is the second most important anthropogenic greenhouse gas (GHG). Its 100-year global warming potential (GWP) is 25 times larger than that for carbon dioxide. The 100-yr integrated GWP of CH4 is sensitive to changes in OH levels. Methane's atmospheric growth rate was estimated to be more than 10 ppb yr(exp -1) in 1998 but less than zero in 2001, 2004 and 2005 (Kirschke et al., 2013). Since 2006, the CH4 is increasing again. This phenomena is yet not well understood. Oxidation of CH4 by OH is the main loss process, thus affecting the oxidizing capacity of the atmosphere and contributing to the global ozone background. Current models typically use an annual cycle of offline OH fields to simulate CH4. The implemented OH fields in these models are typically tuned so that simulated CH4 growth rates match that measured. For future and climate simulations, the OH tuning technique may not be suitable. In addition, running full chemistry, multi-decadal CH4 simulations is a serious challenge and currently, due to computational intensity, almost impossible.

OH↗
NASA NTRSDate not supplied

Modeling Entry Systems to Explore Our Solar System

Exploration of our Solar System is a foundational element of NASA’s identity. Delivering a scientific payload through an atmosphere to the surface of a planetary body requires safely navigating the extreme temperatures and stresses generated by flying many times the speed of sound. An entry system is the outer shell of a vehicle designed to protect the payload and, as a single point of failure for a mission, reliability is paramount. Unfortunately, it is not possible to fully replicate the flight environment in ground test facilities – so how do we confidently design a vehicle that needs to work the first time, every time? Modeling and simulation are critical tools for filling gaps in ground test capability and providing traceability from ground to flight. Modeling an entry system is a truly interdisciplinary effort, requiring knowledge of fluid dynamics, high-temperature chemistry, radiation, materials science, structural dynamics, guidance and control – and, finally, the mathematical and computing capability to pull it all together. In this talk, Dr. Michael Barnhardt will discuss the latest research in modeling entry systems and how it is being used in NASA’s exploration missions.

Entry systems↗
NASA NTRS2021

Petrology, Geochemistry, and Pairing of Lunar Meteorites from the Dominion Range

Introduction: During the 2018-2019 Antarctic Search for Meteorites (ANSMET) field season in the Dominion Range (DOM), 7 lunar meteorite stones were collected: DOM 18242 (15.1 g), DOM 18244 (25.1 g), DOM 18262 (6.8 g), DOM 18509 (16.5 g), DOM 18543 (13.6 g), DOM 18666 (45.9 g), and DOM 18678 (11.6 g). Here we present the initial results of electron microprobe and X-ray computed tomography (XCT) studies of these stones and look at the details of their petrography and mineral chemistry, as well as investigate possible pairing relationships, both with each other and with previously described lunar meteorites. Most of the work presented here is on the DOM 18509, 18543, and 18678 stones; subsamples of the other stones are in hand and similar measurements will be made on them by the time of the meeting. Methods: Textures in DOM 18509, DOM 18543, and DOM 18678 were characterized in 2D by optical microscopy, backscattered electron (BSE) imaging and elemental X-ray images on thin sections, as well as in 3D by X-ray computed tomography (XCT) on sample chips. Mineral compositions were assessed through a combination of wavelength dispersive spectroscopy EPMA (electron probe microanalysis) and x-ray mapping on the JEOL 8530 at NASA JSC. The bulk composition of all three meteorites was determined based on analyses of the fusion crust glass. XCT analyses were done on the Nikon XTH 320 at NASA JSC. ICP-MS data on bulk rock subsamples for each meteorite will be carried out in the near future. Results: The stones are all similar in macroscopic appearance with a dark aphanitic matrix hosting a variety of small- to medium-sized angular mineral and lithic fragments (often light colored in nature) [1,2]. Based on EPMA and XCT results, the three meteorites are polymict regolith breccias comprised of mineral, glass, and lithic clasts ranging up to several mm in length. Melt veins run through all three meteorite samples. Mineral clasts in all 3 stones are dominated by pyroxene and plagioclase (An82-96), with minor amounts of SiO2, olivine (Fo1-52), and FeTi-oxides. Pyroxene grains are mostly Fe-rich pigeonite and augite, and larger clasts are normally zoned and have fine exsolution lamellae. The lithic clasts in all stones consists of: (1) basalt clasts that contain zoned pyroxene, plagioclase laths, and ilmenite, with minor silica and Fe-rich olivine; (2) granulitic clasts; (3) anorthosite clasts; (4) Si-rich clasts that also contain ilmenite, troilite, high-Ca pyroxene, fayalite, and K-feldspar likely mesostasis from late stage basalts). All three meteorites contain glassy fusion crust that is highly vesicular, high in FeO and Al2O3 (15-17 wt% each), ferroan (Mg# of 23-24), and moderately rich in TiO2 (1.4-1.8 wt%). The composition of the fusion crust can serve as a proxy for the bulk meteorite composition and is identical within error for all three meteorites. Implications: The lithic and mineral clasts in all three stones are similar in clast population and assemblages as well as mineral chemistry. In addition, the fusion crust composition, a proxy for bulk composition, is within error of each other for all three stones. Thus DOM 18509, DOM 18543, and DOM 18678 are almost certainly paired. Based on similarities in macroscopic description as well as preliminary classification data [1,2], all 7 lunar stones from DOM are likely paired, though additional quantitative analyses are needed to confirm this. The presence of spherules and vesicular fusion crusts indicates that the DOM pairing group is a regolith breccia. The presence of basaltic and gabbroic clasts as well as more feldspathic materials, suggest that the regolith from which these meteorites formed contained a mixture of feldspathic highland material and mare material, suggesting a possible provenance near a marehighlands boundary. No evidence of KREEPy lithologies have been observed so far in these meteorites, however, future ICP-MS data on bulk rock chips for the stones will reveal any KREEP component if present. The DOM pairing group has many similarities to previously described lunar breccia meteorite MET 01210, however more detailed compositional data will be needed to make a definitive comparison.

R.A. Zeigler↗
NASA NTRSDate not supplied

3D Material Response of the MSL Heatshield Using NuSil-Coated PICA

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield that used NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performance for protecting probes during planetary entry. Data from the Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. A space-grade silicone-based coating commercially known as NuSil CV-1144-0 [3] was applied to the entire MSL heatshield, including the MEDLI plugs, to mitigate the spread of dust from PICA. Modeling the thermal response of PICA-NuSil (PICA-N) system is still an open challenge. Ground testing of PICA-N models exhibited surface temperature jumps of the order of 150 K due to oxide scale formation and sub-sequent NuSil burn-off. It is therefore critical to include a validated model for the material response of the coating in engineering codes. A test campaign has been conducted at the NASA’s Langley HyMETS [4] facility to screen the response of PICA-N and gather detailed data on its behavior [5]. A first model of PICA-N thermal response has been developed using the Hy-METS experiments [6]. The objective of this work is to analyze the material response of the latest PICA-N model compared to the engineering model used to simulate the entry of MSL. The environment and material response around the MSL aeroshell during Mars atmospheric entry is simulated using a collection of tools. The Direct Simulation Monte Carlo SPARTA code [7] is used in the rarefied regime, the Data Parallel Line Relaxation (DPLR) code [8] is used in the continuum regime and radiative heating conditions are provided by the Nonequilibrium air radiation (NEQAIR) code [9] to estimate the environmental conditions. The thermal response inside the material is computed using the Porous material Analysis Toolbox based on Open-FOAM (PATO) [10,11,12]. Thermodynamic and chemistry properties are estimated using the Mutation++ library [13]. The approach implemented in PATO as a first cut PICA-N thermal response model is outlined in Figure 1. While the recession is less than the coating thickness, the Surface mass and energy balance Boundary Condition (SBC) uses the NuSil B’ tables. Once the recession removes the coating, the usual PICA B’ tables are used for the SBC. The B’ tables are computed using an equilibrium solver implemented in Mutation++, given the temperature, pressure, blowing rate, composition of the pyrolysis and environment gases, and the condensed species at the surface. Preliminary results of the 3D material response of the MSL heat-shield at the peak heating (80 sec after Entry Interface) are shown in Figure 2. Current NASA’s mission to Mars, Mars 2020, used the spare heatshield of MSL for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the PATO high-fidelity material response capability was benchmarked against flight data from MEDLI. This effort represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes. This bench-marking is awaiting the final release of the MEDLI-2 data.

Aerospace↗
NASA NTRS2018

Transition Metal Oxides as Cathodes in Li-O2 battery: A First Principles Investigation

Li-O2 batteries have traditionally used carbon based electrodes (graphite, buckypaper) as the cathode of choice due to its good electrical conductivity, stability against non-aqueous electrolytes like Dimethyl ether (DME) and ease of handling. But, the carbon cathode also leads to formation of carbonate by-products that increase overpotentials during charging leading to degradation of cathode and reduction of cyclability. In this work, we investigate some of the well-known oxides as cathodes with focus on the interface between the oxide surfaces and the discharge product: Li2O2, in the Li-O2 battery using first principles computations. Our results show that attention must be paid on choosing the appropriate surface of the oxides. We extend the analysis to suggest other possible oxide chemistries that should be investigated as cathodes in Li-O2 batteries.

Li-O2 battery↗
NASA NTRS2005

Incorporation of a Chemical Equilibrium Equation of State into LOCI-Chem

Renewed interest in development of advanced high-speed transport, reentry vehicles and propulsion systems has led to a resurgence of research into high speed aerodynamics. As this flow regime is typically dominated by hot reacting gaseous flow, efficient models for the characteristic chemical activity are necessary for accurate and cost effective analysis and design of aerodynamic vehicles that transit this regime. The LOCI-Chem code recently developed by Ed Luke at Mississippi State University for NASA/MSFC and used by NASA/MSFC and SSC represents an important step in providing an accurate, efficient computational tool for the simulation of reacting flows through the use of finite-rate kinetics [3]. Finite rate chemistry however, requires the solution of an additional N-1 species mass conservation equations with source terms involving reaction kinetics that are not fully understood. In the equilibrium limit, where the reaction rates approach infinity, these equations become very stiff. Through the use of the assumption of local chemical equilibrium the set of governing equations is reduced back to the usual gas dynamic equations, and thus requires less computation, while still allowing for the inclusion of reacting flow phenomenology. The incorporation of a chemical equilibrium equation of state module into the LOCI-Chem code was the primary objective of the current research. The major goals of the project were: (1) the development of a chemical equilibrium composition solver, and (2) the incorporation of chemical equilibrium solver into LOCI-Chem. Due to time and resource constraints, code optimization was not considered unless it was important to the proper functioning of the code.

Cox, Carey F.↗
NASA NTRS1979

The auroral ionosphere - Comparison of a time-dependent model with composition measurements

A time-dependent model of the auroral ionosphere including the odd nitrogen species, NO, N(D-2), and N(S-4), is used for comparison with data from a coordinated rocket-satellite measurement of an auroral event. The chemical scheme and the adopted rate coefficients have been shown to be compatible with daytime mid-latitude ionospheric chemistry. The electron flux and neutral atmospheric parameters measured on the satellite are used to compute the appropriate ionization and dissociation rates. The calculated NO(plus), O2(plus), O(plus), Ne, and NO densities agree well with the rocket measurements. The calculated N2(plus) densities are larger than the measured densities by a factor of 3 at most altitudes. The calculations show that the nitric oxide content of the aurora (about 1.2 times 10 to the 9th NO molecules/cu cm at 105 km) is below the saturation value.

Gerard, J.-C.↗
NASA NTRS1986

Computer modeling of photodegradation

A computer program to simulate the photodegradation of materials exposed to terrestrial weathering environments is being developed. Input parameters would include the solar spectrum, the daily levels and variations of temperature and relative humidity, and materials such as EVA. A brief description of the program, its operating principles, and how it works was initially described. After that, the presentation focuses on the recent work of simulating aging in a normal, terrestrial day-night cycle. This is significant, as almost all accelerated aging schemes maintain a constant light illumination without a dark cycle, and this may be a critical factor not included in acceleration aging schemes. For outdoor aging, the computer model is indicating that the night dark cycle has a dramatic influence on the chemistry of photothermal degradation, and hints that a dark cycle may be needed in an accelerated aging scheme.

Guillet, J.↗
NASA NTRS1988

Cross-B convection of artificially created, negative-ion clouds and plasma depressions - Low-speed flow

A negative-ion, positive-ion plasma produced by the release of an electron attachment chemical into the F region becomes electrically polarized by collisions with neutrals moving across magnetic field lines. The resulting electric field causes E x B drift of the two ion species and the residual electrons. The cross-field flow of the modified ionosphere is computed using a two-dimensional numerical simulation which includes electron attachment and mutual neutralization chemistry, self-consistent electric fields, and three-species plasma transport. The velocity of the plasma is initially in the direction of the neutral wind because the negative-ion cloud is a Pedersen conductivity enhancement. As the positive and negative ions react, the Pedersen conductivity becomes depressed below the ambient value and the velocity of the plasma reverses direction. A plasma hole remains after the positive and negative ions have mutually neutralized. The E x B gradient drift instability produces irregularities on the upwind edge of the hole.

Bernhardt, Paul A.↗
NASA NTRS1996

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗
NASA NTRS1999

Development of a Computational Chemical Vapor Deposition Model: Applications to Indium Nitride and Dicyanovinylaniline

A computational chemical vapor deposition (CVD) model is presented, that couples chemical reaction mechanisms with fluid dynamic simulations for vapor deposition experiments. The chemical properties of the systems under investigation are evaluated using quantum, molecular and statistical mechanics models. The fluid dynamic computations are performed using the CFD-ACE program, which can simulate multispecies transport, heat and mass transfer, gas phase chemistry, chemistry of adsorbed species, pulsed reactant flow and variable gravity conditions. Two experimental setups are being studied, in order to fabricate films of: (a) indium nitride (InN) from the gas or surface phase reaction of trimethylindium and ammonia; and (b) 4-(1,1)dicyanovinyl-dimethylaminoaniline (DCVA) by vapor deposition. Modeling of these setups requires knowledge of three groups of properties: thermodynamic properties (heat capacity), transport properties (diffusion, viscosity, and thermal conductivity), and kinetic properties (rate constants for all possible elementary chemical reactions). These properties are evaluated using computational methods whenever experimental data is not available for the species or for the elementary reactions. The chemical vapor deposition model is applied to InN and DCVA. Several possible InN mechanisms are proposed and analyzed. The CVD model simulations of InN show that the deposition rate of InN is more efficient when pulsing chemistry is used under conditions of high pressure and microgravity. An analysis of the chemical properties of DCVA show that DCVA dimers may form under certain conditions of physical vapor transport. CVD simulations of the DCVA system suggest that deposition of the DCVA dimer may play a small role in the film and crystal growth processes.

Cardelino, Carlos↗
NASA NTRS1989

Supercomputer requirements for selected disciplines important to aerospace

Speed and memory requirements placed on supercomputers by five different disciplines important to aerospace are discussed and compared with the capabilities of various existing computers and those projected to be available before the end of this century. The disciplines chosen for consideration are turbulence physics, aerodynamics, aerothermodynamics, chemistry, and human vision modeling. Example results for problems illustrative of those currently being solved in each of the disciplines are presented and discussed. Limitations imposed on physical modeling and geometrical complexity by the need to obtain solutions in practical amounts of time are identified. Computational challenges for the future, for which either some or all of the current limitations are removed, are described. Meeting some of the challenges will require computer speeds in excess of exaflop/s (10 to the 18th flop/s) and memories in excess of petawords (10 to the 15th words).

Peterson, Victor L.↗
NASA NTRSDate not supplied

Aerothermal Analysis of the Dragonfly Titan Entry

The Dragonfly mission will send a rotorcraft lander to the surface of Saturn’s moon Titan as part of the New-Frontiers program. This will be the first spacecraft to land on Titan since the Huygens probe’s descent in January of 2005 and only the second spacecraft to enter Titan’s atmosphere. The Dragonfly entry capsule is significantly larger than the Huygens probe and will experience higher aerothermal environments. This poster will provide an overview of the aerothermodynamic models used in the Dragonfly aeroshell design process. Titan provides a unique entry environment that has several fundamental differences from the environments at the more common entry destinations of Earth and Mars. The lack of atmospheric oxygen significantly reduces heatshield recession (expected to be negligible), which simplifies material response analysis but also reduces the efficiency of the thermal protection system. The atmos-phere is composed of mostly nitrogen and trace amounts of methane. During hypersonic entry this methane dissociates and leads to the formation of the molecule CN, which is known to radiate strongly in the shock layer. Current Dragonfly heating predictions estimate that radiation contributes about 50% of the total heat flux along the forebody and up to 90% on the aftbody, making radiation a key element of aerothermal analysis. To properly address the importance of radiation, a unique radiative-heating correlation was developed that incorporates normal-shock equilibrium chemistry predictions based on work from the Mars 2020 program. These new heating correlations agree well with Computational Fluid Dynamics (CFD) predictions and enable large-scale trajectory analyses by providing rapid heating predictions along candidate trajectories. Dragonfly aerothermal environments are generated with the Data Parallel Line Relaxation (DPLR) finite-volume Navier-Stokes solver alongside the Nonequilibrium Radiative Transport and Spectra (NEQAIR) radiation transport code. DPLR produces shock-aligned flowfields that incorporate the effects of both chemical and thermal nonequilibrium while NEQAIR solves the radiation transport equations along lines-of-sight through the flowfield solutions to predict the radiative heating on the aeroshell. Titan’s atmosphere produces complex thermal and chemical nonequilibria and so classical radiation approximations that extrapolate from a single line-of-sight, such as the Tangent-Slab method, were found to overpredict peak radiative heat fluxes by 10-20%. Therefore, all radiative heating is directly calculated using 3-D transport for both the forebody and aftbody to reduce unnecessary conservatism. A 21-species finite-rate chemistry model, including electrons and ionized species, is used in the CFD to ensure accurate vehicle heating predictions (the level of ionization is particularly for accurate radiative flux predictions). This poster will cover: 1) The models and assumptions used to simulate entry into Titan’s atmosphere, 2) Dragonfly aerothermal engineering correlations alongside statistical heating values, and 3) Driving aerothermal considerations for a Titan entry such as CN formation and radiation-flowfield coupling.

EDL↗
NASA NTRS1988

Numerical simulation of hypersonic inlet flows with equilibrium or finite rate chemistry

An efficient numerical program incorporated with comprehensive high temperature gas property models has been developed to simulate hypersonic inlet flows. The computer program employs an implicit lower-upper time marching scheme to solve the two-dimensional Navier-Stokes equations with variable thermodynamic and transport properties. Both finite-rate and local-equilibrium approaches are adopted in the chemical reaction model for dissociation and ionization of the inlet air. In the finite rate approach, eleven species equations coupled with fluid dynamic equations are solved simultaneously. In the local-equilibrium approach, instead of solving species equations, an efficient chemical equilibrium package has been developed and incorporated into the flow code to obtain chemical compositions directly. Gas properties for the reaction products species are calculated by methods of statistical mechanics and fit to a polynomial form for C(p). In the present study, since the chemical reaction time is comparable to the flow residence time, the local-equilibrium model underpredicts the temperature in the shock layer. Significant differences of predicted chemical compositions in shock layer between finite rate and local-equilibrium approaches have been observed.

Yu, Sheng-Tao↗
NASA NTRS1976

Technique for atmospheric rate chemistry calculations

The possibility that predictions of atmospheric photochemistry/transport models are sensitive to uncertainties in reaction rates and other inputs stresses the need for rapid numerical integration schemes in rate photochemistry problems. Reducing the computational burden has a major merit in facilitating sensitivity studies to assess the effect of uncertainties on predicted ozone diminutions from NOx (NO + NO2) in the exhaust plume of SST engines. The paper discusses the validity of an algorithmic approach to integration of rate chemistry problems in combustion, developed by Rubel and Baronti for an approximate calculation of the production rate of the i-th chemical species involved. An analysis of two projected SST engines confirms the validity of the proposed algorithm. Because of the relative arithmetical simplicity, it may be easier to treat diffusion rate chemistry calculations using the Rubel and Baronti approximation than would be possible by other approaches.

Matloff, G. L.↗
NASA NTRS2013

Photosynthesis-dependent Isoprene Emission from Leaf to Planet in a Global Carbon-chemistry-climate Model

We describe the implementation of a biochemical model of isoprene emission that depends on the electron requirement for isoprene synthesis into the FarquharBallBerry leaf model of photosynthesis and stomatal conductance that is embedded within a global chemistry-climate simulation framework. The isoprene production is calculated as a function of electron transport-limited photosynthesis, intercellular and atmospheric carbon dioxide concentration, and canopy temperature. The vegetation biophysics module computes the photosynthetic uptake of carbon dioxide coupled with the transpiration of water vapor and the isoprene emission rate at the 30 min physical integration time step of the global chemistry-climate model. In the model, the rate of carbon assimilation provides the dominant control on isoprene emission variability over canopy temperature. A control simulation representative of the present-day climatic state that uses 8 plant functional types (PFTs), prescribed phenology and generic PFT-specific isoprene emission potentials (fraction of electrons available for isoprene synthesis) reproduces 50 of the variability across different ecosystems and seasons in a global database of 28 measured campaign-average fluxes. Compared to time-varying isoprene flux measurements at 9 select sites, the model authentically captures the observed variability in the 30 min average diurnal cycle (R2 6496) and simulates the flux magnitude to within a factor of 2. The control run yields a global isoprene source strength of 451 TgC yr1 that increases by 30 in the artificial absence of plant water stress and by 55 for potential natural vegetation.

photosynthesis↗
NASA NTRSDate not supplied

Application of a Detached Eddy Simulation Approach with Finite-Rate Chemistry to Mars-Relevant Retropropulsion Operating Environments

Human-scale Mars vehicles will require retropropulsion for descent and landing, replacing heritage supersonic parachute systems with an extended phase of powered flight. Due to the limitations of terrestrial testing in Mars-relevant conditions, design and analysis will increasingly rely on computational modeling and simulation. This paper provides an overview of a computational campaign investigating the aerodynamics of a Mars lander concept along various points on a powered descent trajectory including supersonic, transonic, and subsonic conditions using finite-rate chemistry. Simulations using unstructured grids containing billions of elements are performed at scale using thousands of Graphics Processing Units, enabling run-times of a few days for each simulation presented. At each freestream condition, significant minor species concentrations are observed external to the nozzles in the large mixing region upstream of the vehicle. While the flowfields are highly non-stationary in all cases, the mean integrated forces and moments on the vehicle remain small in comparison to the deceleration provided through retropropulsion.

Ashley Korzun↗