北京航空航天大学学报 ›› 2016, Vol. 42 ›› Issue (10): 2059-2068.doi: 10.13700/j.bh.1001-5965.2015.0667

• 论文 • 上一篇    下一篇

Rayleigh-Taylor不稳定性非线性特性的数值研究

康宁, 黎一锴, 何旭   

  1. 北京理工大学 机械与车辆学院, 北京 100081
  • 收稿日期:2015-10-15 出版日期:2016-10-20 发布日期:2016-04-11
  • 通讯作者: 黎一锴,Tel.:13811862628,E-mail:liyikai@bit.edu.cn E-mail:liyikai@bit.edu.cn
  • 作者简介:康宁,男,博士研究生。主要研究方向:内燃机流动与燃烧。E-mail:46054832@qq.com;黎一锴,男,博士,讲师。主要研究方向:气-液两相流的数值模拟。Tel.:13811862628,E-mail:liyikai@bit.edu.cn
  • 基金资助:
    国家自然科学基金(51476011);北京理工大学青年教师学术启动计划(3030012261599)

Numerical study on nonlinear characteristics of Rayleigh-Taylor instability

KANG Ning, LI Yikai, HE Xu   

  1. School of Mechanical Engineering, Beijing Institute of Technology, Beijing 100081, China
  • Received:2015-10-15 Online:2016-10-20 Published:2016-04-11
  • Supported by:
    National Natural Science Foundation of China (51476011); Beijing Institute of Technology Research Fund Program for Young Scholars (3030012261599)

摘要: 以往对于单模态Rayleigh-Taylor(RT)不稳定性非线性特性的研究主要集中于推导和测量恒定的气泡推进速度上,而缺乏对液态尖钉区域非线性动力学特性的详细分析。采用耦合的Level-Set和Volume-of-Fluid(CLSVOF)界面捕捉方法对单模态RT不稳定性的发展过程进行了精确的数值模拟,并利用模拟得到的压力场和速度场信息对RT不稳定性非线性发展阶段的稳态动力学特性进行了分析。模拟结果表明,在液态尖钉根部由于惯性力作用而引起的水平冲击流会在此处形成一个局部最大压力点,由于此处惯性力与压强梯度的平衡,位于最大压力点附近的流动最终将达到稳态。通过理论分析,确定了此稳态流动中各稳态特征参数与初始扰动波长、惯性加速度之间的关系。这些特征参数的确定有助于将经典低速射流的相关理论扩展应用到RT不稳定性诱导雾化的研究领域。

关键词: Rayleigh-Taylor不稳定性, 非线性, 稳态, CLSVOF方法, 数值模拟

Abstract: The research on the nonlinear dynamics of Rayleigh-Taylor (RT) before mainly focused on deducing and measuring the constant penetration velocity of the bubble and had little detailed analysis of the nonlinear dynamic characteristics in the liquid spike region. An accurate numerical simulation of the single-mode RT instability was carried out based on the coupled Level-Set and Volume-of-Fluid (CLSVOF) interface capturing method. The detailed information on the pressure fields and velocity fields was obtained. In addition, the steady-state dynamic characteristics in the nonlinear development stage were analyzed. Simulation results show that a local maximum pressure point which is caused by the horizontal impinging flow with the action of inertial force appears at the root of the spike. The dependence of the different characteristic parameters of the steady flow on the initial perturbation wavelength and the inertial acceleration is determined. This work may extend the relevant classical theories of the low speed jet to the RT instability inducing atomization field.

Key words: Rayleigh-Taylor instability, nonlinearity, steady state, CLSVOF method, numerical simulation

中图分类号: 


版权所有 © 《北京航空航天大学学报》编辑部
通讯地址:北京市海淀区学院路37号 北京航空航天大学学报编辑部 邮编:100191 E-mail:jbuaa@buaa.edu.cn
本系统由北京玛格泰克科技发展有限公司设计开发