格子Boltzmann方法流体计算的新算法及其可视化研究.pdf

学术研究水平／CAE自主仓U新发展战略 I 格子Botzmann方法流体计算的 新算法及其可视化 钱跃兹 李青 (上海大学力学所) (上海大学计算机工程和科学学院) 1引言 1．1格子Boltzmann方法 1．1．1格子气自动机模型 Gas 格子气自动机(Lattice 的具体化。20世纪70年代，Hardy，de ——．HPP模型。该模型采用二维正方形格子，每个格点上可驻留有流体粒子。在同 一个格点上的流体粒子按照一定的规则进行碰撞，并在一个时间步内沿网格线运动 到相邻的格点。在碰撞和迁移过程中，HPP模型遵循一定的物理规律，即保持质量 守恒、动量守恒——这正是流体宏观运动方程描述的规律。因此，HPP模型能够反 映出流体运动的一些基本特征。但是，由于正方形格子缺乏足够的对称性，HPP模 型对应的宏观方程不能反映正确的非线性性和耗散效应。 LGA还必须满足另一个基本条件：格子需要有足够的对称性。在二维情况下，他们 使用等边三角形网格(或称为正六边形网格)，并由此构造了一个新的LGA模型～ 一FHP模型。可以证明，在～定条件下从FHP模型可以得到宏观尺度下的不可压 Navier—Stokes方程。 由于采用布尔型变量，LGA模型具有绝对的数值稳定性；同时，所有粒子的碰 撞和流动均同时发生，且粒子之间的相互作用上局部的，因而LGA模型方法具有本 质的并行性，非常适合在并行计算机上实现；在固壁边界上，流体粒子反弹回流场 内部，因而LGA在处理复杂边界时十分简便。这些都是LGA模型方法不同于传统数 ⑩ CAE自主创新发展战略＼学术研究水平 值方法的优点。 但是，由于在运行过程中会引入随机噪声，模拟结果往往含有统计噪声。同时， 从LGA导出的宏观方程与标准的流体力学方程之间存在一定的差距，如：不满足伽 利略不变性、压力与速度相关等，这些缺陷是本身不可克服的。 l-1．2格子Boltzmann方法 布尔变量的统计平均量(称为单粒子的分布函数)代替粒子是否存在的布尔量进行 (LatticeBoltzmann Method，LBM)。但是，这一模型仍然采用LGA的碰撞方式， 碰撞操作仍然具有指数复杂性。 1．1．3 LBGK模型 H 1992年，YQian(钱跃兹)等人提出格子Boltzmann方法的单松弛模型，即： 格子BGK(Lattice 撞操作用趋于某一平衡态的松弛过程代替。 LBGK模型极大地简化了模型的计算，并且在一定的条件下，可以从该模型导出 方法的优点。该模型的提出使格子模型方法的研究和应用达到了一个新的水平。目 前，LBGK模型上格子Boltzmann方法中最主要的，也是应用最广泛的模型。 本系统以2维9速LBGK模型(D2Q9)为基础，实现二维流场的计算及其可视 化系统。 1．2科学计算可视化 科学计算数据的高产出与传统数据分析的低效率形成尖锐的矛盾。这一矛盾严 重地影响科学计算的效率和质量，阻碍了超级计算的持续发展。 of is purpose “科学计算的目的是洞察，而不仅是获得数值”(The computing ins numbers)。科学计算可视化基于日趋成熟的计算机图形学和图象处 ight，not 理技术，充分发挥人的视觉潜力，以图形、图象及动画等视觉表现形式展现计算与 数据的本质，并允许与计算和数据进行交互，使人们能够观察到不可见的对象(see ⑧ 学术研究水平／CAE