小雷诺数下颗粒所受Saffman力地研究.pdfVIP

中国工程热物理学会 多相流 第十届年会 编号：016036 小雷诺数下颗粒所受Saffman力的研究 由长福，祁海鹰，徐旭常 清华大学煤的清洁燃烧技术国家重点实验室．北京100084 电话：010- E阳il：youcf@te．tsinghua．edu．cn 摘要：本文采用数值模拟方法对剪切气流中静止颚粒所受到的Saff咖力进行了数值摸拟．控制 方程采用三维稳态N-S方程，有限差分法毗及SII[PI上C算法用于对控制方程的求解．计算结果表明： 颖粒跗蚯的同格尺寸对颗粒受力有较大的影响；当雷诺数小于100时．Saff∞an升力系数与曳力系数 关键词：颗粒受力；两相流动；7升力； 舭萋耋：：颗……赫‰Saffn。n力蒜拟力；’数值模拟 l、引言 颗粒所受各种流体作用力的研究对于准确描述颗粒在气固两相流中的运动具有重要 的意义。一般来讲．这些力主要包括纵向力(如曳力、Basset力、附加值量力)．横向 力在小雷诺数条件下进行了详细的研究．得到了较准确的颗粒受力表达式。本文采用数 值模拟的方法对SMTmaa力在小雷诺数条件下进行研究。 在剪切流场内颗粒所受到的横向升力称为Sa伍n∞力。s绷缸m力的研究始于1965 年s曲_man，EG．例采用角摄动法(singularpe巾l山舶n)得到的描述剪切气流中颗粒所受升 力的表达形式： 刖w扭即料嘲 (1) 在气固两相流动的数值计算中，当考虑Saffman力的影响时，式(】)被广泛地采用。 而其应用条件则被许多研究工作忽略。事实上，式(1)的成立存在严格的条件．即： (地=p一砷 Re，=譬“L ‰=华． (。=踟 ‰=掣 (2) Re．“醚1 为准确描述SMIman力对颗粒运动的影响，在式(2)所规定的条件之外．一些研究人 员也进行了相应的实验与理论研究工作。如最近发表的文献【4，51总结了有关单个颗粒在 剪切气流中的研究并采用数值模拟的方法对不同雷诺数、剪切率条件下的颗粒所受到的 1∞ 升力进行了研究。这些研究均发现，当流动雷诺数逐渐增大时．在某一雷诺数条件下， Saffman力受力方向发生改变：文献[4】的研究结果为Re=50，而文献【5]的研究结果为 Re=60。需要说明的是，以往对Saffman所进行的数值模拟研究多采用容易处理复杂边 界条件的有限元法a而极少有研究者采用有限差分方法对剪切气流中颗粒的受力进行分 析。 本文采用数值模拟方法，对Re100条件下颗粒在剪切气流中所受到的sa丘mⅫ力 进行研究，采用有限差分法对控制方法进行离散，并采用SIMPLE(]算法进行求解。其 中还详细考察了网格尺寸对计算结果的影响，同时对比了s柚陆柚升力系数与曳力系数。 2、数值计算 在许多实际工程问题中，如煤粉颗粒的燃烧，其颗粒雷诺数均小于100。本文对剪 切流动中颗粒所受sa缸an力的研究也在此范围内进行。对此流动进行数值模拟时，控 制方程如下： V·厅=0 甲．(厕：-Vp+卢Vt疗 (3) 本文采用有限差分方法对控制方程进行空间离散，并采用smⅡ，u蛇算法解决速度 与压力的耦合计算问题。计算域如图1所示： 圈I计算域虬及边界条件 计算中进，出口边界条件取为周期性边界务件，即仅仅给定流体的质量流量，而无需给 定速度边界条件。流过颗粒的剪切流动通过在计算域上下表面设为移动边界条件来实 现。通过控制图l中L1与L2的大小使得进E1处速度分布为线性分布。计算中颗粒直