学位论文 > 优秀研究生学位论文题录展示

微纳尺度气体流动和传热的Burnett方程研究

作　者: 包福兵
导　师: 林建忠
学　校: 浙江大学
专　业: 流体力学
关键词: 微纳尺度流动 滑移-过渡流区 Burnett方程 数值模拟方法
分类号: O354
类　型: 博士论文
年　份: 2008年
下　载: 302次
引　用: 1次
阅　读: 论文下载
 

内容摘要


在微纳尺度气体流动中,通常用Knudsen数,Kn,来表征气体的稀薄程度,Kn数定义为气体分子的平均自由程与流场特征尺度之比。微纳尺度气体流动中,流动进入滑移过渡流区,流场中缺乏足够的分子碰撞,流体逐渐偏离热平衡状态。此时,基于线性本构关系的Navier-Stokes(N-S)方程得到的结果与实验结果有较大的偏差,而基于颗粒模型的DSMC方法统计噪声大,收敛慢,需要大量的计算时间。本论文采用Burnett方程来研究微纳尺度通道里的气体流动和传热特性。Burnett方程由Chapman-Enskog展开从Boltzmann方程导出,能够较好地描述稀薄气体的非线性本质,近二十年来,Burnett方程逐渐受到关注。首先采用线性小扰动理论,首次系统分析常规Burnett方程、增广Burnett方程、BGK-Burnett方程以及Woods方程的一维稳定性问题。研究发现这四种方程在没有引入对方程中随体导数的近似时,对于小扰动都是不稳定的,不稳定的临界Kn数分别是0.105,0.105,0.158和0.074。当采用Euler方程来近似这些方程里的随体导数时,增广Burnett方程和BGK-Burnett方程对小扰动是无条件稳定的,而常规Burnett方程和Woods方程是不稳定的,临界Kn数分别是0.499和0.184。当用N-S方程来近似这些方程里的随体导数时,所有四个方程对于小扰动都是无条件稳定的。接着采用增广Burnett方程结合二阶滑移边界条件研究微纳尺度下的平板Couette流动,在边界上采用松弛方法,成功地获得了任意网格尺度时所有Kn数下的收敛解。结果显示Burnett方程的结果与DSMC以及IP方法的结果符合得较好,当Kn数较小时,N-S方程的计算结果与Burnett方程的结果基本符合,但是随着Kn数的增大,N-S方程的结果开始偏离DSMC结果,而Burnett方程仍旧能够较好地符合,Burnett方程更适合模拟高Kn数下的微纳尺度流动。文中还分析比较了不同Kn数和马赫数下的流场特性,给出了壁面上的滑移速度和温度跃变随Kn数和马赫数的变化。最后采用Burnett方程结合高阶通用滑移边界条件研究微纳尺度下的平板Poiseuille流动,首次在Burnett项中引入松弛方法,成功地拓宽了Burnett方程的求解范围,获得了Kn≤0.4时Poiseuille流的收敛结果。比较发现,Burnett方程的结果与理论分析结果、实验结果以及DSMC结果符合得较好。当Kn数比较小时,Burnett方程和N-S方程的结果相符,但是当Kn数变大,达到过渡流区时,两者开始出现差别,N-S方程失效,而Burnett方程仍旧能够符合DSMC结果。文中还分析了入口与壁面具有相同和不同温度时的气体流动和传热特性。本论文首次系统研究了Burnett方程在微纳尺度通道气体流动中的应用,拓宽了Burnett方程的求解范围。把连续性方法延伸到过渡流区,从而提高了计算效率。本文的研究为微纳机电系统的流动和传热研究提出了新的思路,为后续研究奠定了基础。

全文目录


摘要  4-6
Abstract  6-8
符号说明  8-13
第一章 绪论  13-43
  1.1 研究背景  13-20
    1.1.1 微机电系统  13-15
    1.1.2 纳机电系统  15-17
    1.1.3 微纳尺度流动研究的重要性  17-20
  1.2 研究现状  20-35
    1.2.1 微纳尺度流动特点  20-22
    1.2.2 微纳尺度流动的流场参数  22-24
    1.2.3 微纳尺度流动的模拟方法  24-30
    1.2.4 DSMC方法  30-33
    1.2.5 Burnett方程  33-35
  1.3 本文的主要内容和创新点  35-43
第二章 Burnett方程和边界条件  43-65
  2.1 二维增广Burnett方程  43-48
  2.2 三维增广Burnett方程  48-53
  2.3 其他类型Burnett方程  53-54
    2.3.1 原始Burnett方程  53
    2.3.2 Woods方程  53-54
    2.3.3 BGK Burnett方程  54
  2.4 滑移边界条件  54-60
    2.4.1 滑移边界条件表达式  57-58
    2.4.2 切向动量适应系数  58-60
  附录  60-62
  参考文献  62-65
第三章 Burnett方程的稳定性分析  65-83
  3.1 Burnett方程稳定性分析介绍  65-66
  3.2 常规Burnett方程的稳定性分析  66-72
  3.3 增广Burnett方程的稳定性分析  72-75
  3.4 Woods方程的稳定性分析  75-78
  3.5 BGK Burnett方程的稳定性分析  78-81
  3.6 本章小结  81-82
  参考文献  82-83
第四章 Couette流动和传热的模拟  83-105
  4.1 微纳尺度的Couette流动  83-85
  4.2 方程和边界条件  85-94
    4.2.1 控制方程  85-87
    4.2.2 滑移边界条件  87-91
    4.2.3 方程求解  91-94
  4.3 Burnett方程和IP方法结果的比较  94-97
    4.3.1 DSMC方法和IP方法的结果  94-95
    4.3.2 Burnett方程与IP方法结果的比较  95-97
  4.4 结果与讨论  97-102
  4.5 本章小结  102-104
  参考文献  104-105
第五章 Poiseuille流动和传热的模拟  105-131
  5.1 微纳尺度Poiseuille流动  105-106
  5.2 Burnett方程的求解  106-115
    5.2.1 流场求解方法  106-107
    5.2.2 滑移边界条件  107-109
    5.2.3 源项的处理  109-111
    5.2.4 程序的验证  111-115
  5.3 结果与讨论  115-127
    5.3.1 与其他方法结果的比较  115-119
    5.3.2 入口和壁面温度一致时的结果  119-124
    5.3.3 入口和壁面温度不一致时的结果  124-127
  5.4 本章小结  127-129
  参考文献  129-131
第六章 总结和展望  131-135
  6.1 全文总结  131-133
  6.2 下一步工作展望  133-135
攻读博士学位期间发表的论文  135-137
致谢  137

相似论文

  1. 复杂地形风场CFD模拟方法研究,TK81
  2. 电晕枪放电特性的机理研究,O531
  3. 覆冰分裂导线舞动数值模拟方法研究,TM751
  4. 概率理论与随机数值模拟方法在海岸泥沙淤积预测中的应用,TV148
  5. 黑洞的拟正则模,P145.8
  6. 鄂尔多斯盆地上古生界盆地分析模拟,P618.13
  7. 格子Boltzmann方法及其对振荡层流边界层的模拟,O35
  8. 表面界面生长中的两因素的理论探讨,O485
  9. 预应力锚索数值模拟方法及动力承载特性研究,TU452
  10. 数值模拟方法若干问题的研究,O242
  11. 结构损伤数值模拟研究,TU312.3
  12. 自然空泡流数值模拟方法研究,U661.1
  13. 超声速引射器的数值模拟方法及其引射特性研究,TJ95
  14. 方坯连铸宏观传输现象复合数值模拟的研究,TF777
  15. 东中国海海洋动力环境要素的数值模拟方法研究,P731.21
  16. 速度滑移及其对微纳尺度流动影响的分子动力学研究,TB383.1
  17. 含损伤热粘塑性本构数值算法和铝锂合金动态响应研究,O346
  18. 纯引力力相互作用多粒子系统的统计性质研究,P172
  19. 飞秒激光抽运—探测热反射系统建立及金属薄膜热输运过程研究,TK123
  20. 二阶边界滑移气体薄膜挤压特性研究,O354
  21. 大气风场模型研究及应用,O354

中图分类: > 数理科学和化学 > 力学 > 流体力学 > 气体动力学（可压缩流体力学）
© 2012 www.xueweilunwen.com