部分填充多孔介质通道内流动的格子boltzmann研究.pdfVIP

2013年第2期 节 能 (总第365期) ENERGYCONSERVATION 一15一 部分填充多孔介质通道内流动的 格子Boltzmann研究 姬翔，陈宝明，刘芳，雒婉 (山东建筑大学热能工程学院，山东济南250101) 摘要：通过编写格子Boltzmann程序来对部分填充多孔介质的通道内的流体流动进行模拟，在雷诺数 为200的条件下改变了障碍物形状排列，包括由全圆柱以及方柱与圆柱两种障碍物阵列，得到不同工 况下的流场；此外在全圆柱障碍物阵列条件下模拟了雷诺数为200，300，500，600，800时的流动情况， 并取近出口和近入口两个断面，得出断面流速曲线，通过观察曲线变化分析得出在多孔介质层的断面 速度变化情况与纯流体区域变化情况相反，结构复杂的障碍物阵列对于流体的阻力更大。 关键词：格子Boltzmann；通道流动；断面流速 中图分类号：TKl24文献标识码：A文章编号：1004—7948(2013)02—0015—03 doi：1039694．issn．1004—7948．2013．02．004 引言 方程的一种特殊离散形式，这种离散包括速度、时 BoltzmannMeth— 格子Boltzmann方法(Lattice 间和空间离散。模型如下j： od，LBM)为复杂现象的理论探讨和计算机模拟提 供了有效的手段，受到了各界的极大重视。在多孔 OI ‘，， 介质流等相关领域也得到了比较广泛而成功的应 厂9(工，t，V)] (1) 式中：厂4(x，t，1，)—平衡分布函数； J。 用。国内外学者对此进行了大量的研究工作¨。3 格子Boltzmann方法是一种建立在微观模型和细八工，t，v)一粒子的分布函数； 观运动论方法基础上的，宏观上离散，微观上连续 竹一松弛时间。 的直观模拟方法HJ。相对于一般的数值计算方 将式(1)离散处理为： 法，格子Boltzmann方法不需要建立复杂的偏微分 方程，它保留了分子动力学的优点，例如物理条件 清晰，边界条件处理简单，程序易于实施。 (2) f尸(工，f)]=一．EL(工，f)一f，(工，f)] 本文将使用格子Bolmmann方法来构造多个 式中：ei一离散速度； 障碍物组成的障碍物阵列，改变雷诺数和障碍物排 r一松弛因子。 列及形状，取出不同工况同一断面上的速度曲线， 选取分布形式为Maxwell分布的局部平衡分 分析改变参数后流场及断面流速变化规律，通过模 布函数的泰勒级数展开形式‘5|，有： 拟规则多孔介质的流动情况，为今后模拟真实的多 孔介质流动打下基础。 f l数学模型 对式(2)再运用Taylor展开，得 本文采用BGK近似的格子Boltzmann模型， 基金项目：国家自然科学基金资助项目(项目编号；国A(下一寺卜妒(a卿+a以)