用格子Boltzmann方法研究二维Burgers方程.pdfVIP

天津城市建设学院学报 第l8卷 第 1期 2012年3月 JournalofTianjinInstituteofUrbanConstruction Vo1．18 No．1 Mar．2012 用格子Boltzmann方法研究二维Burgers方程 张 伟，李文杰 (天津城市建设学院 能源与机械工程系，天津 300384) 摘要：给出了一个求解二维Burgers方程的格子Boltzmnan方法．分析 了相应格式的单调性和稳 定性，得到了格式的单调性条件 ，并证明了在此务件下，格式是 稳定的．研究表明，通过引 入适当形式的平衡分布函数，采用简单的5速正方格子模型可以恢复宏观Burgers方程．通过与 标准二阶精度的有限差分解的比较，证 明了本文的格子 Boltzmann方法(LBM)是简单而有 效 的． 关 键 词：格子Boltzmann方法；Burgers方程；平衡分布函数；单调性；稳定性分析 中图分类号：O411 文献标志码：A 文章编号：1006—6853(2012)01．0036．05 格子Boltzmann方法 (1atticeboltzmannmethod， 考虑如下形式的二维非线性对流扩散方程，即二 LBM)是近十几年发展起来 的一种新颖 、强劲 维黏性Burgers方程 (robust)和具有广泛应用前景的数值格式，已经被广 a，+ P=v(O +a ) (1) 泛用于模拟各种流体流动问题，这种格式在涉及界面 p(x，0)=Po() (2) 动力学和复杂边界等复杂流动系统中的应用获得了 式中： 是黏性系数． 特别 的成功 1【-4】．与传统数值方法相 比，格子 选择计算量小，边界条件容易实现的5速正方格 Boltzmann方法具有物理图像清晰，算法具有本质并 子构造如上Burgers方程的LB模型，LBE中的碰撞 行性 ，边界条件易处理以及程序简单等独特的优 算子，采用由 Qian[和 Chen等㈣分别独立提出的 点．因此受到了越来越多国内外学者的关注． LBGK模型．通过选择适当的平衡分布函数，利用 对流扩散方程是物理、化学和环境科学等领域中 Chapman．Enskog多尺度展开技术，给出了一个基于 经常用到的一类数学模型．格子Boltzmann方法为求 LBGK模型的LBM． 解对流扩散方程提供 了一个新 的思路 ．在格子 Boltzmann方法中，由于 5速或 (4速)正方形格子缺 乏足够 的对称性而无法由格子Boltzmann方程 (LBE) 1 格子Boltzmann方法 恢复正确的Navier．Stocks方程 J，因此，在 LBM 的 研究中经常采用较为复杂的7速正六边形格子或9速 I．I 平衡态分布函数 正方形格子．但是，一些研究表明，对于某些形式的 首先，将二维空间离散成均匀的正方形格子，5 对流扩散方程，采用简单的格子模型不仅可以恢复宏 个离散速度是c=ce，i：01一，4．其中，c= 为粒 观方程，而且数值模拟结果是正确有效的．Shen[61在 △f