圓柱體繞流環(huán)隙流場軸向流速及水力損失試驗研究

圓柱體繞流環(huán)隙流場軸向流速及水力損失試驗研究

張琪琦1，孫西歡2，李永業(yè)1

(1.太原理工大學(xué) 水利科學(xué)與工程學(xué)院，太原030024；2.山西水利職業(yè)技術(shù)學(xué)院，運城044004)

摘要：為進一步完善圓柱體繞流環(huán)隙流場理論，采用理論分析和模型試驗相結(jié)合的方法，對不同流量及環(huán)隙比條件下圓柱體繞流環(huán)隙流場進行了研究，得出其軸向速度及水力損失隨測點距離管道中心距離的變化關(guān)系曲線，并給出環(huán)隙水流軸向速度及水力損失分布曲線以及隨流量及環(huán)隙比變化的相關(guān)擬合公式。試驗結(jié)果與擬合公式基本一致，相關(guān)系數(shù)大部分在0.95以上，說明擬合公式是可行的。試驗結(jié)果為圓柱體繞流環(huán)隙流場理論研究及實際應(yīng)用提供了支撐。

關(guān)鍵詞：圓柱體；環(huán)隙流場； 軸向流速；水力損失；流量；環(huán)隙比

中圖分類號：TV134文獻標志碼：A

收稿日期：2013-11-01； 修回日期：2014-02-21

作者簡介：張志昌(1954-)，男，陜西西安人，教授，從事水工水力學(xué)方面的研究，(電話)13991942265(電子信箱)zhangzhichang1954@163.com。

DOI:10.3969/j.issn.1001-5485.2015.04.012

收稿日期：2013-11-27; 修回日期：2013-12-31

基金項目：國家水體污染控制與治理科技重大專項(2012ZX07102-001);江蘇省水利廳科技項目(2012086)

作者簡介：解清杰(1973-), 男, 河北獻縣人，教授,博士，主要從事水污染控制與工程方面的研究，(電話)15951289455(電子信箱)xieqingjie73@163.com。

DOI:10.3969/j.issn.1001-5485.2015.04.011

1研究背景

圓柱體繞流環(huán)隙流場問題是一種常見的技術(shù)問題，如流體在活塞與氣缸間或軸與軸承間的流動等，都存在上述問題。目前許賢良等[1]針對動邊界速度較小時沿軸線方向運動產(chǎn)生的環(huán)隙流動作了初步研究；李永業(yè)等[2]針對因圓柱體的速度、縫隙寬度以及流量對環(huán)隙流所造成的分布及大小的影響作了相關(guān)研究；胡仁喜等[3]對沿周線方向產(chǎn)生的動邊界同心環(huán)隙流也作了相關(guān)試驗研究。

本文主要針對不同直徑的圓柱體在不同水流量情況下，對縫隙流場的軸向流速特性及對應(yīng)的水力損失[4]進行試驗研究，并經(jīng)過對比分析，得出圓柱體繞流縫隙流場軸向流速及水力損失隨流量、環(huán)隙比變化的擬合公式，以完善圓柱體繞流環(huán)隙流場理論。

2試驗測試內(nèi)容與測試系統(tǒng)及儀器儀表

2.1　試驗測試系統(tǒng)

圖1　試驗裝置示意圖 Fig.1　Schematic diagram of experimental device

2.2　圓柱體及環(huán)隙流速測定位置

在圓柱體與輸送管內(nèi)壁的縫間，按照極坐標系布置相應(yīng)的測點，每個測點間距為2.5 mm,同時，為了便于比較分析，只選取沿管道中心軸線垂直方向布置的測點進行研究分析，見圖2。圖中1“×”為測點；2為圓柱體，其柱體外徑為d，3為輸送管道內(nèi)壁，其內(nèi)徑為D。

圖2　圓柱體及環(huán)隙流速測點位置示意圖 Fig.2　Section of the cylinder and layout of measurement points of annulus velocity

2.3　試驗測試內(nèi)容

本文共使用4種圓柱體，即長度為100 mm,直徑分別為50，60，70，80 mm，對應(yīng)的相對環(huán)隙比記為W=(D-d)/D，其對應(yīng)值分別為0.5，0.4，0.3，0.2。輸送水流量Q分別為30，40，50，60，70 m3/h共5種。測點x值表示該點距離管道中心的距離，單位mm。

3環(huán)隙水流軸向流速特性

3.1　環(huán)隙流場軸向流速 V Z隨流量 Q的變化關(guān)系

在圓柱體直徑為60 mm的情況下，對流量為40,50,60,70 m3/h的圓柱體中部斷面的環(huán)隙流場進行軸向速度VZ測試,環(huán)隙流場軸向速度值在不同流量下的關(guān)系曲線如圖3所示。圖3中，以管道中心為起點，測點距離管道中心的距離為橫坐標，在實際測試中，起始測點橫坐標為32 mm，即距離管道車外壁2 mm，然后以每2.5 mm依次向外布置測點，最后1個測點橫坐標為49.5 mm，即距離管道內(nèi)壁0.5 mm。圖中軸向速度為無量綱數(shù)，為該測點系列軸向速度值除以該系列最大軸向速度值。

圖3　不同流量下環(huán)隙水流軸向流速測試 Fig.3　Curves of axial velocity of annulus flow field in the presence of different flow rates

從圖3中可以看出，在測試斷面內(nèi)，以中間點為對稱軸的左右兩側(cè)成對數(shù)型分布。這主要是因為在環(huán)隙流段區(qū)間，由于圓柱體處于靜止狀態(tài)，根據(jù)液體的黏滯性，管壁及圓柱體壁面處的液體質(zhì)點依附在壁面處，縫隙中部流速取得最大值。而此時測試圓管段內(nèi)最小雷諾數(shù)Re=80 916>2 000,則環(huán)隙流段內(nèi)處于紊流狀態(tài)，液體質(zhì)點的混滲作用使得質(zhì)點間進行動量交換，流速分布呈現(xiàn)均勻化，所以紊流時過水斷面上的流速是按對數(shù)規(guī)律分布的。因此，為便于分析比較，選取其中一半進行擬合。其擬合關(guān)系式見表1。

通過比較和分析可知，環(huán)隙水流軸向速度隨測點距離管道中心距離的變化而變化關(guān)系基本可以擬合成形如VZ=A1lnx+B1的對數(shù)型公式，其中系數(shù)A1、系數(shù)B1均是與流量有關(guān)的關(guān)系式，根據(jù)軸向速度計算公式V=Q/A,在圓柱體直徑和長度不變的情況下，環(huán)隙流段軸向速度與進入環(huán)隙流段的流量成正比的關(guān)系。再對不同流量條件下的系數(shù)A1、系數(shù)B1進行擬合，得出系數(shù)A1、系數(shù)B1隨流量的擬合值，見表2。

表1　環(huán)隙水流在不同流量下軸向速度隨 測點 x擬合關(guān)系式 Table 1　Fitting formulas of axial flow velocity with measurement point x under different flow rates

表2　系數(shù) A 1、系數(shù) B 1在不同流量下擬合值 Table 2　Fitting results of coefficient A 1 and coefficient B 1 under different flow rates

對表1、表2分析可知，系數(shù)A1、系數(shù)B1在不同流量下擬合關(guān)系式分別為A1=0.041 3Q-0.488 8(相關(guān)系數(shù)R2=0.975 8)和B1=-0.110 8Q+1.891 1(相關(guān)系數(shù)R2=0.954 9)。所以圓柱體繞流環(huán)隙流場軸向速度隨測點距離管道中心距離及流量的變化關(guān)系基本可擬合成式(1),即

VZ=(0.041 3Q-0.488 8)lnx+

(1)

從式(1)可知，環(huán)隙水流軸向速度隨測點距離管道中心距離的變化而變化，呈現(xiàn)對數(shù)規(guī)律分布，且隨著流量的增加，環(huán)隙流段的軸向流速呈線性增加的趨勢。

（69）喜馬拉雅鞭苔 Bazzania himlayana（Mitt.）Schiffn. 熊源新等（2006）；楊志平（2006）

3.2　環(huán)隙流場軸向流速 V Z隨環(huán)隙比的擬合公式

在流量為50 m3/h的情況下，對直徑分別為50，60，70，80 mm的圓柱體中部斷面的環(huán)隙流場進行軸向速度VZ測試，即環(huán)隙比W依次為0.5，0.4，0.3，0.2。環(huán)隙流場軸向速度值在不同環(huán)隙比情況下關(guān)系曲線如圖4所示。圖中軸向速度為無量綱數(shù)，為該測點系列軸向速度值除以該系列最大軸向速度值。

圖4　不同環(huán)隙比下環(huán)隙水流軸向流速測試 Fig.4　Curve of axial velocity of annulus flow field in the presence of different annulus ratios

從上述分析結(jié)合圖4中可以看出，在測試斷面內(nèi)，以中間點為對稱軸的左右兩側(cè)成對數(shù)型分布。因此，為便于分析比較，選取其中一半進行擬合。其擬合曲線及關(guān)系式見表3。

表3　環(huán)隙水流在不同環(huán)隙比下軸向速度隨 測點 x擬合關(guān)系式 Table 3　Fitting formulas of axial flow velocity with measurement point x in the presence of different annulus ratios

同上分析，環(huán)隙水流軸向速度隨測點距離管道中心距離的變化關(guān)系基本可以擬合成形如VZ=A2lnx+B2的對數(shù)型公式，其中系數(shù)A2、系數(shù)B2均是與環(huán)隙比有關(guān)的關(guān)系式。因為環(huán)隙比的增加，環(huán)隙斷面面積減小，在流量一定時，軸向流速增大，根據(jù)軸向速度計算公式V=Q/A，流量一定，環(huán)隙斷面的面積A與環(huán)隙比W的關(guān)系成A=W-2的相關(guān)關(guān)系。再對不同環(huán)隙比條件下的系數(shù)A2、系數(shù)B2進行擬合，得出系數(shù)A2、系數(shù)B2隨環(huán)隙比的擬合值，見表4。

表4　系數(shù) A 2、系數(shù)B 2在不同環(huán)隙比下擬合值 Table 4　Fitting results of coefficient A 2 and coefficient B 2 in the presence of different annulus ratios

根據(jù)表3、表4的分析可知，系數(shù)A2、系數(shù)B2在不同環(huán)隙比下擬合關(guān)系式分別為A2=1.31W-2-1.927 8W-1-1+1.045 8(相關(guān)系數(shù)R2=0.955 6)和B2=2.055W-2-3.858 5W-1+2.391 3(相關(guān)系數(shù)R2=0.892 6)。所以圓柱體繞流環(huán)隙流場軸向速度隨測點距離管道中心距離及環(huán)隙比W的變化關(guān)系基本可以擬合成式(2)，即

(2)

4環(huán)隙水流的主要水力損失研究

環(huán)隙水流的水力損失主要包括圓柱體繞流進口端出現(xiàn)的水流急劇收縮、環(huán)隙段內(nèi)水力損失及圓柱體繞流出口端出現(xiàn)的水流急劇擴散等過程。本文中，在進口端前10 cm及出口端10 cm設(shè)置測壓管，旨在研究圓柱體繞流的局部主要水力損失，并得出環(huán)隙水流主要水力損失與流量Q及環(huán)隙比的擬合關(guān)系式。

4.1　環(huán)隙水流主要水力損失與輸水流量的關(guān)系

二者關(guān)系見圖5。圖5中的公式給出了不同圓柱體在不同流量時的單位水頭損失,圖中曲線從上往下分別代表環(huán)隙比為0.2，0.3，0.4，0.5，圖中公式從上往下分別代表環(huán)隙比為0.2，0.3，0.4，0.5下環(huán)隙水流主要水力損失隨輸水流量變化的數(shù)學(xué)擬合公式，其中橫坐標為流量，縱坐標為單位水頭損失。由圖5可見，總體上水力損失隨流量增加而增大，并呈現(xiàn)乘冪函數(shù)型增長趨勢。這主要是因為，根據(jù)水力學(xué)知識，任何一種局部水頭損失通常都可以用一個系數(shù)和流速水頭的乘積來表示：ΔE=§Vm/(2g),而在環(huán)隙斷面面積不變的前提下，環(huán)隙水流軸向流速與流量成線性相關(guān)的關(guān)系，所以水頭損失與流量成乘冪函數(shù)型相關(guān)增長趨勢。

圖5　環(huán)隙水流水力損失與流量關(guān)系曲線 Fig.5　Curves of hydraulic energy consumption of annulus flow field vs. flow rate

其中，在相同流量條件下，長和直徑為100 mm×60 mm的圓柱體，即環(huán)隙比為0.4時對應(yīng)的單位水頭損失最小。這主要是因為環(huán)隙比的增大，水流的過水斷面面積增大，當環(huán)隙流量的增量大于過水斷面面積的增量時，環(huán)隙流速就隨環(huán)隙比增大而增大，反之則減小。由于環(huán)隙比W<0.4時，環(huán)隙流量的增量大于過水斷面面積增量，環(huán)隙流速隨環(huán)隙比增加呈現(xiàn)上升趨勢，而當環(huán)隙比W>0.4后，環(huán)隙流量的增量小于面積增量，從而導(dǎo)致環(huán)隙流速呈現(xiàn)減小的趨勢,所以環(huán)隙比W=0.4時，同流量下環(huán)隙流軸向速度最大，則水力損失最小。

4.2　環(huán)隙水流主要水力損失與環(huán)隙比 W的關(guān)系

圖6給出了不同環(huán)隙比和不同流量條件下的單位水頭損失及相關(guān)擬合關(guān)系，圖中曲線從上往下分別代表流量為70，60，50，40，30 m3/h下環(huán)隙水流主要水力損失隨環(huán)隙比變化的關(guān)系及數(shù)學(xué)擬合公式。不同環(huán)隙比水力損失亦不同，總的趨勢是環(huán)隙越小，水頭損失越大；在相同環(huán)隙比條件下，流量小，水力損失亦小，這與前文規(guī)律是一致的。

圖6 環(huán)隙水流水力損失與環(huán)隙比的關(guān)系Fig.6 Curvesofhydraulicenergyconsumptionofannulusflowfieldvs.annulusratio

5結(jié)論

(1) 圓柱體繞流環(huán)隙水流軸向速度隨測點距離管道中心距離的變化而變化基本呈現(xiàn)對數(shù)規(guī)律分布，可以擬合成形如VZ=Alnx+B的對數(shù)型公式。其中：軸向速度與流量變化的關(guān)系式為：VZ= (0.041 3Q-0.488 8)lnx+(1.891 1-0.110 8Q)，且隨著流量的增加，環(huán)隙流段的軸向流速呈線性增加的趨勢；軸向速度隨環(huán)隙比變化的關(guān)系式為VZ=(1.31W-2-1.927 8W-1-1+1.045 8)lnx+(2.055W-2-3.858 5W-1+2.391 3)，且隨著環(huán)隙比的增加，環(huán)隙流段的軸向流速呈負冪函數(shù)的趨勢。

(2) 圓柱體繞流環(huán)隙流場的水力損失與環(huán)隙比和輸水流量均有一定的關(guān)系：環(huán)隙比一定時，總體上水力損失隨流量增加而增大，并呈現(xiàn)乘冪函數(shù)型增長趨勢；流量一定時，水力損失隨著環(huán)隙比增加呈現(xiàn)先減小后增加的趨勢：當環(huán)隙比為0.4的圓柱體對應(yīng)的單位水頭損失最小。

參考文獻：

[1]許賢良，朱兵，張軍，等． 同心環(huán)形縫隙流理論研究[J]．安徽理工大學(xué)學(xué)報:自然科學(xué)版，2004，24(3):40-42．(XU Xian-liang，ZHU Bing，ZHANG Jun，etal． The Research of the Concentric Ring Slit Flow[J]. Journal of Anhui University of Science and Technology: Natural Science, 2004, 24( 3) : 40-42．(in Chinese))

[2]李永業(yè)，孫西歡，李飛，等.動邊界同心環(huán)狀縫隙流研究[J].農(nóng)業(yè)機械學(xué)報,2012, 43(3):231-234.(LI Yong-ye, SUN Xi-huan, LI Fei,etal. Cyclical Slit Flow of Concentricity under the Moving Boundary Condition[J]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(3):231-234. (in Chinese))

[3]胡仁喜，苑士華，劉紅寧，等．高壓高速條件下柱塞副泄漏流場分析[J]．農(nóng)業(yè)機械學(xué)報，2009，40(4):221-226． (HU Ren-xi，YUAN Shi-hua，LIU Hong-ning，etal． Analysis on the Leaking Flow Field of the Piston Sector Considering the High Press and High Velocity[J]. Transactions of the Chinese Society for Agricultural Machinery, 2009, 40(4): 221-226．(in Chinese))

[4]林愉.繞靜止圓柱體型料的管道流水力損失[J].農(nóng)業(yè)機械學(xué)報，2008,39(12)：93-96. (LIN Yu. Hydraulic Energy Consumption in Annulus Flow Field of the Flow around Cylinder[J]. Transactions of the Chinese Society for Agricultural Machinery, 2008,39(12)：93-96.(in Chinese))

[5]閻慶紱，張曉東，武利生，煤錠管道輸送的平均能耗[J].水力采煤與管道輸送，2003,29(3)：36-39. (YAN Qing-fu, ZHANG Xiao-dong, WU Li-sheng. The Average Consumption of Coal Log Piped Transportation[J]. Hydraulic Coal Mining & Pipeline Transportation, 2003,29(3)：36-39. (in Chinese))

[6]李永業(yè)，孫西歡，延耀興.管道車不同荷重時筒裝料管道水力輸送特性[J].農(nóng)業(yè)機械學(xué)報，2008,39(12)：93-96. (LI Yong-ye，SUN Xi-huan，YAN Yao-xing． Hydraulic Characteristics of Tube-contained Raw Material Hydraulic Transportation under the Different Loads on the Piped Carriage[J]. Transactions of the Chinese Society for Agricultural Machinery, 2008，39( 12) : 93-96．(in Chinese))

[7]童學(xué)衛(wèi)，李永業(yè)，孫西歡.筒裝料管道水力輸送試驗系統(tǒng)設(shè)計研究[J].山西水利，2008，(2)：42-43. (TONG Xue-wei, LI Yong-ye, SUN Xi-huan． Experimental System Design on the Piped Hydraulic Transportation of Tube-contained Raw Material[J]. Shanxi Water Resources，2008, (2) : 42 -43．(in Chinese))

[8]毛寧. 關(guān)于求解力學(xué)中常用關(guān)系曲線的經(jīng)驗公式[J]. 長江科學(xué)院院報, 2009, 26(10): 54-57. (MAO Ning. How to Find Empirical Formulas of Usual Relation Curves in Mechanics[J]. Journal of Yangtze River Scientific Research Institute, 2009, 26(10): 54-57. (in Chinese))

(編輯：王慰)

Experimental Research on Axial Flow Velocity and HydraulicEnergy Consumption in Annulus Flow Field around Cylinder

ZHANG Qi-qi1,SUN Xi-huan2, LI Yong-ye1

(1.School of Water Conservancy Sciences and Engineering, Taiyuan University of Technology,

Taiyuan030024, China; 2.Shanxi Water Conservancy Technical College, Yuncheng044004, China)

Abstract:To further improve the theory of annulus flow field of the flow around cylinder, we combined theoretical analysis with model test to research the annulus flow field around cylinder in the presence of different flow rates and annulus ratios. The variations of axial velocity and hydraulic loss along with the distance of measurement point to the cylinder center are obtained, and the fitting formulas describing the variations of axial velocity and hydraulic loss with flow rate and annulus ratio are given. The test results are consistent with the fitting formulas, and most of the correlation coefficients are above 0.95, indicating that the fitting formulas are feasible. The result is a test support for the theoretical research and practical application of the annulus flow field around cylinder.

Key words: cylinder; annulus flow field; axial flow velocity; hydraulic energy consumption; flow rate; annulus ratio

2015,32(04):59-64

2015,32(04):55-58