王晶，王亦飛，顏留成，伍子瑋，于廣鎖（華東理工大學(xué)潔凈煤技術(shù)研究所，煤氣化及能源化工教育部重點(diǎn)實(shí)驗(yàn)室，上海 200237）

?

管內(nèi)垂直下降液膜速度與厚度分布特性

王晶，王亦飛，顏留成，伍子瑋，于廣鎖
（華東理工大學(xué)潔凈煤技術(shù)研究所，煤氣化及能源化工教育部重點(diǎn)實(shí)驗(yàn)室，上海 200237）

摘要：對(duì)洗滌冷卻管內(nèi)垂直降膜的流動(dòng)特性進(jìn)行研究，采用超聲波多普勒測(cè)速儀對(duì)管內(nèi)不同周向以及軸向位置的液膜厚度和速度進(jìn)行了無接觸式的測(cè)量，液膜Reynolds數(shù)范圍為1.0×104～3.1×104。結(jié)果表明：在0°周向位置上液膜厚度與速度均達(dá)到最大值，導(dǎo)致該位置局部液膜厚度過大而不能保持穩(wěn)定，部分液體脫離液膜表面，此外還造成了8°和16°位置的液膜厚度激增。在軸向上，當(dāng)Reynolds數(shù)小于2.0×104時(shí)，液膜速度在重力作用下隨流動(dòng)距離增加而增加，反之，液膜速度因?yàn)榱鲃?dòng)阻力會(huì)隨距離增加而減小。隨著Reynolds數(shù)的增大，液膜平均厚度和速度呈增大趨勢(shì)。此外，Reynolds數(shù)的增大還會(huì)使得液膜更加不穩(wěn)定。

關(guān)鍵詞：氣液兩相流；湍動(dòng)；膜；速度測(cè)量；厚度測(cè)量；超聲多普勒測(cè)速儀

2015-12-23收到初稿，2016-02-04收到修改稿。

聯(lián)系人：王亦飛。第一作者：王晶（1990—），女，碩士研究生。

Received date:2015-12-23.

Foundation item: supported by the National High Technology Research and Development Program of China (2012AA053101).

引　言

在多噴嘴對(duì)置式水煤漿氣化技術(shù)中洗滌冷卻室是重要的組成部分，它主要由洗滌冷卻環(huán)、洗滌冷卻管、氣相分離空間、液池等組成[1-2]。在洗滌冷卻室內(nèi)，冷卻水經(jīng)洗滌冷卻環(huán)分布后沿洗滌冷卻管內(nèi)壁呈膜狀向下流動(dòng)，同時(shí)與高溫氣體發(fā)生熱質(zhì)交換，使其冷卻并增濕。洗滌冷卻管內(nèi)垂直降膜的流動(dòng)特性對(duì)其與高溫合成氣、熔渣之間的傳熱傳質(zhì)有著重要的影響。除此之外，氣液兩相降膜流動(dòng)還廣泛用于各種工業(yè)領(lǐng)域，如空調(diào)與制冷、核工業(yè)、機(jī)械加工及石化行業(yè)等[3]。

部分學(xué)者借助高速相機(jī)、電導(dǎo)或電容探針等設(shè)備對(duì)圓管內(nèi)或平面上的氣液兩相降膜流動(dòng)進(jìn)行了研究。宋建等[4]借助高速攝像儀和高頻電容探針測(cè)量了豎直平板表面下降液膜表面波的頻率和周期等參數(shù)。王靈萍等[5]采用雙平行電導(dǎo)探針對(duì)豎直圓管內(nèi)下降液膜的厚度進(jìn)行了測(cè)量。閻維平等[6]借助CCD攝像機(jī)對(duì)豎直平板上下降液膜的厚度進(jìn)行測(cè)量，得到了液膜厚度在時(shí)間和空間上的分布情況。韋勝杰等[7]同樣借助電容法測(cè)量了低Reynolds數(shù)下豎直平板上下降液膜的厚度。Abdulkadir等[8]對(duì)180°豎直彎管內(nèi)氣液兩相流的液膜進(jìn)行了研究，采用平行線狀電導(dǎo)探針測(cè)量液膜較厚的區(qū)域（>2.5 mm），用嵌入式針狀電導(dǎo)探針測(cè)量液膜較薄的區(qū)域（<2.5 mm）。

最早的關(guān)于降膜流動(dòng)的研究可以追溯到1916年提出的Nusselt理論，提出的膜狀層流流動(dòng)模型假設(shè)降膜為無限寬薄液膜、具有光滑的自由剪切氣液界面、物性參數(shù)為常數(shù)、流動(dòng)向無傳熱等。當(dāng)壁面與水平方向夾角為90°時(shí)，液膜速度及厚度的計(jì)算式如下

其中液膜厚度δ的表達(dá)式為

Alekseenko等[9]采用電導(dǎo)探針測(cè)量了內(nèi)徑15 mm的豎直圓管內(nèi)的液膜厚度，液膜Reynolds數(shù)范圍10～135，結(jié)果表明表面波的速度對(duì)逆流氣相速度較為敏感。

Koizumi等[10]對(duì)內(nèi)徑30 mm豎直圓管內(nèi)水和硅油兩種介質(zhì)形成的液膜進(jìn)行了研究，以水為液相的液膜Reynolds數(shù)范圍4200～7700，將熒光染料溶于液相，用氬激光束照射液膜使其發(fā)出熒光，用高速攝像機(jī)拍攝圖像，據(jù)此分析計(jì)算液膜的速度等，推導(dǎo)出了最大液膜厚度與波速的關(guān)聯(lián)式為

此外，也有學(xué)者通過數(shù)值模擬方法對(duì)下降液膜的流動(dòng)特性[11-12]和傳熱傳質(zhì)特性[13-15]進(jìn)行了研究。

目前，對(duì)于降膜流動(dòng)問題的研究多集中于中低Reynolds數(shù)條件下（Re<1.0×104），同時(shí)鮮有學(xué)者通過實(shí)驗(yàn)方法對(duì)液膜的速度進(jìn)行測(cè)量。本工作對(duì)洗滌冷卻管氣液兩相垂直下降液膜的特性進(jìn)行了研究，液膜Reynolds數(shù)最高達(dá)3.1×104。借助超聲多普勒測(cè)速儀對(duì)液膜的厚度和速度進(jìn)行測(cè)量，獲得了不同Reynolds數(shù)條件下液膜的平均厚度和平均速度的分布情況，研究其分布特性，為降膜流動(dòng)在洗滌冷卻室內(nèi)的應(yīng)用提供理論指導(dǎo)。

本工作采用的超聲多普勒測(cè)速儀（ultrasound Doppler velocimetry, UDV）是一款基于超聲波探測(cè)技術(shù)測(cè)量流體速度的儀器，可以用于液態(tài)金屬流動(dòng)測(cè)量[16-17]、固體顆粒運(yùn)動(dòng)測(cè)量[18]、氣泡運(yùn)動(dòng)測(cè)量[19]和固含率測(cè)量[20]等。相比電導(dǎo)法、電容法等接觸式方法，非接觸式測(cè)量具有受限小、不擾亂流場(chǎng)等特點(diǎn)。此外，由于超聲波有較強(qiáng)的穿透性，對(duì)裝置的可視化程度沒有要求，可以穿過不透明壁面進(jìn)行測(cè)量。

1　實(shí)驗(yàn)裝置與方法

圖1　實(shí)驗(yàn)裝置Fig.1　Schematic diagram of experimental apparatus1—UDV; 2—computer; 3—flow-meter; 4—water pump; 5—water tank; 6—scrubbing-cooling ring; 7—water inlet; 8—scrubbing-cooling tube; 9—transducer; 10—chamber

實(shí)驗(yàn)裝置如圖1所示，洗滌冷卻管內(nèi)徑為0.138 m，高度為1 m。實(shí)驗(yàn)介質(zhì)液相為水，氣相為空氣。水由水泵輸送至洗滌冷卻環(huán)，洗滌冷卻環(huán)上共有4個(gè)對(duì)稱的入水口，每個(gè)入水口進(jìn)水流量相同。經(jīng)洗滌冷卻環(huán)分布后，水以液膜的形式沿洗滌冷卻管內(nèi)壁豎直向下流動(dòng)。為了便于闡述，定義周向角度θ，選取入水口所在位置為0°，如圖1所示。定義軸向高度H，洗滌冷卻環(huán)出口處H=0，洗滌冷卻管出口處H=1 m。洗滌冷卻環(huán)出口處液膜平均速度范圍為0.87～2.62 m·s?1，液膜Reynolds數(shù)Re=1.0×104～3.1×104，其計(jì)算式為Re=4Γ/μ，Γ為單位濕周上的質(zhì)量流量。

實(shí)驗(yàn)采用超聲多普勒測(cè)速儀測(cè)量液膜的速度與厚度。因?yàn)閷?shí)驗(yàn)裝置是中心對(duì)稱的，所以測(cè)量時(shí)只在0°～45°范圍內(nèi)選擇測(cè)量點(diǎn)。在進(jìn)行速度測(cè)量時(shí)需要在超聲波探頭外加套管，受裝置尺寸限制，選取了0°、22.5°和45° 3個(gè)周向位置，每個(gè)位置軸向上共取14個(gè)測(cè)量點(diǎn)。厚度測(cè)量選取了0°、8°、16°、24°、34.5°和45° 6個(gè)周向位置，每個(gè)位置軸向上取14個(gè)測(cè)量點(diǎn)。

圖2　速度與厚度的測(cè)量Fig. 2　Measurement of velocity and thickness of liquid film

圖2為超聲多普勒測(cè)速儀測(cè)量速度與厚度的示意圖。測(cè)量流體速度的原理[圖2(a)]基于多普勒效應(yīng)，超聲波探頭向待測(cè)流體發(fā)出一定頻率的超聲波，穿過傳播介質(zhì)以及管壁后散播在流體介質(zhì)中。超聲波傳播方向與流體方向成一定夾角（即多普勒角），當(dāng)遇到介質(zhì)中的固體顆粒、氣泡等示蹤粒子時(shí)反射波的頻率會(huì)發(fā)生改變，探頭接收到不同深度粒子反射波的時(shí)間不同，經(jīng)過信號(hào)動(dòng)態(tài)頻譜分析和計(jì)算便可得到超聲波發(fā)射方向上不同深度的流體速度。在相界面處，因?yàn)槁暡▽挾炔荒芎雎圆挥?jì)，所以測(cè)量得到的速度需要進(jìn)行校正，具體方法可以參考Nowak[21]的工作。圖2(b) 為測(cè)量液膜厚度的原理。超聲波多普勒測(cè)速儀能夠獲得反射波的頻率和強(qiáng)度，而超聲波在穿過兩相界面時(shí)會(huì)有較強(qiáng)的反射，反射波的強(qiáng)度會(huì)出現(xiàn)一個(gè)峰值，根據(jù)發(fā)射與接收聲波信號(hào)的時(shí)間差可以計(jì)算發(fā)射反射的界面所在的位置。本實(shí)驗(yàn)有壁面-液膜和液膜-氣體兩個(gè)相界面，其對(duì)應(yīng)的兩個(gè)較強(qiáng)反射信號(hào)峰之間的距離即為液膜的厚度，并且超聲波自身衰減速度很快，經(jīng)過空氣介質(zhì)會(huì)很快衰減，所以在測(cè)量時(shí)不需考慮超聲波傳遞到對(duì)側(cè)液膜帶來的影響。本實(shí)驗(yàn)所用的傳感器探測(cè)區(qū)域的直徑為5 mm，所發(fā)出的超聲波頻率為4 MHz，采樣速度為每秒100組數(shù)據(jù)。本實(shí)驗(yàn)條件下，超聲多普勒測(cè)速儀厚度測(cè)量的誤差為0.06 mm，速度測(cè)量的誤差為6 mm·s?1。

測(cè)量的相對(duì)標(biāo)準(zhǔn)不確定度為

儀器相對(duì)標(biāo)準(zhǔn)不確定度為

合成不確定度為

液膜厚度的不確定度在0°位置最大，約為2.5%；速度的不確定度最大約1.5%。

2　結(jié)果與討論

2.1液膜速度分布

圖3為不同位置不同Reynolds數(shù)條件下平均液膜厚度分布。實(shí)驗(yàn)測(cè)量了6個(gè)不同周向位置的液膜厚度分布。由圖3可見液膜厚度在0°位置出現(xiàn)最大值。當(dāng)Reynolds數(shù)為1.0×104時(shí)，可以觀察到液膜最厚的位置出現(xiàn)在H≈0.2 m處，隨后液膜厚度在重力作用下開始降低。然而，當(dāng)Reynolds數(shù)增大時(shí)，液膜厚度并沒有明顯的增厚，而是隨流動(dòng)距離增加其大小保持基本不變，這與王靈萍等[5]使用雙平行電導(dǎo)探針測(cè)量得到的結(jié)果不同。其原因主要是當(dāng)液膜Reynolds數(shù)大于10182時(shí)0°位置出現(xiàn)了部分液體脫離液膜進(jìn)入氣相的現(xiàn)象，采用電導(dǎo)探針測(cè)量厚度時(shí)已經(jīng)脫離液膜的液體仍會(huì)經(jīng)過電導(dǎo)探針，因此被當(dāng)成液膜的一部分，從而使得測(cè)量得到的液膜厚度偏大。產(chǎn)生該現(xiàn)象的主要原因是入水口恰好也在0°位置，雖然洗滌冷卻環(huán)已使液體分布盡可能均勻，但0°位置仍然會(huì)有液體積聚，使得該位置液膜較厚。當(dāng)液體流量增大時(shí)，局部液膜越來越厚，液體的內(nèi)聚力不足以維持液膜的穩(wěn)定性，因此部分液體脫離了液膜，在氣相中繼續(xù)下落。

圖3　不同Reynolds數(shù)條件下液膜厚度分布Fig.3　Liquid film thickness distribution at different Reynolds numbers

周向角度為8°位置的液膜與0°位置相比較薄。當(dāng)Reynolds數(shù)為1.0×104時(shí)，液膜厚度隨流動(dòng)距離增加而減小。當(dāng)Reynolds數(shù)較大時(shí)，液膜厚度在H=0.4 m左右出現(xiàn)了突然增加的現(xiàn)象，隨后液膜厚度持續(xù)減小，而且Reynolds數(shù)越大出現(xiàn)拐點(diǎn)的位置越靠后。這主要是受0°位置局部液膜狀態(tài)影響。0°位置局部液膜過厚，因此部分液體會(huì)向兩側(cè)較薄的區(qū)域流動(dòng)。當(dāng)Reynolds數(shù)為10182時(shí)，來自0°位置的液體造成的液膜厚度的增加不足以抵消重力造成的厚度減小，因此并沒有出現(xiàn)明顯的厚度激增情況。相比而言，16°位置的液膜比8°位置的液膜薄得多，因此液膜厚度激增的現(xiàn)象也更加明顯。此外，因?yàn)?6°位置離0°位置更遠(yuǎn)，液膜厚度激增出現(xiàn)的位置也更加靠后。

如圖3所示，24°、34.5°和45°位置的液膜厚度隨周向角度增加而減小。不同液膜Reynolds數(shù)的液膜厚度分布趨勢(shì)大致相同。在管上半部分，液膜厚度隨流動(dòng)距離增加而略微減小，管下半部分趨勢(shì)剛好相反。總體而言，這3個(gè)位置的液膜厚度小于0°、8°和16°位置的液膜厚度。

圖4為不同Reynolds數(shù)條件下0°位置H=0.22 m處的液膜厚度隨時(shí)間的分布，虛線表示液膜的平均厚度。在低Reynolds數(shù)下液膜厚度隨時(shí)間波動(dòng)幅度較小，當(dāng)Reynolds數(shù)逐漸增大時(shí)，較大幅度的波動(dòng)變得更加頻繁，而且波動(dòng)的幅度也越來越大。總體而言，可以發(fā)現(xiàn)液膜厚度波動(dòng)沒有明顯的規(guī)律；隨著Reynolds數(shù)的增大，液膜厚度波動(dòng)的振幅和頻率均有所增加，也就是說液膜變得更加不穩(wěn)定。

2.2液膜速度分布

圖5為不同位置不同Reynolds數(shù)條件下液膜速度分布。在0°位置，與液膜厚度分布類似，液膜速度相比其他位置大。當(dāng)液膜Reynolds數(shù)較?。≧el<2.0×104）時(shí)，隨著流動(dòng)距離的增加，液膜速度增加，同時(shí)如圖3所示液膜的厚度減小。相反地，當(dāng)液膜Reynolds數(shù)較大（Rel>2.0×104）時(shí)，可以發(fā)現(xiàn)，隨著流動(dòng)距離的增加，液膜速度在管上半段迅速減小，在下半段緩慢增加，對(duì)應(yīng)的液膜厚度一直保持穩(wěn)定。這是因?yàn)樵诠軆?nèi)液膜速度主要受重力和阻力兩個(gè)因素影響。當(dāng)Reynolds數(shù)較小時(shí)，液膜的初始速度較小，重力是關(guān)鍵因素，因此液膜的速度隨流動(dòng)距離增加而增加。相反，當(dāng)Reynolds數(shù)較大時(shí)，阻力的影響使得液膜的速度迅速降低。較大Reynolds數(shù)下速度變化趨勢(shì)的轉(zhuǎn)變可能是由0°位置部分液膜的脫離造成的。

與0°位置相比，22.5°與45°位置的液膜速度具有相似的分布趨勢(shì)，當(dāng)Reynolds數(shù)較小時(shí)液膜速度隨流動(dòng)距離增加而增加，當(dāng)Reynolds數(shù)較大時(shí)液膜速度隨流動(dòng)距離增加而一直減小。這兩個(gè)位置的速度值差距較小，而且均略小于0°位置的速度值?？傮w而言，液膜Reynolds數(shù)越大，液膜的速度也越大。

3　結(jié)　論

借助超聲波多普勒測(cè)速儀對(duì)洗滌冷卻管內(nèi)垂直下降液膜的速度和厚度進(jìn)行非接觸式測(cè)量，分析了液膜厚度和速度的分布規(guī)律，得到如下結(jié)論。

圖4　不同Reynolds數(shù)條件下液膜厚度隨時(shí)間的變化Fig. 4　Thickness profiles of liquid film with time at different Reynolds numbers

圖5　不同位置及不同Reynolds數(shù)條件下液膜速度分布Fig.5　Liquid film velocities distribution at different circumferential positions and Reynolds numbers

（1）0°周向位置平均液膜厚度最大。當(dāng)Reynolds更加不穩(wěn)定。

（3）0°周向位置的液膜速度較其他位置大。當(dāng)數(shù)大于1.0×104時(shí)，該位置局部液膜厚度過大而不能保持穩(wěn)定，部分液體脫離液膜，在氣相中繼續(xù)下落。此外，該位置部分液體會(huì)向兩側(cè)液膜較薄的位置流動(dòng)，造成8°和16°位置的液膜厚度激增。其余位置的液膜厚度較薄且變化較緩。

（2）整體液膜厚度和速度均會(huì)隨Reynolds數(shù)增大而增大。此外，Reynolds數(shù)的增大還會(huì)使得液膜Reynolds數(shù)小于2.0×104時(shí)，液膜速度在重力作用下隨流動(dòng)距離增加而增加。與此相反，當(dāng)Reynolds數(shù)大于2.0×104時(shí)，液膜速度因?yàn)樽枇?huì)隨距離增加而減小。不同周向位置的速度分布趨勢(shì)相似。

符號(hào)說明

g——重力加速度，m·s?2

H——軸向高度坐標(biāo)，m

Ql——液體的流量，kg·s?1

Re——Reynolds數(shù)

s——平板面積，m2

t——時(shí)間，s

u——液膜速度，m·s?1

vl——液膜速度，m·s?1

x——以壁面為零坐標(biāo)、垂直指向液體方向?yàn)檎较虻臋M坐標(biāo)，m

?！獑挝粷裰苌系馁|(zhì)量流量，kg·m?1·s?1

δ——液膜厚度，m

θ——周向角度，(°)

μl——液體的黏度，m2·s?1

ρl——液體的密度，kg·m?3

References

[1]ZHENG L, FURINSKY E. Comparison of Shell, Texaco, BGL and KRW gasifiers as part of IGCC plant computer simulation [J]. Energy Conversion and Management, 2005, 46(11/12): 1767-1779.

[2]YU Z H, WANG Y F, ZHOU Z J, et al. A compound bed high temperature gas cooling and washing equipment and its industrial application: 1322795[P]. 2001-11-21.

[3]顏留成, 王亦飛, 王靈萍, 等. 洗滌冷卻環(huán)內(nèi)液體冷態(tài)流動(dòng)行為數(shù)值模擬 [J]. 化工學(xué)報(bào), 2014, 65(8): 2927-2933. DOI: 10.3969/j.issn. 0438-1157.2014.08.011. YAN L C, WANG Y F, WANG L P, et al. Numerical simulation of cold state liquid flow in scrubbing-cooling ring [J]. CIESC Journal, 2014, 65(8): 2927-2933. DOI: 10.3969/j.issn.0438-1157.2014.08.011.

[4]宋建, 胡珀, 韋勝杰, 等. 豎直平板水膜表面波演化行為實(shí)驗(yàn)研究[J]. 原子能科學(xué)技術(shù),2012,46(7)：804-808. SONG J, HU P, WEI S J, et al. Experimental study on surface-wave evolution behavior of water film falling down vertical plate [J]. Atomic Energy Science and Technology, 2012, 46(7): 804-808.

[5]王靈萍, 王亦飛, 郭強(qiáng)強(qiáng), 等. 洗滌冷卻管內(nèi)垂直降膜流動(dòng)特性[J]. 化工學(xué)報(bào), 2013, 64(6): 1959-1968. DOI:10.3969/j.issn.0438-1157. 2013.06.009. WANG L P, WANG Y F, GUO Q Q, et al. Flow characteristics of vertical falling film in scrubbing-cooling pipe [J]. CIESC Journal, 2013, 64(6): 1959-1968. DOI: 10.3969/j.issn.0438-1157.2013.06.009.

[6]閻維平, 李洪濤, 葉學(xué)民, 等. 垂直自由下降液膜厚度的瞬時(shí)無接觸測(cè)量研究 [J]. 熱能動(dòng)力工程, 2007, 22(4): 380-384. DOI:10.3969/j.issn.1001-2060.2007.04.007. YAN W P, LI H T, YE X M, et al. A study of the transient non-contact measurement of vertical free-falling liquid-film thickness [J]. Journal of Engineering for Thermal Energy and Power, 2007, 22(4): 380-384. DOI: 10.3969/j.issn.1001-2060.2007.04.007.

[7]韋勝杰, 宋建, 胡珀, 等. 豎壁冷態(tài)降液膜流動(dòng)統(tǒng)計(jì)特性實(shí)驗(yàn)研究[J]. 原子能科學(xué)技術(shù), 2012, 46(6): 674-678. WEI S J, SONG J, HU P, et al. Statistical characteristics of water film falling down large flat plate [J]. Atomic Energy Science and Technology, 2012, 46(6): 674-678.

[8]ABDULKADIR M, AZZI A, ZHAO D, et al. Liquid film thickness behaviour within a large diameter vertical 180° return bend [J]. Chemical Engineering Science, 2014, 107: 137-148.

[9]ALEKSEENKO S V, AKTERSHEV S P, CHERDANTSEV A V, et al. Primary instabilities of liquid film flow sheared by turbulent gas stream [J]. International Journal of Multiphase Flow, 2009, 35(7): 617-627.

[10]KOIZUMI Y, ENARI R, OHTAKE H. Correlations of wave characteristics for a liquid film falling down along a vertical wall [J]. Journal of Heat Transfer, 2009, 131(8): 082901 1-9.

[11]KIL S H, KIM S S, LEE S K. Wave characteristics of falling liquid film on a vertical circular tube [J]. International Journal of Refrigeration, 2001, 24: 500-509.

[12]MIN J K, PARK I S. Numerical study for laminar wavy motions of liquid film flow on vertical wall [J]. International Journal of Heat and Mass Transfer, 2011, 54(15/16): 3256-3266.

[13]李鐵, 吳晅,袁竹林.氣化爐激冷室下降管內(nèi)氣液兩相熱質(zhì)同時(shí)傳遞過程數(shù)值研究 [J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2008, 28(26): 35-39.DOI:10.3321/j.issn:0258-8013.2008.26.007. LI T, WU X, YUAN Z L. Numerical investigation of gas-liquid heat and mass concurrent transfer process in vertical pipe in gasifier quenching chamber [J]. Proceedings of the CSEE, 2008, 28(26):35-39. DOI:10.3321/j.issn:0258-8013.2008.26.007.

[14]FEDDAOUI M, MEFTAH H, MIR A. The numerical computation of the evaporative cooling of falling water film in turbulent mixed convection inside a vertical tube [J]. International Communications in Heat and Mass Transfer, 2006, 33(7): 917-927.

[15]WANG Y, GUO Q, FU B, et al. Numerical analysis of the flow characteristics and heat and mass transfer of falling-water films in an industrial-scale dip tube of a WSCC in an OMB gasifier [J]. Industrial & Engineering Chemistry Research, 2013, 52(26): 9295-9300.

[16]BRITO D, NATAF H-C, CARDIN P, et al. Ultrasonic Doppler velocimetry in liquid gallium [J]. Experiments in Fluids, 2001, 31(6): 653-663.

[17]ECKERT S, GERBETH G, MELNIKOV Vi. Velocity measurements at high temperatures by ultrasound Doppler velocimetry using an acoustic wave guide [J]. Experiments in Fluids, 2003, 35(5): 381-388.

[18]張歡, 王鐵峰, 王金福. 液固循環(huán)流化床中顆粒軸向速度的實(shí)驗(yàn)研究 [J]. 高?；瘜W(xué)工程學(xué)報(bào), 2002, 16(4): 2-8. DOI:10.3321/j.issn: 1003-9015.2002.04.010. ZHANG H, WANG T F, WANG J F. Study on the axial velocity of particles in a liquid-solid circulating fluidized bed using the ultrasound-Doppler technique [J]. Journal of Chemical Engineering of Chinese Universities, 2002, 16(4): 2-8. DOI:10.3321/j.issn: 1003-9015. 2002.04.010.

[19]LONGO S. The effects of air bubbles on ultrasound velocity measurements [J]. Experiments in Fluids, 2006, 41(4): 593-602.

[20]王鐵峰, 王金福, 張歡, 等. 超聲多普勒測(cè)速儀在液-固和氣-液兩相流測(cè)量中的應(yīng)用 [J]. 化工學(xué)報(bào), 2002, 53(4): 427-432. WANG T F, WANG J F, ZHANG H, et al. Application of ultrasound Doppler velocimetry in multiphase flow measurement [J]. Journal of Chemical Industry and Engineering (China), 2002, 53(4): 427-432.

[21]NOWAK M. Wall shear stress measurement in a turbulent pipe flow using ultrasound Doppler velocimetry [J]. Experiments in Fluids, 2002, 33(2): 249-255.

Characteristics of velocity and thickness distribution of liquid film in vertical falling tube

WANG Jing, WANG Yifei, YAN Liucheng, WU Ziwei, YU Guangsuo
(Key Laboratory of Coal Gasification and Energy Chemical Engineering of Ministry of Education, Institute of Clean Coal Technology, East China University of Science and Technology, Shanghai 200237, China)

Abstract:The flow characteristic of the falling liquid film in the scrubbing-cooling tube was researched in this paper. The ultrasound Doppler velocimetry was adopted to measure the velocity and thickness distribution of the liquid film non-instrusively at different circumferential and axial position. The Reynolds number ranged from 1.0×104to 3.1×104. The results showed that both the velocity and thickness of the liquid film were the largest at the 0° position, leading to unevener thickness distribution and the separation of the local liquid film and the sharp increase of the liquid film thickness at the 8° and 16° positions. In the axial direction, because of the effect of gravity, the liquid film velocity increased with the increase of flow distance when the Reynolds numbers was less than 2.0×104, whereas, it decreased at large Reynolds numbers since that the resistance became a dominant factor. Both of the thickness and velocity of global liquid film increased with increasing Reynolds number. Furthermore, the increase of the Reynolds numbers made the liquid film more unstable.

Key words:gas-liquid flow; turbulence; film; velocity measurement; thickness measurement; ultrasound Doppler velocimetry

中圖分類號(hào)：TQ 021

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：0438—1157（2016）06—2239—07

基金項(xiàng)目：國(guó)家高技術(shù)研究發(fā)展計(jì)劃項(xiàng)目(2012AA053101)。

Corresponding author:Prof. WANG Yifei, wangyf @ ecust. edu. cn

DOI：10.11949/j.issn.0438-1157.20151953