流体输运特性自相关函数积分的稳定性研究(一).pdfVIP

中国工程热物理学会 工程热力学与能源利用 学术会议论文 编号：091082 流体输运特性自相关函数积分的稳定性 研究(一) 陈俊 史琳 安青松 ( 清华大学热科学与动力工程教育部重点实验室 北 京100084 ) Tel: 010E-mail: rnxsl@mail.tsinghua.edu.cn 摘要：纳米制冷剂导热系数、粘度以及扩散系数等输运物性的分子动力学模拟研究，可以从微观 层面对纳米粒子的作用机理进行深层次的理解。本文首先通过采用Green-Kubo 自相关函数和平衡法对 简单流体Ar 的扩散系数、粘度、导热系数等输运参数进行了分子动力学模拟，发现利用平衡分子动力 学模拟方法求解流体输运特性时在一定时长后会出现波动很大的现象，进而对这一现象进行了深入的 探讨，发现这是利用分子动力学模拟方法对自相关函数积分求解流体输运热物性时本身就存在的问题。 本文为系列文章（一）和（二）。（一）主要明确提出稳定性的问题；（二）主要探讨出现波动的原因和 改进措施，本系列文章为纳米制冷剂输运物性较大时长分子动力学模拟奠定基础。 关键词：Green-Kubo 平衡分子动力学 导热系数 粘度 扩散系数 0 引言 纳米粒子由于其具有表面效应、量子尺寸效应、体积效应、宏观量子隧道效应等奇 特的性质，使得其在很多领域得到了应用。有关研究人员[1-3]已经开展了纳米粒子添加 剂在制冷系统中应用的初步研究，取得了很好的效果。导热系数体现了流体的传热特性， 而粘度和扩散系数则与流体的动力和质量传输，所以研究纳米制冷剂的这些输运性质是 纳米制冷剂应用的基础。目前很多学者对纳米制冷剂导热系数等热物性参数进行了深入 研究[4-6]，另外还有学者专门针对碳纳米管制冷剂进行了深入探讨[7]。但总体来说主要 是对宏观规律的掌握，微观的理论研究还需进一步进行。 目前国内外很多学者利用分子动力学对以导热系数为主简单流体的热物性进行了 模拟研究，但由于复杂性，将之应用到纳米流体还刚刚起步[7-11]。因此通过利用线性 响应理论和分子动力学模拟方法对纳米制冷剂的导热系数、粘度以及扩散系数等进行研 究，能从微观层面对纳米粒子的作用机理进行深层次的理解。本文首先对简单流体的输 运性质进行分子动力学模拟，并对在模拟中发现的大波动问题进行了探讨。对于波动性 问题，目前为止大部分研究都是采用收敛后马上截断，从而避免后面波动带来的不必要 的影响；系统收敛后有一段稳定渐变期，此种做法可以避免后面的波动对结果带来的不 必要的影响，由于热物性是对自相关函数的积分过程，而自相关函数自是在无穷远处趋 于0 ，因此在此截断后得到的模拟值会略小于最终模拟值；另外由于固体颗粒的存在， 基金项目：高等学校博士点基金项目(No; 北京市教育委员会科技发展计划重点项目 (No.kz200610016009[0]) 纳米流体自相关函数尤其是更为复杂的纳米制冷剂自相关函数的收敛速率要小于简单 流体，这就要求需要更长的稳定期，因此如何让前期的稳定期延长是个值得研究的课题。 本文从此出发，首先提出需要对稳定性问题进行探讨。 1 分子动力学模拟技术 1.1 扩散系数及相关公式 本文利用 Green-Kubo(G-K)公式对简单流体 Ar 的输运特性进行平衡分子动力学 (EMD)模拟。在EMD 模拟计算中，主要由利用自相关函数计算的G-K 公式和利用均方位 移(MSD)计算的 Einstein 关系式两种等价的公式确定扩散系数。利用速度自相关函数 (VACF)计算的G-K 公式[12]如下： 1∞ D v (t) ⋅v (0) dt (1) 3 ∫ i i 0 相应的归一化速度自相关函数为： v (t) ⋅v (0) i