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