The standard κ-ϵ equations and other turbulence models are evaluated with respect to their applicability in swirling, recirculating flows. The turbulence models are formulated on the basis of two separate viewpoints. The first perspective assumes that an isotropic eddy viscosity and the modified Boussinesq hypothesis adequately describe the stress distributions, and that the source of predictive error is a consequence of the modeled terms in the κ-ϵ equations. Both stabilizing and destabilizing Richardson number corrections are incorporated to investigate this line of reasoning. A second viewpoint proposes that the eddy viscosity approach is inherently inadequate and that a redistribution of the stress magnitudes is necessary. Investigation of higher-order closure is pursued on the level of an algebraic stress closure. Various turbulence model predictions are compared with experimental data from a variety of isothermal, confined studies. Supportive swirl comparisons are also performed for a laminar flow case, as well as reacting flow cases. Parallel predictions or contributions from other sources are also consulted where appropriate. Predictive accuracy was found to be a partial function of inlet boundary conditions and numerical diffusion. Despite prediction sensitivity to inlet conditions and numerics, the data comparisons delineate the relative advantages and disadvantages of the various modifications. Possible research avenues in the area of computational modeling of strongly swirling, recirculating flows are reviewed and discussed.

Bubbles, Drops, and Particles

PDF methods for turbulent reactive flows

Laminar diffusion flamelet models in non-premixed turbulent combustion

The laminar flamelet concept covers a regime in turbulent combustion where chemistry (as compared to transport processes) is fast such that it occurs in asymptotically thin layers—called flamelets—embedded within the turbulent flow field. This situation occurs in most practical combustion systems including reciprocating engines and gas turbine combustors. The inner structure of the flamelets is one-dimensional and time dependent. This is shown by an asymptotic expansion for the Damkohler number of the rate determining reaction which is assumed to be large. Other non-dimensional chemical parameters such as the nondimensional activation energy or Zeldovich number may also be large and may be related to the Damkohler number by a distinguished asymptoiic limit. Examples of the flamelet structure are presented using onestep model kinetics or a reduced four-step quasi-global mechanism for methane flames. For non-premixed combustion a formal coordinate transformation using the mixture fraction Z as independent variable leads to a universal description. The instantaneous scalar dissipation rate χ of the conserved scalar Z is identified to represent the diffusion time scale that is compared with the chemical time scale in the definition of the Damkohler number. Flame stretch increases the scalar dissipation rate in a turbulent flow field. If it exceeds a critical value χ q the diffusion flamelet will extinguish. Considering the probability density distribution of χ , it is shown how local extinction reduces the number of burnable flamelets and thereby the mean reaction rate. Furthermore, local extinction events may interrupt the connection to burnable flamelets which are not yet reached by an ignition source and will therefore not be ignited. This phenomenon, described by percolation theory, is used to derive criteria for the stability of lifted flames. It is shown how values of ∋ q obtained from laminar experiments scale with turbulent residence times to describe lift-off of turbulent jet diffusion flames. For non-premixed combustion it is concluded that the outer mixing field—by imposing the scalar dissipation rate—dominates the flamelet behaviour because the flamelet is attached to the surface of stoichiometric mixture. The flamelet response may be two-fold: burning or non-burning quasi-stationary states. This is the reason why classical turbulence models readily can be used in the flamelet regime of non-premixed combustion. The extent to which burnable yet non-burning flamelets and unsteady transition events contribute to the overall statistics in turbulent non-premixed flames needs still to be explored further. For premixed combustion the interaction between flamelets and the outer flow is much stronger because the flame front can propagate normal to itself. The chemical time scale and the thermal diffusivity determine the flame thickness and the flame velocity. The flamelet concept is valid if the flame thickness is smaller than the smallest length scale in the turbulent flow, the Kolmogorov scale. Also, if the turbulence intensity v′ is larger than the laminar flame velocity, there is a local interaction between the flame front and the turbulent flow which corrugates the front. A new length scale L G =v F 3 /∈ , the Gibson scale, is introduced which describes the smaller size of the burnt gas pockets of the front. Here v F is the laminar flame velocity and ∈ the dissipation of turbulent kinetic energy in the oncoming flow. Eddies smaller than L G cannot corrugate the flame front due to their smaller circumferential velocity while larger eddies up to the macro length scale will only convect the front within the flow field. Flame stretch effects are the most efficient at the smallest scale L G . If stretch combined with differential diffusion of temperature and the deficient reactant, represented by a Lewis number different from unity, is imposed on the flamelet, its inner structure will respond leading to a change in flame velocity and in some cases to extinction. Transient effects of this response are much more important than for diffusion flamelets. A new mechanism of premixed flamelet extinction, based on the diffusion of radicals out of the reaction zone, is described by Rogg. Recent progress in the Bray-Moss-Libby formulation and the pdf-transport equation approach by Pope are presented. Finally, different approaches to predict the turbulent flame velocity including an argument based on the fractal dimension of the flame front are discussed.

Laminar Flamelet Concepts in Turbulent Combustion

A new approach to chemistry modelling for large-eddy simulation of turbulent reacting flows is developed. Instead of solving transport equations for all of the numerous species in a typical chemical mechanism and modelling the unclosed chemical source terms, the present study adopts an indirect mapping approach, whereby all of the detailed chemical processes are mapped to a reduced system of tracking scalars. Here, only two such scalars are considered: a mixture fraction variable, which tracks the mixing of fuel and oxidizer, and a progress variable, which tracks the global extent of reaction of the local mixture. The mapping functions, which describe all of the detailed chemical processes with respect to the tracking variables, are determined by solving quasi-steady diffusion-reaction equations with complex chemical kinetics and multicomponent mass diffusion. The performance of the new model is compared to fast-chemistry and steady-flamelet models for predicting velocity, species concentration, and temperature fields in a methane-fuelled coaxial jet combustor for which experimental data are available. The progress-variable approach is able to capture the unsteady, lifted flame dynamics observed in the experiment, and to obtain good agreement with the experimental data, while the fast-chemistry and steady-flamelet models both predict an attached flame.

Progress-variable approach for large-eddy simulation of non-premixed turbulent combustion

http://staff.mechmining.uq.edu.au/klimenko/pub/pdf/Prog_Energy_Combust_Sci_1999_p.pdf

Conditional moment closure for turbulent combustion

Equations governing the variation through the flow field of averages of quantities such as species mass fractions, conditional on mixture fraction, are derived and modeled. The conditioning variable adds to the independent dimensions of the problem, but it is found that this is offset in some cases by reduction in the spatial dimensionality needed. Predictions are made for the reacting scalar mixing layer, and these show good agreement with experiment. The methodology effectively decouples the kinetics from the large inhomogeneity or macromixing aspects of the flow while preserving the input from the scalar dissipation or micromixing. Arbitrarily complex kinetics may be used within reasonable computational cost.

Conditional moment closure for turbulent reacting flow

Turbulent jet diffusion flames

Flamelet theories are examined in the context of turbulent nonpremixed combustion. The criterion requiring that the reaction zone be thinner than the Kolmogoroff length scale is converted to one in which the range of mixture fraction over which reaction occurs is compared with the scalar scale for the small scale fluctuations. It is found that the criterion is violated in many of the flames of interest, particularly away from the nozzle in jet flames and in recirculating flows. Laser imaging data for scalar mixing is used to illustrate the structure of quasi-equilibrium distributed reaction (QEDR) flames at high Damkohler number. The structure is found to tend toward that of the scalar dissipation. The apparent success of flamelet theories in predicting mean concentrations of CO in hydrocarbon flames is ascribed to the metastable equilibrium of rich mixtures arising from chain terminating reactions of radicals with the fuel. Differences in flame structure for piloted jet diffusion flames near extinction of methane and CO/H 2 /N 2 fuels are discussed in terms of these fluid dynamic and kinetic insights.

The structure of turbulent nonpremixed flames

The theory for the mixing and chemical reaction of two streams of fluid is developed for unsteady laminar and for turbulent flow. Only one conserved scalar is required to fully describe the mixing and various reaction models are considered including the Burke-Schumann flame sheet and shifting equilibrium. A new expression is derived for the instantaneous reaction rate of any species Wi = –pD(▿ξ)2Yi/dξ2 where ξ is a conserved scalar. New insight is obtained into the structure of the instantaneous reaction zone for both laminar and turbulent flames. Diagnostic parameters are presented for determining when various models are applicable and for what properties.

Robert W. Bilger

Papers

Conditional moment closure for turbulent combustion

Conditional moment closure for turbulent reacting flow

Turbulent jet diffusion flames

The structure of turbulent nonpremixed flames

The Structure of Diffusion Flames