应用格子Boltzmann方法模拟街道峡谷内的流场.pdfVIP

第24卷第4期 阜阳师范学院学报(自然科学版) Vo1．24．No．4 2007年12月 Journal of Fuyang Teachers College(Natural Science) Dec．2007 应用格子Boltzmann方法模拟街道峡谷内的流场 樊冬梅 ，姚云飞。 (1．安徽农业大学理学院安徽合肥 230036；2．阜阳师范学院数学与计算科学学院安徽阜阳 236041) 摘 要：用格子Boltzmann方法模拟了入口速度u。一0．1的对称以及不对称街道峡谷内不同时刻的流场，从模拟结果 可以看出街道峡谷内一个逆时针漩涡的形成过程．与传统方法的模拟结果相比，具有很好的一致性，从而说明用格 子玻尔兹曼方法模拟街道峡谷流场是完全可行的． 关键词：Boltzmann方法；街道峡谷流场；一致性 中图分类号：035 文献标识码：A 文章编号：1004—4329(2007)04—0033—03 0 引言 B。ltzmann方法模拟街道峡谷内的流场是可行的· 1 数值算法 格子 Boltzmann方 法 (Lattice Boltzmann Method，LBM)E1]是一种完全离散和并行运行的局 在实际模拟中，我们采用了郭照立等人的不可 部动力学模型，它是使用简单的微观模型来模拟流 压的二维九速模型(D2Q9)，离散速度方向： 体的宏观行为的一种方法，作为一个求解流体系统 r(0，0) =0 偏微分方程的有力工具，是一种与高性能计算机匹 ．{(cos[(／一1)~r／2]，sin[(／一1) 2]) i=1,2，3，4 配的新型流体计算方法．由于格子Boltzmann方法 【^、／ (f。 [( 一1)~r／2+ ／4]，sin[( 一1)~r／2+ ／4]i=5,6，7，8 具有天然的并行性，因此是一种能够与高性能计算 演化方程： 机紧密耦合的完全并行算法．自问世以来，格子 g ( 4-ce At，l 4-At)一 gi( ，1)一 Boltzmann方法已经引起物理学家、计算机科学家 一 {[g ( ，￡)一g／ (z，￡)]， 与数学家们越来越大的兴趣． 由于格子Boltzmann具有天然的完全并行性， 局部平衡态分布函数为g 一aip 4- (“)，其中a。 很好的数值稳定性和易于实现的边界条件，目前，除 一 一 4a／c。，当i一1，2，3，4时a0一一4a／c。，嘶一 了在一般的流体力学问题中得到了成功的应用外， ／c ； 格子Boltzmann方法在湍流Ez]、多孔介质流[33、多相 当 i一5，6，7，8时a —y／c。， 流[4]、反应和燃烧系统[53、磁流体力学[6 等相关领域 其中C—Ax／At为粒子速度，Ax和 是单位网格步 也得到了比较成功的应用． 长和时间步长，压力P—c ID，参数 ， ，7 满足 格子Boltzmann方法在模拟城市街道峡谷问题 + y一 ， + 27一 ， 中也取得了重大进展[7]，但模拟的大多是单个建筑 且 物周围流场，本文用不可压的DZQ9模型来模拟两 个建筑物的街道峡谷内的流场，在给定高宽比和峡 “ 一 [3 5 5 ]， 谷宽度的情况下，实验结果显示了在峡谷内部出现 其中 (DO一 4 了一个逆时针的漩涡的全部过程，与传统方法模拟