格子Boltzmann方法模拟微尺度流动和传热研究.pdf

中国：I=程热物理学会 传热传质学 学术会议论文 编号：063025 格子Boltzmann方法模拟微尺度 流动和传热干．． 田智威，邹春，刘红娟，闫锦霞，郑楚光 (华中科技大学煤燃烧国家重点实验室，湖北武汉，430074) Tel：(027-8411 E-mail-zwtiantian@163．com 得壁面边界的速度滑移和温度阶跃，在含有粘性热耗散的热格予模型的基础上。提出了一种新的直接 基于宏观量的边界处理格式．模拟得到的速度场和温度分布与解析解吻合得相当好，充分说明了本文 采用的模型和边界处理的合理性．同时在模拟中还发现：对于不同的Kn数，均存在使得其上壁面的 温度阶跃为零的临界Ee数，并且其临界值均在3．0附近． 关键词：格子Boltzmann方法；微尺度Couette流；速度滑移I温度阶跃 1、引言 近年来自然科学及工程技术发展的一个重要的趋势是朝微型化迈进，尤其是微电子 机械系统(MEMS)技术的飞速发展，更极大的推动了这一研究热潮【11。而在所有的微 电子机械或纳米器件的设计及应用中，流动和传热又是非常突出而重要的闯题。随着格 el 边界的模拟【3．5】。其中Niea1．首先采用半步长反弹格式来实现固壁面的速度滑移，而 el口，．则提出了将反弹格式和自由滑移的镜面反弹格式按一定比例结合的边界处理 Tang 格式。其次，在模拟中需要解决的另·个问题就是建立Kn数与松弛时问f的关系式。 Lee et el alo[7】 a1．【6】提出了具有一定物理意义的统一的无量纲Kn数定义式。然而Guo 指出其关系式均忽略了微尺度下的壁面限制效应对松弛时问的影响，迸而提出了有界限 el 系统的松弛时问r‘与Kn数的关系式。而对于非等温的微尺度流动的模拟，Sofonea at．is]针对有限差分的多速热格子模型，提出了漫反射边界处理格式以获得滑移速度和温 度阶跃。但是需要指出的是；其需要增加一额外的系数(文中取的是1．15)以获得与理 论值吻合的较好的滑移速度。 本文基于He的含粘性热耗散项的热格二f模型【9】，对非等温的微尺度Couette流动 和传热进行了模拟。．在边界处理上，采用了新的宜接基于宏观量的边界处理格式以更好 模拟速度滑移和温度阶跃。 2、模型及边界处理 2．1、隐格式热格子模型 ’ 为了模拟微尺度的热Coueue流，首先介绍He提出的隐格式热格子LB模型【9】： +基金项tJ：围家蕾点基础研究发胜；}划t973)项目(编譬：2004CB217703) 585 (1) z(x+c，出，H出)一z(x，’)=一赢【肥‘)一∥(x，‘)】 。 其中新分布函数Z和磊的定义如下j6 (3) Z：Z一坐(Z—Zq)，蟊：＆一—0．5—At(＆一∥)+o．5A斫z, 1 1 J