颗粒轨道模型中曳力系数影响.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号：036048 颗粒轨道模型中曳力系数的影响+ 欧阳洁” 崔俊芝 西北工业大学应用数学系，西安，710072 中科院计算数学与科学工程计算所，北京，100080 李静海 孙国刚 中科院过程工程研究所，北京，100080 北京石油大学化工学院，北京， 100083 摘要本文在两相流模拟的颗粒轨道模型中，采用不周的曳力计算公式对密相气固流动系统进行了模 拟。其研究结果表明：对颗粒所受曳力的计算应该以颗粒群的曳力计算为基点，且基于Wen和Yu以 及DiFelice的曳力公式能够较好地模拟密相气周两相流的基本流动特征。而直接从双流体模型中的 气相运动方程出发，对颗粒所受曳力进行计算是否完全合理尚需进一步研究。 关键词两相流，颗粒轨道模型，模拟，曳力 1 引言 近年来，在气固两相流的模拟中，直接对颗粒进行跟踪的颗粒轨道模型由于计算 机硬件的发展而得到关注【”I【6I。这类模型将流体处理为连续介质，将固体处理为离散 颗粒。从理论上说，若不考虑颗粒数目的限制及计算中的误差，颗粒轨道模型对气固 流动系统描述的真实性强于多年来处于主导地位的双流体模型。由于这类模型可以从 介观尺度描述颗粒运动，并可用于认识固相的离散特征。因此对于气固两相流的模拟 而言，颗粒轨道模型是很有前景的模型。 在气固两相流中，颗粒在受气体作用的悬浮过程中，其运动规律可用牛顿第二定 律给出： 碾鲁=(∑n 这儿(’F)．表示颗粒i所受的合力。 作用在颗粒上的合力相当复杂，但在实际计算中，通常忽略影响较小的作用力”J。 在所有相间力中，曳力相对而言是最重要的力。目前对单个颗粒在无界流场中所受曳 力(即曳力系数)的研究已非常充分。但在实际多相流中，颗粒群曳力的大小还受到 许多因素的影响。因此，在不同的情况下，需要对单颗粒的曳力系数进行修正。目前 大部分研究工作是归纳得到各种操作条件下的经验关联式，而不同研究者的关联式计 算结果相差较大，目前尚无很可靠的计算公式【7】。而曳力系数的不同，必然对模拟结果 有影响。本文将针对气固密相流，观察常用曳力系数关联式对模拟结果的影响。 2基本数学模型 散的颗粒相。 · 国家自然科学基金资助项目，陕西省自然科学基金资助项目，开放实验室资助项目 本工作为联系人在中国科学院数学与系统科学研究院访问期间完成。 273 型．t-V．6pgu=08t 丁O(epgu)+V．(瞩Ⅳ砷=一卿一昂+V．(曝)+峨占 其中幽空隙率·风为气体密度，口为气体速度，P为压力·S。为动量交换源项，g为 重力加速度，气相粘性应力为： Tg=(一鲁以)(V，H)n以((VⅣ)+(V甜)7) 2．2颗粒运动方程 在两次碰撞过程之间，颗粒i的运动规律由牛顿动力学方程给出： du i=l，·一，N 7q=÷二=m,g+(乃)，一(0)，VP， df 其中，mI为颗粒质量，H为颗粒速度，(日)。为颗粒所受曳力，(匕)。为颗粒体积，Ⅳ 为系统中的颗粒数。 在稠密气囿两相流中，颗粒问的相互作用不能忽视。目前处理颗粒间相互作用的 模型有两类。第一类是软球模型IlⅢ，这类模型允许多颗粒进行非瞬态接触。第二类是 硬球模型Ij卜16j。这类模型假定颗粒间的碰撞是二体瞬时碰撞。相对硬球模型来说，软 球模型能够处理静态及多颗粒的相互作用。但在软球模型用于气固两相流模拟的研究 中，常常人为地选取小的刚性系数。以减少计算时间。在这些计算中，总是认为人工 调整颗粒的刚性系数是允许的，但从未证明它不会影响模拟结果¨