Microsoft Word - ETASR_V11_N2_pp6943-6948 Engineering, Technology & Applied Science Research Vol. 11, No. 2, 2021, 6943-6948 6943 www.etasr.com Bouziane et al.: Investigation of Swirl Stabilized CH4 Air Flame with Varied Hydrogen Content by using … Investigation of Swirl Stabilized CH4 Air Flame with Varied Hydrogen Content by using Computational Fluid Dynamics (CFD) to Study the Temperature Field and Flame Shape A. Bouziane Laboratory of Materials and Reactive Systems (LMSR) Faculty of Technology University of Djillali Liabes Sidi Bel Abbes, Algeria abouziane@yahoo.fr A. Alami Renewable Energies Department Saad Dahlab Blida University 1 Blida, Algeria alami.ahmed@univ-blida.dz M. Zaitri Laboratory of Gas Combustion and Environment University of Sciences and Technology of Oran Oran, Algeria zaitri_mohammed@yahoo.fr B. Bouchame Laboratory of Mechanics of Structures and Solids (LMSS) Faculty of Technology University of Djillali Liabes Sidi Bel Abbes, Algeria boucham1@yahoo.com M. Bouchetara Laboratory of Gas Combustion and Environment University of Sciences and Technology of Oran Oran, Algeria mbouchetara@hotmail.com Abstract-In the current paper, numerical simulations of the combustion of turbulent CH4-H2 are presented employing the standard k-epsilon and the RNG k-epsilon for turbulence closure. The Fr-ED concept is carried out to account for chemistry/ turbulence interaction. The hydrogen content is varied in the fuel stream from 0% to 100%. The numerical solutions are validated by comparison with corresponding experimental data from the Combustion Laboratory of the University of Milan. The flow is directed radially outward. This method of fuel injection has been already been explored experimentally. The results show that the structure of the flame is described reasonably and both standard k-ɛ and RNG k- ɛ models can predict the flame shape. The general aspect of the temperature profiles is well predicted. The temperature profiles are indicating a different trend between CH4 and CH4/H2 fuel mixtures. Keywords-RNG k-epsilon; swirl; hydrogen; CFD; CH4 I. INTRODUCTION Lean combustion emerges as a potential and an advantageous concept to attain efficiently the combustion process and control NOx emissions in power-generation gas turbines [1–3]. The best approach to obtain a stable lean flame is to add a highly reactive fuel to a traditional hydrocarbon fuel. Due to its high chemical reactivity, hydrogen could be used as a useful fuel addition to enhance the combustion of natural gas. The use of hydrogen in blended fuel can be a good option regarding combustion characteristics, such as extending the lean flammability limits [4, 5] and the environmental advantages by decreasing the emissions [6]. The research on the combustion of NG–H2 fuel mixtures is open and growing. Studies [5, 7-9] on NG–H2 hybrid flame stability behavior are at the center of attention. Authors in [10-11] conducted an experimental investigation on laminar LPG-H2 jet diffusion flame and NG-H2 turbulent flame respectively and affirmed that a very small fluctuation in the visible flame length, in up to 20% of hydrogen addition, enhances the flame length when higher hydrogen concentrations are injected. In the other hand, different investigations [5, 7, 8, 12, 13] were undertaken on the height of the blue cone and pointed out the influence of hydrogen injection in the increment of the combustion rate, the reduction of CO, CO2 and HC emissions, and the extension of the stable lean limit of natural gas burner [1, 9-12, 14-16]. Authors in [17] developed an experimental investigation in optically accessible propane/ hydrogen burner and found that Corresponding author: A. Bouziane Engineering, Technology & Applied Science Research Vol. 11, No. 2, 2021, 6943-6948 6944 www.etasr.com Bouziane et al.: Investigation of Swirl Stabilized CH4 Air Flame with Varied Hydrogen Content by using … the flame pattern of the hydrogen is of a weak turbulent flame [18]. A significant reduction in both CO emissions and soot concentration with hydrogen addition was reported in [4, 11, 12] when turbulent jet diffusion flame was characterized. Authors in [19-21] report an increase in flame adiabatic temperature and a reduction in flame thickness. Authors in [15] investigated the counterflow of NG/H2–air partially premixed flame numerically and experimentally. They stated that the influence of prompt NOx was important in the total NOx formation. One of the purposes of this paper is to support certain theories and assumptions drawn from experiments in swirl flames. The goal was to obtain qualitatively correct predictions. This work focuses on the numerical simulation of turbulent flames utilizing well known turbulence and combustion models [22, 23]. Leaner flames due to the hydrogen addition with different various hydrogen-hydrocarbon mixtures lead to thermo- diffusive instabilities which have an impact on flame propagation and stability [24, 25]. More details are needed about the interaction of these types of flames with turbulence in practical devices (swirl-stabilized burner). A better understanding of the fuel improvement on the fundamental flame properties is critical to the development of combustion systems and the role of hydrogen on these systems appears to be an important parameter. The goal of the current study is the investigation of flame behavior by the use of numerical simulations with the RANS model to treat turbulence coupling (the standard /RNG k-epsilon ) and the Fr-ED model to represent the turbulence/chemistry interactions. For this purpose, the effect of hydrogen addition in methane flames has been examined by modeling a swirl-stabilized combustor in radial injection. The aim is to give more information on the influence of hydrogen enrichment on the CH4 flame properties and a comparison is presented between the experimental results and the numerical solutions by considering a flame shape, temperature field, and velocity vectors. II. NUMERICAL METHODOLOGY The conservation equations governing reactive flows solved by code FLUENT [26] are mass, momentum, species, and energy using several models described with finite volume methods [27]. The schematic of the partially premixed swirl- stabilized burner configuration can be seen in [5]. The grid is shown in Figure 1. The geometry is open to atmospheric pressure and the meshing of all computational domain examined in this work was set up with GAMBIT [24]. The meshing is non-uniform, with a fine meshing in zones of great interest. The numerical simulations carried out in this study were generated using a segregated, implicit solver. The coupled equations were solved with second-order accuracy in time and momentum, continuity, and turbulence parameters. The algorithms PRESTO and SIMPLE were respectively used for the pressure interpolation and for the coupling of pressure and velocity. A uniform velocity profile was assumed for the inlet conditions, and the inlet turbulence intensity was set to 11%. At the inlet, the boundary conditions were validated by the experiments. Theoretically, zero gradient condition is taken at the outlet boundary. These boundary conditions can be listed as follows: inlet, outlet, and wall boundary conditions. The DO model was used to calculate the radiating heat transfer of the unconfined flames [28]. The experimental apparatus is a laboratory scale swirl burner [5]. It consists of two concentric pipes with the annulus (outer radius, Rb = 18mm) supplying the central pipe (rint = 7.5mm) injecting the fuel mixtures (swirled air). The fuel mixture is radially injected into the swirled air stream through 8 circular holes. The holes are distributed on the circumferential section located 3mm upstream from the exit of the burner. Fig. 1. Numerical grid of the burner/furnace system. III. TURBULENCE MODEL To solve the turbulence problem, the two equation-based RANS standard k-ɛ turbulence model was used [29]. It is a semi-empirical turbulence model based on isotropic eddy- viscosity hypotheses. The RNG k-ɛ model is considered as a new version of the k-ɛ family of turbulence models. It is derived by a statistical technique called renormalization group method [30]. IV. RESULTS AND DISCUSSION A. Flame Shape The consequence of hydrogen addition is studied using CFD and compared against measurement data [31, 32]. Streamlines were obtained by two-dimensional (2-D) PIV measurements, in a region located immediately downstream of the exit throat of the burner. Flame imaging of the spontaneous flame emission is used to gain a general understanding (see Figure 2(a)). We can see the simulation plot from the right side. Although these diagnostics provide a somewhat different view of the flame, they provide complementary information. In general, the simulated image of the investigated flame (100% CH4) shows a very similar shape compared to the experimental ones, typical of high swirl processes in transverse injection. The comparison shows that standard k-ɛ is slightly over- predicted at most locations whereas the RNG k-ɛ model results are quite good [33]. Figure 2(a) shows the sudden expansion nozzle at the lower part of the combustion chamber and two locations of high product formation rates are located: one is on centerline in an upside-down V-shape, which can also be clearly seen in the photo of the flame in the experiment and the second in red color indicates a high product formation rate, which means an area with high chemical reaction production rates. The use of finite-rate-Eddy dissipation model generates a reaction at the mixing zone. Engineering, Technology & Applied Science Research Vol. 11, No. 2, 2021, 6943-6948 6945 www.etasr.com Bouziane et al.: Investigation of Swirl Stabilized CH4 Air Flame with Varied Hydrogen Content by using … (a) 100% CH4 40% H2+60% CH4 80% H2 + 20% CH4 100% H2 (b) 100% CH4 40% H2+60% CH4 80% H2 + 20% CH4 100% H2 Fig. 2. Flame shape and position comparison of the simulation (right) and the flame photo from the experiment (left): (a) standard k-ɛ model, (b) RNG k-ɛ model. The experimental reaction takes place at the center of the zone where a strong inner recirculation zone appears, whereas reaction occurs with a certain delay, when simulation is considered. The impact of the different percentage of H2 blending to the CH4–air flame on the global flame characteristics is shown in Figure 2 (different flame shapes and positions). (a) (b) (c) (d) (e) (f) Fig. 3. Radial distribution of mean temperature: (a) RNG k-ɛ 50% CH4 - 50% H2, (b) standard k-ɛ 50% CH4 - 50% H2, (c) RNG k-ɛ 20% CH4 - 80% H2, (d) standard k-ɛ 20% CH4 - 80% H2, (e) RNG k-ɛ 100% CH4, and (f) standard k-ɛ 100% CH4. Engineering, Technology & Applied Science Research Vol. 11, No. 2, 2021, 6943-6948 6946 www.etasr.com Bouziane et al.: Investigation of Swirl Stabilized CH4 Air Flame with Varied Hydrogen Content by using … It can be observed both from the experiments and the numerical simulation that the structure of the flame, i.e. central recirculation zone shape, changes with H2 addition. Enhances in H2 concentration in the mixture increase the axial velocity and decreases the size and strength of the recirculation region. Similar behavior was reported in [34]. B. Temperature Field To provide more information in this section, the computational solutions were compared with the measurements. To examine the role of turbulence closure models, the same configuration was modeled using the standard k-ɛ and RNG k- ɛ turbulence models. The mean temperature profiles are shown in Figure 3, reporting a different trend between NG and H2/NG fuel mixtures. Isothermal 100% CH4 Fig. 4. Experiment (left) and simulation (right) velocity vectors. (a) (b) Fig. 5. Velocity fluctuations. Comparison of the experimental data and simulation results: (a)100 % CH4, (b) 20% CH4 + 80% H2 flame. The calculated temperature profiles are in agreement with the experimental results. When moving on radial direction, the numerical simulation agrees well with the trend of the experimental data. The calculated temperature indicating the existence of the hot area of the flame and then finally to the limit which assigns to the recirculation zone. The temperature profiles can be split in into three separate areas as reported in Figures 5 and 6. The same results were obtained in [35]. The highest temperature shows that the flame takes place in the inner recirculation zone. When burning 100% NG, the peak temperature of the flame is quite close to 1150K, while its experimental value is estimated as 1250K. The peak temperature of the 80% H2 flame increases to about 1800K, as we seen in the experiment [31, 32]. An unexpected temperature value is determined by the CFD when the k-e model is considered. At the periphery, a region characterized by quite uniform temperature values is established. This is related with the formation of a corner recirculating flow induced by the radial expansion of the air stream and by the wall confinement. This model succeeds to predict the hot zone of the flame while over predict the calculated temperature. When r = 0.05m, the model agrees with the experimental results. On the other hand, a good agreement is obtained when the RNG k- ɛ model is used. The RNG k-ɛ model results are in good accordance with the experimental measurements, and are comparable to the results obtained with the k-ɛ model, despite the lack of accounting for non-isotropic turbulent shear stresses expected in the swirling flow. It can be observed that the predicted radial profiles of temperature of the two models are reasonable at most of the axial locations. C. Mean Flow Fields To provide more information on the nature of the flow, we focus our solutions on the behavior of this type of flame by showing the existence of the Central Recirculation Zone (CRZ). The numerical solution shows an open swirl flow core. The Figure points out a more closed form with a small but strong recirculation region. The reaction zone estimated by simulation is not located on the axis as indicated by the measurements. The dynamic flow field is affected greatly by the combustion that redefines a velocity distribution. This discrepancy is mainly the result of shortcoming of the turbulence model which might overestimate the decay of turbulence intensity for these swirling flames. Under- combustion regime (Figure 4), in a region located immediately in the exit of the burner clearly indicates the existence of a CRZ. The width of the CRZ is smaller than under when burning 100% NG. Similar behavior was reported in [29]. As the H2 percentage increases, the above regions reduce in size and move closer to the reactant outlet. When burning 100% NG, the axial flow fields by introducing velocity vectors obtained by the simulation are more large than the 80% H2 flame. Another important parameter describing the flow field, is the velocity fluctuations (Figure 5). Both results show an increase of the intensity a few mm downstream of the respective flame zone. A CRZ can be found, where the axial velocity is smaller than zero. This recirculation increases the flow divergence in the combustion chamber. In general, the numerical solutions agree well with the experimental results Engineering, Technology & Applied Science Research Vol. 11, No. 2, 2021, 6943-6948 6947 www.etasr.com Bouziane et al.: Investigation of Swirl Stabilized CH4 Air Flame with Varied Hydrogen Content by using … (100% methane flame), but we failed to predict the inner recirculation zone with precision. V. CONCLUSION Numerical simulations of a turbulent (CH4-H2)/air flame were conducted in this study. The H2 percentage varied in the mixture from 0% to 100%. The numerical investigation involved simulations employing the standard k-ɛ and the RNG k-ɛ for turbulence closure. The Fr-ED model is used to represent the turbulence/chemistry interactions. Numerical solutions were compared to experimental measurements which were used for model validation. The results show that the structure of the flame is described reasonably in both standard k-ɛ and RNG k-ɛ models which can predict satisfactorily the flame shape with the RNG k-ɛ model providing better performance than the k-e model. From the flame simulations, the hydrogen addition alters considerably the global flame characteristics by different flame shapes and flame positions. It can be observed from the experiments and the simulations that the structure of the flame varies with hydrogen blending. In general, the temperature data were well represented. The RNG k-ɛ model consistently demonstrated superiority over the standard k-ɛ turbulence model. The dynamic flow field is affected greatly by the combustion. A region located immediately downstream of the exit throat of the burner clearly indicates the existence of a CRZ. Under combustion conditions the width of the CRZ is smaller than under when burning 100% NG. It was concluded that the use of a 2D axial symmetric grid is a promising toolbox for modeling turbulent combustion and can be used for predicting the operating conditions of complex industrial furnaces. NOMENCLATURE �� Mass fraction of a particular reactant � ��,� � Stoichiometric coefficient for reactant in reaction ��,�: �� Stoichiometric coefficient for product in reaction N Number of chemical species in the system � ,� Molecular weight of species � ,� Molecular weight of species� � ,� Molecular weight of a particular reactant � k/ε Large-eddy mixing time scale � Empirical constant equal to 4.0 � Empirical constant equal to 0.5. α Thermal diffusivity [m 2 . s- 1 ] β Thermal expansion coefficient [K -1 ] φ Solid volume fraction Ɵ Dimensionless temperature µ Dynamic viscosity [kg. m -1 .s -1 ] REFERENCES [1] R. W. Schefer, "Hydrogen enrichment for improved lean flame stability," International Journal of Hydrogen Energy, vol. 28, no. 10, pp. 1131–1141, Oct. 2003, https://doi.org/10.1016/S0360-3199(02)00199-4. [2] S. O. Bade Shrestha and G. A. Karim, "Hydrogen as an additive to methane for spark ignition engine applications," Hydrogen as an additive to methane for spark ignition engine applications, vol. 24, no. 6, pp. 577–586, 1999. [3] A. A. Konnov, G. P. Alvarez, I. V. Rybitskaya, and J. D. Ruyck, "The Effects of Enrichment by Carbon Monoxide on Adiabatic Burning Velocity and Nitric Oxide Formation in Methane Flames," Combustion Science and Technology, vol. 181, no. 1, pp. 117–135, Dec. 2008, https://doi.org/10.1080/00102200802380173. [4] A. R. Choudhuri and S. R. Gollahalli, "Combustion characteristics of hydrogen–hydrocarbon hybrid fuels," International Journal of Hydrogen Energy, vol. 25, no. 5, pp. 451–462, May 2000, https://doi.org/10.1016/ S0360-3199(99)00027-0. [5] F. Cozzi and A. Coghe, "Behavior of hydrogen-enriched non-premixed swirled natural gas flames," International Journal of Hydrogen Energy, vol. 31, no. 6, pp. 669–677, May 2006, https://doi.org/10.1016/ j.ijhydene.2005.05.013. [6] A. J. Akande, "Production of hydrogen by reforming of crude ethanol," Ph.D. dissertation, University of Saskatchewan, Saskatoon, Saskatchewan, Canada, 2005. [7] S. S. Shy, Y. C. Chen, C. H. Yang, C. C. Liu, and C. M. Huang, "Effects of H2 or CO2 addition, equivalence ratio, and turbulent straining on turbulent burning velocities for lean premixed methane combustion," Combustion and Flame, vol. 153, no. 4, pp. 510–524, Jun. 2008, https://doi.org/10.1016/j.combustflame.2008.03.014. [8] A. R. Choudhuri and S. R. Gollahalli, "Effects of Isothermal Mixing Process on the Stability of Hydrogen-hydrocarbon Blend Turbulent Jet Diffusion Flames," in International Joint Power Generation Conference, Florida, FL, USA, Jul. 2000. [9] M. Karbasi and I. Wierzba, "The effects of hydrogen addition on the stability limits of methane jet diffusion flames," International Journal of Hydrogen Energy, vol. 23, no. 2, pp. 123–129, Feb. 1998, https://doi.org/10.1016/S0360-3199(97)00031-1. [10] P. Kumar and D. P. Mishra, "Experimental investigation of laminar LPG–H2 jet diffusion flame," International Journal of Hydrogen Energy, vol. 33, no. 1, pp. 225–231, Jan. 2008, https://doi.org/ 10.1016/j.ijhydene.2007.09.023. [11] A. R. Choudhuri and S. R. Gollahalli, "Characteristics of hydrogen– hydrocarbon composite fuel turbulent jet flames," International Journal of Hydrogen Energy, vol. 28, no. 4, pp. 445–454, Apr. 2003, https://doi.org/10.1016/S0360-3199(02)00063-0. [12] H. J. Burbano, A. A. Amell, and J. M. Garcia, "Effects of hydrogen addition to methane on the flame structure and CO emissions in atmospheric burners," International Journal of Hydrogen Energy, vol. 33, no. 13, pp. 3410–3415, Jul. 2008, https://doi.org/10.1016/j.ijhydene. 2008.04.020. [13] C. J. Tseng, "Premixed combustion of hydrogen/methane mixtures in a porous medium burner," in 36th Intersociety Energy Conversion Engineering Conference, Savannah, GE, USA, Aug. 2001. [14] S. Naha and S. K. Aggarwal, "Fuel effects on NOx emissions in partially premixed flames," Combustion and Flame, vol. 139, no. 1, pp. 90–105, Oct. 2004, https://doi.org/10.1016/j.combustflame.2004.07.006. [15] G. J. Rortveit and J. E. Hustad, "NOx Formations in Diluted CH4/H2 Counterflow Diffusion Flames," International Journal of Energy for a Clean Environment, vol. 4, no. 4, 2003, https://doi.org/10.1615/ InterJEnerCleanEnv.v4.i4.20. [16] S. A. El-Ghafour, "Pollutant emissions control utilizing hybrid fuel and air nozzle-cascading techniques," M.S. thesis, Port Said University, Port- Said, Egypt, 2011. [17] J. B. Heywood and F. R. Vilchis, "Comparison of Flame Development in a Spark-Ignition Engine Fueled with Propane and Hydrogen," Combustion Science and Technology, vol. 38, no. 5–6, pp. 313–324, Jul. 1984, https://doi.org/10.1080/00102208408923778. [18] C. K. Law, Combustion Physics. Cambridge, UK: Cambridge University Press, 2006. [19] Y. Lafay, B. Renou, G. Cabot, and M. Boukhalfa, "Experimental and numerical investigation of the effect of H2 enrichment on laminar methane–air flame thickness," Combustion and Flame, vol. 153, no. 4, pp. 540–561, Jun. 2008, https://doi.org/10.1016/j.combustflame. 2007.10.002. [20] M. Fukuda, K. Korematsu, and M. Sakamoto, "On Quenching Distance of Mixtures of Methane and Hydrogen with Air," Bulletin of JSME, vol. 24, no. 193, pp. 1192–1197, 1981, https://doi.org/10.1299/jsme1958. 24.1192. [21] Z. Dulger, "Adiabatic flame temperature and product composition for lean combustion of hydrogen-methane combination," The American Society of Mechanical Engineers, vol. 2, pp. 707–712, 1994. Engineering, Technology & Applied Science Research Vol. 11, No. 2, 2021, 6943-6948 6948 www.etasr.com Bouziane et al.: Investigation of Swirl Stabilized CH4 Air Flame with Varied Hydrogen Content by using … [22] A. C. Mangra, "Design and Numerical Analysis of a Micro Gas Turbine Combustion Chamber," Engineering, Technology & Applied Science Research, vol. 10, no. 6, pp. 6422–6426, Dec. 2020, https://doi.org/ 10.48084/etasr.3835. [23] M. W. Khalid and M. Ahsan, "Computational Fluid Dynamics Analysis of Compressible Flow Through a Converging-Diverging Nozzle using the k-ε Turbulence Model," Engineering, Technology & Applied Science Research, vol. 10, no. 1, pp. 5180–5185, Feb. 2020, https://doi.org/ 10.48084/etasr.3140. [24] J. B. Bell, R. K. Cheng, M. S. Day, and I. G. Shepherd, "Numerical simulation of Lewis number effects on lean premixed turbulent flames," Proceedings of the Combustion Institute, vol. 31, no. 1, pp. 1309–1317, Jan. 2007, https://doi.org/10.1016/j.proci.2006.07.216. [25] M. Day, J. Bell, P.-T. Bremer, V. Pascucci, V. Beckner, and M. Lijewski, "Turbulence effects on cellular burning structures in lean premixed hydrogen flames," Combustion and Flame, vol. 156, no. 5, pp. 1035–1045, May 2009, https://doi.org/10.1016/j.combustflame. 2008.10.029. [26] G. Solero, A. Coghe, and G. Scribano, "Influence of natural gas injection procedure in a swirl Burner," in 7th Mediterranean Combustion Symposium, Sardinia, Italy, Sep. 2011. [27] M. R. Lecic, A. S. Cocic, and S. M. Cantrak, "Original measuring and calibration equipment for investigation of turbulent swirling flow in circular pipe," Experimental Techniques, vol. 38, no. 3, pp. 54–62, May 2014, https://doi.org/10.1111/j.1747-1567.2012.00812.x. [28] B. E. Launder, Lectures in Mathematical Models of Turbulence. London, New York: Academic Press, 1972. [29] V. Yakhot and S. A. Orszag, "Renormalization group analysis of turbulence. I. Basic theory," Journal of Scientific Computing, vol. 1, no. 1, pp. 3–51, Mar. 1986, https://doi.org/10.1007/BF01061452. [30] F. Cozzi, A. Olivani, and A. Coghe, "Combustion of NG + H2 fuel mixtures in a partially premixed swirlburner," presented at the Third European Combustion Meeting ECM 2007, 2007. [31] A. Olivani, F. Cozzi, G. Solero and A. Coghe, "Analysis of the Environmental Impact of a Swirl Burner:Hydrogen and Natural Gas Mixture Feeding," presented at the 28th Meeting of the Italian Section of The Combustion Institute, 2005 [32] J. F. Widmann, S. R. Charagundla, and C. Presser, Benchmark Experimental Database for Multiphase Combustion Model Input and Validation: Characterization of the Inlet Combustion Air. Gaithersburg, MD, USA: National Institute of Standards and Technology, 1999. [33] M. W. Eldeen, M. A. Taemah, and W. M. El-Maghalany, "Influence of Swirl Intensity and Blended Methane-Hydrogen on Gas Turbine Combustion Characteristics," International Journal of Engineering Research & Technology, vol. 7, no. 7, pp. 405–412, 2018. [34] A. Frassoldati, S. Frigerio, E. Colombo, F. Inzoli, and T. Faravelli, "Determination of NOx emissions from strong swirling confined flames with an integrated CFD-based procedure," Chemical Engineering Science, vol. 60, no. 11, pp. 2851–2869, Jun. 2005, https://doi.org/10.1016/j.ces.2004.12.038.