过冷水-壁面接触面积对冰成核行为影响实验研究
Ice Nucleation Behavior in Supercooled Water with Varying Wall Contact Area
通讯作者: 王福新,研究员,博士生导师;E-mail:fuxinwang@sjtu.edu.cn.
责任编辑: 王一凡
收稿日期: 2023-11-10 修回日期: 2023-12-26 接受日期: 2023-12-29
| 基金资助: |
|
Received: 2023-11-10 Revised: 2023-12-26 Accepted: 2023-12-29
作者简介 About authors
陆晨律(1999—),硕士生,从事飞机结冰研究.
认知过冷水与固体表面接触对冰成核过程的影响对过冷物质储运以及防冰表面设计等工程应用非常关键,而接触面积变化对冰成核的影响尚不明确,尤其是当接触面积较大时.通过改变硅胶管和聚氯乙烯(PVC)纤维增强软管的长度改变过冷水与壁面接触面积的方式,研究了接触面积对结冰温度的影响.实验结果表明,接触面积变化将显著影响过冷水的结冰温度.根据所得变化规律建立了接触面积-过冷度预测模型,并利用实验结果对经典成核理论中的面积项进行了修正.修正后的成核率预测结果与实验中面积增长对成核的影响基本一致,且接触面积对成核率的影响并非目前认为的线性,而呈非线性.研究方法和所提出的模型为未来工程中过冷水的利用以及接触面积变化影响成核过程的进一步研究奠定了良好的基础.
关键词:
Understanding the influence of the contact area between supercooled water and surface on the nucleation is critical for engineering applications such as anti-icing surface design and the supercooling preservation. However, the effect of contact area between supercooled water and surface on the nucleation still remains unclear, especially when the contact area is large. Therefore, the influence of the contact area on the freezing temperature of supercooled water inside the hose was experimentally studied by changing the length of silicone and polyvinyl chloride (PVC)-reinforced hoses. The results show that variations in contact area significantly affect the freezing temperature of supercooled water. Based on experimentally observed trends, a predictive model for contact area-supercooling relationship was developed. Furthermore, the area term in the classical nucleation theory (CNT) was modified using the experimental data. The modified nucleation rate predictions align well with the experimental findings regarding the influence of increasing contact area on nucleation. Notably, the effect of contact area on nucleation rate was found to be nonlinear, contrary to the commonly assumed linear relationship. The research methodology and the proposed model in this paper provide a solid foundation for future applications of supercooled water in engineering, as well as for further investigation into the role of contact area in nucleation phenomena.
Keywords:
本文引用格式
陆晨律, 王利平, 孟航飞, 刘洪, 王福新.
LU Chenlü, WANG Liping, MENG Hangfei, LIU Hong, WANG Fuxin.
过冷是指液体温度低于熔点仍保持液态的一种亚稳定状态,其普遍存在于人类生活与自然界中.研究表明在一些温度低至-37.5 ℃的云层中还存在过冷的液态水[1].水在过冷状态发生凝固需首先经历冰成核过程,成核又分为均质成核和异质成核[2].均质成核是在深过冷度下热力学波动的结果,其他情况下则存在壁面及杂质等作为凝结核诱发异质成核,这是由于凝结核的存在降低了成核过程的能垒[2-3].过冷和成核与很多应用领域相关,比如过冷能够延长生物细胞与器官的保存时间[4],过冷保存能维持果蔬与肉制品等食品的品质[5],过冷水能够用于快速制冰[6],输电线、飞机及风力机等表面需要避免成核发生来实现防冰[7],等等.同时有些应用领域需要避免出现过冷甚至主动触发成核,比如应用相变材料进行热能储存[8].过冷导致物质性质的异常,且成核机理非常复杂[9],目前还未被完全认识.因此理解存在不同外界扰动时的成核行为是促进相关应用与技术发展的重要基础.
在经典成核理论中,异质成核率与接触面积成正比[17].成核率可以通过改变接触面积来改变,比如超疏水表面接触角越大,水滴与表面间的接触面积就越小,从而可以在不改变成核能垒的条件下降低成核率[18].根据经典成核理论,通常条件下温度每变化1 ℃,成核率会呈现数个量级的变化[3,19-20].因此,试验条件相同时,理论上接触面积应存在量级上的改变才能使结冰温度变化1 ℃.然而,根据Zhang等[21]对铝表面单液滴冰成核特性实验研究结果,等效直径2.12 mm和2.67 mm水滴平均成核温度分别为-17.04 ℃ 和-15.82 ℃,此时水滴与表面的接触面积分别为6.11 mm2和9.66 mm2,比值为1∶1.59,并未出现量级上的变化.另外,Inada等[22]实验研究了玻璃表面水滴的冰成核规律,认为当水与玻璃间的接触面积大于0.01 mm2后,经典成核理论不再适用于冰成核预测.在一些实际应用中,过冷物质与表面的接触面积会远大于单个水滴与壁面的接触面积,比如过冷保存生物样品容积可达数百毫升[4],有些用于储能的相变材料胶囊直径在10 cm左右[23],未来结冰风洞模拟过冷大水滴云雾结冰条件也会需要大量过冷水在管道内进行输运[24].但目前,过冷水体积较大时过冷水-壁面接触面积对结冰的影响规律尚不清晰.
为了理解过冷水-壁面接触面积对过冷水结冰温度的影响,本文开展了过冷水-壁面接触面积较大条件下其对冰成核过程影响的实验研究,建立接触面积与成核温度间的关系,结合实验统计和经典成核理论,对过冷水与不同面积的不同材料表面接触时的异质成核行为进行分析,并根据实验结果对经典成核理论中成核率预测公式进行修正,以实现实际应用系统中过冷水与壁面接触面积较大时壁面诱导冰成核的成核率与成核温度的预测.
1 实验方法
图1
制冷设备使用的是一台高低温交变湿热试验箱(LRHS-504B-LJS, Shanghai Linpin Instrument, Ltd.).其温度控制范围为-40~150 ℃,控温精度±0.5 ℃,箱内温度均匀度小于2 ℃,最大降温速率约为3.0 ℃/min.
实验测试部分由实验样本(注入超纯水的圆管)与相应的支撑架组成.实验所用的管道为工程中常用的硅胶管与聚氯乙烯(PVC)纤维增强软管,内径分别为4 mm和8 mm,壁厚分别为1 mm和2 mm,长度范围分别为0.5~8 m和0.2~4 m.为了尽可能全面地反映过冷水结冰温度随接触面积的变化趋势,硅胶管的长度选取为0.5、0.65、0.8、1、2、5和8 m,PVC管的长度为0.2、0.25、0.5、1、2.5和 4 m.对应的过冷水-壁面接触面积范围分别为6.28×10-3~1.01×10-1 m2和 5.03×10-3~1.01×10-1 m2,过冷水体积范围分别为6.28~100.53 mL(硅胶管)和10.05~201.06 mL(PVC管).为了验证相同接触面积条件下过冷水体积变化对结冰温度的影响,选取内径2、4、6和8 mm的硅胶管进行了实验,对应的管道长度分别为2、1、0.67和0.5 m.受试验箱工作室尺寸(700 mm×800 mm×900 mm)限制,支撑架上每次最多可安装3条实验管道同时进行实验.为了降低冰成核的随机性对实验结果的影响,每个实验条件至少测试20个实验样本.需注意的是在将管路安装在支撑架上时应避免管道弯折.图2展示了实验所用的两种管道,图2(a)为实物图.利用原子力显微镜(atomic force microscopy, AFM)对两种管道内壁面的表面形貌与粗糙度进行了测量,结果如图2(b)、2(c)所示,测得硅胶管与PVC管内壁面的表面粗糙度Ra分别为17 nm和62 nm.实验中所用的超纯水由Rephile Bioscience, Ltd.的RS2200QUV制备,其电阻率为18.2 MΩ·cm@25 ℃.
图2
测温系统由一台温度测量仪(THM200K, TENGHUI)和PT100温度传感器组成,测温精度为±0.1 ℃,测量频率为1 Hz.为标定管道内外壁面温差,实验前将两个PT100分别粘贴在管道同一位置的内外壁面上,在降温过程中对两种管道的内外壁面温度分别进行了测量,温度曲线如图3所示.降温过程中,管道内壁面温度高于外壁面,管道内外壁面温差Td基本稳定,硅胶管为0.1 ℃,管壁较厚的PVC管为0.8 ℃.结冰温度测量实验中每条管道外壁面至少粘贴两个温度传感器,同一管道上所有传感器温度的平均值记为该实验样本外壁面的温度,再根据实验前标定的内外壁面温差进行修正,得到管道内过冷水的平均温度.过冷水结冰过程中会释放潜热使温度迅速上升至熔点,因此温度上升表明过冷水内已经发生冰成核并开始结冰,温度上升前过冷水的温度则视为结冰温度,亦为冰成核温度TN,如图3所示.图中:Tsi-in和Tsi-out分别表示硅胶管内壁面与外壁面的温度;TPVC-in和TPVC-out分别表示PVC管内壁面与外壁面的温度.
图3
为了降低三相界面对过冷水成核的影响,本文通过对管道两端进行加热提高端口处温度来降低此处成核发生的概率,故设计了端口控温系统.该系统由直流电源、缠绕在管道外侧的电阻丝以及外部包裹的保温棉组成.在实验开始降温时对电阻丝施加电流进行加热,使端口处的温度始终高于管道温度至少3 ℃.对两种管道端口加热装置的影响范围分别进行了测量,结果显示:对于硅胶管,距电阻丝 5 cm 处与管路中段处温差小于0.1 ℃;对于PVC软管,这个距离是2 cm.通过这种方法,实验中未在端口处发生成核.
对于每次实验,首先将超纯水注入管内,将端口弯折并用鱼尾夹夹住防止水外流.需要注意的是,在往管内注水过程中要避免出现气泡.之后将实验样本放入高低温试验箱并开始降温,同时开启直流电源控制端口温度.实验样本从室温开始降温直到所有样本发生结冰,为节约单次实验时间,当温度高于5 ℃时,降温速率设定为1.5 ℃/min,当温度降低到5 ℃后,降温速率设定为0.3 ℃/min直到本次实验结束.实验结束后,需等待管路中的冰融化后将水排出,再重新注入超纯水重复上述步骤进行下次实验.
2 结果与讨论
2.1 成核温度分布
过冷水在不同长度硅胶管和PVC管内的结冰温度如图4所示.图中:l为各组实验中管道长度.在实验条件完全相同的情况下,各次实验的结冰温度具有随机性,但仍在一定范围内.一般而言,管道材质相同时,管道长度越短,结冰温度就越低.在本文实验的长度范围内,硅胶管内过冷水的结冰温度分布为-12~-21 ℃,PVC管内的过冷水结冰温度分布为-12~-20 ℃.
图4
图4
过冷水在不同长度管内形成不同接触面积条件下的TN
Fig.4
TN of the supercooled water under different contact areas formed by different lengths
图5给出了硅胶管与PVC管内过冷水平均结冰温度TAVG、中值结冰温度TMED与l之间的变化趋势.图中TAVG的误差线为平均绝对偏差.硅胶管与PVC管的TAVG、TMED的变化趋势较为相似,都大致呈现随l增加先迅速增长,随后增长速度放缓并在实验范围内逐步趋向于定值.对于硅胶管,以2 m为分界点,当l<2 m,l的增加会导致TAVG、TMED的迅速上升,当l>2 m后,l的进一步增加不再使TAVG、TMED发生明显的变化,TAVG、TMED在实验所测试的长度范围内逐步趋向于稳定.对于PVC管而言,这个分界点为1 m.平均绝对偏差的变化总体上呈现一个随l增加而略微减小的趋势,这可能是l增大后成核的随机性减弱的结果.
图5
当l发生变化时,A相应改变,同时管道内过冷水的体积V也将发生变化.根据经典成核理论,温度为260 K时,发生在水内部的均质成核率仅有约10-150 cm-3·s-1[29],因此,理论上本文实验中水体积的改变并不会对结冰温度造成影响.为验证,本文进行了保持A不变而改变V的实验.图6展示了实验中TN的散点图,图7展示了不同体积过冷水的TAVG.可以看出,当A相同,V从6.28 mL变化到25.13 mL时,TAVG的变化很小.V相同,A不同(内径din=4 mm、l=0.5 m与din=2 mm、l=2 m)的硅胶管内的过冷水,TAVG相差了4.1 ℃.这表明在本文所进行的实验中,A的变化相比于V的变化对结冰温度有更显著的影响.因此,尽管管道长度的变化会同时导致A与V的变化,本文仅考虑A变化对成核带来的影响.
图6
图6
过冷水在A相同(0.012 6 m2) V不同时的TN
Fig.6
TN of the supercooled water with the same contact area of 0.012 6 m2 and different volumes
图7
2.2 过冷度与接触面积间的关系
后续将以A代替l分析成核过程.由图5可见,当过冷水A达到一定大小后,A的进一步增加对冰成核的促进作用将不再明显,从而使TAVG趋向于一个定值.另外,管道材质的差异使A的增长对成核的促进程度存在一定的差异,在建立A与平均过冷度ΔTAVG之间关系时需要合理设置参数以体现材料之间的差异.已有研究表明,当撞击物自由下落撞击壁面触发过冷水结冰时,随着下落高度的增加,过冷水结冰温度迅速上升,且上升速度逐渐放缓,结冰温度在下落高度达到一定值后不再发生变化[16,30].本文中A的变化对ΔTAVG的影响与上述研究中结冰温度的变化趋势有一定的相似之处,因此本文参考以上研究中的关系模型建立过冷水-壁面接触面积A与过冷水结冰平均过冷度ΔTAVG两者之间的关系:
图8
表1 平均成核过冷度与接触面积关系式拟合参数
Tab.1
| 管道材质 | a | b | c |
|---|---|---|---|
| 硅胶 | -13.525 11 | -0.002 82 | 0.020 78 |
| PVC | -13.368 30 | 0.002 38 | 0.028 76 |
需要指出的是,理想情况下,当过冷水与外界物质接触面积为0时,结冰始于均质成核,根据过冷水的体积,成核温度在-35~-40 ℃[31].因此,当过冷水与壁面接触面积趋近于0时,结冰温度也应当趋近一个有限的值.根据式(1)计算,在A趋向于0时,PVC管的ΔTAVG最终会达到25.0 ℃,硅胶管的拟合结果则出现了奇点,这是由于实验中在硅胶管A接近0的过程中TN迅速下降,进而导致b为负值.由于本文实验方法的限制,未能实现更小接触面积下过冷水结冰温度的测量实验,所以根据实验结果建立的预测模型的适用范围应大于本文实验最小面积,即约0.01 m2.当面积小于0.01 m2时,本文所建立的模型需进一步验证或修正.
2.3 基于经典成核理论考虑接触面积变化的成核行为分析
2.3.1 基于成核温度分布的成核率估算
图9
图9
成核次数与成核的累积概率随温度变化关系
Fig.9
Frequency of nucleation and cumulative probability
图10
γiw为冰-水界面能,单位为J/m2.在237.15~273.15 K范围内可以表示为[19]
ΔHV为冰的体积熔解焓,为306 MJ/m3[26].
表2 成核率拟合参数
Tab.2
| 管道材质 | lg J0 | f |
|---|---|---|
| 硅胶 | 0.441 52 | 0.014 73 |
| PVC | 0.771 53 | 0.017 85 |
2.3.2 基于接触面积修正的成核率预测
在2.3.1节的分析中,可以看到当过冷水与壁面间的接触面积大范围变化时,成核率拟合结果仍具备较好的符合性,但由于在分析时将面积标准化,最终主要体现是管道材质不同对过冷水内冰成核的影响,未能体现过冷水与壁面接触面积变化带来的影响.为了进一步分析接触面积变化对成核的影响,本节将通过具体计算成核率的变化来将其影响定量化,并据此对经典成核理论成核率公式中的面积项做出一定的修正,以使本文实验的结果与其结果一致.
E与TR为常数,对于150~273 K温度范围内的液态水,E=892 K,TR=118 K.
异质成核是一个随机事件,可以使用随机过程理论的方法进行研究,当温度为T,持续时间为δt时,样本发生结冰的概率[2]为
参考Shardt等[31]的参数设定方法,根据本文实验中的降温速率,取δt=10 s.在结冰概率已知的条件下,结合式(8)和(11)以及Fletcher因子f,可以建立成核温度与接触面积之间的一一对应关系.但现有技术条件测量冰核-水-凝结核间杨氏接触角的值较为困难,目前一般认为f的值与表面材质有关[36].因此,本文分别将硅胶管和PVC软管实验中最小接触面积时的TMED与A代入式(8)和(11),来计算两种管路材质表面的f值.计算得到f值后,可以根据式(8)和(11)分别建立硅胶管与PVC管对应的TMED与A的关系.这里计算得到的面积值与同一温度下实验中的面积大小不一定相等,因此将根据式(8)和(11)计算出的面积定义为虚拟面积AV.表3和4分别展示了硅胶管和PVC管实验中的真实接触面积Aexp、该面积对应的TMED与该温度下AV的值.
表3 硅胶管的Aexp、TMED与AV
Tab.3
| Aexp/m2 | TMED/℃ | AV/m2 | ln AV |
|---|---|---|---|
| 6.28×10-3 | -19.20 | 6.28×10-3 | -5.07 |
| 8.17×10-3 | -18.43 | 3.63×100 | 1.29 |
| 1.01×10-2 | -15.75 | 5.30×1013 | 31.60 |
| 1.26×10-2 | -15.65 | 2.32×1014 | 33.08 |
| 2.51×10-2 | -13.90 | 1.34×1028 | 64.77 |
| 6.28×10-2 | -13.67 | 2.36×1030 | 69.94 |
| 1.01×10-1 | -13.93 | 6.99×1027 | 64.11 |
表4 PVC管的Aexp、TMED与AV
Tab.4
| Aexp/m2 | TMED/℃ | AV/m2 | ln AV |
|---|---|---|---|
| 5.03×10-3 | -17.50 | 5.03×10-3 | -5.29 |
| 6.28×10-3 | -18.00 | 1.02×10-4 | -9.19 |
| 1.26×10-2 | -15.55 | 1.32×106 | 14.10 |
| 2.51×10-2 | -13.95 | 1.30×1016 | 37.10 |
| 6.28×10-2 | -13.78 | 3.00×1017 | 40.24 |
| 1.01×10-1 | -13.90 | 7.31×1016 | 38.83 |
图11
表5 面积修正模型的拟合参数
Tab.5
| 管道材质 | u | v |
|---|---|---|
| 硅胶 | 77.13 | -0.53 |
| PVC | 42.62 | -0.28 |
图11中的点表示根据经典成核理论计算得到的值,实线为拟合曲线.随着实验中Aexp的增大,ln AV迅速增长,且增长速度逐渐放缓.当Aexp趋向0时,ln AV趋向负无穷,即AV趋向0,拟合模型的结果与实际相符.
将式(12)代入式(8),并对其中的指数项进行整理,可以得到修正后的异质成核率预测公式:
为了说明修正面积项对成核率量级的影响,本文选择PVC管在温度为-15 ℃条件下,计算了A变化时式(13)与式(8)计算值的差异,如图12所示.黑色实线与红色虚线分别代表式(8)与式(13)计算得到的结果.A从0.01 m2变化到0.20 m2过程中,根据式(8),A增长对成核率的影响很小,成核率的量级从10-14 s-1增长到10-11 s-1,只改变了3个量级;而根据面积项修正后的计算结果,A的增加会令成核率在量级上迅速增长,导致成核发生并结冰,与本文实验情况相符.
图12
图12
根据经典成核理论应用虚拟接触面积和真实接触面积预测的成核率随A的变化(PVC,-15 ℃)
Fig.12
Variation of nucleation rate with contact area predicted by CNT and modified CNT using virtual and actual contact areas (PVC, -15 ℃)
对成核率分析的结果表明,接触面积变化对成核所造成的影响十分明显.Niedermeier等[37]通过假设凝结核表面存在多个成核区域,各区域内成核性质符合特定的分布函数,利用统计模拟建立了对应的成核模型,初步解释了造成成核行为随机性与奇异性的原因.同样的,实际中与过冷水接触的材料表面上各个区域的性质(如杨氏接触角)并不一定相同,从而会影响不同区域成核发生的概率.根据现有技术水平,还不能准确测量得到接触表面上成核点的分布及数量,另外成核点位的性质是否会随着温度发生改变也无法确定,这些性质的差异最终有可能体现为随着材料表面积增加成核率并非线性增加,从而影响接触面积增加时过冷水可以达到的最低温度.目前,对于与壁面接触的过冷水,基本上只能通过数值模拟或者实际试验积累数据后,再根据拟合或数据分析的方法对过冷水的成核行为进行预测[38-39].未来需要发展纳微尺度成核探测技术,再结合先进的表面表征技术对表面成核开展更深入的研究.
3 结论
本文通过实验研究了不同长度硅胶管和PVC纤维增强软管中过冷水与壁面接触面积对过冷水结冰温度影响,得出以下主要结论:
(1) 在本文实验接触面积范围内,当A较小时,过冷水TAVG会随着A增加迅速上升,且上升速度逐渐放缓,当A增加到一定值后,TAVG基本不再随A增加而变化.尽管管道长度的变化会同时导致过冷水体积的变化,但在本文的实验条件下,实验证明其对冰成核的影响很小.
(2) 通过对实验数据拟合,建立了A与ΔTAVG间的关系模型.基于经典成核理论开展了理论分析,使用经典成核理论中的成核率预测模型对实验获得的成核率进行拟合,发现即使实验样本A变化范围较大,面积均一化后的成核率仍具备较好的拟合结果.
(3) 为理解A变化对实验中过冷水内冰成核的影响,本文定义了过冷水与壁面间虚拟接触面积AV,根据经典成核理论建立了Aexp与AV间的关系,发现本文实验中冰成核率随接触面积的变化并不是通常认为的线性,而是非线性的,并建立了基于面积修正的成核率预测方法.
(4) 无论是A-ΔT还是A-AV关系模型,总体变化趋势与表面材质无关,但关系模型中的拟合参数与表面材质有关.
本文结果能够为表面接触面积变化对成核的影响提供新的见解,同时也能够为大体积过冷相关设施设备的设计提供一定指导.
参考文献
Deep convective clouds with sustained supercooled liquid water down to -37.5 ℃
[J].
A review of classical and nonclassical nucleation theories
[J].
Long-term deep-supercooling of large-volume water and red cell suspensions via surface sealing with immiscible liquids
[J].
DOI:10.1038/s41467-018-05636-0
PMID:30097570
[本文引用: 3]
Supercooling of aqueous solutions is a fundamentally and practically important physical phenomenon with numerous applications in biopreservation and beyond. Under normal conditions, heterogeneous nucleation mechanisms critically prohibit the simultaneous long-term (> 1 week), large volume (> 1 ml), and low temperatures (< - 10 degrees C) supercooling of aqueous solutions. Here, we report on the use of surface sealing of water by an oil phase to significantly diminish the primary heterogeneous nucleation at the water/air interface. We achieve deep supercooling (down to - 20 degrees C) of large volumes of water (up to 100 ml) for long periods (up to 100 days) simultaneously via this approach. Since oils are mixtures of various hydrocarbons we also report on the use of pure alkanes and primary alcohols of various lengths to achieve the same. Furthermore, we demonstrate the utility of deep supercooling via preliminary studies on extended (100 days) preservation of human red blood cells.
Sub-zero temperature preservation of fruits and vegetables: A review
[J].
Research status of supercooled water ice making: A review
[J].
Icephobic surfaces: Definition and figures of merit
[J].
DOI:S0001-8686(18)30331-2
PMID:31096074
[本文引用: 1]
Icephobic surfaces have a critical footprint on human daily lives ranging from aviation systems and infrastructures to energy systems, but creation of these surfaces for low-temperature applications remains elusive. Non-wetting, liquid-infused and hydrated surfaces have inspired routes for development of icephobic surfaces. However, high freezing temperature, high ice adhesion strength and subsequent ice accretion, low mechanical durability, and high production cost have restricted their practical applications. In this review, we provide a comprehensive definition for icephobicity through thermodynamics, heat transfer and mechanics of ice/water-material interface and elucidate physic-based routes through which nano-scale could help to achieve exceptional icephobic surfaces. Based on conservation laws, mathematical models are developed that accurately predict ice growth rate on various substrates and wind conditions. Through physics of fracture at ice-icephobic material interface, we cast a standard method for ice adhesion measurement that has the potential to eliminate discrepancies between reported ice adhesion from different laboratories. To assure long-time performance of icephobic surfaces, durability metrics need to be defined. We provide standard methods to examine mechanical, chemical, and environmental durability of icephobic surfaces. In the developed comprehensive framework on icephobicity in this review, performance of state-of-the-art icephobic surfaces are compared and main deficiencies in this field are highlighted.Copyright © 2019 Elsevier B.V. All rights reserved.
Supercooling of phase change materials: A review
[J].
Understanding crystal nucleation mechanisms: Where do we stand? General discussion
[J].DOI:10.1039/d2fd90021a PMID:35789238 [本文引用: 1]
Isochoric conditions enhance stability of metastable supercooled water
[J].
Temperature-gradient effects on heterogeneous ice nucleation from supercooled water
[J].
Suppression of ice nucleation in supercooled water under temperature gradients
[J].
Fundamental research on external factors affecting the freezing of supercooled water
[J].
Review on the control of ice nucleation by ultrasound waves, electric and magnetic fields
[J].
Enhanced control over ice nucleation stochasticity using a carbohydrate polymer cryoprotectant
[J].
Nucleation in supercooled water triggered by mechanical impact: Experimental and theoretical analyses
[J].
Control of ice nucleation: Freezing and antifreeze strategies
[J].
DOI:10.1039/c8cs00626a
PMID:30137078
[本文引用: 1]
Water freezing remains a perennial topic of great relevance to many important aspects of our lives; from the climate to human society and from economics to medicine, frozen water profoundly influences our living environment and life activities. There have been numerous publications on water freezing; however, confusion regarding the process of freezing remains. In this review, we mainly focused on the nucleation aspects of water freezing; in particular, we focused on the effect of the surface morphology and nanostructure of foreign bodies. This review covers the recent progress in ice nucleation and anti-freezing strategies within the framework of nucleation principles. In this regard, we first summarize the crystal nucleation theories. Due to high interfacial energy, ice crystallization is primarily controlled by heterogeneous nucleation events, because the homogeneous nucleation barrier of ice is extremely high. In addition to the interfacial energy, the interfacial morphology or nanostructure of foreign bodies plays a diverse role under different supercooling regimes due to the Gibbs-Thomson effect. This effect gives rise to the inverse homogeneous-like nucleation phenomenon, in which foreign bodies have little influence on the nucleation barrier. This ensures the accurate measurement of the nucleation barrier, critical size, and water-ice interfacial energy, in agreement with the latest studies based on a microemulsions approach, metadynamics, the mW model, etc. As a consequence, anti-freezing strategies can be implemented by reducing the nucleation rate through restriction of the contact area of the water/substrate interface, by increasing the heterogeneous nucleation barrier through modification of the interfacial properties of foreign particles, including the interfacial structure and the interaction between the water and foreign particles and by kink kinetics. Within this context, the anti-freezing mechanism of superhydrophobic substrates was reviewed. Therefore, it follows that by significantly reducing the contact area between the water and substrate, superhydrophobic materials can effectively reduce the heterogeneous nucleation rate. We hope that this review will provide a unified picture and guidance for future work on water freezing.
Pancake bouncing on superhydrophobic surfaces
[J].
DOI:10.1038/nphys2980
PMID:28553363
[本文引用: 1]
Engineering surfaces that promote rapid drop detachment1,2 is of importance to a wide range of applications including anti-icing3-5, dropwise condensation6, and self-cleaning7-9. Here we show how superhydrophobic surfaces patterned with lattices of submillimetre-scale posts decorated with nano-textures can generate a counter-intuitive bouncing regime: drops spread on impact and then leave the surface in a flattened, pancake shape without retracting. This allows for a four-fold reduction in contact time compared to conventional complete rebound1,10-13. We demonstrate that the pancake bouncing results from the rectification of capillary energy stored in the penetrated liquid into upward motion adequate to lift the drop. Moreover, the timescales for lateral drop spreading over the surface and for vertical motion must be comparable. In particular, by designing surfaces with tapered micro/nanotextures which behave as harmonic springs, the timescales become independent of the impact velocity, allowing the occurrence of pancake bouncing and rapid drop detachment over a wide range of impact velocities.
Classical nucleation theory of homogeneous freezing of water: Thermodynamic and kinetic parameters
[J].
DOI:10.1039/c4cp04184d
PMID:25627933
[本文引用: 1]
The probability of homogeneous ice nucleation under a set of ambient conditions can be described by nucleation rates using the theoretical framework of Classical Nucleation Theory (CNT). This framework consists of kinetic and thermodynamic parameters, of which three are not well-defined (namely the interfacial tension between ice and water, the activation energy and the prefactor), so that any CNT-based parameterization of homogeneous ice formation is less well-constrained than desired for modeling applications. Different approaches to estimate the thermodynamic and kinetic parameters of CNT are reviewed in this paper and the sensitivity of the calculated nucleation rate to the choice of parameters is investigated. We show that nucleation rates are very sensitive to this choice. The sensitivity is governed by one parameter - the interfacial tension between ice and water, which determines the energetic barrier of the nucleation process. The calculated nucleation rate can differ by more than 25 orders of magnitude depending on the choice of parameterization for this parameter. The second most important parameter is the activation energy of the nucleation process. It can lead to a variation of 16 orders of magnitude. By estimating the nucleation rate from a collection of droplet freezing experiments from the literature, the dependence of these two parameters on temperature is narrowed down. It can be seen that the temperature behavior of these two parameters assumed in the literature does not match with the predicted nucleation rates from the fit in most cases. Moreover a comparison of all possible combinations of theoretical parameterizations of the dominant two free parameters shows that one combination fits the fitted nucleation rates best, which is a description of the interfacial tension coming from a molecular model [Reinhardt and Doye, J. Chem. Phys., 2013, 139, 096102] in combination with the activation energy derived from self-diffusion measurements [Zobrist et al., J. Phys. Chem. C, 2007, 111, 2149]. However, some fundamental understanding of the processes is still missing. Further research in future might help to tackle this problem. The most important questions, which need to be answered to constrain CNT, are raised in this study.
Experimental investigation and statistical analysis of icing nucleation characteristics of sessile water droplets
[J].
Ice nucleation in water droplets on glass surfaces: From micro-to macro-scale
[J].
Encapsulation surface roughness effect on the performance of cool storage systems
[J].
Investigating the effects of solid surfaces on ice nucleation
[J].
DOI:10.1021/la3014915
PMID:22741592
[本文引用: 4]
Understanding the role played by solid surfaces in ice nucleation is a significant step toward designing anti-icing surfaces. However, the uncontrollable impurities in water and surface heterogeneities remain a great challenge for elucidating the effects of surfaces on ice nucleation. Via a designed process of evaporation, condensation, and subsequent ice formation in a closed cell, we investigate the ice nucleation of ensembles of condensed water microdroplets on flat, solid surfaces with completely different wettabilities. The water microdroplets formed on flat, solid surfaces by an evaporation and condensation process exclude the uncontrollable impurities in water, and the effects of surface heterogeneities can be minimized through studying the freezing of ensembles of separate and independent water microdroplets. It is found that the normalized surface ice nucleation rate on a hydrophilic surface is about 1 order of magnitude lower than that on a hydrophobic surface. This is ascribed to the difference in the viscosity of interfacial water and the surface roughness.
Ice nucleation behaviour on sol-gel coatings with different surface energy and roughness
[J].
DOI:10.1039/c5cp03243a
PMID:26220055
[本文引用: 5]
In this paper, the ice nucleation temperatures of 10 μL water droplets on a series of sol-gel coatings with different roughness and surface energies were obtained using a customized automatic measurement system. Classical nucleation theory was then employed to explain the different icing behaviour on the coatings. It was found that the wetting mode at low temperatures is strongly correlated with the icing behavior of the droplets on the surfaces. Ice-phobic coatings can lower the icing temperature of the droplet on the surface by up to 6.9 °C compared with non-icephobic ones. Using classical nucleation theory, our results support some recent observations that the dominant nucleation sites are along the substrate-water-vapour three-phase contact line rather than at the substrate-water interface.
Fast imaging of freezing drops: No preference for nucleation at the contact line
[J].
Faster nucleation of ice at the three-phase contact line: Influence of interfacial chemistry
[J].
Nucleation rate and supercooling degree of water-based graphene oxide nanofluids
[J].
The mechanical stimulus to crystallization
[J].
Homogeneous freezing of water droplets for different volumes and cooling rates
[J].
DOI:10.1039/d2cp03896j
PMID:36413087
[本文引用: 2]
To understand the crystallization of aqueous solutions in the atmosphere, biological specimens, or pharmaceutical formulations, the rate at which ice nucleates from pure liquid water must be quantified. There is still an orders-of-magnitude spread in the homogeneous nucleation rate of water measured using different instruments, with the most important source of uncertainty being that of the measured temperature. Microfluidic platforms can generate hundreds to thousands of monodisperse water-in-oil droplets, unachievable by most other techniques. However, most microfluidic devices previously used to quantify homogeneous ice nucleation rates have reported high temperature uncertainties, between ±0.3 and ±0.7 K. We use the recently developed Microfluidic Ice Nuclei Counter Zurich (MINCZ) to observe the freezing of spherical water droplets with two diameters (75 and 100 μm) at two cooling rates (1 and 0.1 K min). By varying both droplet volume and cooling rate, we were able to probe a temperature range of 236.5-239.3 K with an accuracy of ±0.2 K, providing reliable data where previously determined nucleation rates suffered from large uncertainties and inconsistencies, especially at temperatures above 238 K. From these data and from Monte Carlo simulations, we demonstrate the importance of obtaining a sufficiently large dataset so that underlying nucleation rates are not overestimated at higher temperatures. Finally, we obtain new parameters for a previous parameterisation by fitting to our newly measured nucleation rates, enabling its use in applications where ice formation needs to be predicted.
Two-dimensional nucleation of ice from supercooled water
[J].
Rational nanostructuring of surfaces for extraordinary icephobicity
[J].
DOI:10.1039/c3nr06644d
PMID:24667802
[本文引用: 5]
Icing of surfaces is commonplace in nature, technology and everyday life, bringing with it sometimes catastrophic consequences. A rational methodology for designing materials with extraordinary resistance to ice formation and adhesion remains however elusive. We show that ultrafine roughnesses can be fabricated, so that the ice nucleation-promoting effect of nanopits on surfaces is effectively counteracted in the presence of an interfacial quasiliquid layer. The ensuing interface confinement strongly suppresses the stable formation of ice nuclei. We explain why such nanostructuring leads to the same extremely low, robust nucleation temperature of ∼-24 °C for over three orders of magnitude change in RMS size (∼0.1 to ∼100 nm). Overlaying such roughnesses on pillar-microtextures harvests the additional benefits of liquid repellency and low ice adhesion. When tested at a temperature of -21 °C, such surfaces delayed the freezing of a sessile supercooled water droplet at the same temperature by a remarkable 25 hours.
Mechanism of supercooled droplet freezing on surfaces
[J].
Heterogeneous ice nucleation rate coefficient of water droplets coated by a nonadecanol monolayer
[J].
A simple classical model for predicting onset crystallization temperatures on curved substrates and its implications for phase transitions in confined volumes
[J].
Heterogeneous ice nucleation: Exploring the transition from stochastic to singular freezing behavior
[J].
Predicting heterogeneous ice nucleation with a data-driven approach
[J].
DOI:10.1038/s41467-020-18605-3
PMID:32963232
[本文引用: 1]
Water in nature predominantly freezes with the help of foreign materials through a process known as heterogeneous ice nucleation. Although this effect was exploited more than seven decades ago in Vonnegut's pioneering cloud seeding experiments, it remains unclear what makes a material a good ice former. Here, we show through a machine learning analysis of nucleation simulations on a database of diverse model substrates that a set of physical descriptors for heterogeneous ice nucleation can be identified. Our results reveal that, beyond Vonnegut's connection with the lattice match to ice, three new microscopic factors help to predict the ice nucleating ability. These are: local ordering induced in liquid water, density reduction of liquid water near the surface and corrugation of the adsorption energy landscape felt by water. With this we take a step towards quantitative understanding of heterogeneous ice nucleation and the in silico design of materials to control ice formation.
Crystal nucleation in liquids: Open questions and future challenges in molecular dynamics simulations
[J].
DOI:10.1021/acs.chemrev.5b00744
PMID:27228560
[本文引用: 1]
The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insights into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that, in the past few decades, have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state-of-the-art computational methods by reviewing simulations of such processes as ice nucleation and the crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insights into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that, as a result, the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that, by improving (i) existing interatomic potentials and (ii) currently available enhanced sampling methods, the community can move toward accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.
/
| 〈 |
|
〉 |
