李奧典 韋 衡

（河海大學(xué)水利水電學(xué)院 江蘇南京 210098）

0 引言

在現(xiàn)行的河流流速測量規(guī)范中,通常采用流速儀和浮標(biāo)測量河流流速[1]，前者費時費工，極不方便。后者受外界的影響較大，精度較低，同樣不是十分的理想。為此，有學(xué)者提出了采用浮子式積深法測量河流流量的設(shè)想[2]，在河底布置氣泡發(fā)生裝置，利用氣泡發(fā)生裝置產(chǎn)生氣泡作為浮子，通過測得浮子從河底升至水面的時間和順流而下的水平距離，計算出河流垂線上的平均流速，從而達到測量流量的目的。楊玉閣根據(jù)實測資料擬合出氣泡冒出時間和水深之間的關(guān)系曲線[3]，論證了氣泡作為浮子在小水深、低流速河流中測量流量的可行性，但缺乏理論依據(jù)。在氣泡上升速度的研究中，程文給出了圓形氣泡上升過程中的運動學(xué)方程[4]，并計算最終速度；繳健研究了氣泡上升過程中加速階段的情況[2],指出了氣泡上升過程中加速階段可以忽略；徐炯討論了氣泡的動力理論特性[5]，得到穩(wěn)定狀態(tài)下不同尺寸單個氣泡的速度變化規(guī)律、運動軌跡規(guī)律等各種運動特性,指出了當(dāng)氣泡大于4 mm時，氣泡在運動過程中開始變形，變成扁橢球形氣泡。對于氣泡在液體中上升運動規(guī)律問題，在石油化工、船舶制造等領(lǐng)域應(yīng)用比較廣泛，并由此建立了氣泡動力學(xué)[6]。筆者在前人研究基礎(chǔ)上分析了水溫對不同形態(tài)氣泡運動速度的影響，并自行設(shè)計了一個實驗裝置，驗證了理論上的球體形態(tài)氣泡在不同溫度中運動的可靠性，同時給出了流速測量所用的氣泡發(fā)生器的氣流量大小和針頭大小。

1 水溫對氣泡運動最終速度的影響

1.1 氣泡冒出最終速度

氣泡在水中的形狀因直徑的大小而不同，主要有三種形狀：球狀、長軸水平的橢圓球形和球帽狀[7]，本文只討論直徑在0～6 mm氣泡的運動情況，直徑在0～4 mm氣泡按球狀氣泡處理，直徑在4～6 mm氣泡按橢圓球形氣泡處理。不同形狀氣泡的最終速度計算公式不盡相同。

1.1.1 直徑在4 mm內(nèi)氣泡的最終速度

根據(jù)Stoke理論分析受力情況，可以將氣泡近似地看作剛性球體，直徑在4 mm內(nèi)的氣泡形態(tài)近似為球體。本文以過渡區(qū)的直徑為1 mm和2 mm氣泡以及紊流區(qū)的直徑3 mm和4 mm氣泡為例，設(shè)定氣泡運動的溫度范圍在0℃至40℃之間,以20℃為基準(zhǔn)溫度，每隔5℃(即0℃、5℃、10℃、15℃、20℃、25℃、30℃、35℃、40℃)作為一個推算點（共9個）。繳健給出了這種情況下水中氣泡上升的最終速度公式[2],如公式(1)所示，式中ρ1為液體密度 (即水的密度)、ρg為空氣密度、dp為氣泡直徑、g為重力加速度（取9.81 m/s2）、CD為阻力系數(shù)。其中，ρ1在各個溫度推算點的對應(yīng)值[8]分別是 999.8、1000.0、999.7、999.1、998.2、997.0、995.7、994.1、992.2 kg/m3；ρg可以利用氣態(tài)方程求出，各個溫度推算點的值分別是 1.293、1.270、1.247、1.225、1.205、1.184、1.165、1.146、1.128 kg/m3。

公式(1)中CD可以通過公式求解，但不同流態(tài)的計算公式不同，繳健采用試算的方法給出了不同直徑的氣泡所處的最終流態(tài)[2]，在基準(zhǔn)溫度下，經(jīng)試算得出：直徑為0.5～2.0 mm的氣泡流態(tài)處于過渡區(qū)，此時，CD的計算公式如公式(2)所示，式中Re為顆粒雷諾數(shù)、v為水的動力黏度。流體力學(xué)顆粒雷諾數(shù)可以用公式(3)計算[9]。直徑為2.0～4.0 mm的氣泡的流態(tài)處于紊流區(qū)，此時的CD≈0.44[9]。對于非球形顆粒其顆粒雷諾數(shù)取等體積球形顆粒的當(dāng)量直徑[10]。

利用公式(1)、(2)、(3)及相關(guān)參數(shù)值可以計算出直徑為1 mm和2 mm氣泡在各個溫度推算點的最終速度，利用式(1)、(2)及相關(guān)參數(shù)值及CD≈0.44可以計算出直徑為3 mm和4 mm氣泡在各個溫度推算點的最終速度。計算結(jié)果如表1所示。

表1 不同溫度下不同直徑氣泡冒出最終速度

1.1.2 直徑大于4 mm氣泡的最終速度

當(dāng)氣泡直徑大于4 mm時，氣泡形狀發(fā)生變形，徐炯認(rèn)為此時氣泡的形態(tài)可按照扁柱體來處理。本文以直徑為5 mm和6 mm氣泡為例，設(shè)定氣泡運動的溫度范圍也在0℃至40℃之間,也以20℃為基準(zhǔn)溫度，每隔5℃（同上）作為一個推算點。徐炯給出了這種情況下氣泡上升的最終速度求解公式[5]，如公式(4)所示，式中σ為水的表面張力系數(shù)，水的表面張力系數(shù)在推算點對應(yīng)值[8]分別為 0.0762、0.0754、0.0748、0.0741、0.0731、0.0726、0.0718、0.071、0.0701 N/m，其他參數(shù)含義同公式(1)。

利用公式(4)和上述各個相關(guān)參數(shù)值，可以計算出在各個溫度推算點的氣泡最終速度，結(jié)果如表1所示。

1.2 敏感指數(shù)與敏感性分析

水溫對不同直徑氣泡影響程度可以用敏感指數(shù)表示，Ng.E和R.S.Loomis給出了敏感指數(shù)S(sensitivity index)公式[11]，如公式(5)所示，式中n為自變量個數(shù)、Xni為各自變量變化后求得的因變量值、Xci為標(biāo)準(zhǔn)自變量下求得的因變量值、δ為自變量的變化率。取常溫20℃下氣泡上升速度為基準(zhǔn)值，利用此公式可以計算出不同直徑氣泡對溫度的敏感指數(shù)，計算結(jié)果如表1所示。

從表1結(jié)果可以看出：水溫的變化對不同直徑氣泡的最終速度有一定的影響；而對處于紊流區(qū)的氣泡影響較小；直徑小于4 mm的氣泡在過渡區(qū)和紊流區(qū)對水溫的敏感指數(shù)隨直徑的增加而增加，當(dāng)直徑大于4 mm時，氣泡對水溫的敏感指數(shù)為一個相對定值。

2 靜水實驗材料與測驗方法

2.1 試驗材料與裝置

本試驗水箱材料選用6 mm厚的有機玻璃，規(guī)格為30 cm(L)×30 cm(W)×120 cm(H）；在底部安置氣泡發(fā)生器；在頂端開有直徑2 cm的孔（目的是既便于水箱內(nèi)的水接觸空氣，也便于向水箱注水）；頂端孔中懸掛一支溫度計（測定水溫）；水箱的背板固定有標(biāo)尺和黑紙（便于觀察和拍攝水箱內(nèi)的氣泡大?。Ｕ麄€裝置如圖1(右)所示。

圖1 氣泡發(fā)生器實物圖（左）及實驗裝置示意圖（右）

圖2 四種規(guī)格針頭產(chǎn)生的氣泡直徑（左到右依次是0.5 mm、0.7 mm、0.9 mm、1.2 mm）

氣泡發(fā)生器由針頭、護套、底座、進氣連接頭等四個部分組成，針頭豎直向上，固定在有機玻璃支架上，如圖1(左)所示。針頭材料為不銹鋼，有四種規(guī)格，內(nèi)徑/外徑分別為0.3 mm/0.5 mm、0.5 mm/0.7 mm、0.7 mm/0.9 mm、1.0 mm/1.2 mm，長度均為4 cm。為了便于更換不同規(guī)格的針頭，水箱設(shè)計成上、下可分離的兩個部分（高度分別為90 cm、30 cm），用法蘭連接。氣泡發(fā)生器的氣流由氣泵產(chǎn)生，通過橡膠連接管接入，如圖1(右)所示。

2.2 氣泡發(fā)生觀測方法

為了觀測氣泡發(fā)生的條件，在針頭進氣管上安裝了氣體流量計和氣壓表測試適宜的流量與壓力，如圖1所示。因氣體流量較小，本實驗采用靈敏度較高的LZB-3WB型轉(zhuǎn)子流量計進行測量氣流?？紤]在室內(nèi)試驗，氣體壓力差不大，故采用U形壓差式氣壓計測定壓差。

運用高速拍攝的方法記錄氣泡形狀及大小，由于在拍攝中需要捕捉到運動中的氣泡，相機的快門速度不宜過低，本實驗的快門速度設(shè)定在1/3200～1/2500 s范圍內(nèi)，光圈和ISO根據(jù)光源強度調(diào)整，鏡頭中心距離水箱底板的高度為70 cm。

2.3 氣泡上升速度的測量

試驗在水深120 cm、100 cm、80 cm等三種環(huán)境下進行。分別試驗四種規(guī)格的針頭及其不同出氣流量、氣水壓差所產(chǎn)生的氣泡形態(tài)。通過實測數(shù)據(jù)推算出氣泡的最終速度，結(jié)果如表2所示。

2.4 氣泡上升形態(tài)大小的測量與分析

四種不同規(guī)格針頭產(chǎn)生的氣泡所拍攝的照片如圖2所示，從照片中可以看出外徑為0.7 mm和0.5 mm的針頭所產(chǎn)生的氣泡大小在3 mm左右，形態(tài)呈圓形或略扁的橢圓；外徑為0.9 mm和1.2 mm的針頭所產(chǎn)生的氣泡直徑在5 mm左右，形態(tài)呈較扁的橢圓，有少部分呈現(xiàn)上突下凹的球帽狀。

對比上述理論計算出的氣泡上升最終速度，當(dāng)氣泡形態(tài)呈球體（針頭的外徑為0.5 mm和0.7 mm）時理論值與實際值吻合程度非常高，誤差只在6%左右，當(dāng)氣泡的形態(tài)呈非球體（針頭的外徑為0.9 mm和 1.2 mm）時，采用扁柱體求得最終速度的理論值與實測數(shù)據(jù)相差較大，可見氣泡在大直徑下的實際運動狀況較為復(fù)雜。

2.5 氣壓、流量對氣泡上升形態(tài)的影響

針對外徑分別為0.7 mm、0.9 mm、1.2 mm等三種規(guī)格的針頭，試驗不同的出氣流量、氣水壓差所產(chǎn)生的氣泡形態(tài)。在實驗中發(fā)現(xiàn)：

（1）外徑0.7 mm針對的適宜氣水壓差可取7～25 cmH2O；外徑0.9 mm針對的適宜氣水壓差可取5～14 cmH2O；外徑1.2 mm針對的適宜氣水壓差可取2～3 cmH2O。可見，對于同一種規(guī)格針頭不同水深的情況，適宜氣水壓差之間的差別不大；但是對于同一水深不同規(guī)格的針頭，適宜的氣水壓差之間的差別卻很大。

（2）外徑0.7 mm針頭的適宜出氣流量可取0.035～0.100 L/min；外徑0.9 mm針頭的適宜出氣流量可取0.035～0.180 L/min；外徑1.2 mm針頭的適宜出氣流量可取0.045～0.170 L/min?？梢?，對于不同規(guī)格的針頭、不同的水深，適宜出氣流量之間的差別不大

表2 氣泡流速計算表

3 結(jié)論

（1）在氣泡上升過程中，溫度對于氣泡的影響較小，在實際測量過程中，溫度對氣泡上升的最終速度的影響基本可以忽略。

（2）采用氣泡浮子法測河流流量的過程中，氣泡發(fā)生器上所選用的針頭內(nèi)徑需小于0.5 mm，內(nèi)徑大于0.5 mm時，由于氣泡實際運動性質(zhì)發(fā)生了變化，運動狀態(tài)較為復(fù)雜，與現(xiàn)有模型[5]計算結(jié)果出入較大，不宜采用。

（3）供氣裝置提供的氣壓對產(chǎn)生的氣泡浮子形態(tài)的影響較為敏感，但水深對氣水壓差的影響不大，實際的出氣量大小以現(xiàn)場測定為準(zhǔn)。

［1］中華人民共和國水利部.GB 50179-1993,河流流量測驗規(guī)范［S］.北京：中國計劃出版社，1993.

［2］繳健，等.基于水動力學(xué)計算的氣泡上升規(guī)律研究［J］.水利信息化,2011(3):41-45.

［3］楊玉閣,陳祖輝.氣泡上升規(guī)律的研究與應(yīng)用［J］.東北水利水電,2010(10):47-48.

［4］程文,周孝德,郭瑾瓏,等.水中氣泡上升速度的實驗研究［J］.西安理工大學(xué)學(xué)報,2000(01):57-60.

［5］徐炯,王彤,楊波,等.靜止水下氣泡運動特性的測試與分析［J］.水動力學(xué)研究與進展A輯,2008(06):709-714.

［6］A.Shima.Studies on bubble dynamics［J］.Shock Waves,1997(1):33-42.

［7］王獻孚.高等流體力學(xué)(第一版)［M］.武漢:華中科技大學(xué)出版社,2003:209-215.

［8］謝振華.工程流體力學(xué)(第三版)［M］.北京:冶金工業(yè)出版社,2007:5.

［9］謝振華.工程流體力學(xué)(第三版)［M］.北京:冶金工業(yè)出版社,2007:170-173.

［10］鐘秦等.化工原理［M］.北京:國防工業(yè)出版社,2009:106.

［11］LiangB.C，Gregorich E.G，Mackenzie A.F.Modelingthe effects ofinorganic and organic amendments on organic matter in a Quebec soil［J］.Soil Science,1996，161(2):109-114.