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

气固两相流的高效格子Boltzmann算法设计及其直接数值模拟

作 者: 韩海锋
导 师: 郑楚光
学 校: 华中科技大学
专 业: 热能工程
关键词: 气固两相流 格子Bolztmann方法 高性能计算 非牛顿流体 热传导 直接数值模拟
分类号: TK12
类 型: 博士论文
年 份: 2013年
下 载: 41次
引 用: 0次
阅 读: 论文下载
 

内容摘要


气固两相流在自然界和许多工业过程中都扮演着十分重要的角色。颗粒和流体的相互作用对两相流理论具有重要意义,对此过程的内在机理的认识、探讨和模型预报是气固两相流体力学的核心问题。通过实验来研究颗粒对流场的影响,会有相当大的局限性。作为一种新的计算流体动力学技术,格子Boltzmann方法(LBM)已经成为研究气固两相流的有力工具。但是,用格子Boltzmann方法直接模拟气固两相流时,由于颗粒占据一定的流场体积,这样造成计算时间和内存消耗问题特别突出,LBM代码的优化成为一个重要的课题。因此,我们提出了一个索引算法来处理这个问题。在此算法的基础上,我们对气固两相流的一些前沿问题进行了一些有益的探索。本文的研究工作有如下几个方面:首先,为改善现有算法的缺点,我们设计了一个新的LBM优化算法索引算法。该算法通过索引(指针)的轮换,可以把二维或者三维问题转化为一维问题,经过此步骤后,通过碰撞和流动的耦合,一个计算格点仅需要两个临近的格点即可完成一次迭代。索引算法的内存消耗只是Two-lattice算法的一半。与其它最快算法设计相比,索引算法在二维问题方面的加速比达到2.3X,三维情形的加速比也达到1.3X,而且,它在缓存性能方面综合性能也是最优的。此外,该算法可以很容易地用于具有复杂边界的气固两相流的直接模拟,通过对1024个颗粒的沉降问题模拟,发现使用Intel(R) Xeon(R) E5645CPU,计算速度为20.19MLUPS。其次,为了消减颗粒表面由于网格离散引起的离散误差,我们利用LBM多松弛模型在边界处理方面的优点,通过对比单颗粒在静止状态和运动状态的计算结果,给出了使用MRT模拟颗粒运动的最佳参数取值。再次,研究了单颗粒在牛顿流体中的沉降,分析了单颗粒最终沉降雷诺数和颗粒运动状态的关系,并对管道宽度对颗粒运动的影响进行了分析。然后,分析了双颗粒在牛顿流体中沉降的相互关系,给出了双颗粒三种运动关系以及它们速度,角速度和位置变化。接着,研究了单颗粒和双颗粒在幂律非牛顿流体中的运动,分析了不同非牛顿指数对颗粒运动轨迹的影响,并进一步分析了非牛顿指数引起的流场涡量场的变化。最后,我们对单个恒温颗粒在非等温流场中的沉降进行了研究,根据管道宽度的不同,分两部分进行研究。首先对管道宽度为颗粒直径4倍的情形,根据颗粒运动的状态、平衡位置以及尾涡状态等对Grashof数进行分类;其次对管道宽度为颗粒直径1.5倍的情形,管道壁面对颗粒运动的影响加强,我们发现了一个新现象:颗粒的沉降速度随着Grashof数的增大而变小,当Grashof数为4083时,颗粒的最终沉降速度为0,达到了热悬浮状态。当Grashof数继续增大时,颗粒先向下沉降,然后改变速度方向,向上运动。

全文目录


摘要  4-6
Abstract  6-10
1 绪论  10-21
  1.1 研究背景以及现状  10-12
  1.2 格子 Boltzmann 方法在气固两相流中的应用  12-18
  1.3 本文研究目的与内容  18-19
  1.4 全文安排  19-21
2 气固两相流的格子 Boltzmann 方法  21-50
  2.1 格子 Boltzmann 方法基本模型  21-32
  2.2 格子 Boltzmann 方法外力项的处理  32-34
  2.3 气固两相流模拟的若干关键问题  34-49
  2.4 本章小结  49-50
3 格子 Boltzmann 方法高效算法设计  50-68
  3.1 格子 Boltzmann 方法简介及其程序设计  52-54
  3.2 索引算法  54-60
  3.3 计算性能和缓存表现  60-64
  3.4 索引算法模拟颗粒群的性能  64-66
  3.5 本章小结  66-68
4 颗粒在牛顿流体中的运动  68-85
  4.1 格子量纲与物理量纲转换  68
  4.2 网格离散误差处理  68-71
  4.3 Segré-Silberberg 效应  71
  4.4 牛顿流体中单颗粒的沉降  71-76
  4.5 牛顿流体中双颗粒的沉降  76-82
  4.6 本章小结  82-85
5 颗粒在非牛顿流体中的运动  85-95
  5.1 非牛顿流体的格子 Boltzmann 模型  85-88
  5.2 程序验证  88-89
  5.3 单颗粒在非牛顿流体中的沉降  89-91
  5.4 双颗粒在非牛顿流体中的沉降  91-94
  5.5 本章小结  94-95
6 非等温气固两相流直接数值模拟  95-113
  6.1 控制方程  95-96
  6.2 格子 Boltzmann 非耦合双分布热模型  96-98
  6.3 程序验证  98-103
  6.4 位于方腔偏心位置具有温度的颗粒产生的自然对流  103-105
  6.5 恒温单颗粒的自由沉降  105-111
  6.6 本章小结  111-113
7 总结与展望  113-117
  7.1 全文总结  113-115
  7.2 研究展望  115-117
致谢  117-118
参考文献  118-128
附录 1 攻读学位期间发表的学术论文  128-130
附录 2 攻读学位期间参加的学术会议  130-131
附录 3 攻读学位期间参与的科研项目  131-132
附录 4 Ladd 模型的宏观模型推导  132-140

相似论文

  1. 喷动床内气固两相流动特性的研究,TK173
  2. 三维非稳态热传导边界元方法研究及数值系统开发,O241.82
  3. 打散回转筒内流场和温度场数值模拟,TK173
  4. 平房仓储粮温度变化规律的研究,S379
  5. 基于诱导风流机制下的粉尘运移规律研究,X513
  6. 基于精细积分算法的瞬态热传导分析的IMLPG方法,TK124
  7. 方腔内气固两相流动的格子Boltzmann模拟,O359
  8. 细粉加料装置料腔内气固两相流场的数值模拟分析,TK124
  9. CFB锅炉炉内流化与防磨的数值模拟与实验研究,TK229.6
  10. 基于多核SMP集群环境的光线追踪模拟卫星成像并行研究与实现,TP391.41
  11. 管道内颗粒悬浮流中大颗粒和湍流相互作用的直接数值模拟研究,O359
  12. 肺动脉及其分支血液动力学的数值模拟研究,R312
  13. Cu/DT9/Cu叠层材料的轧制复合与数值模拟研究,TG335
  14. 基于热传导原理的气化炉内高温测量技术研究,TQ546
  15. 基于DEM方法的散体颗粒热传导模拟,TK124
  16. 流化床气固两相流数值模拟及实验研究,TK221
  17. 热传导方程反问题的若干方法研究,O551.3
  18. 两种正则化方法在二维逆时热传导问题中应用,O241.82
  19. CFD在工程中的一些应用,O35
  20. 二维方腔内非牛顿流体的格子Boltzmann模拟,O373
  21. 热传导方程的一种自适应有限元算法,O241.82

中图分类: > 工业技术 > 能源与动力工程 > 热力工程、热机 > 热力工程理论
© 2012 www.xueweilunwen.com