欢迎访问过程工程学报, 今天是

过程工程学报 ›› 2022, Vol. 22 ›› Issue (1): 41-49.DOI: 10.12034/j.issn.1009-606X.221057

• 研究论文 • 上一篇    下一篇

平面激波诱导斜界面失稳的大涡模拟

管永康,刘仪*,章利特   

  1. 浙江理工大学机械与自动控制学院, 浙江 杭州 310018
  • 收稿日期:2021-02-22 修回日期:2021-03-18 出版日期:2022-01-28 发布日期:2022-01-28
  • 通讯作者: 刘仪 Yi.liu@zstu.edu.cn
  • 作者简介:管永康(1993-),男,安徽省广德市人,硕士研究生,流体机械及工程专业,E-mail: davidgyk1128@163.com;刘仪,通讯联系人,E-mail: Yi.liu@zstu.edu.cn
  • 基金资助:
    浙江理工大学科研启动基金

Large eddy simulation of inclined interface instability induced by planar shock wave

Yongkang GUAN,  Yi LIU*,  Lite ZHANG   

  1. Faculty of Mechanical Engineering & Automation, Zhejiang Sci-Tech University, Hangzhou, Zhejiang 310018, China
  • Received:2021-02-22 Revised:2021-03-18 Online:2022-01-28 Published:2022-01-28

摘要: 近年来国内外学者从多个角度深入地研究了经典平面激波与相界面之间的相互作用规律,但针对平面激波与气-液斜界面的研究仍不够系统和全面,探究平面激波驱动气-液斜界面失稳过程的演化规律在超声速燃烧和惯性约束聚变等工程应用中具有重要的实际意义。本工作结合VOF模型和大涡数值模拟方法,研究了平面激波诱导气-液两相斜界面变形和湍流混合的现象,分析了入射激波马赫数、入射激波初始倾角和斜界面倾角等因素对界面失稳的影响规律。结果表明,入射激波马赫数大小是界面失稳现象发展的主要影响因素,增大激波马赫数可明显地增强界面变形和湍流混合现象的发展程度。入射激波初始倾角和斜界面倾角对界面失稳发展也有一定影响。随着入射激波和斜界面倾角增大,相同时刻时的湍流混合区宽度相应增加。

关键词: VOF模型, 大涡模拟, 界面变形, 湍流混合, 气-液两相斜界面

Abstract: Interfacial instability and its induced interfacial deformation and turbulent instability are widely used. In medium characteristic scale supersonic combustion, the interfacial instability enhances the mixing of fuel and combustion aids. In small characteristic scale inertial confinement fusion, the turbulent mixing induced by interfacial instability dilutes and cools the fuel twice, thus reducing the reaction rate and even causing ignition failure. In recent years, scholars at home and abroad have deeply studied the interaction between the classical planar shock wave and the phase interface from many angles, but the research on the interaction between the shock wave and the gas-liquid oblique interface is still not systematic and comprehensive. It is of great significance to explore the evolution of shock driving gas-liquid interface instability in practical engineering applications such as supersonic combustion and inertial confinement fusion. In the process of numerical simulation, VOF model is used, combined with large eddy simulation method and appropriate boundary conditions to study the evolution process of shock induced gas-liquid two-phase oblique interface deformation and turbulent mixing in two-dimensional plane. The influence of Mach number of incident shock, initial inclination angle of shock and inclined interface inclination angle on interface instability is analyzed. The results show that the Mach number of incident shock is the most important factor in the development of interface instability. Under the same conditions, compared with the other two factors, increasing the shock Mach number can significantly increase the interface deformation and the development of turbulent mixing. At the same time, the number of serrated structures on the interface increases, the forming speed is obviously accelerated, and the width of turbulent mixing zone increases obviously. In addition, with the increase of the initial angle of shock wave and inclined interface angle, the number and forming speed of serrated structure on the interface increase at the same time, and the width of turbulent mixing zone also increases.

Key words: VOF model, large eddy simulation, interface deformation, turbulent mixing, gas-liquid two-phase oblique interface