Search NASA⌕ Search

Engineering topics

Cline, M. C.

Publications and source records attributed to Cline, M. C..

Analysis of Fuel Vaporization, Fuel-Air Mixing, and Combustion in Integrated Mixer-Flame Holders

Requirements to limit pollutant emissions from the gas turbine engines for the future High-Speed Civil Transport (HSCT) have led to consideration of various low-emission combustor concepts. One such concept is the Integrated Mixer-Flame Holder (IMFH). This report describes a series of IMFH analyses performed with KIVA-II, a multi-dimensional CFD code for problems involving sprays, turbulence, and combustion. To meet the needs of this study, KIVA-II's boundary condition and chemistry treatments are modified. The study itself examines the relationships between fuel vaporization, fuel-air mixing, and combustion. Parameters being considered include: mixer tube diameter, mixer tube length, mixer tube geometry (converging-diverging versus straight walls), air inlet velocity, air inlet swirl angle, secondary air injection (dilution holes), fuel injection velocity, fuel injection angle, number of fuel injection ports, fuel spray cone angle, and fuel droplet size. Cases are run with and without combustion to examine the variations in fuel-air mixing and potential for flashback due to the above parameters. The degree of fuel-air mixing is judged by comparing average, minimum, and maximum fuel/air ratios at the exit of the mixer tube, while flame stability is monitored by following the location of the flame front as the solution progresses from ignition to steady state. Results indicate that fuel-air mixing can be enhanced by a variety of means, the best being a combination of air inlet swirl and a converging-diverging mixer tube geometry. With the IMFH configuration utilized in the present study, flashback becomes more common as the mixer tube diameter is increased and is instigated by disturbances associated with the dilution hole flow.

Deur, J. M.↗

Analysis of Lean Premixed/Prevaporized Combustion with KIVA-2

Requirements to reduce the emissions of pollutants from gas turbines used in aircraft propulsion and ground based power generation have led to consideration of lean premixed/prevaporized (LPP) combustion concept. This paper describes some of the LPP flame tube analyses performed at the NASA Research Center with KIVA-2, a well-known multi-dimensional CFD code for problems including sprays, turbulence, and combustion. Modifications to KIVA-2's boundary condition and chemistry treatments have been made to meet the needs of the present study. The study itself focuses on two key aspects of the LPP concept, low emissions and flame stability (including flashback and lean blowoff.

Deur, J. M.↗

A three-dimensional algebraic grid generation scheme for gas turbine combustors with inclined slots

A 3D algebraic grid generation scheme is presented for generating the grid points inside gas turbine combustors with inclined slots. The scheme is based on the 2D transfinite interpolation method. Since the scheme is a 2D approach, it is very efficient and can easily be extended to gas turbine combustors with either dilution hole or slot configurations. To demonstrate the feasibility and the usefulness of the technique, a numerical study of the quick-quench/lean-combustion (QQ/LC) zones of a staged turbine combustor is given. Preliminary results illustrate some of the major features of the flow and temperature fields in the QQ/LC zones. Formation of co- and counter-rotating bulk flow and shape temperature fields can be observed clearly, and the resulting patterns are consistent with experimental observations typical of the confined slanted jet-in-cross flow. Numerical solutions show the method to be an efficient and reliable tool for generating computational grids for analyzing gas turbine combustors with slanted slots.

Yang, S. L.↗

Numerical simulation of a low-emission gas turbine combustor using KIVA-II

A modified version of the KIVA-II code was used to obtain a multidimensional numerical solution for the turbulent two-phase chemically reacting flows inside a staged turbine combustor (STC). The STC under consideration is equipped with an advanced airblast fuel nozzle and encompasses a fuel nozzle (FN), a rich-burn (RB) zone, a converging connecting pipe, a quick-quench (QQ) zone, a diverging connecting pipe, and a lean-combustion (LC) zone. The STC was divided into two subsystems, namely, FN/RB zone and QQ/LC zones, and the numerical solutions were obtained separately for each subsystem. Preliminary data characterize the major features of the flow and temperature fields inside the STC. Information on velocity, temperature, and some critical species in the FN/RB zone is presented. In the QQ/LC zones, formation of the co- and counter-rotating bulk flow and the sandwiched-ring-shaped temperature field can be clearly seen. The calculations of the mass-weighted standard deviation and the pattern factor of temperature indicated that the mixing performance of the STC is very promising.

Yang, S. L.↗

Computation of the Space Shuttle solid rocket booster nozzle start-up transient flow

The first NASA Space Shuttle flight (STS-1) produced an overpressure wave that exceeded preflight predictions by as much as 5 to 1. This second overpressure wave occurred just after the solid rocket booster (SRB) igniter wave. To understand this overpressure phenomenon, a numerical simulation effort was undertaken. Both the SRB static firing test and STS-1 geometries were studied for two-dimensional, inviscid and viscous flow. The inviscid calculations did not produce significant second overpressure waves. However, the viscous calculations did produce second overpressure waves that qualitatively agree with experiment. These overpressure waves were present in both the static firing test and STS-1 geometries. This second overpressure wave is generated by the motion of the boundary layer separation point and the subsequent radial motion of the exhaust jet during the start-up of the SRB nozzle flow. The presence of the mobile launch platform exhaust hole wale amplifies this wave, but does not appear to be the source of any additional overpressure waves. The lack of good quantitative agreement between theory and experiment indicates that other overpressure sources, not accounted for by this simulation, may be present.

Cline, M. C.↗

Computation of high Reynolds number internal/external flows

A general, user oriented computer program, called VNAF2, developed to calculate high Reynolds number internal/external flows is described. The program solves the two dimensional, time dependent Navier-Stokes equations. Turbulence is modeled with either a mixing length, a one transport equation, or a two transport equation model. Interior grid points are computed using the explicit MacCormack scheme with special procedures to speed up the calculation in the fine grid. All boundary conditions are calculated using a reference plane characteristic scheme with the viscous terms treated as source terms. Several internal, external, and internal/external flow calculations are presented.

Cline, M. C.↗

VNAP2: A Computer Program for Computation of Two-dimensional, Time-dependent, Compressible, Turbulent Flow

A computer program, VNAP2, for calculating turbulent (as well as laminar and inviscid), steady, and unsteady flow is presented. It solves the two dimensional, time dependent, compressible Navier-Stokes equations. The turbulence is modeled with either an algebraic mixing length model, a one equation model, or the Jones-Launder two equation model. The geometry may be a single or a dual flowing stream. The interior grid points are computed using the unsplit MacCormack scheme. Two options to speed up the calculations for high Reynolds number flows are included. The boundary grid points are computed using a reference plane characteristic scheme with the viscous terms treated as source functions. An explicit artificial viscosity is included for shock computations. The fluid is assumed to be a perfect gas. The flow boundaries may be arbitrary curved solid walls, inflow/outflow boundaries, or free jet envelopes. Typical problems that can be solved concern nozzles, inlets, jet powered afterbodies, airfoils, and free jet expansions. The accuracy and efficiency of the program are shown by calculations of several inviscid and turbulent flows. The program and its use are described completely, and six sample cases and a code listing are included.

Cline, M. C.↗

Computation of high Reynolds number internal/external flows

A general, user oriented computer program, called VNAP2, has been developed to calculate high Reynolds number, internal/external flows. VNAP2 solves the two-dimensional, time-dependent Navier-Stokes equations. The turbulence is modeled with either a mixing-length, a one transport equation, or a two transport equation model. Interior grid points are computed using the explicit MacCormack scheme with special procedures to speed up the calculation in the fine grid. All boundary conditions are calculated using a reference plane characteristic scheme with the viscous terms treated as source terms. Several internal, and internal/external flow calculations are presented.

Cline, M. C.↗

Computation of high Reynolds number internal/external flows

A general, user oriented computer program, called VNAP2, was developed to calculate high Reynolds number, internal/ external flows. The VNAP2 program solves the two dimensional, time dependent Navier-Stokes equations. The turbulence is modeled with either a mixing-length, a one transport equation, or a two transport equation model. Interior grid points are computed using the explicit MacCormack Scheme with special procedures to speed up the calculation in the fine grid. All boundary conditions are calculated using a reference plane characteristic scheme with the viscous terms treated as source terms. Several internal, external, and internal/external flow calculations are presented.

Cline, M. C.↗

Stability aspects of diverging subsonic flow

The paper describes the computation of two-dimensional, subsonic, diverging internal flows and how they differ from the corresponding converging flows. Such diverging or decelerating flows occur in such obvious places as subsonic diffusers and inlets; however, such flows also occur in supersonic nozzles in the presence of a normal shock. The flow instability and its relation to the numerical method used, boundary conditions, and viscous effects are assessed both analytically and numerically. The inviscid flow is shown to be physically unstable and a poor representation of the true viscous flow.

Cline, M. C.↗

The computation of steady nozzle flow by a time-dependent method

The steady flow in two-dimensional and axisymmetric nozzles was computed using a time-dependent method. In this method the interior mesh points were computed using the MacCormack finite-difference scheme, while a characteristic scheme was used to calculate the boundary mesh points. No explicit artificial viscosity term was included. The fluid was assumed to be a perfect gas. This method was used to compute the flow in a 45 deg - 15 deg conical, converging-diverging nozzle, a 15 deg conical, converging nozzle, and a 10 deg conical, plug nozzle. Good agreement between the numerical solution and experimental data was found. In contrast to previous time-dependent methods, the computational times were less than one minute on a CDC 6600 computer.

Cline, M. C.↗

Computation of steady nozzle flow by a time-dependent method

The equations of motion governing steady, inviscid flow are of a mixed type, that is, hyperbolic in the supersonic region and elliptic in the subsonic region. These mathematical difficulties may be removed by using the so-called time-dependent method, where the governing equations become hyperbolic everywhere. The steady-state solution may be obtained as the asymptotic solution for large time. The object of this research was to develop a production type computer program capable of solving converging, converging-diverging, and plug two-dimensional nozzle flows in computational times of 1 min or less on a CDC 6600 computer.

Cline, M. C.↗

The analysis of nonequilibrium, chemically reacting, supersonic flow in three dimensions using a bicharacteristics method.

A second-order numerical method of characteristics based on a bicharacteristics scheme was developed for the evaluation of steady, supersonic, nonequilibrium, chemically reacting flows. The absolute accuracy and order of accuracy of the method was shown by comparisons with spherical source flow and axisymmetric, nonequilibrium flows. Numerical results were obtained for elliptical and super-elliptical nozzles. These results illustrate the complex nature of three-dimensional flows and the inadequacy of quasi-three-dimensional methods which neglect cross flows.

Cline, M. C.↗

Comparison of characteristic schemes for three-dimensional, steady, isentropic flow.

The accuracy and order of three bicharacteristic schemes based on the numerical method of characteristics for three-dimensional, steady, isentropic, supersonic flow were compared based on numerical calculations. Comparisons were made for source flows and for flows in contoured axisymmetric nozzles. The first two methods were tetrahedral line networks, differing only in the manner in which the cross derivatives were evaluated. The third method was the pentahedral line network proposed by Butler. All three methods were numerically stable and produced efficient and accurate solutions. A discussion of the relative merits and disadvantages of each method is presented.

Cline, M. C.↗