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

摘要

采用介观格子Boltzmann方法模拟界面张力作用下三维流体界面的Rayleigh-Taylor(RT)不稳定性的增长过程,主要分析表面张力对流体界面动力学行为及尖钉和气泡后期增长的影响机制.首先发现三维RT不稳定性的发生存在临界表面张力(σc),其值随着流体Atwood数的增大而增大,且数值预测值与理论分析结果σc=(ρh−ρl)g/k^(2)一致.另外,随着表面张力的增大,不稳定性演化过程中界面卷吸程度和结构复杂性逐渐减弱,系统中界面破裂形成离散液滴的数目也显著减少.相界面的后期动力学行为也从非对称发展转向始终保持关于中轴线对称.尖钉与气泡振幅在表面张力较小时对其变化不显著,当表面张力增大到一定值后,可以有效地抑制尖钉与气泡振幅的增长.进一步发现,高雷诺数三维RT不稳定性在不同表面张力下均经历4个不同的发展阶段:线性阶段、饱和速度阶段、重加速和混沌混合阶段.尖钉与气泡在饱和速度阶段以近似恒定的速度增长,其渐进速度的值与修正的势流理论模型结果一致.受非线性Kelvin-Helmholtz旋涡的剪切作用,尖钉与气泡随后的增长被加速,导致在重加速阶段的演化速度超过势流模型的解析解.重加速阶段不能持续发展下去,尖钉与气泡在不稳定性后期的增长速度会随时间上下波动,这表明不稳定性的演化进入了混沌混合阶段.通过数值分析,证实了三维RT不稳定性在后期的混沌混合阶段具有二次增长的规律,并且尖钉与气泡增长率总体上随着表面张力的增大而逐渐减少.