基于旋转爆震三维流场结构分析的计算模型对比研究

doi:10.13675/j.cnki.tjjs.190536

推进技术 ›› 2020, Vol. 41 ›› Issue (12): 2757-2765.DOI: 10.13675/j.cnki.tjjs.190536

基于旋转爆震三维流场结构分析的计算模型对比研究

刘朋欣^1,2，郭启龙^1,2，赵炜¹，李辰^1,2，李沁³，张涵信^1,2

1.中国空气动力研究与发展中心计算空气动力研究所，四川绵阳 621000;2.中国空气动力研究与发展中心空气动力学国家重点实验室，四川绵阳 621000;3.厦门大学航空航天学院，福建厦门 361102

发布日期:2021-08-15
基金资助:
国家自然科学基金（11902345；91541105）；四川省科技计划（2019YJ0272）。

Computational Models Based on Analysis of Three Dimensional Flow Field Structures in Rotating Detonation

1.Computational Aerodynamics Institute，China Aerodynamics Research and Development Center，Mianyang 621000，China;2.State Key Laboratory of Aerodynamics，China Aerodynamics Research and Development Center，Mianyang 621000，China;3.School of Aerospace Engineering，Xiamen University，Xiamen 361102，China

Published:2021-08-15

摘要/Abstract

摘要： 为了研究不同计算模型对三维旋转爆震数值模拟的影响，分别基于Euler方程、N-S方程、RANS方法和IDDES方法并结合滑移、无滑移壁面边界条件，耦合氢气/空气的有限化学反应速率模型（7组分8基元反应），采用高分辨率的五阶有限差分格式WENO-PPM5离散对流项，对圆环筒型燃烧室内的三维旋转爆震波进行数值模拟，得到了不同模型下旋转爆震波的流场结构、传播特性和推力性能。对比了不同计算模型对流场结构的影响。当采用滑移壁面边界条件时，Euler方程和N-S方程的计算结果较为一致。当使用无滑移壁面边界条件时，边界层的存在会导致可燃混气与燃烧产物之间接触面上的爆燃燃烧区域沿壁面向上游渗透，增大爆燃区域的范围；且爆震波锋面非受限侧变窄、严重变形，不同计算方法计算的变形程度有所不同。

关键词: 燃烧室;旋转爆震;三维;计算模型;流场结构;传播特性;推力性能

Abstract: The effects of numerical model on three-dimensional rotating detonation in an annular chamber are investigated by using different models including Euler， Navier-Stokes（N-S） and Reynolds-Averaged N-S（RANS）， Improved Delayed Detached-Eddy- Simulation（IDDES） and different wall boundary conditions including slip wall and non-slip wall. The chemical sources are evaluated by 7-species-and-8-reaction finite-rate chemical reaction model of hydrogen/air mixtures， and the convection terms are discretized by a fifth order WENO-type scheme with high resolution， namely WENO-PPM5. The three-dimensional rotating detonation wave in an annular chamber is simulated. And the flow structure， propagation property and propulsion performance are obtained. Comparisons of flow structure are made among the results of different numerical models. When using slip wall condition， the results of Euler and N-S are very similar. While considering non-slip wall condition， the existence of boundary layer leads the deflagration at the interface between reactant and product extending to upstream along the wall， and enlarges the range of deflagration zone. Besides， the unbounded side of detonation front becomes narrow and appears serve deformation. The deformation degrees show some difference in different computational methods.

Key words: Combustor;Rotating detonation;Three-dimensional;Numerical model;Flow structures;Propagation characteristics;Thrust performance

刘朋欣，郭启龙，赵炜，李辰，李沁，张涵信. 基于旋转爆震三维流场结构分析的计算模型对比研究[J]. 推进技术, 2020, 41(12): 2757-2765.

LIU Peng-xin^1,2, GUO Qi-long^1,2, ZHAO Wei¹, LI Chen^1,2, LI Qin³, ZHANG Han-xin^1,2. Computational Models Based on Analysis of Three Dimensional Flow Field Structures in Rotating Detonation[J]. Journal of Propulsion Technology, 2020, 41(12): 2757-2765.

参考文献

[1] Wolanski Piotr. Detonation Propulsion［J］. Proceedings of the Combustion Institute， 2013， 34： 125-158.
[2] Frank K Lu， Braun Eric M. Rotating Detonation Wave Propulsion： Experimental Challenges， Modeling， and Engine Concepts［J］. Journal of Propulsion and Power， 2014， 30（5）： 1125-1142.
[3] Bykovskii F A， Zhdan S A. Current Status of Research of Continuous Detonation in Fuel-Air Mixtures （Review）［J］. Combustion， Explosion， and Shock Waves， 2015， 51（1）： 21-35.
[4] Kailasanath K. Recent Developments in the Research on Rotating Detonation Wave Engines［R］. AIAA 2017-0784.
[5] Zhou Rui， Wu Dan， Wang Jian-ping. Progress of Continuously Rotating Detonation Engines［J］. Chinese Journal of Aeronautics， 2016， 29（1）： 15-29.
[6] 刘世杰，刘卫东，林志勇，等. 连续旋转爆震波传播过程研究（I）：同向传播模式［J］. 推进技术， 2014， 35（1）： 138-144.
[7] 刘世杰，刘卫东，林志勇，等. 连续旋转爆震波传播过程研究（II）：双波对撞传播模式［J］. 推进技术， 2014， 35（2）： 269-275.
[8] Zhang Li-feng， Zhang Shu-jie， Wang Jian-ping. A Direct Simulation of Continuous Detonation Engine with the Naiver-Stokes Equations［R］. AIAA 2017-2323.
[9] Frolov S M， Dubrovskii A V， Ivanov V S. Three Dimensional Numerical Simulation of the Operation of the Rotating Detonation Chamber［J］. Combustion， Explosion and Shock Waves， 2012， 31（3）： 32-45.
[10] Frolov S M， Dubrovskii A V， Ivanov V S. Three Dimensional Numerical Simulation of the Operation of a Rotating Detonation Chamber with Separate Supply of Fuel and Oxidizer［J］. Combustion， Explosion and Shock Waves， 2013， 32（2）： 56-65.
[11] Strakey P A， Ferguson D， Sisler A， et al. Computationally Quantifying Loss Mechanisms in a Rotating Detonation Engine［R］. AIAA 2016-0900.
[12] Cocks P A， Holley A T， Greene C B， et al. Development of a High Fidelity RDE Simulation Capability［R］. AIAA 2015-1823.
[13] Cocks P A， Holley A T， Rankin B A， et al. High Fidelity Simulations of a Non-Premixed Rotating Detonation Engine［R］. AIAA 2016-0125.
[14] Gaillard T， Davidenko D， Dupoirieux F. Numerical Simulation of a Rotating Detonation with a Realistic Injector Designed for Separate Supply of Gaseous Hydrogen and Oxygen［J］. Acta Astronautica， 2017， 141： 64-78.
[15] 刘朋欣，李沁，张涵信. 初始点火条件引起的三维旋转爆震波单/双波传播模式［J］. 气体物理， 2018， 3（1）： 20-27.
[16] Li Qin， Liu Peng-xin， Zhang Han-xin. Further Investigations on the Interface Instability between Fresh Injections and Burnt Products in 2-D Rotating Detonation［J］. Computers and Fluids， 2018， 170： 261-272.
[17] Liu Pengxin， Li Qin， Huang Zhangfeng， et al. Interpretation of Wake Instability at Slip Line in Rotating Detonation［J］. International Journal of Computational Fluid Dynamics， 2018， 32（8-9）： 379-394.
[18] Menter F R. Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications［J］. AIAA Journal， 1994， 32（8）： 1598-1605.
[19] Shur M L， Spalart P R， Strelets M K， et al. A Hybrid RANS-LES Model with Delayed DES and Wall-Modeled LES Capabilities［J］. International Journal of Heat and Fluid Flow， 2008， 29： 1638-1649.
[20] Evans J S， Schexnayder C J. Influence of Chemical Kinetics and Unmixedness on Burning in Supersonic Hydrogen Flames［J］. AIAA Journal， 1980， 18（2）： 188-193.
[21] Jiang G S， Shu C W. Efficient Implementation of Weighted ENO Schemes［J］. Journal of Computational Physics， 1996， 126： 202-28.
[22] Li Qin， Liu Peng-xin， Zhang Han-xin. Piecewise Polynomial Mapping Method and Corresponding WENO Scheme with Improved Resolution［J］. Communications in Computational Physics， 2015， 18（5）： 1417-1444.
[23] Pan Z H， Fan B C. Wavelet Pattern and Self-Sustained Mechanism of Gaseous Detonation Rotating in a Coaxial Cylinder［J］. Combustion and Flame， 2011， 158： 2220-2228.

基于旋转爆震三维流场结构分析的计算模型对比研究

Computational Models Based on Analysis of Three Dimensional Flow Field Structures in Rotating Detonation

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 0

编辑推荐

Metrics

本文评价