基于OpenFOAM的三维H2/Air连续旋转爆轰流场数值模拟

doi:10.13675/j.cnki.tjjs.190370

推进技术 ›› 2020, Vol. 41 ›› Issue (6): 1351-1360.DOI: 10.13675/j.cnki.tjjs.190370

基于OpenFOAM的三维H₂/Air连续旋转爆轰流场数值模拟

孟豪龙¹，翁春生¹，武郁文¹，郑权¹，魏万里¹

南京理工大学瞬态物理国家重点实验室，江苏南京 210094

发布日期:2021-08-15
作者简介:孟豪龙，博士生，研究领域为爆轰推进。E-mail：18362916356@163.com
基金资助:
国家自然科学基金（11802137；11702143）；中央高校基本科研业务专项资金（30919011259；30918011343）。

Three-Dimensional Numerical Simulation of H2/Air Continuous Rotating Detonation Flow Field Based on OpenFOAM

National Key Lab of Transient Physics,Nanjing University of Science and Technology,Nanjing 210094,China

Published:2021-08-15

摘要/Abstract

摘要： 为进一步研究旋转爆轰流场特征，基于开源计算流体动力学软件OpenFOAM，采用9组分19步的基元化学反应模型，对H₂/Air连续旋转爆轰流场进行了三维数值模拟，得到了旋转爆轰波稳定传播时燃烧室内部流场的详细结构，研究了燃烧室头部激波的传播特性，分析了旋转爆轰燃烧室的压力增益性能。结果表明：旋转爆轰波后的第一道反射激波在由燃烧室外壁面向内壁面传播过程中反射激波的高度增加并在靠近内壁面附近与滑移线交汇形成局部高温高压区域；旋转爆轰波在外壁面位置处相位约落后于内壁面0.003rad~0.15rad，其相位差随燃烧室曲率差的增加而增大；燃烧室头部反射激波数目受到曲率差和进气总压的影响，燃烧室曲率差增大，反射激波数目减少，进气总压增大，反射激波数目增多；燃烧室压力增益保持在0.3以上，在进气总压一定的条件下，压力增益随着燃烧室曲率差的增大有增加的趋势。研究结果揭示了三维旋转爆轰流场的精细结构和燃烧室头部激波的传播规律。

关键词: 旋转爆轰流场;反射激波;传播特性;数值模拟;压力增益

Abstract: To further study the flow field characteristic of rotating detonation, based on the open source computational fluid dynamics (CFD) software OpenFOAM, 19 species and 9-step reaction kinetics model was used to simulate the three-dimensional flow field of H₂/Air continuous rotating detonation. The detailed structure of the flow field in the combustor was obtained when the rotating detonation wave propagated steadily. The propagation characteristics of shock wave near the combustor head were studied. The pressure gain performance of rotating detonation combustor was analyzed. The results show that the first reflected shock wave after rotating detonation wave propagates from the outer wall of combustion chamber to the inner wall, and the height of reflected shock wave increases, and a local high temperature and pressure region is formed near the inner wall and intersects with the slip line. The phase of the rotating detonation wave on the outer wall is about 0.003rad to 0.15rad behind that on the inner wall, and the phase difference increases with the increase of the difference of combustion chamber curvature. The number of reflected shocks near the head of combustion chamber is affected by the stagnation pressure and the difference of curvature. The number of reflected shocks decreases with the increase of the difference of combustion chamber curvature, while the number of reflected shocks increases with the increase of stagnation pressure. The pressure gain of combustor keeps above 0.3. When the stagnation pressure is constant, the pressure gain increases with the increase of the curvature difference of the combustor. The results reveal the fine structure of the three-dimensional rotating detonation flow field and the propagation law of shock wave near the combustor head.

Key words: Rotating detonation flow field;Reflected shock wave;Propagation characteristics;Numerical simulation;Pressure gain

孟豪龙，翁春生，武郁文，郑权，魏万里. 基于OpenFOAM的三维H₂/Air连续旋转爆轰流场数值模拟[J]. 推进技术, 2020, 41(6): 1351-1360.

MENG Hao-long¹, WENG Chun-sheng¹, WU Yu-wen¹, ZHENG Quan¹, WEI Wan-li¹. Three-Dimensional Numerical Simulation of H2/Air Continuous Rotating Detonation Flow Field Based on OpenFOAM[J]. Journal of Propulsion Technology, 2020, 41(6): 1351-1360.

参考文献

[1] Kailasanath K. Review of Propulsion Applications of Detonation Waves[J]. AIAA Journal， 2000， 38(9): 1698-1708.
[2] 王健平，周蕊，武丹. 连续旋转爆轰发动机的研究进展[J]. 实验流体力学， 2015， 29(4): 12-25.
[3] Bykovskii F A， Zhdan S A， Vedernikov E F， et al. Effect of Combustor Geometry on Continuous Spin Detonation in Syngas-Air Mixtures[J]. Combustion Explosion and Shock Waves， 2015， 51(6): 688-699.
[4] Bykovskii F A， Zhdan S A. Current Status of Research of Continuous Detonation in Fuel-Air Mixtures[J]. Combustion Explosion and Shock Waves， 2015， 51(1): 21-35.
[5] Kindracki J， Wolanski P， Gut Z. Experimental Research on the Rotating Detonation in Gaseous Detonation in Gaseous Fuels-Oxygen Mixtures[J]. Shock Wave， 2011， 21（2）: 75-84.
[6] Wolański P. Detonative Propulsion[J]. Proceedings of the Combustion Institute， 2013， 34(1): 125-158.
[7] 刘世杰，刘卫东，林志勇，等. 连续旋转爆震波传播过程研究（I）：同向传播模式[J]. 推进技术， 2014， 35(1): 138-144.
[8] 刘世杰，林志勇，刘卫东，等. 连续旋转爆震波传播过程研究（II）：双波对撞传播模式[J]. 推进技术， 2014， 35(2): 269-275.
[9] Lin W， Zhou J， Liu S， et al. An Experimental Study on CH₄/O₂ Continuously Rotating Detonation Wave in a Hollow Combustion Chamber[J]. Experimental Thermal and Fluid Science， 2015， 62: 122-130.
[10] Wang Y H， Wang J P. Coexistence of Detonation with Deflagration in Rotating Detonation Engines[J]. International Journal of Hydrogen Energy， 2016， 41(32):14302-14309.
[11] Wang Y H， Wang J P. Effect of Equivalence Ratio on the Velocity of Rotating Detonation[J]. International Journal of Hydrogen Energy， 2015， 40（25）: 7949-7955.
[12] 郑权，翁春生，白桥栋. 当量比对液体燃料旋转爆轰发动机爆轰影响实验研究[J]. 推进技术， 2015， 36(6): 947-952.
[13] Zhdan S A， Bykovskii F A， Vedernikov E F. Mathematical Modeling of a Rotating Detonation Wave in a Hydrogen-Oxygen Mixture[J]. Combustion Explosion and Shock Waves， 2007， 43(4): 449-459.
[14] Davidenko D M， Gokalp I， Kudryavtsev A N. Numerical Simulation of the Continuous Rotating Hydrogen-Oxygen Detonation with a Detailed Chemical Mechanism[C]. Moscow: Proceedings of the West-East High Speed Flow Field Conference， 2007.
[15] 潘振华，范宝春，张旭东，等. 连续旋转爆轰三维流场的数值模拟[J]. 兵工学报， 2012， 33(5): 594-599.
[16] Hishida M， Fujiwara T， Wolanski P. Fundamentals of Rotating Detonations[J]. Shock Waves， 2009， 19(1): 1-10.
[17] Schwer D， Kailasanath K. Numerical Study of the Effects of Engine Size on Rotating Detonation Engines[C]. Florida： 49th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition， 2011.
[18] Eude Y， Davidenko D， Falempin F， et al. Use of the Adaptive Mesh Refinement for 3D Simulations of a CDWRE (Continuous Detonation Wave Rocket Engine) [C]. California：AIAA International Space Planes & Hypersonic Systems & Technologies Conference， 2011.
[19] Zhou R， Wang J P. Numerical Investigation of Shock Wave Reflections near the Head Ends of Rotating Detonation Engines[J]. Shock Waves， 2013， 23(5): 461-472.
[20] Braun J， Saracoglu B H， Paniagua G. Unsteady Performance of Rotating Detonation Engines with Different Exhaust Nozzles[J]. Journal of Propulsion and Power， 2016， 33(1): 1-10.
[21] Braun J， Saracoglu B H， Magin T E. et al. One-Dimensional Analysis of the Magnetohydro-Dynamic Effect in Rotating Detonation Combustors[J]. AIAA Journal， 2016, 54(12): 1-7.
[22] Sousa J， Braun J， Paniagua G. Development of a Fast Evaluation Tool for Rotating Detonation Combustors[J]. Applied Mathematical Modelling， 2017， 52: 42-52.
[23] Weller H G， Tabor G， Jasak H C， et al. A Tensorial Approach to Computational Continuum Mechanics Using Object Orientated Techniques[J]. Computers in Physics， 1998， 12(6): 620-631.
[24] Florian E， Vollmer K G， Thomas S. Numerical Simulation of the Deflagration-to-Detonation Transition in Inhomogeneous Mixtures[J]. Journal of Combustion， 2014，（5）: 1-15.
[25] Yamada T， Hayashi A K， Yamada E， et al. Detonation Limit Thresholds in H₂/O₂ Rotating Detonation Engine[J]. Combustion Science and Technology， 2010， 182，(11-12).

基于OpenFOAM的三维H₂/Air连续旋转爆轰流场数值模拟

Three-Dimensional Numerical Simulation of H2/Air Continuous Rotating Detonation Flow Field Based on OpenFOAM

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 0

编辑推荐

Metrics

本文评价