弯管内水煤浆流动阻力的数值实验研究

第35卷算10期鼻束電力Vol 35 No 102007年10月Fast china Electric Poweroct.2007弯管内水煤浆流动阻力的数值实验研究赵国华,陈良勇,段钰锋(东南大学洁净煤发电及燃烧技术教育部重点实验室江苏南京21006)摘要;结合管流法和旋转粘度计法对神华煤水煤浆进行了流变特性测量得出符合屈服一幂律流体模型的水煤浆流变方程。利用CFD软件平台建立了水煤浆∞η°弯管的流动模型,通过计算值与实验值的比较验证了计算模型的正确性。运用正确的计算模型对水煤浆流经弯管的局部阻力系数进行了多变量数值模拟实验,得出了弯管工程设计中最佳弯径比与流速的关系,总结出适合工程应用的水煤浆流经较小弯径比时的局部阻力系数的经验公式。且弯管内最大速度点随弯径比的增大往弯管外侧移动;随弯径比或弯曲角度的增大弯管内的速度场趋于一致。关键词:水煤浆;弯管流动;数值试验;阻力损失系数基金项目:国家重点基础研究发展计划(973计划)资助项目(2004cB217701)作者简介:赵国华(1983-),男,硕土研究生研究方向为洁净煤发电及燃烧技术。中图分类号:TQ53文献标识码;A文章编号:1001-9529(2007)10004Numerical experimentation on flow resistance of coal-water slurry in curved pipesZHAO Guo-hua, CHEN Liang-feng, DUAN Yu-fengKey Laboratory of Clean Coal Power Generation and Combustion Technology of Ministry of EducationSoutheast Univ., Nanjing 210096, ChinaAbstract: Based on the pipe flow method and the rotational viscometer method, the rheological behavior of Shenhuacoal-water slurry was measured and the rheological equation in conformity with the model of Herschel-Bulkley waser slurry now in a curved 90oby using CFD software platformand the computational model was proved correct through comparison between caleulated values and experimental vaues. With the correct computational model, numerical experiments on local resistance coefficient waB conducted,pro-ding the interrelations between the best bend diameter ratio and flow velocity in engineering design as well as theempirical equation of local resistance coefficient for practical engineering. It is also concluded that with increasing ofthe bend diameter ratio the dot of maximum velocity is near to the wall outside of a pipe bend, and the distribution ofthe velocity feld in the pipe tends to become the same with increasing of the bend diameter ratio or the bend angleKey words: coal-water slurry; pipe bend flow; numerical experimentation; frictional-loss coefficient水煤浆管道输送是发展大规模高效气流床煤出非牛顿流体流经弯管的局部阻力系数公式。对气化技术的基础性研究并应用于新型高效的洁具有非牛顿流体性质的水煤浆运用数值试验的方净煤发电技术ICC中。弯曲管道是水煤浆管道法就流经小弯径比的弯管局部阻力系数进行研究输送过程中常见的部件。当流体流过弯管时,由尚未见报道。于向心力及壁面压力梯度的作用,使得流体的运水煤浆是一种非牛顿流体,流经弯管时的局动状态变得极为复杂2),不仅在同一横截面上部阻力系数与牛顿流体不同,研究中发现影响水存在二次流动而且在弯头的进出口还存在流体煤浆流经弯管的阻力系数的因素很多包含水煤脱离壁面的现象。 Hwang and Pal1、 Deshpande浆的流变模型弯管的几何形状和浆体的流速等。and barigou4对非牛顿流体的突扩突缩管件进行结合数值模拟实验的方法对水煤浆流经弯管进行了实验研究, Jure Marr(运用二次方程流变模型数值试验对局部阻力系数进行研究可以减少实对剪切变稠的高浓度灰水混合物流经9°弯管进验周期、工作量和实验耗资。行过数值模拟, Sing and Mishra6在大弯径比下提本文从实验出发,总结出关于弯径比和雷诺赵国华等弯管内水煤菜流动阻力的数值实验研究7(总904)数的经验公式为工业设计和设备选型提供依据。-P+2B(4)1试验研究(5)1.1水煤浆流变方程的确定-P+2B水煤浆的流变模型很多,影响因素也很复杂对水煤浆在管道内的恒定剪切流动通过管流法Tr =Tr =B(-2测量其压降和流量,可推算剪切应力与应变速度的关系,从而确定流变方程。管流法一般可测量BC中、高剪切速率下的水煤浆流变方程;旋转粘度计B C法测量低剪切速率下的水煤浆流变方程。在25℃恒温下,对质量浓度为57.08%的神华水煤浆运用管流法和旋转粘度计法相结合的方法,测B=[出水煤浆在较大范围内的剪切应力与剪切速率的关系,运用屈服一幂律模型得出水煤浆的真实流A=[2()2+2()2+2(弘+)2+变方程为:r=845+(0.068·y)(10)1.2实验弯管式中T坐标系3个方向上的法向应力;Tn试验弯管管径为50mm,R/R分别取3、8、Tnnr—坐标系3个方向上的切向应力10.0为弯管部分对应的圆心角,R为弯管半3网格划分和数值计算条件径与管道半径的比值(弯径比),是衡量局部管件阻力损失系数的一个重要参数。试验分别测试了计算中网格采用四面体和六面体非结构网不同工况下水煤浆流经弯管时上、下游的压差格,节点距离在2~4mm之间。人口边界条件:水煤浆的密度为11982计算模型kg/m3,采用平均速度进口条件。由于水煤浆在2.1控制方程弯管中平均流速比较小,表观粘度较大,Re<水煤浆在管内的流动是一个等温过程,不考2200,流动可认为是层流流动。差分格式采用一虑温度变化,直角坐标系中定常条件、不可压缩流阶迎风格式,算法采用 SIMPLE算法收敛误差控体演流动通用形式方程为制在0.1%。(pd)+dv(pU)= div(grade)+S(1)4数值计算结果与试验比较式中φ为自变量,p为流体的密度,U为速度矢量,F为水煤浆是不透明、粘稠的非牛顿流体,在试验扩散系数,S为源项量厂与S、中值有关。当中=1时,通测量中只能对宏观参数进行测量(例如弯管上下用微分方程变为质量守恒方程;当φ=L2、B,时分别为坐标系3个方向上的动量方程。“为坐标系3个游的压差)。受测量水平的限制对速度场和表方向上的速度。观粘度场难以测量。在相同试验段计算出试验管2.2非牛顿流体运动方程件弯管上下游的压差,并与实测值进行对比。描述非牛顿流体的流变方程,应力与剪切速图1、图2和图3是不同弯径比下计算的弯率和流体速度的关系表达式为管上、下游压差与实测值的对比。可以看出计算式中y剪切速率;剪切应加命(2)值与实验值的结果比较接近,相对误差小于r=古+(·y)30%,说明运用非牛顿流体模型建立的数学模型服剪切应力;mψ、c—常系数。进行计算是可行的。基于此模型进行数值试验探讨局部损失系数与流速、弯管几何形状等的关=-p+2B(3)系为工程应用提供可以借鉴参考的设计参数。8(总905)束電力2007,35(10)R/D=1.5一试验值f=2-P18一-计算值式中p一人口端的压力,Pa;P2出口端的压力Paip水煤浆的密度,kg/m3;流速。5.2雷诺数在非牛顿流体流动中,雷诺数的定义很多,根据水煤浆的性质运用基于壁面表观粘度的雷诺数速度/mg比较合适,定义为图1弯径比为15时试验测量压差与计算值比较式中p—水煤浆的密度,kg/m3;a——流速,m/s;R/D=1.5试验值d—管径,m;—表观粘度,Pa·ss→计算值5.3数值实验工况影响水煤浆流经弯管的局部阻力系数的因素很多,例如弯管的几何形状、水煤浆的流变特性和浆体的流速等。本文主要考虑弯管的几何形状和浆体的流速对弯管的阻力系数的影响。对管径为50mm的弯管,弯径比按照工程上常见的标准取速度/mst弯径比为1.2、1.6、23、4、6、8、10、12时,进囗流图2弯径比为4时试验测量压差与计算值比较速分别为0.25、0.5、0.75、1、1.25、1.5、1.752m/s进行模拟计算,研究局部阻力系数与弯径比、试验值流速的关系0-计算值5.4数值实验结果与分析54.1局部阻力系数根据 Mishra- Gupta对大弯径比下总结出的非牛顿流体流经弯管与牛顿流体流经弯管局部阻力系数的关系,运用此方法对小弯径比下水煤浆流经90弯管的局部阻力系数进行研究,其关系如速度sg图4。由图4可见,雷诺数越大,水煤浆的局部阻力系数越偏离牛顿流体的局部阻力系数,且随着图3弯径比为6时试验测量压差与计算值比较弯径比的增大,变化明显。5弯管内水煤浆流动的数值实验研究5.1局部阻力系数对非牛顿流体流经弯径比不同的弯管局部阻力系数研究, Mishra-Gupta对弯径比大于12的管道进行了研究,提出了经验公式:f=1+0.033( logUe)(11)式中f——流经弯管的局部阻力系数;f=16/Re;De=工程上,局部损失系数与流速的关系通常通图4各弯径比下的局部阻力系数过局部损失系数与雷诺数的关系表现出来。局部从以上实验数据总结出水煤浆流经小弯径比损失系数定义为:时,局部损失系数与弯径比及雷诺数的关系如下。赵国华,等弯营内水煤浆流动阻力的数值实脸研究9(总906)18=-1.18+6.84P+(0.53-0.07P)(lgDe)(B41P6式中f—流经弯管的局部阻力系数=16/Re;De=Re√D2P;P=2R!/R。5.42最佳弯径比最佳弯径比是相同流速下水煤浆流经弯管产生的局部阻力系数最小时的弯径比。对弯径比从12~12各个流速下进行数值实验,得出各个流速下产生最小局部阻力系数的弯径比。最佳弯径图6=45°截面速度分布比出现在弯径比1.2~2.4之间,在弯径比1.8附近比较集中,实际工程90°弯管设计时弯径比取18左右较为合适。最佳弯径比随水煤浆流速的增加而增加,基本成线性变化。5.4.3速度分布沿着管轴线方向,在弯曲不同部位截取面图5、图6和图7分别是=30°45°60°面无量纲速度随管半径的变化分布图。由于离心力的作用,最大的速度不像流体流经直管一样出现在管径的中心处。弯径比为1.2和2时,弯曲管道中,内侧图7=60°截面速度分布的速度高外侧速度低,最大速度出现在靠近弯管生的局部阻力系数进行研究,总结出适合工程应内侧。随着水煤浆流经管道截面(6分别为30°、用的经验公式。(3)对弯管截面的速度场分析,45°、60°时),最大速度的偏离更加明显速度场由速度场随截面夹角θ的增大趋于集中最大速度点分散逐渐趋于集中化,最后几乎重合。随着弯径随弯径比的增加往弯管外侧移动最终趋于相同。比的增大,最大速度往弯曲管道外侧偏移在靠近参考文献:弯曲管道外侧处产生很大的速度梯度,弯曲管径外侧的剪切速率比预期的大,而内侧比预期的小。[1] Chisholm D. Brief communications of two-phase iow in bends[J]. Multiphase Flow,1980,6(2):363367R/R[2] Churchil S W.all no[冂]Chem. Engng Japan, 1977(84): 91-92.-6[3] HwangC Y J, Pal R. Flow of two-phase oil/wster mixturesEg,Japn,1997(68):57-63[4]ls N S, Bargou M. Foam flow phenomena in suddenexpansions and contractions[ J]. Multiphase Flow, 2001[5] Jure Marn, Primoz Temik Laminar flow of shear-hickeningfluid im a 90 pipe bend[ n]. Fluid Dynamics Researh, 2006图58=30°截面速度分布[6] Singh R P, Mishra P, Friction factor for Newtonian and non6结论Newtonian fuid flow in curved pipes[ J].Chem Eng.,Japan,80(13):275-280.(1)运用计算模型比较好地模拟出了水煤浆[7]岑可法,姚强曹欣玉,等煤浆燃烧流动、传热和气化流经弯管的流动状况。(2)通过数值实验得出最的理论与应用技术[M]杭州浙江大学出版社,1997佳弯径比与速度的关系,实际工程90°弯管设计收稿日期:2000747时弯径比取1.8比较合适;对水煤浆流经弯管产本文编辑:王延婷