搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

高雷诺数下非混相Rayleigh-Taylor不稳定性的格子Boltzmann方法模拟

胡晓亮 梁宏 王会利

引用本文:
Citation:

高雷诺数下非混相Rayleigh-Taylor不稳定性的格子Boltzmann方法模拟

胡晓亮, 梁宏, 王会利

Lattice Boltzmann method simulations of the immiscible Rayleigh-Taylor instability with high Reynolds numbers

Hu Xiao-Liang, Liang Hong, Wang Hui-Li
PDF
HTML
导出引用
  • 本文采用相场格子Boltzmann方法研究了竖直微通道内中等Atwoods数流体的单模Rayleigh-Taylor不稳定性问题, 系统分析了雷诺数对相界面动力学行为以及扰动在各发展阶段演化规律的影响. 数值结果表明高雷诺数条件下, 不稳定性界面扰动的增长经历了四个不同的发展阶段, 包括线性增长阶段、饱和速度阶段、重加速阶段及混沌混合阶段. 在线性增长阶段, 我们计算获得的气泡与尖钉振幅符合线性稳定性理论, 并且线性增长率随着雷诺数的增加而增大. 在第二个阶段, 我们观察到气泡与尖钉将以恒定的速度增长, 获得的尖钉饱和速度略高于Goncharov经典势能模型的解析解[Phys. Rev. Lett. 2002 88 134502], 这归因于系统中产生了多个尺度的旋涡, 而涡之间的相互作用促进了尖钉的增长. 随着横向速度和纵向速度的差异扩大, 气泡和尖钉界面演化诱导产生的Kelvin–Helmholtz不稳定性逐渐增强, 从而流体混合区域出现许多不同层次的涡结构, 加速了气泡与尖钉振幅的演化速度, 并在演化后期阶段, 导致界面发生多层次卷起、剧烈变形、混沌破裂等行为, 最终形成了非常复杂的拓扑结构. 此外, 我们还统计了演化后期气泡与尖钉的无量纲加速度, 发现气泡和尖钉的振幅在后期呈现二次增长规律, 其增长率系数分别为0.045与0.233. 而在低雷诺条件下, 重流体在不稳定性后期以尖钉的形式向下运动而轻流体以气泡的形式向上升起. 在整个演化过程中, 界面变得足够光滑, 气泡与尖钉在后期的演化速度接近于常数, 未观察到后期的重加速与混沌混合阶段.
    In this paper, an advanced phase-field lattice Boltzmann method based on the multiple-relaxation-time collision model is used to simulate the immiscible single-mode Rayleigh-Taylor instability with a moderate Atwoods number in a long tube, and we systematically analyze the effect of the Reynolds number on the interfacial dynamics and the late-time development stages of interface disturbance. The highest Reynolds number in the current simulation reaches up to 10000. The numerical results show that the Reynolds number significantly affects the development of the instability. For high Reynolds numbers, the instability undergoes a sequence of different growth stages, which include the linear growth, saturated velocity growth, reacceleration, and chaotic mixing stages. In the linear growth stage, the developments of the bubble and spike conform to the classical linear growth theory, and it is shown that the growth rate increases with the Reynolds number. In the second stage, the bubble and spike evolve with the constant velocities, and the numerical prediction for spike velocity is found to be slightly larger than the solution of the potential flow theory proposed by Goncharov [Phys. Rev. Lett. 2002 88 134502], which can be attributed to the formation of vortices in the proximity of the spike tip. In addition, it is found that increasing the Reynolds number reduces the bubble saturated velocity, which then is smaller than the solution of the potential model. The nonlinear evolutions of the bubble and spike induce the Kelvin–Helmholtz instability, producing many vortex structures with different scales. Due to the interactions among the vortices, the instability eventually enters into the chaotic mixing stage, where the interfaces undergo the roll-up at multiple layers, sharp deformation, and chaotic breakup, forming a very complicated topology structure. Furthermore, we also measured the bubble and spike accelerations and find that the dimensionless values fluctuates around the constants of 0.045 and 0.233, indicating a mean quadratic growth. And for low Reynolds numbers, the heavy fluid will fall down in the form of the spike, and the interface in the whole process becomes very smooth without the appearances of the roll-up and vortices. The late-time evolutional stages such as the reacceleration and chaotic mixing cannot also be observed.
      通信作者: 梁宏, lianghongstefanie@163.com
    • 基金项目: 省部级-浙江省自然科学基金一般项目(LY19A020007)
      Corresponding author: Liang Hong, lianghongstefanie@163.com
    [1]

    Remington B A, Drake R P, Ryutov D D 2006 Rev. Mod. Phys. 78 755Google Scholar

    [2]

    Whitehead J A, Luther D S 1975 J. Geophys. Res. 80 705Google Scholar

    [3]

    Lindl J D, Amendt P, Berger R L, et al. 2004 Phys. Plasmas 11 339Google Scholar

    [4]

    Zhou C T, Yu M Y, He X T 2007 J. Appl. Phys. 101 103302Google Scholar

    [5]

    Rayleigh L 1883 Proc. London Math. Soc. 14 170

    [6]

    Taylor G I 1950 Proc. R. Soc. London 201 192Google Scholar

    [7]

    Mitchner M, Landshoff R K M 1964 Phys. Fluids 7 862Google Scholar

    [8]

    Chandrasekhar S 1961 Hydrodynamic and Hydromagnetic Stability (Oxford: Oxford University Press) pp1−653

    [9]

    Menikoff R, Mjolsness R C, Sharp D H, Zemach C 1977 Phys. Fluids 20 2000Google Scholar

    [10]

    Lewis D J 1950 Proc. R. Soc. London Ser. A 202 81Google Scholar

    [11]

    Sharp D H 1984 Physica D 12 3Google Scholar

    [12]

    Zhou Y 2017 Phys. Rep. 720−722 1Google Scholar

    [13]

    Zhou Y 2017 Phys. Rep. 723−725 1Google Scholar

    [14]

    Wei Y K, Wang Z D, Dou H S, Qian Y H 2017 Comput. Fluids 156 97Google Scholar

    [15]

    李德梅, 赖惠林, 许爱国, 等 2018 物理学报 67 080501Google Scholar

    Li D, Lai H, Xu A, et al. 2018 Acta Phys. Sin. 67 080501Google Scholar

    [16]

    Waddell J T, Niederhaus C E, Jacobs J W 2001 Phys. Fluids 13 1263Google Scholar

    [17]

    Goncharov V N 2002 Phys. Rev. Lett. 88 134502Google Scholar

    [18]

    Wilkinson J P, Jacobs J W 2007 Phys. Fluids 19 124102Google Scholar

    [19]

    He X Y, Zhang R Y, Chen S Y, Doolen G D 1999 Phys. Fluids 11 1143Google Scholar

    [20]

    Glimm J, Li X L, Lin A D 2002 Acta Math. Appl. Sin. 18 1

    [21]

    Celani A, Mazzino A, Ginanneschi P M, Vozella L 2009 J. Fluid Mech. 622 115Google Scholar

    [22]

    Ramaprabhu P, Dimonte G, Woodward P, et al. 2012 Phys. Fluids 24 074107Google Scholar

    [23]

    Wei T, Livescu D 2012 Phys. Rev. E 86 046405Google Scholar

    [24]

    Liang H, Shi B C, Guo Z L, Chai Z H 2014 Phys. Rev. E 89 053320Google Scholar

    [25]

    Liang H, Li Q X, Shi B C, Chai Z H 2016 Phys. Rev. E 93 033113Google Scholar

    [26]

    Hu Z X, Zhang Y S, Tian B L, He Z W, Li L 2019 Phys. Fluids 31 104108Google Scholar

    [27]

    郭照立, 郑楚光, 格子Boltzmann方法的原理及应用(北京: 科学出版社) 第1—250页

    Guo Z L, Zheng C G 2009 Theory and Applications of Lattice Boltzmann Method (Beijing: Science Press) pp1−250 (in Chinese)

    [28]

    Liang H, Chai Z H, Shi B C, Guo Z L, Zhang T 2014 Phys. Rev. E 90 063311Google Scholar

    [29]

    Liang H, Li Y, Chen J X, Xu J R 2019 Int. J. Heat Mass Tranfer 130 1189Google Scholar

    [30]

    Liang H, Shi B C, Chai Z H 2016 Phys. Rev. E 93 013308Google Scholar

    [31]

    Liang H, Xu J R, Chen J X, Chai Z H, Shi B C 2019 Appl. Math. Model. 73 487Google Scholar

    [32]

    娄钦, 李涛, 杨茉 2018 物理学报 67 234701Google Scholar

    Lou Q, Li T, Yang M 2018 Acta Phys. Sin. 67 234701Google Scholar

    [33]

    臧晨强, 娄钦 2017 物理学报 66 134701Google Scholar

    Zang C Q, Lou Q 2017 Acta Phys. Sin. 66 134701Google Scholar

    [34]

    梁宏, 柴振华, 施保昌 2016 物理学报 65 204701Google Scholar

    Liang H, Chai Z H, Shi B C 2016 Acta Phys. Sin. 65 204701Google Scholar

    [35]

    Liang H, Liu H H, Chai Z H, Shi B C 2019 Phys. Rev. E 99 063306Google Scholar

    [36]

    Lallemand P, Luo LS 2000 Phys. Rev. E 61 6546

    [37]

    Wei Y K, Wang Z D, Yang J F, Dou H S, Qian Y H 2015 Comput. Fluids 118 167Google Scholar

    [38]

    Wei Y K, Yang H, Lin Z, Wang Z D, Qian Y H 2018 Appl. Math. Comput. 339 556

    [39]

    Liang H, Xu J R, Chen J X, Wang H L, Chai Z H, Shi B C 2018 Phys. Rev. E 97 033309Google Scholar

    [40]

    Abarzhi S I, Gorobets A, Sreenivasan K R 2005 Phys. Fluids 17 081705Google Scholar

    [41]

    Sreenivasan K R 1984 Phys. Fluids 27 1048Google Scholar

  • 图 1  雷诺数对非混相RT不稳定性中相界面演化图案的影响 (a) $ Re = 10000 $; (b) $ Re = 2048 $; (c) $ Re = 50 $; (d) $ Re = 5 $

    Fig. 1.  The effect of the Reynolds number on the evolution of interfacial patterns in the immiscible RT instability: (a) $ Re = 10000 $; (b) $ Re = 2048 $; (c) $ Re = 50 $; (d) $ Re = 5 $.

    图 2  雷诺数对无量纲化的气泡与尖钉随时间演化振幅的影响

    Fig. 2.  The effect of the Reynolds number on the dimensionless bubble and spike amplitudes.

    图 3  雷诺数对无量纲化的气泡和尖钉演化速度的影响

    Fig. 3.  The effect of the Reynolds number on the dimensionless bubble and spike velocities.

    图 4  不同雷诺数下, 气泡和尖钉振幅在初始阶段的演化曲线, 其中数据点是统计结果, 实线则是拟合结果

    Fig. 4.  The curves of the early-time bubble and spike amplitudes with different Reynolds numbers, where the data points and solid lines are the statistical and fitting results.

    图 5  高雷诺数下, 气泡和尖钉的无量纲加速度演化曲线, 红色和蓝色实线分别为0.045和0.233

    Fig. 5.  The curves of dimensionless bubble and spike accelerations at a high Reynolds number, and the red and blue solid lines are 0.045 and 0.233.

  • [1]

    Remington B A, Drake R P, Ryutov D D 2006 Rev. Mod. Phys. 78 755Google Scholar

    [2]

    Whitehead J A, Luther D S 1975 J. Geophys. Res. 80 705Google Scholar

    [3]

    Lindl J D, Amendt P, Berger R L, et al. 2004 Phys. Plasmas 11 339Google Scholar

    [4]

    Zhou C T, Yu M Y, He X T 2007 J. Appl. Phys. 101 103302Google Scholar

    [5]

    Rayleigh L 1883 Proc. London Math. Soc. 14 170

    [6]

    Taylor G I 1950 Proc. R. Soc. London 201 192Google Scholar

    [7]

    Mitchner M, Landshoff R K M 1964 Phys. Fluids 7 862Google Scholar

    [8]

    Chandrasekhar S 1961 Hydrodynamic and Hydromagnetic Stability (Oxford: Oxford University Press) pp1−653

    [9]

    Menikoff R, Mjolsness R C, Sharp D H, Zemach C 1977 Phys. Fluids 20 2000Google Scholar

    [10]

    Lewis D J 1950 Proc. R. Soc. London Ser. A 202 81Google Scholar

    [11]

    Sharp D H 1984 Physica D 12 3Google Scholar

    [12]

    Zhou Y 2017 Phys. Rep. 720−722 1Google Scholar

    [13]

    Zhou Y 2017 Phys. Rep. 723−725 1Google Scholar

    [14]

    Wei Y K, Wang Z D, Dou H S, Qian Y H 2017 Comput. Fluids 156 97Google Scholar

    [15]

    李德梅, 赖惠林, 许爱国, 等 2018 物理学报 67 080501Google Scholar

    Li D, Lai H, Xu A, et al. 2018 Acta Phys. Sin. 67 080501Google Scholar

    [16]

    Waddell J T, Niederhaus C E, Jacobs J W 2001 Phys. Fluids 13 1263Google Scholar

    [17]

    Goncharov V N 2002 Phys. Rev. Lett. 88 134502Google Scholar

    [18]

    Wilkinson J P, Jacobs J W 2007 Phys. Fluids 19 124102Google Scholar

    [19]

    He X Y, Zhang R Y, Chen S Y, Doolen G D 1999 Phys. Fluids 11 1143Google Scholar

    [20]

    Glimm J, Li X L, Lin A D 2002 Acta Math. Appl. Sin. 18 1

    [21]

    Celani A, Mazzino A, Ginanneschi P M, Vozella L 2009 J. Fluid Mech. 622 115Google Scholar

    [22]

    Ramaprabhu P, Dimonte G, Woodward P, et al. 2012 Phys. Fluids 24 074107Google Scholar

    [23]

    Wei T, Livescu D 2012 Phys. Rev. E 86 046405Google Scholar

    [24]

    Liang H, Shi B C, Guo Z L, Chai Z H 2014 Phys. Rev. E 89 053320Google Scholar

    [25]

    Liang H, Li Q X, Shi B C, Chai Z H 2016 Phys. Rev. E 93 033113Google Scholar

    [26]

    Hu Z X, Zhang Y S, Tian B L, He Z W, Li L 2019 Phys. Fluids 31 104108Google Scholar

    [27]

    郭照立, 郑楚光, 格子Boltzmann方法的原理及应用(北京: 科学出版社) 第1—250页

    Guo Z L, Zheng C G 2009 Theory and Applications of Lattice Boltzmann Method (Beijing: Science Press) pp1−250 (in Chinese)

    [28]

    Liang H, Chai Z H, Shi B C, Guo Z L, Zhang T 2014 Phys. Rev. E 90 063311Google Scholar

    [29]

    Liang H, Li Y, Chen J X, Xu J R 2019 Int. J. Heat Mass Tranfer 130 1189Google Scholar

    [30]

    Liang H, Shi B C, Chai Z H 2016 Phys. Rev. E 93 013308Google Scholar

    [31]

    Liang H, Xu J R, Chen J X, Chai Z H, Shi B C 2019 Appl. Math. Model. 73 487Google Scholar

    [32]

    娄钦, 李涛, 杨茉 2018 物理学报 67 234701Google Scholar

    Lou Q, Li T, Yang M 2018 Acta Phys. Sin. 67 234701Google Scholar

    [33]

    臧晨强, 娄钦 2017 物理学报 66 134701Google Scholar

    Zang C Q, Lou Q 2017 Acta Phys. Sin. 66 134701Google Scholar

    [34]

    梁宏, 柴振华, 施保昌 2016 物理学报 65 204701Google Scholar

    Liang H, Chai Z H, Shi B C 2016 Acta Phys. Sin. 65 204701Google Scholar

    [35]

    Liang H, Liu H H, Chai Z H, Shi B C 2019 Phys. Rev. E 99 063306Google Scholar

    [36]

    Lallemand P, Luo LS 2000 Phys. Rev. E 61 6546

    [37]

    Wei Y K, Wang Z D, Yang J F, Dou H S, Qian Y H 2015 Comput. Fluids 118 167Google Scholar

    [38]

    Wei Y K, Yang H, Lin Z, Wang Z D, Qian Y H 2018 Appl. Math. Comput. 339 556

    [39]

    Liang H, Xu J R, Chen J X, Wang H L, Chai Z H, Shi B C 2018 Phys. Rev. E 97 033309Google Scholar

    [40]

    Abarzhi S I, Gorobets A, Sreenivasan K R 2005 Phys. Fluids 17 081705Google Scholar

    [41]

    Sreenivasan K R 1984 Phys. Fluids 27 1048Google Scholar

  • [1] 赖瑶瑶, 陈鑫梦, 柴振华, 施保昌. 基于格子Boltzmann方法的钉扎螺旋波反馈控制. 物理学报, 2024, 73(4): 040502. doi: 10.7498/aps.73.20231549
    [2] 刘程, 梁宏. 三相流体的轴对称格子 Boltzmann 模型及其在 Rayleigh-Plateau 不稳定性的应用. 物理学报, 2023, 72(4): 044701. doi: 10.7498/aps.72.20221967
    [3] 陈效鹏, 冯君鹏, 胡海豹, 杜鹏, 王体康. 基于格子Boltzmann方法的二维气泡群熟化过程模拟. 物理学报, 2022, 71(11): 110504. doi: 10.7498/aps.70.20212183
    [4] 陈效鹏, 冯君鹏, 胡海豹, 杜鹏, 王体康. 基于格子Boltzmann方法的二维汽泡群熟化过程模拟. 物理学报, 2022, (): . doi: 10.7498/aps.71.20212183
    [5] 马聪, 刘斌, 梁宏. 耦合界面张力的三维流体界面不稳定性的格子Boltzmann模拟. 物理学报, 2022, 71(4): 044701. doi: 10.7498/aps.71.20212061
    [6] 张恒, 任峰, 胡海豹. 基于格子Boltzmann方法的幂律流体二维顶盖驱动流转捩研究. 物理学报, 2021, 70(18): 184703. doi: 10.7498/aps.70.20210451
    [7] 黄皓伟, 梁宏, 徐江荣. 表面张力对高雷诺数Rayleigh-Taylor不稳定性后期增长的影响. 物理学报, 2021, 70(11): 114701. doi: 10.7498/aps.70.20201960
    [8] 李碧勇, 彭建祥, 谷岩, 贺红亮. 爆轰加载下高纯铜界面Rayleigh-Taylor不稳定性实验研究. 物理学报, 2020, 69(9): 094701. doi: 10.7498/aps.69.20191999
    [9] 李德梅, 赖惠林, 许爱国, 张广财, 林传栋, 甘延标. 可压流体Rayleigh-Taylor不稳定性的离散Boltzmann模拟. 物理学报, 2018, 67(8): 080501. doi: 10.7498/aps.67.20171952
    [10] 解文军, 滕鹏飞. 声悬浮过程的格子Boltzmann方法研究. 物理学报, 2014, 63(16): 164301. doi: 10.7498/aps.63.164301
    [11] 史冬岩, 王志凯, 张阿漫. 任意复杂流-固边界的格子Boltzmann处理方法. 物理学报, 2014, 63(7): 074703. doi: 10.7498/aps.63.074703
    [12] 袁永腾, 王立峰, 涂绍勇, 吴俊峰, 曹柱荣, 詹夏宇, 叶文华, 刘慎业, 江少恩, 丁永坤, 缪文勇. 掺杂对CH样品Rayleigh-Taylor不稳定性增长的影响. 物理学报, 2014, 63(23): 235203. doi: 10.7498/aps.63.235203
    [13] 夏同军, 董永强, 曹义刚. 界面张力对Rayleigh-Taylor不稳定性的影响. 物理学报, 2013, 62(21): 214702. doi: 10.7498/aps.62.214702
    [14] 曾建邦, 李隆键, 蒋方明. 气泡成核过程的格子Boltzmann方法模拟. 物理学报, 2013, 62(17): 176401. doi: 10.7498/aps.62.176401
    [15] 陶烨晟, 王立锋, 叶文华, 张广财, 张建成, 李英骏. 任意Atwood数Rayleigh-Taylor和 Richtmyer-Meshkov 不稳定性气泡速度研究. 物理学报, 2012, 61(7): 075207. doi: 10.7498/aps.61.075207
    [16] 曾建邦, 李隆键, 廖全, 陈清华, 崔文智, 潘良明. 格子Boltzmann方法在相变过程中的应用. 物理学报, 2010, 59(1): 178-185. doi: 10.7498/aps.59.178
    [17] 朱昌盛, 王智平, 荆 涛, 肖荣振. 二元合金微观偏析的相场法数值模拟. 物理学报, 2006, 55(3): 1502-1507. doi: 10.7498/aps.55.1502
    [18] 赵代平, 荆 涛, 柳百成. 相场方法模拟铝合金三维枝晶生长. 物理学报, 2003, 52(7): 1737-1742. doi: 10.7498/aps.52.1737
    [19] 李华兵, 黄乒花, 刘慕仁, 孔令江. 用格子Boltzmann方法模拟MKDV方程. 物理学报, 2001, 50(5): 837-840. doi: 10.7498/aps.50.837
    [20] 吕晓阳, 李华兵. 用格子Boltzmann方法模拟高雷诺数下的热空腔黏性流. 物理学报, 2001, 50(3): 422-427. doi: 10.7498/aps.50.422
计量
  • 文章访问数:  6411
  • PDF下载量:  114
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-10-05
  • 修回日期:  2019-10-29
  • 刊出日期:  2020-02-20

/

返回文章
返回