水导热系数的分子动力学模拟研究.pdfVIP

中周工程热物理学会 工程热力学1』能源利用 5 学术会议论文 编4，：06111 水导热系数的分子动力学模拟 刘娟芳，曾丹苓，刘朝 重庆大学动力工程学院(400044) E-mail：vacoljf@yahoo．com．ca，Tel：023 摘要，奉文选用平衡分子动力学模拟方法-应用不『一的势能模型对强极性分子水的导热系数进行 了模拟计算。结果表明，用分子动力学模拟方法计算水的导热系数时·模拟结果对模型的选取极为 敏感．在目前应用较为广泛的几种模型中咀TlP4P模型较为适用，所得的导热系数与实验值较为接 近． 关键词；分子动力学模拟方法，微观热流，Green-Kubo公式，导热系数 1引言 物质的物性参数一般通过实验测定、模型计算和计算机分子模拟三种方法得到， 导热系数是物质的一种很重要的输运系数。近年来，随着计算机技术的迅速发展，分 子模拟已成为在微观尺度下研究流体性质的一种强有力的工具，越来越被广泛地应用 于各个领域。 分子动力学模拟方法不仅能获得物质的平衡性质，也能得到其输运特性。它通过 求解体系中粒子的运动方程获得粒子的运动速度和运动轨迹，并经统计平均求得物质 的宏观特性。本文的目的旨在以．T：程常用的』=质水为例，用分子动力学模拟方法计算 其导热系数．但文中所提供的方法也可适用其它类似的流体· 2势能模型及其势能函数 分子动力学模拟结果的准确性主要取决于是否选取了能够准确地描述分子微观结 构的势能模型。对于强极性水分子．由于其结构和运动较为复杂-许多学者提出了不 模型，还有晟近提出的TIP5P、SPC柔性模型等·这些模型各有优劣·本文选取较为 常用的SPC、SPCE和TIP4P三种模型计算水的扩散系数-并对模拟结果进行了比较 和评价。 SPC、SPCE和TIP4P三种模型均为刚体模型，其势能函数形式如下： 嘛最，置)=；菩等慨陋]12-刚 下标WW表示水分子之间的相互作用，OO表示水分子中氧原子间的相互作_}jfj；％为分 W W 子j和分子j质心之间的距离；Ql和Q2分别为描述分子i和j的方位角；矿为分子i 的口作用鹰所带的电量：以为模型中一个水分子的作用座数t对于SPC和SPCE模型 ^，。=3，TIP4P模型Na=4。，尹为分子i的口作用座和分子j的∥作用座之间的距离： 占和盯分别为LJ势能的能量和尺寸参数。三种模型的位能参数如表1所示。 表ITIP4P，SPC和SPCE三种模型的位能参数 Z ’钿 ，．OM 0ID Fo／k 模型 HOH ％(e)钆(e) r) (局 (胸 (南 (K) TIP4P104．5209575O．15 3．154 78．020．1．04 0．52 SPC109．471．0 3．165678．182-o．82 041 SPCE109．47 3．166 78．182-08476 1．O 0．4238 3模拟方法 模拟体系由256个水分子组成。初始时刻分子按FCC构型排布。水分子起始平动 速度按Maxwell分布随机取值，初始转动速度均取为零，初始方位为随机取向。然后 对速度进行标定，使得体系的总动量为零。为消除界面效应，模拟盒的三个方向均采 用周期性边界条件． 式(”中长程静电作用项衰减得较为缓慢，在模拟中采用反应场法对其进行处 理。为此势能函数变为如下形式t 式中g，为环境介电常数．计算中取为无穷。式(1)中的短程u作用，采用球形切 断法，截断半径为模拟盒边长的一半，截断半径以外的分子间