基于反转法的二氧化碳势能计算及迁移性质预测.pdfVIP

中国工程热物理学会 工程热力学与能源利用 学术会议论文 编号：091058 基于反转法的二氧化碳势能计算及迁移 性质预测 王晓坡 宋渤 吴江涛 刘志刚 （西安交通大学动力工程多相流国家重点实验室，710049，西安） 联系电话：029 E-mail: wangxp@mail.xjtu.edu.cn 摘要：以常见气体二氧化碳为计算对象，采用反转法迭代得到了二氧化碳的新势能参数。在此基础上， 根据分子动力理论，计算了二氧化碳在零密度下的迁移性质，包括黏度η、自扩散系数D 和热扩散因 子α ，计算的温度范围为300-3500 K 。与文献比较表明计算结果可以满足实际工程应用。 0 关键词：二氧化碳；热物理性质，反转法，新势能 Calculation of the New Potential Energy Function of Carbon Dioxide by Inversion Method Wang Xiao-Po, Song Bo, Wu Jiang-Tao, Liu Zhi-Gang (State Key Laboratory of Multiphase Flow in Power Engineering, Xi’an Jiaotong University, 710049, Xi’an) Abstract: A new potential energy function of carbon dioxide was obtained by means of inversion method. According to the kinetic theory of gas, the transport properties of carbon dioxide, including viscosity η, self-diffusion coefficient D and isotropic thermal diffusion factor α , were obtained in temperature between 0 300 K and 3500 K at zero-density. Keywords: Carbon dioxide; Transport properties; Inversion method, Potential energy function 0 前言 流体的热物性数据在工程应用和科学研究中具有重要作用，通过实验测试能够得到 比较精确的热物性数据。但是在实际应用过程中，经常会用到高温下（大于1000 K）的 热物性数据，这时实验测试就变得非常困难。利用分子热力学方法推算高温流体热物性 数据越来越受到广泛关注，其最关键的问题就是需要可靠的分子间作用力描述。势能参 资助项目：国家自然科学基金重点项目资助）；全国优秀博士学位论文作者专项资金 （200540 ） 数主要是通过精确的实验数据（比如维里系数、气相黏度、音速等）回归得到。但目前 常用的势能模型一般都是两个或三个可调参数且精度不高，同时，由于模型的近似性， 通过不同性质回归得到的势能参数往往也不相同。因此，Smith 等人[1]提出了反转法来 计算势能参数，该方法不依赖于特定的势能模型，只需假定一初始参数，通过迭代即可 得到物质的新势能。近年来，Najafi 和Papari 等学者[2,3]对这种方法进行了研究改进，提 高了其计算精度。 本文以常见气体CO2 为对象，通过反转法得到了CO2 的新势能，在此基础上根据分 子动力学理论，计算了二氧化碳在零密度下的迁移性质，包括黏度η、自扩散系数D 和 热扩散因子α0 ，计算的温度范围为300-3500 K 。 1 理论模型 [4]