粒子在锥形管中运动的晶格玻尔兹曼方法研究.pdf

物理学报 Acta Phys. Sin. Vol. 62, No. 8 (2013) 084702 粒子在锥形管中运动的晶格玻尔兹曼方法研究* 1 1 1 1 2 1† 薛泽 施娟 王立龙 周锦阳 谭惠丽 李华兵 1) ( 桂林电子科技大学材料科学与工程学院, 桂林 541004 ) 2) ( 广西师范大学物理科学与技术学院, 桂林 541004 ) ( 2012年9月28 日收到; 2012年12月22 日收到修改稿) 运用晶格玻尔兹曼方法对单个悬浮粒子在锥形管中的运动进行了数值计算, 给出了锥形管流体的速度分布和 压力分布等. 粒子所受的流体作用力分别用动量交换法、改进的动量交换法和压力张量积分法进行计算. 分析了在 不同初始位置释放的粒子的运动轨迹和速度变化情况, 结果表明压力张量积分法和改进的动量交换法的计算结果 一致, 而没有改进的动量交换法的计算结果和前两者略有不同. 关键词: 晶格玻尔兹曼方法, 锥形管, 悬浮粒子, 改进的动量交换法 PACS: 47.11.−j DOI: 10.7498/aps.62.084702 可以忽略它的影响, 这就是所谓锥度角效应. 有必 1 引言 要对有锥度的血管的流动进行深入的研究1011 . 本文主要基于LBM, 建立了以黏性不可压缩 微小粒子悬浮运动现象广泛存在于自然界和 牛顿流体在轴对称二维非弹性锥形管的流动模型, 工业生产中, 如花粉粒子的输运与传播、血管中 以单个刚性粒子在锥形管中的运动来模拟血液在 红细胞随血浆的运动、废水中带粒流的分离与处 小血管中的运动(这里红细胞是不可形变的粒子), 理等. 了解悬浮体的运动规律对于改善生态系统 获得了流体的速度分布和压力分布等. 流体对粒子 的循环、了解血液流变的奥秘、提高化学工业的 12 边界的作用力分别采用动量交换法 、压力张量 生产效率有重要的意义. 晶格玻尔兹曼方法(lattice 13 14 积分法 和改进的动量交换法 来计算, 分析并 Boltzmann method, LBM) 是在格子气(lattice gas) 模 比较了三种方法在不同初始位置释放的粒子的轨 型的基础上发展起来的一种新的探索复杂系统的 迹和速度变化情况. 数值计算方法和模拟方法. LBM 已成为目前流体 系统的研究中最有发展前景的方法之一, 被广泛用 2 D2Q9 LBM 于研究各种复杂流体的运动1−4 , 以及流体与粒子 的相互作用5−7 . 由于LBM 中流体格点间的相互 为简单起见, 我们采用二维九速四方格子 作用是局域的, 特别适合动边界和变形边界问题的 (D2Q9) 的LBM 模型来计算粒子在锥形管中的 研究, 而血液流系统是有代表性的边界特别复杂的 运动. 流体格点上分布函数的演化方程为15 生物系统. 近几十年来, 在国内外见之于公开文献 1 eq