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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Prediction of the Thrust Performance and the Flowfield of Liquid Rocket Engines

In an effort to improve the current solutions in the design and analysis of liquid propulsive engines, a computational fluid dynamics (CFD) model capable of calculating the reacting flows from the combustion chamber, through the nozzle to the external plume, was developed. The Space Shuttle Main Engine (SSME) fired at sea level, was investigated as a sample case. The CFD model, FDNS, is a pressure based, non-staggered grid, viscous/inviscid, ideal gas/real gas, reactive code. An adaptive upwinding differencing scheme is employed for the spatial discretization. The upwind scheme is based on fourth order central differencing with fourth order damping for smooth regions, and second order central differencing with second order damping for shock capturing. It is equipped with a CHMQGM equilibrium chemistry algorithm and a PARASOL finite rate chemistry algorithm using the point implicit method. The computed flow results and performance compared well with those of other standard codes and engine hot fire test data. In addition, the transient nozzle flowfield calculation was also performed to demonstrate the ability of FDNS in capturing the flow separation during the startup process.

Wang, T.-S.↗
NASA NTRS1996

Aerothermodynamic Analysis of Commercial Experiment Transporter (COMET) Reentry Capsule

An aerothermodynamic analysis of the Commercial Experiment Transporter (COMET) reentry capsule has been performed using the laminar thin-layer Navier-Stokes solver Langley Aerothermodynamic Upwind Relaxation Algorithm. Flowfield solutions were obtained at Mach numbers 1.5, 2, 5, 10, 15, 20, 25, and 27.5. Axisymmetric and 5, 10, and 20 degree angles of attack were considered across the Mach-number range, with the Mach 25 conditions taken to 90 degrees angle of attack and the Mach 27.5 cases taken to 60 degrees angle of attack. Detailed surface heat-transfer rates were computed at Mach 20 and 25, revealing that heating rates on the heat-shield shoulder ,can exceed the stagnation-point heating by 230 percent. Finite-rate chemistry solutions were performed above Mach 10, otherwise perfect gas computations were made. Drag, lift, and pitching moment coefficients are computed and details of a wake flow are presented. The effect of including the wake in the solution domain was investigated and base pressure corrections to forebody drag coefficients were numerically determined for the lower Mach numbers. Pitching moment comparisons are made with direct simulation Monte Carlo results in the more rarefied flow at the highest Mach numbers, showing agreement within two-percent. Thin-layer Navier-Stokes computations of the axial force are found to be 15 percent higher across the speed range than the empirical/Newtonian based results used during the initial trajectory analyses.

Wood, William A.↗
NASA NTRS1996

A-priori testing of sub-grid models for chemically reacting nonpremixed turbulent shear flows

The beta-assumed-pdf approximation of (Cook & Riley 1994) is tested as a subgrid model for the LES computation of nonpremixed turbulent reacting flows, in the limit of cold infinitely fast chemistry, for two plane turbulent mixing layers with different degrees of intermittency. Excellent results are obtained for the computation of integrals properties such as product mass fraction, and the model is applied to other quantities such as powers of the temperature and the pdf of the scalar itself. Even in these cases the errors are small enough to be useful in practical applications. The analysis is extended to slightly out of equilibrium problems such as the generation of radicals, and formulated in terms of the pdf of the scalar gradients. It is shown that the conditional gradient distribution is universal in a wide range of cases whose limits are established. Within those limits, engineering approximations to the radical concentration are also possible. It is argued that the experiments in this paper are essentially in the limit of infinite Reynolds number.

Jimenez, J.↗
NASA NTRS2005

Applications of Quantum Chemistry to the Study of Carbon Nanotubes

For several years, scientists at NASA Ames have been studying the properties of carbon nanotubes using various experimental and computational methods. In this talk, I will compare different strategies for using quantum chemistry calculations to describe the electronic structure, deformation and chemical functionalization of single wall carbon nanotubes (SWNT) and the physisorption of small molecules on nanotube surfaces. The SWNT can be treated as an infinite (periodic) or finite length carbon cylinder or as a polycyclic aromatic hydrocarbon (PAH) molecule with an imposed curvature maintained by external constraints (as if it were cut out of the SWNT surface). Calculations are carried out using DFT and MP2 methods and a variety of atomic orbital basis sets from minimal (STO-3G) to valence triple zeta. The optimal approach is based on the particular SWNT property of interest. Examples to be discussed include: nanotube fluorination and other functionalization reactions; coating of nanotubes by water vapor and low-molecular weight organic molecules; and the nature of the interface between SWNT and liquids such as water and amines. In many cases, the quantum chemistry calculations are used to parameterize or validate force fields for molecular dynamics simulations. The results of these calculations have helped explain experimental data and contributed to the design of novel materials and sensors based on carbon nanotubes. Some of this research is described in the following papers:

Jaffe, Richard L.↗
NASA NTRS2013

Multi-model Mean Nitrogen and Sulfur Deposition from the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP): Evaluation Historical and Projected Changes

We present multi-model global datasets of nitrogen and sulfate deposition covering time periods from 1850 to 2100, calculated within the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP). The computed deposition fluxes are compared to surface wet deposition and ice-core measurements. We use a new dataset of wet deposition for 2000-2002 based on critical assessment of the quality of existing regional network data. We show that for present-day (year 2000 ACCMIP time-slice), the ACCMIP results perform similarly to previously published multi-model assessments. For this time slice, we find a multi-model mean deposition of 50 Tg(N) yr1 from nitrogen oxide emissions, 60 Tg(N) yr1 from ammonia emissions, and 83 Tg(S) yr1 from sulfur emissions. The analysis of changes between 1980 and 2000 indicates significant differences between model and measurements over the United States but less so over Europe. This difference points towards misrepresentation of 1980 NH3 emissions over North America. Based on ice-core records, the 1850 deposition fluxes agree well with Greenland ice cores but the change between 1850 and 2000 seems to be overestimated in the Northern Hemisphere for both nitrogen and sulfur species. Using the Representative Concentration Pathways to define the projected climate and atmospheric chemistry related emissions and concentrations, we find large regional nitrogen deposition increases in 2100 in Latin America, Africa and parts of Asia under some of the scenarios considered. Increases in South Asia are especially large, and are seen in all scenarios, with 2100 values more than double 2000 in some scenarios and reaching 1300 mg(N) m2 yr1 averaged over regional to continental scale regions in RCP 2.6 and 8.5, 3050 larger than the values in any region currently (2000). The new ACCMIP deposition dataset provides novel, consistent and evaluated global gridded deposition fields for use in a wide range of climate and ecological studies.

Lamarque, J.-F.↗
NASA NTRS2013

Multi-Model Mean Nitrogen and Sulfur Deposition from the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP): Evaluation of Historical and Projected Future Changes

We present multi-model global datasets of nitrogen and sulfate deposition covering time periods from 1850 to 2100, calculated within the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP). The computed deposition fluxes are compared to surface wet deposition and ice core measurements.We use a new dataset of wet deposition for 2000-2002 based on critical assessment of the quality of existing regional network data. We show that for present day (year 2000 ACCMIP time slice), the ACCMIP results perform similarly to previously published multi-model assessments. For this time slice, we find a multimodel mean deposition of approximately 50 Tg(N) yr−1 from nitrogen oxide emissions, 60 Tg(N) yr−1 from ammonia emissions, and 83 Tg(S) yr−1 from sulfur emissions. The analysis of changes between 1980 and 2000 indicates significant differences between model and measurements over the United States but less so over Europe. This difference points towards a potential misrepresentation of 1980 NH3 emissions over North America. Based on ice core records, the 1850 deposition fluxes agree well with Greenland ice cores, but the change between 1850 and 2000 seems to be overestimated in the Northern Hemisphere for both nitrogen and sulfur species. Using the Representative Concentration Pathways (RCPs) to define the projected climate and atmospheric chemistry related emissions and concentrations, we find large regional nitrogen deposition increases in 2100 in Latin America, Africa and parts of Asia under some of the scenarios considered. Increases in South Asia are especially large, and are seen in all scenarios, with 2100 values more than double their 2000 counterpart in some scenarios and reaching >1300 mg(N)m−2 yr−1 averaged over regional to continental-scale regions in RCP 2.6 and 8.5, 30-50% larger than the values in any region currently (circa 2000). However, sulfur deposition rates in 2100 are in all regions lower than in 2000 in all the RCPs. The new ACCMIP multimodel deposition dataset provides state-of-the-science, consistent and evaluated time slice (spanning 1850-2100) global gridded deposition fields for use in a wide range of climate and ecological studies.

Lamarque, J.-F.↗
NASA NTRS2012

Effect of Non-Equilibrium Surface Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using non-equilibrium gas/surface interaction model provides time-accurate solutions for the multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and AblatioN Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas mass conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between the material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of nonequilibrium gas/surface interaction chemistry between air and the carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kanq↗
NASA NTRS2001

Chemistry-Climate Models of the Stratosphere

Over the last decade, improved computer power has allowed three-dimensional models of the stratosphere to be developed that can be used to simulate polar ozone levels over long periods. This paper compares the meteorology between these models, and discusses the future of polar ozone levels over the next 50 years.

Austin, J.↗
NASA NTRS1987

Cometary MHD and chemistry

An MHD and chemical comet-coma model was developed, applying the computer program of Huebner (1985) for the detailed chemical evolution of a spherically expanding coma and the program of Schmidt and Wegman (1982) and Wegman (1987) for the MHD flow of plasma and magnetic field in a comet to the Giotto-mission data on the ion abundances measured by the HIS ion mass spectrometer. The physics and chemistry of the coma are modeled in great detail, including photoprocesses, gas-phase chemical kinetics, energy balance with a separate electron temperature, multifluid hydrodynamics with a transition to free molecular flow, fast-streaming atomic and molecular hydrogen, counter and cross streaming of the ionized species relative to the neutral species in the coma-solar wind interaction region with momentum exchange by elastic collisions, mass-loading through ion pick-up, and Lorentz forces of the advected magnetic field. The results, both inside and outside of the contact surface, are discussed and compared with the relevant HIS ion mass spectra.

Wegmann, R.↗
NASA NTRS1994

Navier-Stokes simulations of Orbiter aerodynamic characteristics including pitch trim and bodyflap

An analysis of the longitudinal aerodynamics of the shuttle orbiter in the hypersonic flight regime is made through the use of computational fluid dynamics. Particular attention is given to establishing the cause of the 'pitching moment anomaly,' which occurred on the orbiter's first flight, and to computing the aerodynamics of a complete orbiter configuration at flight conditions. Data from ground-based facilities as well as orbiter flight data are used to validate the computed results. Analysis shows that the pitching moment anomaly is a real-gas chemistry effect that was not simulated in ground-based facilities, which used air as a test gas. Computed flight aerodynamics for the orbiter are within 5% of the measured flight values and trim bodyflap deflections are predicted to within 10%.

Weilmuenster, K. James↗
NASA NTRS1976

The production of nitric oxide in the troposphere as a result of solid-rocket-motor afterburning

As part of an ongoing assessment of the environmental effects of solid-rocket-motor operations in the troposphere, estimates were made of the nitric oxide produced in the troposphere by the space shuttle and Titan 3-C boosters. Calculations were made with the low-altitude plume computer program and included the effects of coupled finite-rate chemistry and turbulent mixing. A recent measurement of nitric oxide taken in the effluent cloud of a Titan 3-C booster is compared with calculations made with this computer code. The various chemical reactions of the exhaust gases are listed in tabular form.

Stewart, R. B.↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗
NASA NTRS1977

Silicon material task. Part 3: Low-cost silicon solar array project

The feasibility of a process for carbon reduction of low impurity silica in a plasma heat source was investigated to produce low-cost solar-grade silicon. Theoretical aspects of the reaction chemistry were studied with the aid of a computer program using iterative free energy minimization. These calculations indicate a threshold temperature exists at 2400 K below which no silicon is formed. The computer simulation technique of molecular dynamics was used to study the quenching of product species.

Roques, R. A.↗
NASA NTRS1992

Analysis of a hydrocarbon scramjet with augmented preburning

This paper presents the results of a feasibility study of a hydrocarbon scramjet design utilizing an augmented preburner upstream of the main fuel injector locations. The combustor design evaluated here is for a small hypersonic research vehicle. It consists of a preburner into which a small amount of fuel is burned with on-board liquid oxygen and injected into the main airflow, upstream of the main fuel injector locations, thus ensuring that combustion is present and uninterrupted. Two degrees of analysis are presented including a one-dimensional cycle analysis and a complete computational fluid dynamic analysis with finite-rate chemistry and a two-equation turbulence model. Comparison of these analyses show good agreement when the CFD-predicted fuel consumption schedule is used in the cycle analysis.

Molvik, Gregory A.↗
NASA NTRS1993

Axisymmetric computational fluid dynamics analysis of a film/dump-cooled rocket nozzle plume

Prediction of convective base heating rates for a new launch vehicle presents significant challenges to analysts concerned with base environments. The present effort seeks to augment classical base heating scaling techniques via a detailed investigation of the exhaust plume shear layer of a single H2/O2 Space Transportation Main Engine (STME). Use of fuel-rich turbine exhaust to cool the STME nozzle presented concerns regarding potential recirculation of these gases to the base region with attendant increase in the base heating rate. A pressure-based full Navier-Stokes computational fluid dynamics (CFD) code with finite rate chemistry is used to predict plumes for vehicle altitudes of 10 kft and 50 kft. Levels of combustible species within the plume shear layers are calculated in order to assess assumptions made in the base heating analysis.

Tucker, P. K.↗
NASA NTRS1999

The Influence of Extremely Large Solar Proton Events in a Changing Stratosphere. Stratospheric Influence of Solar Proton Events

Two periods of extremely large solar proton events (SPEs) occurred in the past thirty years, which forced significant long-term polar stratospheric changes. The August 2-10, 1972 and October 19-27, 1989 SPEs happened in stratospheres that were quite different chemically. The stratospheric chlorine levels were relatively small in 1972 (approximately 1.2 ppbv) and were fairly substantial in 1989 at about (approximately 3 ppbv). Although these SPEs produced both HO(x) and NO(y) constituents in the mesosphere and stratosphere, only the NO(y) constituents had lifetimes long enough to affect ozone for several months to years past the events. Our recently improved two-dimensional chemistry and transport atmospheric model was used to compute the effects of these gigantic SPEs in a changing stratosphere. Significant upper stratospheric ozone depletions > 10% are computed to last for a few months past these SPEs. The long-lived SPE-produced NO(y) constituents were transported to lower levels during winter after these huge SPEs and caused impacts in the middle and lower stratosphere. During periods of high halogen loading these impacts resulted in interference with the chlorine and bromine loss cycles for ozone destruction. The chemical state of the atmosphere, including the stratospheric sulfate aerosol density, substantially affected the predicted stratospheric influence of these extremely large SPEs.

Jackman, Charles H.↗
NASA NTRS2000

Satellite Measured Solar Protons from 1963-1993 and Their Influence on Ozone in a Changing Stratosphere

The fluxes of solar protons have been measured by a series of Interplanetary Monitoring Platform (IMP) satellites since 1963. Eight IMP satellites have been launched since 1963 providing excellent coverage of solar protons with energies from a few MeV to several hundred MeV. IMP 8, launched in October 1973, continues to provide proton measurements, twenty-six and a half (26 1/2) years after launch. These high energy solar protons rain down on the earth's polar atmosphere sporadically, primarily during solar proton events (SPEs) which typically last a few days. Solar protons with energies of 30 MeV or greater are capable of reaching the stratosphere and causing increases in odd nitrogen (NO(y)) constituents at polar latitudes ($>$ 60 degrees geomagnetic) which last for several months to years past the events. These enhanced NO(y) species can lead to significant upper stratospheric ozone depletions $ greater than $10\% during the gigantic SPEs that occurred in August 1972 and October 1989. We studied the effects of SPEs on the stratosphere from 1963 to the mid 1990s. During this time period, the stratospheric chlorine levels changed from relatively small in 1963 ($\sim$ 1 ppbv) to fairly substantial amounts in the mid-1990s ($\sim$ 3.3 ppbv). Our recently improved two-dimensional chemistry and transport atmospheric model was used to compute the effects of SPEs in this changing stratosphere. The long-lived SPE-produced NO(y) constituents were transported to lower stratospheric levels during winter after the events and caused impacts in the middle and lower stratosphere. Generally, the SPEs resulted in a decrease in ozone. However, during periods of high halogen loading these impacts resulted in interference with the chlorine and bromine loss cycles for ozone destruction. This interference actually led to a predicted total ozone increase that was especially notable in the time period 1992-4, a few years after the October 1989 SPE.

Jackman, Charles H.↗