微流动的同伦分析法解析求解研究.pdfVIP

中国T程热物理学会 传热传质学 学术会议论文 编学：073234 微流动的同伦分析法解析求解 张田田， 贾 力， 王志成 (北京交通大学机械与电予控制工程学院动力系，北京100044) (联系电话：010 Email：1{ia@center．nltu．edLl．cn) 摘要：引入相似变换将微尺度滑移流动的控制方程变换为非线性常微分方程。采用通论分析法解析求 解f该方程，通过’j文献结果对比验证r本文求解结果的正确性。 关键词：微流动，速度滑移，同伦分析法，相似变换 1前言 当前微米纳米尺度传热学已经成为热科学领域中最为激动入心的前沿之一，人们认 识到，当微尺度和瞬态作用时间小于一定数值时，传统热和流体理论将不雨适用于描述 或偏离所观测到的现象…。7在O．001砌0．1范围内，稀薄效应造成在壁面处流体分予 的速度与壁面边界条件不同，此种现蒙称1为速度滑移。在此范嗣内，仍可利用 针对滑移流动与换热进行了大麓的实验和数值研究。但是，由于具有滑移边界的控制方 程是高度非线性的偏微分方程，求解十分困难。因此，解析解仍十分匮乏。 廖世俊教授1992年提出了求解高度非线性方程的同伦分析法。…。它是建立在拓扑 理论的同伦方法的基础上，其有效性与所研究的方程是否含有小参数无关。这个方法对 许多类型的非线形常、偏微分方程均有效。同伦分析法本文采用同伦分析法对微尺度滑 移流动进行了解析求解。 2滑移流动的数学描述 本文对二维无限大平行平板间滑移流动进行了研究。假定平行平板间高为2鼠给定 进口速度，出口压力设为大气压，流体与擘面温度一致，不考虑换热。壁面速度滑移。 由于流场关丁-平板中心线对称，冈此，为简化计算将计算区域取对称的一!k此外，假 定为二维稳态层流流动；忽略体力、可压缩性和粘性加热的影响：动餐渊借系数取为l。 根据上面的物理模型，建立对于不可压缩气体的控制方程。表达如F： 连续性方程 盟+盟：o ax a)) 动餐方挥 V——二 Lz』(2) “ “塑+v卫：～上盟+u堡～———o+U—__ 一+ a 8x av D x 8V! 边界条什： 壁面 y=O 掰=力a掰／ay (3) 中心线 y=∥ a“／ay=O v=O (4) 其中，旯为分子平均自由程， Ⅳ、盼别为置y方向的速度，以7运动黏度，p为密度， p为压强。 3相似变换 引入相似变换变量： 77=吉 哟)=M+《抓列 (5) 由上述相似变换可将控制方程及相应的边界条件转化为如下常微分方程 。／…+够“一召(厂’)2：墨日。 (6) ，p 边界条件： 壁面上 K聆·厂(o)一／1(o)=o‘ (7) 中心线 ／”(1)=o厂(1)=O (8) 4同伦分析方法求解 考虑非线性方