直管内水煤浆流动阻力的数值实验

第39卷第3期2008年5月BOILER TECHNOLOGYMay.,2008文章编号:CN31-1508(2008)03-0001-04直管内水煤浆流动阻力的数值实验赵国华,陈良勇,刘猛,段钰锋(东南大学洁净煤发电及燃烧技术教育部重点实验室,江苏南京210096)关键词:水煤浆;非牛顿流体模型;阻力系数;数值实验摘要:运用非牛顿流体模型对水煤浆流经直管进行数值模拟,通过压差实验值与计算值的比较验证计算模型的正确性。基于正确的计算模型,进行数值实验研究沿程阻力系数。总结得出水煤浆流经直管时沿程阻力系数的经验公式。通过对速度场的分析,得出水煤浆流动的核心区。中图分类号:TQ53文献标识码:A沿程阻力系数定义1前(1)水煤浆管道输送比其它运输方式如卡车或火车具有明显的优势,同样具有经济性。管道输送是可靠的不受恶劣气候如暴风雪或低温的影式中:△p—测量段压差,Pa;密度,kg/m2;响,并且由于一定程度的自动操作,人工劳动相V—流速,m/s;对减少,也开辟了输送能源新技术的开端,那就d一直径,m;是远距离高容量煤浆管道输送。水煤浆在直管中输送,存在沿程阻力损失,对水煤浆流经直管L—测量长度,m。为了研究沿程阻力系数对非牛顿流体在圆时沿程阻力的数值实验研究具有重要的意义管中的流体力学特性,许多研究者2针对不同本文采用管流法和旋转粘度计法相结合,得出水的研究对象和研究目的,提出不同的雷诺数定煤浆的流变方程,运用非牛顿流体模型对浓度分义,本文采用基于壁面表观粘度的雷诺数别为57.08%和59.10%的神华水煤浆进行数3n+1、8值模拟。通过实验值与计算值的比较,验证计算Re,=_evd模型的正确性;基于正确的计算模型对管径分别式中P—密度kg/m;为25mm32mm、40mm和50mm的直管进行V—流速,m/s;数值实验,研究沿程阻力系数。d—直径,m;2沿程阻力系数v—表现粘度流变指数。1933年尼古拉兹采对一系列管道进行实验,得出了管道沿程阻力系数变化的一般规律:3水煤浆的流变方程(1)雷诺数Re<2000时,水流为层流,λ与运用管流法和旋转粘度计法得出2个浓度Re呈倒数关系,且A=64/Re下水煤浆的流变方程:(2)20004000时,水流处于湍流状式中国煤化工态,受到的影响因素较多。CNMHG收稿日期:2007-07-03基金项目:国家重点基础研究发展计划(973计划)资助项日(2004CB2l7701)作者简介:赵国华(1983-),男,东南大学能源与环境学院硕士研究生,主要从事洁净煤技术研究锅炉技术第39卷计算模型进行数值模拟。通过计算得到实验压差与计算压差对比如下图:4.1控制方程水煤浆在管内流动是一个等温过程,不考虑管径25mm质量浓度温度变化,直角坐标系中定常条件、不可压缩流实验值体流动通用形式方程为:t ot pg)+div(pup)=div(rgradp)+S(3)中:卓—自变量;流体的密度U—速度矢量;r—扩散系数;流速/m-1s—源项量厂和S与ψ值有关图1管径为25mm中计算值与实验值对比当∮=1时通用微分方程变为质量守恒方程;当=L、乱、时分别为坐标系3个方向上的动管径32mm质囊浓度实验值量方程。M、42、为坐标系3个方向上的速度。4.2非牛顿流体运动方程描述非牛顿流体的流变方程,建立应力与剪切速率的关系,表达式为:r2=2+(g÷)式中:y—剪切速率,1/s;r——剪切应力,Pa;3.75流速/m1—屈服剪切应力,Pa图2管径为32mm中计算值与实验值对比m、g、C—常系数(5)管径40m实验值0计算值au,⊥aur(8)rx=rn=Bu4⊥a(9)au, au淀速加y24图3管径为40mm中计算值与实验值对++A(11)A-2y+2y+22y++2y+实验度5浓度管径50mm质量375}0-计算)2+rx、r、r为坐标系3个方向上的法向应力rx、rs、r、ra、r为坐标系3个方向上的切中国煤化工向应力CNMHG5数值模拟结果与实验结果对比对2个浓度的水煤浆运用非牛顿流体模块图4管径为50mm中计算值与实验值对比赵国华等直管内水蝶浆流动阻力的数值实验10管径25mm质量浓度在20%以内,从而运用非牛顿流体模型中的一实验值Herschel- Bulkley模型是正确的,运用此模型进行数值实验,探讨沿程阻力损失系数与流速、管径的关系是可靠、有工程应用的实际价值的6数值实验结果6.1沿程阻力系数对管径为25mm、32mm、40mm和50mm在各个流速下进行数值实验,探讨沿程阻力系数图5管径为25mm中计算值与实验值对比与雷诺数的关系。图9和图10是2个浓度下各个管径中的沿管径32mm质量浓度实验值59.10%程阻力系数与雷诺数的关系,并且与牛顿流体中沿程阻力系数与雷诺数的关系相比较。雷诺数采用基于表观粘度的雷诺数,从图中可以看出,沿程阻力系数随雷诺数的增加而减小,并趋于相同;随煤浆的浓度增加而增加;相同雷诺数下,水煤浆的沿程阻力系数比牛顿流体的沿程阻力系数大,随雷诺数的增加,水煤浆沿程阻力系数与流速/m-1牛顿流体的沿程阻力系数的差值越小。图6管径为32mm中计算值与实验值对比质量浓度质量浓度额12牛顿流体0400800120016002000240000051015202.530雷诺数流速m-4图9沿程阻力系数与雷诺数的关系图7管径为40mm中计算值与实验值对比都1240080011600200020030.6091.21.518雷诺数流速/m图10沿程阻力系数与雷诺数的关系图8管径为50mm中计算值与实验值对比中国煤化工力系数与雷诺从图1至图8中显示单位长度压差随流速的数的关HCNMHG增加而增加,基本呈线性变化。通过实验值与计浓度57.08%时:A=100/Re142算值比较,计算值与实验值相当接近,相对误差浓度59.10%时:A=100/Re15锅炉技术第39卷经验公式的形式与牛顿流体的沿程阻力系径中的无量纲速度分布应该是完全重合的,从图数公式A=64/Re相比还是相当接近的。l1图122个图可以看出无量纲速度分布重合度6.2速度分布非常好。速度从管壁到管中心逐渐增加,速度梯速度分布是反映水煤浆在管道中的流动状度越来越小,中间产生一个流动核心区,从图中况。图11和图12是2个浓度下水煤浆在各个管可以大概看出核心区为之间径中的无量纲速度分布,横坐标是无量纲半径,为各点半径与管道半径的比值纵坐标是无量纲7结论速度,为不同半径处的速度与管道平均速度的比运用非牛顿流体模型结合进口条件运用外值。理论上对于同一个浓度的水煤浆在不同管接函数处理,对2个浓度的水煤浆流经4个管径的直管进行数值计算,通过实验值和计算值的压差比较,验证计算模型的正确性。基于正确的计算模型对水煤浆流经直管进行数值实验,沿程阻力系数随雷诺数的增加而减小,并趋于相同;随煤浆的浓度增加而增加;相同雷诺数下,水煤浆的沿程阻力系数比牛顿流体的沿程阻力系数大,32随雷诺数的增加,水煤浆沿程阻力系数与牛顿流体的沿程阻力系数的差值越小。总结得出水煤浆流经直管时沿程阻力系数的经验公式,并与牛图11浓度57.08%无量纲速度分布顿流体流经直管时沿程阻力系数相比较,形式比较接近。水煤浆在直管流动过程中,产生的流动核心区大概为一3R-3R之间参考文献[]孔现.工程流体力学[M门.北京;中国电力出版社,199225[2] Cohen Y, Metner A. B. Apparent slip flow of polymer solutions[J].J. Sheol,1985,34(29):67-102[3]岑可法,姚强,曹欣玉,等.煤浆燃烧、流动、传热和气化的理论与应用技术[M]杭州,浙江大学出版社,1997[4]沈崇棠刘鹤年非牛顿流体力学及其应用[M]北京:高等教图12浓度59.10%无量纲速度分布育出版社,1989Numerical Experiment of Flow Resistance ofCoal-water Slurry in Straight TubeZHAO Guo-hua, CHEN Liang-yong, LIU Meng, DUAN Yu-feng(Key laboratory of Clean Coal Power Generation and Combustion Technologyof the Ministry of Education, Southeast University, Nanjing 210096, China)Key words: coal-water slurry; non-newtonian model; friction coefficient: numerical simulationAbstract:Numerical simulation of coal-water slurry flow in straight tube by Non-Newtoni-an model, the correctness of computational model was中国煤化工n of pres-sure drop between experimental values and calculatedmerical experiments were carried, the distance frictionCNMH GThe expe-rienced formula of distance friction coefficient was gained. By the study on velocity fieldcore region of flow in straight tube was given out