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 343 records · Page 19

Investigation of chemically-reacting supersonic internal flows

This report covers work done on the research project Analysis and Computation of Internal Flow Field in a Scramjet Engine. The work is supported by the NASA Langley Research Center (Computational Methods Branch of the High-Speed Aerodynamics Division) through research grant NAG1-423. The governing equations of two-dimensional chemically-reacting flows are presented together with the global two-step chemistry model. The finite-difference algorithm used is illustrated and the method of circumventing the stiffness is discussed. The computer program developed is used to solve two model problems of a premixed chemically-reacting flow. The results obtained are physically reasonable.

Chitsomboon, T.↗

Re-Entry Aeroheating Analysis of Tile-Repair Augers for the Shuttle Orbiter

Computational re-entry aerothermodynamic analysis of the Space Shuttle Orbiter s tile overlay repair (TOR) sub-assembly is presented. Entry aeroheating analyses are conducted to characterize the aerothermodynamic environment of the TOR and to provide necessary inputs for future TOR thermal and structural analyses. The TOR sub-assembly consists of a thin plate and several augers and spacers that serve as the TOR fasteners. For the computational analysis, the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) is used. A 5-species non-equilibrium chemistry model with a finite rate catalytic recombination model and a radiation equilibrium wall condition are used. It is assumed that wall properties are the same as reaction cured glass (RCG) properties with a surface emissivity of epsilon = 0.89. Surface heat transfer rates for the TOR and tile repair augers (TRA) are computed at a STS-107 trajectory point corresponding to Mach 18 free stream conditions. Computational results show that the average heating bump factor (BF), which is a ratio of local heat transfer rate to a design reference point located at the damage site, for the auger head alone is about 1.9. It is also shown that the average BF for the combined auger and washer heads is about 2.0.

Mazaheri, Ali R.↗

Computation of wake/exhaust mixing downstream of advanced transport aircraft

The mixing of engine exhaust with the vortical wake of high speed aircraft operating in the stratosphere can play an important role in the formation of chemical products that deplete atmospheric ozone. An accurate analysis of this type of interaction is therefore necessary as a part of the assessment of the impact of proposed High Speed Civil Transport (HSCT) designs on atmospheric chemistry. This paper describes modifications to the parabolic Navier-Stokes flow field analysis in the UNIWAKE unified aircraft wake model to accommodate the computation of wake/exhaust mixing and the simulation of reacting flow. The present implementation uses a passive chemistry model in which the reacting species are convected and diffused by the fluid dynamic solution but in which the evolution of the species does not affect the flow field. The resulting analysis, UNIWAKE/PCHEM (Passive CHEMistry) has been applied to the analysis of wake/exhaust flows downstream of representative HSCT configurations. The major elements of the flow field model are described, as are the results of sample calculations illustrating the behavior of the thermal exhaust plume and the production of species important to the modeling of condensation in the wake. Appropriate steps for further development of the UNIWAKE/PCHEM model are also outlined.

Quackenbush, Todd R.↗

Scramjet combustor and nozzle computations

The SHIP3D PNS code is used to conduct a 3D parametric study of a scramjet combustor and nozzle configuration at flight Mach numbers of 8, 10, 12, 16, and 20. The study demonstrated the powerful computational strategy of using a highly efficient, modified PNS code and a nonkinetic chemistry model to analyze test data and conduct parametric studies of scramjet combustors and nozzles. The mixing efficiency was found to decrease by 0.2 when the flight Mach number increased from 8 to 12 and showed little change at higher Mach numbers. The protruding ramp injector generated an oblique shock that impinged on the cowl. Ramp-slot injection produced better mixing, both because of the interaction of this shock with the slot fuel, and because of better fuel distribution in the duct. The former effect was enhanced when the ramp angle was increased. Mixing was also enhanced when the ramp fuel injection angle and sweep angle were increased.

Kamath, Pradeep S.↗

Numerical modeling of wall-injected scramjet experiments

A wall-injected, hydrogen-fueled scramjet is modeled numerically using a parabolic Navier-Stokes computer code with a k-epsilon turbulence model and finite-rate chemistry. The numerical results are compared to experimental scramjet data taken in a shock tunnel and are found to be in reasonable agreement. Under the conditions studied, the numerical results show that combustion in the scramjet is mainly mixing limited.

Brescianini, C. P.↗

Simulation of summertime ozone over North America

The concentrations of O3 and its precursors over North America are simulated for three summer months with a 3D, continental-scale photochemical model using meteorological input from the Goddard Institute for Space Studies (GISS) GCM. The model has 4 x 5 deg grid resolution and represents nonlinear chemistry in urban and industrial plumes with a subgrid nested scheme. Simulated median afternoon O3 concentrations at rural U.S. sites are within 5 ppb of observations in most cases, except in the south central U.S., where concentrations are overpredicted by 15-20 ppb. The model captures successfully the development of regional high-O3 episodes over the northeastern United States on the back side of weak, warm, stagnant anticyclones. Simulated concentrations of CO and nonmethane hydrocarbons are generally in good agreement with observations, concentrations of NO(x) are underpredicted by 10-30 percent, and concentrations of PANs are overpredicted by a factor of 2 to 3. The overprediction of PANs is attributed to flaws in the photochemical mechanism, including excessive production from oxidation of isoprene, and may also reflect an underestimate of PANs deposition. Subgrid nonlinear chemistry as captured by the nested plumes scheme decreases the net O3 production computed in the U.S. boundary layer by 8 percent on average.

Jacob, Daniel J.↗

Turbulent Flame Processes Via Diffusion Flame-Vortex Ring Interactions

Flame-vortex interactions are canonical configurations that can be used to study the underlying processes occurring in turbulent reacting flows. This configuration contains many of the fundamental aspects of the coupling between fluid dynamics and combustion that could be investigated with more controllable conditions than are possible under direct investigations of turbulent flames. Diffusion flame-vortex ring interaction contains many of the fundamental elements of flow, transport, combustion, and soot processes found in turbulent diffusion flames. Some of these elements include concentrated vorticity, entrainment and mixing, strain and nonequilibrium phenomena, diffusion and differential diffusion, partial premixing and diluent effects, soot formation and oxidation, and heat release effects. Such simplified flowfield allows the complex processes to be examined more closely and yet preserving the physical processes present in turbulent reacting flows. Furthermore, experimental results from the study of flame-vortex interactions are useful for the validation of numerical simulations and more importantly to deepen our understanding of the fundamental processes present in reacting flows. Experimental and numerical results obtained under microgravity conditions of the diffusion flame-vortex ring interaction are summarized in this paper. Results are obtained using techniques that include Flame Luminosity Imaging (FLI), Laser Soot-Mie Scattering (LSMS), Computational Fluid Dynamics and Combustion (CFDC), and Diode Laser Spectroscopy/Iterative Temperature with Assumed Chemistry (DLS/ITAC).

Dahm, Werner J. A.↗

Studies of Trace Gas Chemical Cycles Using Inverse Methods and Global Chemical Transport Models

We report progress in the first year, and summarize proposed work for the second year of the three-year dynamical-chemical modeling project devoted to: (a) development, testing, and refining of inverse methods for determining regional and global transient source and sink strengths for long lived gases important in ozone depletion and climate forcing, (b) utilization of inverse methods to determine these source/sink strengths using either MATCH (Model for Atmospheric Transport and Chemistry) which is based on analyzed observed wind fields or back-trajectories computed from these wind fields, (c) determination of global (and perhaps regional) average hydroxyl radical concentrations using inverse methods with multiple titrating gases, and (d) computation of the lifetimes and spatially resolved destruction rates of trace gases using 3D models. Important goals include determination of regional source strengths of methane, nitrous oxide, methyl bromide, and other climatically and chemically important biogenic/anthropogenic trace gases and also of halocarbons restricted by the Montreal protocol and its follow-on agreements and hydrohalocarbons now used as alternatives to the restricted halocarbons.

Prinn, Ronald G.↗

TPSAS-NF1676L-21180-DND

Ablative materials are used as the Thermal Protection System (TPS) for space vehicles that enter a planetary atmosphere. During the entry phase, the vehicle experiences high rates of aerodynamic heating. These ablative materials, acting as a heat shield, protect the vehicle from the harsh environment and ensure a safe entry. Ablative materials reduce the incoming heat flux through various phenomena such as near-surface oxidation, pyrolysis chemical reaction, vaporization, and other erosive processes, collectively known as ablation. Spallation is one such ablative process, in which solid particles are ejected from the material into the flow field. The many conditions leading to spallation, however, are not well understood and cannot be described accurately by simple phenomenological models. However, the spallation phenomenon might impact the aerodynamic heating rates of re-entry vehicles using ablative heat shields. For example, the solid particles ejected from the surface can modify the near surface chemical processes and flow field. To investigate spallation effects, a code was developed to compute the dynamics of spalled particles. The code uses a finite-rate chemistry model to study the chemical interactions of the particles with the flow field. The spallation code is dynamically coupled to a CFD solver that models the time accurate hypersonic flow field around an ablative sample. An experimental investigation was also performed at the NASA Langley HYMETS facility. HYMETS is a 400 kW arcjet that can simulate hypersonic flight and Earth or Martian entry conditions. Using high speed cameras, an IR camera, NIR and UV-VIS spectrometers, and other diagnostic instruments, test data on various graphite, LI-900, FiberForm, and PICA test samples were obtained, using target heat flux conditions of 100, 200, and 400 W/cm 2 . From the test campaign, measurements of particle velocity, acceleration and trajectory are obtained.

A Martin↗

A Look at the Truths and Misconceptions of the Variational Quantum Eigensolver and the Implications of Overparameterization

In this work, we investigate loss landscapes of the variational quantum eigensolver (VQE) by quantifying the number of local minima through empirical analyses. We focus on minimal models in chemistry and physics so that we can do a complete analysis using more computationally expensive tools. We employ Hessian eigenvalue calculations and the nudged elastic band algorithm to characterize these landscapes. Our results expand upon the existing literature by highlighting the optimization challenges faced by VQE. We find that, as the number of parameters in our ansatz increases, the number of basins increases while the corresponding loss function values converge toward the global minimum value. This observation implies that overparameterization may lead to an ``effective convexity'' in VQE loss landscapes, a phenomenon supported by theoretical and numerical work in classical machine learning.

quantum computing↗

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Comparisons are made between potential energy surfaces (PES) for N2 + N and N2 + N2 collisions and between rate coefficients for N2 dissociation that were computed using the quasiclassical trajectory method (QCT) on these PESs. For N2 + N we compare the Laganas empirical LEPS surface with one from NASA Ames Research Center based on ab initio quantum chemistry calculations. For N2 + N2 we compare two ab initio PESs (from NASA Ames and from the University of Minnesota). These use different methods for computing the ground state electronic energy for N4, but give similar results. Thermal N2 dissociation rate coefficients, for the 10,000K-30,000K temperature range, have been computed using each PES and the results are in excellent agreement. Quasi-stationary state (QSS) rate coefficients using both PESs have been computed at these temperatures using the Direct Molecular Simulation of Schwartentruber and coworkers. The QSS rate coefficients are up to a factor of 5 lower than the thermal ones and the thermal and QSS values bracket the results of shock-tube experiments. We conclude that the combination of ab initio quantum chemistry PESs and QCT calculations provides an attractive approach for the determination of accurate high-temperature rate coefficients for use in aerothermodynamics modeling.

quasiclassical trajectory method↗

Numerical simulations of non-equilibrium shock layers with efficient implicit schemes

Current and future calculations of nonequilibrium shock layers require the use of a very large number of equations, due to a multiplicity of chemical species, excited states, and internal energy modes. The computational cost associated with the use of standard implicit methods becomes prohibitive; it is therefore desirable to examine the potential of several methods and determine if any can be projected to be more efficient and accurate for large systems of equations. Here, the performance of several implicit schemes on several simple practical examples of reacting flows is examined. The Euler equations are solved by three different implicit methods, and two methods of coupling between the fluid dynamics and the chemistry are studied. Several cases of stiffness are considered and both 1D and 2D examples are computed.

Cambier, Jean-Luc↗

Numerical solutions of the complete Navier-Stokes equations

The physical phenomena within supersonic flows that sustain chemical reactions are investigated. An earlier study to develop accurate physical models for supersonic reacting flowfields focused on 2-D laminar shear layers. The objective is to examine the mixing and subsequent combustion within turbulent reacting shear layers. To conduct this study, a computer program has been written to solve the axisymmetric Reynolds averaged Navier-Stokes equations. The numerical method uses a cell-centered finite volume approach and a Runge Kutta time stepping scheme. The Reynolds averaged equations are closed using the eddy viscosity concept. Several zero-equation models have been tested by making calculations for an H2-air nonreacting coaxial jet flow. Comparisons made with experimental data show that Cohen's eddy viscosity model provides best agreement. The finite rate chemistry model used in the study of 2-D laminar shear layers is incorporated into the computer program and data is compared from a recent experiment performed at NASA Langley.

Hassan, H. A.↗

Bohm's Quantum Potential and the Visualization of Molecular Structure

David Bohm's ontological interpretation of quantum theory can shed light on otherwise counter-intuitive quantum mechanical phenomena including chemical bonding. In the field of quantum chemistry, Richard Bader has shown that the topology of the Laplacian of the electronic charge density characterizes many features of molecular structure and reactivity. Visual and computational examination suggests that the Laplacian of Bader and the quantum potential of Bohm are morphologically equivalent. It appears that Bohmian mechanics and the quantum potential can make chemistry as clear as they makes physics.

Levit, Creon↗

Overview of High-Fidelity Modeling Activities in the Numerical Propulsion System Simulations (NPSS) Project

A high-fidelity simulation of a commercial turbofan engine has been created as part of the Numerical Propulsion System Simulation Project. The high-fidelity computer simulation utilizes computer models that were developed at NASA Glenn Research Center in cooperation with turbofan engine manufacturers. The average-passage (APNASA) Navier-Stokes based viscous flow computer code is used to simulate the 3D flow in the compressors and turbines of the advanced commercial turbofan engine. The 3D National Combustion Code (NCC) is used to simulate the flow and chemistry in the advanced aircraft combustor. The APNASA turbomachinery code and the NCC combustor code exchange boundary conditions at the interface planes at the combustor inlet and exit. This computer simulation technique can evaluate engine performance at steady operating conditions. The 3D flow models provide detailed knowledge of the airflow within the fan and compressor, the high and low pressure turbines, and the flow and chemistry within the combustor. The models simulate the performance of the engine at operating conditions that include sea level takeoff and the altitude cruise condition.

Veres, Joseph P.↗

A numerical method for predicting hypersonic flowfields

The flow about a body traveling at hypersonic speed is energetic enough to cause the atmospheric gases to react chemically and reach states in thermal nonequilibrium. In this paper, a new procedure based on Gauss-Seidel line relaxation is shown to solve the equations of hypersonic flow fields containing finite reaction rate chemistry and thermal nonequilibrium. The method requires a few hundred time steps and small computer times for axisymmetric flows about simple body shapes. The extension to more complex two-dimensional body geometries appears straightforward.

Maccormack, Robert W.↗

The computation of non-equilibrium chemically-reacting flows

Numerical techniques for the analysis of flows with nonequilibrium thermodynamics and chemistry are developed and demonstrated. Particular attention is given to the formulation of the thermodynamic and chemical gas models, the derivation of the governing equations, and the solution procedures. Results for typical problems involving shock tubes and supersonic diffusers are presented in extensive graphs and briefly characterized.

Grossman, B.↗

Radiative Extinction of Gaseous Spherical Diffusion Flames in Microgravity

Radiative extinction of spherical diffusion flames was investigated experimentally and numerically. The experiments involved microgravity spherical diffusion flames burning ethylene and propane at 0.98 bar. Both normal (fuel flowing into oxidizer) and inverse (oxidizer flowing into fuel) flames were studied, with nitrogen supplied to either the fuel or the oxygen. Flame conditions were chosen to ensure that the flames extinguished within the 2.2 s of available test time; thus extinction occurred during unsteady flame conditions. Diagnostics included color video and thin-filament pyrometry. The computations, which simulated flow from a porous sphere into a quiescent environment, included detailed chemistry, transport and radiation, and yielded transient results. Radiative extinction was observed experimentally and simulated numerically. Extinction time, peak temperature, and radiative loss fraction were found to be independent of flow rate except at very low flow rates. Radiative heat loss was dominated by the combustion products downstream of the flame and was found to scale with flame surface area, not volume. For large transient flames the heat release rate also scaled with surface area and thus the radiative loss fraction was largely independent of flow rate. Peak temperatures at extinction onset were about 1100 K, which is significantly lower than for kinetic extinction. One observation of this work is that while radiative heat losses can drive transient extinction, this is not because radiative losses are increasing with time (flame size) but rather because the heat release rate is falling off as the temperature drops.

Santa, K. J.↗