应用分子动力学方法研究纳米流道不同固–液间相互作用势能指数对液体流动特性的影响。研究结果显示,纳米尺度下流道出口高度及固液间相互作用势能不同会引起流体压力的改变,并出现射流现象。出口端距离较小时,流道内液体流动规律不再服从NS方程——按照分子动力学模拟与按NS或Renolds方程求解得到的压力分布结果存在较大的差异,此时射流现象比较明显,射流速度随着固液间势能指数的增强呈线性增大趋势。随着出口端高度的增加,纳米尺度效应减弱,压力分布逐渐趋近宏观流动下的压力分布曲线,固液间势能指数对于射流速度的大小不再产生明显的影响。 Molecular dynamics method is applied to study the influence of potential interaction strength be-tween the liquid and the solid on the properties of fluid film in nanochannels. The results indicate: the pressure distribution of fluid film is changed in Nanochannels with the change of the height of channel in the outlet end and the liquid-solid potential interaction strength. At the same time, the jet phenomenon can occur in the outlet end. The difference of the pressure distribution between the results obtained by molecular dynamics simulation and that by NS or Renolds equation is much bigger. At this point, the jet phenomenon is more obvious. The jet velocity increases linearly as the liquid-solid potential rises. With the increasing height of the outlet, the nanoscale effect becomes weaker and weaker. The pressure profile gradually approaches to that of the macro-flow. The liquid-solid potential has no significant effect on the jet velocity.
贾 妍,薛晔
西安工程大学,陕西 西安
Email: jiayan05@126.com
收稿日期:2015年4月28日;录用日期:2015年5月12日;发布日期:2015年5月18日
应用分子动力学方法研究纳米流道不同固–液间相互作用势能指数对液体流动特性的影响。研究结果显示,纳米尺度下流道出口高度及固液间相互作用势能不同会引起流体压力的改变,并出现射流现象。出口端距离较小时,流道内液体流动规律不再服从NS方程——按照分子动力学模拟与按NS或Renolds方程求解得到的压力分布结果存在较大的差异,此时射流现象比较明显,射流速度随着固液间势能指数的增强呈线性增大趋势。随着出口端高度的增加,纳米尺度效应减弱,压力分布逐渐趋近宏观流动下的压力分布曲线,固液间势能指数对于射流速度的大小不再产生明显的影响。
关键词 :分子动力学,势能指数,液体薄膜,射流
宏观状态下当运动表面正常运转时,两个固体壁面间液体的流动特性可以通过NS方程或Reynolds方程研究。但随着微纳技术在众多科学领域的广泛应用,纳米尺度下流体的流动特性是否依然符合经典的流体动力学理论已经受到人们越来越广泛的重视。当物体尺寸达到纳米级以后其性质会发生重大的改变,此时用传统的连续介质理论研究纳米尺寸下流体的流动特性可能已不再有效,而从统计物理学中衍生出来的分子动力学模拟,由于不需要引入经典流体动力学的连续介质假设,以微观粒子间的相互作用势为基础可以研究尺寸极小的系统,因此在探索纳米尺寸下物质特性方面有着独到的优势。
Shin-Pon Ju [
然而,目前无论对于宏观领域内的研究还是微观领域内液体薄膜的研究都是针对两固体壁面平行条件下液体薄膜行程以后,而对于存在于一般机械系统中的两接触壁面分离过程中流体薄膜在纳米尺度下流动演化特性的研究很少涉及。本文采用分子动力学的方法主要研究两个固体壁面分离过程中形成的非平行纳米流道内不同固液间相互作用势指数对液体薄膜特性的影响。
模拟系统由固体壁面和壁面间的液体分子构成,如图1所示。为了研究方便模拟过程中保持上壁面静止不动,让下壁面沿x方向产生相对运动。
模拟过程中边界条件和分子插入原则按照ITSUO [
图1. 模拟系统示意图。○表示流体分子,■表示固体壁面分子
纳米管内液体的流动特性和水的射流现象时选取的边界条件和分子插入原则。仅z方向采用周期性边界条件。插入分子的速度按照初始温度根据麦克斯韦–玻尔兹曼分布随机确定;插入分子的位置在采取随机确定的同时还要遵循以下的原则:
首先,当初始条件下模拟区域内的分子数确定时,靠近入口端处单层分子晶格中所包含的分子的数目即为初始条件下尝试插入的分子数目。并且尝试插入的数目随着模拟区域中分子数目和壁面运动速度的变化而变化。当模拟区域的分子数增加时,认为尝试插入的分子数减小;当模拟区域内的分子数小于初始条件对应的分子数时,尝试插入的分子数增加。每时间步尝试插入的分子数遵循以下的方式:
其中:Nat表示每时间步尝试插入的分子数;v0表示下壁面运动速度;Nadm表示初始状态外部区间每时间步长尝试插入分子数;Nmax表示模拟区间最大容纳分子数;Nt表示模拟区间瞬时分子数;Nint表示初始状态模拟区间内能容纳分子数。
模拟区域最大容纳的分子数Nmax根据YAO [
其次,分子的尝试插入过程中还应满足MIHALY [
为了研究模拟区域内液体沿x方向不同层面特性的变化,将模拟区域沿x方向划分成若干个小的子区域。由于模拟中两固体壁面处于非平行状态,为了方便各小区域内压强的比较,小区域的划分遵循等体积原则。
其中:Vl表示小层面的体积;dx,dy,dz表示模拟区域沿x,y和z方向尺寸;h表示两固体壁面间出口端的距离;L表示划分小区域的个数。
各个小层面压力计算公式 [
压力计算公式中,第一项表示分子动能对压力作用效果,第二项表示粒子间相互作用势能对压力的作用效果。
其中:α,β表示分别表示x,y和z方向;viα表示第i个粒子沿α方向的速度分量;viβ表示第i个粒子沿β方向的速度分量;mi表示第i个粒子的质量;rijα表示第i与第j个粒子间距离沿α方向分量;fijβ表示第i与第j个粒子间相互作用力沿β方向分量。
模拟系统中的液体采用液态氩,液态氩分子间的相互作用势采用Lennard-Jones (LJ) [
其中:rij表示第i个分子与第j个分子之间的距离;σ表示分子特征长度,也是分子的直径;ε表示分子间相互作用强度的能量特征值。
整个模拟过程采用无量纲参数,具体约化单位处理见表1。
为反映液体与固体之间不同作用势能对液体润滑膜的影响效果,分子动力学模拟选取的液体和固体壁面之间的LJ相互作用势为 [
其中:usl表示液体分子和固体分子间的LJ势能;εsl表示液体分子和固体分子之间的作用势能;σsl表示液体分子和固体分子之间的距离参数;c表示用来调整液体分子和固体壁面之间的相互作用效果,即浸润性,其意义是指液体分子与固体壁面分子之间的作用力相对于距离相等的液体分子和固体分子作用力的倍数,c越大则表示液体分子和壁面固体分子之间的作用势能越强,c = 1则表示液体分子和固体分子之间的势能作用。本文中为了方便计算选取固液间分子的作用势能和距离参数等于液体分子的作用势能和距离参数。
物理量 | 表达式 | 数值 | 国际单位 | |||
---|---|---|---|---|---|---|
长度 | σ | 3.4 × 10−10 | m | |||
能量 | ε | 1.65 × 10−21 | J | |||
质量 | m | 6.69 × 10−26 | Kg | |||
时间 | 2.16 × 10−12 | s | ||||
玻尔兹曼常数 | kB | 1.381 × 10−23 | J/K | |||
温度 | ε/kB | 119 | K | |||
速度 | σ/τ | 157 | m/s | |||
加速度 | σ/τ2 | 7.27 × 10−13 | m/s2 | |||
压力 | ε/σ3 | 4.22 × 107 | pa | |||
数密度 | 1/σ3 | 2.54 × 1028 | 1/m3 | |||
表1. 分子动力学模拟中物理量的约化单位
对液体进行分子动力学模拟研究时液体分子的初始位置不会改变模拟系统最终的特性,原则上对于液体任意的初始分布都是合理的 [
模拟过程中纳米流道入口端的宽度保持不变均为hy = 9.35,系统在x和z方向的计算域尺寸分别为34.8和5.1。h为纳米流道出口端不同分离阶段两个固体壁面间的距离。
当粒子初始位置,速度和粒子上的作用力确定后,根据牛顿运动定律,液体分子的运动方程为:
为研究液体与固体壁面间不同势能作用指数对模拟区域内液体流动特性的影响。模拟中保持下壁面以稳定的速度v0 = 1.0运行,液体与固体壁面间势能作用的强弱通过势能作用指数c进行调节。本文采用的几种典型固液间势能作用指数分别为0.2、0.5、0.7、0.9、1.0、2.0、3.0和4.0。
在保持纳米流道入口端距离不变的情况下改变出口端距离,获得出口端不同状态下所对应的沿x方向液体压力随固液间相互作用势能指数的变化曲线,如图2所示。
当出口端的距离较小时,模拟获得的液体的压力分布曲线与宏观状态下通过NS方程求得的液体压力分布曲线存在较大差异,随着纳米流道出口端距离的增加,模拟区域内液体压力的最大值不断下降,并随固体壁面间的距离的增加越来越偏离出口端,同时液体压力的变化曲线越来越平缓,当h > 5.1时,模拟区域内液体压力分布曲线与宏观状态下的压力分布曲线基本相符合。
纳米流道内,固液间相互作用势能指数较弱时,其对液体压力的影响较大,当c < 1.0时,随着固液间相互作用势能指数的增加,模拟区域内液体的压力不断升高;当c > 1.0时,固液间势能作用指数较强时模拟区域内液体压力分布基本不受固液间相互作用势能指数增加的影响。引起这种现象的主要原因是由于当c < 1.0时,固液间相互作用势能指数较弱,液体在固体表面表现为疏水性,固体壁面运动运动过程中壁面的固体分子始终对液体分子的运动产生影响;当c > 1.0时,固液间相互作用势能指数较强,液体在固体表面表现为亲水性,固体壁面运动过程中由于固体分子与液体分子的强烈作用会使固体表面附着较多的液体分子,在分子相互间力的作用范围内实际表现为液体分子之间的相互作用,此时固液间相互作用势能指数的增加对液体压力的影响基本保持不变。
纳米流道出口端不同状态下,模拟区域内液体最大压力随固液间势能指数变化曲线,如图3所示。
当h = 1.7时,随着固液间势能指数的增加,液体压力的最大值不断增加;当h > 3.4,c < 1.0时,对于固液间相互作用势能指数较弱的情况,随着固液间势能作用指数的增强不断增加;但当c > 1.0时,固液间相互作用势能指数改变对液体最大压力的影响基本可以忽然不计。
纳米流道出口端距离改变时,射流速度随固液间相互作用势能指数的变化曲线,如图4所示。
当h = 1.7时,由于出口端两固体壁面间的距离较小,固液间势能作用强弱对射流速度的影响较为明显,随着势能指数从0.2增加到4.0,出口端的射流速度呈线性的增加。
[
[
[
[
[
[
[
[
[
[
[