管排绕流的格子Boltzmann模拟.pdfVIP

中国工程热物理学会 燃烧学学术会议 编号：984032 管排绕流的格子Boltzmann模拟+ 阮剑郑楚光邢华伟柳朝晖 (华中理工大学煤燃烧国家重点实验室武汉430074) 熊圣武 (武汉大学软件工程国家重点实验室，武汉430072) 摘 要 作为格子气方法的改进，格子Boltzmann方法在保留了格子气模型的～些优点的同时，克 服了格子气模型的不足之处。应用格子Bohzmann方法，可以对更复杂的流体进行模拟。本文 阶精度上表现为标准的Navier．Stokes方程。用该模型对圆柱绕流和管排流动进行模拟的结果显 示，该模型能较好地模拟复杂流动现象，并具有一定工程应用背景． 关键词： 格子Boltzmann方法，Navier—Stokers方程，数值计算，管排，圆柱绕流 0．引言 对流动过程的准确描述和模拟，是揭示和预报燃烧过程乃至许多工程问题的基础 和关键。长期以来，人们一直在寻找着对流动有效的模拟方法。由于计算机技术的飞 速发展，一些数值计算方法(如差分法、有限体积法、有限元法等)已经在工程上得 到了较好的应用。但由于流动本身的复杂性、计算能力的限制和各种模型固有的局限 性在很大程度上阻碍了这些方法对某些复杂流动模拟的准确性和有效性。与传统的计 算流体方法不同，八十年代后期形成的格子气(LG)方法基于原始的物理背景，从运 动层次论上对流体进行描述，用若干粒子表示流体微元，并满足有关物理关系(如守 恒定律、对称性等)。这些粒子在宏观系综平均上较好的反映了一些流体运动现象。LG 方法具有数值稳定、边界条件容易引入、程序简单，可以大规模并行实现、以及可方 便运用可视化技术等优点。但LG方法亦表现出某些问题：(1)需要用大量的空间和 时间平均来消去其统计噪声，这极大的影响了LG方法的计算效率；(2)对流项不是 p的函数：(3)格子气状态方程中的压力项依赖于速度，这将产生一些非本质的流体 的上述三个缺点，并保留了LG方法的一些优点。 作了一定的改善，而最有影响的是sY．Chent”和Y．H．QianE4】等几乎同时提出的单一时间 TimeRelaxation 松弛方法(theSingle Method)，并应用它广泛地研究了二维、三维流 体力学。华中理工大学煤燃烧国家重点实验室采用基于单一时间松弛方法的FHP模型， 在对分离板后的流场等流动的模拟中取得了有意义的成果【51。 本文中，我们提出了更灵活、更自然的九点正方形格子Boltzmann模型，并用此 模型模拟了流体力学中典型的圆柱绕流以及更具有实际意义且边界也更加复杂的管排 绕流，得到了较好的结果。 1．九点正方形格子Boltzmann模型 格子Boltzmann模型的动力学演化方程是 ，，(x+e。，，十lJ=／，(x，rJ+s2t (1) 国家自然科学基金资助项目 IV--1 Q．为碰撞函数，表示由于碰撞过程而引起^(x，r)的改变率。格点处的局部宏观流 体密度、速度由下列各式确定 p=∑工，pu=∑厂P， (2) 动力学演化方程(1)把流体运动分成两部分：迁移过程和松弛过程。在迁移过程中， 厂，(z，，)沿8．方向运动并与其邻居交换质量：在松弛过程中，每个格点处的速度分 布函数向它的局部平衡状态松弛。为简单起见，松弛过程采用BGK单松弛模型。 。 1 Q，=一。(f，一f，”) (3) f