有限差分格子Boltzmann方法模拟稀疏气固两相流研究.pdfVIP

中国工程热物理学会 学科： 多相流 第十一届学术会议 编号： 056085 l 有限差分格子Bo 固两相流· 陈胜，柳朝晖，贺铸，郑楚光 (华中科技大学煤燃烧国家重点实验室，武汉，430074) 1H：027 Email：isdon影lle@163．com 摘要 本文采用一种新的耦合有限差分格子Boltz毗11Ⅱ模型，模拟了气固后台阶下的气固两相流。其 中的颗粒相使用一种新的指教Lagr柚gi姐方法进行跟踪．由于考虑了气固两相耦合作用，所得结果与 实验结果和其他的数值模||；【结果非常吻合．并且与以往的格子Boltz毗nn模型相比，计算量大大减少。 关键词有限差分格子B01tzmann，双向耦合，气固两相流，后台阶 气固两相流是日常生活和工业生产中非常常见的一种流动现象，如在垃圾燃烧及气 力管道运输中。在这类流动中有两种特别值得关注的现象：颗粒弥散和颗粒对流场的改 由于颗粒的存在，对流体分离点和重新汇合点的改变非常明显。而流体分离和重新汇合 所形成的回流区对工业设备的运行性能有很大影响。为了提高生产效益，虽然已开展了 大量的工作来研究这类气固两相流动的特性，但由于该流动本身的复杂性，目前相关的 工业装置的设计和控制工作主要还是依靠经验模型。随着计算机技术的发展，数值模拟 方法在气固两相流研究中起的作用越来越明显[1]。 法。欧拉一拉格朗日方法指把流体(气相)当作连续介质而将颗粒视为离散体系，流场 (气相)使用欧拉方法求解而颗粒(固相)使用拉格朗日方法求解。而在欧拉一欧拉方 法中，颗粒与流体都被看成是共同存在且相互渗透的连续介质。由于把颗粒作为拟流体， 所以无论是气相还是周相都是使用欧拉方法求解，如常用的双流体模型(t釉fluid model)【2—4]。 有十几年的发展历史，因此不可避免的会有一些局限性。其中最大的问题是格子 主要原因有二：一是格子BoltzⅡIann方法本身收敛速度较慢； 法一般采用均匀网格对流场进行剖分。尤其是第二种原因，虽然均匀格子有利于格子 物理量在靠近流场边界的区域变化一般比较剧烈，对这些区域需要使用较细的网格离 散：而在远离边界的区域物理量变化相对比较平缓，使用较大的网格就可以描述流场了。 ·教育部新世纪创新人才计划项目资助项目 277 如果对正规区域采用同样的网格离散，自然就会降低算法的效率，造成资源浪费，或者 不能细致的描述流动。虽然有很多学者开始研究基于非均匀网格的格子Boltzmann模型 拟两相流的工作依然很少。因此，在本文中我们提出了一种新的在非均匀网格中处理相 模型在工程领域的拓展提供新的尝试。 1．数值方法 有限差分格子Boltzmann模型的演化方程如下： a，Z+e『a。Z=Q(Z) (1．1) 函数。f是松弛时间。为了体系出两相间的作用力，我们将平衡态分布函数由标准的 z=删+警+型拿手盟】 nz， 改为： 抽p[1+等+竖％盟圳n，， 离散速度。矿的定义为： 屯=％+加。 (1．4) 其中n为待定常数，而口。=疋／p。p是流体密度，，是气固两相间的作用力。式(1．3) 中晟后一项^f表达式为 肌堂等盟 ms， 其中参数ck同n一样，由宏观方程确定。通过标准的多尺度展开，我们可以得到，为 了还原出正确的Navier—stokes方程，有