交变电场作用下单液滴蒸发的分子动力学模拟1)

2021-05-31富庆飞

力学学报 2021年5期

王群富庆飞

(北京航空航天大学宇航学院,北京 100190)

引言

液滴的蒸发是自然界的一种常见现象,在许多领域中都起着重要的作用[1-4].1976 年,日本学者浅川勇吉发现了“浅川效应”,即把水放在高压静电场中,而电场能够显著促进水蒸发的一种现象[5].因此,很多学者的研究集中在电场对传热传质、流体流动、液滴蒸发等方面的影响上.外加电场作为流体的驱动力已经被应用在很多技术上,如静电纺丝[6]、喷墨打印[7-8]、喷雾干燥[9]、静电雾化[10-11]等.

在电场力作用下,水分子的排列结构和动力学都会发生改变,从而影响液滴的蒸发过程[12].到目前为止,外加电场对液滴蒸发过程的影响一直是研究的热点.Kiyoshi 等[13]、Hashinaga 等[14]和Li 等[15]通过实验发现,当电场垂直作用于高温壁面上的液滴时,电场对液滴蒸发有着较强的促进作用,且交变电场比匀强电场更显著地促进液滴蒸发.Heidarinejad 等[16]对直流电场作用下通道内的水蒸发进行了数值研究.结果表明,在电场作用下,水的蒸发速率增大,但在高雷诺数下,电场的作用减弱.常景彩等[17]通过数值模拟方法得到了高压静电场下液膜的蒸发速率会升高.Jaiswal 等[18]利用PIV 技术观察到在电场作用下的三维悬浮液滴,其内部的热对流被大大加强了,从而导致液滴蒸发速率加快.但Zuo 等[19]利用数值模拟方法发现,当外部静电场平行于高温壁面时,会抑制高温壁面上的液滴内部的流动,从而降低了液滴的蒸发速率,这与Almohammadi 等[20]通过实验得到的结论定性上保持一致.

分子动力学(molecular dynamics,MD)模拟是基于粒子模型的微观尺度研究方法,是从微观上分析液滴蒸发行为的有力工具.Wu 等[21]利用分子动力学方法研究了二氟甲烷纳米液滴的蒸发特性,包括环境温度、液滴尺寸、液滴初始温度等因素的影响.Vaitheeswaran 等[22]、Peng 等[23]采用MD 方法研究了电场作用下板之间的水或水模的蒸发行为,发现在外加电场作用下,会增强水的蒸发.Okuno 等[24]利用分子动力学方法研究了不同方向的匀强电场对水滴蒸发过程的影响,结果表明,电场方向垂直于气液界面时,电场会促进水滴蒸发;电场方向平行于气液界面时,电场会抑制水滴的蒸发[25].Hens 等[26]利用分子动力学方法模拟了水平和垂直方向的电场对固体壁面上液滴的扩散和蒸发行为,结论与文献[24] 类似.Wang 等[27]利用MD 方法研究了电场对悬浮三维微纳米水滴以及盐溶液液滴蒸发过程特性的影响,主要探讨了在电场下不同的盐溶液以及盐溶液浓度对蒸发特性的影响.

综上所述,学者研究的焦点主要集中在匀强电场的影响上,而匀强电场对液滴或水膜蒸发的影响与电场的方向有关.交变电场的方向总是随时间的变化而不断变化的,目前交变电场对蒸发影响的研究相对较少.张立等[28]通过实验研究了均匀和不均匀电场对液滴蒸发的影响,结果表明,当电场强度足够大时,不均匀电场相比于匀强电场更显著地促进液滴的蒸发.Zhang 等[29]利用MD 方法研究了方波电场下水膜和盐溶液膜的蒸发特性,表明方波电场比匀强电场更能够促进水薄膜的蒸发,且频率越大,蒸发越快,但在交变电场作用下液滴蒸发过程中参数的变化和机理尚未详细讨论.本文利用分子动力学模拟方法,研究正弦交变电场对三维悬浮水滴在超临界氮气环境下蒸发过程的影响,系统地讨论了电场幅值和频率对液滴蒸发特性的影响,包括液滴的蒸发寿命和速率、瞬时蒸发速率、液滴温度变化等方面.并从水分子的排列结构、能量变化和所加交变电场的关系等角度解释了交变电场作用下液滴蒸发特性的不同.

1 模型建立及验证

1.1 模型建立

建立了如图1 所示的液滴蒸发模型.其中,模拟区域为正六面体,其边长l=24 nm,分为气态和液态两个区域.其中气态的成分为27 000 个氮气分子,均匀地分布在整个模拟区域内,位于模拟区域中心位置的球形液滴由8000 个水分子组成.

图1 液滴蒸发模型Fig.1 Droplet evaporation model

在模拟过程中,蒸发模型采用的是周期性边界条件,即原子从一个边界出去,会从另一个边界进来,保证了蒸发过程中原子数的恒定,同时克服了模拟体系较小的问题.

在上述建立的模型中,水分子选取了刚性四点(TIP4P/2005) 水模型[30],由3 个固定点电荷和一个Lennard-Jones 中心组成,如图2 所示.其中,m点为引入的无质量电荷点.这里采用简单的L-J 势能模型描述水分子和氮气分子中原子间的相互作用力,如式(1)所示

图2 TIP4P/2005 水模型示意图Fig.2 TIP4P/2005 water model diagram

式中,ε 为势阱深度;σ 为势能为零的点距原子核心的距离;r为原子间的距离,rc为截断距离;exp(-κr)为阻尼因子;C为库伦常数,C=9 GN·m2/c2;εr为相对介电常数.其中,水分子和氮分子的L-J 势能参数设置如表1 所示[30].表中1 kcal=4.18 kJ.

表1 水和氮气分子的L-J 势能参数Table 1 L-J potential energy parameters of water and nitrogen molecules

对于不同原子间的相互作用的L-J 势能参数,可按以下公式选取

此外,水、氮气的成键作用力或键角由谐振函数控制

式中,Kl为拉伸强度,l为瞬时键长,l0为平衡键长;Kθ为旋转强度,θ 为瞬时键角,θ0为平衡键角.表2 中给出了水和氮气分子的结构参数.

表2 水和氮气分子的结构参数Table 2 Structural parameters of water and nitrogen molecules

1.2 模型验证

为了验证上述建立的分子模型以及蒸发模型的正确性,下面分别计算了水的部分物性参数和液滴在亚临界环境下的蒸发特性,并分别与理论数据或文献中的结果进行对比.

首先模拟计算了水在不同状态下的热导率.在本文中,热导率是通过Green-Kubo 公式计算求得,如式(5)和式(6),表明输运系数等于自关联函数对时间的积分[31-32],其中,热导率是热流J的自关联函数的积分[33]

式中,V表示模拟区域体积,KB表示玻尔兹曼常数,T为模拟体系的温度大小.

利用式(5) 和式(6),得到了不同状态下水的热导率模拟结果,如图3 所示.可以看出,无论是在液态、气态或是超临界状态下,模拟计算得到的热导率和黏性系数结果与IAPWS 理论数据均符合良好,分子模型可以正确反映水的宏观性质,从而验证了分子模型的正确性,同时也验证了微观小尺度的模拟可以反映物质的宏观现象规律的结论.

图3 水的热导率模拟结果Fig.3 Simulation results of thermal conductivity of water

接着,为了验证蒸发模型的正确性,选取了与文献[27] 中一致的参数,模拟了水滴在亚临界氮气环境下和匀强电场作用下的蒸发过程.首先使用NVT系综使环境和液滴达到平衡状态,即使氮气分子均匀地充满整个空间,液滴在表面张力作用下变成球形,并使液滴温度保持为298 K.接着提高环境气体的温度如1000 K,然后去掉NVT 系综并利用NVE系综,使系统中能量自由交换,能量从高温的氮气流向低温的水分子[34-35],从而使液滴的温度升高,不断蒸发.

在模拟过程中,匀强电场的施加是通过对体系中每一个带电的原子施加一个电场力F=qE来实现的,其中F的正负表示电场力的方向.需要注意的是,电场的施加是在体系达到平衡后,即在NVE 系综下施加.

对模拟得到的结果与文献中结果进行对比,如图4 所示.可以看出,模拟得到液滴分子数变化曲线与文献中的数据几乎重合,表明蒸发模型以及蒸发流程是正确的.

2 结果及讨论

2.1 电场幅值的影响

选定电场频率f=5 GHz,分别设置电场幅值E0=0.4,1.2,4,8 以及12 V/nm,其电场形式如式(7)所示,其中电场在t=80 ps 时刻施加到模拟区域中.

接着,在初始环境温度T=1000 K,初始环境压力P=26.961 MPa 工况下,耦合上述电场,模拟了水滴在氮气环境中的蒸发过程,得到了不同电场幅值下液滴内的水分子数随时间的变化关系,如图5 所示.对图5 中的数据进行处理,得到了液滴的蒸发寿命,并与匀强电场下液滴蒸发的寿命进行对比,如图6 所示.可以看出,在交变电场下,随着电场幅值的增大,液滴蒸发的寿命是不断下降的,且相比于匀强电场,相同幅值的交变电场更能够促进液滴的蒸发,这与文献[14]中的结果在定性上保持一致.

图5 不同电场幅值下液滴内分子数随时间的变化Fig.5 Change of molecular number in droplet with different electric field amplitudes

图6 交变电场和匀强电场下的液滴蒸发寿命对比Fig.6 Comparison of droplet evaporation lifetime under alternating and uniform electric field

相比于匀强电场,交变电场的方向和大小总是随时间周期性变化,这会导致液滴内水分子所受到的电场力也是周期性变化的.图7 给出了液滴的瞬时蒸发速率随时间的变化.在本文中,液滴瞬时蒸发速率表示在一个步长的时间内,液滴内水分子数的变化,用来表示液滴的瞬时质量蒸发速率.

图7 不同电场形式下液滴的瞬时蒸发速率Fig.7 Instantaneous evaporation rate of droplet under different electric field

图7 中分别给出了无电场、匀强电场以及交变电场作用下的液滴瞬时蒸发速率(电场强度或幅值均为4 V/nm).可以看出,当电场为匀强电场或不存在电场时,瞬时蒸发速率曲线表现出无规则的振荡特性,主要表现为小振荡幅值的无规则的振荡,这是由分子动力学本身的热力学涨落和统计误差造成的波动;而当电场为交变电场时(红色曲线),曲线振荡变得更加有规律,主要表现为更显著的振荡幅值以及两个振荡峰值之间的间距几乎保持一致.为了进一步探究交变电场下液滴瞬时蒸发速率的振荡特性与所加电场的关系,统计了蒸发速率的峰值对应的时刻,并与电场的节点时刻进行对比,如图8 所示.可以看出,液滴瞬时蒸发速率的振荡峰值总是发生在电场的节点所对应的时刻,这表明液滴的瞬时蒸发速率的振荡频率是所加电场频率的二倍,且总是在每一个电场的节点处,液滴瞬时蒸发速率达到最大.

图8 蒸发速率振荡峰值与交变电场的对应关系Fig.8 Relationship between evaporation rate frequency and electric field frequency

为了解释交变电场对液滴蒸发寿命和瞬时蒸发速率的影响,下面主要从液滴的能量和水分子的排列结构两个方面进行分析.不同幅值的交变电场作用下蒸发过程中液滴的温度随时间的变化如图9 所示.当不存在电场或电场强度较小,A=0.12 时,液滴温度的变化曲线几乎重合,且液滴温度会达到一个稳定温度.随着电场幅值的继续增大,液滴温度变化曲线表现出明显的振荡特性,通过分析,可以得到其振荡频率仍然是是所加电场频率的二倍,而可以明显看出振荡幅值随着电场幅值的增大而不断增大.同时,液滴温度随时间是不断增大的,没有表现出与无电场下液滴温度相对稳定类似的性质,且电场幅值越大,液滴温度上升地越快,在同一时刻液滴温度越大.这就导致了电场幅值越大,液滴蒸发的寿命越小.

图9 液滴温度随时间的变化Fig.9 Change of droplet temperature with time

接着,从水分子排列结构方面分析交变电场对液滴蒸发的影响.水分子间O-H 的归一化径向分布函数的定义为:系统中一个氧原子周围的氢原子的局域密度与系统中氢原子的平均密度之比(不包括与该氧原子组成水分子的氢原子),可表示为

式中,dM为距离一个氧原子r+dr内的氢原子数,ρ为氢原子的平均数密度.图10 给出了不同时刻下体系中水分子间O-H 的归一化径向分布函数.

图10 E=8 V/nm 水分子的O-H 归一化径向分布函数Fig.10 O-H normalized radial distribution function of water molecules at E=8 V/nm

图10 中给出了在电场一个周期内的O-H 的归一化径向分布函数变化,可以看出,在不同的时刻,第一峰值的位置几乎不发生变化,而第一峰值的大小并不是单调变化的,而是随着时间时而下降时而升高.有学者认为0.18 nm 是形成氢键的氧原子和氢原子间最可能的距离[31],而径向分布函数中第一峰值的位置为0.173 5 nm,在一定程度上可认为是形成氢键的最可能的长度.第一峰值的大小表示一个水分子最可能的归一化氢键密度.第一峰值位置和大小可以表示一个水分子受到的作用力大小.

为了更清楚地研究水分子排列结构在交变电场下的变化规律,下面给出了电场幅值E=8 V/nm 下的第一峰值大小随时间的变化关系,如图11 所示.可以看出,第一峰值的大小也会产生与所加正弦电场类似的振荡,且第一峰值振荡的峰值发生在所加电场对应的结点时刻附近.这表明体系中水分子间的相互作用力也是周期性变化的,同时观察图11 可以看出,在施加电场后,第一峰值的振荡波谷是与电场的节点相对应的,表明在该时刻水分子间的相互作用力相对最小,使水分子更容易从液滴内扩散到环境中,这导致了液滴的瞬时蒸发速率在电场的节点时刻会产生一个峰值.

图11 E=8 V/nm 下径向分布函数第一峰值随时间的变化Fig.11 Change of the first peak value with time at E=8 V/nm

2.2 电场频率的影响

本节中研究了电场频率对液滴蒸发过程的影响.选定电场幅值E=4 V/nm,通过改变式(7)中的t0大小,来改变电场的频率大小,表3 给出了部分选定的t0以及对应的周期和频率大小.

表3 选择的部分交变电场频率Table 3 Frequency of alternating electric field

图12 给出了不同频率的交变电场作用下,液滴蒸发寿命的变化规律.可以看出,蒸发寿命随着频率的增加并不是单调变化的,在f=5 GHz 时,液滴的蒸发寿命达到一个极大值.当f<5 GHz 时,随着频率的增大,液滴蒸发寿命不断升高,但通过对数据进行分析,可以得出在f=3.13～5 GHz 内,液滴蒸发寿命的变化较为不明显,这与文献中频率对液滴蒸发的影响较小的结论保持一致.当频率继续增大,f>5 GHz 时,液滴蒸发寿命开始逐渐下降,且随着频率的大幅增大,液滴的蒸发寿命也会急剧下降,这与文献[29] 中的结论一致.但作者仅对0.5 GHz,100 GHz,200 GHz,500 GHz 的方波对液膜蒸发的影响进行了模拟,缺少频率较小时的数据.

图12 液滴蒸发寿命和频率的关系Fig.12 Relationship between droplet evaporation life and frequency

为了进一步研究交变电场对液滴蒸发过程的影响,图13 给出了液滴的瞬时蒸发速率.

图13 f <6.25 GHz 电场下的液滴瞬时蒸发速率Fig.13 Droplet evaporation rate at f <6.25 GHz

可以看出,在不同频率的电场作用下液滴的瞬时蒸发速率都存在着明显的规律,如显著的振荡幅值、几乎一致的振荡时间间隔等.在f=5 GHz 时,液滴的瞬时蒸发速率会产生所加电场频率二倍频率的振荡特性,下面给出了不同频率下液滴的瞬时蒸发速率振荡峰值对应的时刻与电场波节时刻的关系,如图14 所示,其中空心符号表示液滴瞬时蒸发速率的峰值对应的时刻,实心符号表示电场节点的时刻.可以看出,随着电场频率的增大,蒸发速率峰值的时刻与电场波节对应的时刻符合得更好,当电场较小时,二者存在一些偏差,而当频率较大时,蒸发速率峰值出现的频率明显下降,在某些电场波节时刻,蒸发速率并不产生峰值,表示不再满足二倍频率的关系.这是由于当频率较大时,所加电场的变化较快,电场力与粒子本身的热运动相互耗散,导致液滴瞬时蒸发速率的振荡特性逐渐消失.

图14 不同频率下蒸发速率峰值和电场波节对应关系Fig.14 Relationship between evaporation rate frequency and electric field frequency

为了研究上述液滴蒸发的寿命以及瞬时蒸发速率变化趋势的产生原因,同样从液滴的能量变化和水分子的排列结构变化两个方面进行分析.由2.1 节中内容可知,交变电场对水分子结构的影响主要体现在第一壳层水分子的密度,即径向分布函数的第一峰值大小,而对第一峰值的位置几乎没有影响.

因此,图15 和图16 给出了不同频率的电场作用下归一化O-H 径向分布函数的第一峰值大小随时间的变化关系.可以看到,当f<5 GHz 时,第一峰值的振荡峰值和波谷几乎不发生变化,而频率继续增大时,峰值同样不发生变化,而波谷显著下降,如图17所示,表明分子间作用力显著下降,从而造成了当频率较大时,液滴的蒸发寿命不断下降.

图15 f ≥6.25 GHz 下O-H 径向分布函数第一峰值大小随时间的变化Fig.15 The first peak value of O-H radial distribution function at f ≥6.25 GHz

图16 f ≤5 GHz 下O-H 径向分布函数第一峰值大小随时间的变化Fig.16 The first peak value of O-H radial distribution function at f ≤5 GHz

图17 第一峰值振荡的波谷值与频率的关系Fig.17 Relationship between the trough value and the frequency of the first peak oscillation

接着对图15 和图16 中的径向分布函数的第一峰值的振荡频率进行了分析,如图18 所示.可以看出,第一峰值的振荡频率和交变电场的频率有着良好的线性关系,通过拟合处理可以得到斜率k=2.001 61,表明第一峰值的振荡频率为所加交变电场频率的二倍.

图18 第一峰值振荡频率与所加电场频率的关系Fig.18 Relationship between the oscillation frequency and the frequency of the electric field

上面的结论可以解释在高频电场下(f>6.25 GHz)液滴蒸发寿命下降的原因,以及液滴的瞬时蒸发速率产生所加电场频率二倍的振荡特性的原因.但无法解释在f=5 GHz 时,液滴的蒸发寿命产生极值点.在f<6.25 GHz 时径向分布函数的第一峰值的振荡幅值几乎保持一致,对水分子结构的影响也应是一致的,但考虑到水分子排列结构变化频率的不同,造成在一个电场周期内的一定时间内,电场对液滴做的功是不一致的,这会导致水分子能量变化的不一致,从而导致液滴蒸发寿命产生差异.

接着,给出了不同频率的交变电场下液滴的动能随时间的变化,如图19 所示.可以看到,液滴的动能同样会产生与所加电场频率的二倍频率的振荡.当频率较大f>6.25 GHz 时,随着频率的增大,液滴动能的振荡幅值也在逐渐增大,同时液滴所达到的稳定的动能也更大,这也会导致液滴蒸发的加快;而当频率较小(f<5 GHz)时,动能的振荡幅值和液滴最终达到的平均温度几乎保持一致,且f=5 GHz时所达到的稳定动能要稍微小一些,f=6.25 GHz所达到的稳定动能要稍微大一些,此外,在初始阶段,f=5 GHz 比f<5 GHz 以及f=6.25 GHz 时的动能变化更慢,且随着频率升高,动能在初始阶段的变化更慢,而在f=6.25 GHz 时,动能在初始阶段变化又加快,这些都会导致液滴的蒸发寿命在f=5 GHz 时达到极大值.