采用双原子分子旋转计算及动量控制方法的空气分子动力学计算.pdfVIP

第27卷第3期 V01．27 分子科学学报 No．3 2011年6月 OFMo】皿CI『I．AR J叫RNAL SCⅡNCE 2011 June [文章编号]1000-9035[2011)03．0203-05 采用双原子分子旋转计算及动量控制方法 的空气分子动力学计算 张志杰，张文飞。，张钧惠，孙虹宇，李小龙 (燕山大学建筑工程与力学学院，河北秦皇岛066004) [摘要]用分子动力学方法，采用Lennard．Jones势函数对空气的物理性质进行了模拟计算． 在模拟过程中，针对气体的特点，应用了对系统动能和动量的控制，以及双原子分子的旋转计 算．模拟了不同高度下空气压强的变化，并计算了在标准大气压下空气温度和密度不同时黏度 的变化．结果表明，用分子动力学方法模拟得到的空气物理参数与实测数据吻合较好． [关键词]分子动力学模拟；Lennard—Jones势函数；压强；黏度 641 150·30 [中图分类号]0 [学科代码] [文献标识码]A 0引言 目前，关于气体性质的分子动力学模拟已经引起了广泛关注．刘艳杰等采用巨正则蒙特卡罗方法研 究了c02分子在H．sⅡ分子筛中的吸附u 对氮气的对势进行了研究匕J，得出了对势的第二维里系数、黏度、导热性、自扩散系数、热扩散因子等参 K一1600 数；Garruchet等应用分子动力学方法对在950K时0原子在Ni中的扩散进行了研究[3]3，发现 聚苯乙烯中的渗透和扩散进行了分子动力学模拟[4]，计算得到的扩散系数与实验数据和以前的模拟结 果吻合很好． 空气在不同情况下的物理性质计算，特别在超高压和超低压等条件下的空气的物理性质研究在很 多工程问题中都有重要的应用．因此，如果能用计算机模拟的方法来算出各种不同条件下空气的物理性 质是有意义的．本文运用分子动力学模拟的方法，计算气体的压强、黏度等参数．在计算的过程中根据气 体分子动力学的特点，提出了一些新的计算技巧．对计算结果与宏观理论研究和实验观测的结果进行了 比较，并分析总结了空气状态对空气物理性质变化的影响． 1分子动力学模拟技术 1．1势函数的选取 原子间的相互作用势是用分子动力学模拟材料性质的基础和关键．作为双原子分子的势函数，常用 该势函数计算简单，结果满足精度要求． 收稿日期：2010-07．14 联系人简介：张文飞(1954．)，男，教授，主要从事分子动力学研究 F_，mall：zwenfei@gmail．Ⅻ1 分子科学学报 第27卷 Lennard．Jones(12—6)势函数的表达式为 u(r)：4￡【(詈)12一(詈)6】 (1) 式中长度参数仃和能量参数￡要根据发生作用的原子的不同有所区别，当0原子和0原子作用时000= 0．3365nnl，eoo=1．597J；当N原子和N原子作用时ONN=0．363nln，￡NN=1．394Jmo；N原子和0原子 25 21．584 J．憎。 GNO=0．5×(ONN+000)=0．363nln，eNo2~／ENNgeOO 1．2数值计算方法 本文采用正则系综分子动力学算法【t03(N、呵算法)进行模拟．NVT算法是指对粒子数Ⅳ、体积y和