孙友勋,谭章荣,范瑾初? (同济大学 环境科学与工程学院,上海 200092)
摘 要:异波折板絮凝池的水头损失在设计中均采用明渠渐放和渐缩公式来计算,本研究在模型试验中发现利用给排水设计手册计算值与实测值相差较大。根据 模型试验数据,通过量纲分析,用最小二乘法拟合得出异波折板单元的半经验半理论水头损失计算公式。 关键词:异波折板絮凝池;水头损失;絮凝 中图分类号:TU991.22 文献标识码:B 文章编号 :1000-4602(1999)10-0029-04 异波折板絮凝池是在隔板絮凝工艺的基础上由我国研究发展起来的净水新工艺,具有结构紧凑、絮凝时间短、效果好、维护管理方便等优点,自80年代以来在新建及改造水 厂中得到了广泛应用,然而折板单元水头损失计算还存在许多问题。 1 量纲分析 异波折板单元如图1所示。 假定折板单元是光滑的,即忽略粗糙度的影响,则折板单元的压力损失(ΔP)是表1 中各变量的函数。 写成函数形式: ΔP=f(μ,L,B,R,ρ,v,α)? (1)? 式中 ΔP——因变量? ? μ、L、B、R、ρ、v、α——自变量? 根据π定理[1],Z=8,r=3 ?N(π)=Z-r=5 (2) 式中?Z——所有变量的总数 ?r——基本量纲的数目 ? N(π)——独立的无量纲数的最大数目? 表1 压力损失变量表表1 压力损失变量表变量名称 | 符号 | 量纲 | 流体粘度 | μ | ML-1T-1 | 流体密度 | ρ | ML-3 | 波峰流速 | υ | LT-1 | 折板板长 | L | L | 波峰间距 | B | L | 波峰水力半径 | R | L | 转折角度 | α | 无 | 压力损失 | △P | ML-1T-2 | 无量纲数π的一般表达式可以写为: ? π=(ΔP)aμbLcBdReρfvgαh?(3)? 式中?a、b、c、d、e、f、g、h?——各变量的量纲幂 它们的组合应使π成为无量纲数(由于α本身无量纲,可将α视作无量纲数)。? 将各变量的相应量纲代入并整理得: ? π=Ma+b+f L-a-b+c+d+e-3f+gT-2a-b-g? (4)?? 因变量是7个,方程只有3个,所以有4个变量可以任意选择。选择4组a、b、c、d的值代入方程组 ① a=1,b=0,c=0,d=0,得e=0,f=-1,g=-2,于是π1=ΔP/ρυ2 (7) ② a=0,b=1,c=0,d=0,得e=-1,f=-1,g=-1,于是π2=μ/Rρυ (8) ③ a=0,b=0,c=1,d=0,得e=-1,f=0,g=0,于是π3=L/R (9) ④ a=0,b=0,c=0,d=1,得e=-1,f=0,g=0,于是π4=B/R (10) 再加上五量纲数π5=ah,因此存在一函数关系f(π1,π2,π3,π4,π5)=0,即 ΔP/ρυ2=f(μ/Rρυ,L/R,B/R,α) (11) 这样,问题便简化为求解五个无量纲数之间的关系。 参照管道水头损失公式的形式,假定折板单元压力损失公式是以下形式: ΔP/ρυ2=K(μ/Rρυ)α1(L/R)α2(B/R)α3αα4 (12) 式中 K——无量纲常数 α1,α2,α3,α4——幂次,待定 无量纲数α也可能是sinα、cosα的形式。 两边取对数: ln(ΔP/ρυ2)=lnK+α1ln(μ/Rρυ)+α2ln(L/R)+α3ln(B/R)+α4lnα 令y=ln(ΔP/ρυ2),α0=lnK,x1=ln(μ/Rρυ),x2=ln(L/R),x3=ln(B/R),x4=lnα 则公式变为: y=α0+α1x1+α2x2+α3x3+α4x4 式(15)是多元变量线形方程。这样,可以根据水力试验数据所得的N组数据(yi,X1i,X2i,X3i,X4i)(i=1,……,N)利用最小二乘法拟合解出α0,α1,α2,α3,α4得到水头损失计算公式。 2 试验结果 2.1 模型试验设计? 由于进行生产性试验比较困难,故设计了一组异波折板单元模型,材料采用平板玻璃,符合 壁面光滑的假定。设计模型参数见表2。 表2 模型参数模型编号 | 板长L(mm) | 夹角α | 波峰间距B(mm) | M1 | 100 | 90° | 10 | M2 | 100 | 120° | 10 | M3 | 150 | 90° | 20 | M4 | 150 | 120° | 10 | 2.2 实测结果与给排水设计手册计算结果对比 模型单元水头损失用特制的测压板测量,实测水头损失与给水排水设计手册计算值的比较见表3。 表3 实测水头损失与手机计算值比较 波峰流速(m/s) | 实测水头损失 | 手册计算值 | 渐扩段(mm) | 渐缩段(mm) | 总水头损失(mm) | 渐扩段(mm) | 渐缩段(mm) | 总水头损失(mm) | 0.3 | 4 | 1 | 5 | 2.3 | 5.0 | 7.3 | 0.4 | 8 | 2 | 10 | 4.1 | 8.9 | 13.0 | 0.445 | 11 | 2 | 13 | 5.0 | 11.0 | 16.0 | 0.556 | 17 | 2 | 19 | 7.8 | 17.2 | 25.0 | 由表3可见,折板单元实测水头损失总是渐扩段大于渐缩段,这与文献[2]中的阐述相符。而手册计算则正相反。关于总水头损失,手册计算所得值略大于实测值,这可能是因壁面粗糙及水头损失系数取值太大所致。从模型试验数据来看,手册中渐扩与渐缩段水头损失计算公式值得商榷。实测渐扩与渐缩水头损失值与手册计算值相差较大可能是以下原因:? ①手册计算系将明渠渐扩及渐缩段简单相加。实际上折板单元内渐扩和渐缩段组成一个整体,其中流态较为复杂(见图2),与两者简单相加时流态不同。? ②波谷流速的计算取值影响。手册计算波谷流速采用流量除以面积来确定,而折板单元中由于存在较强的紊流涡旋,故以断面平均流速来计算流速水头便存在一定的误差。折板单元波谷断面流速分布如图3所示。? 以上分析可见,目前异波折板絮凝单元的水头损失计算公式在理论和实际中都存在重大的缺陷。? 根据量纲分析,我们将渐扩与渐缩段视为一个单元,由于两端流速相同,由流速计算引起的误差可以抵消,通过对模型水力试验数据的拟合,可以得出一个针对异波折板絮凝单元的水头损失计算公式。 2.3 最小二乘法拟合水头损失计算公式 四组模型各进行四种流速的试验,共得16组数据如表4。 表4 试验结果 模型编号 | 流量(m3/h) | 折板长度L(mm) | 折板夹角α | 粘度μ(Pa.s) | 波峰流速υ(m/s) | 波峰间距B(mm) | 水的密度ρ(kg/m3) | 波峰水力半径R(mm) | 总水头损失△P(mm) | M1 | 0.54 | 100 | 45° | 0.00113 | 0.3 | 10 | 1000 | 4.17 | 5 | 0.72 | 100 | 45° | 000113 | 0.4 | 10 | 1000 | 4.17 | 10 | 0.8 | 100 | 45° | 0.00113 | 0.445 | 10 | 1000 | 4.17 | 13 | 1 | 100 | 45° | 0.00113 | 0.556 | 10 | 1000 | 4.17 | 19 | M2 | 0.54 | 100 | 60° | 0.00113 | 0.3 | 10 | 1000 | 4.17 | 6.5 | 0.72 | 100 | 60° | 0.00113 | 0.4 | 10 | 1000 | 4.17 | 11.5 | 0.8 | 100 | 60° | 0.00113 | 0.445 | 10 | 1000 | 4.17 | 15 | 1 | 100 | 60° | 0.00113 | 0.556 | 10 | 1000 | 4.17 | 23 | M3 | 0.54 | 150 | 45° | 0.00113 | 0.15 | 20 | 1000 | 7.14 | 1.5 | 0.72 | 150 | 45° | 0.00113 | 0.2 | 20 | 1000 | 7.14 | 2.5 | 0.8 | 150 | 45° | 0.00113 | 0.222 | 20 | 1000 | 7.14 | 3.3 | 1 | 150 | 45° | 0.00113 | 0.278 | 20 | 1000 | 7.14 | 4.5 | M4 | 0.54 | 150 | 60° | 0.00113 | 0.3 | 10 | 1000 | 4.17 | 7.5 | 0.72 | 150 | 60° | 0.00113 | 0.4 | 10 | 1000 | 4.17 | 13 | 0.8 | 150 | 60° | 0.00113 | 0.445 | 10 | 1000 | 4.17 | 17 | 1 | 150 | 60° | 0.00113 | 0.556 | 10 | 1000 | 4.17 | 24 | 2.4 公式的推广? 上述水头损失公式是在假定壁面光滑的情况下得出的,而实际生产中,折板单元壁面粗糙度 不应忽略。故而可考虑加上一个修正系数K,即公式变为: 依据公式 (14)计算出16组(X1i,X2i,X3i,X4i,yi)。利用计算程序计算各系数α0,α1,α2,α3,α4得(当角度以 αh形式出现时): α0=-1.50264,α1=0.0073,α2=0.27,α3=0.46,α4=0.64,代入式(12)得: ΔP/ρυ2=0.223(μ/Rρυ)0.0073(L/R)0.27(B/R)0.46α0.64 (16) 以ΔP=ρgΔh、ρ=1000代入,Δh为水头损失,得: Δh=0.424[μ0.0073L0.27B0.46α/R0.733]·[v1.99/2g] (17) 由式(17)可以看出,异波折板单元水头损失受L、μ、B、R、α及υ的共同影响,且水头损失为流速的1.99次方。对于一定的折板单元,受μ的影响较小。以水温20℃时粘度μ=1.01×10-3代入式(17)得: Δh=0.400[L0.27B0.46α0.64/R0.733]·[v1.99/2g] (18) 同理,当角度以(sinα)h形式出现时,计算可得:a0=1.34253,a1=0.0073,a2=0.27,a3=0.46,a4=-0.53,代入式(12)得: Δh=0.472[L0.27B0.46(sinα)0.91/R0.733]·[v1.99/2g] (19) 当角度以(cosα)h形式出现时,计算可得:a0=-1.84128,a1=0.0076,a2=0.46,a4=-0.53,代入式(12)得: Δh=0.286[L0.27B0.46/(R0.733(cosα)0.53)]·[v1.99/2g] (20) 式(18)-(20)即为异波折板单元的水头损失计算公式。 2.4 公式的推广 上述水头损失公式是在假定壁面光滑的情况下得出的,而实际生产中,折板单元壁面粗糙度不应忽略。故而可考虑加上一个修正系数K,即公式变为: 式中 K——由试验确定的系数? 由于折板单元高度与角度是sin α关系,而波谷宽度与角度是cos α关系,因此式 (20)、(22)都是合理的,建议实际计算中采用。因 它们是利用试验数据解出的理论公式,故而可称之为半经验半理论公式。在絮凝试验中用此 公式计算的水头损失与实测值符合得较好。? 3 结论? ①通过模型水力试验,得出了异波折板絮凝单元的水头损失计算公式。公式虽然在模型试验中得到了较好的验证,但由于壁面粗糙度等因素的影响,应用于生产实践中尚待进一步研究。? ②试验中发现异波折板絮凝单元的水头损失渐扩段大于渐缩段,与手册计算公式相反。因此,研究异波折板絮凝单元的水头损失计算公式,具有理论及实践的双重意义。?? 参考文献:? [1]李士豪.流体力学[M].高等教育出版社,1990. [2]康士坦丁诺夫 M.水力学[M].钟用升译.江西高校出版社,1990.
电 话:(0)13901720712 (021)65017139? 收稿日期:1999-02-11 |