CMC simulations of lifted turbulent jet flame in a vitiated coflow

Lifted turbulent jet diffusion flame is simulated using Conditional Moment Closure (CMC). Specifically, the burner configuration of Cabra et al. [R. Cabra, T. Myhrvold, J.Y. Chen. R.W. Dibble, A.N. Karpetis, R.S. Barlow, Proc. Combust. Inst. 29 (2002) 1881-1887] is chosen to investigate H-2/N-2 jet flame supported by a vitiated coflow of products of lean H-2/air combustion. A 2D, axisymmetric flow-model fully coupled with the scalar fields, is employed. A detailed chemical kinetic scheme is included, and first order CIVIC is applied. Simulations are carried out for different jet velocities and coflow temperatures (T-c) The predicted liftoff generally agrees with experimental data, as well as joint-PDF results. Profiles of mean scalar fluxes in the mixture fraction space, for T-c = 1025 and 1080 K reveal that (1) Inside the flame zone, the chemical term balances the molecular diffusion term, and hence the Structure is of a diffusion flamelet for both cases. (2) In the pre-flame zone, the structure depends on the coflow temperature: for the 1025 K case, the chemical term being small, the advective term balances the axial turbulent diffusion term. However, for the 1080 K case. the chemical term is large and balances the advective term, the axial turbulent diffusion term being small. It is concluded that, lift-off is controlled (a) by turbulent premixed flame propagation for low coflow temperature while (b) by autoignition for high coflow temperature. (C) 2009 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Validation of a Modified Eddy Dissipation Concept Model for Stationary and Nonstationary Turbulent Diffusion Flames

A transient flame simulation tool based on unsteady Reynolds average Navier Stokes (RANS) is characterized for stationary and nonstationary flame applications with a motivation of performing computationally affordable flame stability studies. Specifically, the KIVA-3V code is utilized with incorporation of a recently proposed modified eddy dissipation concept for simulating turbulence-chemistry interaction along with a model for radiation loss. Detailed comparison of velocities, turbulent kinetic energies, temperature, and species are made with the experimental data of the turbulent, non-premixed DLR_A CH4/H-2/N-2 jet flame. The comparison shows that the model is able to predict flame structure very well. The effect of some of the modeling assumptions is assessed, and strategies to model a stationary diffusion flame are recommended. Unsteady flame simulation capabilities of the numerical model are assessed by simulating an acoustically excited, experimental, oscillatory H-2-air diffusion flame. Comparisons are made with oscillatory velocity field and OH plots, and the numerical code is observed to predict transient flame structure well.

Forcing of self-excited round jet diffusion flames

In this experimental and numerical study, two types of round jet are examined under acoustic forcing. The first is a non-reacting low density jet (density ratio 0.14). The second is a buoyant jet diffusion flame at a Reynolds number of 1100 (density ratio of unburnt fluids 0.5). Both jets have regions of strong absolute instability at their base and this causes them to exhibit strong self-excited bulging oscillations at welldefined natural frequencies. This study particularly focuses on the heat release of the jet diffusion flame, which oscillates at the same natural frequency as the bulging mode, due to the absolutely unstable shear layer just outside the flame. The jets are forced at several amplitudes around their natural frequencies. In the non-reacting jet, the frequency of the bulging oscillation locks into the forcing frequency relatively easily. In the jet diffusion flame, however, very large forcing amplitudes are required to make the heat release lock into the forcing frequency. Even at these high forcing amplitudes, the natural mode takes over again from the forced mode in the downstream region of the flow, where the perturbation is beginning to saturate non-linearly and where the heat release is high. This raises the possibility that, in a flame with large regions of absolute instability, the strong natural mode could saturate before the forced mode, weakening the coupling between heat release and incident pressure perturbations, hence weakening the feedback loop that causes combustion instability. © 2009 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Scalar and its dissipation in the near field of turbulent lifted jet flame

In this study various scalar dissipation rates and their modelling in the context of partially premixed flame are investigated. A DNS dataset of the near field of a turbulent hydrogen lifted jet flame is processed to analyse the mixture fraction and progress variable dissipation rates and their cross dissipation rate at several axial positions. It is found that the classical model for the passive scalar dissipation rate ε{lunate}̃ZZ gives good agreement with the DNS, while models developed based on premixed flames for the reactive scalar dissipation rate ε{lunate}̃cc only qualitatively capture the correct trend. The cross dissipation rate ε{lunate}̃cZ is mostly negative and can be reasonably approximated at downstream positions once ε{lunate}̃ZZ and ε{lunate}̃cc are known, although the sign cannot be determined. This approach gives better results than one employing a constant ratio of turbulent timescale and the scalar covariance c'Z'̃. The statistics of scalar gradients are further examined and lognormal distributions are shown to be very good approximations for the passive scalar and acceptable for the reactive scalar. The correlation between the two gradients increases downstream as the partially premixed flame in the near field evolves ultimately to a diffusion flame in the far field. A bivariate lognormal distribution is tested and found to be a reasonable approximation for the joint PDF of the two scalar gradients. © 2011 The Combustion Institute.

Numerical study on NOx/CO emissions in the diffusion flames of high-temperature off-gas of steelmaking converter

The combustion of high-temperature off-gas of steelmaking converter with periodical change of temperature and CO concentration always leads to CO and NOx over-standard emissions. In the paper, high-temperature off-gas combustion is simulated by adopting counterflow diffusion flame model, and some influencing factors of CO and NOx emissions are investigated by adopting a detailed chemistry GRI 3.0 mechanism. The emission index of NOx (EINOx) decreases 1.7–4.6% when air stoichiometric ratio (SR) increase from 0.6 to 1.4, and it dramatically increases with off-gas temperature at a given SR when the off-gas temperature is above 1500 K. High-concentration CO in off-gas can result in high NOx emissions, and NOx levels increase dramatically with CO concentration when off-gas temperature is above 1700 K. Both SR and off-gas temperature are important for the increase of CO burnout index (BICO) when SR is less than 1.0, but BICO increase about 1% when off-gas temperature increases from 1100 K to 1900 K at SR > 1.0. BICO increases with CO concentration in off-gas, and the influence of off-gas temperature on BICO is marginal. BICO increases with the relative humidity (RH) in air supplied, but it increases about 0.5% when RH is larger than 30%.

Soot Measurements in Steady and Pulsed Ethylene/Air Diffusion Flames Using Laser-Induced Incandescence

Combustion-generated carbon black nano particles, or soot, have both positive and negative effects depending on the application. From a positive point of view, it is used as a reinforcing agent in tires, black pigment in inks, and surface coatings. From a negative point of view, it affects performance and durability of many combustion systems, it is a major contributor of global warming, and it is linked to respiratory illness and cancer. Laser-Induced Incandescence (LII) was used in this study to measure soot volume fractions in four steady and twenty-eight pulsed ethylene diffusion flames burning at atmospheric pressure. A laminar coflow diffusion burner combined with a very-high-speed solenoid valve and control circuit provided unsteady flows by forcing the fuel flow with frequencies between 10 Hz and 200 Hz. Periodic flame oscillations were captured by two-dimensional phase-locked LII images and broadband luminosity images for eight phases (0°- 360°) covering each period. A comparison between the steady and pulsed flames and the effect of the pulsation frequency on soot volume fraction in the flame region and the post flame region are presented. The most significant effect of pulsing frequency was observed at 10 Hz. At this frequency, the flame with the lowest mean flow rate had 1.77 times enhancement in peak soot volume fraction and 1.2 times enhancement in total soot volume fraction; whereas the flame with the highest mean flow rate had no significant change in the peak soot volume fraction and 1.4 times reduction in the total soot volume fraction. A correlation (ƒv Reˉ1 = a+b· Str) for the total soot volume fraction in the flame region for the unsteady laminar ethylene flames was obtained for the pulsation frequency between 10 Hz and 200 Hz, and the Reynolds number between 37 and 55. The soot primary particle size in steady and unsteady flames was measured using the Time-Resolved Laser-Induced Incandescence (TIRE-LII) and the double-exponential fit method. At maximum frequency (200 Hz), the soot particles were smaller in size by 15% compared to the steady case in the flame with the highest mean flow rate.

On the ignition and extinction problems in forced convection systems

In this paper the problem of ignition and extinction has been formulated for the flow of a compressible fluid with Prandtl and Schmidt numbers taken as unity. In particular, the problems of (i) a jet impinging on a wall of combustible material and (ii) the opposed jet diffusion flame have been studied. In the wall jet case, three approximations in the momentum equation namely, (i) potential flow, (ii) viscous flow, (ii) viscous incompressible with k = 1 and (iii) Lees' approximation (taking pressure gradient terms zero) are studied. It is shown that the predictions of the mass flow rates at extinction are not very sensitive to the approximations made in the momentum equation. The effects of varying the wall temperature in the case (i) and the jet temperature in the case (ii) on the extinction speeds have been studied. The effects of varying the activation energy and the free stream oxidant concentration in case (ii), have also been investigated.

Sandwich propellant combustion: Modeling and experimental comparison

This paper reports reacting fluid dynamics calculations for an ammonium percholrate binder sandwich and extracts experimentally observed features including surface profiles and maximum regression rates as a function of pressure and binder thickness. These studies have been carried out by solving the two-dimensional unsteady Navier-Stokes equations with energy and species conservation equations and a kinetic model of three reaction steps (ammonium perchlorate decomposition flame, primary diffusion flame, and final diffusion flame) in the gas phase. The unsteady two-dimensional conduction equation is solved in the condensed phase. The regressing surface is unsteady and two dimensional. Computations have been carried out for a binder thickness range of 25-125 mum and a pressure range of 1.4 to 6.9 MPa. Good comparisons at several levels of detail are used to demonstrate the need for condensed-phase two-dimensional unsteady conduction and three-step gas-phase reactions. The choice of kinetic and thermodynamic parameters is crucial to good comparison with experiments. The choice of activation energy parameters for ammonium percholrate combustion has been made with stability of combustion in addition to experimentally determined values reported in literature. The choice of gas-phase parameters for the diffusion flames are made considering that (a) primary diffusion flame affects the low-pressure behavior and (b) final diffusion flame affects high-pressure behavior. The predictions include the low-pressure deflagration limit of the sandwich apart from others noted above. Finally, this study demonstrates the possibility of making meaningful comparisons with experimental observations on sandwich propellant combustion.

Progress in numerical modeling of non-premixed combustion

Progress is made on the numerical modeling of both laminar and turbulent non-premixed flames. Instead of solving the transport equations for the numerous species involved in the combustion process, the present study proposes reduced-order combustion models based on local flame structures.

For laminar non-premixed flames, curvature and multi-dimensional diffusion effects are found critical for the accurate prediction of sooting tendencies. A new numerical model based on modified flamelet equations is proposed. Sooting tendencies are calculated numerically using the proposed model for a wide range of species. These first numerically-computed sooting tendencies are in good agreement with experimental data. To further quantify curvature and multi-dimensional effects, a general flamelet formulation is derived mathematically. A budget analysis of the general flamelet equations is performed on an axisymmetric laminar diffusion flame. A new chemistry tabulation method based on the general flamelet formulation is proposed. This new tabulation method is applied to the same flame and demonstrates significant improvement compared to previous techniques.

For turbulent non-premixed flames, a new model to account for chemistry-turbulence interactions is proposed. %It is found that these interactions are not important for radicals and small species, but substantial for aromatic species. The validity of various existing flamelet-based chemistry tabulation methods is examined, and a new linear relaxation model is proposed for aromatic species. The proposed relaxation model is validated against full chemistry calculations. To further quantify the importance of aromatic chemistry-turbulence interactions, Large-Eddy Simulations (LES) have been performed on a turbulent sooting jet flame. %The aforementioned relaxation model is used to provide closure for the chemical source terms of transported aromatic species. The effects of turbulent unsteadiness on soot are highlighted by comparing the LES results with a separate LES using fully-tabulated chemistry. It is shown that turbulent unsteady effects are of critical importance for the accurate prediction of not only the inception locations, but also the magnitude and fluctuations of soot.

A novel theoretical approach to passive control of thermo-acoustic oscillations: Application to ducted heat sources

In this paper, we develop a linear technique that predicts how the stability of a thermo-acoustic system changes due to the action of a generic passive feedback device or a generic change in the base state. From this, one can calculate the passive device or base state change that most stabilizes the system. This theoretical framework, based on adjoint equations, is applied to two types of Rijke tube. The first contains an electrically-heated hot wire and the second contains a diffusion flame. Both heat sources are assumed to be compact so that the acoustic and heat release models can be decoupled. We find that the most effective passive control device is an adiabatic mesh placed at the downstream end of the Rijke tube. We also investigate the effects of a second hot wire and a local variation of the cross-sectional area but find that both affect the frequency more than the growth rate. This application of adjoint sensitivity analysis opens up new possibilities for the passive control of thermo-acoustic oscillations. For example, the influence of base state changes can be combined with other constraints, such as that the total heat release rate remains constant, in order to show how an unstable thermo-acoustic system should be changed in order to make it stable. Copyright © 2013 by ASME.

Estudo experimental de chamas difusivas livres turbulentas de gás natural submetidas a oscilações acústicas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

The internal structure of hydrogen-air difussion flames

The Internal Structure of Hydrogen-Air Diffusion Flames. Tho purpose of this paper is to study finite rate chemistry effects in diffusion controlled hydrogenair flames undor conditions appearing in some cases in a supersonic combustor. Since for large reaction rates the flame is close to chemical equilibrium, the reaction takes place in a very thin region, so thata "singular perturbation "treatment" of the problem seems appropriate. It has been shown previously that, within the inner or reaction zone, convection effects may be neglocted, the temperature is constant across the flame, and tho mass fraction distributions are given by ordinary differential equations, whore tho only independent variable involved is tho coordinate normal to the flame surface. Tho solution of the outer problom, which is a pure mixing problem with the additional condition that fuol and oxidizer do not coexist in any zone, provides t h e following information: tho flame position, rates of fuel consumption, temperature, concentrators of species, fluid velocity outside of tho flame, and the boundary conditions required to solve the "inner problem." The main contribution of this paper consists in the introduction of a fairly complicated chemical kinetic scheme representing hydrogen-oxygen reaction. The nonlinear equations expressing the conservation of chemical species are approximately integrated by means of an integral method. It has boen found that, in the case considered of a near-equilibrium diffusion flame, tho role played by the dissociation-recombination reactions is purely marginal, and that somo of the second order "shuffling" reactions are close to equilibrium. The method shown here may be applied to compute the distanco from the injector corresponding to a given separation from equilibrium, say ten to twenty percent. For the casos whore this length is a small fraction of the combustion zone length, the equilibrium treatment describes properly tho flame behavior.

Dynamics of thermal ignition of spray flames in mixing layers

Conditions are identified under which analyses of laminar mixing layers can shed light on aspects of turbulent spray combustion. With this in mind, laminar spray-combustion models are formulated for both non-premixed and partially premixed systems. The laminar mixing layer separating a hot-air stream from a monodisperse spray carried by either an inert gas or air is investigated numerically and analytically in an effort to increase understanding of the ignition process leading to stabilization of high-speed spray combustion. The problem is formulated in an Eulerian framework, with the conservation equations written in the boundary-layer approximation and with a one-step Arrhenius model adopted for the chemistry description. The numerical integrations unveil two different types of ignition behaviour depending on the fuel availability in the reaction kernel, which in turn depends on the rates of droplet vaporization and fuel-vapour diffusion. When sufficient fuel is available near the hot boundary, as occurs when the thermochemical properties of heptane are employed for the fuel in the integrations, combustion is established through a precipitous temperature increase at a well-defined thermal-runaway location, a phenomenon that is amenable to a theoretical analysis based on activation-energy asymptotics, presented here, following earlier ideas developed in describing unsteady gaseous ignition in mixing layers. By way of contrast, when the amount of fuel vapour reaching the hot boundary is small, as is observed in the computations employing the thermochemical properties of methanol, the incipient chemical reaction gives rise to a slowly developing lean deflagration that consumes the available fuel as it propagates across the mixing layer towards the spray. The flame structure that develops downstream from the ignition point depends on the fuel considered and also on the spray carrier gas, with fuel sprays carried by air displaying either a lean deflagration bounding a region of distributed reaction or a distinct double-flame structure with a rich premixed flame on the spray side and a diffusion flame on the air side. Results are calculated for the distributions of mixture fraction and scalar dissipation rate across the mixing layer that reveal complexities that serve to identify differences between spray-flamelet and gaseous-flamelet problems.

Measurement and Simulation of Suppression Effects in a Buoyant Turbulent Line Fire

An experimental and numerical study of turbulent fire suppression is presented. For this work, a novel and canonical facility has been developed, featuring a buoyant, turbulent, methane or propane-fueled diffusion flame suppressed via either nitrogen dilution of the oxidizer or application of a fine water mist. Flames are stabilized on a slot burner surrounded by a co-flowing oxidizer, which allows controlled delivery of either suppressant to achieve a range of conditions from complete combustion through partial and total flame quenching. A minimal supply of pure oxygen is optionally applied along the burner to provide a strengthened flame base that resists liftoff extinction and permits the study of substantially weakened turbulent flames. The carefully designed facility features well-characterized inlet and boundary conditions that are especially amenable to numerical simulation. Non-intrusive diagnostics provide detailed measurements of suppression behavior, yielding insight into the governing suppression processes, and aiding the development and validation of advanced suppression models. Diagnostics include oxidizer composition analysis to determine suppression potential, flame imaging to quantify visible flame structure, luminous and radiative emissions measurements to assess sooting propensity and heat losses, and species-based calorimetry to evaluate global heat release and combustion efficiency. The studied flames experience notable suppression effects, including transition in color from bright yellow to dim blue, expansion in flame height and structural intermittency, and reduction in radiative heat emissions. Still, measurements indicate that the combustion efficiency remains close to unity, and only near the extinction limit do the flames experience an abrupt transition from nearly complete combustion to total extinguishment. Measurements are compared with large eddy simulation results obtained using the Fire Dynamics Simulator, an open-source computational fluid dynamics software package. Comparisons of experimental and simulated results are used to evaluate the performance of available models in predicting fire suppression. Simulations in the present configuration highlight the issue of spurious reignition that is permitted by the classical eddy-dissipation concept for modeling turbulent combustion. To address this issue, simple treatments to prevent spurious reignition are developed and implemented. Simulations incorporating these treatments are shown to produce excellent agreement with the experimentally measured data, including the global combustion efficiency.

Regimes of spray vaporization and combustion in counterflow configurations

This article addresses the problem of spray vaporization and combustion in axisymmetric opposed-jet configurations involving a stream of hot air counterflowing against a stream of nitrogen carrying a spray of fuel droplets. The Reynolds numbers of the jets are assumed to be large, so that mixing of the two streams is restricted to a thin mixing layer that separates the counterflowing streams. The evolution of the droplets in their feed stream from the injection location is seen to depend fundamentally on the value of the droplet Stokes number, St, defined as the ratio of the droplet acceleration time to the mixing layer strain time close to the stagnation point. Two different regimes of spray vaporization and combustion can be identified depending on the value of St. For values of St below a critical value, equal to 1/4 for dilute sprays with small values of the spray liquid mass loading ratio, the droplets decelerate to approach the gas stagnation plane with a vanishing axial velocity. In this case, the droplets located initially near the axis reach the mixing layer, where they can vaporize due to the heat received from the hot air, producing fuel vapor that can burn with the oxygen in a diffusion flame located on the air side of the mixing layer. The character of the spray combustion is different for values of St of order unity, because the droplets cross the stagnation plane and move into the opposing air stream, reaching distances that are much larger than the mixing layer thickness before they turn around. The vaporization of these crossing droplets, and also the combustion of the fuel vapor generated by them, occur in the hot air stream, without significant effects of molecular diffusion, generating a vaporization-assisted nonpremixed flame that stands on the air side outside the mixing layer. Separate formulations will be given below for these two regimes of combustion, with attention restricted to the near-stagnation-point region, where the solution is self-similar and all variables are only dependent on the distance to the stagnation plane. The resulting formulations display a reduced number of controlling parameters that effectively embody dependences of the structure of the spray flame on spray dilution, droplet inertia, and fuel preferential diffusion. Sample solutions are given for the limiting cases of pure vaporization and of infinitely fast chemistry, with the latter limit formulated in terms of chemistry-free coupling functions that allow for general nonunity Lewis numbers of the fuel vapor.