气体微通道流动的传递特性.pdfVIP

论文 第51卷第22期2006年11月 讲学屯柱 气体微通道流动的传递特性 王利民①② 葛蔚①4 陈飞国①② (①中国科学院过程工程研究所多相复杂系统国家重点实验室，北京100080；②中国科学院研究生院，北京100049． 4联系人，E—mail：wge@home．ipe．ac．cn) 摘要采用微观的拟颗粒模型研究了努森数在o．01～0．20范围内气体微通道流动的速度和温度分布．发 现速度分布主要受努森数和所加外场力的影响，随努森数增加，壁面滑移速度先增后降，而温度分布也 受外力场的显著影响．达-酉(Darcy)摩擦因子随努森数增大而增大，而它随马赫数的变化与所谓层流向 湍流转捩提前的趋势相似． 关键词 微通道 拟颗粒模拟 努森数 马赫数 微流动研究起源于对紧凑、高效的微尺度换热器 对应关系．因此，在一种简化的物理模型下考察微观 的开发…．微流体控制元件具有不受电磁场干扰的特 流动的表现，为建立实际流动的完整模型作某一方 性，已应用于精密控制系统和医学的药物注射装置【2J． 面的探讨也是有益的． 微反应系统可以强化反应过程的混合、传热和传质， 对此，我们注意到，粒子模拟方法可以不受很多 提高反应产率，也已应用于硝化反应、氢气催化燃烧 实验条件的限制，建立相对简单的理想化模型，直接 J． 从离散的物理图景就问题的不同方面分别深入探讨， 等的研究p 在微流动中，常规尺度下可忽略的力和现象却 促进这方面的研究【111．目前，对于过渡流区和自由分 可能起着重要的作用．如液体会呈现颗粒性状，气体 子流区的传递行为主要采用直接模拟蒙特卡罗 会呈现稀薄效应，而边界也会“移动”等[4】．因此传 统的流动和传递理论需要不同程度的修正才能适用． 拟，对稠密和高度非平衡的气体流动的模拟还存在 目前这方面的研究还不太成熟，阻碍了微化工技术 不少困难．而更严格的分子动力学(MD)方法则会因 的进一步发展． 为引入了很多分子及其与边界间作用的细节而使问 J 由于现象的复杂性，微流动的实验研究也存在 题复杂化．本文将采用微观的拟颗粒模型(PPM)upl7 诸多困难．比如气体微通道流动的阻力特性就是一 从物理上直接模拟简单的硬球气体的微通道流动， 个颇有争议的问题．Wu和Little【56]的测量发现气体流 过微小槽道时阻力系数远远高于层流理论计算值． 和导热的基本规律，为更复杂条件下的实验与模拟 李志信等【_7]的试验研究表明，微细圆管气体流动范宁 提供相比较的基准． 摩擦系数大于层流理论计算值．而Pfahler等人【81在 1模拟方法及条件 测量氮气和硅油流过硅蚀刻微通道时，流体阻力系 数都小于层流理论计算值．Choi等[9]发现氮气在微细 拟颗粒模型是结合分子动力学模拟【18】和直接蒙 圆管中的阻力系数明显小于常规尺度下的理论值． 特卡罗模拟方法提出的一种粒子方法．如图1所示， 拟颗粒没有内部结构，每个拟颗粒有4个属性，即质 此外，Wu和Little[5，6】在气体流过微小槽道的实验中还 发现层流向湍流转捩的临界雷诺数大大减小．而对 于转捩提前，现在一般认为是由于微槽中的较大相 系统，质量和半径是定值．模拟的每个时间步中，粒子 对粗糙度造成的u…． 从本质上说，影响强非平衡的微观流动的因素 并且进一步靠拢的各组拟颗粒按硬球方式处理其碰 要远多于宏观和近平衡的流动．如果只按少数几个 撞，即碰撞后的速度为： 准则将一种模型对应于某一试验或另一模型，而不 _l，1，：¨+2坠盟兰掣(B一只)，‘ 考虑许多目前还难以定量的分子尺度的细节对流动