三维颗粒群沉降的格子Boltzmann模拟.PDF

中国工程热物理学会 多相流 学术会议论文 编号：146091 三维颗粒群沉降的格子 Boltzmann 模拟 查露，曹传胜，栗晶，柳朝晖*，郑楚光 （华中科技大学煤燃烧国家重点实验室，武汉 430074） (Tel: 02787542417, Email: jingli@hust.edu.cn) 摘 要 ： 颗粒悬浮运动研究一直是多相流领域中的热点难点问题，而格子 Boltzmann （LB ）方法是一 种天然并行、易于处理复杂边界的介观方法。本文采用 LB-EBF 方法对三维悬浮颗粒群进行了全解析 直接数值模拟，详细分析了 1152 颗粒（1000 数量级）的沉降规律，具体分为三个方面：三维多颗粒沉 降速度及统计量分析，三维多颗粒沉降位置及统计量分析，三维多颗粒沉降颗粒流场相互作用分析。 关键词：颗粒群沉降；格子 Boltzmann 方法；统计分析；相互作用 0 前言 悬浮颗粒与流体之间的相互作用机理研究一直是多相流领域中的热点问题，流体颗 粒系统广泛存在于化工、环境、能源、生命科学等领域，其与流化床、燃烧、颗粒悬浮 等应用密切相关。研究悬浮颗粒运动中颗粒聚集及扩散机理对于解决实际应用中的问题 起到很重要的作用。无论是从应用角度还是从基础研究角度，研究流体与悬浮颗粒之间 相互作用都非常重要，因此吸引了众多研究者进行实验以及模拟方面的研究[1-3]。 在过去的几十年中，研究颗粒沉降的过程发现颗粒尾涡对颗粒的运动产生重要的影 响作用，其中最重要的是由两颗粒间相互作用而发生的“拖曳-碰撞-翻滚”（DKT）运动。 这种相互作用最早由Fortes et al.(1987)[4]试验发现，在被发现之后，DKT运动被众 多的研究复现，包括实验研究[5]，直接数值模拟（DNS）[6,7]以及格子Boltzmann模拟 （LB）[8]。 实际的颗粒悬浮系统包含了数量非常巨大的有限体积的颗粒，因此研究实际悬浮颗 粒与流体的相互作用是非常复杂及困难的，并且计算效率很低。但是，在模拟中全解析 颗粒周围的流场，并准确捕捉颗粒的运动规律对于研究颗粒与流场及颗粒间相互作用非 常重要。JohnsonTezduyar (1997)[9]利用全解析数值模拟方法模拟了100个颗粒在三 维通道中的运动。Pan et al.(2002)[10]模拟研究了1204个圆球颗粒的流化运动规律。 Kajishima(2004) [11]采用DNS研究了2048个颗粒悬浮运动中颗粒群与湍流场的相互作 用。Xiaolong Yin and Donald L.Koch（2007） [12]研究了悬浮固体颗粒在雷诺数为 1-20范围内的颗粒群受阻沉降速度规律及相应的微观结构。 从上面的研究进展中可以看到颗粒悬浮运动中颗粒与流场相互作用起到了重要作 用，并且多颗粒悬浮运动的研究还很有限，因此，在本文中详细研究了多颗粒悬浮运动 *基金项目：国家自然科学基金 51306065） 的规律，并试图了解颗粒与流场的相互作用对颗粒运动产生的影响。本文采用格子 Boltzmann 方法及 EBF 曲边界处理格式研究了三维大流场（全周期边界）中多颗粒的悬 浮运动过程，解析了颗粒周围的流场，并且捕捉了颗粒的运动规律，详细分析了颗粒群 运动过程中相关统计量的规律。 1 格子Boltzmann方法简介 格子 Boltzmann（LB）方法是以编程简单、天然并行、易于处理复杂边界为特色的 介观方法，数十年来，格子Boltzmann方法在复杂流动等领域有了很大发展[1,13]。本 文采用 LBGK 模型进行模拟，速度模型是 D3Q15，详细介绍见文献[1,13,14]。对于边界 处理部分，考虑到Wu Aidun (2010)[6]提出的附加边界应力方法非常适合于处理运动 边界，本文球形曲面边界处理均使用该方法，详见文献[1,15]。颗粒碰撞部分的处理采 用 Glowinski et al.(2001)[7]提出的处理方法，该方法假设当两颗粒间的距离小于给 定的阀值时，颗粒间将产生排斥作用力，避免颗粒边界相互