气固两相流中气体湍流强度的研究

《气固两相流中气体湍流强度的研究》由会员分享，可在线阅读，更多相关《气固两相流中气体湍流强度的研究（4页珍藏版）》请在金锄头文库上搜索。

1、1998 年 1 月J OU RNAL O F XIAN J IAO TON G U N IV ERSI T YJ an . 1998气固两相流中气体湍流强度的研究魏进家姜培正胡春波(西安交通大学 ,710049 ,西安)摘要在气固两相流的流体相湍流模型中计入颗粒尾迹影响 ,采用低雷诺数模型 ,对竖直上升管中气固两相流动进行了数值模拟 ,较好地揭示了流体相湍流强度随颗粒大小 、体积浓度变化的规律.关键词气固两相流颗粒尾迹湍流强度中国图书资料分类法分类号O359Study on Ga s Turbulence Intensity in Ga s2Sol id Two2Pha se Fl o ws

2、W ei J i n j i aJ i a n g Pei z hen gH u Ch u nbo( Xian J iaoto ng U niversit y , 710049 , Xian)This paper perfo r med a numerical simulatio n fo r gas2solid t wo2p hase upflow in a verticalAbstractpipe wit h low2Reynolds number mo del by inco rpo rating t he effect s of wakes behind t he particles

3、intot he gas t urbulent mo del . The variatio n charader of gas t urbulence intensit y wit h t he size and co ncent ratio n of particles were well revealed.Key wordsgas2sol i d t w o2p hase f low st he p a rt icles w akesgas t u rbulence i n tensi t y气固两相流动广泛存在于自然界和许多工业部门 ,如含尘风机 、管道输送 、旋风分离器 、流化床等 ,

4、因而研究气固两相流动规律具有重要的意义 . 对气 固两相流动建立适当的流动模型进行数值模拟 ,以预测两相流动规律 ,已是两相流研究的一个热门课题 . 目前的两相流数值研究普遍采用双流体模型 ,该 模型中的流体相湍流模型是通过借鉴单相流体而推 得的 ,从而未能计入大颗粒尾迹的影响 ,使流体相湍流强度的数值预测值始终因颗粒存在而降低 ,这与 Tsuji 等1 所得的颗粒较大时流体相湍流强度反而 增强的实验结果不符. 因此 ,不计入大颗粒尾迹对湍 流强度的增强效应这一物理事实 ,是不能较正确地 预测流体相湍流强度的 . Go re 等2 通过总结分析大量实验数据 ,提出用颗粒直径与流体湍流涡长度尺度

5、的比值来判断流体相湍流强度是否因颗粒存在而增强 ,当该比值大于临界值 011 时 ,湍流强度增强 ,反之则减小 ; Het sro ni 3 则 提 出 当 颗 粒 雷 诺 数 大 于400 时 ,流体相湍流强度因颗粒尾迹而增强. 但这些都不能定量地表示湍流强度因颗粒存在的改变值.Yarin 等4 研究了颗粒尾迹对湍流强度的影响 , 得到了定量的表达式 ,但只是通过对实验数据的分析来研究由于颗粒尾迹而引起的湍流强度的经验关系式 ,实际上湍流强度还与流体平均速度梯度有关 ,因此未能全面考虑流体相的湍流强度. 本文则是将文献 4 的研究结果应用于流体相湍流模型中 ,并考虑颗粒对流体湍动能的耗散作用

6、 ,进行数值模拟 ,较全地研究了颗粒大小 、浓度高低对流体相湍流强度将 ( 5) 式代入 ( 4) 式得 Kf = 35 . 035(p d g) / (f CD) (p / f C3/ 2 ) 8/ 9影响规律 ,从而为完善流体相湍流模型打下基础 ,D气固两相流动数值模拟更接近于工程应用实际 .( 6)将该项加入 ( 2) 、( 3) 式得控制方程S K= G - -D + K( 7)ffk f f feffS= C1 ( Gk + Kf )- C2 (fff + De)f本文采用文献 5 中的两相双流体湍流模型 ,流KfKf相采用 Kf 2f 模型 , 颗粒相采用局部追随模式 , 颗( 8

7、)由于颗粒的存在使流体相在壁面的速度分布律 已大为改变 6 , 故本文采用低雷诺数模型 , 以解决 壁面函数不再适用的困难 . 本文借用了文献 7 中处 理单相流体的观点 , 在控制方程扩散系数中计入分子扩散系数 , 同时 ( 7) 、( 8) 式变为间的相互作用通过引入颗粒温度 Tp 方程来加以1虑 , Tp 表示单颗粒脉动的强弱 , Tp =u2 , 其中p3为单颗粒脉动速度.任意曲线坐标系下 , 稳态的气固两相流动方程表示为 5 S K1/ 22Gk - fff -De + Kf -2f ( 5 Kf / 5 y)( 9)=flm 5k 5 l5 X l kkJ ( U kk - g5

8、X m ) = J Skkf= C1 f 1 ( Gk + Kf ) - C2 f 2 (fff +S f( 1)Kf中 k 、k 分别为相 k 的体积浓度和密度 ; k = f 、pDe)f / Kf + 2ft ( 52 U f / 5 y 2 )( 10)别表示流体相和颗粒相 ; J 为 J aco b 行列式 ; glm为式 ( 9) 、( 10) 分别为计入了颗粒尾迹对流体相湍流强变度量张量 ; U l 为相 k 速度的逆变分量 ; 为度影响的流体相 Kf 、f 方程的源项. 其中右边最后kk一项为考虑壁面处湍流耗散各向异性的附加项 ;t 、的扩散系数 ; k 可分别代表 1 、U

9、ki ( i = 1 , 2 , 3) 、f 1 、f 2 的具体表达式同文献 7 .、f 及 Tp ; U ki 为相 k 速度的协变分量 ; Kf 为流体湍动能 ;f 为流体相湍动能耗散率 ; Tp 为颗粒相2数值解法度 ; S 为方程源项.k对于 Kf 方程 , 其源项 S K 的表达式为f8Spalding借用求解单相流体的 S IM PL E 算法的思想 ,提出了 IPSA 算法 ,该算法通过由压力修正 方程得到的压力修正值来进行速度和压力的修正.当用压力修正量修正速度时 ,可以使速度场满足连 续方程 ;但用于修正压力时 ,则会导致一个相当强烈的压力修正场 ,不易收敛 ,尤其对于两相

10、流场 ,由于相间作用力的非线性耦合作用导致各求解变量源项 的复杂化 ,从而使收敛性更差. 为了克服这一缺点 , 本文借用 S IM PL ER 算法的思想 ,将压力修正量仅用于修正速度场 ,而压力场则通过求解压力修正方 程得到. 为此 ,本文算法迭代步骤为 :(1) 由颗粒浓度方程计算颗粒相体积浓度 ,再由 总守恒方程计算流体相体积浓度 ;(2) 由流体相压力方程计算流体相压力场 ,根据 气固两相间压力关系计算颗粒相压力 ;(3) 根据气固两相压力场 ,计算两相动量方程 ,得出两相速度场 ;(4) 计算气固两相浓度加权压力修正方程 ,得出S K= Gk - fff - De( 2)f中Gk 为

11、湍动能产生项 ; De 为因颗粒与流体湍作用而引起的流体相湍流强度的耗散.对于 f 方程 , 其源项 S的表达式为fffS- C2 (fff + De)= C1 Gk( 3)fKKff中C1 、C2 为经验常数.Yarin 等的研究结果表明 , 由颗 粒 尾 迹 使 流 体湍流强度的增加量为 Kf = 26 . 275(p / f C3/ 2 ) 8/ 9 ( u-u) 2Df p( 4)中 为经验系数 ; CD 为两相间阻力系数 ; u、f分别为流体和颗粒相脉动速度.由于 当 颗 粒 与 流 体 的 密 度 比 值 大 于 100 时 ,f 主要与重力有关 , 由文献 4 得= 4 p d

12、g( u-u) 2( 5)f p3 f CD(5) 由气固两相浓度加权压力修正值修正两相的速度场 ;(6) 计算流体相 Kf2f 双方程及颗粒相温度方 程 ,求出流体湍流扩散系数 、颗粒相湍流扩散系数和 层流粘性系数 ;(7) 重复以上过程直至达到收敛指标.边界条件处理方式如下.(1) 进口边界 : 给定流体相 、颗粒相速度分布及 相应的体积浓度分布 ( 本文取均布) ,流体相进口湍 动能按进口动能的 1 %给出 ,进口 f 按下式给出动均方值. 由图可见 ,随颗粒浓度增加 ,流体湍流度降低 ,在管道中心区出现了极端扁平现象. 考虑颗粒尾迹影响所得的结果更接近实验值 ,因此 ,即使小颗粒情况也

13、应加入颗粒尾迹影响. 另外 ,本文为了探讨密相情况下颗粒对流体湍流度的影响 ,对实验没测试的 2 %浓度情况下的流体相湍流强度也进行了计算 .C3/ 43/ 2 / 0 . 03 lfin=Kf( 11)其中l 为进口特征长度.进口 Tp 为9= Cpc l( 12)Tpin其中C 为常数 ;pc为颗粒相进口平均体积浓度.(2) 出口条件 :按局部单向化处理.(3) 壁 面 条 件 : 流 体 相 速 度 满 足 无 滑 移 条 件 ,Kf 、f 均为 0 ,颗粒相速度及温度同文献 10 .(4) 轴对称条件 :5k5 n( 13)| sym= 0计算结果与分析3文献 1 对竖直上升管中的气固

14、两相流动进行了实验研究. 实验圆管内径 D = 3015 mm ,介质为空 气和聚苯乙烯球 ,气流雷诺数为 23 000 ,颗粒密度为1 020 kg/ m3 ,体积浓度从 0 %014 %变化.本文对文献 1中竖直上升管的气固两相流动 进行了数值计算 ,认为颗粒尾迹较明显并出现流体湍流强度增强现象的颗粒为大颗粒 ,反之则为小颗粒 . 因此 ,本计算选取平均直径 012 mm 和 015 mm的聚苯乙烯球分别作为小颗粒和大颗粒的典型代表 进行研究 ,取管道轴对称的一半 ,长度为 2 m ,以满足充分发展的条件 . 划分网格为 50 80 , 其中粘性 支层区均布 40 个节点 ,粘性支层厚度取法是先利用单相流体经验公式算出一粘性支层厚度值 ,然后将 其放大 112115 倍 . 计算表明 ,该方法能有效地布置壁面节点 ,进