Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chemically reacting flows”

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 541 records · Page 30

Multigrid Method for Modeling Multi-Dimensional Combustion with Detailed Chemistry

A highly accurate and efficient numerical method is developed for modeling 3-D reacting flows with detailed chemistry. A contravariant velocity-based governing system is developed for general curvilinear coordinates to maintain simplicity of the continuity equation and compactness of the discretization stencil. A fully-implicit backward Euler technique and a third-order monotone upwind-biased scheme on a staggered grid are used for the respective temporal and spatial terms. An efficient semi-coarsening multigrid method based on line-distributive relaxation is used as the flow solver. The species equations are solved in a fully coupled way and the chemical reaction source terms are treated implicitly. Example results are shown for a 3-D gas turbine combustor with strong swirling inflows.

Zheng, Xiaoqing↗

Theoretical kinetic computations in complex reacting systems

Nasa Lewis' studies of complex reacting systems at high temperature are discussed. The changes which occur are the result of many different chemical reactions occurring at the same time. Both an experimental and a theoretical approach are needed to fully understand what happens in these systems. The latter approach is discussed. The differential equations which describe the chemical and thermodynamic changes are given. Their solution by numerical techniques using a detailed chemical mechanism is described. Several different comparisons of computed results with experimental measurements are also given. These include the computation of (1) species concentration profiles in batch and flow reactions, (2) rocket performance in nozzle expansions, and (3) pressure versus time profiles in hydrocarbon ignition processes. The examples illustrate the use of detailed kinetic computations to elucidate a chemical mechanism and to compute practical quantities such as rocket performance, ignition delay times, and ignition lengths in flow processes.

Bittker, David A.↗

Application of a Reynolds Stress turbulence model to a supersonic hydrogen-air diffusion flame

A second-order differential Reynolds Stress turbulence model has been applied to the Favre-averaged Navier-Stokes equations for the study of supersonic flows undergoing hydrogen-air chemical reactions. An assumed Beta Probability Density Function is applied to account for the chemical source terms in the conservation equations. An algebraic Reynolds Flux model is used for the fluctuating density-velocity as well as the species mass fraction-velocity correlations. The variances of temperature and species fluctuations are also modelled using an algebraic flux technique. A seven-species, seven-reaction finite rate chemistry mechanism is used to simulate the combustion processes. The resulting formulation is validated by comparison with experimental data on reacting supersonic axisymmetric jets. Results obtained for specific conditions indicate that the effect of chemical reaction on the turbulence is significant.

Chandrasekhar, R.↗

Ignition and structure of a laminar diffusion flame in the field of a vortex

The distortion of flames in flows with vortical motion is examined via asymptotic analysis and numerical simulation. The model consists of a constant density, one step, irreversible Arrhenius reaction between initially unmixed species occupying adjacent half-planes which are then allowed to mix and react in the presence of a vortex. The evolution in time of the temperature and mass fraction fields is followed. Emphasis is placed on the ignition time and location as a function of vortex Reynolds number and initial temperature differences of the reacting species. The study brings out the influence of the vortex on the chemical reaction. In all phases, good agreement is observed between asymptotic analysis and the full numerical solution of the model equations.

Macaraeg, Michele G.↗

Evaluation of Finite-Rate Gas/Surface Interaction Models for a Carbon Based Ablator

Two sets of finite-rate gas-surface interaction model between air and the carbon surface are studied. The first set is an engineering model with one-way chemical reactions, and the second set is a more detailed model with two-way chemical reactions. These two proposed models intend to cover the carbon surface ablation conditions including the low temperature rate-controlled oxidation, the mid-temperature diffusion-controlled oxidation, and the high temperature sublimation. The prediction of carbon surface recession is achieved by coupling a material thermal response code and a Navier-Stokes flow code. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and Ablation Program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting full Navier-Stokes equations using Data Parallel Line Relaxation method. Recession analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities with heat fluxes ranging from 45 to 1100 wcm2 are performed and compared with data for model validation. The ablating material used in these arc-jet tests is Phenolic Impregnated Carbon Ablator. Additionally, computational predictions of surface recession and shape change are in good agreement with measurement for arc-jet conditions of Small Probe Reentry Investigation for Thermal Protection System Engineering.

Ablation↗

Validation and application of a multiphase CFD model for hydrodynamics, temperature field and RTD simulation in a pilot-scale biomass pyrolysis vapor phase upgrading reactor

Accurate prediction of transport phenomena is critical for VPU reactor design, optimization, and scale-up. The current study focused on the validation and application of a multiphase CFD model within an open-source code MFiX for hydrodynamics, temperature field, and residence time distribution (RTD) simulation in a non-reacting circulating fluidized bed riser for biomass pyrolysis vapor phase upgrading (VPU). First, an Eulerian-Eulerian approach three-dimensional CFD model was employed to simulate the pilot-scale VPU riser on the supercomputer Joule. Excellent quantitative agreement between experimental and simulated results was achieved for pressure drops and temperature field in a range of operating conditions. Then the validated multiphase CFD model was applied to predict gas and solid residence time distributions (RTDs) since prediction and analysis of RTD is an important tool to study the complex multiphase flow behavior and mixing inside chemical reactors. The predictions show that solid mean residence time is 3.5 times the gas residence time; the solid RTD is more sensitive to the process gas flow rate than the solids circulation rate.

09 BIOMASS FUELS↗

Hypersonic Chemically Reacting Boundary-Layer Stability using LASTRAC

Hypersonic boundary layer transition is critical to the design of all hypersonic vehicles due to its effect on the heat transfer into the vehicle surface and potential drag enhancement or reduction during reentry. Boundary layer transition and boundary layer stability analysis under hypersonic conditions has been studied for decades, yet there is ample room for improved accuracy and further investigations into the relevant phenomena. In this work, we present a recent implementation of chemical equilibrium, finite-rate chemistry, and thermochemical nonequilibrium capabilities into LASTRAC, an existing well-established boundary-layer stability analysis code. Verification against existing numerical results in the literature are presented. LASTRAC was previously able to address calorically perfect flows. By using solutions of the Parabolized Stability Equations (PSE) with chemical and thermal nonequilibrium, we are able to investigate the effects of chemical and thermal nonequilibrium on a variety of phenomena including stationary crossflow instability on a swept wing and 2nd mode instabilities over a wedge.

Kline, H. L.↗

The Reduction of Lunar Regolith by Carbothermal Processing Using Methane

The processing of lunar regolith for the production of oxygen is a key component of the In-Situ Resource Utilization plans currently being developed by NASA. In the carbothermal process, a portion of the surface of the regolith in a container is heated by exposure to a heat source so that a small zone of molten regolith is established. A continuous flow of methane is maintained over the molten regolith zone. In this paper, we discuss the development of a chemical conversion model of the carbothermal process to predict the rate of production of carbon monoxide. Our model is based on a mechanism where methane pyrolyzes when it comes in contact with the surface of the hot molten regolith to form solid carbon and hydrogen gas. Carbon is deposited on the surface of the melt, and hydrogen is released into the gas stream above the melt surface. We assume that the deposited carbon mixes in the molten regolith and reacts with metal oxides in a reduction reaction by which gaseous carbon monoxide is liberated. Carbon monoxide bubbles through the melt and is released into the gas stream. It is further processed downstream to ultimately produce oxygen.

Balasubramaniam, R.↗

The Reduction of Lunar Regolith by Carbothermal Processing Using Methane

The processing of lunar regolith for the production of oxygen is a key component of the In-Situ Resource Utilization plans Currently being developed by NASA. In the carbothermal process, a portion of the surface of the regolith in a container is heated by exposure to a heat source so that a small zone of molten regolith is established. A continuous flow of methane is maintained over the molten regolith zone. In this paper, we discuss the development of a chemical conversion model of the carbothermal process to predict the rate of production of carbon monoxide. Our model is based on a mechanism where methane pyrolyzes when it comes in contact with the surface of the hot molten regolith to form solid carbon and hydrogen gas. Carbon is deposited on the surface of the melt, and hydrogen is released into the gas stream above the melt surface. We assume that the deposited carbon mixes in the molten regolith and reacts with metal oxides in a reduction reaction by which gaseous carbon monoxide is liberated. Carbon monoxide bubbles through the melt and is released into the gas stream. It is further processed downstream to ultimately produce oxygen.

Balasubramaniam, R.↗

Application of a Reynolds stress turbulence model to a supersonic radiating hydrogen-air diffusion flame

A second-order differential Reynolds Stress turbulence model has been applied to the Favre-averaged Navier-Stokes equations for the study of supersonic flows with finite-rate chemistry and radiation. An assumed Beta Probability Density Function is applied to account for the chemical source terms and the radiative flux terms in the conservation equations. A seven-species, seven-reaction finite rate chemistry mechanism is used to simulate the combustion process. The tangent slab approximation is used in radiative flux formulation. A pseudo-gray gas model is used to represent the absorption-emission characteristics of the participating species. The turbulence/radiation interaction is achieved via a new formulation. The resulting formulation is validated by comparison with experimental data on reacting supersonic axisymmetric jets. Results obtained for specific conditions indicate that the effect of chemical reaction on the turbulence is significant. Also, the radiative heat transfer is enhanced by the turbulence.

Chandrasekhar, R.↗

Base Heating Sensitivity Study for a 4-Cluster Rocket Motor Configuration in Supersonic Freestream

In support of launch vehicle base heating and pressure prediction efforts using the Loci-CHEM Navier-Stokes computational fluid dynamics solver, 35 numerical simulations of the NASA TND-1093 wind tunnel test have been modeled and analyzed. This test article is composed of four JP-4/LOX 500 lbf rocket motors exhausting into a Mach 2 - 3.5 wind tunnel at various ambient pressure conditions. These water-cooled motors are attached to a base plate of a standard missile forebody. We explore the base heating profiles for fully coupled finite-rate chemistry simulations, one-way coupled RAMP (Reacting And Multiphase Program using Method of Characteristics)-BLIMPJ (Boundary Layer Integral Matrix Program - Jet Version) derived solutions and variable and constant specific heat ratio frozen flow simulations. Variations in turbulence models, temperature boundary conditions and thermodynamic properties of the plume have been investigated at two ambient pressure conditions: 255 lb/sq ft (simulated low altitude) and 35 lb/sq ft (simulated high altitude). It is observed that the convective base heat flux and base temperature are most sensitive to the nozzle inner wall thermal boundary layer profile which is dependent on the wall temperature, boundary layer s specific energy and chemical reactions. Recovery shock dynamics and afterburning significantly influences convective base heating. Turbulence models and external nozzle wall thermal boundary layer profiles show less sensitivity to base heating characteristics. Base heating rates are validated for the highest fidelity solutions which show an agreement within +/-10% with respect to test data.

Mehta, Manish↗

Reacting viscous-shock-layer solutions with multicomponent diffusion and mass injection

Numerical solutions of the viscous-shock-layer equation where the chemistry is treated as being either frozen, equilibrium, or nonequilibrium are presented. Also the effects of the diffusion model, surface catalysis, and mass injection on surface transport and flow parameters are considered. The flow is treated as a mixture of five inert and thermally perfect species. The viscous-shock-layer equations are solved by using an implicit-difference scheme. All calculations are for hyperboloids with included angles of 20 and 45. The flight conditions are those for various altitudes and velocities in the earth's atmosphere. Data are presented to show the effects of the chemical models; diffusion models; surface catalysis; and mass injection of air on heat transfer; skin friction; shock standoff distance; wall pressure distribution; and tangential velocity, temperature, and species profiles. The results show that an equilibrium analysis can substantially overpredict the heat-transfer rates for flow conditions experienced by earth-orbital entry vehicles. Moreover, at such conditions surface catalysis significantly influences heat-transfer and flow-field properties. If a binary rather than a multicomponent diffusion model is assumed, negligible errors in most flow properties result. Quantitative results are presented that show the effect of mass injection on flow properties within and downstream of the injection region.

Moss, J. N.↗

Effects of unburnt reaction progress on stretch flame dynamics under elevated temperatures

Flames in practical combustors and engines are inevitably subject to the effects of upstream chemical reaction progress and stretch, induced by elevated thermodynamic conditions and flow non-uniformities, respectively. Recent shock tube experiments and simulation studies on flame propagation have shown that flame speed under engine-relevant conditions can be enhanced with non-negligible upstream chemical reaction progress, especially when low-temperature heat release is involved in the unburnt mixture. On the other hand, depending on the mixture equivalence ratio and diluents, nonequidiffusion (including the non-unity Lewis number effect and the preferential diffusion effect) can couple with the flame stretch to fundamentally affect the flame propagation, which can manifest as either facilitation or suppression. Depending on the transport property of a reacting mixture, there hence can be either inhibition or promotion from unburnt reaction progress and stretch effects on flame propagation. In the current work, through one-dimensional numerical simulations of transient planar and spherical flames of n-heptane/air under elevated thermodynamic conditions, the combined effects of upstream chemical reaction progress and stretch on flame propagation are investigated. Results show that for both lean and rich n-heptane/air mixtures, flame speed can be substantially promoted with reaction progress, while the rich mixtures can exhibit opposite stretch dependence after first-stage ignition. Different definitions of effective Lewis number are adopted to explain the change in Markstein length for spherical flame in lean and rich mixtures reformed by the low temperature chemistry. Further, this work fills an important gap in laminar premixed flame research relevant to practical combustion systems and can provide useful insight into local turbulent flame behaviors and phi-sensitivity in engine combustion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics and structure of turbulent premixed flames

In earlier work (Mantel & Bilger, 1994) the structure of the turbulent premixed flame was investigated using statistics based on conditional averaging with the reaction progress variable as the conditioning variable. The DNS data base of Trouve and Poinsot (1994) was used in this investigation. Attention was focused on the conditional dissipation and conditional axial velocity in the flame with a view to modeling these quantities for use in the conditional moment closure (CMC) approach to analysis of kinetics in premixed flames (Bilger, 1993). Two remarkable findings were made: there was almost no acceleration of the axial velocity in the flame front itself; and the conditional scalar dissipation remained as high, or higher, than that found in laminar premixed flames. The first finding was surprising since in laminar flames all the fluid acceleration occurs through the flame front, and this could be expected also for turbulent premixed flames at the flamelet limit. The finding gave hope of inventing a new approach to the dynamics of turbulent premixed flames through use of rapid distortion theory or an unsteady Bernoulli equation. This could lead to a new second order closure for turbulent premixed flames. The second finding was contrary to our measurements with laser diagnostics in lean hydrocarbon flames where it is found that conditional scalar dissipation drops dramatically below that for laminar flamelets when the turbulence intensity becomes high. Such behavior was not explainable with a one-step kinetic model, even at non-unity Lewis number. It could be due to depletion of H2 from the reaction zone by preferential diffusion. The capacity of the flame to generate radicals is critically dependent on the levels of H2 present (Bilger, et al., 1991). It seemed that a DNS computation with a multistep reduced mechanism would be worthwhile if a way could be found to make this feasible. Truly innovative approaches to complex problems often come only when there is the opportunity to work close at hand with the (in this case numerical) experimental data. Not only can one spot patterns and relationships in the data which could be important, but one can also get to know the limitations of the technique being used, so that when the next experiment is being designed it will address resolvable questions. A three-year grant from the Australian Research Council has enabled us to develop a small capability at the University of Sydney to work on DNS of turbulent reacting flow, and to analyze data bases generated at CTR. Collaboration between the University of Sydney and CTR is essential to this project and finding a workable modus operandum for this collaboration, given the constraints involved, has been a major objective of the past year's effort. The overall objectives of the project are: (1) to obtain a quantitative understanding of the dynamics of turbulent premixed flames at high turbulence levels with a view to developing improved second order closure models; and (2) to carry out new DNS experiments on turbulent premixed flames using a carefully chosen multistep reduced mechanism for the chemical kinetics, with a view to elucidating the laser diagnostic findings that are contrary to the findings for DNS using one-step kinetics. In this first year the objectives have been to make the existing CTR data base more accessible to coworkers at the University of Sydney, to make progress on understanding the dynamics of the flame in this existing CTR data base, and to carefully construct a suitable multistep reduced mechanism for use in a new set of DNS experiments on turbulent premixed flames.

Bilger, R. W.↗

Intumescent Coatings as Fire Retardants

The development of fire-retardant coatings to protect surfaces which may be exposed to fire or extreme heat is a subject of intense interest to many industries. A fire-retardant paint has been developed which represents a new chemical approach for preparing intumescent coatings, and potentially, is very important to fire-prevention authorities. The requirements for a superior coating include ease of application, suitability to a wide variety of surfaces and finishes, and stability over an extended period of time within a broad range of ambient temperature and humidity conditions. These innovative coatings, when activated by the heat of a fire, react to form a thick, low-density, polymeric coating or char layer. Water vapor and sulphur dioxide are released during the intumescent reaction. Two fire-protection mechanisms thus become available: (1) the char layer retards the flow of heat, due to the extremely low thermal conductivity; and (2) water vapor and sulfur dioxide are released, providing fire quenching properties. Still another mechanism functions in cases where the char, by virtue of its high oxidation resistance and low thermal conductivity, reaches a sufficiently high temperature to re-radiate much of the incident heat load. The coatings consist of dispersions of selective salts of a nitro-amino-arornatic compound. Specifically, para-nitroaniline bisulfate and the ammonium salt of para-nitroaniline-ortho sulphuric acid (2-amino-5-nitrobenzenesulphuric acid) are used. Suitable vehicles are cellulose nitrate of lacquer grade, a nitrite-phenolic modified rubber, or epoxy-polysulfide copolymer. Three separate formulations have been developed. A solvent is usually employed, such as methylethyl ketone, butyl acetate, or toluene, which renders the coatings suitably thin and which evaporates after the coatings are applied. Generally, the intumescent material is treated as insoluble in the vehicle, and is ground and dispersed in the vehicle and solvent like an ordinary coating pigment. The char found on intumescence is better in terms of yield and physical properties than chars obtained from many previously known intumescent materials. Prior to intumescence, the coating has a density of 85 pounds per cubic foot. After intumescence, the density is approximately 0.3 pounds per cubic loot. The linear expansion of the coatings ranges from 70 to 200 times the applied coating thickness.

Parker, J. A.↗

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

Surface modeling and grid generation for aeropropulsion CFD

The efforts in geometry modeling and grid generation at the NASA Lewis Research Center, as applied to the computational fluid dynamic (CFD) analysis of aeropropulsion systems, are presented. The efforts are mainly characterized by a focus on the analysis of components of an aeropropulsion system, which involve turbulent viscous flow with heat transfer and chemistry. Thus, this discussion will follow that characterization and will sequence through the components of typical propulsion systems consisting of inlets, compressors, combustors, turbines, and nozzles. For each component, some applications of CFD analysis will be presented to show how CFD is used to compute the desired performance information, how geometry modeling and grid generation are performed, and what issues have developed related to geometry modeling and grid generation. The discussion will illustrate the following needs related to geometry modeling and grid generation as observed in aeropropulsion analysis: (1) accurate and efficient resolution of turbulent viscous and chemically-reacting flowfields; (2) easy-to-use interfaces with CAD data for automated grid generation about complex geometries; and (3) automated batch grid generation software for use with design and optimization software.

Choo, Yung K.↗

Passive components for high temperature operation

The development of resistors, capacitors, and interconnect metalizations utilizing a low pressure chemical vapor deposition (LPCVD) process is described. The LPCVD process allows the passive components to be fabricated at temperatures higher than their highest operating temperatures. The deposition of thin films by LPCVD is accomplished by reacting one or more gases on the surface of a heated substrate. The substrates to be coated are placed on a graphite susceptor and then loaded into the center of a quartz reaction tube. The RF power is applied to the coil on the outside of the reaction tube which in turn is coupled into the graphite susceptor causing it to heat. Pressures of several torr or less are typical, with carrier flow rates of 0.1 to 2.0 liters/min. Nitrogen, hydrogen, and helium are typical carrier gases. These are controlled with mass flow controllers and the pressure is continuously monitored with a capacitive manometer. Materials selection and device testing are also discussed.

Raymond, L. S.↗