二维明渠非恒定流的格子Boltzmann模拟-水科学进展.PDF

第 14 卷 第 1 期 水 科 学 进 展 Vol14 ,No1 2003 年 1 月 　　 ADVANCES IN WATER SCIENCE 　　J an . ,2003 　 二维明渠非恒定流的格子 Boltzmann 模拟 1 2 程永光 , 索丽生 ( 1 武汉大学水利水电学院 , 湖北 武汉 　430072 ; 2 河海大学水利水电学院 , 江苏 南京 　210098) 摘要 : 根据二维浅水波方程的特点 , 应用由格子 Boltzmann 方程展开的多尺度方程 , 建立了一个能模拟二维明渠非 恒定流的格子 Boltzmann 模型 。通过对一维溃坝波和二维圆形溃坝波的模拟分析 , 验证了该模型的有效性 。以某水 电站尾水渠为例 , 模拟分析了由机组甩负荷和增负荷引起的水位波动规律 , 为尾水渠的设计提供了依据 。 关 　键 　词 : 格子 Boltzmann 方法 ; 二维明渠非恒定流 ; 二维浅水波方程 ; 尾水明渠 ( ) 中图分类号 : TV 1332 　　　文献标识码 : A 　　　文章编号 : 2003 0100906 ( ) 近年来 , 作为一种新的数值计算方法 , 格子 Boltzmann 方法 Lattice Boltzmann Method , LBM 由于具有算法简 单 、边界条件易处理 、压力直接计算 、适宜并行计算等优点 , 引起人们的广泛重视 , 并在偏微分方程求解 、多 相流 、多组分流 、多孔介质流 、反应扩散流 、生物流场 , 尤其多维不可压流场的模拟方面得到应用[ 1] 。笔者曾 尝试将LBM 应用于一维 、二维水力过渡过程[2 ,3 ] 和一维明渠非恒定流[4 ] , 初步的研究证明该方法有效可行 。本 文旨在发展LBM , 建立其求解二维明渠非恒定流模型 , 并尝试模拟实际的二维非恒定流问题 。 1 　LB 方程和多尺度方程 LBM 是以运动论作为背景的一种数值方法 , 它通过粒子分布函数按 Lattice Boltzmann (LB) 方程的演化来模 拟宏观流动[ 1] 。常用的LB 方程是带B GK 碰撞项的LB GK 方程[ 1] 1 (