FTM Numerical Simulation of Point Heat Source in Marangoni Without Gravity

doi:10.19596/j.cnki.1001-246x.8395

Abstract

Abstract:

Front-tracking method (FTM) is used to study droplet motion caused by Marangoni effect with a point heat source in the diagonal direction. Droplet motion at different Marangoni (Ma) numbers was simulated. It was found that the droplet velocity which decreases as Ma number increasing, increases rapidly to a stable migration velocity then decreases. However, it reverses at t=1.2, which increases as Ma number increasing. Whereupon, the vortex of Hill sphere split and the vortex center moved slightly. The same reflux exists in the droplet as in the Hill sphere vortex. With the increasing of Ma, end topology of temperature field bifurcates twice: The lowest temperature point jumps into the droplet from the stagnation point when occurs firstly in the lower critical Ma number; The inner shell cooling zone breaks from the central point when occurs secondly at upper, then forming toroidal cooling zone.

Key words: Marangoni, point heat source, FTM, numerical modeling

Mengying SUN, Ming MA, Hailong GUO, Mengjun YAO, Meng XU, Ying ZHANG. FTM Numerical Simulation of Point Heat Source in Marangoni Without Gravity[J]. Chinese Journal of Computational Physics, 2022, 39(2): 191-200.

Figures/Tables 14

Fig.1 Lagrangian grid and Eulerian grid

Fig.2 Weight coefficient map (indicated by area fractions)

Fig.3 Physical model of the simulation

Fig.4 Droplet migration velocity and theoretical prediction

Fig.5 Grid independence test

Fig.6 Aspect ratio as a function of Ma

Fig.7 Droplet trajectories

Fig.8 Droplet migration velocity as functions of time at different Ma

Fig.9 Dimensionless migration velocity of droplets as a function of Ma

Fig.10 Influence of Ma number on (A) streamline and (B) pressure distributions in and around a droplet

Fig.11 Droplet surface temperature difference as functions of Ma

Fig.12 Position of vortex center relative to droplet center

Fig.13 Isotherms at different Ma

Fig.14 Movement of the lowest temperature in the drop difference at different Ma numbers

References 25

1	刘会敏. 喷墨打印电致发光薄膜及器件制备方法的研究[D]. 广州: 华南理工大学, 2016.
2	田野, 宋义虎, 吴子良, 等. 基于马兰戈尼效应制备物理水凝胶薄膜[C]//中国化学会高分子学科委员会. 中国化学会2017全国高分子学术论文报告会摘要集——主题C: 高分子物理与软物质, 2017: 1.
3	LE C , LIU L , LI Z . Numerical investigation of the effect of rotation on the oscillatory thermocapillary convection and dopant transport in a silicon liquid bridge[J]. Journal of Crystal Growth, 2019, 523, 125149. DOI
4	谢雄. AZ31镁合金NSA-TIG焊接溶池运动规律及接头性能强化机理研究[D]. 重庆: 重庆大学, 2016.
5	方肇伦, 方群. 微流控芯片发展与展望[J]. 现代科学仪器, 2001, 18 (4): 3- 6. DOI
6	汪洋, 陈杰, 王智慧, 等. 多相传质过程中的Marangoni效应[J]. 化工学报, 2013, 64 (1): 124- 132. DOI
7	YOUNG N O , GOLDSTEIN J S , BLOCK M J . The motion of bubbles in a vertical temperature gradient[J]. Journal of Fluid Mechanics, 1959, 6 (3): 350- 356. DOI
8	BRETUKHIN Y K . Thermocapillary drift of a droplet of viscous liquid[J]. Fluid Dynamics, 1975, 10 (5): 833- 837.
9	BALASUBRAMANIAM R , CHAI A T . Thermocapillary migration of droplets: An exact solution for small Marangoni numbers[J]. Journal of Colloid and Interface Science, 1987, 119 (2): 531- 538. DOI
10	XIE J C , LIN H , HAN J H , et al. Experimental investigation of thermocapillary migration of isolated drops[J]. Advances in Space Research, 1999, 24 (10): 1409- 1415. DOI
11	XIE J C , LIN H , HAN J H , et al. Experimental investigation on Marangoni drop migrations using drop shaft facility[J]. International Journal of Heat and Mass Transfer, 1998, 41 (14): 2077- 2081. DOI
12	MEYYAPPAN M , SUBRAMANIAN R S , WILCOX W R , et al. Bubble behavior in molten glass in a temperature gradient[J]. Mrs Proceedings, 1981, 9, 311.
13	TREUNER M , GALINDO V , GERBETH G , et al. Thermocapillary bubble migration at high Reynolds and Marangoni numbers under low gravity[J]. Journal of Colloid and Interface Science, 1996, 179 (1): 114- 127. DOI
14	HADLAND P H , BALASUBRAMANIAM R , WOZNIAK G , et al. Thermocapillary migration of bubbles and drops at moderate to large Marangoni number and moderate Reynolds number in reduced gravity[J]. Experiments in Fluids, 1999, 26 (3): 240- 248. DOI
15	ZHAO J F , ZHANG L , LI Z D , et al. Topological structure evolvement of flow and temperature fields in deformable drop Marangoni migration in microgravity[J]. International Journal of Heat and Mass Transfer, 2011, 54 (21-22): 4655- 4663. DOI
16	LU M , DONG B H , ZHANG Y , et al. Numerical simulation of free-rising of bubbles in a tube with fins using front tracking method[J]. Chinese J Comput Phys, 2018, 35 (1): 47- 54.
17	YAO M J , SHANG W Q , ZHANG Y , et al. Numerical analysis of Kelvin-Helmholtz instability in inclined walls[J]. Chinese J Comput Phys, 2019, 36 (4): 403- 412.
18	ZHANG Y , YAO M , SHANG W , et al. Numerical analysis of coupled Kelvin-Helmholtz and Rayleigh-Taylor instability on inclined walls[J]. Canadian Journal of Physics, 2020, 98 (8): 790- 800. DOI
19	GAO H , GAO R F , YAO M J , et al. Numerical simulation of fluid-thermal coupled in droplet liquid film: Front tracking method[J]. Chinese J Comput Phys, 2020, 37 (4): 422- 430.
20	MA M , LU J , TRYGGVASON G . Using statistical learning to close two-fluid multiphase flow equations for a simple bubbly system[J]. Physics of Fluids, 2015, 27 (9): 092101. DOI
21	雷杰, 王昱, 马明, 等. 基于FTM方法的双气泡融合特性模拟[J]. 过程工程学报, 2019, 19 (2): 263- 270.
22	TRYGGVASON G , SCARDOVELLI R , ZALESKI S . Direct numerical simulations of gas-liquid multiphase flows[M]. Cambridge: Cambridge University Press, 2011.
23	张莹, 许术方, 李文彬, 等. 液滴撞击亲疏水间隔条纹表面[J]. 化工进展, 2020, 39 (2): 461- 467.
24	LI P , PENG C , DU P , et al. The investigation of the viscous fingering phenomenon of immiscible fluids displacement by the lattice Boltzmann method[J]. Canadian Journal of Physics, 2020, 98 (7): 650- 659.
25	潘雁妮, 吴忱韩, 张莹. 基于界面追踪法的不同初始形状同轴双气泡上升直接数值模拟[J]. 科学技术与工程, 2018, 18 (4): 316- 320. DOI

[1]	Han XIAO, Chunxi LI, Haozhe SU, Xuemin YE. Numerical Simulation on Instability of Vertical Liquid Drainage [J]. Chinese Journal of Computational Physics, 2021, 38(6): 661-671.
[2]	YAO Mengjun, SHANG Wenqiang, ZHANG Ying, GAO Hui, ZHANG Daoxu, LIU Peiyao. Numerical Analysis of Kelvin-Helmholtz Instability in Inclined Walls [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2019, 36(4): 403-412.
[3]	LU Min, DONG Boheng, ZHANG Ying, ZENG Liang, LIN Xiangquan, DU Peng. Numerical Simulation of Free-Rising of Bubbles in a Tube with Fins Using Front Tracking Method [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2018, 35(1): 47-54.
[4]	YE Xuemin, JIANG Kai, LI Chunxi. Marangoni Effect on Soluble Surfactant Droplet Spreading in Different Spreading Systems [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2013, 30(5): 683-691.
[5]	PENG Weicai, HE Zeng, LI Peng. Coupled Finite Element/Wave Based Method for Acoustic Analysis [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2009, 26(3): 396-402.
[6]	SHEN Lin-fang, CHEN Deng-peng, QIAN Jing-ren. THE NUMERICAL MODELING OF ERBIUM-DOPED SUPERFLUORESCENT FIBER SOURCE [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2001, 18(1): 1-6.