Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive flow”

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 37 records · Page 2

Stability and Interaction of Coherent Structure in Supersonic Reactive Wakes

A theoretical formulation and analysis is presented for a study of the stability and interaction of coherent structure in reacting free shear layers. The physical problem under investigation is a premixed hydrogen-oxygen reacting shear layer in the wake of a thin flat plate. The coherent structure is modeled as a periodic disturbance and its stability is determined by the application of linearized hydrodynamic stability theory which results in a generalized eigenvalue problem for reactive flows. Detailed stability analysis of the reactive wake for neutral, symmetrical and antisymmetrical disturbance is presented. Reactive stability criteria is shown to be quite different from classical non-reactive stability. The interaction between the mean flow, coherent structure and fine-scale turbulence is theoretically formulated using the von-Kaman integral technique. Both time-averaging and conditional phase averaging are necessary to separate the three types of motion. The resulting integro-differential equations can then be solved subject to initial conditions with appropriate shape functions. In the laminar flow transition region of interest, the spatial interaction between the mean motion and coherent structure is calculated for both non-reactive and reactive conditions and compared with experimental data wherever available. The fine-scale turbulent motion determined by the application of integral analysis to the fluctuation equations. Since at present this turbulence model is still untested, turbulence is modeled in the interaction problem by a simple algebraic eddy viscosity model. The applicability of the integral turbulence model formulated here is studied parametrically by integrating these equations for the simple case of self-similar mean motion with assumed shape functions. The effect of the motion of the coherent structure is studied and very good agreement is obtained with previous experimental and theoretical works for non-reactive flow. For the reactive case, lack of experimental data made direct comparison difficult. It was determined that the growth rate of the disturbance amplitude is lower for reactive case. The results indicate that the reactive flow stability is in qualitative agreement with experimental observation.

Menon, Suresh↗

LES, DNS, and RANS for the Analysis of High-Speed Turbulent Reacting Flows

A filtered density function (FDF) method suitable for chemically reactive flows is developed in the context of large eddy simulation. The advantage of the FDF methodology is its inherent ability to resolve subgrid scales (SGS) scalar correlations that otherwise have to be modeled. Because of the lack of robust models to accurately predict these correlations in turbulent reactive flows, simulations involving turbulent combustion are often met with a degree of skepticism. The FDF methodology avoids the closure problem associated with these terms and treats the reaction in an exact manner. The scalar FDF approach is particularly attractive since it can be coupled with existing hydrodynamic computational fluid dynamics (CFD) codes.

Colucci, P. J.↗

Comparison of PDF and Moment Closure Methods in the Modeling of Turbulent Reacting Flows

In modeling turbulent reactive flows, Probability Density Function (PDF) methods have an advantage over the more traditional moment closure schemes in that the PDF formulation treats the chemical reaction source terms exactly, while moment closure methods are required to model the mean reaction rate. The common model used is the laminar chemistry approximation, where the effects of turbulence on the reaction are assumed negligible. For flows with low turbulence levels and fast chemistry, the difference between the two methods can be expected to be small. However for flows with finite rate chemistry and high turbulence levels, significant errors can be expected in the moment closure method. In this paper, the ability of the PDF method and the moment closure scheme to accurately model a turbulent reacting flow is tested. To accomplish this, both schemes were used to model a CO/H2/N2- air piloted diffusion flame near extinction. Identical thermochemistry, turbulence models, initial conditions and boundary conditions are employed to ensure a consistent comparison can be made. The results of the two methods are compared to experimental data as well as to each other. The comparison reveals that the PDF method provides good agreement with the experimental data, while the moment closure scheme incorrectly shows a broad, laminar-like flame structure.

Norris, Andrew T.↗

A Gas-Kinetic Scheme for Multimaterial Flows and Its Application in Chemical Reaction

This paper concerns the extension of the multicomponent gas-kinetic BGK-type scheme to multidimensional chemical reactive flow calculations. In the kinetic model, each component satisfies its individual gas-kinetic BGK equation and the equilibrium states of both components are coupled in space and time due to the momentum and energy exchange in the course of particle collisions. At the same time, according to the chemical reaction rule one component can be changed into another component with the release of energy, where the reactant and product could have different gamma. Many numerical test cases are included in this paper, which show the robustness and accuracy of kinetic approach in the description of multicomponent reactive flows.

Lian, Yongsheng↗

A fast sampling device for the mass spectrometric analysis of liquid rocket engine exhaust

The design of a device to obtain compositional data on rocket exhaust by direct sampling of reactive flow exhausts into a mass spectrometer is presented. Sampling at three stages differing in pressure and orifice angle and diameter is possible. Results of calibration with pure gases and gas mixtures are erratic and of unknown accuracy for H2, limiting the usefulness of the apparatus for determining oxidizer/fuel ratios from combustion product analysis. Deposition effects are discussed, and data obtained from rocket exhaust spectra are analyzed to give O/F ratios and mixture ratio distribution. The O/F ratio determined spectrometrically is insufficiently accurate for quantitative comparison with cold flow data. However, a criterion for operating conditions with improved mixing of fuel and oxidizer which is consistent with cold flow results may be obtained by inspection of contour plots. A chemical inefficiency in the combustion process when oxidizer is in excess is observed from reactive flow measurements. Present results were obtained with N2O4/N2H4 propellants.

Ryason, P. R.↗

Systematic development of reduced reaction mechanisms for dynamic modeling

A method for systematically developing a reduced chemical reaction mechanism for dynamic modeling of chemically reactive flows is presented. The method is based on the postulate that if a reduced reaction mechanism faithfully describes the time evolution of both thermal and chain reaction processes characteristic of a more complete mechanism, then the reduced mechanism will describe the chemical processes in a chemically reacting flow with approximately the same degree of accuracy. Here this postulate is tested by producing a series of mechanisms of reduced accuracy, which are derived from a full detailed mechanism for methane-oxygen combustion. These mechanisms were then tested in a series of reactive flow calculations in which a large-amplitude sinusoidal perturbation is applied to a system that is initially quiescent and whose temperature is high enough to start ignition processes. Comparison of the results for systems with and without convective flow show that this approach produces reduced mechanisms that are useful for calculations of explosions and detonations. Extensions and applicability to flames are discussed.

Frenklach, M.↗

OMV/VTE variable thrust engine analysis

The objective of the present work is to develop a predictive CFD based analytical tool for the variable thrust engine (VTE) in the orbital maneuvering vehicle (OMV). This objective is being accomplished within the framework of the Los Alamos KIVA computer code for chemically reactive flows with sprays. For the OMV application, the main structure of KIVA is to be retained while reformulating many of the phenomenological submodels, enhancing some of the numerics, and adding more features. The analytical model consists of the general conservation equations for two phase reactive flows and of submodels for turbulence, chemical reactions, and bipropellant sprays. Tailoring this model to the OMV engine brings about the added complexities of combustion and flow processes that occur in a liquid hypergolic propellant rocket chamber. This report exposes the foundation upon which the analytical tool is being constructed and developed. Results from a cursory computational exercise involving the simulation of the flow and combustion processes in a hypothetical N2H4/N2O4 rocket engine thrust chamber is presented and discussed.

Larosillere, Louis↗

OMV/VTE variable thrust engine analysis

The objective of the present work is to develop a predictive CFD based analytical tool for the Variable Thrust Engine (VTE) in the Orbital Maneuvering Vehicle (OMV). This objective is being accomplished within the framework of the Los Alamos KIVA computer code for chemically reactive flows with sprays. For the OMV application, the main structure of KIVA is to be retained while reformulating many of the phenomenological submodels, enhancing some of the numerics, and adding more features. The analytical model consists of the general conservation equations for two-phase reactive flows and of submodels for turbulence, chemical reactions, and bipropellant sprays. Tailoring this model to the OMV engine brings about the added complexities of combustion and flow processes that occur in a liquid hypergolic propellant rocket chamber. This report exposes the foundation upon which the analytical tool is being constructed and developed. Results from a cursory computational exercise involving the simulation of the flow and combustion processes in a hypothetical N2H4/N204 rocket engine thrust chamber is presented and discussed.

Larosillere, L.↗

Modeling of Alkane Oxidation Using Constituents and Species

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

Bellan, Jasette↗

Numerical solutions of reactive fluid flows during postignition transients in hybrid rocket systems.

A computational method has been developed for the study of the post-ignition transients in hybrid rocket systems. The particular system chosen consisted of a gaseous oxidizer flowing within a tube of solid fuel, resulting in heterogeneous combustion. With the appropriate assumptions, two-dimensional, time-dependent conservation equations were derived for the reacting gas phase, and for the solid phase, in a cylindrical coordinate system. These were then programmed for numerical computation, using two implicit finite-difference schemes, the Lax-Wendroff scheme for the gas phase, and the Crank-Nicolson scheme for the solid phase. Appropriate initial and boundary conditions were represented, including heat and mass conservation at the interface between gas and solid. Initially, no attempt was made to relate the recession rate at the surface to the surface temperature, or to include heat transfer by radiation. A simple case was selected for preliminary calculations, with aluminum and oxygen as fuel and oxidizer, and aluminum oxide as the product.

Hung, W. S. Y.↗

Transport phenomena of reactive fluid flow in heterogeneous combustion processes.

A previously developed computer program was used to model two transient hybrid combustion processes involving tubes of solid Plexiglas. In the first study, representing combustion of a hybrid rocket, the oxidizing gas was oxygen, and calculations were continued sufficiently long to obtain steady-state values. Systematic variations were made in reaction rate constant, mass flow rate, and pressure, alternatively using constant and temperature dependent regression rate models for the fuel surface. Consistent results were obtained, as is evidenced by the values for the mass function of the reaction product and the flame temperature, for which plots are supplied. In the second study, fire initiation in a duct was studied, with an air mixture as the oxidizing gas. It was demonstrated that a satisfactory flame spread mechanism could be reproduced on the computer. In both of the above applications, the general, transient, two-dimensional conservation equations were represented, together with chemical reactions, solid-fuel interface conditions, and heat conduction in the solid fuel.

Hung, W. S. Y.↗

Modeling of dense reactive granular flows

Because of the lack of fundamental knowledge of both the fluid mechanical and thermal behavior of multiphase flow systems, the design of these processes has been necessarily based on empirical correlations and experiments performed in laboratory or pilot scale units. Detailed models and simulations of these systems would allow for the optimization of the processes and their scaling, without the need for expensive and time consuming testing. In this paper, we present the proposed hydrodynamic model in detail.

granular↗

Modeling of dense reactive granular flows

A general model of dense granular flows derived elsewhere is here assessed to understand the validity of some of the assumptions embedded in the derivation.

granular↗

Modelling of jet- and swirl-stabilized reacting flows in axisymmetric combustors

Turbulent nonreactive and reactive flows with and without swirl are analyzed, with particular attention given to the flow fields of a gas-fueled nonpremixed swirl-stabilized combustor and a premixed opposed-jet combustor. Local mean flow properties, including velocity, temperature, and major species concentrations, are calculated by solving numerically the governing partial differential equations with associated submodels for turbulence and combustion. The results of the study indicate that the constant-density k-epsilon turbulence model provides a satisfactory representation of the aerodynamics in most practical combustor flows. The exception is the case of jet-stabilized combustor flow, due to the fact that the k-epsilon model cannot replicate the highly dissipative phenomenon found in such flows.

Nikjooy, M.↗

Generation of boundary and boundary-layer fitting grids

The details of extended physical processes, such as the gas dynamic flow over an airfoil, the reactive flow through a combustor, or the electric field in a multi-contact transistor, are understood by solving the differential equations of a mathematical model of the process. The accuracy of finite difference methods for the numerical solution of the equations is increased if the underlying mesh fits the region boundaries and is closely spaced in regions where the solution is rapidly varying. Automatic methods for producing a satisfactorily adjusted mesh were developed for one dimensional problems. In one simple, effective scheme of this kind the unknown function and the distribution of mesh modes are found simultaneously, the nodes being placed so that they correspond to points uniformly spaced on the solution curve. In a two dimensional generalization, the nodes correspond to points equally spaced on the solution surface in two directions that are as nearly orthogonal as possible. Examples of such meshes are shown.

Ablow, C. M.↗

Forebody and base region real gas flow in severe planetary entry by a factored implicit numerical method. II - Equilibrium reactive gas

The factored-implicit, gas-dynamic algorithm has been adapted to the numerical simulation of equilibrium reactive flows. Changes required in the perfect gas version of the algorithm are developed, and the method of coupling gas-dynamic and chemistry variables is discussed. A flow-field solution that approximates a Jovian entry case was obtained by this method and compared with the same solution obtained by HYVIS, a computer program much used for the study of planetary entry. Comparison of surface pressure distribution and stagnation line shock-layer profiles indicates that the two solutions agree well.

Davy, W. C.↗

Selective Two-Photon Absorptive Resonance Femtosecond-Laser Electronic-Excitation Tagging (STARFLEET) Velocimetry in Flow and Combustion Diagnostics

Selective Two-Photon Absorptive Resonance Femtosecond-Laser Electronic-Excitation Tagging (STARFLEET), a non-seeded ultrafast-laser-based velocimetry technique, is demonstrated in reactive and non-reactive flows. STARFLEET is pumped via a two-photon resonance in N2 using 202.25-nm 100-fs light. STARFLEET greatly reduces the per-pulse energy required (30 μJ/pulse) to generate the signature FLEET emission compared to the conventional FLEET technique (1.1 mJ/pulse). This reduction in laser energy results in less energy deposited in the flow, which allows for reduced flow perturbations (reactive and non-reactive), increased thermometric accuracy, and less severe damage to materials. Velocity measurements conducted in a free jet of N2 and in a premixed flame show good agreement with theoretical velocities and further demonstrate the significantly less-intrusive nature of STARFLEET.

Jiang, Naibo↗

On the modelling of non-reactive and reactive turbulent combustor flows

A study of non-reactive and reactive axisymmetric combustor flows with and without swirl is presented. Closure of the Reynolds equations is achieved by three models: kappa-epsilon, algebraic stress and Reynolds stress closure. Performance of two locally nonequilibrium and one equilibrium algebraic stress models is analyzed assuming four pressure strain models. A comparison is also made of the performance of a high and a low Reynolds number model for combustor flow calculations using Reynolds stress closures. Effects of diffusion and pressure-strain models on these closures are also investigated. Two models for the scalar transport are presented. One employs the second-moment closure which solves the transport equations for the scalar fluxes, while the other solves the algebraic equations for the scalar fluxes. In addition, two cases of non-premixed and one case of premixed combustion are considered. Fast- and finite-rate chemistry models are applied to non-premixed combustion. Both show promise for application in gas turbine combustors. However, finite rate chemistry models need to be examined to establish a suitable coupling of the heat release effects on turbulence field and rate constants.

Nikjooy, Mohammad↗