耦合界面张力的三维流体界面不稳定性的格子Boltzmann模拟  

Lattice Boltzmann simulation of three-dimensional fluid interfacial instability coupled with surface tension

作　　者：马聪 刘斌 梁宏 Ma Cong;Liu Bin;Liang Hong(School of Science,Hangzhou Dianzi University,Hangzhou 310018,China)

机构地区：[1]杭州电子科技大学理学院,杭州310018

出　　处：《物理学报》2022年第4期147-157,共11页Acta Physica Sinica

基　　金：国家自然科学基金(批准号:11972142)资助的课题.

摘　　要：采用介观格子Boltzmann方法模拟界面张力作用下三维流体界面的Rayleigh-Taylor(RT)不稳定性的增长过程,主要分析表面张力对流体界面动力学行为及尖钉和气泡后期增长的影响机制.首先发现三维RT不稳定性的发生存在临界表面张力(σc),其值随着流体Atwood数的增大而增大,且数值预测值与理论分析结果σc=(ρh−ρl)g/k^(2)一致.另外,随着表面张力的增大,不稳定性演化过程中界面卷吸程度和结构复杂性逐渐减弱,系统中界面破裂形成离散液滴的数目也显著减少.相界面的后期动力学行为也从非对称发展转向始终保持关于中轴线对称.尖钉与气泡振幅在表面张力较小时对其变化不显著,当表面张力增大到一定值后,可以有效地抑制尖钉与气泡振幅的增长.进一步发现,高雷诺数三维RT不稳定性在不同表面张力下均经历4个不同的发展阶段:线性阶段、饱和速度阶段、重加速和混沌混合阶段.尖钉与气泡在饱和速度阶段以近似恒定的速度增长,其渐进速度的值与修正的势流理论模型结果一致.受非线性Kelvin-Helmholtz旋涡的剪切作用,尖钉与气泡随后的增长被加速,导致在重加速阶段的演化速度超过势流模型的解析解.重加速阶段不能持续发展下去,尖钉与气泡在不稳定性后期的增长速度会随时间上下波动,这表明不稳定性的演化进入了混沌混合阶段.通过数值分析,证实了三维RT不稳定性在后期的混沌混合阶段具有二次增长的规律,并且尖钉与气泡增长率总体上随着表面张力的增大而逐渐减少.In this paper,the development of three-dimensional fluid interfacial Rayleigh-Taylor(RT)instability coupled with the surface tension was numerically studied using the mesoscopic lattice Boltzmann method.We mainly analyzed the influence of surface tension on fluid interfacial dynamics and spike/bubble late-time growth.The numerical experiments show that there exists the critical surface tension(σc)in the three-dimensional RT instability,above which the RT phenomenon does not appear and below which it would take place.It is found that the critical surface tension increases with the fluid Atwood number and the corresponding numerical predictions show good agreements with those of the theoretical analysisσc=(ρh−ρl)g/k^(2).In addition,we can find that increasing surface tension reduces the roll-up of the interface and the complexity of interfacial structure,also preventing the breakup of the interface into the individual droplets.The late-time dynamics of phase interface change from the asymmetric development to the symmetry with respect to the middle axis.When the surface tension is sufficiently low,the spike and bubble amplitudes almost no longer change with it,and further increasing the surface tension can slow down the growth of the spike and bubble amplitudes.Furthermore,we can observe that the development of the high-Reynolds-number RT instability under different surface tensions can also be divided into four distinct stages,including the linear growth,saturated velocity growth,reacceleration,and chaotic mixing.The spike and bubble grow with approximately constant velocities at the saturated stage and their asymptotic values are consistent with those of the modified potential flow theory.In the following,the spike and bubble driven by the increasing Kelvin-Helmholtz vortices are accelerated such that their evolutional velocities exceed the solutions of the potential flow model at the reacceleration stage.The reacceleration stage cannot last infinitely and the spike and bubble velocities at the late time fl

关 键 词：RAYLEIGH-TAYLOR不稳定性 格子BOLTZMANN方法 界面张力 雷诺数 湍流混合 

分 类 号：O35[理学—流体力学]