江 山,張海金,齊 鵬
(綏化水文局,黑龍江綏化152000)
在管道水力系統(tǒng)中,沿程水頭損失計(jì)算是水力計(jì)算的一個(gè)重要環(huán)節(jié)。在實(shí)際工程中,為方便計(jì)算,不管是層流和紊流,管線的沿程水頭損失都是采用層流時(shí)所求的沿程阻力系數(shù)來(lái)進(jìn)行計(jì)算,而大多數(shù)情況下管道中水流流態(tài)都是紊流的,紊流狀態(tài)下的沿程阻力系數(shù)與層流狀態(tài)下的沿程阻力系數(shù)是有所差別的,其公式也不一樣。因此,本文分別分析了各種水流狀態(tài)下沿程阻力系數(shù)求解公式,并利用Excel軟件介紹了一種求解柯列勃洛克公式的簡(jiǎn)便方法。
層流狀態(tài)下,沿程阻力系數(shù)只與雷諾數(shù)有關(guān),與管壁粗糙程度無(wú)關(guān),其計(jì)算公式一般由謝才公式得來(lái):
式中:g為重力加速度;C為謝才系數(shù)。
謝才系數(shù)公式可由曼寧公式算出:
式中:R為水力半徑;n為管道的糙率。
紊流中沿程阻力系數(shù)與雷諾數(shù)及粗糙度之間的理論,在理論上至今沒有完全解決,1933年,尼古拉茲根據(jù)大量人工粗糙管試驗(yàn)發(fā)表了反映圓管流動(dòng)情況的實(shí)驗(yàn)結(jié)果。描繪出了尼古拉茲實(shí)驗(yàn)曲線圖。人們?cè)诤髞?lái)的研究中依據(jù)普朗特混合長(zhǎng)度理論和尼古拉茲試驗(yàn)數(shù)據(jù)提出了紊流沿程阻力系數(shù)的半經(jīng)驗(yàn)公式。
1.2.1 紊流光滑管區(qū)
在Re=5×10 ~3×10情況下,沿程阻力系數(shù)滿足尼古拉茲光滑管公式:
1.2.2 紊流粗糙管區(qū)
此區(qū)黏性底層失去意義,絕對(duì)粗糙度對(duì)水頭損失起決定作用。在此區(qū)沿程阻力系數(shù)滿足尼古拉茲粗糙管公式:
式中:r0為管道半徑;△為管道絕對(duì)粗糙度。
而在工業(yè)管道水力計(jì)算應(yīng)用較為廣泛的是柯列勃洛克公式,該公式是柯列勃洛克根據(jù)大量工業(yè)管道試驗(yàn)資料提出來(lái)的,實(shí)際上是尼古拉茲光滑區(qū)和粗糙區(qū)公式的結(jié)合。其公式為:
式中:λ為沿程阻力損失系數(shù);d為管道直徑;△為管道絕對(duì)粗糙度;Re為雷諾數(shù)。其工業(yè)管道絕對(duì)粗糙度,見表1。
表1 常見工業(yè)管道絕對(duì)粗糙度
因該公式應(yīng)用比較麻煩,需要編制程序來(lái)進(jìn)行計(jì)算,在實(shí)際工程中應(yīng)用很少。本文現(xiàn)介紹一種在Excel中求解柯列勃洛克公式的簡(jiǎn)便方法。
現(xiàn)將式(1)中式子右邊部分移到左邊則變?yōu)?
根據(jù)給定管道參數(shù)、流量的情況下,雷諾數(shù)Re可通過(guò)公式算出,公式(2)中只有λ為未知數(shù)。在某一單元格1中輸入λ值,而在另外一單元格2中輸入公式2左邊部分,然后逐步調(diào)整λ值,直至單元格2中數(shù)值在滿足精度要求的情況下近似等于0時(shí),則此時(shí)的λ即為通過(guò)柯列勃洛克公式求解的沿程阻力系數(shù)值。
某一水管長(zhǎng)500 m,直徑為200 mm,管壁粗糙高度△為0.1 mm,輸送流量為10 L/s水的運(yùn)動(dòng)黏性系數(shù)為0.01310 cm2/s。經(jīng)計(jì)算雷諾數(shù)Re=48595。
利用上述求解柯列勃洛克公式方法求得λ=0.0225,查莫迪圖查得λ =0.0208,由此可以看出,該計(jì)算方法應(yīng)該說(shuō)是更為準(zhǔn)確的。
對(duì)于管道紊流狀態(tài)情況,利用Excel軟件來(lái)求解柯列勃洛克公式中的沿程阻力系數(shù),避免了編制繁瑣的迭代程序,簡(jiǎn)便易行,在工程實(shí)際中易于實(shí)現(xiàn),且計(jì)算結(jié)果與實(shí)際情況是比較吻合的。該方法同樣適用于其他流體的管道計(jì)算中。
[1]李大美.水力學(xué)[M].武漢:武漢大學(xué)出版社,2003.
[2]劉福祥.關(guān)于工業(yè)管道紊流阻力系數(shù)的柯列勃洛克公式迭代算法的收斂區(qū)間研究[J].大連大學(xué)學(xué)報(bào),2000(4):35-39.