MonteCarlo方法研究气液相转变区热力学参数研究.pdfVIP

| ‘ 中国工程热物理学会 991042 工程热力攀与能润利用学术会议 1 MonteCarlo方法研究气液相转变区热力学参犷 1 | 冯晓利 李志信 过增元 (清华大学工程力学系 北京 100084) 1 | | 摘 共:本文用MonteCarlo方法棋拟了二维扭分子系统的热力学乡数。模拟范围包括气、液单相区和 | 气液两相区，侧重于研究气液两相区。模拟在正则 A(4的系综中进行，分子间作用势采用截断的 | ﹄ Leonard-Jones12-6势能。计算得到的等温线显示出与vanderWaals曲线相同的性状，表明所棋拟的二维 分子系统的等温相转变具有连续渐变的特征;并能够以稳定的状态存在于vanderWaals方程指出的宏观 热力学状态的不德定区。依据理论分析得到的二维vanderWaals方程对模拟等温线进行非线性拟合，确 定了该二维灿 阮Waals方程中与温度有关的参数。在同样的正则系综中沿等容线的计算显示了定容比 卜 热的奇异特性. 1 W, MrMCarl.方法 热力学参教 Leonard-Jon.系统 vanderWaals理论 1 | | 1.前言 近十余年来，热科学中的空间/时间微小尺度问题愈来愈受到研究者的关注，如徽细 s 空间的流动与换热、MEMS领域中涌现出的热现象与热问题、微小时间尺度热传导、相 1 | 变与界面动力学等0.=1。围绕这些问题所进行的探索既开拓出了新的研究领域，也引发了 | 对于已有理论的重新思考。半个世纪以来，随着计算机科学技术的进步，MonteCarlo方 法和分子动力学 (MolecularDynamics)模拟方法逐渐发展为现代物理和化学中较为成熟 而有效的工具n[,+1。在热科学中引入分子模拟技术，细致地研究微小尺度(微米一纳米量级 的空间尺度，皮秒一飞秒量级的时间尺度)的热现象11[2,1，不仅有助于揭示现象的物理机制， 而且是联结微观与宏观的有效途径。 MonteCarlo方法由vonNeumann,Ulam和Mertopolis等在二次世界大战末期发展 起来，并由于在计算中大量使用随机数而得名4(1。早期的开创性研究为计算机模拟由分子 /原子组成的系统的热力学行为莫定了基础。基于统计热力学理论和分子间作用势模型， MonteCarlo方法构造随机过程来模拟粒子系统的热力学演进。在对统计系综模拟中，根 据Metropolis重要性抽样原则，由对系统位形进行随机抽样所得到的相空间中的状态点 组成一条Markov链，进而由Markov链中的状态点的统计平均计算系统的热力学函数13,51. MonteCarl。方法的原理决定了其适宜于模拟平衡系统的特性，因此广泛地应用于模拟多 相、多组分流体的相平衡。近来，A.Z.Panagiotopouloslsl综述了针对于各种具体问题的分 衬 . 子间作用势能模型和MonteCarlo模拟技术的最新进展。 二维Leonard-Jones系统是物理吸附的典型模型，具有重要的理论和实验意义161，在 物理化学领域内针对该系统(氢为其典型代表)己经开展了不少研究工作。在对二维 Leonard门】ones系统油固相转变的持续的研究中，J.A.Barker等门将微扰理论(用于液态〕 和自洽胞腔理论(用于固态)的分析与MonteCarlo模拟相结合，得到了形状与三维情形相 *国家自然科学基金资助项 目 I一 55 蠢舞”一 ‘妻 一 -- 一