伍亡格子气自动机模型的并行运行.PDFVIP

2O0O年3月 重庆走学学报 (自然科学版) VOl 23 弟23卷弟2期 ~,mmt of Chongqing University(Natural Science Edition) Mar 20o0 文章编号：1000—582x【2000)02—0092—04 1 c【 13一Bit格子气自动机模型的并行运行 c I 墨 龟—篓，高小强 (I 重庆走学热能工程学院，重庆 400044；2 重庆大学材料科学与工程学院，重庆 4O0044) 摘 要：采用按“住”鳊码乖计算 把流体拉子的运动扣碰撞过程处理为运动状态的转贮和节点信 息的回迭两个过程，得 实现并行运算．建立了计算速度·睫、可包含节点量较大的13．Bit格子气自动机 模型的运行机制。焦炭多孔介质流体流场模拟结果表明该机制是正确有效的。 墼 中凰分类号：o 357； TP 311 ； 文献标识码：标 A l L J I J 1972年，法国J．Hardy，Y．Pomeau和0．de Pazzis为 得出流动的微观细节。由它所得的流动的详尽信息， 解决二维流体力学中输运系数发散问题，提出第一个 是其它方法望尘莫及的 。为了能较充分描述流体流 完全离散的格子气自动机模型(Lattice gas automata，缩 动的复杂性和模拟传热过程，文献[3，4]提出和研究了 写为LGA)模型HPP以来… 格子气自动机模型得到 一 种新型的多速模型——13．Bit格子气自动机模型。 迅猛发展。先后产生了单速模型(其中包括FlipⅡ模 1 13．Bit格子气自动机模型 型)、多速模型、高维模型等。LGA模型处理流体流动 问题的特点是将流体及其存在的时间、空间完全离散， 该模型在不改变FHPⅡ模型正六边形网格(如图 把流体想象成由大量只有质量没有体积的流体粒子组 1所示)划分的基础上．只增加一种能量较1速粒子 成，所有这些粒子同步地随着离散的整时间步长，依据 (在一个时间步运动一个空间步长)高的2速粒子(在 给定的简单碰撞规则在网格节点上相互作用，沿网格 一 个时间步运动两个空间步长)，并设计相应的碰撞规 线在节点之间运动。它认为流体运动的宏观特缸是由 则。大量的数值实验证明．13 Bit LGA模型能够有效 微观流体粒子相互碰撞而在整体上表现出来的统计规 模拟典型绕流等复杂几何边界流体流动．模型的通用 律。 性较强” 。初始条件和边界条件易于调整。 LGA模型是微观模型，模拟空间内的微观流体粒 子运动状态瞬息万变。但每个粒子的运动和碰撞须遵 循质量守恒定律、动量守恒定律和能量守恒定律，粒子 运动时还遵循PaⅡ【i不相容原理。模拟一个运动过程 时．预先将大量粒子垒方向或仅在流动方向上随机投 图 I I3-Bit LcA模型的网格划分 放于模拟区域内，格子气自动机经过一段时间的运行， 13．Bit LGA模型源于FHPⅡ模型，但它并不是FHP 粒子的分布状态最终在整体上会达到一种动态平衡， Ⅱ模型的简单推广。引人2速粒子后，粒子的运动和 经过进一步的统计分析．可以获得流体流场及相关宏 碰撞——格子气自动机的运行相对复杂。首先，粒子 观统计量，如雷诺数、速度等，即宏观统计规律。它能 种类增加，不但涉及同类粒子的碰撞，还包括不同类粒 以数值实验方法处理理论分析、传统数值模拟方法难 子的碰