Second-order modeling of non-premixed turbulent methane-air combustion

Abstract

The main purpose of this research is the second-order modeling of flow and turbulent heat flux in non-premixed methane-air combustion. A turbulent stream of non-premixed combustion in a stoichiometric condition, is numerically analyzed through the Reynolds averaged Navier-Stokes (RANS) equations. For modeling radiation and combustion, the discrete ordinates (DO) and eddy dissipation concept model have been applied. The Reynolds stress transport model (RSM) also was used for turbulence modeling. For THF in the energy equation, the GGDH model and high order algebraic model of HOGGDH with simple eddy diffusivity model have been applied. Comparing the numerical results of the SED model (with the turbulent Prandtl 0.85) and the second-order heat flux models with available experimental data follows that applying the second-order models significantly led to the modification of predicting temperature distribution and species mass fraction distribution in the combustion chamber. Calculation of turbulent Prandtl number in the combustion chamber shows that the assumption of Pr_t of 0.85 is far from reality and Pr_t in different areas varies from 0.4 to 1.2.

Keywords

combustion modeling / turbulent Prandtl number / second-order models / Reynolds stress transport model

References

[1] TurnsS R, MantelS JAn introduction to combustion [M], 20002nd eds.New York, McGraw Hill

[2] WisemanS, RirthM, GruberA, DawsonJ, ChenJ. A comparison of the blow-out behavior of turbulent premixed ammonia/ hydrogen/nitrogen-air and methane — air flames [J]. Proceedings of the Combustion Institute, 2021, 38(2): 2869-2876

[3] LeeJ, MuellerM. Closure modeling for the conditional Reynolds stresses in turbulent premixed combustion [J]. Proceedings of the Combustion Institute, 2021, 38(2): 3031-3038

[4] YuS, LiuX, BaiX S, ElbazA M, RobertsW L. LES/PDF modeling of swirl-stabilized non-premixed methane/air flames with local extinction and re-ignition [J]. Combustion and Flame, 2020, 219: 102-119

[5] ShinJ, GeY, LampmannA, PfitzerM. Adata-driven subgrid scale model in large eddy simulation of turbulent premixed combustion [J]. Combustion and Flame, 2021, 231: 111486

[6] LuoK, YangJ, BaiY, FanJ. Large eddy simulation of turbulent combustion by a dynamic Second-order moment closure model [J]. Fuel, 2017, 187457-467

[7] ZianiL, ChakerA, ChetehounaK, MalekA, MahmahB. Numerical simulations of non-premixed turbulent combustion of CH₄ — H₂ mixtures using the PDF approach [J]. Int J Hydrogen Energy, 2013, 38(20): 8597-8603

[8] RoyR N, SreedharaS. Modelling of methanol and H2/CO bluffbody flames using RANS based turbulence models with conditional moment closure model [J]. Appl Therm Eng, 2016, 93: 61-70

[9] WEN X, WANG P, MAN S. Numerical study of performance of reverse flow combustor [C]// IEEE Aerospace Conference, 2016: 1–13. DOI: https://doi.org/10.1109/AERO.2016.7500596.

[10] NORWAZAN A, JAAFAR M. Studies of isothermal swirling flows with different RANS models in unconfined burner [C]// IEEE, World Congress on Sustainable Technologies (WCST-2014), 2014: 48–53. DOI: https://doi.org/10.1109/WCST.2014.7030095.

[11] AlemiE, Rajabi-ZargarabadiM. Effects of jet characteristics on NO formation in a jet-stabilized combustor [J]. Int Journal of Thermal Science, 2017, 112: 55-67

[12] YilmazH, CamO, TangozS, YilmazI. Effect of different turbulence models on combustion and emission characteristics of hydrogen/air flames [J]. Int Journal of Hydrogen Energy, 2017, 42(40): 25744-25755

[13] HossainM, JonesJC, MalalasekeraW. Modelling of a bluff-body non-premixed flame using a coupled radiation/flamelet combustion model [J]. Flow Turbulence and Combustion, 2001, 67: 217-234

[14] SilvaC V, FrancaF H, VielmoH A. Analysis of the turbulent, non-premixed combustion of natural gas in a cylindrical chamber with and without thermal radiation [J]. Combust Sci and Tech, 2007, 179: 1605-1630

[15] Garre’tonD, SimoninO. Aerodynamics of steady state combustion chambers and furnaces. Final results [C]. Proceedings of the First Workshop, 1994, Chatou, France, EDF, 2935

[16] LaunderB E. On the computation of convective heat transfer in complex turbulent flows [J]. Journal of Heat Transfer, 1986, 110: 1112-1128

[17] WikstromP M. Derivation and investigation of a new explicit algebraic model for the passive scalar flux [J]. Physics of Fluids, 2000, 12: 688-702

[18] AbeK, SugaK. Towards the development of a Reynolds-averaged algebraic turbulent scalar flux model [J]. Int Journal of Heat and Fluid Flow, 2001, 22: 19-29

[19] RogersM, MansourN, ReynoldsW C. An algebraic model for the turbulent flux of a passive scalar [J]. Journal of Fluid Mechanics, 1989, 203: 77-101

[20] NaganoY, ShimadaM. Development of a two-equation heat transfer model based on direct simulations of turbulent flows with different Prandtl numbers [J]. Physics of Fluids, 1996, 8: 3379-3402

[21] YangW-h, JiangS-j, HsiaoT-c, YangL-xing. Numerical simulation of high temperature air combustion flames properties [J]. Journal of Central South University, 2000, 7(3): 156-158

[22] ChenH, LiuW-qiang. Numerical study of effect of front cavity on hydrogen/air premixed combustion in a microcombustion chamber [J]. Journal of Central South University, 2019, 26(8): 2259-2271

[23] DalyB J, HarlowF H. Transport Equation in Turbulence [J]. Phys. Fluids, 1970, 13: 2634-2649

[24] SugaK, AbeK. Nonlinear eddy viscosity modeling for turbulence and heat transfer near wall and shear-free boundaries [J]. Int Journal of Heat and Fluid Flow, 2000, 21: 37-48

[25] MagelH, SchnellU, HeinK R G. Simulation of detailed chemistry in a turbulent combustor flow [J]. Symposium (International) on Combustion, 1669, 26(1): 67-74

[26] CentenoF R, Da-SilvaC V, BrittesR, FrancaFHR. Numerical simulations of the radiative transfer in a 2D axisymmetric turbulent non-premixed methane–air flame using up-to-date WSGG and gray-gas models [J]. J Braz Soc Mech Sci Eng, 2015, 37(6): 1839-1850

[27] Da-SilvaC V, DeonD L, CentenoF R, FrancaF H R, PereiraFM. Assessment of combustion models for numerical simulations of a turbulent non-premixed natural gas flame inside a cylindrical chamber [J]. Combustion Science and Technology, 2018, 190(9): 1528-1556

[28] CentenoF R, BrittesR, FrancaF H R, Da-SilvaC V. Application of the WSGG model for the calculation of gas–soot radiation in a turbulent non-premixed methane-air flame inside a cylindrical combustion chamber [J]. Int J of Heat and Mass Transfer, 2016, 93: 742-753

[29] ChouaiebS, KriaaW, MhiriH, BournotP. Presumed PDF modeling of micro jet assisted CH₄-H₂/air turbulent flames [J]. Energy Conversion and Management, 2016, 120: 412-421

[30] LuoX, LiuP-l, LuoH-an. Improvement of Prandtl mixing length theory and application in modeling of turbulent flow in circular tubes [J]. Journal of Central South University, 2008, 15: 774-778

[31] JiangG-b, HeY-s, ShuS, XiaoY-xiong. Numerical prediction of inner turbulent flow in conical diffuser by using a new five-point scheme and DLR k — e turbulence model [J]. Journal of Central South University, 2008, 15(S1): 181-186

[32] HinzeJ OTurbulence [M], 1975, New York, McGraw-Hill Publishing Co.

[33] ANSYS. FLUENT user’s manual, 2016, Version 17 [M]. ANSYS Inc, 2016.

[34] LiX, RenJ, JiangH. Application of algebraic anisotropic turbulence models to film cooling flows [J]. Int J of Heat and Mass Transfer, 2015, 91: 7-17

[35] SugaK, NagaokaM, HorinouchiN. Application of a higher order GGDH heat flux model to three-dimensional turbulent U-bend duct heat transfer [J]. Journal of Heat Transfer, Technical Notes, 2003, 125: 200-203

[36] Jia-qiangE, WuJ-h, LiuT, ChenJ-w, DengY-w, PengQ-g. Effects analysis on catalytic combustion characteristic of hydrogen/air in the micro turbine engine by fuzzy grey relation method [J]. Journal of Central South University, 2019, 26(8): 2214-2223

[37] HishidaM, NaganoY. Structure of turbulent velocity and temperature fluctuations in fully developed pipe flow [J]. J Heat Transfer, 1979, 101: 15-22

[38] Rajabi-ZargarabadiM, Bazdidi-TehraniF. Implicit algebraic model for predicting THF in film cooling flow [J]. Int J Numer Meth Fluids, 2010, 64: 517-531

[39] Bazdidi-TehraniF, Rajabi-ZargarabadiM. Application of second moment closure and higher order generalized gradient diffusion hypothesis to impingement heat transfer [J]. Transactions of the Canadian Society for Mechanical Engineering, 2008, 32(1): 91-106

[40] BaeY Y. A new formulation of variable turbulent Prandtl number for heat transfer to supercritical fluids [J]. Int J of Heat and Mass Transfer, 2016, 92: 792-806