高浓度颗粒气固两相流动的二阶矩模型的数值模拟研究.pdfVIP

中国工程热物理学会 多相流 第十一届学术会议 编号：056012 高浓度颗粒气固两相流动的二阶矩模型 的数值模拟 陆慧林，何玉荣，刘欢鹏，刘文铁，李响 哈尔滨工业大学能源科学与丁程学院，哈尔滨，15000l E·mail=h尘曲鼬i!：鲤丛：￡n，Tcl：0451 摘要：被广泛应用于气蚓两相流动研究的颗粒动理学是借鉴非均匀稠密气体分子运动 论，将单颗粒随机运动比拟为气体分子热运动，假定颗粒速度分布服从各向同性 M“well分布，建立颗粒相输运方程，以及与颗粒温度相关联的颗粒相物性参数计算模 型。然而，颗粒相互碰撞造成能量的耗散，导致颗粒速度分布呈现各向异性。采用动 程，分别取零次矩、～次矩和二次矩得到颗粒相连续方程、动量方程和二阶矩方程。 模拟计算得到上升管内颗粒流场分布和脉动速度分布，并与实验结果进行比较。结果 表明各向异性颗粒相湍流模型能全面反映颗粒相Reynolds应力、湍流动能等湍流特征 量，表明各向异性颗粒动理学能够反映颗粒运动的微观特性，从而能够更好地模拟和 认识气固两相流动。 关键词：各向异性，湍流模型，二阶矩模型，颗粒动理学。 一、前言 高颗粒浓度气固两相流动广泛存在各种工业应用领域。如化工与石油丁业中的裂 化催化反应器、洁净煤燃烧技术中的流化床燃烧技术，等。在高颗粒浓度气固两相流 动研究中，颗粒动理学方法被广泛应用于颗粒流动特性的预测如提升管反应器、流化 床锅炉等。但是在高颗粒浓度气固两相流动中颗粒呈较强的脉动流动特性。颗粒脉动 往往在某个主导方向上最强，而在其他方向上较弱，颗粒流动出现强烈的非均匀性和 各向异性。颗粒动理学是基于分子动理学理论，将固体颗粒类比为气体分子，提出反 映颗粒脉动大小的“颗粒温度”概念，即日=CC，3，其中，c为颗粒脉动速度， 并借助颗粒速度分布函数的BoItzmallll微积分方程推导出颗粒压力、颗粒粘度等物性 特征量的表达式。由分子动理学可知，由M“we¨速度分布函数获得的气体各速度分 量的分布函数具有对称性，表明气体处于平衡态时，气体的性质与方向无关，分子向 各个方向运动的几率均等。所以分子在三个速度分量的方差必然相等。由此可见颗粒 动理学方法是以各向同性为基础，应用M“wen速度分布函数，确定颗粒相输运参 数。因此在使用颗粒动理学方法模拟高颗粒浓度气固两相流动内各向异性的颗粒湍流 时有时会出现较火的误差 二、气固两相流动模型及颗粒相矩方法 气固两相流动应满足质量、动量和能量守恒。对于气相和颗粒相的质最和动量守 恒方程分别如F。 2．1气相守恒方程 气相连续性方程： 昙k以)+vk掣。)=o 基金项目：国家自然科学基金瓷助项目(No 188 气相动量方程 (2) 勺V名+％艮g+V-k一声09一％) 昙k以％)+V如。艮％‰)= 五堪gPgHgpVugpgHgHg卜 式中，名为气相压力，g为重力加速度， ∥为气I古】相问曳力系数，79为气相应力张 姑。 fg=∥g[V“占+(V“占)’卜詈∥g(V‘Ⅳg)I (3) 2．2颗粒相守恒方程 假设，为颗粒速度分布函数，它是时间r、速度Ⅳ和坐标，的函数，且满足统一微 分方程 鲁+c善+磋=c考k“西 加 ＆、西…“ ca，～ 方群式(4)表示了颗粒碰撞引起的相空间体元(△，△c)中的颗粒数的净增率 【钟／a)，△f△，△c等于相空间中对流引起的颗粒数净流出与相空间体元(△，△c)中的 颗粒数的净增量之和。假设颗粒是二体碰撞。则因碰撞而导致的生成源项可表示为