稠密气固两相流动数值模拟的最新进展.pdfVIP

中国工程热物理学会 编号：986007 多相流学术会议 稠密气固两相流动数值模拟的最新进展 祁海鹰李字红徐翘常 (清华大学热能工程系，北京100084) 摘 要 本文对气固两相流动的“双流体”假设和在此基础上发展而成的欧拉(Euler)数值模拟方法加以 介绍，阐述了在稠密条件下所应用的理论和部分物理模型的最新进展，包括描述固体颗粒相的离散质 动力理论、边界条件和气固夹带阻力模型．文中还给出了部分数值模按结果． 关键词：双流体，欧拉方法，离散介质动力理论。边界条件，气圈阻力 一、引言 稠密多相流动是普遍存在的物理现象，其中以燃煤循环流化床内的稠密气固两 相湍流流动最具代表性．此外，多相流研究在探讨诸如管道气力输送、旋风分离器 近壁浓相区特性、河流泥沙沉积规律、粮食堆放及选煤过程、燃煤烟气脱硫除尘以 及沙漠变迁都显示出广阔的和潜在的应用前景． 较之稀疏流而言，稠密气固两相流的特殊性在于：一，颗粒相的体积百分比浓 度高，故它对流动的影响已占有相当的比重，而不能加以忽略，它在流体中的运动 也不能再简化地视作完全对流体的跟随．二、颗粒之间的相互作用大大强化，并产 生团聚现象，这同样是不能加以忽略的．为了合理描述稠密气固两相流动。 “双流 体”的概念和假设应运而生c】】． 所谓“双流体”是指由气相和被处理成“拟流体”的固体颗粒相组成的气固两 相流．在一定的浓度条件下，可将颗粒与气相质点相比拟，而将其看成是如气相一 样的连续介质，使之具备连续介质的动力学特性，从而能够通过形式上与气相相同 的流体力学守恒方程描述离颗粒这种离散介质．运用“双流体”假设意味着从整体 上观察颗粒相，因此，由此发展而成的“歇拉”数值模拟方法就能够模拟稠密的气 固两相流动．这里所谓的“欧拉”方法不是指一般流体力学中不考虑粘性的欧拉方 程，而是指在多相流数值模拟中，与追踪单颗粒轨迹的拉格朗日(Lagrange)方法 相对的、对颗粒相作宏观的空间观察的一种方法．该方法计及颗粒的体积，通过丽 相各自的体积百分比浓度对当地的诸流体力学参数进行空间上的修正，并耦合气固 两相的相互作用．其它各种影响因素的作用以及湍流效应都能够通过模型封闭在这 种方法中得以进一步的考虑．因此，与拉格朗日方法相比，欧拉方法在研究稠密气 固两相流问题上具有独特的优越性。自s．LSodll于1967年提出“双流体”的概念 之后，M．Ishii 化床【4J． 在稠密条件下，颗粒的离散介质特性，即由颗粒碰撞所引发的一系列颗粒内部 应力不容忽视，为此，作为应用欧拉方法的重点之一，人们着重对颗粒的脉动特性 迸行了深入研究，发展了“离散介质动力理论BP，使颗粒特性的模型描述更加完 善合理。过去，人们主要采用“固体弹性模量”(solid elasticitymodulus)G(￡0来描 写颗粒应力【6】．它是通过变换颗粒应力张量而导出的称之为“颗粒压力”的一个经 Ⅱ一58 验表达式，并被用在动量方程中．如今，它已被离散介质动力理论所取代． 在稠密条件下的颗粒相边界条件变得比稀疏流动要复杂得多．这是因为颗粒内 部应力及其与壁面的碰撞引发一系列能量传输．在壁面处的能量生成与耗散决定了 必须采用不同的边界条件，否则，数值模拟结果会受到直接的影响． 鉴于循环流化床内流动的特殊性和复杂性，描写气固间相互作用、耦合两相微 分方程的气固夹带阻力模型成为数值模拟中的关键环节．研究表明，只有考虑了颗 粒的团聚效应，才可能实现循环其气固两相流的数值模拟并获得合理的结果[7】． =、离散介质动力理论(ThekineticTheo吖GranularMedien) 式(1)和式(2)中分别是描写气固两相流的质量和动量通用守恒方程： 警+毒(帆)；峨 (1) 擀uk