顯式Colebrook-White 摩阻系數方程

苑偉民

中國石化天然氣分公司廣西液化天然氣工程項目部，廣西 北海 536000

1 Colebrook-White 方程

在計算摩阻系數的眾多方程中，Colebrook-White方程較為常用［1-3］：

式中： f 為為水力摩阻系數；ε 為為管道內壁粗糙度，m；D 為為管道內徑，m；Re 為為雷諾數，無量綱量。

2 Colebrook-White 方程求解

隱式Colebrook-White 方程［4-6］，早期借助圖表進行估算，誤差較大；隨著數值算法及編程軟件的發展，可運用牛頓法、弦截法、拋物線法等迭代方法［7-9］或者使用數學軟件進行求解［10-13］；但迭代法計算用時較長，應用數學軟件有可能不能直接得出數值解。

2.1 Lambert W 格式的方程

使用MATLAB 求解式（1）得到如下最簡形式：

Lambert W 函數定義為：W=Lambert W（x），表達的數學意義為：w*exp（w）=x。 LambertW 函數是一個超越方程的形式，不能直接應用于工程計算。

2.2 WrightOmega 格式的方程

使用Maple 軟件求解式（1）得到如下最簡形式：

2.3 Product Log 格式的方程

使用Mathematica 軟件求解式 （1） 得到帶有Product Log 函數的方程， 形式更為復雜。 Product Log函數與Wright Omega 函數定義相同，仍然不能直接應用于工程計算。

3 Colebrook-White 顯式方程

3.1 顯式方程的推導

在數值方法研究的基礎上，綜合牛頓法和多重迭代法［7］，提出以下迭代法：

使用式（5）代替求導：

式（4）具有4 階以上收斂階。

3.2 顯示方程的提出

將式（1）代入式（4），經過數學推導，簡化得到顯式Colebrook-White 方程：

4 計算實例

4.1 方程使用

4.1.1 已知參數

雷諾數Re=5 813 924， 粗糙度k=0.02 mm， 管徑D=600 mm。

4.1.2 計算過程

將已知數值帶入b、c、x0計算， 得到表1 數值，進一步帶入式（6），即可解得摩阻系數。

4.2 結果分析

將式（6）與文獻［10］中的方程（見表2）作對比見表3［10］。

表1 數值表

表2 顯式Colebrook-White方程［10］

表3 摩阻系數計算參數對比［10］

表3 中， 采用的精確值為1.050 336 403 193 8×10-2，該值是采用拋物線法和弦截法分別求解式（1）所得計算值的平均值。

生成100 組實驗數據進行計算，式（6）計算結果的平均相對誤差為（1.5×10-7）%，絕對誤差為1.5×10-11；使用計算機編程求解， 計算所用時長約為拋物線法的1/32，約為弦截法的1/9。

5 結論

顯式Colebrook-White 方程具有以下特點：

a）適用范圍廣，適用于所有雷諾數范圍；

b）計算結果精度高，絕對誤差數量級10-11；

c）計算時間短，遠小于迭代算法；

d）形式簡單，既適合手算也適合計算機計算。

［1］ Brown O G. The History of The Darcy-Weisbach Equation for Pipe Flow Resistance，In Environmental and Water Resources History［J］.ASCE Conference Proceedings 126 2002，126（4）：34-43.

［2］ Coelho M，Carlos Pinho. Considerations about Equations for Steady State Flow in Natural Gas Pipelines［J］.Braz.Soc.of Mech.Sci.&Eng.Vol.2007：XXIX，No.（3），262-273.

［3］［S.N.］［EB/OL］.http：//www.colebrookequation.com/，2008-05-06/2011-07-17.

［4］ Thomas G，Lester P E. Solving the Colebrook Equation for Friction Factors Part 1 of 3 in a Series of Articles［EB/OL］.http：//www.cheresources.com/colebrook1.shtml，2011-12-30/2012-01-05.

［5］ Thomas G，Lester P E. Solving the Colebrook Equation for Friction Factors Part 2 of 3 in a Series of Articles［EB/OL］.http：//www.cheresources.com/colebrook2.shtml，2011-12-30/2012-01-05.

［6］ Thomas G，Lester P E. Solving the Colebrook Equation for Friction Factors Part 3 of 3 in a Series of Articles［EB/OL］.http：//www.cheresources.com/colebrook3.shtml，2011-12-30/2012-01-05.

［7］ 鐘爾杰， 黃延祝. 數值分析 ［M］. 北京： 高等教育出版社，2004：1-3.

［8］ 苑偉民，賀 三，袁宗明，等. 求解BWRS 方程中密度根的數值方法［J］. 天然氣與石油，2009，27（1）：4-6.

［9］ 苑偉民，賀 三，袁宗明，等. 求解BWRS 方程中壓縮因子的數值方法［J］. 管道技術與設備，2009，16（3）：14-16.

［10］ 苑偉民，青 青，袁宗明，等. Colebrook-White 方程顯式公式對比研究［J］. 天然氣與石油，2010，28（4）：5-7.

［11］ 苑偉民，賀 三，袁宗明，等.VB 與MATLAB 混合編程在求解天然氣物性參數中的運用［J］. 中國科學論壇，2008，8（9）：21-23.

［12］［S.N.］.Colebrook Equation［EB/OL］.http：//www.engineeringtoolbox. com/colebrook -equation -d_1031.html，2011 -11-03/2011-08-05.

［13］ Wikipedia. org. Darcy Friction Factor formulae ［EB/OL］.http：//en.wikipedia.org/wiki/Colebrook -White_equation，2011-01-26/2011-07-05.

［14］ Corless R M，Jeffrey D J. Wright Omega Function［EB/OL］.http：//www.maplesoft.com/support/help/Maple/view.aspx?path=Wrightomega，2011-11-03/2011-08-06.

［15］ The Math works. Lambert W Function ［EB/OL］. http：//www.mathworks.com/help/toolbox/mupad/stdlib/wright.html，2011-11-03/2011-08-06.