用格子Boltzmann方程模拟浅水波问题力学学报.PDFVIP

第　29 　卷　第　5 　期 力 　学 　学 　报 　 Vol. 29 , No. 5 1997 　年　9 　月 ACTA MECHANICA SINICA Sept . , 1997 用格子 Boltzmann 方程模拟浅水波问题 1) 施卫平　　胡守信　　阎广武 (吉林大学数学系力学教研室 , 长春 130023) ( ) 摘要 　提出了用格子 Boltzmann 方程 LBE 模拟浅水波问题的方法. 通过无粘气体运动方 程与浅水波方程的比较 , 确定了 LBE 方法中平衡态的形式 , 使宏观方程与浅水波方程一致. 计算了二维浅水波的一个问题 , 数值结果与精确解做了比较. 关键词 　格子 Boltzmann 方程 , 浅水波 引　言 　　1986 年 , U. Frish , B. Hasslacher 和 U. Pomeau[1 ] 提出了一种六角形的 FHP 格子气自动 机 (简称 L GA) 模型. 早期的计算表明 , 无论是运算速度 , 还是计算精度都有很大的提高. L GA 模型用布尔量表示粒子的分布函数 , 每一节点最多允许 7 个粒子存在 , 并且满足泡利不 相容原理. 粒子在节点处发生碰撞 , 碰撞过程满足质量 , 动量守恒. 宏观量通过统计平均获 得[2 ] . L GA 模型模拟 NavierStokes 方程存在几个问题 : 低速度限制 ; 宏观流体方程中对流项 系数一般不等于 1 ; 计算结果存在噪声. 格子 Boltzmann 方程[3 ,4 ] (Lattice Boltzmann Equation , 简称 LBE) 是在 L GA 方法基础上 提出来的 , 它取实数型的格子气 , 即将粒子分布函数由布尔量变为 [0 , 1 ] 范围内的实数 , 碰 撞过程可通过选择特定形式的平衡态分布函数来计算. 平衡态的选取要求满足 : ①各向同性 , ②质量、动量守恒 , ③宏观量满足流体力学方程. LBE 方法很好地解决了 L GA 方法的困难. LBE 方法对应的宏观方程与 NavierStokes 方程一致 , 可计算高 Re (雷诺) 数问题 , 计算结果 噪声小. 碰撞过程是局部的 , 适合并行计算. LBE 方法用于模拟不可压缩流体问题已取得一 些进展. 本文对 LBE 方法进行了改动 , 修改了平衡态分布函数 , 使宏观方程与浅水长波方程 一致. 作为数值试验 , 计算了“抛物面水波”问题 , 并与精确解进行了比较. 1 　BGK 型 FHP 模型 如图 1 , FHP 单松弛模型将二维空间离散化成正六角形 , 每一个节点通过 6 条格线与相邻 的 6 个节点相连 , 粒子沿此 6 条格线方向运动. 7bit 的 FHP 模型中 , 粒子可沿 e1 , …, e6 6 α α 个速度方向运动, 其速度的大小| e | = C , = 1 , …, 6 ; = 0 表示静止粒子, | e | = 0. 粒 α 0 子分布函数的演化方程为 1 ( ) ( ) ( ) eq ( )