涨落现象对非均匀系分子模拟结果的影响.pdfVIP

中国工程热物理学会 T程热力学及能源利用 第卜一届学术会议 编l{}：051120 涨落现象对非均匀系 分子模拟结果的影响 蔡治勇 曾丹苓 刘娟芳 (重庆大学动力T程学院，400044) Tel：02365112184；E—mail：dlzeng@cqu．edu．Cll 摘要；本文提出了一种新的计算方法，该方法考虑了分子涨落现象对汽一液界面层的影 相中饱和气体密度、饱和液体密度以及非均匀相中表面张力变化，与传统的统计方法相 比，该方法计算非均匀相中表面张力的误差更小，更加节省模拟时间。 关键词：涨落现象， 非均匀系， 分子模拟 气液界面层是一个随时间变化而涨落起伏的介于气相区与液相区之间的过渡区域， 这个过渡区域并不是一个稳定的区域，在气液界面层内部密度、温度、压力张量、表面 张力等以及由此算得的热力学性质只是统计的结果。实验证明，在平衡态附近，气液界 面层内的微观结构、界面层厚度、表面张力等仍随着分子的涨落而波动。对于气、液体 相而言，由于其具有各向同性，涨落并不起重要作用。然而对于汽一液共存的过渡区而 言。由于本身的各向异性，且其包括的分子数较少，分子的位置、速度、能量的涨落将 不可忽略。这在计算表面张力是不能不考虑的。为此本文提出了一种新的表面张力的计 相平衡进行了分子动力学模拟，得出了均匀相中饱和气体密度、饱和液体密度以及非均 匀相中表面张力随模拟时间变化的分布情况，与传统的统计方法相比，非均匀相中表面 张力的计算结果更加接近实验值．且更加1，省模拟时间。 1分子动力学模拟方法 本文以C凡为例对气液平衡态进行了分子动 力学模拟．分子间的相互作用是基于L-J12—6模 |’ |。 }’ |’ 型，分子i和j相互间作州的势函数表示为： vapor liquid vapor 嘶)-4sk，。y2一-／勺6J (1) 其中参数a=3．7275x10…m，：E=2．057x10“J。模拟．| | | { 体系由分子质量m=2．655×10“kg的CH-分子构成， 采用直角坐标系Oxyz，分子在初始时刻为饱平¨液 幽1．模拟盒尺寸幽 136 相分子以边长为L的正方体FaG结构分布于模拟盒中部，其质心在原点(0，0，0)处。 模拟过程中，系统采用正则系综．即分子数N、容积V和温度T保持不变，并不断调整 体系质心使之处于原点。模拟盒在x、Y方向的K度Lx=Ly=L，为了保证足够的气相空间， 取在z方向的长度Lz=3L(如图1)。本模拟在三个方向上都采用周期性边界条件。 在计算机分子动力学模拟过程中，为确保模拟体系各物理量之间的数量级较为接 近，减小模拟过程中出现的误筹，将所用到的备物理量进行无量纲化。采用无量纲物理 量进行分子动力学模拟，分子的各物理量利用卜J12—6势幽数中的参数s和仃，以及体 系中分子质量m对模拟流体再物理量进行无量纲化处理如r： 长度r=L佑时间 t’=t(8／m一)”2温度T’=啪、／6表面张力 y+=Y62／￡ 能最U==U／g密度p’=p∥／m 压力P*=Po’／￡ 2计算公式 在进行分子动力学模拟时，质量为m的分子i的运动规律由联立求解各分子的运动 方程 弘≯b)叫争川：“ cz， 得出，其中F。是作川在分子i上的总作_}}；{力，0是i分子与j分子之间的距离，分子i 和j闻的相互作用力，?k，}是势函数∥m∥是竹，的导数，对L-J流体，如果取截断半径为 R，则作用力F“可表示为：