非常规颗粒（颗粒群）曳力的格子Boltzmann模拟-热能工程专业论文.docx

独创性声明 本人声明所呈交的学位论文是我个人在导师的指导下进行的研究工 作及取得的研究成果。尽我所知，除文中已标明引用的内容外，本论文 不包含任何其他人或集体已经发表或撰写过的研究成果。对本文的研究 做出贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到 本声明的法律结果由本人承担。 学位论文作者签名： 日期： 年 月 日 学位论文版权使用授权书 本学位论文作者完全了解学校有关保留、使用学位论文的规定， 即：学校有权保留并向国家有关部门或机构送交论文的复印件和电子 版，允许论文被查阅和借阅。本人授权华中科技大学可以将本学位论文 的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫 描等复制手段保存和汇编本学位论文。 保 密□，在 年解密后适用本授权书。 本论文属于  不保密□。 （请在以上方框内打“√”） 学位论文作者签名： 指导教师签名： 日期： 年 月 日 日期： 年 月 日 华中 科技 大学 华 中 科 技 大 学 硕 士 学 位 论 文 I I 摘 要 气固两相流广泛存在于工业生产中，特别是燃煤电厂锅炉中存在大量复杂的非规 则颗粒（颗粒群）。本文采用格子波尔兹曼方法对非常规颗粒及颗粒群在不可压粘性流 体中的性质进行了研究，重点研究了颗粒所受的曳力。具体而言，本文的主要工作和 得到的结论有： 研究了圆柱颗粒和圆球颗粒在均匀流中的流场性质及受力特征，计算了不同 Re 数 时颗粒的曳力系数、升力系数及斯特劳哈尔数等无量纲特征数，并和文献结果进行了 比较，证实了 LBM 模拟的可靠性。 研究了四种非规则截面颗粒（三角形、十字形、三叶形及矩形）在斯托克斯流中 的产生的压降。三角形颗粒和三叶形颗粒具有相似的几何结构，放置角度为 0°时颗粒 压降最大，随着放置角度增大，压降依次减小。十字形颗粒压降对放置角度最不敏感， 而矩形颗粒对放置角度最为敏感。对比四种非规则颗粒压降可知，矩形在 0°放置时压 降最小，三叶形（ε=2）在 0°放置时压降最大。 研究了简单颗粒阵列、八颗粒阵列及随机分布颗粒群的流场性质及曳力。使用硬 球蒙特卡洛（HSMC）方法对颗粒群进行初始化。定量验证了简单颗粒阵列所受曳力 随体积分数的变化规律。对于八颗粒阵列，中心位置处的主流方向速度随着 Re 数的增 大而增大，且与实验结果吻合较好。对于随机颗粒群，当体积分数较大时，空隙处速 度概率密度峰值出现在接近 vz=0 处；当体积分数减小至 0.35 时，概率密度出现双峰分 布的趋势；当体积分数继续减小，概率密度变为单峰分布。对颗粒群中颗粒所受平均 曳力进行定量分析并比较不同理论公式，发现 Koch Sangani 公式更加精确，但需修 正其分段函数转折点为 0.2。此外，颗粒群中的单个颗粒无量纲曳力的概率密度符合平 均值为 1、标准差为 0.518 的高斯分布。 研究了二维点源颗粒团聚体和三维有限体积颗粒团聚体在斯托克斯流体中所受曳 力。得到了以下结论：对于二维颗粒团聚体，当颗粒数目较少时，团聚体各向同性特 II II 征不显著，不同入流角度下团聚体所受曳力差别较大。此外，对不同体积分数下的团 聚体所受曳力进行计算，然后利用已有的二维圆柱颗粒曳力理论公式进行数据拟合， 得到了团聚体曳力的计算公式。对不同颗粒数目的三维颗粒团聚体进行模拟，拟合得 到了团聚体曳力和颗粒数目以及投影面积之间的关系式。 初步研究了两种特殊工况下颗粒在流场中的受力性质，即颗粒壁面滑移效应及颗 粒运动。模拟了存在壁面滑移效应的二维圆柱绕流，计算了不同 Re 数和 Kn 数下的圆 柱颗粒曳力。模拟了单颗粒在二维通道内的沉降过程，对比颗粒最终沉降速度和文献 相吻合，并探讨了不同颗粒-流体密度比时颗粒所受曳力及曳力系数的变化情况。 关键词：格子 Boltzmann 方法；曳力；非规则颗粒；颗粒群；颗粒团聚体；滑移边界； 颗粒运动 PAGE IV PAGE IV Abstract Gas-solid tow-phase flow is very common in various industrial practice. In particular, there exists a lot of complex irregular particles and multi-particles in coal-fired power plant boilers. The properties of irregular particle and multi-particles in impressible viscous flow were investigated in this paper, especially the drag force of particles. To