颗粒数量平衡方程的直接矩积分求解.pdf

中国工程热物理学会 多相流 学术会议论文 编号：066094 颗粒数量平衡方程的直接矩积分求解 苏军伟顾兆林1李云 西安交通大学动力工程多相流国家重点实验宣．710049，西安1联系人．通讯地址：西安交通大学 人居学院环境科学与技术系电话：029E-mail：g筵h艘!娶@理巫墨j丛￡垂婪￡堕 摘要：双流体模型将各相看成连续的相．要想通过该模型描述颗粒微观行为t如颗粒的聚并和破碎等， 必须引入数量平衡方程。关于数量平衡方程的数值求解方法有三类：分组方法(ClassMethod)．Monte 的直接矩积分方法(DOMOM)，并对各向同性颗粒系统的颗粒长大及其颗粒的聚并和破碎等微观现象进 行模拟，验证直接矩积分方法的合理性． 关键词：直接矩积分方法：数量平衡方程；颗粒系统 1引言 工业上的很多处理过程(比如：制药，结晶)都有颗粒参与，目前描述气固相互作 用的多流体模型将各相看成连续的耜，很难微观描述实际操作过程中粒子的动态行为 (比如颗粒合并和破碎，粒子的长大)，因此需要引入新的模型——数量平衡模型。数 量平衡可以看成跟踪任何实体的数量，这类实体可以是固体颗粒、液滴、气泡或单元等， 更抽象地来说．是能够描述被研究系统行为的任何存在或者发生的事件。在颗粒系统中． 颗粒的分布不仅和空间位置有关系还与颗粒本身的属性有关，方程呈现内部属性坐标和 外部空间时间坐标的双重特性【1]。由于其内在的复杂性，目前，对于～般化的数量平 衡方程得不到其分析解。就其数值求解方法而言，主要分为三类：分组方法(Class №thod)、Monte of Canto方法(MCM)和矩方法(MethodMoment)。分组方法基于对颗 粒的内部属性进行分组，得到一系列离散的偏微分方程，然后对偏微分方程进行求解。 要得到满足工程需求的解必须将数量平衡方程分成很多组(20’30)，计算量相当大，该 方法很难和CFD耦合进行实际的工程模拟[2]。Monto Canto方法基于统计的思想，为了 减少统计误差，需要大量的统计数据，如果考虑到颗粒分布的空间异性的话，对一个真 537 实的系统进行模拟该方法是不可取的[3]。MOM并不跟踪颗粒分布的变化，而是直接跟踪 颗粒的矩的变化，该方法计算量小．能够很好的和CFD进行耦合．但是由于标准矩方法 的封闭问题，限制了其使用，McGraw对颗粒分布引入了高斯积分从而使该方法得以封闭， 近几年该方法得到了广泛应用[3—5]，并且该方法在已经向多变量属性坐标扩展[6-8]． 基于矩积分方法Marchisio＆＆Fox已经导出了直接矩积分方法(DirectQuadrature MethodOf 得到了广泛的应用[9-11]。 目前国内化工领域对数量平衡方程的研究相对较少，而且主要是针对汽液两相流 动，采用的分组方法，基于矩方法的气固两相流动的数量平衡方程的研究缺乏。 本文引入数量平衡方程的一般形式，对数量平平衡方程的标准矩方法作简单介绍， 引出直接矩积分方法，最后对各向同性颗粒系统的颗粒长大及其颗粒的聚并和破碎等微 观现象进行了模拟，得到的结果与分析解以及文献中的结果比较吻合。说明直接矩积分 方法的特点。 2数量平衡方程及其求解 2．1数量平衡方程 气固多分散系统通常采用颗粒分布函数月色“jt∥来描述，其中￡表示颗粒的特 征长度。u。表示颗粒的速度。该函数表示某时刻t，空间位置为x，特征长度为L，特征 速度为tI．的颗粒的数目，数量平衡方程用来描述分布函数的随上述两套坐标(内部属性 坐标和外部时间和空间坐标)的变化关系，对于气固系统的单属性数量平衡方程可以用 下面来描述’ (1) 兰!!!!!：：：姜三兰王!+V．[“，n(上，虬；x，f)】+V也·【Fn(L，虬；x，f)】=S(L，“；；x，t) 其中，X表示空间坐标：t表示时间：S(L乩j鼻0代表源项，与颗粒系统的微观行