纳米流道内液体特性分子动力学研究.docVIP

纳米流道内液体特性分子动力学研究 　　摘　要：以非平衡分子动力学方法研究了纳米流道?纫禾宓牧鞫?特性，重点探讨上壁面以不同剪切速度运动，下壁面保持静止且在其附近的温度保持相对稳定时，两壁面间液体分子的运动特性及作用机理，模拟获得了不同剪切速度下流体的速度分布、密度分布及温度分布，并分析了速度滑移率、密度分布和最高温度与剪切速度之间的关系，研究结果表明，随着剪切速度的增加，两固体壁面间液体靠近运动平板处的温度呈线性增加，在距离运动壁面0.8倍分子直径范围之内存在着明显的滑移现象，而在这个区域之外，液体分子的速度基本呈线性变化，纳米流道内液体分子密度的非均匀分布和有序化排列结构不受剪切速度的影响，平均剪切应力随着剪切速度的增加而不断增加。 　　关键词：分子动力学；剪切速度；纳米流道；模拟 　　中图分类号：TH117　文献标志码：A 文章编号：0253-987X(2008)01-0009-04 　　 　　目前，应用分子动力学模拟方法研究纳米级流体的流动特性已取得了许多成果，Thompson和Robbins采用Lennard-Jones(L-J)液体(液态氩)模拟了2个平板间的剪切流动，着重研究了液体在壁面附近形成的有序结构及平板间流体的速度剖面。Thompson和Troian模拟了纳米流道内液体的无滑移流动和滑移流动，揭示了滑移长度和速度剪切率之间的关系，Bitsanis和Magda则采用L-J流体对固体边界附近的牛顿流体进行了详细的研究，得出了在剪切作用下液体的速度剖面、密度剖面和温度剖面，同时计算出液体的剪切应力和等效黏度，Cieplak等研究了流体在纳米流道内的密度分布现象，研究结果表明，滑移和壁面附近流体的组织结构有着密切的关系。 　　本文在模拟开始之前，首先假设在固体壁面吸附一层液体分子，流体中的液体分子与固体壁面的作用可以看作液体分子与附着在壁面上的液体分子之间的相互作用，然后采用非平衡分子动力学方法研究在不同剪切速度下，纳米流道内流体的运动特性。重点研究上壁面以不同剪切速度运动时，两壁面间液体分子的运动特性及作用机理，与以往处理方法不同之处在于，首先当上壁面开始以速度。v沿x方向运动时，如果有液体分子穿过上壁面，则此时液体分子速度在x方向的分量还应加上v，其次，以前关于流体在固体壁面之间的运动特性的研究，主要是建立在系统处于恒温或者恒压的基础之上的，但在实际应用中的有些情况下很难使整个系统的温度保持稳定，比如轴承转子的运动，就很难保证运动过程中润滑油膜整体温度稳定不变，因此，在模拟过程中，仅考虑静止壁面附近的液体分子温度处于稳定状态，即在实际应用中相当于给静止壁面增加了一个冷却作用，使靠近静止壁面的液体层的温度保持稳定。 　　 　　1 基本原理与模拟系统 　　 　　模拟系统由2个平行的固体壁面和壁面间的流体分子构成，如图1所示，液体采用液态氩，液态氩分子间的相互作用势采用L-J势，即 　　φ(rij)=4ε((σ/rij)12-(σ/rij)6)　(1)式中：rij为第i个分子与第j个分子之间的距离；σ为分子特征长度，也可以理解为分子的直径；ε为分子间相互作用强度的特征能量值。 　　设系统中液体分子与附着在固体壁面上液态氩分子之间的作用势也采用L-J势，但是此时势函数中的特征能量和特征长度分别为εwf和σwf，并且εwf/ε=1，σw/σ=1，计算中使用的模拟参数分别为：液态氩的特征长度σ=0.34mm，特征能量ε=1.65×10-21J，玻尔兹曼常数kB=1.381×10-23/K，系统温度r=132K，氩原子的质量m=6.96×10-24kg，数密度p=0.8，为了方便计算，整个模拟过程中数据采用归一化处理：时间r=(mσ/ε)1/2；速度v=σ/t；加速度a=σ/t2；温度T=ε/kg。 　　系统在xyz方向的尺寸为20.4σ×9.35σ×8.5σ，该系统包含1200个液态氩分子和240个附着在固体壁面上的液态氩分子，在初始时刻，流体分子按照面心立方品格结构排列，分子间的距离按照给定的密度值确定，根据已知温度确定分子的初始速度v0，方向则随机给出，当粒子初始位置、速度和粒子上的作用力确定后，根据牛顿运动定律，液体分子的运动方程为Fi=ma(i=1，2，…，N)，通过积分运动方程，可以得出任意时刻液体分子的速度和位置，积分过程选用蛙跳法，积分过程中的时间步长△h=0.005r，为节约计算时间，在计算分子间作用力时，取截断半径rc=2.5σ，即当2个粒子间的距离大于截断半径时。其相互作用可以忽略不计。 　　 　　在模拟过程中，为减少有限尺寸对结果的影响，系统在x和z方向采用周期性边界条件，即沿x和z方向模拟系统被其本身的映像系统所包围，而在y方向由于2个固体壁面的存在，因此采用Ashu