水平管入口段内水煤浆流动特性数值模拟

第40卷第2期东南火学学报(自然科学版)Vol 40 No. 22010年3月JOURNAL OF SOUTHEAST UNIVERSITY(Natural Science Edition)Mar.2010doi:10.3969/jiss.1001-05052010.02.036水平管入口段内水煤浆流动特性数值模拟陈良勇12段钰锋!刘猛!蒲文灏'赵长遂(东南大学能源与环境学院南京210096)(2南京工业大学能源学院南京210000摘要:以 Eulerian多相流方法为基础建立了水煤浆流动的液固两相流体动力学模型,并对水平管入口段内水煤浆流动进行了模拟.模型采用颗粒动理学理论求解固相本构方程,釆用RNGk-ε湍流混合模型描述颗粒间具有强烈作用的两相湍流流动.针对水煤浆中煤粉顆粒的双峰分布特性,将煤粉看作2种大小不同的固相,同时考虑固相与液相、固相与固相之间的动量交换.模型有效性通过 Kaushal等试验结果和水煤浆的压降试验验证.通过模拟考察了入口段内速度及浓度分布过程和入口段长度的变化规律.结果表明:重力和颗粒间的相互作用对浓度和速度分布过程具有重要影响;固相体积分数在30%~49.5%范围内,入口段长度随浓度的增加而减小,随平均流速(0.2~5.0m/s)先增加后减小关键词:水煤浆;入口段;多相流模型;颗粒动理学中图分类号:O373文献标志码:A文章编号:1001-0505(2010)02040207Numerical study of developing coal-water slurry flowin entrance region of horizontal pipeChen Liangyong, 2 Duan Yufeng Liu Meng Pu Wenhao Zhao Changsui( School of Energy and Environment, Southeast University, Nanjing 210096, China)(College of Energy, Nanjing University of Technology, Nanjing 210009, ChinaAbstract: An Eulerian multi-fluid flow model is developed to carry out a numerical study on coal-water slurry flow in the entrance region of horizontal pipe. The kinetic theory of granular flow ised to calculate constitutive equations of the solid-phase components and the rng k-e turbulentmodel is incorporated into the governing equation to model turbulent two phase flow with strong par-ticle-particle interactions. In this model, the coal particles with bimodal distribution are consideredtwo solid-phase components and the moment exchange between solid and liquid phase as well as thatbetween solid and solid phase are taken into account. The model was validated with Kaushal expenibutions in the entrance region as well as the entry length. The results show that the gravitationalforce and the strong particle-particle interaction have significant effects on concentration and velocitydistribution profile. The entry length decreases as solid concentration increases from 30% to 49. 5%and for the same slurry, the entry length firstly increases and then decreases when the inlet velocityincreases from 0.2 to 5.0 msKey words: coal water slurry; entrance region; multi- flI: kinetic theory of granular flow中国煤化工收稿日期:200909-30.作者简介:陈良勇(1977—),男博土生;段钰锋(CNMHGuan@seu.edu.cn基金项目:国家重点基础研究发展计划(973计划)资助项目(2010cB227001)引文格式:陈良勇段钰锋刘猛等.水平管人口段内水煤浆流动特性数值模拟[J]东南大学学报:自然科学版,2010,40(2):402-408第2期陈良勇,等:水平管入口段内水煤浆流动特性数值模拟水煤浆(CWs)管道输送技术是实现大规模高双颗粒组分浆体流动的试验结果21,但所提供效煤气化技术或CWS燃烧技术集成的关键组成的流动信息仅限于管道阻力损失.另外,尚无文献部分之一.目前,由于测量技术的限制和液-固两相报道双颗粒成分浆体在管道入口段内流动规律基流动的复杂性,除阻力损失外,对CwS这类高浓于以上分析,本文以 Eulerian方法和KTGF理论为度浆体的流动测量十分困难流体动力学建模方法基础建立CwS流动的液-固两相流体动力学模是获得浆体内部全面流动信息的有效替代方法,可型,采用RNGkε湍流模型描述液-固两相湍流流为管道输送系统的设计和运行参数的优化提供动;将煤粉(双峰分布)看作大小不同的2种固相同时考虑固相与液相、固相与固相之间的动量交在高浓度浆体流动的研究中, Eulerian方法为换模型的有效性通过文献[l3]的测量数据和描述液-固两相流动特性提供了合理的框架,不但cws压降试验数据验证;基于CWS实际流动工况可以全面地考虑液相颗粒相和壁面之间的相互作进行模拟,重点考察入口段内速度及浓度分布过程用,而且可以在动力学规律上描述各相内部的动力和入口段长度的变化规律学性质,固相还可以考虑碰撞与摩擦.基于 Eulerian方法,H"建立了适用于非胶体特性浆体的多1数学模型相流模型,其中考虑了固相对液相湍流的影响.As-1.1基本守恒方程sar2基于各相的动量守恒方程和浓度的对流扩散连续方程方程,建立了高浓度液固两相流动模型.作为Eule(1a)rian模型的简化形式,混合模型( mixture model)V·(ap,4,)=0i=1,2(1b)泛用于预测浆体流动压降、浓度和速度分布3动量方程其中引入代数滑移方程描述液固两相的相对运动,模拟计算只需求解混合相的连续方程、动量方程以及颗粒相的体积分数方程,可直接引入湍流模型进∑Bn(a,-H1)+apg行计算Hsu和Asar的模型取得了一定的成功,V·(ap,·")=-awVp-卩ps+V·T+但封闭控制方程需大量的经验式或需通过试验确Bn(u-M)+6:(4-)+aJPg(2b)定相关系数,在CWS这类高浓度浆体流动的数值式中,ar,a1和P分别为液相的体积分数速度和密模拟研究中受到较大限制;虽然混合模型计算量较度;an,"和p分别为第i固相的体积分数速度少,但只在中等浓度以下获得合理结果最近,在和密度.在任意空间点满足Eulerian多相流模型的框架下采用颗粒动理学理(3)论( kinetic theory of granular flow,KTGF)封闭控制方程成为高浓度液-固两相流动模拟的一个新的r=aμa[Vu1+(Vu1)(V·)l发展方向.KTGF理论能够合理描述颗粒间的摩擦4a与碰撞,赋予固相特性参数和输运系数明确的物理意义,避免过多的经验处理,使 Eulerian模型可以T =au,[Vu +(Vu)+求解复杂的高浓度浆体流动问题.应用该理论框2(V·4)l(4b)架文献[6-8]通过数值模拟研究了流化床内高式中,7和r为液固两相的应力张量;Vp和Vp浓度浆体的流动规律文献[9-11]应用 Eulerian为各相共享压力及颗粒碰撞产生的固相压力B方法和KTGF理论模拟直管和弯管内液固两相为液固相间的动量传递系数,计算中仅考虑曳力作用,当a≤0.8时B采用 Ergun公式计算,当ar水煤浆的煤粉颗粒通常呈现双峰分布粗、细>0.8时采用Wen和u公式1):为两固相间的颗粒的粒径差别较大,在CWS流动问题的研究中若将煤粉看作单一粒径的颗粒会产生较大偏差合动量传递系数理的方法是按颗粒分布特性将其看作2种固相,在(d+d)gJ模型中考虑不同固相之间的相互作用和影响E月,=中国煤化工lerian方法与KTGF理论相结合为这类问题的研CNMHG(5)究提供了良好的理论基础2).目前,有文献报道了式中,e为碰撞恢复系数;CB,为摩擦系数,取值404东南大学学报(自然科学版)第40卷0.15;d和d为颗粒直径1.2固相本构方程S.。sd8oa+d8o业KTGF理论认为流体-颗粒两相流动中的颗式中,a,m为固相混合物的最大填充分数粒相在微观尺度上的特性与分子尺度下的气体分1.3湍流模型子特性类似,颗粒的脉动速度遵从 Maxwellian分液固两相湍流流动采用 RNG k-e模型描述布通过类比定义颗粒的温度和建立颗粒拟温度液固两相具有相同的湍动能及湍动能耗散率为方程表示颗粒相的脉动和相应的传递特性V·(pn"nk)=V·( a, k)+颗粒拟温度方程13v·(p.a,nn)=(-p,I+x):V.+v·(pnne)=V·(a4mmVe)+R(16b·(knV6)-ys++Ds(6)式中,6,=1(2)为颗粒拟温度,反映颗粒间碰式中,n和n为各相的平均密度和平均速度μmcm和μn:为各相平均有效黏度和湍流黏度撞产生的小尺度脉动强度,(c2)为颗粒脉动速度的均方根;k为颗粒的相扩散系数;y。为颗粒间2模拟对象和计算方法的非弹性碰撞引起的脉动动能耗散率;dw为颗粒21模拟对象相的脉动动能在连续相中的耗散;D为颗粒的能模拟计算以水平管内CWS流动阻力试验为量耗散率,分别定义如下:基础试验装置、测量方法见文献[19].模拟工况和基本参数如下:煤粉真密度1.465;煤粉平均直=34(1+6)861+a8(4+en)径d=134.5μm;浆体温度T=20℃;试验管径D=32mm试验测量了3种体积分数和不同流速下2a p,d, gou(1 +eu)(7)直管稳定段上的流动阻力浆体的总体积分数q分别为30%,41.7%和49.5%.流速在1~5m/s煤.=3(1-6)8p, d- 4-w粉粒径分布见图1,按双峰分布特性,将煤粉看作(8直径由65,345μm的2种球形颗粒按体积比3:1混dn=-3B76(9)合而成3.0da1-un|2(10峰值1固相压力p包含动力学和颗粒碰撞2部2.0分m,即P,=Pa, 0 +223(1+eu)P, 0, 8o.@w ay固相剪切黏度μ和体积黏度A。1m为煤粉粒径/m图1煤粉粒度分布a(1+e)2.2计算方法801+3a计算采用管道长度为200D,基于对称性取管道一半作为计算区域;端面网格采用O型辐射网3u.00+y(1)格计算区域的离散化采用有限容积方法方程的离散采用一阶迎风格式,同时求解流动方程体积分数中国煤化工方程压力和速辐射分布函数g0n和801为度的CNMHGE方法在壁面,-[1-(a)1+05处液相来用无移週齐承仟,相采用 Johnsond等提出的部分滑移边界条件.在管道入口,假设(14)固相的体积分数一定,各相混合均匀,人口速度均第2期陈良勇,等:水平管入口段内水妹菜流动特性数值模叔匀,径向和切线速度为零;湍动能为来流动能的1.5%,湍动能耗散率为0.1k2.管道出口采用充分发展管流条件;颗粒间及颗粒与壁面间的碰撞恢复系数取0.95~0.993结果与讨论3.1模拟与试验结果的对比0.5Kaushal等实验测量了双颗粒组分浆体流动在稳定直管段上的压降和垂直径向上的浓度分J-1布浆体由水和2种直径的玻璃珠混合构成玻璃0.20.40.60,8珠密度为2470kg/m3;玻璃珠平均直径分别为125(a)I m/s(b)2 m/s和440μm,均为窄筛分分布;试验管道直径54.9mm,水平布置;管道入口浆体总体积分数为40%,粗细颗粒体积比为1:1,流速v=1~5m/s.图2给出了用本文所建模型得到的压降模拟值与试验值(L为管长;△P为管长L上的压差)图3给出了不同平均流速下垂直径向上(R为管道半径;r为离管道中心的径向位置)总浓度分布的模拟值与0.5测量值图中还包括粗细玻璃珠的分浓度φL的模拟值由图可见,在不同流速下,总浓度分布的模拟1.0值与测量值基本吻合.分浓度的计算结果显示,在(e)3 m/s(d)5 m/s低流速下,粗玻璃珠在重力作用下向管道底部沉4/m:25模拟);a40模拟);-04+25(模拟);降,细玻璃珠在垂直径向上分布相对均匀,这与文■440+125(测量)献[13]的试验结论一致图3垂直径向上玻璃珠的浓度分布35试验值模拟值5V/(m1)图2125和40μm玻璃珠混合浆体的压降图4cwS的压降模拟值与试验值图4给出了稳定直管段上Cws流动的压降管道断面的大部分区域;在管道底部,少量细颗粒模拟值与试验值,其中,离散点为试验值对应的连开始被增加的粗颗粒“挤出”,细颗粒体积分数由续曲线为模拟值由图可见,各浓度下的模拟值与31.3%减小为30%.随x的增加,粗颗粒体积分数试验值基本吻合,两者平均误差小于15%9%-12%的中间区域面积不断缩小,体积分数低3.2入口段内的浓度分布发展过程于9%和高于12%的上、下两区域相应地增大,与图5列出了φ=41.7%,V=0.5m/s时入口此同时粗颗粒平均体积分数在上部区域逐步减段内各相体积分数的分布发展过程在轴向距离x小在下部区域逐步增加如图5(a)中x=10D5D时,各相体积分数的显著变化仅限于靠近上、30D中国煤化工显示,随x的增下管壁的较小区域内此时粗颗粒在重力作用下加,CNMHG管道下部较小区开始向管道底部聚集,管道上部体积分数开始减域内;在中、上部区域细颗粒体积分数仅有少量增小底部体积分数开始增加,9%-12%的范围占据加且分布始终较为均匀对比粗细颗粒的分布过406东南大学学报(自然科学版)第40卷I=5D-00. 0401A(a)粗颗粒16、x=5Dx=60024-y312-0103T0.306(b)细颗粒-064-064---{02Y0.60]0图5入口段内CWS各相体积分数的变化过程(g=4.7%,v=0.5m/s)程可以看出,在管道底部粗颗粒沉积越多,被“挤1.0出”的细颗粒也越多.可见,除粗颗粒的强沉降特性外粗、细颗粒间的相互作用也是导致各相重新分布的主要动力对比图5(a)~(c)中x=30D和x=60D可见,x=30D时各相的分布已经完成,颗粒体积分数分布不再变化,水分由最初的均匀分布达到管道上部体积分数增多、下部体积分数减少的分布随轴向距离的变化过程在入口段的初娘1分布格局图6给出了入口段内垂直径向上粗颗粒浓度浓度分布(q中国煤化工浓度均匀分布的中间区域迅速缩小;在x=30D左右,浓度分布不再变化流体CNMHG沉降性的浆体3.3入口段内轴向速度分布发展过程入口段内速度分布的发展过程还与固相浓度的分对单相流体,人口段内的速度分布发展过程与布过程有关图7给出了q=4.7%,V=0.5m/s第2期陈良勇,等:水平管入口段内水煤浆流动特性数值模拟时入口段内轴向速度分布的变化过程由于液固两域的速度分布均匀而且垂直和水平径向上的速度相间的相对速度较小,此处速度值采用各相的平均分布较为相似随x的增加,壁面对速度分布的影速度在x=5D处,由于浆体黏性和壁面影响,壁响不断向中心扩展的同时,颗粒浓度分布不均的影面附近的流速迅速降低,为保持流动的连续性,中响逐步显现,2个方向上的速度分布差异显著增心区域的流速相应地增加,但此时壁面的影响仅限加,即在靠近管道中间区域,垂直径向上的速度梯于靠近壁面的有限区域内在x=5D~30D的范围度显著增大,水平径向上的速度梯度变化相对较内,随x的增加壁面附近的流速进一步降低,壁面小,这导致端面内的等速度线呈椭圆形分布阻滞的影响逐步向管道中心扩展;当x=30D时壁面影响已扩展到管中心区域对比速度分布和浓度分布的发展过程看出,x=30D时浓度分布过程已基本完成,速度分布在x从30D增加到60D时仍然有较大变化,并在x=60D时才达到稳定.可见,浓度分布过程完成时速度分布尚未完成,速度的稳定长度要比浓度的稳定长度更大,大约是浓度0.5稳定长度的2倍因此,人口段长度应以速度的稳定长度计算由图7可见,在x<15D的范围内等速度线基V/(m·s-1)()垂直径向本呈圆形分布,这说明颗粒分布不均的影响并不显著;当x>15D时,等速度线开始偏离圆形分布,由重力和颗粒间相互作用导致的总浓度分布不均以0.5及粗、细颗粒分浓度的分布不均开始对速度分布产生显著的影响,随x的增加,等速度线偏离圆形的三品程度显著增加图8给出了入口段内垂直和水平径向上轴向速度分布的变化过程.在人口段的初始段内,壁面对速度的影响仅限于靠近壁面的较小区域,中间06v/(m·8-1)=10Df0x4=15D(b)水平径向图8入口段内轴向速度分布(q=41.7%,v=0.5m/s)画、3.4入口段长度人口段长度L是入口段流动分析的重要参数,在经过L。长度后,管道截面上的流动参数不再受进口处速度分布的影响图9给出了L随浓度和速度的变化规律各浓度下,L。均先随V的增加而迅速增加,然后随V的增加而缓慢降低,这是因为在较髙的流速下,浆体湍动能和湍流黏度迅速增x=30D〓60D么、良中国煤化工CNMHG图7入口段内的速度分布过程图9入口段长度408东南大学学报(自然科学版)第40卷大导致固相浓度分布和速度分布更快地达到了稳dict the velocity and concentration profile for solid-liquid定状态另外,L随φ的增加而减小,这是因为浓lurry flow in pipelines C]//The 17th intemational度增高导致浆体的整体黏度增加,从而使浆体流动Conference on the Hydraulic Transport of Solids. Capetown Southen African 2007: 149-176在较短的距离内达到稳定状态[10]刘永兵陈纪忠阳永荣管道内液固浆液输送的数4结论值模拟[浙江大学学报,2006,40(5):858-863Liu Yongbing, Chen Jizhong, Yang Yongrong Numeri-1)当流速和固相浓度在较大范围内变化时,cal simulation of liquid-solid two-Phase flow in slurry所建模型都能很好地预测水煤浆管内流动,获得各pipeline transportation[ J] Journal of Zhejiang Univer相浓度分布和速度分布及压降等重要流动信息iy,2006,40(5):858-863.( in Chinese)2)粗颗粒的强沉降特性和粗、细颗粒间的强11hiho, David S, Koshi o. Numerical study on de-烈作用是导致各相重新分布的主要动力人口段relopment of particle concentration profiles in a curvedmicrochannel [J]. Chemical Engineering Sciend内,浆体由均匀混合状态逐步过渡到粗颗粒沿重力2006,61(11):3714-3724.方向上具有较大浓度梯度细颗粒整体分布相对均[12]HeYR, Chen H s, Ding YL,etal. Solids motion and匀的状态segregation of binary mixtures in a rotating drum mixer3)速度稳定长度约是浓度稳定长度的2倍[J. Chemical Engineering Research and Design颗粒浓度分布不均对速度分布产生显著影响主要007,85(7):963-973发生在入口段的后段[13]Kaushal D R, Kimihiko S, Takeshi T, et al. Effect ofparticle size distribution on pressure drop and concer4)入口段长度随浓度的增加而减小,随平均tration profile in pipeline flow of highly concentrated流速先增加后减小slurry[ J]. Intemational Journal of Multiphase Flow参考文献( References)05,31(7):809-823[14] Skudarmov P V, Lin CX. Ebadian M A Double-species[1 Hsu F L. flow of non-colloidal slurries in pipelines, ex-slurry flow in a horizontal pipeline [J] Journal of FIperiment and numerical simulations[ D]. Chicago, USAids engineering,2004,126(1):125-132University of Illinois, 1987[15] Gidaspow D. Multiphase flow and fluidization: continu-[2]Assar M H. Theoretical and experimental study of slurryum and kinetic theory descriptions [M]. New York:flow in pipe[ D]. Ohio, USA: Case Westem Reserve UNew York Press. 1994.[16]Ding J, Gidaspow D. A bubbling fluidization model u-[3] Ling J, Skudamov P V, Lin CX, et al. Numerical inves-sing kinetic theory of granular flow[J].AIChE Jourtigations of liquid-solid slury flows in a fully developedna,1990,36(4):523-538.urbulent flow region [J]. Intemarional Joumal of Heat [17] Lun C KK, Savage S B, Jeffrey D J, et al. Kinetic theand Fluid Flow,2003,24(3):389-398ories for granular flow inelastic particle in Couette flow[43 Lin C X, Ebadian M A. A numerical study of developingand slightly inelastic particles in a general flow fieldslurry flow in the entrance region of a horizontal pi[]. Journal of Fluid Mechanics, 1984, 140(2): 223-[J]. Computers&Fwds,2008,37(8):965-974.[5] Xu Rouelle A, Smith K M, et al. Two-phase flow of [18]Ogawa S, Umemura A, Oshima N On the equations ofsolid hydrogen particles and liquid helium[ J].Cryogen-ully fluidized granular materials [J]. Z Angew Mathics,2004,44(6/7/8):459-466.,31(4):4836] Jack t c, Fariborz t, Renaud e,ta, CFD modeling o[19]陈良勇段钰锤刘猛等温度和固相粒径与浓度对a liquid-solid fluidized bed[ J]. Chemical Engineering水煤浆管内流动壁面滑移的影响[].化工学报Science,2007,62(22):6334-63482008,59(9):2206-2213[7 Paola L, Renzo D F, Roberta P, et al. CFD modelling ofChen Liangyong, Duan Yufeng, Liu Meng, et al. Effectsliquid fluidized beds in slugging mode [J].Powderof temperature, solid particle size and concentration onTechnology,2006,167(2):94-103[8] Panneerselvam R. Savithri S, Surender G D CFD basedJournal of Chemical Industry and Engineering, 2008中国煤化工due to solid motion in liquid fluidised bed J]. Chemical [20]Engineering Joumal, 2007, 132(1/2/3): 159-1CNMHG IJ可[9] Lahiri S K, Ghanta K C Computational technique to pre-Fluid Mechanics, 1987, 176(2): 67-93