颗粒绕流的格子boltzmann法模拟.pdfVIP

中国T程热物理学会 燃烧学 学术会议论文 编’寸：074158 11tzmann法模拟 颗粒绕流的格子bo 张健 刘向军 北京科技人学热能工程系，100083，北京 术电话：0l0 E-mail：Liu蜀@me．ustb．edu．cn 摘要本文采用格I了．b01tzmann方法对气体绕流颗粒的流动与受力规律进行了研究，计算了颗粒粒径为 100um，水流速度为0．-t5ln／s时双颗粒，多颗粒绕流的流动状况，分析I：：r颗粒周嘲流体的速度场，腻强场的 变化及颗粒的受力情况，得i{=“各个颗粒详细的受力大小及{i!{!律，为分析稠密气周两相流中颗粒的复杂崭e动现 象提供__r的理论甚础。 关键词 格了boltz魄nn：稠密气固两相流：颗粒绕流 、 1．前言 气I矧两相流广泛存在下．=f：程实际中，在化’l：及能源动力等领域有着非常J1‘泛的应用。 气矧两相反应，尤其是稠密气同两相反应流内颗粒一颗粒，颗粒一气体的相互作用是十分剧 烈、复杂的。恒温恒压。卜．的自发过程总是向降低自由能的方向进行。气吲丽相流中同体不能 像液体那样可以改变自己的形状来降低表面积，唯一的途径就是与周围的颗粒相互粘附减小 比表面积达剑较为稳定的状态。冈此，细颗粒存在的颗粒体系中总有粘附现象发生…。颗粒 团是稠密气俐两相流的一个主要的基本特征。此时颗粒团中颗粒所受气动力与单颗粒状态时 的受力特征有了明显的改变，分析研究多颗粒及颗粒团详细的绕流规律，对稠密气俐两相流 动的深入研究具有重要意义。 数值研究颗粒的绕流规律目前主要还是采川传统的基1iN—S方程的数值模拟。而最近发 展起来的格子boltzmann方法，有着传统数值模拟无法比拟的优势Ⅲ‘： (1)边界条件易丁．处理； (2)群序编写简单，适合丁大规模的并行运算； (3)能很容易地应刚于微尺度流体流动问题上。 u 本文采用格子boltzmann方法，对颗粒粒径为100m，米流速度为0．5m／s时双颗粒， 多颗粒绕流的流动状况进行了计算，分析了颗粒周同流体的速度场，压力场的变化及颗粒的 受办睛况，得出各个颗粒详细的受力人小及规律，为分析稠密气蚓两相流中颗粒的复药}流动 现象提供了的理论基础。 2．D2Q9模型 格子boftzma彻方法是一种不同丁：传统数值方法的计算方法，它源丁格子气自动机。 其理论整础是，流体由流体分f构成，流体的宏观运动是分子微观运动的平均结果，宏观动 基金项目：圈咏自然科学綦会顾日 Natj()naJNnturalScienceFoundation I)r【)lccI ofChina(NSFC) supportcdb 460 力学特征对流体微观运动的细节并不敏感。它的基本思想是构造一个简化的动力学模型，、该 模型能够反映流体运动的微观的或细观的物理本质，并且在宏观上遵循流体运动的规律。 本文采刚D2Q9模型。’I。旧1，其分布函数Z的动力方程如下： f：o，l，…’8 (1) ／：(x+万；，，f+万)一Z(x，f)：一三【Z(x，f)一，钾(x，f)】， 9速方格BGK局部平衡分布函数，钾具有如。f：_．形式： ．f州=w，p[1+3(虿-二)+兰(虿·；)2一吾磊·云】， (2) ；，=压(c。s叭f一4一i)／2)sin㈣f一4一：)／2))，，=5，67’8