Search NASA⌕ Search

SEARCH · Search NASA

Results for “flow 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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Development of Comprehensive Reduced Kinetic Models for Supersonic Reacting Shear Layer Simulations

Large-scale simulations of multi-dimensional unsteady turbulent reacting flows with detailed chemistry and transport can be computationally extremely intensive even on distributed computing architectures. With the development of suitable reduced chemical kinetic models, the number of scalar variables to be integrated can be decreased, leading to a significant reduction in the computational time required for the simulation with limited loss of accuracy in the results. A general MATLAB-based automated mechanism reduction procedure is presented to reduce any complex starting mechanism (detailed or skeletal) with minimal human intervention. Based on the application of the quasi steady-state (QSS) approximation for certain chemical species and on the elimination of the fast reaction rates in the mechanism, several comprehensive reduced models, capable of handling different fuels such as C2H4, CH4 and H2, have been developed and thoroughly tested for several combustion problems (ignition, propagation and extinction) and physical conditions (reactant compositions, temperatures, and pressures). A key feature of the present reduction procedure is the explicit solution of the concentrations of the QSS species, needed for the evaluation of the elementary reaction rates. In contrast, previous approaches relied on an implicit solution due to the strong coupling between QSS species, requiring computationally expensive inner iterations. A novel algorithm, based on the definition of a QSS species coupling matrix, is presented to (i) introduce appropriate truncations to the QSS algebraic relations and (ii) identify the optimal sequence for the explicit solution of the concentration of the QSS species. With the automatic generation of the relevant source code, the resulting reduced models can be readily implemented into numerical codes.

Zambon, A. C.↗
NASA NTRS2018

One-Dimensional Modeling Methodology for Shock Tubes: Application to the EAST Facility

In this work, a one-dimensional methodology for simulating shock tubes is developed. The model accounts for the viscous interactions of the shock with the shock tube wall by adding an area change source term in the 1-D conservation equations corresponding to the boundary layer growth. This source term corresponds to the mass and energy going into the boundary layer. The boundary layer growth is computed using a simple model with a scaling factor. This scale factor is used to tailor a solution to match the deceleration profile of a shock tube test. In doing so, not only will the source term take into account boundary layer losses, it will also cover any effect due to radiative cooling loses from the gas. For this study, the Electric Arc Shock Tube(EAST) facility at NASA Ames Research Center is modeled for Earth reentry conditions. The purpose of this paper is to investigate if anomalies identified for certain conditions in the EAST data are due to shock deceleration. These anomalies include measuring electron number density above equilibrium predictions and observing that radiance profiles can continually increase behind the shock, never reach steady state, for certain shots (typically those less than 10 km/s). An eleven species air mixture is chosen to study the chemistry of the flow. Comparisons of the simulations to the experimental results are presented. Good agreement with the shock deceleration profiles was achieved by tuning in the boundary layer scale factor. The temperature as well as electron number density increases behind the shock, as has also been observed in the experiments. Finally, radiance comparisons between results from NEQAIR and experiments also show good agreement for some shots, but significant discrepancies are still observed for others.

Sharma, Maitreyee↗
NASA NTRS2018

One-Dimensional Modeling Methodology for Shock Tubes: Application to the EAST Facility

In this work, a one-dimensional methodology for simulating shock tubes is developed. The model accounts for the viscous interactions of the shock with the shock tube wall by adding an area change source term in the 1-D conservation equations corresponding to the boundary layer growth. This source term corresponds to the mass and energy going into the boundary layer. The boundary layer growth is computed using a simple model with a scaling factor. This scale factor is used to tailor a solution to match the deceleration profile of a shock tube test. In doing so, not only will the source term take into account boundary layer losses, it will also cover any effect due to radiative cooling loses from the gas. For this study, the Electric Arc Shock Tube(EAST) facility at NASA Ames Research Center is modeled for Earth reentry conditions. The purpose of this paper is to investigate if anomalies identified for certain conditions in the EAST data are due to shock deceleration. These anomalies include measuring electron number density above equilibrium predictions and observing that radiance profiles can continually increase behind the shock, never reach steady state, for certain shots (typically those less than 10 km/s). An eleven species air mixture is chosen to study the chemistry of the flow. Comparisons of the simulations to the experimental results are presented. Good agreement with the shock deceleration profiles was achieved by tuning in the boundary layer scale factor. The temperature as well as electron number density increases behind the shock, as has also been observed in the experiments. Finally, radiance comparisons between results from NEQAIR and experiments also show good agreement for some shots, but significant discrepancies are still observed for others.

Sharma, Maitreyee↗
NASA NTRS1991

Monte Carlo analysis of dissociation and recombination behind strong shock waves in nitrogen

Computations are presented for the relaxation zone behind strong, 1D shock waves in nitrogen. The analysis is performed with the direct simulation Monte Carlo method (DSMC). The DSMC code is vectorized for efficient use on a supercomputer. The code simulates translational, rotational and vibrational energy exchange and dissociative and recombinative chemical reactions. A model is proposed for the treatment of three body-recombination collisions in the DSMC technique which usually simulates binary collision events. The model improves previous models because it can be employed with a large range of chemical-rate data, does not introduce into the flow field troublesome pairs of atoms which may recombine upon further collision (pseudoparticles) and is compatible with the vectorized code. The computational results are compared with existing experimental data. It is shown that the derivation of chemical-rate coefficients must account for the degree of vibrational nonequilibrium in the flow. A nonequilibrium-chemistry model is employed together with equilibrium-rate data to compute the flow in several different nitrogen shock waves.

Boyd, I. D.↗
NASA NTRS1986

Time-dependent solution for axisymmetric flow over a blunt body with ideal gas, CF4, or equilibrium air chemistry

A time-asymptotic method has been used to obtain steady-flow solutions for axisymmetric inviscid flow over several blunt bodies including spheres, paraboloids, ellipsoids, and spherically blunted cones. Comparisons with experimental data and results of other computational methods have demonstrated that accurate solutions can be obtained using this approach. The method should prove useful as an analysis tool for comparing with experimental data and for making engineering calculations for blunt reentry vehicles.

Hamilton, H. H., II↗
NASA NTRS1994

Modified SEAGULL

Original version of program incorporated into program SRGULL (LEW-15093) for use on National Aero-Space Plane project, its duty being to model forebody, inlet, and nozzle portions of vehicle. However, real-gas chemistry effects in hypersonic flow fields limited accuracy of that version, because it assumed perfect-gas properties. As a result, SEAGULL modified according to real-gas equilibrium-chemistry methodology. This program analyzes two-dimensional, hypersonic flows of real gases. Modified version of SEAGULL maintains as much of original program as possible, and retains ability to execute original perfect-gas version.

Salas, M. D.↗
NASA NTRS1987

Chemistry-split techniques for viscous reactive blunt body flow computations

The weak-coupling structure between the fluid and species equations has been exploited and resulted in three, closely related, time-iterative implicit techniques. While the primitive variables are solved in two separated groups and each by an Alternating Direction Implicit (ADI) factorization scheme, the rate-species Jacobian can be treated in either full or diagonal matrix form, or simply ignored. The latter two versions render the split technique to solving for species as scalar rather than vector variables. The solution is completed at the end of each iteration after determining temperature and pressure from the flow density, energy and species concentrations. Numerical experimentation has shown that the split scalar technique, using partial rate Jacobian, yields the best overall stability and consistency. Satisfactory viscous solutions were obtained for an ellipsoidal body of axis ratio 3:1 at Mach 35 and an angle of attack of 20 degrees.

Li, C. P.↗
NASA NTRS1974

A source flow characteristic technique for the analysis of scramjet exhaust flow field

The factors which influence the design and selection of a nozzle for a hypersonic scramjet are described. A two dimensional second-order characteristic procedure capable of analyzing the aerodynamic performance of typical nozzle configurations is presented. Equations of motion governing the two dimensional, axisymmetric, or axially expanding inviscid flow of a gas mixture, with frozen chemistry, are provided. Diagrams of the flow conditions for various configurations are included.

Delguidice, P.↗
NASA NTRS1992

A code validation strategy and facility for nonequilibrium, reacting flows

The rotational and vibrational population distributions in a hypersonic flow of hot iodine vapor were investigated using a pilot wind tunnel. The tunnel provides run times of about 20 min, with a scale-up to a larger facility appearing to be possible. It will provide a 20-cm-diameter, Mach 9 flow of iodine vapor with only small amounts of free stream nonequilibrium. An analysis of the flow response to a normal shock wave indicates that significant chemistry will occur in model flow fields in the larger facility.

Muntz, E. P.↗
NASA NTRS1989

A detailed analysis of inviscid flux splitting algorithms for real gases with equilibrium or finite-rate chemistry

The extension of the known flux-vector and flux-difference splittings to real gases via rigorous mathematical procedures is demonstrated. Formulations of both equilibrium and finite-rate chemistry for real-gas flows are described, with emphasis on derivations of finite-rate chemistry. Split-flux formulas from other authors are examined. A second-order upwind-based TVD scheme is adopted to eliminate oscillations and to obtain a sharp representation of discontinuities.

Shuen, Jian-Shun↗
NASA NTRS1993

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.↗
NASA NTRS1979

Combustor modelling for scramjet engines

A system of computer programs is being developed to analyse and predict the complex flow fields found in hydrogen-fueled scramjet combustors. Each program is designed to solve the governing equation system for the type of flow present in a particular combustor region. A two-dimensional parabolic program has been found to be valuable in the development and experimental evaluation of turbulence and chemistry models for supersonic flow, and in the development of a program to model supersonic flow downstream of the fuel injection struts by means of solutions to the three-dimensional parabolic Navier-Stokes equations and species equations. A partially elliptic code has been derived to account for local subsonic flow regions, and fully elliptic programs have been developed by the consideration of streamwise diffusion effects for the recirculating flow fields near transverse fuel injectors. The programs are currently being applied to problems of scramjet engine development.

Drummond, J. P.↗
NASA NTRS1989

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 chemically react and reach states in thermal nonequilibrium. The prediction of hypersonic flowfields requires a numerical method capable of solving the conservation equations of fluid flow, the chemical rate equations for specie formation and dissociation, and the transfer of energy relations between translational and vibrational temperature states. Because the number of equations to be solved is large, the numerical method should also be as efficient as possible. The proposed paper presents a fully implicit method that fully couples the solution of the fluid flow equations with the gas physics and chemistry relations. The method flux splits the inviscid flow terms, central differences of the viscous terms, preserves element conservation in the strong chemistry source terms, and solves the resulting block matrix equation by Gauss Seidel line relaxation.

Maccormack, Robert W.↗
NASA NTRS1992

Finite-Difference Solution for Laminar or Turbulent Boundary Layer Flow over Axisymmetric Bodies with Ideal Gas, CF4, or Equilibrium Air Chemistry

A computer code was developed that uses an implicit finite-difference technique to solve nonsimilar, axisymmetric boundary layer equations for both laminar and turbulent flow. The code can treat ideal gases, air in chemical equilibrium, and carbon tetrafluoride (CF4), which is a useful gas for hypersonic blunt-body simulations. This is the only known boundary layer code that can treat CF4. Comparisons with experimental data have demonstrated that accurate solutions are obtained. The method should prove useful as an analysis tool for comparing calculations with wind tunnel experiments and for making calculations about flight vehicles where equilibrium air chemistry assumptions are valid.

Hamilton, H. Harris, II↗
NASA NTRS1987

An analysis of artificial viscosity effects on reacting flows using a spectral multi-domain technique

Analytical and experimental techniques for modeling the aerothermodynamics of hypersonic flight are assessed, together with the problems which will be encountered in developing reusable hypersonic vehicles. Emphasis is placed on a numerical coupling between a nonequilibrium chemistry model and hypersonic flow kinematics. Finite difference and finite element descriptions of flow fields in which molecules encounter a shock wave and undergo various motion (and thereby energy) transformations are discussed. The effects of artificial smearing of the shock wave are considered in terms of the resulting effects on the distribution of the energies and chemical composition of the transition region. Results are provided from schlieren photographs of shock-tube experiments, Navier-Stokes calculations of axisymmetric flow over a conical body, and calculations using a spectral multidomain approach for chemically reacting flows.

Macaraeg, Michele G.↗
NASA NTRS2020

Automated Identification and Differentiation of Spectrally Similar Hydrothermal Minerals on Mars

Early telescopic observations corroborated hydration related absorptions on Mars in the infrared. Images from the Viking missions led to speculation of hydrothermal alteration and were followed by two missions which mapped the spatial variability of the ~ 3 μm hydration feature. Since then, the Compact Reconnaissance Imager for Mars (CRISM) has provided high spatial resolution (up to 18m) spectral identification of a suite of hydrothermal and diagenetic minerals which have illuminated a range of formation mechanisms. Presence/absence and spatial segregation or mixing of minerals like prehnite, epidote, chlorite amphiboles, and mixed-layer Fe/Mg smectite-chlorite provide valuable evidence for the geologic setting of deposits on Earth, and these phases are often used as temperature and aqueous chemistry indicators in terrestrial systems. Mapping the distribution of these phases will help to answer whether Mars had widespread conditions favorable for low-grade metamorphism and diagenesis, or only focused hydrothermal systems in areas of high heat flow. Further characterizing the chemistry and structure of these phases will then help to answer how most of the widespread Fe/Mg phyllosilicates formed, further defining early geochemical cycling and climate. A fully automated approach for accurate mapping of important hydrothermal mineral phases on Mars has been a challenge. Due to overlapping features in the M-OH region (~2.2-2.4 μm), the strongest absorption features of chlorite, prehnite, and epidote in the short-wave infrared are difficult to distinguish from one another and from the most commonly occurring hydrated silicates on Mars, Fe/Mg smectites. Weaker absorptions are present in both prehnite and epidote which help to distinguish them from chlorite and smectites, but their relative strength in the presence of noise and spatial mixing is often too low to confidently identify them without the noise suppression and feature enhancement methods described here. The spectral signatures of mixed-layer Fe/Mg smectite-chlorite and partially chloritized Fe/Mg smectites have not yet been adequately assessed. Here we evaluate the effectiveness of two empirical and statistical methods for identifying and differentiating these phases using CRISM data.

Rasmussen, B. P.↗
NASA NTRS2017

Fully Coupled Aero-Thermochemical-Elastic Simulations of an Eroding Graphite Nozzle

A multiphysics simulation capability has been developed that incorporates mutual interactions between aerodynamics, structural response from aero/thermal loading, ablation/pyrolysis, heating, and surface-to-surface radiation to perform high-fidelity, fully coupled aerothermoelastic ablation simulations, which to date had been unattainable. The multiphysics framework couples CHAR (a 3-D implicit charring ablator solver), Loci/CHEM (a computational fluid dynamics solver for high-speed chemically reacting flows), and Abaqus (a nonlinear structural dynamics solver) to create a fully coupled aerothermoelastic charring ablative solver. The solvers are tightly coupled in a fully integrated fashion to resolve the effects of the ablation pyrolysis and charring process and chemistry products upon the flow field, the changes in surface geometry due to recession upon the flow field, and thermal-structural analysis of the body from the induced aerodynamic heating from the flow field. The multiphysics framework was successfully demonstrated on a solid rocket motor graphite nozzle erosion application. Comparisons were made with available experimental data that measured the throat erosion during the motor firing. The erosion data is well characterized, as the test rig was equipped with a windowed nozzle section for real-time X-ray radiography diagnostics of the instantaneous throat variations for deducing the instantaneous erosion rates. The nozzle initially undergoes a nozzle contraction due to thermal expansion before ablation effects are able to widen the throat. A series of parameters studies were conducted using the coupled simulation capability to determine the sensitivity of the nozzle erosion to different parameters. The parameter studies included the shape of the nozzle throat (flat versus rounded), the material properties, the effect of the choice of turbulence model, and the inclusion or exclusion of the mechanical thermal expansion. Overall, the predicted results match the experiment very well, and the predictions were able to bound the data within acceptable limits.

Blades, E. L.↗
NASA NTRS1993

Axisymmetric computational fluid dynamics analysis of Saturn V/S1-C/F1 nozzle and plume

An axisymmetric single engine Computational Fluid Dynamics calculation of the Saturn V/S 1-C vehicle base region and F1 engine plume is described. There were two objectives of this work, the first was to calculate an axisymmetric approximation of the nozzle, plume and base region flow fields of S1-C/F1, relate/scale this to flight data and apply this scaling factor to a NLS/STME axisymmetric calculations from a parallel effort. The second was to assess the differences in F1 and STME plume shear layer development and concentration of combustible gases. This second piece of information was to be input/supporting data for assumptions made in NLS2 base temperature scaling methodology from which the vehicle base thermal environments were being generated. The F1 calculations started at the main combustion chamber faceplate and incorporated the turbine exhaust dump/nozzle film coolant. The plume and base region calculations were made for ten thousand feet and 57 thousand feet altitude at vehicle flight velocity and in stagnant freestream. FDNS was implemented with a 14 species, 28 reaction finite rate chemistry model plus a soot burning model for the RP-1/LOX chemistry. Nozzle and plume flow fields are shown, the plume shear layer constituents are compared to a STME plume. Conclusions are made about the validity and status of the analysis and NLS2 vehicle base thermal environment definition methodology.

Ruf, Joseph H.↗