粗糙度对微细圆管内流动特性影响的摄动分析研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：073316 粗糙度对微细圆管内流动特性影响的 摄动分析 张春平唐大伟 韩鹏祝捷 中国科学院丁程热物理研究所北京 100080 联系电话：010 Email：z垫塑g￡迎韭卫ingj!@bQ地鱼!!￡Q啦 摘要：本文以瓤糙度为小参数，采用正则摄动法求解了微细圆管内的层流流动。计算结果表明，壁面 相对粗糙度和粗糙元的无量纲间距是影响微细通道流场的主要阕素。壁面相对粗糙度仅影响扰动流函 数的峰值但不影响流场受扰区域的大小：I『l：『粗l糙曲线窄问波数对扰动流函数的峰值和扰动区域的人小 都有影响。当裥糙度相同减小窄问波数时，流函数的扰动峰值和受扰区域都将增人，这为合理设计拳fi 糙元同时抑制摩擦阻力增大提供了理论依据。 关键词：粗糙度；丁F则摄动；微细通道 1引言 微通道内的流动与常规尺度一卜．的流动有较大的区别，从众多文献的研究来看，一个 很重要的原冈来自于壁面粗糙度对流动的影响。而且近年来随着计算技术和模拟技术的 发展，使用数值模拟研究微尺度条件‘卜．各种新的物理现象已发展成为一类重要的方法。 李志信、于玮等川数值模拟了交错布置和对称布置的规则突起粗糙元对微尺度平板间层 流流动的影响，发现粗糙元大小、分布密度及布置形式均影响流动阻力。Koo和 Kleistreuer【玉3 J用PML(PorousMedja Layer)模型研究了壁面粗糙度对微通道流动阻力的 的影响。 相对粗糙度可以定义为粗糙元高度与流道定性尺寸的比值，微细通道的相对粗糙度 基本在O．0l一一O．05范围内，与定性尺寸相比为小参数。文献[5】提出壁面有微小变形可用 正则摄动法进行求解。众多文献也在不同的研究中采用了这种方法[6~8]。基下上述思想， 本文采用止则掇动分析方法结合数值模拟，对二维长商微细圆管内肇面粗糙度对层流流 动影响的规律进行了研究。 2基本控制方程 幽l围管微通道内粗糙j÷分布物理模型 ·微细圆管内粗糙元分布物理模型如图l所示。假发管道充分艮，直径为D，且流动 为定常流动。微细圆管内表面粗糙度高度为七俐，假设粗糙元上卜_分布对称且为周期函 数形式分布，则相对粗糙度七俐慨可以设为函数／阳与一个小参数f的乘积，其中函数 ／M是周期函数。则圆管上表面点的坐标为，／2+f俐。 项目名称：中科院知识创新丁程系要方向性项日(KGcxl一sw．12)——微系统荇十前沿技术的研究 881 由物理模型可知圆管微通道内的流动为典型的二维流动问题。取特征长度为微管直 径D，特征速度为截面平均速度瓦特征压强为p矿2对二维圆柱坐标下的N—S方程进行 无量纲处理，圆管上半区域的无量纲方程和边界条件建立如‘卜．： 矗——+1，——=一上+——f—彳+一——f，．——?|| c，，、‘／ 舐 毋 缸 Re l舐2 r西L毋／l 矗罢+y雾=一塞+去[窘+吉昙(，．考)] 甜～+v一=一二+～l一+_卜_I厂V川 (2) 苏 加 加Re l叙2升Lr西、’／l 甜妻+V≯争去眵托釉)l 丝+三盟：o (3) 反，咖 zf=v=O，，．=l／2+矿(x)；甜，1，∞，r=O (4) 式中Re为雷诺数，甜，v，p分别为x、r方向的速度和压力。