-
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.
-
Keywords:
- Lattice Boltzmann method /
- Rayleigh-Taylor instability /
- phase field /
- Reynolds number /
- quadratic growth
[1] Remington B A, Drake R P, Ryutov D D 2006 Rev. Mod. Phys. 78 755
Google Scholar
[2] Whitehead J A, Luther D S 1975 J. Geophys. Res. 80 705
Google Scholar
[3] Lindl J D, Amendt P, Berger R L, et al. 2004 Phys. Plasmas 11 339
Google Scholar
[4] Zhou C T, Yu M Y, He X T 2007 J. Appl. Phys. 101 103302
Google Scholar
[5] Rayleigh L 1883 Proc. London Math. Soc. 14 170
[6] Taylor G I 1950 Proc. R. Soc. London 201 192
Google Scholar
[7] Mitchner M, Landshoff R K M 1964 Phys. Fluids 7 862
Google 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 2000
Google Scholar
[10] Lewis D J 1950 Proc. R. Soc. London Ser. A 202 81
Google Scholar
[11] Sharp D H 1984 Physica D 12 3
Google Scholar
[12] Zhou Y 2017 Phys. Rep. 720−722 1
Google Scholar
[13] Zhou Y 2017 Phys. Rep. 723−725 1
Google Scholar
[14] Wei Y K, Wang Z D, Dou H S, Qian Y H 2017 Comput. Fluids 156 97
Google Scholar
[15] 李德梅, 赖惠林, 许爱国, 等 2018 物理学报 67 080501
Google Scholar
Li D, Lai H, Xu A, et al. 2018 Acta Phys. Sin. 67 080501
Google Scholar
[16] Waddell J T, Niederhaus C E, Jacobs J W 2001 Phys. Fluids 13 1263
Google Scholar
[17] Goncharov V N 2002 Phys. Rev. Lett. 88 134502
Google Scholar
[18] Wilkinson J P, Jacobs J W 2007 Phys. Fluids 19 124102
Google Scholar
[19] He X Y, Zhang R Y, Chen S Y, Doolen G D 1999 Phys. Fluids 11 1143
Google 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 115
Google Scholar
[22] Ramaprabhu P, Dimonte G, Woodward P, et al. 2012 Phys. Fluids 24 074107
Google Scholar
[23] Wei T, Livescu D 2012 Phys. Rev. E 86 046405
Google Scholar
[24] Liang H, Shi B C, Guo Z L, Chai Z H 2014 Phys. Rev. E 89 053320
Google Scholar
[25] Liang H, Li Q X, Shi B C, Chai Z H 2016 Phys. Rev. E 93 033113
Google Scholar
[26] Hu Z X, Zhang Y S, Tian B L, He Z W, Li L 2019 Phys. Fluids 31 104108
Google 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 063311
Google Scholar
[29] Liang H, Li Y, Chen J X, Xu J R 2019 Int. J. Heat Mass Tranfer 130 1189
Google Scholar
[30] Liang H, Shi B C, Chai Z H 2016 Phys. Rev. E 93 013308
Google Scholar
[31] Liang H, Xu J R, Chen J X, Chai Z H, Shi B C 2019 Appl. Math. Model. 73 487
Google Scholar
[32] 娄钦, 李涛, 杨茉 2018 物理学报 67 234701
Google Scholar
Lou Q, Li T, Yang M 2018 Acta Phys. Sin. 67 234701
Google Scholar
[33] 臧晨强, 娄钦 2017 物理学报 66 134701
Google Scholar
Zang C Q, Lou Q 2017 Acta Phys. Sin. 66 134701
Google Scholar
[34] 梁宏, 柴振华, 施保昌 2016 物理学报 65 204701
Google Scholar
Liang H, Chai Z H, Shi B C 2016 Acta Phys. Sin. 65 204701
Google Scholar
[35] Liang H, Liu H H, Chai Z H, Shi B C 2019 Phys. Rev. E 99 063306
Google 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 167
Google 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 033309
Google Scholar
[40] Abarzhi S I, Gorobets A, Sreenivasan K R 2005 Phys. Fluids 17 081705
Google Scholar
[41] Sreenivasan K R 1984 Phys. Fluids 27 1048
Google Scholar
-
图 1 雷诺数对非混相RT不稳定性中相界面演化图案的影响 (a)
$ Re = 10000 $ ; (b)$ Re = 2048 $ ; (c)$ Re = 50 $ ; (d)$ Re = 5 $ Figure 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 雷诺数对无量纲化的气泡与尖钉随时间演化振幅的影响
Figure 2. The effect of the Reynolds number on the dimensionless bubble and spike amplitudes.
图 3 雷诺数对无量纲化的气泡和尖钉演化速度的影响
Figure 3. The effect of the Reynolds number on the dimensionless bubble and spike velocities.
图 4 不同雷诺数下, 气泡和尖钉振幅在初始阶段的演化曲线, 其中数据点是统计结果, 实线则是拟合结果
Figure 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
Figure 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 755
Google Scholar
[2] Whitehead J A, Luther D S 1975 J. Geophys. Res. 80 705
Google Scholar
[3] Lindl J D, Amendt P, Berger R L, et al. 2004 Phys. Plasmas 11 339
Google Scholar
[4] Zhou C T, Yu M Y, He X T 2007 J. Appl. Phys. 101 103302
Google Scholar
[5] Rayleigh L 1883 Proc. London Math. Soc. 14 170
[6] Taylor G I 1950 Proc. R. Soc. London 201 192
Google Scholar
[7] Mitchner M, Landshoff R K M 1964 Phys. Fluids 7 862
Google 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 2000
Google Scholar
[10] Lewis D J 1950 Proc. R. Soc. London Ser. A 202 81
Google Scholar
[11] Sharp D H 1984 Physica D 12 3
Google Scholar
[12] Zhou Y 2017 Phys. Rep. 720−722 1
Google Scholar
[13] Zhou Y 2017 Phys. Rep. 723−725 1
Google Scholar
[14] Wei Y K, Wang Z D, Dou H S, Qian Y H 2017 Comput. Fluids 156 97
Google Scholar
[15] 李德梅, 赖惠林, 许爱国, 等 2018 物理学报 67 080501
Google Scholar
Li D, Lai H, Xu A, et al. 2018 Acta Phys. Sin. 67 080501
Google Scholar
[16] Waddell J T, Niederhaus C E, Jacobs J W 2001 Phys. Fluids 13 1263
Google Scholar
[17] Goncharov V N 2002 Phys. Rev. Lett. 88 134502
Google Scholar
[18] Wilkinson J P, Jacobs J W 2007 Phys. Fluids 19 124102
Google Scholar
[19] He X Y, Zhang R Y, Chen S Y, Doolen G D 1999 Phys. Fluids 11 1143
Google 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 115
Google Scholar
[22] Ramaprabhu P, Dimonte G, Woodward P, et al. 2012 Phys. Fluids 24 074107
Google Scholar
[23] Wei T, Livescu D 2012 Phys. Rev. E 86 046405
Google Scholar
[24] Liang H, Shi B C, Guo Z L, Chai Z H 2014 Phys. Rev. E 89 053320
Google Scholar
[25] Liang H, Li Q X, Shi B C, Chai Z H 2016 Phys. Rev. E 93 033113
Google Scholar
[26] Hu Z X, Zhang Y S, Tian B L, He Z W, Li L 2019 Phys. Fluids 31 104108
Google 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 063311
Google Scholar
[29] Liang H, Li Y, Chen J X, Xu J R 2019 Int. J. Heat Mass Tranfer 130 1189
Google Scholar
[30] Liang H, Shi B C, Chai Z H 2016 Phys. Rev. E 93 013308
Google Scholar
[31] Liang H, Xu J R, Chen J X, Chai Z H, Shi B C 2019 Appl. Math. Model. 73 487
Google Scholar
[32] 娄钦, 李涛, 杨茉 2018 物理学报 67 234701
Google Scholar
Lou Q, Li T, Yang M 2018 Acta Phys. Sin. 67 234701
Google Scholar
[33] 臧晨强, 娄钦 2017 物理学报 66 134701
Google Scholar
Zang C Q, Lou Q 2017 Acta Phys. Sin. 66 134701
Google Scholar
[34] 梁宏, 柴振华, 施保昌 2016 物理学报 65 204701
Google Scholar
Liang H, Chai Z H, Shi B C 2016 Acta Phys. Sin. 65 204701
Google Scholar
[35] Liang H, Liu H H, Chai Z H, Shi B C 2019 Phys. Rev. E 99 063306
Google 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 167
Google 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 033309
Google Scholar
[40] Abarzhi S I, Gorobets A, Sreenivasan K R 2005 Phys. Fluids 17 081705
Google Scholar
[41] Sreenivasan K R 1984 Phys. Fluids 27 1048
Google Scholar
DownLoad:
Catalog
Metrics
- Abstract views: 8049
- PDF Downloads: 126
- Cited By: 0
