Numerical Investigation of an Impinging Diffusion Flames-Effects of Fuel Variability

In our study, we investigate the differences between the combustion of different hydrocarbon fuels CH4, C3H8, C4H10. A numerical simulation of an impinging jet diffusion flames is used. The jet injector has a 10 mm in diameter and the distance between the jet flame and the vertical wall is 2 time half diameter. The fuel jet velocity was fixed for 11.8 m/s, corresponding to a Reynolds number of 6881. The flame characteristics varied from hydrocarbon to another for the same Reynolds number. The combustion products of CO, CO2, NO, OH, are depending on the methane and propane and butane flames for the same conditions. The temperature of the flame was varied from hydrocarbon to another the same as for the chemical species production rate. The concentration of the thermal and prompt NO pollutant depends on the temperature flow field and on the thermochemical characteristics of the hydrocarbon fuels.


Introduction
A diffusion flame is a very common practical system, such as a match flame or a candle flame, and is fundamental to many more complex systems.However, in a diffusion flame there is a relatively long time scale, allowing sufficient recombination reactions to take place such that it involves complex interaction of flow, transport and chemistry.There is more chemical and physical interaction in diffusion flames than in a premixed flame.The first systematic analysis of a confined, jet diffusion flame dates back to 1928 by Burke and Schumann (published in 1948) [1].Due to environmental concerns, the amount of NO x emission in combustion processed and its formation mechanism have always been of interest.Being of particular interest is source of nitric oxide, which is a major pollutant to the atmosphere.It has been postulated that there are two major mechanisms Zeldovich thermal and Fennimore prompt NO in the production of total NO.Barlow and Carter (1993) [2] reported a rather complete simultaneous measurement on temperature and species concentrations (NO, OH, and major species) in hydrogen flames.However, the amount of naturally generated NO in a gas phase hydrocarbon flame is very small, and most measurements were done with a certain amount of NO addition.Dong et al. [3] studied the characteristics of an impinging inverse diffusion flame jet the length of the jet impingement region, the characteristic of the wall gauge static pressure and the heat flux in the impingement region they provide a correlation between heat transfer performance and the hydrodynamic behavior, which facilitates the optimization of the jet impingement system design and operation.The aim is to understand the characterization of the extinction limits of fuel-air mixtures from low extinction strain rate methane-air flames to high extinction strain rate ethylene-air flame.Sarnacki et al. [4] reported an experimental and computational work, and they found that the variation of local extinction strain rate with changes in separation distance was within uncertainty of the experimental data.Zhen et al. [5] investigated the thermal and heat transfer behaviors of multifuel jet inverse diffusion flame with induce swirl and non-swirling under identical air/fuel rates.Flame appearances, temperature fields and wall static pressures were examined, and they found that the main reaction zone in the swirling flame is closer to the burner exit and the flame length is much shorter than the non-swirling flame and the wall static pressure and radial heat flux both are influenced by the swirl effect.A large eddy simulation technique was used to study the fuel variability on the dynamics of hydrogen and syngas impinging flames by Mira Martinez et al. [6] including a pure H 2 and (20% CO, 80% H 2 ), (40% CO, 60% H 2 ), (20% CO, 20% CO 2 , 60% H 2 ) the results of their study show that the flames develop vertical structures in the primary jet associated with the buoyancy and shear layer instability, and the wall jet progresses parallel to the impinging plate forming large scale vortex rings at different locations and strengths as a consequence of the fuel compositions.
Choe McDaid et al. [7] studied the effect of ignition location on the propagation of premixed and diffusion flames of hydrogen and mixtures of hydrogen and carbon dioxide towards.They found that the fuel velocity and Reynolds had a large effect on the observed flame velocities, and the type of the fuel affected the velocities and accelerations of the flame front.Experimental and numerical study had been conducted by Subhash et al. [8] to investigate the occurrence of off stagnation peak for laminar/air flame impinging on a flat surface of this off stagnation peak using a commercial CFD code Fluent.They found that this off stagnation in heat flux is primarily due to the peak in the axial velocity close to the impingement surface.Jaramillo et al. [9] applied the DNS and RANS techniques to study the fluid flow and heat transfer in plane impinging jets.The DNS results have been used as reference solution to assess the performance of several Reynolds averaged Navier Stokes (RANS) models.However all the models predict correctly the local Nusselt number at the stagnation region to investigate the complex flow field of an impinging jet Naseem et al. [10] used an LES simulation with a dynamic smagorinsky model.The large eddy simulation helps to understand the reason of occurrence of second peak in the radial distribution of Nusselt number at the target wall.They found that the LES simulation of a complex flow impinging jet is highly sensitive to the quality of the grid.
Zhen et al. [11] studied the emission of CO and NO x from swirling and non-swirling of an impinging inverse diffusion flames and they found that the parameters of air jet Reynolds numbers, overall equivalence ration and Nozzle to plate distance have significant influence on the overall pollutants emission.To reduce the NO x production a new down fired combustion technology based on multiple injection and multiple staging was developed by Min Kuang et al. [12].An experimental study was performed by Zhen et al. [13] to investigate the effects of the nozzle length on the air pollutant emission and noise radiation.They found that the noise radiation from the inner reaction cone of the flame is stronger than that from the lower and upper parts of the flame for the stoichiometric air/fuel ration.Zhen et al. [14] performed an experimental work to compare between the emission and impingement heat transfer behaviors liquefied petro-leum gas added of hydrogen and methane air flames (LPG-H 2 -air) and (CH 4 -H 2 -air) comparison shows a more significant change in the laminar burning speed, temperature and CO/NO x emissions in the CH 4 flames.Gurpreet et al. [15] investigated the heat transfer charactiristics of natural gas/air swirling flame impinging on a flat surface the dip in heat flux at and around stagnation point was observed in almost all cases which could be the main cause of non uniformity even in case of heating with swirling impinging flames.
Nadjib et al. [16] investigated the influence of the nitrogen dilution on the extinction of methane impinging diffusion flame, and they found that when the dilution rate increased the extinction of the diffusion flame increased.
In the current study, we investigated an impinging diffusion flame with three fuels, Methane, Propane and Butane with fixed fuel jet velocity.We also reported on the temperature response to the increase of the CH atom mole fraction, on the other hand we studied the relation between NO fraction production and heating high temperature for a non-premixed turbulent hydrocarbons jet flame situation.

Turbulent Governing Equations
In the present study, Fluent [20] (commercial CFD software) was used to model the flow field and heat transfer for diffusion turbulent methane/air flame impinging vertically on a flat surface with a reduced reaction mechanism (8 species are considered).
The general form of transport equations for two dimensional stationary turbulent reactive flows can be written as: where  denotes 1, , , , Fluent uses a control-volume-based technique to convert the governing equations to algebraic equations that can be solved numerically.
In order to resolve the turbulent flow problem we are used K-Epsilon turbulent RNG based model.

 
    Copyright © 2013 SciRes.OJFD The model constants appearing in the above equations are The effects of the mean strain rate and mean rotation on turbulent diffusion have been affected by using the renormalized RNG k- model Yakhot et al., 1992 [21], which employs equations of the same form as the standard k- model.The RNG k- model assumes different model coefficients evaluated by the renormalization group theory which vary with the ration of the turbulent to the mean strain, n, as described below: Where S is the strain rate modulus and the ratio k  is the turbulence time scale.

Modeling Non-Premixed Combustion
In non-premixed combustion, fuel and oxidizer enter the reaction zone in distinct streams.This is in contrast to premixed systems, in which reactants are mixed at the molecular level before burning.Examples of non-premixed combustion include methane combustion, pulverized coal furnaces, and diesel (compression) internalcombustion engines.
In order to resolve the turbulent chemistry interaction we are focused to use Pre PDF model based on the resolution of mean transport species equation and its variance 2 f  [20].
The species fractions considered in this investigation are   H . he impingement surface is parallel to ding vertically on the the fuel jet is 10 mm Fenimore "prompt" NO.In the numerical work, both Zeldovich thermal and Fenimore prompt formation of

Computational Domaine
Using a 2D model, t the fuel jet; and the jet was sprea impinging plate.The diameter of Figure 1 a total of seven transport equations (continuity, axial, and radial momentums, turbulence kinetic energy and its dissipation rate, energy, and radiative intensity), are solved using the commercial CFD package FLUENT [20].A second-order discretization scheme was used to solve all governing equations.Solution convergence was determined by two criteria.First is ensuring that the residuals of the solved equations drop below specified thresholds set at 10 −3 for all variables, while a residual of 10 −6 was used for the energy equation.The second convergence criterion is ensuring that the value of a sensitive property (e.g., concentration of a radical species) at a critical spatial location has stabilized and is no longer changing with iterations.

Numerical Results
In this simulation the temperature and concentration of major species C 4 H 10 , C 3 H 8 , CH 4 , H 2 , H 2 O, CO 2 , N 2 , and O 2 and minor species NO, CO, and OH was performed using Fluent software code.CH 4 , C 3 H 8 , C 4 H 10 , was used as the fuel in the impinging jet.Velocity magnitude, flame temperature and the structure of the turbulent flame region were compared with different fuel jet flames.The effect of N 2 dilution level, in the fuel stream, on the flame temperature is quite substantial with a drop from 1700 K in the case propane flame, to 1400 K in case of butane flame.The temperature at the centerline is 420 K and is the same for all cases.The mean temperature in the jet vicinity of the flow is 1200 K and is also consistent between all flames Figure 2 this result is compared with the experimental data in references [16][17][18][19] for the same authors.When we look in Figure 3 we observe that the mass fraction of CO 2 is maximal for the butane jet flame.The same results are obtained in Figure 4.
Figures 7 and 9 show how the calculated NO and CO 2 concentrations vary with the three jet fuels for the same simulation conditions.axial position (100 e and the mean peak NO is almost three times that of other flames.
The NO formation rate is significantly higher in the case of butane fuel compared with methane and propane fuels Figure 9

Conclusions
In this study, we investigated the characteristics of three hydrocarbons methane, propane and butane fuel flames.The temperature and species and NO production of turbulent jet flames under non-premixed conditions are investigated.
Though the modification of the above constants of the model, it is intended to simulate and control the modelling of the energy dissipation.The RNG k- model is a modification version of the standard k- turbulent model.It adopts a non-equilibrium strain parameter n kS   , C H ,CH ,O ,N ,H O,NO,CO ,CO,O In order to model nitric oxide formation in a flame, the ch ken into account.The reactions emical reactions involving nitrogen compounds must be ta

Figure 1 .
Figure 1.Schematic of mesh domain and boundary conditions for jet diffusion flame.

Figures 2 - 6
show contours plots for temperature and mass fractions of OH and CO 2 and OH and NO and production rate plotted versus mixture fraction for flames C 4 H 10 , C 3 H 8 , and CH 4 at the same axial location.

Figures 7 -
10 show radial profiles of temperature and mass fractions of OH and CO and NO and turbulent intensity for turbulent flames C 4 H 10 , C 3 H 8 , and CH 4 at an axial position of 100 mm above the jet exit.
due to the higher O atom and CH concentrations promoted by the fuel oxidation.Also shown in Figures 11-13 is the concentration of the NO production in the three flames zone on the same horizontal positions calculated from the turbulent flames modeling for the propane and butane flames more NO concentration near the impinging wall due to the effect of molecular weight on the mixing with the surrounding air.The mean NO mole fractions for flame CH 4 at the positions of 100 and 200 and 300 mm and 500 mm have the same peak of 2.5 ppm.For flames C H 8 and C 4 H 10 the of CH 4 flame.A n zone location appears for both flames.At locations 100 and 200 mm, flame el for flame C 4 H 10 , while at location 200 m, the peak level for both flames is similar.It is worth noting that the peak temperature is roughly the same for all locations and that the peak NO levels reported here for all flames are very small when compared with standard diffusion flames with similar Reynolds numbers.This indicates that dilution effects are presented in Figures 12 and 13.

Figure 7 .
Figure 7. Dioxide mass fraction for methane, propane and butane flames.

Figure 9 .
Figure 9. NO concentration profiles for methane, propane and butane flames.

Figure 11 .
Figure 11.Radial profiles of mean NO mass fraction at different radial locations for methane flame.

Figure 12 .
Figure 12.Radial profiles of mean NO mass fraction at different radial locations for propane flame.

-
The flame temperature varied with the fuel type with the increasing of CH atom.We obtain an increase in flame temperature level.From the numerical re of simultaneous Methane Propane Butane calculations, -We concluded that the turbulent flame propagation velocity varied from a fuel to another an ing in the fuel type will influence on the NO masse fraction produced.The same conclusion is obtained for the other species CO, CO 2 .