用格子Boltzmann方法模拟平板层流边界层流动.PDFVIP

　 辑第 16 卷第 3 期　　　　　 水 动 力 学 研 究 与 进 展 　　　 . , . 16, . 3 A Ser A V o l N o 　　2001 年 9 月　　　　　　　JOU RNAL O F H YDROD YNAM ICS　　　　　Sep. , 2001 文章编号:(2001) 用格子Bo ltzm ann 方法模拟平板层流 边界层流动 ① 闫广武1, 2 ( 1. 中国科学院计算数学与科学工程计算研究所, 北京 100080; 2. 吉林大学数学系, 吉林长春 130023) 　　摘　要: 　本文应用格子Bo ltzm ann 方法模拟了平板层流边界层流动。考虑到平板边界的影 响, 我们对模型进行了修正, 在标准格子Bo ltzm ann 方程中引入压力修正项。该项满足一定的假定 条件, 从而消除了压缩性引起的偏差。作为算例, 我们用修正的模型计算了平板层流边界层流动, 并与经典的B lasius 解及 L avallee 的格子气自动机结果进行了比较, 结果是可以接受的。 关　键　词: 　格子Bo ltzm ann 方法; 单松弛BGK 模型; 平板层流边界层 中图分类号: 　 357. 4　　　文献标识码: O A 1　前言 格子Bo ltzm ann 方法 (LBM ) 是近几年发展起来的一种模拟复杂系统的新方法[ 13 ] 。这种 方法的基本思想是利用离散的BGK 型Bo ltzm ann 的模型方程式, 通常是将时间、空间与粒子 速度一并离散, 称为格子 Bo ltzm ann 方程, 通过选择平衡态分布函数, 利用多重尺度技术和 展开, 可以得到相应的宏观方程。最早提出和应用单松弛的 模型的是 Chapm an En skog BGK . . [ 1 ] 与 . . [ 2 ] 等, 他们放弃精确的碰撞算子, 引用一种可以控制的碰撞项, 即 S Y Chen Y H Q ian “希望”一种平衡态, 这个平衡态是局部的, 完全依赖于局部的宏观物理量, 并且由此能够得到 宏观的动力学方程, 从而保证了伽利略不变性并去掉了非物理压力。该模型在低速度限制下, 可以得到不可压缩流体的 方程。这完全可以模拟一些不可压流动, 但我们注意 N avier Stokes 到其应用到不可压缩流动中方程有一个偏差, 当去掉这个偏差之后, 才得到不可压缩流动的 [ 4 ] N avier Stokes 方程。在低速时这个偏差不是很明显, 但在可压缩流动中变得严重起来 。本文 在文 [ 1 ] 的基础上, 通过引入附加的碰撞项, 消除了这个偏差, 然后应用该方法模拟了平板层 流边界层流动, 并与经典结果进行了比较。 ① 收稿日期: ( ) 基金项目: 　国家重点基础发展规划研究基金资助项目 G 1999032801 和国家自然科学基金资助 ( ) 项目 (