界面微观特性分子动力学模拟研究.pdf

Study on Microscope Characteristics of Interface by Molecular Dynamics Simulation A Thesis Submitted to Chongqing University in Partial Fulfillment of the Requirement for the Degree of Doctor of Engineering by Cai Zhiyong Supervisor: Prof. Zeng Danling Major: Power Engineering EngineeringThermophysics College of Power Engineering of Chongqing University, Chongqing, China April 2008 中文摘要 摘 要 本 论 文的主要目的是利用分子动力学模拟方法（ Molecular Dynamic Simulation，简称MD ）从微观上系统地研究相变及界面现象。 相变是工程中常见的物理现象，相界面是不同相间的边界面。由于界面非常 薄，在许多场合下仅有几个分子直径厚，因而常被忽略而将其想像为一个没有厚 度的几何面。但正是这一薄薄的界面的存在，对界面区的热力学特性及其附近的 流动和输运、化学过程产生了重要的影响。近几年来，国内外对界面性质及相间 输运现象的研究非常活跃，界面现象的探索长期以来都是科学研究的一个重点领 域。 本论文借助分子动力学模拟方法，从微尺度上全面系统地研究界面层的物理 特性和流体宏观特性与微观特性之间的关系，找出了宏观现象对应的微观图象， 使我们能更好地把握现象的本质，以丰富理论、指导实践。 本论文所作的主要工作如下： ①详细阐述了分子动力学模拟的原理和关键技术，构造了分子动力学模拟方 法的计算程序。以LJ 流体氩、甲烷和强极性物质水为例，计算了流体的饱和特性， 确定了饱和压力与温度的相互关系。并通过与实际气体状态方程的计算结果及相 应的实验值进行了比较，证明了分子动力学模拟方法在模拟气体物性上是成功的， 并对现有实际气体状态的正确性利用 MD 方法进行了对比、验证。另外，运用曾 丹苓教授提出的用分数布朗函数作为描写实际流体分子无规则运动的概率密度函 数的理论，获取了分子动力学模拟中大量的分子运动的微观信息，定量描述了实 际气体偏离理想气体的程度。 ②用 MD 方法研究了气-液相平衡。发现此时除气、液两相外，还存在一个各 向异性的界面相。本文仔细地研究了界面相中流体粒子的数密度分布、温度分布 及法向切向应力分布，计算了表面张力。在气-液界面特性的分子动力学模拟方法 研究中，由于界面层存在各向异性，本文特别分析了势能截断半径的选取对均匀 相中饱和气体和饱和液体的密度以及非均匀相中表面张力的影响，提出了一种能 够同时满足均匀相中饱和气体、饱和液体密度和非均匀相中表面张力计算要求的 截断半径选取方法-变截断半径算法。此算法克服了传统模拟方法中统一的截断半 径不可能同时满足均匀相中和非均匀相中相关参数的计算精度要求的缺点，提高 了计算精度，节约了计算时间；并将该算法运用在分子动力学模拟中，进一步分 析了涨落现象对非均匀相的影响。 ③在相变过程界面现象的研究中，本论文从平界面的相平衡发展到球界面的相 I 重庆大学博士学位论文 平衡研究，从微观尺度上分析研究了纳米液滴蒸发、冷却凝固的微观过程。以径 向分布函数理论为基础，更为详细地探究了纳米液滴在蒸发和冷却凝固过程中界 面现象及过程的微观机制。并从微尺度液滴的模拟中获取了相应的宏观性质，这 将为探索宏观系统尚未知或未确定的物理性质打下基础。 ④本论文研究了分子动力学模拟中固壁对