-
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 (
$\sigma_{\rm{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${\sigma_{\rm{c}}}= {{({{\rho_{\rm{h}}}-{\rho_{\rm{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 fluctuate with time, implying that the growth of the RT instability enters into the chaotic mixing stage. By numerical analysis, we demonstrate that the three-dimensional RT instability at the chaotic mixing stage has a quadratic growth and also report that the spike and bubble growth rates decrease with the surface tension in general.-
Keywords:
- Rayleigh-Taylor instability /
- lattice Boltzmann method /
- surface tension /
- Reynolds number /
- turbulent mixing
[1] Zhou Y 2017 Phys. Rep. 720
[2] Boffetta G, Mazzino A 2017 Annu. Rev. Fluid Mech. 49 119
Google Scholar
[3] Rayleigh L 1883 Proc. London Math. Soc. 14 170
[4] Taylor G I 1950 Proc. R. Soc. London Ser. A 201 192
[5] Lewis D J 1950 Proc. R. Soc. London Ser. A 202 81
Google Scholar
[6] Waddell J T, Niederhaus C E, Jacobs J W 2001 Phys. Fluids 13 1263
Google Scholar
[7] Goncharov V N 2002 Phys. Rev. Lett. 88 134502
Google Scholar
[8] Sohn S I 2009 Phys. Rev. E 80 055302(R
Google Scholar
[9] Glimm J, Li X L, Lin A D 2002 Acta Math. Appl. Sin. 18 1
[10] Wilkinson J P, Jacobs J W 2007 Phys. Fluids 19 124102
Google Scholar
[11] Ramaprabhu P, Dimonte G, Young Y N, Calder A C, Fryxell B 2006 Phys. Rev. E 74 066308
Google Scholar
[12] Bian X, Aluie H, Zhao D X, Zhang H S, Livescu D 2020 Physica D 403 132250
Google Scholar
[13] Ramaprabhu P, Dimonte G, Woodward P, Fryer C, Rochefeller G, Muthuraman K, Lin P H, Jayaraj J 2012 Phys. Fluids 24 074107
Google Scholar
[14] Wei T, Livescu D 2012 Phys. Rev. E 86 046405
Google Scholar
[15] Hu Z X, Zhang Y S, Tian B L, He Z W, Li L 2019 Phys. Fluids 31 104108
Google Scholar
[16] 李德梅, 赖惠林, 许爱国, 张广财, 林传栋, 甘延标 2018 67 080501
Google Scholar
Li D M, Lai H L, Xu A G, Zhang G C, Lin C D, Gan Y B 2018 Acta Phys. Sin. 67 080501
Google Scholar
[17] Liang H, Li Q X, Shi B C, Chai Z H 2016 Phys. Rev. E 93 033113
Google Scholar
[18] 胡晓亮, 梁宏, 王会利 2020 69 044701
Google Scholar
Hu X L, Liang H, Wang H L 2020 Acta Phys. Sin. 69 044701
Google Scholar
[19] Liang H, Xia Z H, Huang H W 2021 Phys. Fluids 33 082103
Google Scholar
[20] Daly B J 1969 Phys. Fluids 12 1340
Google Scholar
[21] Zhang R Y, He X Y, Chen S Y 2000 Comput. Phys. Commun. 129 121
Google Scholar
[22] Young Y N, Ham F E 2006 J. Turbul. 7 71
Google Scholar
[23] Matsuoka C 2009 Phys. Fluids 21 092107
Google Scholar
[24] Sohn S I 2009 Physical Review E 80 055302(R
[25] 夏同军, 董永强, 曹义刚 2013 62 214702
Google Scholar
Xia T J, Dong Y Q, Cao Y G 2013 Acta. Phys. Sin. 62 214702
Google Scholar
[26] 黄皓伟, 梁宏, 徐江荣 2021 70 114701
Google Scholar
Huang H W, Liang H, Xu J R 2021 Acta Phys. Sin. 70 114701
Google Scholar
[27] Guo Z L, Shu C 2013 Lattice Boltzmann Method and its Applications in Engineering (Singapore: World Scientific Publishing Company) pp239–284
[28] Wang H L, Yuan X L, Liang H, Chai Z H, Shi B C 2019 Capillarity 2 33
Google Scholar
[29] Liang H, Hu X L, Huang X F, Xu J R 2019 Phys. Fluids 31 112104
Google Scholar
[30] Liang H, Shi B C, Guo Z L, Chai Z H 2014 Phys. Rev. E 89 053320
Google Scholar
[31] Liang H, Shi B C, Chai Z H 2017 Comput. Math. Appl. 73 1524
Google Scholar
[32] d’Humieres D, Ginzburg I, Krafczyk M, Lallemand P, Luo L S 2002 Philos. Trans. R. Soc. London, Ser. A 360 437
Google Scholar
[33] He X Y, Chen S Y, Zhang R Y 1999 J. Comput. Phys. 152 642
Google Scholar
[34] Abarzhi S I, Gorobets A, Sreenivasan K R 2005 Phys. Fluids 17 081705
Google Scholar
[35] Olson D H, Jacobs J W 2009 Phys. Fluids 21 034103
Google Scholar
-
图 1 三维单模RT不稳定性问题的示意图
Figure 1. Schematic of three-dimensional single-mode RT instability problem
图 2 不同Atwood数下三维RT不稳定性的临界表面张力
Figure 2. Critical surface tensions of three-dimensional RT instability at various Atwood numbers
图 3 表面张力对高雷诺数三维单模RT不稳定性相界面演化的影响,
$ Re = 5000 $ ,$ A_t = 0.1 $ (a)$ \sigma = 1\times10^{-6} $ ; (b)$\sigma = 1\times $ $ 10^{ - 5}$ ; (c)$ \sigma = 2\times10^{ - 4} $ ; (d)$ \sigma = 5\times10^{ - 4} $ Figure 3. Effect of the surface tension on the evolution of phase interface in three-dimensional single-mode RT instability with high Reynolds number,
$ Re = 5000 $ ,$ A_t = 0.1 $ : (a)$ \sigma = 1\times10^{-6} $ ; (b)$ \sigma = 1\times10^{ - 5} $ ; (c)$ \sigma = 2\times10^{ - 4} $ ; (d)$ \sigma = 5\times10^{ - 4} $ .
图 4 表面张力对高雷诺数三维单模RT不稳定性的对角平面(x = y平面)的相界面演化的影响,
$ Re = 5000 $ ,$ A_t = 0.1 $ (a)$\sigma = $ $ 1\times10^{ - 6}$ ; (b)$ \sigma = 1\times10^{ - 5} $ ; (c)$ \sigma = 2\times10^{ - 4} $ ; (d)$ \sigma = 5\times10^{ - 4} $ Figure 4. Effect of the surface tension on the evolution of phase interface in the diagonal plane of three-dimensional single-mode RT instability with high Reynolds number,
$ Re = 5000 $ ,$ A_t = 0.1 $ : (a)$ \sigma = 1\times10^{ - 6} $ ; (b)$ \sigma = 1\times10^{ - 5} $ ; (c)$ \sigma = 2\times10^{ - 4} $ ; (d)$\sigma = $ $ 5\times10^{ - 4}$ .
图 5 表面张力对随时间演化的尖钉(左)和气泡(右)振幅的影响
Figure 5. Influence of surface tension on the time evolutions of spike (left) and bubble (right) amplitudes.
图 6 表面张力对随时间演化的尖钉(左)和气泡(右)增长速度的影响. 黑色和蓝色虚线分别表示修正的势能模型所预测尖钉与气泡渐进速度在
$ \sigma = 1\times10^{ - 5} $ $ \sigma = 5\times10^{ - 4} $ 时的解析解Figure 6. Influence of surface tension on the time evolutions of spike (left) and bubble (right) growth velocities. The black and blue dotted lines represent the analytical solutions of the spike and bubble asymptotic velocities from the modified potential flow model at
$ \sigma = 1\times10^{ - 5} $ and$ \sigma = 5\times10^{ - 4} $ .
图 7 不同界面张力下三维RT不稳定性中尖钉(左)与气泡(右)振幅的开方
$ \sqrt{h_{{\rm{s, b}}}} $ 与演化时间$ \sqrt{gA_t}t $ 的关系曲线, 实线表示在后期混沌混合阶段的线性拟合结果Figure 7. Relations between the square roots of the spike and bubble amplitudes
$ \sqrt{h_{{\rm{s, b}}}} $ and the time$ \sqrt{g{A_t}}t $ in the three-dimensional RT instability with different surface tensions. The solid lines represent the linear fitting results at the late-time chaotic mixing stage. -
[1] Zhou Y 2017 Phys. Rep. 720
[2] Boffetta G, Mazzino A 2017 Annu. Rev. Fluid Mech. 49 119
Google Scholar
[3] Rayleigh L 1883 Proc. London Math. Soc. 14 170
[4] Taylor G I 1950 Proc. R. Soc. London Ser. A 201 192
[5] Lewis D J 1950 Proc. R. Soc. London Ser. A 202 81
Google Scholar
[6] Waddell J T, Niederhaus C E, Jacobs J W 2001 Phys. Fluids 13 1263
Google Scholar
[7] Goncharov V N 2002 Phys. Rev. Lett. 88 134502
Google Scholar
[8] Sohn S I 2009 Phys. Rev. E 80 055302(R
Google Scholar
[9] Glimm J, Li X L, Lin A D 2002 Acta Math. Appl. Sin. 18 1
[10] Wilkinson J P, Jacobs J W 2007 Phys. Fluids 19 124102
Google Scholar
[11] Ramaprabhu P, Dimonte G, Young Y N, Calder A C, Fryxell B 2006 Phys. Rev. E 74 066308
Google Scholar
[12] Bian X, Aluie H, Zhao D X, Zhang H S, Livescu D 2020 Physica D 403 132250
Google Scholar
[13] Ramaprabhu P, Dimonte G, Woodward P, Fryer C, Rochefeller G, Muthuraman K, Lin P H, Jayaraj J 2012 Phys. Fluids 24 074107
Google Scholar
[14] Wei T, Livescu D 2012 Phys. Rev. E 86 046405
Google Scholar
[15] Hu Z X, Zhang Y S, Tian B L, He Z W, Li L 2019 Phys. Fluids 31 104108
Google Scholar
[16] 李德梅, 赖惠林, 许爱国, 张广财, 林传栋, 甘延标 2018 67 080501
Google Scholar
Li D M, Lai H L, Xu A G, Zhang G C, Lin C D, Gan Y B 2018 Acta Phys. Sin. 67 080501
Google Scholar
[17] Liang H, Li Q X, Shi B C, Chai Z H 2016 Phys. Rev. E 93 033113
Google Scholar
[18] 胡晓亮, 梁宏, 王会利 2020 69 044701
Google Scholar
Hu X L, Liang H, Wang H L 2020 Acta Phys. Sin. 69 044701
Google Scholar
[19] Liang H, Xia Z H, Huang H W 2021 Phys. Fluids 33 082103
Google Scholar
[20] Daly B J 1969 Phys. Fluids 12 1340
Google Scholar
[21] Zhang R Y, He X Y, Chen S Y 2000 Comput. Phys. Commun. 129 121
Google Scholar
[22] Young Y N, Ham F E 2006 J. Turbul. 7 71
Google Scholar
[23] Matsuoka C 2009 Phys. Fluids 21 092107
Google Scholar
[24] Sohn S I 2009 Physical Review E 80 055302(R
[25] 夏同军, 董永强, 曹义刚 2013 62 214702
Google Scholar
Xia T J, Dong Y Q, Cao Y G 2013 Acta. Phys. Sin. 62 214702
Google Scholar
[26] 黄皓伟, 梁宏, 徐江荣 2021 70 114701
Google Scholar
Huang H W, Liang H, Xu J R 2021 Acta Phys. Sin. 70 114701
Google Scholar
[27] Guo Z L, Shu C 2013 Lattice Boltzmann Method and its Applications in Engineering (Singapore: World Scientific Publishing Company) pp239–284
[28] Wang H L, Yuan X L, Liang H, Chai Z H, Shi B C 2019 Capillarity 2 33
Google Scholar
[29] Liang H, Hu X L, Huang X F, Xu J R 2019 Phys. Fluids 31 112104
Google Scholar
[30] Liang H, Shi B C, Guo Z L, Chai Z H 2014 Phys. Rev. E 89 053320
Google Scholar
[31] Liang H, Shi B C, Chai Z H 2017 Comput. Math. Appl. 73 1524
Google Scholar
[32] d’Humieres D, Ginzburg I, Krafczyk M, Lallemand P, Luo L S 2002 Philos. Trans. R. Soc. London, Ser. A 360 437
Google Scholar
[33] He X Y, Chen S Y, Zhang R Y 1999 J. Comput. Phys. 152 642
Google Scholar
[34] Abarzhi S I, Gorobets A, Sreenivasan K R 2005 Phys. Fluids 17 081705
Google Scholar
[35] Olson D H, Jacobs J W 2009 Phys. Fluids 21 034103
Google Scholar
DownLoad:
Catalog
Metrics
- Abstract views: 3772
- PDF Downloads: 71
- Cited By: 0
