喷动床内颗粒流动特性的数值模拟研讨.pdf

喷动床内颗粒流动特性的数值模拟研究 何玉荣*，詹文博，门玉宾，刘文铁，陆慧林 哈尔滨工业大学，能源科学与工程学院，黑龙江 哈尔滨150001 关键词：喷动床，欧拉-欧拉模型，摩擦应力，高颗粒浓度 NUMERICAL SIMULATION OF BEHAVIOR OF PARTICLES IN SPOUTED BEDS HE Yurong, ZHAN Wenbo, MEN Yubin, LIU Wentie and LU Huilin ( School of Energy Science and Engineering , Harbin Institute of Technology , Harbin 150001 , Heilongjiang , China) Abstract An improved Euler-Euler two-fluid model is applied in this work to predict dense gas-solid flow in spouted beds. Constitutive equations describing the particle interactions and friction of particles and viscosity were incorporated into a hydrodynamic simulation computer program. A body fitted coordinate is used in the simulation to make the boundary fit well with the real object. Three zones of spouted beds: spout, annulus and fountain were numerically studied. Numerical results are compared with experimental results of He et al. [1, 2] with good agreement. Keywords spouted bed, Euler-Euler two-fluid model, frictional stress, dense gas-solid flow 前言 喷动床被广泛的应用于燃料燃烧、气化、物料颗粒干燥及催化聚合等工业过程 中。随着喷动床在工业过程中应用范围的不断扩大和其优点的不断体现，喷动床内 的流动行为的研究愈显重要。 国内外学者对喷动床内气固两相运动进行了部分数值模拟研究工作[3-7] 。在上述 喷动床气固两相流动模型中忽略了颗粒间半接触状态下相互作用力，大部分研究者 没有考虑倒锥体内颗粒和气体的流动，仅考虑圆柱体内气体和颗粒的流动，在数值 模拟中多采用直角坐标系、用折线代替锥体倾斜面。 本文提出气固两相流动模型预测喷动床内流体动力行为。模型中引入了颗粒相 应力动力-摩擦封闭模型。动力模型基于颗粒动理学理论。应用Johnson 等[8-9]提出的 摩擦应力模型和由Syamlal 等[10]提出的修正摩擦剪切粘度模型，研究颗粒相互作用， 数值分析喷动床内气固两相流动。同时应用贴体坐标系以使网格很好的符合喷动床 倒锥体边界。 1.数学模型 1.1 质量守恒方程 ∂ (ε ρ ) + ∇ ⋅(ε ρ v ) 0 (1) i i i i i ∂t 其中，i 表示气体或者颗粒相，ε 是气体或颗粒浓度, v 是气体或颗粒的速度矢量, ρ 是气体或颗粒的密度。 1.2 动量守恒方程 气体动量守恒方程可表示为： ∂ (ε ρ v ) + ∇ ⋅(ε ρ v v ) ε ∇ ⋅τ + ε ρ g − ε ∇p − β (v −