超声速弹丸诱导爆轰波的无网格数值模拟

推进技术 ›› 2015, Vol. 36 ›› Issue (5): 664-670.

超声速弹丸诱导爆轰波的无网格数值模拟

吴　伟 ,许厚谦,王　亮,薛　锐,蔡峰峰

南京理工大学能源与动力工程学院，江苏南京 210094,南京理工大学能源与动力工程学院，江苏南京 210094,南京理工大学能源与动力工程学院，江苏南京 210094,南京理工大学能源与动力工程学院，江苏南京 210094,南京理工大学能源与动力工程学院，江苏南京 210094

发布日期:2021-08-15
作者简介:吴　伟(1985—)，男，博士生，研究领域为计算流体力学中的无网格方法。

Numerical Simulation of Denotation Wave Induced

School of Power Engineering，Nanjing University of Science & Technology，Nanjing 210094，China,School of Power Engineering，Nanjing University of Science & Technology，Nanjing 210094，China,School of Power Engineering，Nanjing University of Science & Technology，Nanjing 210094，China,School of Power Engineering，Nanjing University of Science & Technology，Nanjing 210094，China and School of Power Engineering，Nanjing University of Science & Technology，Nanjing 210094，China

Published:2021-08-15

摘要/Abstract

摘要： 为研究锥角对超声速弹丸诱导爆轰波形态的影响，发展了耦合有限速率化学反应模型的最小二乘显式无网格算法，其流体动力学采用含化学反应源项的多组分Euler方程建模，对流项和时间项分别采用多组分HLLC(Harten-Lax-van Leer-Contact)格式和四阶Runge-Kutta法计算。对尖劈诱导斜爆轰以及激波诱导燃烧流场进行了模拟，验证了算法的有效性。最后对等当量比甲烷/空气预混气体中，不同锥角弹丸诱导爆轰波流场进行了模拟，云图同实验阴影照片吻合较好，结果表明当锥角处于70°～100°时，易形成驻定斜爆轰波；锥角较小不利于可燃混合气体的点燃，仅能形成驻定斜激波；过大的锥角将导致爆轰波的脱体。

关键词: 无网格方；HLLC格式；有限速率反应模型；爆轰波

Abstract: The gridless method coupled with finite rate chemical model was developed to study the effects of conical angle on the structure of denotation shock induced by supersonic projectile. The fluid dynamics was modeled by multi-component Euler equations with chemical sources. The numerical flux was calculated by the multi-component HLLC (Harten-Lax-van Leer-Contact) scheme. A four-stage Runge-Kutta algorithm was employed to advance the equations in time. The flows of wedge-induced oblique denotation and shock-induced combustion were simulated to validate in respect of accuracy firstly. The denotation wave phenomena induced by supersonic projectile with different conical angles in the stoichiometric methane/air mixture were simulated subsequently. The contours agree well with the experimental shadowgraphs，and the results show that the oblique detonation wave is prone to be standing with conical angle at the range of 70°～100°，when the conical angle is not large enough ，it does not avail to ignite combustible mixture，only oblique shock is found in the flow field，contrarily，the denotation wave would be detached.

Key words: Gridless method；HLLC scheme；Finite rate chemical model；Detonation wave

吴　伟,许厚谦,王　亮,薛　锐,蔡峰峰. 超声速弹丸诱导爆轰波的无网格数值模拟[J]. 推进技术, 2015, 36(5): 664-670.

WU Wei,XU Hou-qian,WANG Liang,XUE Rui,CAI Feng-feng. Numerical Simulation of Denotation Wave Induced[J]. Journal of Propulsion Technology, 2015, 36(5): 664-670.

参考文献

[1] Nicholls J A, Dabora E K. Recent Results on Standing Detonation Waves[J]. Symposium (International) on Combustion, 1961, 8(1): 644-655.
[2] Menees. Analytical and Experimental Investigations of the Oblique Detonation Wave Engine Concept [R]. AGARR-CP- 479, 1991.
[3] 崔东明, 范宝春, 邢晓江. 驻定在高速弹丸上的斜爆轰波[J]. 爆炸与冲击, 2002, 22(3): 263-266.
[4] Green Brent S, Nelson H F. Effect of Inflow Pressure and Wall Heat Transfer on Shock-Induced Combustion / Detonation [R]. AIAA 96-0111.
[5] Figueira da Silva L F, Azevedo Joao L F, Heidi Korzenowski. Unstructured Adaptive Grid Flow Simulations of Inert and Reactive Gas Mixtures[J]. Journal of Computational Physics, 2000, 160: 522-540.
[6] HuiYuan Fan, Frank K Lu. Numerical Study of Reactive Flow past a Wedge in a Channel[R]. AIAA 2005-1168.
[7] 代淑兰. 复杂化学反应流并行数值模拟[D]. 南京: 南京理工大学, 2008.
[8] 代淑兰. 驻定斜爆轰波并行数值模拟[J]. 爆炸与冲击, 2006, 26(6): 572-576.
[9] 归明月, 范宝春. 尖劈诱导的斜爆轰波的精细结构及其影响因素[J]. 推进技术, 2012, 33(3): 490-494. (GUI Ming-yue, FAN Bao-chun. Fine Structure and Its Influence Factor of Wedge-Induced Oblique Detonation Waves[J]. Journal of Propulsion Technology, 2012, 33(3): 490-494.)
[10] 周平, 范宝春, 归明月. 可燃介质中飞行圆球诱导斜爆轰的流场结构[J]. 爆炸与冲击, 2012, 32(3): 278-282.
[11] 周星. 含动边界复杂非定常流动的无网格算法研究[D]. 南京: 南京理工大学, 2012.
[12] 孙迎丹, 王刚, 叶正寅. AUSM+-up格式在无网格算法中的推广[J]. 空气动力学学报, 2005, 23(4): 511-515.
[13] 吴伟, 许厚谦, 王亮, 等. AUFS格式在无网格方法中的应用[J].力学与实践, 2014, 36(1): 76-80.
[14] 马志华.自适应无网格及网格和无网格混合算法研究[D].南京: 南京航空航天大学, 2008.
[15] 马新建. 最小二乘无网格法的改进及其拓展应用研究[D]. 南京: 南京理工大学, 2012.
[16] Xiaowei Cai, Junjie Tan, Xinjian Ma, et al. Application of Hybrid Cartesian Grid and Gridless Approach to Moving Boundary Flow Problems[J]. International Journal for Numerical Method in Fluid, 2013,72:994-1013.
[17] Koh E P C, Tsai H M. Euler Solution Using Cartesian Grid with Least Squares Technique[R]. AIAA 2003-1120.
[18] 吴伟, 许厚谦. 激波诱导燃烧流场模拟的无网格算法[J]. 力学与实践, 2013, 35(6): 19-23.
[19] 卓长飞, 武晓松, 封锋, 等. 基于无网格算法的化学非平衡流数值模拟[J]. 推进技术, 2014, 35(2): 244-250. (ZHUO Chang-fei, WU Xiao-song, FENG Feng, et al. Numerical Simulation of Chemical Non-Equilibrium Flow Based on Gridless Algorithm[J]. Journal of Propulsion Technology, 2014, 35(2): 244-250.)
[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] 归明月, 范宝春. 尖劈诱导斜爆轰的胞格结构的数值研究[J]. 弹道学报, 2012, 24(2): 83-87.
[22] Moeljo Soetrisno, Scott T Imlay. Simulation of the Flow Field of a Ram Accelerator[R].AIAA 91-1915.(编辑:史亚红)　　　　　　　　　　　　　*　收稿日期:2014-03-04；修订日期:2014-05-13。作者简介:吴伟(1985—),男,博士生,研究领域为计算流体力学中的无网格方法。E-mail:wuwei8363@163.com