Numerical study on turbulent combustion in high Mach number scram-jet engine considering thermal non-equilibrium effect
-
摘要: 高超声速飞行过程中,随着马赫数升高,边界层内温度快速上升,分子振动能、电子能被激发,气体分子发生离解甚至电离,使得量热完全气体假设失效。针对高温热化学非平衡气体,采用热力学平衡模型及Park的双温度非平衡模型对飞行马赫数12条件下的超燃冲压发动机进行了数值模拟。研究结果表明:与平衡态计算结果相比,非平衡效应使得波系位置前移,激波间峰值压力升高,对冷场的影响更加明显;非平衡态平动温度场与平衡态差别不大,热力学非平衡效应使平动温度略微升高;非平衡态燃烧室出口截面燃烧效率更低,热力学非平衡效应会降低反应程度。Abstract: During hypersonic flight, the temperature inside the boundary layer rapidly increases with the increase of Mach number, leading to the excitation of molecular vibration and electron energy. Gaseous molecules undergo dissociation or even ionization, making the assumption of complete gas calorimetry invalid, and thereby affecting the characteristics of the scram-jet engine. The turbulent combustion in the scram-jet engine under the Ma 12 flight condition is numerically studied via a thermodynamic equilibrium model and the Park's dual temperature non-equilibrium model where the high-temperature thermal non-equilibrium effects are considered. The results indicate that compared to the equilibrium case, the non-equilibrium effect causes the position of the shock wave trains to shift forward, and the increase of the peak pressure between shock waves. This is more significant for the frozen flow field. The temperature field Ttr in the non-equilibrium case is not significantly different from that of the equilibrium case, and the thermodynamic non-equilibrium effect slightly increases Ttr. The combustion efficiency at the outlet section is lower in the non-equilibrium case, and the thermodynamic non-equilibrium effect slightly weakens the intensity of the reactions.
-
表 1 组分特征振动温度
Table 1. The characteristic vibrational temperature of the species
组分 θv,s/K N2 3371 O2 2256 H2 6215 H2O 2294、5180、5400 OH 5375 HO2 1577、2059 H2O2 1250、1970、2030、2070、4130、4910 表 2 空气和燃料边界条件
Table 2. Boundary conditions for air and fuel stream
马赫数Ma 流速u∞/(m·s−1) 静压p∞/Pa 平动温度Ttr/K 空气入口 10 3000 1450 227 燃料入口 1 1208 541778 249 表 3 初始条件
Table 3. Initial conditions
Ttr/K Tv/K p∞/Pa XN2 XO2 XH2 2000 900 70000 0.69 0.21 0.10 -
[1] 陈贤亮, 符松. 热化学非平衡高超声速平板边界层线性稳定性分析[J]. 空气动力学学报, 2020, 38(2): 316–325. doi: 10.7638/koqdIlxxl-2019.0156CHEN X L, FU S. Linear stability analysis of hypersonic boundary layer on a flat-plate with thermal-chemical non-equilibrium effects[J]. Acta Aerodynamica Sinica, 2020, 38(2): 316–325. doi: 10.7638/koqdIlxxl-2019.0156 [2] 袁军娅, 任翔, 蔡国飙, 等. 双锥/双楔流动中的高温气体效应仿真模拟[J]. 气体物理, 2022, 7(4): 10–18. doi: 10.19527/j.cnki.2096-1642.0927YUAN J Y, REN X, CAI G B, et al. Simulation of high temperature gas effects in high enthalpy double cone/wedge flows[J]. Physics of Gases, 2022, 7(4): 10–18. doi: 10.19527/j.cnki.2096-1642.0927 [3] 傅杨奥骁, 董维中, 丁明松, 等. 高焓电弧风洞试验热化学非平衡流场数值模拟[J]. 实验流体力学, 2019, 33(3): 1–12. doi: 10.11729/syltlx20180138FU-Y A X, DONG W Z, DING M S, et al. Numerical simulation of thermochemical non-equilibrium flow field in arc-jet tunnel[J]. Journal of Experiments in Fluid Mechanics, 2019, 33(3): 1–12. doi: 10.11729/syltlx20180138 [4] 岳连捷, 张旭, 张启帆, 等. 高马赫数超燃冲压发动机技术研究进展[J]. 力学学报, 2022, 54(2): 263–288. doi: 10.6052/0459-1879-21-547YUE L J, ZHANG X, ZHANG Q F, et al. Research progress on high-mach-number scramjet engine technologies[J]. Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(2): 263–288. doi: 10.6052/0459-1879-21-547 [5] LONGO J M A, HANNEMANN K, HANNEMANN V. The challenge of modeling high speed flows[C]//Proceedings of the 6th EUROSIM Congress on Modelling and Simulation. 2007. [6] ANNA A. Numerical modeling of surface chemistry processes for hypersonic entry environments[D]. Ann Arbor, MI, US: The University of Michigan, 2013. [7] 熊文韬. 高温非平衡效应下双楔绕流中激波干扰研究[D]. 合肥: 中国科学技术大学, 2017.XIONG W T. On shock-shock interaction in double-wedge flow with high temperature non-equilibrium effects[D]. Hefei: University of Science and Technology of China, 2017. [8] BRICALLI M G, BROWN L, BOYCE R R, et al. Scramjet performance with nonuniform flow and swept nozzles[J]. AIAA Journal, 2018, 56(10): 3988–4003. doi: 10.2514/1.j056963 [9] CASSEAU V, ESPINOZA D E R, SCANLON T J, et al. A two-temperature open-source CFD model for hypersonic reacting flows, part two: multi-dimensional analysis[J]. Aerospace, 2016, 3(4): 45. doi: 10.3390/aerospace3040045 [10] CASSEAU V. An open-source CFD solver for planetary entry[D]. Glasgow, United Kingdom, University of Strathclyde, 2017. doi: 10.48730/fc4y-6m27. [11] LANDAU L D, TELLER E. On the theory of sound dispersion[J]. Physikalische Zeitschrift der Sowjetunion, 1936, 10(34): 1. doi: 10.1016/B978-0-08-010586-4.50027-4 [12] JOHNSON R D Ⅲ. NIST Computational chemistry comparison and benchmark database[R]. NIST Standard Reference Database Number 101, 2016. [13] MENTER F R, KUNTZ M, LANGTRY R. Ten years of industrial experience with the SST turbulence model[J]. Turbulence, heat and mass transfer, 2003, 4(1): 625–632. [14] CONAIRE M Ó, CURRAN H J, SIMMIE J M, et al. A comprehensive modeling study of hydrogen oxidation[J]. International Journal of Chemical Kinetics, 2004, 36(11): 603–622. doi: 10.1002/kin.20036 [15] AO Y, WU K, LU H B, et al. Combustion dynamics of high Mach number scramjet under different inflow thermal nonequilibrium conditions[J]. Acta Astronautica, 2023, 208: 281–295. doi: 10.1016/j.actaastro.2023.04.020 [16] BRICALLI M G, BROWN L M, BOYCE R R. Numerical investigation into the combustion behavior of an inlet-fueled thermal-compression-like scramjet[J]. AIAA Journal, 2015, 53(7): 1740–1760. doi: 10.2514/1.j053513 [17] 白贺之, 陈兵, 徐旭. 基于PNS方程的高超声速热化学非平衡喷管流动数值模拟[J]. 战术导弹技术, 2018(4): 25–31. doi: 10.16358/j.issn.1009-1300.2018.7.197BAI H Z, CHEN B, XU X. Numerical simulation of hypersonic nozzle flow in thermochemical nonequilibrium based on parabolized Navier–Stokes equation[J]. Tactical Missile Technology, 2018(4): 25–31. doi: 10.16358/j.issn.1009-1300.2018.7.197 [18] VATANSEVER D, ÇELİK B. An open-source hypersonic solver for non-equilibrium flows[J]. Journal of Aeronautics and Space Technologies (Havacilik ve Uzay Teknolojileri Dergisi), 2021, 14(1): 35–52.