微通道内流动沸腾的耗散动力学模拟.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号: 083157 微通道内流动沸腾的耗散动力学模拟 * 谢荣建，马虎根 ,黄万鹏 (上海理工大学 热工程研究所，上海 200093) 联系电话：021 E-mail:mahugen406@ 摘要：耗散动力学模拟是一种介观尺度的数值模拟方法，它适合模拟复杂的流体运动。本文使用了一 种两相的耗散动力学模型来模拟微通道内流体沸腾，统计得到了微通道内流体的平均速度分布图，在 粒子分布较均匀的情况下，速度分布呈现 Poiseuille 流，之后逐渐上移； DPD 粒子流动过程的形状 趋势图，出现部分变浓的间隙流动，最后稀化区占主导的现象。 关键词：耗散动力学 微通道 流动沸腾 1. 引言 微通道内两相流动沸腾高的热换系数使得冷却和加热更加迅速, 因此已经被广泛地 应用于微电子机械(MEMS 和 MST) ，计算机 CPU 的冷却，和燃料电池等对空间要求 较为苛刻的行业领域。迄今为止人们对于微通道的传热机理、干涸(dry - out) 现象等的 认识还相当有限,许多有关的研究尚处于起步阶段,并且对同类现象的实验结果及分析还 没有统一的认识，这使得微通道内沸腾换热问题成为近年来热物理学科的又一研究热 点。 在微观尺度上, 分子动力学通过追踪每个流体分子的位置和动量从而对流体的流动 行为进行模拟。虽然分子动力学在模拟较少数量的分子系统方面取得了相当的成功, 但 是由于计算机速度的限制, 目前分子动力学所能计算的时空尺度还限于纳秒、纳米级。 因此，发展介观尺度的流体模拟技术是近年来国际计算机模拟领域的发展趋势。介观尺 度的流体模拟技术试图克服传统方法运用到复杂流体系统时所面临的上述固有困难。介 观尺度的流体模拟方法中利用系统的特征时间尺度和空间尺度来简化整个模拟过程。最 早的介观尺度流体模拟技术是气体格子法(lattice gas approach), 随后的 lattice-Boltzmann 方法由于采用了分布函数而有较大的改进。 耗散粒子动力学(dissipative particle dynamics,简称 DPD) 是一门新兴的介观尺度数 值模拟技术, 于 1992 年由 Hoogerbrugge 和 Koelman 提出[1], DPD 结合了分子动力学 和气体格子法的优点, Espanol、Warren 和 Marsh 等人的工作[2,3]奠定了耗散粒子动力 学在统计力学方面的基础。在 DPD 系统中, 基本的单元是一些离散的被称为“粒子” 的动量载体, 这些粒子在连续的空间和离散的时间上运动, 每个粒子在一定范围内和周 围粒子发生相互作用, 粒子之间的相互作用力包括保守力、耗散力和随机力三种。DPD 系统中所谓的“粒子”具有粗粒化(coarse-graining)的概念, 每颗粒子的运动代表的是大 量分子( 即所谓的流体粒子) 的集体行为。与前述的介观尺度的流体模拟技术不同,DPD 系统中的粒子在连续的空间而非离散的格子上运动, 因此遵循Galilean 不变性。另一方 面, DPD 又可以理解为宏观的微分流动控制方程在小尺度上的随机描述[4], 从这个意 义上说, DPD 方法是连接微观分子动力学方法和宏观流体力学方法的一座桥梁, 是一种 真正的介观尺度的模拟技术。 DPD 方法其优越性在于能够比较方便地实现对复杂流体系统的模拟。近年来，有 学者开始使用 DPD 模拟相变的过程。2000 年，Pagonabarraga et al[5]