-
Based on the multi-component phase field theory, in this paper we propose an axisymmetric lattice Boltzmann model for three-phase fluids. The proposed model takes advantage of two particle distribution functions for capturing phase interface among three different fluids, and another particle distribution function for solving the hydrodynamic equations for flow field. In order to describe the axisymmetric effect arising from the coordinate transformation, we elaborately design the equilibrium distribution function and forcing distribution function in the evolution equation, which ensures that the model can accurately recover the macroscopic governing equation for three-phase fluids. Also, the introduced source terms accounting for the axisymmetric effect contain no additional gradient term, which makes it be simpler than the existing lattice Boltzmann model for axisymmetric three-phase fluids. To validate the proposed model, a series of axisymmetric multiphase benchmark examples are performed, including the static double droplets, the spreading of liquid lens, and the binary-fluid Rayleigh-Plateau instability. It is reported that the present model can accurately capture the phase interface, and the predicted steady shapes of the liquid lens agree well with the analytical profiles. Then, the proposed model is used to study the three-phase Rayleigh-Plateau instability and the effects of the wavenumber and the radius ratio of liquid column on the interfacial dynamic behaviour, the breakup time of liquid threads and the size of daughter droplet are investigated in detail. It can be found that the compound liquid thread at a high wavenumber could break up into one main droplet and one satellite droplet, but the multiple satellite droplets can be produced at a low wavenumber, which leads to that the sizes of main and satellite droplets increase with the wavenumber at first and then decrease with it. Besides, we can observe that the inner fluid undergoes the breakup at earlier time than the middle fluid, and the breakup time for both inner and middle fluids increases with the decrease of the wavenumber. Finally, we can find that increasing the radius ratio of liquid column accelerates the breakup of inner-fluid thread, but prevents the breakup of the middle-fluid thread. In addition, the size of the compound main droplet increases with the radius ratio of liquid column, while the size of the compound satellite droplet doest not change much with it.
[1] Smith K A, Solis F J, Chopp D L 2002 Interfaces Free Bound. 4 263
[2] Bonhomme R, Magnaudet J, Duval F, Piar B 2012 J. Fluid Mech. 707 405
[3] Kim J 2007 Comput. Methods Appl. Mech. Eng. 196 4779
Google Scholar
[4] Guo Z L, Shu C 2013 Lattice Boltzmann Method and Its Applications in Engineering (Singapore: World Scientific) pp10–32
[5] 胡晓亮, 梁宏, 王会利 2020 物理学报 69 044701
Google Scholar
Hu X L, Liang H, Wang H L 2020 Acta Phys. Sin. 69 044701
Google Scholar
[6] 贺传晖, 刘高洁, 娄钦 2021 物理学报 70 244701
Google Scholar
He C H, Liu G J, Lou Q 2021 Acta Phys. Sin. 70 244701
Google Scholar
[7] Spencer T J, Halliday I, Care M C 2010 Phys. Rev. E 82 066701
Google Scholar
[8] Leclaire S, Reggio M, Trepanier J 2013 J. Comput. Phys. 246 318
Google Scholar
[9] Yu Y, Liu H H, Liang D, Zhang Y H 2019 Phys. Fluids 31 012108
Google Scholar
[10] Semprebon C, Kruger T, Kusumaatmaja H 2016 Phys. Rev. E 93 033305
Google Scholar
[11] Wöhrwag M, Semprebon C, Moqaddam M A, Karlin I 2018 Phys. Rev. Lett. 120 234501
Google Scholar
[12] Liang H, Shi B C, Chai Z H 2016 Phys. Rev. E 93 013308
Google Scholar
[13] Liang H, Xu J R, Chen J X, Chai Z H, Shi B C 2019 Appl. Math. Model. 73 487
Google Scholar
[14] Liang H, Chai Z H, Shi B C, Guo Z L, Zhang T 2014 Phys. Rev. E 90 063311
Google Scholar
[15] Haghani-Hassan-Abadi R, Rahimian M 2020 Acta Mech. 231 2323
Google Scholar
[16] Yang J X, Li Y B, Lee C Y, Kim J 2021 Eur. J. Mech.-B/Fluid 89 203
Google Scholar
[17] Boyer F, Lapuerta C 2006 ESAIM: Math. Model. Numer. Anal. 40 653
Google Scholar
[18] Boyer F, Lapuerta C, Minjeaud S, Piar B, Quintard M 2010 Transp. Porous Media 82 463
Google Scholar
[19] Jacqmin D 1999 J. Comput. Phys. 155 96
Google Scholar
[20] Qian Y, d’Humieres D, Lallemand P 1992 Europhys. Lett. 17 479
Google Scholar
[21] Liang H, Shi B C, Guo Z L, Chai Z H 2014 Phys. Rev. E 89 053320
Google Scholar
[22] Liang H, Li Y, Chen J X, Xu J R 2019 Int. J. Heat Mass Transfer 130 1189
Google Scholar
[23] Rawlinson J S, Widom B 1982 The Molecular Theory of Capillarity (Oxford: Clarendon Press) pp208–214
[24] Jiang F, Tsuji T 2017 Water Resour. Res. 53 11
Google Scholar
[25] Zheng L, Zheng S, Zhai Q L 2021 Comput. Fluid 218 104857
Google Scholar
[26] Tjahjadi M, Stone H A, Ottino J M 1992 J. Fluid Mech. 243 297
Google Scholar
-
图 1 静态双液滴测试中三个序参数的稳态分布 (a)
$ \phi_1 $ 的分布; (b)$ \phi_2 $ 的分布; (c)$ \phi_3 $ 的分布Figure 1. Steady-state distributions of three order parameters in the static double-droplet test: (a) Distribution of
$ \phi_1 $ ; (b) distribution of$ \phi_2 $ ; (c) distribution of$ \phi_3 $ 
图 2 静态双液滴测试中序参数沿水平线的稳态分布
Figure 2. Steady distributions of the order parameters across the horizontal line in static double-droplet test
图 3 液体透镜相关几何参数示意图
Figure 3. Schematic diagram of relevant geometric parameters of liquid lens
图 4 不同界面张力比液体透镜的平衡形态 (a)
$\sigma_{12}:\sigma_{13}: $ $ \sigma_{23}=0.6:0.6:1$ ; (b)$ \sigma_{12}:\sigma_{13}:\sigma_{23}=1:1:1 $ ; (c)$\sigma_{12}: $ $ \sigma_{13}:\sigma_{23}=1:4/3:1$ .Figure 4. Equilibrium morphology of liquid lens at different interfacial tension ratios: (a)
$ \sigma_{12}:\sigma_{13}:\sigma_{23}=0.6:0.6:1 $ ; (b)$ \sigma_{12}:\sigma_{13}:\sigma_{23}=1:1:1 $ ; (c)$\sigma_{12}:\sigma_{13}:\sigma_{23}= 1:4/ $ $ 3:1$ .
图 5 不同波数下两相Rayleigh-Plateau不稳定性的界面演化过程 (a)
$ k=0.63 $ ; (b)$ k=0.31 $ ; (c)$ k=0.21 $ .Figure 5. Interfacial evolution processes of two-phase Rayleigh-Plateau instability with different wavenumbers: (a)
$ k=0.63 $ ; (b)$ k=0.31 $ ; (c)$ k=0.21 $ 
图 6 不同波数下三相Rayleigh-Plateau不稳定性的界面演化过程 (a)
$ k=0.42 $ ; (b)$ k=0.21 $ ; (c)$ k=0.14 $ Figure 6. Interfacial evolution processes of three-phase Rayleigh-Plateau instability with different wavenumbers: (a)
$ k=0.42 $ ; (b)$ k=0.21 $ ; (c)$ k=0.14 $ 
图 7 不同波数下复合液体线程破裂生成主液滴和卫星液滴的尺寸
Figure 7. The sizes of main and satellite droplets originating from the breakup of compound liquid threads with different wavenumbers
图 8 不同内外液柱半径比下三相Rayleigh-Plateau不稳定性的界面演化过程 (a)
$ RM/RI=1.5 $ ; (b)$RM/RI= $ $ 2$ ; (c)$ RM/RI=2.5 $ Figure 8. Interfacial evolution of three-phase Rayleigh-Plateau instability with different radius ratios of the inner and outer liquid columns: (a)
$ RM/RI=1.5 $ ; (b)$ RM/RI=2 $ ; (c)$ RM/RI=2.5 $ 
图 9 不同液柱半径比下复合液体线程破裂生成复合主液滴和卫星液滴的尺寸
Figure 9. The sizes of the compound main and satellite droplets originating from the breakup of compound liquid threads with different radius ratios
表 1 不同表面张力比的液体透镜的平衡接触角
$\theta_1$ 和$\theta_2$ Table 1. Equilibrium contact angles
$\theta_1$ and$\theta_2$ of liquid lens at different surface tension ratios表面张力比
$(\sigma_{12}:\sigma_{13}:\sigma_{23})$数值解 解析解 相对误差 $\theta_1$ $\theta_2$ $\theta_1$ $\theta_2$ $\theta_1$ $\theta_2$ $0.6:0.6:1$ 34.0 34.0 33.6 33.6 1.2% 1.2% $1:1:1$ 60.1 60.1 60 60 0.2% 0.2% $1:4/3:1$ 84.4 49.0 83.6 48.2 1.0% 1.7% 
DownLoad: CSV
表 2 不同表面张力比液体透镜的长度d和高度
$h_1$ ,$h_2$ Table 2. Length d and heights
$h_1$ ,$h_2$ of liquid lens at different surface tensions表面张力比$(\sigma_{12}:\sigma_{13}:\sigma_{23})$ 数值解 解析解 相对误差 d $h_1$ $h_2$ d $h_1$ $h_2$ d $h_1$ $h_2$ $0.6:0.6:1$ 322.1 49.2 49.2 391.8 59.1 59.1 17.8% 16.8% 16.8% $1:1:1$ 253.0 73.2 73.2 277.0 80.0 80.0 8.7% 8.5% 8.5% $1:4/3:1$ 234.5 106.1 53.4 251.6 112.5 56.3 6.8% 5.7% 5.2%
DownLoad: CSV
-
[1] Smith K A, Solis F J, Chopp D L 2002 Interfaces Free Bound. 4 263
[2] Bonhomme R, Magnaudet J, Duval F, Piar B 2012 J. Fluid Mech. 707 405
[3] Kim J 2007 Comput. Methods Appl. Mech. Eng. 196 4779
Google Scholar
[4] Guo Z L, Shu C 2013 Lattice Boltzmann Method and Its Applications in Engineering (Singapore: World Scientific) pp10–32
[5] 胡晓亮, 梁宏, 王会利 2020 物理学报 69 044701
Google Scholar
Hu X L, Liang H, Wang H L 2020 Acta Phys. Sin. 69 044701
Google Scholar
[6] 贺传晖, 刘高洁, 娄钦 2021 物理学报 70 244701
Google Scholar
He C H, Liu G J, Lou Q 2021 Acta Phys. Sin. 70 244701
Google Scholar
[7] Spencer T J, Halliday I, Care M C 2010 Phys. Rev. E 82 066701
Google Scholar
[8] Leclaire S, Reggio M, Trepanier J 2013 J. Comput. Phys. 246 318
Google Scholar
[9] Yu Y, Liu H H, Liang D, Zhang Y H 2019 Phys. Fluids 31 012108
Google Scholar
[10] Semprebon C, Kruger T, Kusumaatmaja H 2016 Phys. Rev. E 93 033305
Google Scholar
[11] Wöhrwag M, Semprebon C, Moqaddam M A, Karlin I 2018 Phys. Rev. Lett. 120 234501
Google Scholar
[12] Liang H, Shi B C, Chai Z H 2016 Phys. Rev. E 93 013308
Google Scholar
[13] Liang H, Xu J R, Chen J X, Chai Z H, Shi B C 2019 Appl. Math. Model. 73 487
Google Scholar
[14] Liang H, Chai Z H, Shi B C, Guo Z L, Zhang T 2014 Phys. Rev. E 90 063311
Google Scholar
[15] Haghani-Hassan-Abadi R, Rahimian M 2020 Acta Mech. 231 2323
Google Scholar
[16] Yang J X, Li Y B, Lee C Y, Kim J 2021 Eur. J. Mech.-B/Fluid 89 203
Google Scholar
[17] Boyer F, Lapuerta C 2006 ESAIM: Math. Model. Numer. Anal. 40 653
Google Scholar
[18] Boyer F, Lapuerta C, Minjeaud S, Piar B, Quintard M 2010 Transp. Porous Media 82 463
Google Scholar
[19] Jacqmin D 1999 J. Comput. Phys. 155 96
Google Scholar
[20] Qian Y, d’Humieres D, Lallemand P 1992 Europhys. Lett. 17 479
Google Scholar
[21] Liang H, Shi B C, Guo Z L, Chai Z H 2014 Phys. Rev. E 89 053320
Google Scholar
[22] Liang H, Li Y, Chen J X, Xu J R 2019 Int. J. Heat Mass Transfer 130 1189
Google Scholar
[23] Rawlinson J S, Widom B 1982 The Molecular Theory of Capillarity (Oxford: Clarendon Press) pp208–214
[24] Jiang F, Tsuji T 2017 Water Resour. Res. 53 11
Google Scholar
[25] Zheng L, Zheng S, Zhai Q L 2021 Comput. Fluid 218 104857
Google Scholar
[26] Tjahjadi M, Stone H A, Ottino J M 1992 J. Fluid Mech. 243 297
Google Scholar
DownLoad:
Catalog
Metrics
- Abstract views: 7293
- PDF Downloads: 86
- Cited By: 0
