微气体表面粗糙效应的格子Boltzmann模拟

论文摘要

随着微机电系统工业和微制造业的发展,对微流体的研究受到越来越多的关注。对于微通道中的液体流动的研究主要集中在表面张力、壁面可湿度等因素的影响,在理论或实验上的研究也取得了不少的进展。对于微通道中的气体流动的研究,情况又有些不同。因为在微通道的气体流动中,由于连续介质假设不成立导致基于Navier-Stokes方程（NSE）的解的可靠性受到了质疑。在微气体流动或稀薄气体流动的数值模拟中,我们需要求解一些基于分子（伪粒子）的模型,包括分子动力学（MD）模拟,Boltzmann方程（BE）,直接Monte-Carlo模拟（DSMC）、信息保存（IP）方法等。近年来,格子Boltzmann方法（LBM）提供了一个有效的研究微流体系统的新工具。关于格子Boltzmann在微流体或稀薄气体上的研究已经取得了一些很好的成果。与分子动力学、直接Monte-Carlo模拟或其它同类型的方法相比,格子Boltzmann方法在计算效率上,在处理复杂边界问题上具有相当的优势。然而,尽管在小尺度中粗糙效应重要性增加,但由于缺乏有效的模拟工具,到目前为止对高Knudsen数的微气体流动中的粗糙效应的研究成果仍相对较少而且不能统一。为了研究氮气和氦气在微通道中的流动,本文设计了一个微通道模型,其高度为H=1.2μm,长度为L=5H或L=10H,对应氮气和氦气的Knudsen数分别为KnN2=0.055和KnHe=0.16。微通道粗糙表面的几何形状被分别设计成方形、Sin函数形、三角形和分形形状。利用这个模型,我们做了下面主要工作:1)在微流体格子Boltzmann模型的状态方程方面,我们从气体运动论出发,重新推导松弛时间因子以及它和Knudsen数的关系。通过设置具有物理意义的网格时间步长,使得Knudsen数的表达式与网格精度无关。同时,通过与实验比较我们验证了这一修正的松弛因子跟原模型相比可以消除高Knudsen数的微气体的人工粘性,得到和实验相比更吻合的而且具有明确物理意义的结果。2)在微流体格子Boltzmann模型的边界处理方面,我们考虑在微通道中由于摩擦或者流体杂质沉淀而形成的表面粗糙具有分形性质。所以我们在微流体格子Boltzmann方法中引进分形粗糙面并使用具有物理意义的参数如分维来控制粗糙面,我们认为它比其他学者使用的规则型,如矩形或三角形粗糙面更接近真实情况。此外,我们还计算并讨论了漫反射边界条件中粗糙表面的切向动量容许系数的影响,得到其在格子Boltzmann模型中对应的取值。3)利用修正的微流体格子Boltzmann模型,我们模拟了微通道中的氮气和氦气的流动,捕捉到了最小Knudsen数现象。我们计算并分析了固壁表面粗糙效应和粗糙面形状对流体速度、压力、质量流量以及摩阻系数的影响,讨论了微通道中粗糙效应与稀薄效应与可压缩效应的耦合作用,得到一些有意义的结论。对氮气和氦气的模拟表明表面粗糙不但影响微通道中的速度分布,而且影响其压力分布。随着相对粗糙度的增加质量流量快速下降,摩阻系数快速增加。微通道中稀薄性效应、可压缩性效应和表面粗糙效应的作用强烈耦合,综合影响宏观物理量的变化趋势。在微通道中氮气和氦气的稀薄性、可压缩性和表面粗糙效应强烈耦合。在模拟分形粗糙度时,我们也分析了分形层次和网格精度的影响。除了相对粗糙度之外,所使用的粗糙几何形状对结果的影响也很大。这些不同于宏观流动的特性对了解微流气体流动有着重要的意义。

论文目录

• 摘要
• ABSTRACT
• 主要符号说明
• 第一章 绪论
• §1.1 研究背景
• §1.2 微通道气体流动研究概况
• §1.3 格子Boltzmann方法在Knudsen流体中的应用
• §1.4 论文提出与模型建立
• §1.5 本文的主要工作和创新点
• 第二章 稀薄气体动力学
• §2.1 基本介绍
• §2.2 Knudsen数和其它基础概念
• §2.3 Boltzmann方程
• §2.4 气体和固壁的作用
• §2.5 微气体流动与稀薄气体流动的相似性条件
• §2.6 影响微通道气体流动的主要因素
• 第三章 修正的格子Boltzmann模型
• §3.1 格子Boltzmann方法与理论
• §3.2 平衡分布函数的推导
• §3.3 Knudsen流体的格子Boltzmann模型
• §3.4 光滑边界微通道中的气体流动特征
• §3.5 修正的格子Boltzmann模型的松弛时间
• §3.6 修正模型的验证
• §3.7 小结
• 第四章 粗糙边界条件的处理
• §4.1 压力边界条件
• §4.2 分形边界的生成
• 4.2.1 迭代函数系统法
• 4.2.2 独立分割法
• §4.3 格子Boltzmann模型的固壁边界的处理
• 4.3.1 向后反弹法
• 4.3.2 向前向后反弹结合法
• 4.3.3 散射核边界条件
• §4.4 切向动量容许系数的测量
• §4.5 不规则边界的处理
• §4.6 分形层次和网格精度的影响
• §4.7 小结
• 第五章 微气流的粗糙效应的数值模拟
• §5.1 气体参数和计算条件
• §5.2 对速度的表面粗糙效应
• §5.3 对压力分布的分形表面粗糙效应
• §5.4 对质量流量的表面粗糙效应
• §5.5 对摩阻系数的表面粗糙效应
• §5.6 粗糙表面几何形状对速度剖面的影响
• §5.7 粗糙表面几何形状对摩阻系数的影响
• §5.8 粗糙效应与稀薄性和可压缩性的耦合作用
• §5.9 小结
• 第六章 总结
• §6.1 本文工作与结论
• §6.2 后续研究工作展望
• 参考文献
• 攻读博士学位期间成果
• 致谢