搜索

x

留言板

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

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

接触角迟滞和气-液界面张力温度敏感性对液滴蒸发动态特性的影响

李春曦 程冉 叶学民

引用本文:
Citation:

接触角迟滞和气-液界面张力温度敏感性对液滴蒸发动态特性的影响

李春曦, 程冉, 叶学民

Effect of contact angle hysteresis and sensitivity of gas-liquid interfacial tension to temperature of a sessile-drop on evaporation dynamics

Li Chun-Xi, Cheng Ran, Ye Xue-Min
PDF
HTML
导出引用
  • 基于润滑理论和滑移边界条件, 构建了考虑接触角迟滞时液滴在均匀加热固体壁面上蒸发的数学模型, 探讨了接触角迟滞对蒸发液滴动力学的影响, 并分析了不同气-液界面张力温度敏感性的液滴蒸发特征. 研究结果表明: 接触角迟滞对液滴蒸发过程有较大影响, 随接触角迟滞增大, 液滴的接触线钉扎时间延长, 铺展阶段和去钉扎阶段时长缩短, 液滴蒸发显著加快; 迟滞角的增大使前进接触角增加, 后退接触角减小, 且后退接触角减小的幅度大于前进接触角增加的幅度. 提高气-液界面张力对温度的敏感性系数可通过减小后退接触角, 改善液滴在壁面上的润湿性, 从而加强液滴传热, 致使加快液滴蒸干. 因此, 改变接触角迟滞和气-液界面张力对温度的敏感性均可实现对液滴运动的调控, 从而控制其蒸发进程.
    The evaporation process of drops on a solid surface is widely applied to daily life and industrial fields. Both contact angle hysteresis and the sensitivity of gas-liquid interfacial tension to temperature are important factors affecting the drop evaporation reflected in the contact line and contact angle. To investigate the internal mechanism, according to the lubrication theory and slip boundary conditions, we establish a mathematical model of the drop evaporation on a uniformly heated solid wall with considering the effect of contact angle hysteresis. This model is numerically solved by using a coordinate transformation method and Freefem++14.3, a highly efficient solver. The accuracy of the numerical calculation method is verified by comparing the numerical results with experimental results, and the grid independence is validated. The effect of contact angle hysteresis on the dynamics of evaporating drops is discussed, and the evaporation characteristics of drops with different tension sensitivities of the air-liquid interface to temperature are further investigated. The results show that the contact angle hysteresis has an apparent influence on the drop evaporation process which includes drop spreading stage, contact line pinning stage, and depinning stage. In the drop spreading stage, the increase in the hysteresis angle shortens the spreading time, and reduces the spreading velocity and radius, while in the contact line pinning stage, the pinning time is prolonged and the reduction of drop mass is significantly increased with hysteresis angle increasing. In the contact line depinning stage, the contact angle hysteresis reduces the contact angle, and a flatter shape emerges, thereby enhancing the ability to transfer heat and accelerating evaporation as well as shortening the depinning time. In addition, the large hysteresis angle leads to a large advancing contact angle and a small receding contact angle. The reduction in receding contact angle is more notable than the increment of advancing contact angle. The temperature sensitivity coefficient of the gas-liquid interfacial tension can be increased by reducing the receding contact angle, thereby improving the wettability of the drops on the wall enhancing the heat transfer and accelerating the evaporation. Regulating the contact angle hysteresis and the sensitivity of the interfacial tension to temperature can realize the manipulation of the drop movement, thus controlling the evaporation process.
      通信作者: 叶学民, yexuemin@163.com
    • 基金项目: 国家自然科学基金(批准号: 51876065)资助的课题
      Corresponding author: Ye Xue-Min, yexuemin@163.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 51876065)
    [1]

    Kavehpour P, Ovryn B, McKinley G H 2002 Colloids Surf., A 206 409Google Scholar

    [2]

    Lee K S, Ivanova N, Starov V M, Hilal N, Dutschk V 2008 Adv. Colloid Interface Sci. 144 54Google Scholar

    [3]

    Wijshoff H 2010 Phys. Rep. 491 77Google Scholar

    [4]

    Putnam S A, Briones A M, Byrd L W, Ervin J S, Hanchak M S, White A, Jones J G 2012 Int. J. Heat Mass Transfer 55 5793Google Scholar

    [5]

    Ye X M, Zhang X S, Li M L, Li C X 2018 Phys. Fluids 30 112103Google Scholar

    [6]

    Lopes M C, Bonaccurso E 2012 Soft Matter 8 7875Google Scholar

    [7]

    Kiper I, Fulcrand R, Pirat C, Simon G, Stutz B, Ramos S M M 2015 Colloids Surf., A 482 617Google Scholar

    [8]

    Eral H B, Mannetje D, Oh J M 2013 Colloid. Polym. Sci. 291 247Google Scholar

    [9]

    Neumann A W, Good R J 1972 J. Colloid Interface Sci. 38 341Google Scholar

    [10]

    Johnson R E, Dettre R H 1964 Contact Angle, Wettability, and Adhesion (Washington: American Chemical Society) pp112−135

    [11]

    Dettre R H, Johnson R E 1964 Contact Angle, Wettability, and Adhesion (Washington: American Chemical Society) pp136−144

    [12]

    Yu H Z, Soolaman D M, Rowe A W, Banks J T 2004 ChemPhysChem 5 1035Google Scholar

    [13]

    Li Y F, Sheng Y J, Tsao H K 2013 Langmuir 29 7802Google Scholar

    [14]

    Trybala A, Okoye A, Semenov S, Agogo H, Rubio R G, Ortega F, Starov V M 2013 J. Colloid Interface Sci. 403 49Google Scholar

    [15]

    Kulinich S A, Farzaneh M 2009 Appl. Surf. Sci. 255 4056Google Scholar

    [16]

    Lin T S, Zeng Y H, Tsay R Y, Lin S Y 2016 J. Taiwan Inst. Chem. Eng. 62 54Google Scholar

    [17]

    Kuznetsov G V, Feoktistov D V, Orlova E G, Batishcheva K A 2016 Colloid J. 78 335Google Scholar

    [18]

    Ajaev V S 2005 J. Fluid Mech. 528 279Google Scholar

    [19]

    Masoud H, Felske J D 2009 Phys. Rev. E 79 016301Google Scholar

    [20]

    Semenov S, Starov V M, Rubio R G, Agogo H, Velarde M G 2012 Math. Modell. Nat. Phenom. 7 82Google Scholar

    [21]

    Karapetsas G, Sahu K C, Matar O K 2013 Langmuir 29 8892Google Scholar

    [22]

    Ye X M, Zhang X S, Li M L, Li C X, Dong S 2019 Int. J. Heat Mass Transfer 128 1263Google Scholar

    [23]

    叶学民, 李永康, 李春曦 2016 物理学报 65 104704Google Scholar

    Ye X M, Li Y K, Li C X 2016 Acta Phys. Sin. 65 104704Google Scholar

    [24]

    Karapetsas G, Craster R V, Matar O K 2011 J. Fluid Mech. 670 5Google Scholar

    [25]

    焦云龙, 刘小君, 逄明华, 刘焜 2015 物理学报 65 016801Google Scholar

    Jiao Y L, Liu X J, Pang M H, Liu K 2015 Acta Phys. Sin. 65 016801Google Scholar

    [26]

    Nagy M, Škvarla J 2013 Acta Montan. Slovaca 18 125

    [27]

    王晓东, 彭晓峰, 陆建峰, 王补宣 2003 应用基础与工程科学学报 11 174Google Scholar

    Wang X D, Peng X F, Wang B X 2003 Chin. J. Chem. Eng. 11 174Google Scholar

    [28]

    Smith M K 1995 J. Fluid Mech. 294 209Google Scholar

    [29]

    Craster R V, Matar O 2000 J. Fluid Mech. 425 235Google Scholar

    [30]

    Ma C, Liu J, Xie S, Liu Y 2020 Chem. Eng. Sci. 214 115418Google Scholar

    [31]

    Vasu B, Dubey A, Bég O A, Gorla R S R 2020 Comput. Biol. Med. 126 104025Google Scholar

    [32]

    Jiménez Bolaños S, Vernescu B 2017 Phys. Fluids 29 057103Google Scholar

    [33]

    朱仙仙, 闵春华, 郭宇虹, 王坤, 解立垚 2021 热科学与技术 20 28

    Zhu X X, Min C H, Guo Y H, Wang K, Xie L Y 2021 J. Therm. Sci. Tech. 20 28

    [34]

    Hu H, Larson R G 2005 Langmuir 21 3963Google Scholar

    [35]

    Gatapova E Y, Semenov A A, Zaitsev D V, Kabov O A 2014 Colloids Surf., A 441 776Google Scholar

    [36]

    Chu F, Wu X, Zhu Y, Yuan Z 2017 Int. J. Heat Mass Transfer 111 836Google Scholar

    [37]

    Bormashenko E, Musin A, Zinigrad M 2011 Colloids Surf., A 385 235Google Scholar

    [38]

    Brutin D, Sobac B (Brutin D Ed.) 2015 Droplet Wetting and Evaporation (Oxford: Academic Press) pp25−30

  • 图 1  液滴在加热表面铺展示意图

    Fig. 1.  Schematic of a drop disposed on a horizontally heated substrate.

    图 2  计算模型的验证 (a) 接触角的演化过程; (b) 接触线的演化过程; (c) 网格无关性验证

    Fig. 2.  Verification of the computational model: (a) Evolution of contact angle; (b) evolution of contact line; (c) validation of grid independence.

    图 3  不同Δθ时液滴演化过程中的特征参数变化 (a) 接触角; (b) 接触半径; (c) 接触线移动速度; (d) 蒸发剩余质量

    Fig. 3.  Characteristic parameters of drop evolution with time at different Δθ: (a) Dynamic contact angle; (b) contact line; (c) contact line velocity; (d) remaining mass.

    图 5  动态接触角与接触线移动速度关系图

    Fig. 5.  Relation between the dynamic contact angle and velocity of contact line.

    图 4  液滴演化过程轮廓图 (a) t = 10000; (b) t = 50000; (c) t = 70000; (d) t = 74000

    Fig. 4.  Drop profile during the drop evolution process: (a) t = 10000; (b) t = 50000; (c) t = 70000; (d) t = 74000.

    图 6  液滴演化过程中的特征参数变化 (a) 接触角; (b) 接触半径; (c) 接触线移动速度; (d) 蒸发剩余质量

    Fig. 6.  Characteristic parameters of drop evolution with time: (a) Dynamic contact angle; (b) contact line; (c) contact line velocity; (d) remaining mass.

    表 1  不同Δθ对液滴演化时间的影响

    Table 1.  Effect of varied Δθ on the drop evolution time.

    Δθ钉扎时刻去钉扎时刻蒸干时刻铺展阶段钉扎阶段去钉扎阶段
    0.271891039750780000—1891018910—3975039750—78000
    0.351653043250762000—1653016530—4325043250—76200
    0.401462045100746000—1462014620—4510045100—74600
    下载: 导出CSV

    表 2  不同迟滞角下θaθr 变化

    Table 2.  Variation in θa and θr of drop with different Δθs

    ΔθθaθrΔθaΔθr
    0.270.9660.6960.106–0.164
    0.350.9970.6470.031–0.049
    0.401.0160.6160.019–0.031
    下载: 导出CSV

    表 3  不同Ωlg时液滴演化过程中参数变化

    Table 3.  Varieties of parameters during drop evolution with different Ωlg.

    Ωlg钉扎时刻去钉扎时刻蒸干时刻θeθaθr
    0.0031891039750780000.8600.9660.686
    0.0052028040140767000.8370.9420.672
    0.012220040260710000.7750.8800.610
    下载: 导出CSV
  • [1]

    Kavehpour P, Ovryn B, McKinley G H 2002 Colloids Surf., A 206 409Google Scholar

    [2]

    Lee K S, Ivanova N, Starov V M, Hilal N, Dutschk V 2008 Adv. Colloid Interface Sci. 144 54Google Scholar

    [3]

    Wijshoff H 2010 Phys. Rep. 491 77Google Scholar

    [4]

    Putnam S A, Briones A M, Byrd L W, Ervin J S, Hanchak M S, White A, Jones J G 2012 Int. J. Heat Mass Transfer 55 5793Google Scholar

    [5]

    Ye X M, Zhang X S, Li M L, Li C X 2018 Phys. Fluids 30 112103Google Scholar

    [6]

    Lopes M C, Bonaccurso E 2012 Soft Matter 8 7875Google Scholar

    [7]

    Kiper I, Fulcrand R, Pirat C, Simon G, Stutz B, Ramos S M M 2015 Colloids Surf., A 482 617Google Scholar

    [8]

    Eral H B, Mannetje D, Oh J M 2013 Colloid. Polym. Sci. 291 247Google Scholar

    [9]

    Neumann A W, Good R J 1972 J. Colloid Interface Sci. 38 341Google Scholar

    [10]

    Johnson R E, Dettre R H 1964 Contact Angle, Wettability, and Adhesion (Washington: American Chemical Society) pp112−135

    [11]

    Dettre R H, Johnson R E 1964 Contact Angle, Wettability, and Adhesion (Washington: American Chemical Society) pp136−144

    [12]

    Yu H Z, Soolaman D M, Rowe A W, Banks J T 2004 ChemPhysChem 5 1035Google Scholar

    [13]

    Li Y F, Sheng Y J, Tsao H K 2013 Langmuir 29 7802Google Scholar

    [14]

    Trybala A, Okoye A, Semenov S, Agogo H, Rubio R G, Ortega F, Starov V M 2013 J. Colloid Interface Sci. 403 49Google Scholar

    [15]

    Kulinich S A, Farzaneh M 2009 Appl. Surf. Sci. 255 4056Google Scholar

    [16]

    Lin T S, Zeng Y H, Tsay R Y, Lin S Y 2016 J. Taiwan Inst. Chem. Eng. 62 54Google Scholar

    [17]

    Kuznetsov G V, Feoktistov D V, Orlova E G, Batishcheva K A 2016 Colloid J. 78 335Google Scholar

    [18]

    Ajaev V S 2005 J. Fluid Mech. 528 279Google Scholar

    [19]

    Masoud H, Felske J D 2009 Phys. Rev. E 79 016301Google Scholar

    [20]

    Semenov S, Starov V M, Rubio R G, Agogo H, Velarde M G 2012 Math. Modell. Nat. Phenom. 7 82Google Scholar

    [21]

    Karapetsas G, Sahu K C, Matar O K 2013 Langmuir 29 8892Google Scholar

    [22]

    Ye X M, Zhang X S, Li M L, Li C X, Dong S 2019 Int. J. Heat Mass Transfer 128 1263Google Scholar

    [23]

    叶学民, 李永康, 李春曦 2016 物理学报 65 104704Google Scholar

    Ye X M, Li Y K, Li C X 2016 Acta Phys. Sin. 65 104704Google Scholar

    [24]

    Karapetsas G, Craster R V, Matar O K 2011 J. Fluid Mech. 670 5Google Scholar

    [25]

    焦云龙, 刘小君, 逄明华, 刘焜 2015 物理学报 65 016801Google Scholar

    Jiao Y L, Liu X J, Pang M H, Liu K 2015 Acta Phys. Sin. 65 016801Google Scholar

    [26]

    Nagy M, Škvarla J 2013 Acta Montan. Slovaca 18 125

    [27]

    王晓东, 彭晓峰, 陆建峰, 王补宣 2003 应用基础与工程科学学报 11 174Google Scholar

    Wang X D, Peng X F, Wang B X 2003 Chin. J. Chem. Eng. 11 174Google Scholar

    [28]

    Smith M K 1995 J. Fluid Mech. 294 209Google Scholar

    [29]

    Craster R V, Matar O 2000 J. Fluid Mech. 425 235Google Scholar

    [30]

    Ma C, Liu J, Xie S, Liu Y 2020 Chem. Eng. Sci. 214 115418Google Scholar

    [31]

    Vasu B, Dubey A, Bég O A, Gorla R S R 2020 Comput. Biol. Med. 126 104025Google Scholar

    [32]

    Jiménez Bolaños S, Vernescu B 2017 Phys. Fluids 29 057103Google Scholar

    [33]

    朱仙仙, 闵春华, 郭宇虹, 王坤, 解立垚 2021 热科学与技术 20 28

    Zhu X X, Min C H, Guo Y H, Wang K, Xie L Y 2021 J. Therm. Sci. Tech. 20 28

    [34]

    Hu H, Larson R G 2005 Langmuir 21 3963Google Scholar

    [35]

    Gatapova E Y, Semenov A A, Zaitsev D V, Kabov O A 2014 Colloids Surf., A 441 776Google Scholar

    [36]

    Chu F, Wu X, Zhu Y, Yuan Z 2017 Int. J. Heat Mass Transfer 111 836Google Scholar

    [37]

    Bormashenko E, Musin A, Zinigrad M 2011 Colloids Surf., A 385 235Google Scholar

    [38]

    Brutin D, Sobac B (Brutin D Ed.) 2015 Droplet Wetting and Evaporation (Oxford: Academic Press) pp25−30

  • [1] 扶龙香, 贺少波, 王会海, 孙克辉. 离散忆阻混沌系统的Simulink建模及其动力学特性分析. 物理学报, 2022, 71(3): 030501. doi: 10.7498/aps.71.20211549
    [2] 刘春杰, 赵新军, 高志福, 蒋中英. 高分子混合刷吸附/脱附蛋白质的模型化研究. 物理学报, 2021, 70(22): 224701. doi: 10.7498/aps.70.20211219
    [3] 扶龙香, 贺少波, 王会海, 孙克辉. 离散忆阻混沌系统的Simulink建模及其动力学特性分析. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211549
    [4] 徐珂, 许龙, 周光平. 考虑水蒸气蒸发和冷凝的球状泡群中泡的动力学特性. 物理学报, 2021, 70(19): 194301. doi: 10.7498/aps.70.20210045
    [5] 宁利中, 张珂, 宁碧波, 刘爽, 田伟利. 倾斜Poiseuille-Rayleigh-Bénard流动的对流分区与动力学特性. 物理学报, 2020, 69(12): 124401. doi: 10.7498/aps.69.20191941
    [6] 肖利全, 段书凯, 王丽丹. 基于Julia分形的多涡卷忆阻混沌系统. 物理学报, 2018, 67(9): 090502. doi: 10.7498/aps.67.20172761
    [7] 叶学民, 张湘珊, 李明兰, 李春曦. 液滴在不同润湿性表面上蒸发时的动力学特性. 物理学报, 2018, 67(11): 114702. doi: 10.7498/aps.67.20180159
    [8] 包涵, 包伯成, 林毅, 王将, 武花干. 忆阻自激振荡系统的隐藏吸引子及其动力学特性. 物理学报, 2016, 65(18): 180501. doi: 10.7498/aps.65.180501
    [9] 吴赛, 李伟斌, 石峰, 蒋世春, 蓝鼎, 王育人. 受限胶体液滴蒸发过程中胶体颗粒沉积过程观察. 物理学报, 2015, 64(9): 096101. doi: 10.7498/aps.64.096101
    [10] 王小发, 李骏. 短外腔偏振旋转光反馈下1550 nm垂直腔面发射激光器的动力学特性研究. 物理学报, 2014, 63(1): 014203. doi: 10.7498/aps.63.014203
    [11] 葛宋, 陈民. 接触角与液固界面热阻关系的分子动力学模拟. 物理学报, 2013, 62(11): 110204. doi: 10.7498/aps.62.110204
    [12] 蒋树庆, 甯家敏, 陈法新, 叶繁, 薛飞彪, 李林波, 杨建伦, 陈进川, 周林, 秦义, 李正宏, 徐荣昆, 许泽平. Z箍缩动态黑腔动力学及辐射特性初步实验研究. 物理学报, 2013, 62(15): 155203. doi: 10.7498/aps.62.155203
    [13] 李洪伟, 周云龙, 刘旭, 孙斌. 基于随机子空间结合稳定图的气液两相流型分析. 物理学报, 2012, 61(3): 030508. doi: 10.7498/aps.61.030508
    [14] 黄丽莲, 辛方, 王霖郁. 新分数阶超混沌系统的研究与控制及其电路实现. 物理学报, 2011, 60(1): 010505. doi: 10.7498/aps.60.010505
    [15] 吴渝, 杨晶晶, 王利. Swarm模型突现行为的动力学特性分析. 物理学报, 2011, 60(10): 108902. doi: 10.7498/aps.60.108902
    [16] 和兴锁, 李雪华, 邓峰岩. 平面柔性梁的刚-柔耦合动力学特性分析与仿真. 物理学报, 2011, 60(2): 024502. doi: 10.7498/aps.60.024502
    [17] 包伯成, 周国华, 许建平, 刘中. 斜坡补偿电流模式控制开关变换器的动力学建模与分析. 物理学报, 2010, 59(6): 3769-3777. doi: 10.7498/aps.59.3769
    [18] 蒋飞, 刘中, 胡文, 包伯成. 连续混沌调频信号的动力学设计与分析. 物理学报, 2010, 59(1): 116-122. doi: 10.7498/aps.59.116
    [19] 郑桂波, 金宁德. 两相流流型多尺度熵及动力学特性分析. 物理学报, 2009, 58(7): 4485-4492. doi: 10.7498/aps.58.4485
    [20] 蔡国梁, 谭振梅, 周维怀, 涂文桃. 一个新的混沌系统的动力学分析及混沌控制. 物理学报, 2007, 56(11): 6230-6237. doi: 10.7498/aps.56.6230
计量
  • 文章访问数:  2011
  • PDF下载量:  49
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-02-08
  • 修回日期:  2021-05-16
  • 上网日期:  2021-09-30
  • 刊出日期:  2021-10-20

接触角迟滞和气-液界面张力温度敏感性对液滴蒸发动态特性的影响

  • 华北电力大学动力工程系, 保定 071003
  • 通信作者: 叶学民, yexuemin@163.com
    基金项目: 国家自然科学基金(批准号: 51876065)资助的课题

摘要: 基于润滑理论和滑移边界条件, 构建了考虑接触角迟滞时液滴在均匀加热固体壁面上蒸发的数学模型, 探讨了接触角迟滞对蒸发液滴动力学的影响, 并分析了不同气-液界面张力温度敏感性的液滴蒸发特征. 研究结果表明: 接触角迟滞对液滴蒸发过程有较大影响, 随接触角迟滞增大, 液滴的接触线钉扎时间延长, 铺展阶段和去钉扎阶段时长缩短, 液滴蒸发显著加快; 迟滞角的增大使前进接触角增加, 后退接触角减小, 且后退接触角减小的幅度大于前进接触角增加的幅度. 提高气-液界面张力对温度的敏感性系数可通过减小后退接触角, 改善液滴在壁面上的润湿性, 从而加强液滴传热, 致使加快液滴蒸干. 因此, 改变接触角迟滞和气-液界面张力对温度的敏感性均可实现对液滴运动的调控, 从而控制其蒸发进程.

English Abstract

    • 固体表面上液滴的蒸发过程不仅在日常生活中随处可见, 而且在涂覆工艺、焊接技术、油墨印刷、DNA分析等工业和医学领域也有广泛应用[1,2]. 掌控液滴在固体基底上的运动特性对于加快微电子元件冷却速率、提高喷涂质量及提升医疗诊断准确性等诸多实际应用至关重要[3].

      近年来的研究表明, 液滴蒸发过程的动态特性与固体表面温度[4]、液体性质[5]、固体材质[6]和环境条件[7]密切相关. 除此之外, 存在于实际表面上的接触角迟滞[8]也是影响液滴运动的重要因素. 理论上讲, 对于任意一个理想表面, 都可以由Young方程, $ \cos {\theta _Y}{\text{ = (}}{\sigma _{{\text{sg}}}} - {\sigma _{{\text{ls}}}})/{\sigma _{{\text{lg}}}} $, 唯一地确定其表面上所附着液滴的接触角. 但实际上很难通过实验来测出这一角度, 因为附着液滴在固体表面上的接触角值并不唯一, 而是在相对稳定的两个角度之间变化, 这种现象被称为接触角迟滞. 其中上限和下限分别为前进接触角θa和后退接触角θr, 二者的差值为迟滞角Δθ. 对于接触角迟滞, Neumann和Good[9]认为接触角迟滞现象归因于表面粗糙[10]、化学异质性[11]和亚稳表面能量态, 而Eral等[8]则发现表面变形、液体吸附和滞留、润湿时分子重排和相互扩散也是影响接触角迟滞的重要因素.

      目前, 已有实验研究了液滴蒸发过程中接触角迟滞的影响. Yu等[12]通过在自组装分子层上微液滴的蒸发实验, 发现蒸发模式的转变源于接触角迟滞的存在. Li等[13]发现接触线钉扎与接触角迟滞有关, 而接触线钉扎是咖啡环沉积形成的必要条件, 接触角迟滞直接影响液滴蒸发模式, 从而影响沉积模式, 液滴长时间附着在固体表面是由接触角迟滞引起的. Trybala等[14]通过实验探究了含无机纳米颗粒的悬浮液的固着液滴在三种不同润湿性基材(PE, PTFE和硅片)上的蒸发动力学行为, 指出当未达到静态后退接触角时, 三种表面上的液滴蒸发过程包括定接触线(CCR)和定接触角(CCA)两个阶段, 或仅呈现CCR阶段; 与纯水蒸发的普遍特征相比, 在所使用基底上观察到了良好一致性, 仅前进接触角和后退接触角数值有差异. Kulinich和Farzaneh[15]观察了具有相同接触角(≥ 150°)、不同钝化程度的超疏水聚合物表面上的蒸发模式, 发现在高迟滞表面上液滴蒸发遵循CCR模式, 而在低迟滞表面上则遵循CCA模式; 固着水滴的蒸发行为不是由高接触角决定, 而是受润湿滞后(或后退接触角)控制. Lin等[16]实验研究了6种不同表面粗糙度的聚甲基丙烯酸甲酯(PMMA)基底上水滴蒸发的动力学行为, 同样发现当壁面粗糙度增加时, 从CCR模式向CCA模式的转变延迟, 后退接触角也减小. Kuznetsov等[17]研究了不同粗糙度的铜表面上蒸馏水的蒸发特性, 发现液滴接触半径和钉扎持续时间随粗糙度增加而增大, 可通过改变基底粗糙度来控制液滴的铺展过程.

      关于液滴蒸发过程的理论研究, 通常针对接触线动力学的不同阶段建立不同的数学模型[18-20]进行描述. 其中, Semenov等[20]通过建立CCR和CCA模型描述接触角迟滞条件下固着液滴蒸发过程的两个阶段: 在CCR模型中, 假设液滴为球形, 接触半径始终等于其初始值L0, 仅研究接触角的演化; 在CCA模型中, 接触角被设定为常数并等于θr, 仅研究接触半径的演化, 即从L0减小到几乎为零. 虽然分开描述这两个阶段可以清晰得到各阶段的运动特征, 但却不能捕捉从CCR过渡到CCA时的重要动态特征, 因此不能全面刻画液滴蒸发过程的运动学特征. 此外, 还有学者利用润滑近似理论研究液滴在加热表面的动力学过程. Karapetsas等[21]通过将接触线速度和动态接触角与平衡接触角之间的差耦合成一个关系式来模拟接触线运动, 并引入不同界面张力对温度的敏感性, 研究了非等温倾斜固体基底上液滴的二维动力学, 证明了温度引起的平衡接触角变化会引起复杂的动力学过程. Ye等[22]进一步研究了均匀加热基底上固着液滴蒸发时的动力学特性, 分析了气-液、固-气和固-液界面张力温度敏感性对壁面润湿性对液滴动态特性的影响, 表明通过改变固-气界面张力温度敏感系数调控处于蒸发状态下的液滴运动更加有效. 此外, 叶学民等[23]研究了不同加热条件下温度影响壁面润湿性及铺展过程的内部机理, 发现随温度梯度增大, 液滴所受Marangoni效应增强, 从而加快液滴向低温区的铺展速率.

      上述实验研究表明, 接触角迟滞是影响液滴蒸发动力学特征的重要因素. 针对不同阶段建立的不同蒸发模型易忽略过渡阶段间的重要动态过程, 且模型不具有整体性, 而基于润滑理论建立的已有模型虽能描述整个液滴蒸发过程, 但考虑的均为理想表面, 未涉及接触角迟滞对实际表面上液滴润湿性的影响, 因此也不能完整刻画液滴的实际蒸发动力学特征. 此外, 在加热表面上的液滴界面张力受温度影响, 但至今尚未有能够准确测量固-液、固-气界面张力的可靠方法, 同时固-液、固-气界面张力温度敏感性系数也不易调控, 一般采用气-液界面张力的温度敏感性描述壁面润湿性对液滴运动特性的影响. 为此, 本文依据相关实验结果, 通过引入接触角迟滞模型, 基于润滑理论和滑移边界条件建立二维液滴厚度的演化方程, 采用数值计算方法模拟液滴蒸发全过程的动力学特征, 分析接触角迟滞对接触角、接触线位置和移动速度、蒸干时间等典型特征量的影响, 阐明接触角迟滞影响液滴蒸发动态过程的内在机理, 并分析考虑接触角迟滞时气-液界面张力的温度敏感性对液滴蒸发过程的影响.

    • 假设一不可压缩液滴在均匀加热的固体表面上不断蒸发, 如图1所示, 液滴初始最大厚度为$ \tilde H $, 流动方向尺度为$ \tilde L $, 壁面温度为$ \tilde T_{\text{w}}^{} $. 由于$\varepsilon = \tilde H/\tilde L $$ \ll 1$ (图1中的$ \tilde H $, $ \tilde L $仅为表示清楚, 不代表真实比例), 因此可应用润滑理论(上标~表示有量纲量, 下同).

      图  1  液滴在加热表面铺展示意图

      Figure 1.  Schematic of a drop disposed on a horizontally heated substrate.

      液滴运动过程的控制方程包括连续性方程、动量方程和能量方程:

      $ \tilde \nabla \cdot\tilde {\boldsymbol{u}} = 0, $

      $ {\tilde {\boldsymbol{u}}_{\tilde t}} + \tilde {\boldsymbol{u}}\cdot\tilde \nabla \tilde {\boldsymbol{u}} = \tilde g-\frac{1}{{\tilde \rho }}\tilde \nabla \tilde p + \frac{{\tilde \mu }}{{\tilde \rho }}{\tilde \nabla ^2}\tilde {\boldsymbol{u}}, $

      $ \tilde \rho {\tilde C_{\text{p}}}({\tilde T_{\tilde t}} + \tilde {\boldsymbol{u}}\cdot\tilde \nabla \tilde T) = \tilde \lambda {\tilde \nabla ^2}\tilde T, $

      式中, $ \tilde p $, $ \tilde T $$ \tilde t $分别表示压强、温度和时间, $ \tilde \mu $$ \tilde \rho $分别为液体的动力黏度与密度, $ \tilde \lambda $$ {\tilde C_{\text{p}}} $为液体的导热系数和定压比热; $ \tilde {\boldsymbol{u}} $= ($ \tilde u $, $ \tilde w $), $ \tilde u $$ \tilde w $分别为水平和垂直方向上的速度分量.

      在气-液界面$ \tilde z = \tilde h $处, 满足切向和法向应力平衡:

      $ \tilde{{\boldsymbol{n}}}\cdot \tilde{{\boldsymbol{\tau}} }\cdot \tilde{{\boldsymbol{t}}}=\tilde{{\boldsymbol{t}}}\cdot {\tilde{\nabla }}_{\text{s}}{\tilde{\sigma }}_{\text{lg}}, $

      $ \tilde {\boldsymbol{n}}\cdot\tilde {\boldsymbol{\tau}} \cdot\tilde {\boldsymbol{n}} = \tilde \kappa {\tilde \sigma _{{\text{lg}}}}, $

      式中, $ \tilde {\boldsymbol{n}} $, $ \tilde {\boldsymbol{t}} $, $ \tilde {\boldsymbol{\tau}} $分别为界面处的单位法向和切向向量及总应力张量, $ \tilde{\kappa } $为气-液界面处的平均曲率, $ {\tilde \sigma _{{\text{lg}}}} $为液-气界面张力.

      在气-液界面$ \tilde z = \tilde h $处, 运动边界及热边界条件为

      $ {\tilde h_{\tilde t}} + \tilde u{\tilde h_{\tilde x}} = \tilde w - \frac{{\tilde J}}{{\tilde \rho }}, $

      $ \tilde \lambda \tilde {\boldsymbol{n}} \cdot \tilde \nabla \tilde T + \tilde J{\tilde L_{\text{a}}} = 0, $

      式中, $ {\tilde L_{\text{a}}} $为汽化潜热, $ \tilde J $为界面质量迁移的质量跃变流率.

      界面处的蒸发流量与界面温度有关, 可表示为

      $ \tilde K\tilde J = {\tilde T_{\text{i}}} - {\tilde T_{\text{s}}}, $

      式中, $ \tilde{K} $为界面热阻; $ {\tilde T_{\text{i}}} $为界面温度, $ {\tilde T_{\text{s}}} $为液体的饱和温度.

      在壁面$ \tilde z = 0 $处, 垂直方向上满足无渗透条件:

      $ \tilde w = 0 . $

      为避免接触线处出现应力奇点, 在水平方向施加Navier滑移边界条件[24,25]:

      $ \tilde u = \tilde \beta {\tilde u_{\tilde z}}, $

      式中$ \tilde \beta $为滑移长度.

      壁面温度满足

      $ {\tilde T_{\text{w}}} = {\tilde T_{{\text{w0}}}}. $

      由于固-液、固-气界面张力目前仍难以通过实验直接测量, 假设固-液、固-气界面张力不随温度改变, 而气-液界面张力与温度满足线性关系[21]:

      $ \begin{split} &{{\tilde \sigma }_{{\text{sg}}}} = {{\tilde \sigma }_{{\text{s}}g,{{\tilde T}_{\text{s}}}}}, ~~ {{\tilde \sigma }_{{\text{ls}}}} = {{\tilde \sigma }_{{\text{ls,}}{{\tilde T}_{\text{s}}}}}, \hfill \\ &{{\tilde \sigma }_{{\text{lg}}}} = \tilde \sigma {}_{{\text{lg}},{{\tilde T}_{\text{s}}}} - {{\tilde \xi }_{1{{\tilde T}_{\lg }}}}({{\tilde T}_{{\text{lg}}}} - {{\tilde T}_{\text{s}}}), \hfill \end{split} $

      式中, ${\tilde \sigma _{\lg, {{\tilde T}_{\text{s}}}}}$为气-液界面在${\tilde T_{\text{s}}}$下的界面张力, ${\tilde \xi _{1{{\tilde T}_{\lg }}}}$为温度引起的表面张力梯度.

      将控制方程和边界条件采用下式进行无量纲化:

      $\begin{split} &(\tilde x,\tilde z,\tilde h) = \tilde L(x,\varepsilon z,\varepsilon h),\\ &\tilde t = \frac{{\tilde L}}{{\tilde U}}t, ~~ \left( {\tilde u,\tilde w} \right) = \tilde U\left( {u,\varepsilon w} \right), \\ &\tilde p = \frac{{\tilde S}}{{\tilde H}}p, ~~ \tilde T = T({\tilde T_m} - {\tilde T_{\text{s}}}) + {\tilde T_{\text{s}}}, \\ &{\varOmega _{\lg }} = \frac{{{{\tilde \xi }_{1{{\tilde T}_{\lg }}}}({{\tilde T}_{\text{m}}} - {{\tilde T}_{\text{s}}})}}{{{{\tilde \sigma }_{\lg,{{\tilde T}_{\text{s}}}}} - {{\tilde \sigma }_{\lg,{{\tilde T}_{\text{m}}}}}}}, \\ &{\tilde{\sigma }}_{\text{i}}=\tilde{S}{\sigma }_{\text{i}} +{\tilde{\sigma }}_{\text{lg},{\tilde{T}}_{\text{m}}}, \\ &\tilde{J}={J}_{0}J, ~~ \tilde{K}=\frac{K\epsilon \tilde{L}{\tilde{L}}_{\text{a}}^{}}{\tilde{\lambda }}, \end{split}$

      式中, ${J_0} = \dfrac{{\tilde \lambda ({{\tilde T}_{\text{m}}} - {{\tilde T}_{\text{s}}})}}{{\tilde H\tilde L_{\text{a}}^{}}}$, $\tilde U = \dfrac{{\tilde S\tilde H}}{{\tilde \mu \tilde L}}$为特征速度, $\tilde S = $$ {\tilde \sigma _{{\text{lg}}, {{\tilde T}_{\text{s}}}}} - {\tilde \sigma _{{\text{lg}}, {{\tilde T}_{\text{m}}}}}$为铺展系数, Ωlg为气-液界面张力对温度的敏感系数.

      通过无量纲变换和保留数量级 ≥ Ο(ε), 可得无量纲控制方程组为

      $ {u_x} + {w_z} = 0, $

      $ {p_x} = {u_{zz}}, $

      $ {p_z} = - \varepsilon Bo, $

      $ {T_{zz}} = 0, $

      式中, $Bo = \dfrac{{\tilde \rho \tilde g{{\tilde H}^2}}}{{\tilde \mu \tilde U}}$为邦德数.

      无量纲边界条件为

      $ z = 0, \;\;\; w = 0, $

      $ u = \beta {u_z}, $

      $ T = {\gamma _0} ; $

      $ z = h, \;\;\; {u_z} = {\sigma _{{\text{lg}}x}}, $

      $ p = - Ca{h_{xx}}, $

      $ {h_t} + u{h_x} = w - EJ, $

      $ {T_z} = - J, $

      $ KJ = {T_{}}, $

      式中, $\beta = \dfrac{{\tilde \beta }}{{\tilde H}}$, ${\gamma _0} = \dfrac{{{{\tilde T}_{\text{w}}} - {{\tilde T}_{\text{s}}}}}{{{{\tilde T}_{\text{m}}} - {{\tilde T}_{\text{s}}}}}$, $Ca = \dfrac{{{\varepsilon ^2}{{\tilde \sigma }_{{\text{lg}}, {{\tilde T}_{\text{m}}}}}}}{{\tilde S}}$$E = \dfrac{{\tilde \lambda ({{\tilde T}_{\text{m}}} - {{\tilde T}_{\text{s}}})}}{{\varepsilon \tilde H\tilde U\tilde \rho \tilde L_{\text{a}}^{}}}$分别表示滑移系数、加热系数、毛细数和蒸发数.

      由(19)式和(24)式可得$J = \dfrac{{{\gamma _0}}}{{K + h}}$.

      无量纲形式的界面张力与温度关系为

      $ \begin{split} &{\sigma _{{\text{sg}}}} = {\sigma _{{\text{sg}},{T_{\text{s}}}}}, \hfill \\ &{\sigma _{{l} {\text{s}}}} = {\sigma _{{l} {\text{s}},{T_{\text{s}}}}}, \hfill \\ &{\sigma _{\lg }} = {\sigma _{\lg,{T_{\text{s}}}}}(1 - {\varOmega _{\lg }}T). \hfill \end{split} $

      液滴表面温度为

      $ T = - \frac{{{\gamma _0}}}{{K + h}}h + {\gamma _0} . $

      采用积分法可得液滴厚度的演化方程为

      $\begin{split} {h_{\rm{t}}} =\;& \left[ - Ca{h_{xxx}}\left(\frac{{{h^3}}}{3} + \beta {h^2}\right) + \varepsilon Bo\left(\frac{{{h^3}}}{3} + \beta {h^2}\right){h_x} \right.\\ &\left.- \frac{{{\varOmega _{\lg }}h{\gamma _0}}}{{{{(K + h)}^2}}}{h_x}\left(\frac{{{h^2}}}{2} + \beta h\right) \right]_x - \frac{{E{\gamma _0}}}{{K + h}} . \\[-15pt] \end{split}$

      (27)式右侧依次代表毛细力、重力、热毛细力和蒸发在液滴运动过程中的作用.

      液滴在理想光滑表面上, 在三相接触线处, Young方程[26]

      $ {\tilde \sigma _{{\text{sg}}}}({x_{\text{c}}}) = {\tilde \sigma _{\lg }}({x_{\text{c}}})\cos {\tilde \theta _{\text{e}}} + {\tilde \sigma _{{\text{ls}}}}({x_{\text{c}}}), $

      式中, $ \tilde \theta _{\text{e}}^{} $为液滴在理想光滑表面上的有量纲的平衡接触角.

      由于在液膜厚度演化模型推导过程中, 采用润滑理论对液滴纵横比进行了放大, 因此, 液滴平衡接触角也需进行相应放大, 即cos$ \tilde \theta _{\text{e}} $→cos(εθe), 又有$ 1 - ({\varepsilon ^2}\theta _{\text{e}}^{{\kern 1 pt} {\text{2}}})/2 \approx \cos (\varepsilon {\theta _{\text{e}}}) $. 结合(25)式和(28)式可得:

      $ \theta _{_{\text{e}}}^{{\kern 1pt} {\kern 1pt} {\kern 1pt} 2} = - \frac{2}{{{\varepsilon ^2}}}\left( {\frac{{C + {\varOmega _{\lg }}{\gamma _0}}}{{1 - {\varOmega _{\lg }}{\gamma _0}}}} \right), $

      式中, $ C = {\sigma _{{\text{sg}}, {T_{\text{s}}}}} - {\sigma _{{\text{ls}}, {T_{\text{s}}}}} - 1 $.

      为描述接触角迟滞影响, 采用以下线性经验关联式描述前进接触角和后退接触角与迟滞角之间的关系[27]:

      $ {\theta _{\text{a}}} = {\theta _{\text{e}}} + A\Delta \theta, $

      $ {\theta _{\text{r}}} = {\theta _{\text{e}}} + (A - 1)\Delta \theta, $

      式中, A ≥ 0是与粗糙度无关, 与液体性质有关的常数; 迟滞角$ \Delta \theta $与表面粗糙度和化学异质性有关.

      接触线处的液滴高度为零, 即h(x = xcl = xcr, t) = 0, 其中xclxcr分别代表左侧和右侧接触线位置. 采用下式描述接触线移动速率与接触角间的关系[28](下标c代表接触线):

      $ {\upsilon _{\text{c}}} = \frac{{{\text{d}}{x_{\text{c}}}}}{{{\text{d}}t}} = \left\{ \begin{aligned} &B{(\theta - {\theta _{\text{a}}})^m},\;{\text{ }}\theta > {\theta _{\text{a}}}, \hfill \\ & - B{({\theta _{\text{r}}} - \theta )^m},\;\theta < {\theta _{\text{r}}}, \hfill \end{aligned} \right. $

      式中, Bm为常数; θ为动态接触角; υc > 0代表接触线向外铺展, υc < 0代表向内收缩.

      计算初始条件为

      $ h(x,t = 0) = (1 - {x^2})[F(1 - x) - F( - 1 - x)], $

      $ {x_{{\text{cl}}}}(t = 0) = - 1,~~ {x_{{\text{cr}}}}(t = 0) = 1, $

      式中, F(x) = 0.5[1 + tanh(20x)][29].

      基于Karapetsas等[21]提出的坐标变换法, 将瞬态物理域(x, t)映射到固定计算域(x', t'),

      $ x' = (x - x_{{\text{cl}}}^{})\frac{{x_{{\text{cr}}}^{}(t = 0) - x_{{\text{cl}}}^{}(t = 0)}}{{x_{{\text{cr}}}^{} - x_{{\text{cl}}}^{}}} . $

      液滴运动区域设定为–1 ≤ x' ≤ 1, 时间导数表示如下:

      $ \frac{\partial }{{\partial t}} = \frac{\partial }{{\partial t'}} - \frac{{{\text{d}}x'}}{{{\text{d}}t}}\frac{\partial }{{\partial x'}}, $

      (36)式用于替换控制方程(27)的相应项.

      数值求解采用Freefem++14.3, 该软件是求解偏微分方程的强大工具, 是一款免费的、开放源代码的有限元计算软件, 目前已经比较成熟地运用在各类数值模拟中, 对于液滴、液膜的流体数值计算[30-32]也很常见. 该软件采用Delaunay算法生成离散偏微分方程所需网格, 具有网格自适应和移动网格生成的功能; 采用稀疏矩阵存储格式, 内存需求少, 计算速度快, 是求解复杂区域问题的高效计算软件.

    • 参照朱仙仙等[33]用4 μL去离子水在加热温度为60°的玻璃表面上的蒸发实验数据, 换算得到无量纲参数, 各参数取值为: ε = 0.1, A = 0.39, B = 0.007, m = 3, C = –0.004, Bo = 0.007, Ca = 0.5, K = 10, E = 1 × 10-3, Ωlg = 0.003. 首先, 为了验证数值模拟结果的准确性, 将接触角和接触线随时间变化的实验结果无量纲化后与模拟值进行对比, 结果如图2(a)图2(b)所示. 可以看出, Δθ = 0.27时的模拟值与实验结果总体符合良好, 即采用考虑接触角迟滞的理论模型可更准确地反映液滴的运动过程. 值得注意的是, 朱仙仙等[33]所给实验结果并未呈现液滴铺展阶段, 可能是由此阶段变化较快、不易观测记录所致; 模拟结果和实验结果有偏差的原因可能是实验时基底加热时受热不均匀或者基底传热能力不一致, 也有可能是基于润滑理论建立模型时, 未考虑液滴蒸气扩散[34]所致.

      图  2  计算模型的验证 (a) 接触角的演化过程; (b) 接触线的演化过程; (c) 网格无关性验证

      Figure 2.  Verification of the computational model: (a) Evolution of contact angle; (b) evolution of contact line; (c) validation of grid independence.

      此外, 模拟时进行了网格无关性验证, 在x方向上对比了网格数为300, 500和800的结果, 发现增加网格数使计算时长增大, 而各特征参数计算结果数值差距不大, 以图2(c)的接触线演化过程为例, 网格数为500时可同时满足精确性与高效性的要求. 因此, 计算中时间步长设为10–6, 将[0, 2] × [0, 1]的计算域划分为 500 × 2 个均匀网格. 因方程(27)是一维的, 液滴厚度hy方向无关, 所以为节约计算时间, 在y方向仅划分2层网格.

    • 综合朱仙仙等[33]和Gatapova等[35]的实验数据, 换算得到无量纲参数, 如无特殊说明, 各参数取值如下: ε = 0.1, A = 0.39, B = 0.001, m = 3, $ {\sigma _{{\text{ls}}, 0}} $=1, C = –0.004, Bo = 0.5, β = 1 × 10–5, γ0 = 0.1, Ca = 0.1, K = 10, E = 1 × 10–3, Ωlg = 0.003.

    • Gatapova等[35]和Chu等[36]通过测量不同表面上多种液滴的接触角迟滞特征发现, 迟滞角Δθ范围为0—90°. 为保证θa > 0和θr > 0, 代入(30)式和(31)式, 可得0 < Δθ < 1.4; 由于本文初始液滴的表观接触角为1.1, 为保证蒸发动态过程的完整呈现, 则有0 < Δθ < 0.6. 为探究不同迟滞角对液滴运动过程的影响, 下面选取Δθ = 0.27, 0.35和0.40开展分析, 其他参数保持不变, 因液滴运动过程中始终保持在x = 0两侧对称, 下面以右侧特征参数来阐明其动态过程.

      图3图5为液滴演化过程的动态特征. 图3表明, 三种情形下的接触角、接触线及其移动速度变化趋势大致相同, 液滴演化过程均可分为液滴铺展、接触线钉扎和去钉扎三个阶段(图3(b)图3(c)), 但各阶段的参数变化和持续时间却有明显不同, 各阶段开始时刻及持续时间如表1所列, 图3中以Δθ = 0.27为例用虚线划分各个阶段.

      图  3  不同Δθ时液滴演化过程中的特征参数变化 (a) 接触角; (b) 接触半径; (c) 接触线移动速度; (d) 蒸发剩余质量

      Figure 3.  Characteristic parameters of drop evolution with time at different Δθ: (a) Dynamic contact angle; (b) contact line; (c) contact line velocity; (d) remaining mass.

      图  5  动态接触角与接触线移动速度关系图

      Figure 5.  Relation between the dynamic contact angle and velocity of contact line.

      Δθ钉扎时刻去钉扎时刻蒸干时刻铺展阶段钉扎阶段去钉扎阶段
      0.271891039750780000—1891018910—3975039750—78000
      0.351653043250762000—1653016530—4325043250—76200
      0.401462045100746000—1462014620—4510045100—74600

      表 1  不同Δθ对液滴演化时间的影响

      Table 1.  Effect of varied Δθ on the drop evolution time.

      表1可知, Δθ = 0.27下液滴开始钉扎和去钉扎的时刻分别为t = 18910和t = 39750, Δθ = 0.35下开始钉扎和去钉扎的时刻为t = 16530和t = 43250, 而Δθ = 0.40下则为t = 14620和t = 45100, 即接触线钉扎开始时刻提前, 去钉扎时刻延后. 由此可知, 增大迟滞角, 延长了接触线钉扎时间, 缩短了铺展阶段和去钉扎阶段时长, 同时加快了液滴蒸干历程.

      在铺展阶段, 三种迟滞角下的液滴均表现出接触角不断减小至前进接触角(图3(a)), 接触线向外铺展(图3(b)), 且铺展速度逐渐减慢(图3(c)). 不同的是, 随着迟滞角的增大, 液滴铺展半径和铺展速度减小, Δθ = 0.40时的最大铺展半径仅为1.00417, 远小于Δθ = 0.27和Δθ = 0.35情形. 该阶段液滴铺展直径的增长率通常可表示为(xcr xcl)与t之间的幂律关系, 三种情形下分别为~t 0.004, ~t 0.002和~t 0.001, 即接触角迟滞降低了铺展速率, 抑制了铺展过程. 另外, Δθ = 0.4情形下的接触角始终大于其他两种情形, 液滴整体的铺展情况受到了显著抑制.

      在钉扎阶段, 三种迟滞角下液滴的接触线均钉扎在最大接触半径处不动, 接触角从前进接触角减小至后退接触角, Δθ = 0.4情形下的液滴接触角仍高于其他两种情形下的动态接触角, 这是受到铺展阶段的接触角迟滞影响. Bormashenko等[37]通过实验发现, 水滴在迟滞角度约为50°的铝表面的蒸发, 主要以钉扎方式进行. 在本文中, 钉扎阶段的液滴质量快速减少(图3(d)), 且Δθ = 0.4情形下液滴质量减少了0.59, 占初始液滴总质量的44.36%, 高于Δθ = 0.27和Δθ = 0.35情形的30.83%和39.10%, 即增大接触角迟滞后, 钉扎阶段的液滴蒸发在整个液滴演化过程中起重要作用, 这一结果与Bormashenko等[37]的实验结论在定性上一致, 且迟滞角越大其影响越明显.

      在去钉扎阶段, 三种情形下的接触角呈现不断减小的趋势, 液滴不断收缩, 接触线收缩速率均不断增加. 通过计算拟合得到Δθ = 0.27, 0.35和0.40情形下的液滴收缩直径(xcrxcl)与t的幂律关系分别为~t –9.06, t –8.61t –7.14, 即较大的接触角迟滞加快了液滴收缩速率, 进一步缩短此阶段时长.

      图3(d)表明, 接触角迟滞加快了液滴的蒸干过程. 当t = 74000时, Δθ = 0.40情形下的液滴接近蒸干, 此时液滴质量为0.00479; 而Δθ = 0.27和0.35时的液滴质量为0.01805和0.01009, 分别约是Δθ = 0.27情形下的3.77倍和2.11倍, 即随着迟滞角的增大, 液滴剩余质量减小(图3(d)), 液滴蒸干过程加快.

      Brutin等[38]认为液滴蒸发速率与润湿半径和接触角密切相关, 即润湿半径越大、接触角越小, 其蒸发速率越快. 接触角迟滞加快蒸干的原因可归结为: 增大迟滞角使去钉扎阶段的接触角减小, 接触线附近的空气流通阻碍减弱, 同时接触半径增大, 固-液之间接触面积增大, 加强传热, 进而使得液滴蒸干更快.

      液滴演化过程如图4所示, 由图4(a)可知, 在铺展阶段, 随着迟滞角的增大, 液滴铺展半径缩小, 液滴厚度增大, 但因该阶段液滴表面温度较低, 蒸发速率较慢, 不同迟滞角下的液滴轮廓差异并不显著. 图4(b)图4(d)表明, 在去钉扎阶段, 当迟滞角增大时, 液滴收缩速度加快, 液滴厚度减小. 随着时间的推移, 液滴表面温度提高, 蒸发速率也不同程度地增大, 液滴轮廓逐渐显出差异; 由图4(d)可知, t = 74000时, Δθ = 0.27, 0.35和0.40下的液滴收缩半径分别为0.2798, 0.2579和0.2352, 液滴高度分别为0.04825, 0.0291和0.01496, 即迟滞角越大, 液滴收缩半径和高度都越小.

      图  4  液滴演化过程轮廓图 (a) t = 10000; (b) t = 50000; (c) t = 70000; (d) t = 74000

      Figure 4.  Drop profile during the drop evolution process: (a) t = 10000; (b) t = 50000; (c) t = 70000; (d) t = 74000.

      图5比较了三种迟滞角下动态接触角与接触线移动速度间的对应关系. 可以看出, 迟滞角对铺展、钉扎和去钉扎阶段等全过程均有不同程度的影响, 在钉扎阶段, Δθ = 0.27下的液滴前进接触角和后退接触角分别为θa = 0.966和θr = 0.696, Δθ = 0.35和0.40下液滴对应于θa = 0.997, θr = 0.647和θa = 1.016, θr = 0.616. 如表2所列, 随着迟滞角的增大, 前进接触角增大, 后退接触角减小, 后退接触角减小的幅度大于前进接触角增大的幅度, 这与Karapetsas等[21]的研究结果一致. 因此, 可通过改变迟滞角来改变前进、后退接触角, 控制接触线钉扎时间, 从而调控液滴蒸发.

      ΔθθaθrΔθaΔθr
      0.270.9660.6960.106–0.164
      0.350.9970.6470.031–0.049
      0.401.0160.6160.019–0.031

      表 2  不同迟滞角下θaθr 变化

      Table 2.  Variation in θa and θr of drop with different Δθs

    • 为探讨气-液界面张力对温度的敏感性(以下简称温度敏感性)对液滴蒸发的影响, 将Ωlg由0.003增至0.005和0.01, 其他参数不变, Δθ分别设置为0和0.27. Δθ = 0表示不考虑接触角迟滞的情形. 图6给出了有无迟滞两种情况下, 不同温度敏感性系数对液滴蒸发时的特征参数变化的影响. 图6中以Ωlg = 0.003, Δθ = 0.27为代表用虚线划分各个阶段, 可以清晰地看出, 考虑接触角迟滞与否对液滴蒸发是有影响的, 以Ωlg = 0.003的情况为例, 在液滴铺展阶段, Δθ = 0时的液滴铺展半径变化速率明显高于Δθ = 0.27的情况, 接触角迟滞显著降低了铺展速度, 抑制了铺展过程, 此时Δθ = 0.27时的液滴铺展半径相对Δθ = 0.27情形较小, 而Δθ = 0.27情形下的接触角始终大于Δθ = 0的接触角, 液滴整体的铺展情况受到了显著抑制. 受接触角迟滞的影响, 液滴钉扎时间显著增大, 接触角仍然较大. 在去钉扎阶段, 接触角迟滞加快了液滴收缩速率和接触角减小速率, 从而加快蒸发.

      图  6  液滴演化过程中的特征参数变化 (a) 接触角; (b) 接触半径; (c) 接触线移动速度; (d) 蒸发剩余质量

      Figure 6.  Characteristic parameters of drop evolution with time: (a) Dynamic contact angle; (b) contact line; (c) contact line velocity; (d) remaining mass.

      三种温度敏感性下的液滴演化历程如图6所示, 接触角、接触线及其移动速度变化趋势大致相同. 图6(a)表明, 提高Ωlg后液滴在三个阶段的动态接触角均减小, 而液滴铺展半径(图6(b))和铺展速度(图6(c))增大, 致使钉扎时Ωlg = 0.01的液滴接触线更大, 液滴蒸发更快. 在去钉扎阶段, 液滴的收缩速度随Ωlg的增加而加快, 接触线半径快速减小, 但此阶段Ωlg = 0.01的液滴接触角远小于其他两种情形, 液滴质量(图6(d))仍快速减小, 直至t = 71000完全蒸干, 比Ωlg = 0.003和Ωlg = 0.005的蒸干时间分别快约10%和7%.

      液滴铺展、钉扎和去钉扎三个阶段的重要时间参数及角度参数如表3所列. 与增加接触角迟滞影响不同的是, Ωlg增加后, 液滴的θe, θaθr均减小, 强化液滴在壁面上的润湿效果, 液滴传热也得以加强, 液滴蒸发加快. 接触线开始钉扎和去钉扎时刻延后, 这是由于达到较小的θaθr需要一定时间. 在3.1节中, 增加Δθ时, θa增大, θr减小, 液滴蒸发同样加快, 因此通过增加接触角迟滞, 亦或提高气-液界面张力对温度的敏感性, 均可以减小后退接触角, 从而实现液滴的快速蒸干. 这一结论与Kulinich和Farzaneh[15]所提出的固着水滴的蒸发行为受润湿滞后(或后退接触角)控制是相符的.

      Ωlg钉扎时刻去钉扎时刻蒸干时刻θeθaθr
      0.0031891039750780000.8600.9660.686
      0.0052028040140767000.8370.9420.672
      0.012220040260710000.7750.8800.610

      表 3  不同Ωlg时液滴演化过程中参数变化

      Table 3.  Varieties of parameters during drop evolution with different Ωlg.

    • 接触角迟滞对液滴蒸发有重要影响, 可加快液滴蒸发进程. 接触角迟滞对液滴蒸发的三个阶段(液滴铺展、接触线钉扎和去钉扎阶段)均有不同程度的影响, 考虑接触角迟滞后, 接触线钉扎时间延长, 铺展阶段和去钉扎阶段时长缩短. 在液滴铺展阶段, 迟滞角的增大使液滴铺展速度降低, 铺展半径减小, 铺展受到抑制; 接触角迟滞可延长接触线钉扎时间, 维持一定的固-液接触面积, 增强传热, 从而加快蒸发. 此外, 接触角迟滞使去钉扎阶段的动态接触角减小, 液滴整体呈扁平状, 气-液界面接触面积大, 加快蒸发进行, 由此液滴收缩速率增加, 进一步缩短去钉扎阶段的时长.

      接触线开始钉扎和去钉扎时液滴分别达到前进接触角和后退接触角, 随着迟滞角的增加, 前进接触角增大, 后退接触角减小, 且后退接触角减小的幅度大于前进接触角增大的幅度. 通过改变接触角迟滞有利于调控液滴运动, 从而控制蒸发.

      与增加接触角迟滞不同的是, 气-液界面张力敏感性系数增加后, 液滴的平衡接触角、前进接触角和后退接触角均减小, 强化了液滴在壁面上的润湿效果, 液滴传热也得以加强, 液滴蒸发加快. 结合液滴蒸发过程的动态特性, 通过增加接触角迟滞或气-液界面张力对温度的敏感性将减小后退接触角, 进而加快液滴蒸干进程.

参考文献 (38)

目录

    /

    返回文章
    返回