任意复杂流固边界的格子Boltzmann处理方法.pdfVIP

物理 学报 ActaPhys．Sin．Vo1．63，No．7(2014)074703 任意复杂流．固边界的格子Boltzmann处理方法冰 史冬岩 ) 王志凯 ) 张阿漫2)十 1)(哈尔滨工程大学机电工程学院，哈尔滨 150001) 2)(哈尔滨工程大学船舶工程学院，哈尔滨 150001) (2013年10月17日收到；2014年1月6日收到修改稿) 本文提出了一种适用于流固耦合领域中任意复杂边界条件的latticeBoltzmann处理方法．该方法基于 half-way反弹模型，在流固耦合处构建了一层虚拟边界，并结合有限差分的方法，获取虚拟边界上的变量值． 改进后的方法确保了粒子反弹位置与宏观速度采集点的位置相同，计入了实际物理边界与网格线不重合时， 偏移量对计算结果的准确影响，而且其适用范围被扩展到了任意静止或运动、平直或弯曲的复杂边界．文中研 究了该方法在Poiseuille流、圆柱绕流和 Couette流等经典条件下的边界处理能力，结果表明该方法与理论值 符合 良好，且当实际物理边界与网格线不重合时，与已发表文献中的结果相 比，具有更高的精度． 关键词：格子Boltzmann方法，边界处理，流固耦合，虚拟边界 PACS：47．11．Qr，47．10．ad，47．15．Rq DOI：10．7498／aps．63．074703 领域 中对气液边界、液液边界、液固边界处理上所 1 引 言 表现出的优势，己得到相关领域学者们的广泛关 注 [15一l71． 在 处 理 复 杂 流 体 动 力 学 问题 时格 子 到 目前为止，已有多种基于LBM 的边界处理 Boltzmann方法 (LBM)作为一种相对崭新的数值 模型 [18,19]．从LGA一直被沿用到LBM的标准反 方法正得到迅速的发展 [112】．相对于传统的数值方 弹模型，虽然能够在保证质量守恒的情况下，在边 法，起源于格子气 自动机 (LGA)和连续Boltzmann 界处获得非滑移边界条件．但1995年Noble等发现 方程的LBM，采用 自下而上的方式，通过首先确定 该模型只有一阶精度，不满足原latticeBoltzmann 微观层面上粒子的分布，然后积分运算获得宏观 方程在低mach数条件下，时间和空间上的二阶精 控制方程 [0 】．1992年，Qian等提出了BGK模型， 度要求 ．half-way反弹模型是标准反弹模型的 通过多尺度展开，获得了具有二阶精度 的不可压流 一 种变形，该模型在实际物理边界与其 内侧第一 场的Navier—Stokes方程3[]_拓展了LBM在计算流 层流体粒子之间的1／2处构造了一层虚拟边界．流 体力学领域 中的应用． 体粒子在运动 了0．5 后到达该虚拟边界，然后发 LBM在模拟复杂流场时具有三点比较明显的 生反弹，经过下一个0．5 返回到原来的位置．1997 优势．第一，天然的可并行性；第二，易于理解和编 年，He等验证了half-way模型具有二阶精度 [21]． 程实现 [6】．第三，灵活的多相边界可处理性 [7】_相 修正反弹模型是标准反弹模型的另一种变形，虽然 对于有限体积法对于液固边界的处理、边界元方法 同样具有二阶精度，但不同于前两种模型不允许在 对于气液边界的处理 [8-11】，LBM能够从微观层面 边界格点上执行碰撞操作，该模型在保证碰撞前指 上，对表面张力、黏度系数和滑移速度等物理参量 向流场内部的分布函数已知的前提下，允许边界格 进行灵活、有效的控制 [12-14]．特别