颗粒绕流流动与受力特性的格子Boltzmann模拟研究-热能工程专业论文.docxVIP

内蒙 内蒙古科技大学硕士学位论文 I I 摘 要 颗粒流体系统广泛存在于工业中，例如流体流经固定颗粒床、固体颗粒流态化 及密闭管道中的气力输送等。两相流动往往伴随传热、传质和化学反应等过程，致 使颗粒与颗粒、颗粒与流体之间的动力学特性相当复杂，需要对颗粒流体复杂的局 部环境及颗粒流体的作用机理进行研究以便更有效地设计、改进和控制生产过程。 本文以颗粒流体系统的基本单元——双颗粒、四颗粒绕流现象为研究对象，数 学模型采用格子 Boltzmann 方法 D2Q9 模型，计算域出、入口及上、下壁面的边界处 理采用非平衡外推法，颗粒曲表面的处理采用结合空间插值的非平衡外推法，流固 表面受力的计算采用动量交换法。基于 C++平台自编程序，研究了不同颗粒间距、 颗粒排布方式和不同雷诺数下的流动模式与受力特性的变化规律。研究结果表明： （1）Re=200 时，并列排布下，流动模式按间距比可分为单体模式、间隙流模 式、平行涡街模式，临界间距比分别为 1.2 和 2.2。不同流动模式内，两颗粒的平均 曳力系数（CD）均随着间距比的增大而减小。平均升力系数（CL）随间距比的增大 而减小，表明间距比越大，上下排颗粒互斥作用越弱。串列排布时，流动模式按 间距比可分为单体模式、再附着模式、双脱落模式，临界间距比分别为 2.0 和 4.1。 在单体模式和再附着模式内，上游颗粒的 CD 值随着 L/D 增大而减小，下游颗粒 的 CD 值为负值，此时两颗粒相互吸引。在双脱落模式内，上、下游颗粒的 CD 值 均随着 L/D 增大而缓慢增大，表明两颗粒间引力减弱。 （2）四颗粒绕流流动模式及受力特性与Re、间距比及排布方式有关。正形 排布时，流动模式分为单体模式、稳定屏蔽模式、摆动屏蔽模式及涡脱落模式四 种，临界间距的划分受Re的影响，Re越小越滞后。流动模式的转变是CD、曳力系 数均方根、升力系数均方根发生突变的主要原因。随间距增大，上游颗粒CD值渐 小，而下游颗粒CD值渐大，表明上游颗粒所受曳力的一部分转移到下游颗粒上。 菱形排布的流动模式随颗粒间距比的增大可以分为单体模式、狭窄间隙流模式、 涡冲击模式三种，临界间距的划分同样受Re的影响，Re=200时临界间距为1.2和 2.5，Re=100时为1.2和3.5。流动模式的过渡会导致各颗粒受力系数的显著跳跃。 本文的研究结果为进一步研究大规模、形状复杂、排布方式混乱、高雷诺数下 的颗粒群的输运特性提供了参考。 关键词：格子 Boltzmann 方法；双颗粒；四颗粒；流动模式；受力特性 II II Abstract The particle-fluid system is widespread in the industry, such as fluid flows through a static bed filled with particles, the fluidization of solid particles, and pneumatic conveying in closed pipe. The process of two-phase flow in particle-fluid system is often accompanied with heat transfer, mass transfer, and chemical reactions. As the complexity of kinetics characteristics between particles and particle-fluid, it is necessary to research on the local environment between particles and the mechanism between particles and fluid, which will help to facilitate the design, improvement and control of the production process. In this paper, cross-flow over the basic unit of particle-fluid system―double particles and four particles were researched. Lattice Boltzmann method with D2Q9 model was employed as the mathematical model. The non-equilibrium extrapolation was used to process the boundary of the computa