摘要:應用粒子成像測速技術(shù)獲得了柴油計量表葉片入口流場(chǎng)的速度分布信息，并基于該測量結果，運用T-G模型理論得出流量計的響應。通過(guò)與以往所采用的幾種典型的入口速度分布計算得到的結果比較分析表明，基于PIV測量的結果更接近于柴油計量表的真實(shí)響應。還比較分析了渦輪入口速度分布對柴油計量表響應的影響機理，相關(guān)結果可望為改進(jìn)柴油計量表響應的計算分析方法以及優(yōu)化設計提供有價(jià)值的參考。
1引言
柴油計量表作為一種速度式儀表，因其具有諸多優(yōu)點(diǎn)被廣泛應用到工業(yè)生產(chǎn)以及實(shí)際生活中。在柴油計量表的實(shí)際使用中，一般需要結合流量計本身的響應曲線(xiàn)來(lái)計算被測管流的實(shí)際流量。獲得柴油計量表響應曲線(xiàn)的方法主要有2種，一是通過(guò)標準流量平臺標定，二是通過(guò)理論模型計算獲得其響應。其中標定方法在實(shí)際生產(chǎn)中應用更廣泛，不過(guò)特定的標定曲線(xiàn)往往僅適用于某些單一工況下的響應，局限性較大。因此通過(guò)對柴油計量表理論模型的探索和改進(jìn)從而更準確地預測流量計的響應曲線(xiàn)具有重要意義。
1970年，Thompson和Grey基于葉柵理論和不可壓縮勢流提出了較為系統的計算柴油計量表響應的理論模型(以下簡(jiǎn)稱(chēng)TG模型)。該模型可以將速度入口信息以及柴油計量表各部件的幾何和運動(dòng)參數均納入考慮，因而被此后的研究廣泛采用。流量測量設備的內流場(chǎng)對其響應有著(zhù)重要影響，目前部分研究采用數值模擬手段對其進(jìn)行計算進(jìn)而分析儀表的響應情況。對于渦輪內流場(chǎng)的實(shí)際流動(dòng)情況Xu采用了激光多普勒(LDA)技術(shù)對口徑100mm的柴油計量表輪轂與管壁間的12個(gè)不同位置的速度進(jìn)行了測量并代入理論模型進(jìn)行計算，理論計算結果與實(shí)驗結果比較吻合。
以上研究都表明，獲得準確的柴油計量表入口速度分布，結合TG模型可大大提高計算的準確性。由于實(shí)際渦輪上游情況比較復雜，不同的導流葉片、輪轂前緣設計等因素都對入口速度分布有著(zhù)重要的影響，因而實(shí)際的渦輪入口速度分布，并非均勻分布或充分發(fā)展的環(huán)空分布，難以通過(guò)簡(jiǎn)單的黏性流理論獲得通用的速度分布計算方法。因而采用實(shí)驗的手段，獲得能夠反映切合實(shí)際的渦輪入口速度分布具有重要的應用價(jià)值?？v觀(guān)以往的研究，對于入口速度分布的獲得，多是采用理論計算或是數值模擬的方法，僅有Xu采用了LDA技術(shù)觀(guān)察了渦輪環(huán)空的速度分布。LDA技術(shù)雖然計量精確，但其多光束匯聚和單點(diǎn)測量的性質(zhì)決定了它測點(diǎn)較少，難以同時(shí)獲得全流場(chǎng)信息的缺陷，因而其僅能用于口徑較大的柴油計量表流場(chǎng)測量。隨著(zhù)粒子圖像測速(PIV)技術(shù)的發(fā)展，其瞬時(shí)獲得全場(chǎng)信息的能力也被應用到流量測量的研究中，也可以用于柴油計量表葉輪輪轂與管壁的研究中?；谝陨峡紤]，本文應用粒子成像測速技術(shù)(PIV)來(lái)獲得流場(chǎng)的流速信息，以便通過(guò)更準確和全面的入口速度分布進(jìn)而對流量計的響應獲得認識上的深入。
2柴油計量表理論模型
在柴油計量表處于穩定響應的狀態(tài)下，角加速度為零，此時(shí)作用在葉輪上的各力矩(見(jiàn)圖1)需滿(mǎn)足力矩平衡方程:
Td－Th－Tt－Tw－Tb－Tm=0(1)
式中:Td為葉片驅動(dòng)力矩;Th為輪轂周邊黏性阻力矩;Tt為葉片頂隙黏性阻力矩;Tw為輪轂端面黏性阻力矩;Tb為軸承黏性阻力矩;Tm為軸與軸尖機械阻力矩

渦輪轉速ω為所求的量。獲得各不同力矩，通過(guò)求解力矩平衡方程可得對應工況的渦輪轉速ω。
下面分別對各力矩計算方法進(jìn)行介紹。
2．1葉片驅動(dòng)力矩
本文采用的理論模型葉片驅動(dòng)力矩類(lèi)似Xu文章中的處理方法，模型假設在葉輪輪轂和管道內壁之間的環(huán)空空間內不存在沿半徑方向的流動(dòng)，因而可以將三維的渦輪葉片離散成有限個(gè)二維葉柵計算不同葉柵上葉片的受力。針對半徑為r處的葉柵，可計算其所受的驅動(dòng)力系數Cdr(見(jiàn)圖1)：
Cdr=Clcosβ－Cdsinβ(2)
式中:Cl和Cd為葉柵的升力系數和阻力系數，二者可通過(guò)不可壓縮勢流的方法計算，與葉片翼型、葉輪轉速、來(lái)流速度和半徑位置等參數相關(guān)。
通過(guò)對輪轂半徑Ｒh到葉頂半徑Ｒt的驅動(dòng)力矩進(jìn)行積分可以獲得作用在整個(gè)葉片上的驅動(dòng)力矩Tdr(不含黏性影響):

式中:ρ為流體密度;N為葉片數量;C為葉片弦長(cháng);Ur(r)為葉柵處速度，與葉柵位置相關(guān)，在本文中通過(guò)PIV測量結果插值獲得。
實(shí)際流動(dòng)中受黏性影響，流體還會(huì )在葉柵表面產(chǎn)生黏性力，采用黏性流體力學(xué)中二維渠道流平板黏性力的計算方法計算黏性阻力Fv

式中:t為葉柵柵距，ν為流體的運動(dòng)黏性系數。
由葉柵黏性阻力Fv可求得葉片整體所受黏性阻力矩Tv:

進(jìn)而可以獲得葉片上所受的整體驅動(dòng)力矩Td:
Td=Tdr－Tv(7)
2．2輪轂周面黏性阻力矩
本文理論模型中輪轂周面黏性阻力矩由2部分構成:葉片部分輪轂和葉片上游輪轂。
依據Tsukamoto的計算，葉片部分輪轂黏性阻力矩Thb的計算式為：

式中:Bt為葉片厚度。
2．4輪轂端面黏性阻力矩
依據Tsukamoto的計算，輪轂端面黏性阻力矩Tw的計算式為:

2．5軸承黏性阻力矩和機械摩擦阻力矩
根據同軸圓筒黏性阻力矩計算方法可得軸承黏性阻力矩Tb

式中:Ｒb和Ｒbo分別為軸和軸承半徑，lb為軸的等效長(cháng)度。
機械摩擦阻力矩基本不受轉速影響可設置為定值，本文中渦輪機械摩擦阻力矩取為5×10－7N·m。
2．6理論模型綜合分析
當渦輪進(jìn)入線(xiàn)性響應區間后，起主要作用的是葉片驅動(dòng)力矩和葉片頂隙阻力矩之間的平衡，其他各阻力矩相對較小。葉片頂隙阻力矩與渦輪軌速矩近似成正比關(guān)系，驅動(dòng)力矩則主要受入口速度分布Ur(r)影響，獲得準確的入口速度分布可以使理論模型的計算結果與實(shí)際更為符合，傳統的理論模型中入口速度分布多采用均勻分布假設(即各不同半徑入口速度相等)或充分發(fā)展的環(huán)空空間速度分布，本文則通過(guò)PIV技術(shù)測量了實(shí)驗使用渦輪的入口速度分布并代入理論模型進(jìn)行計算。
3實(shí)驗系統
實(shí)驗在中國科大水平循環(huán)水實(shí)驗平臺上進(jìn)行，實(shí)驗采用20mm口徑的柴油計量表，量程范圍是1～80方/天，其中較好線(xiàn)性段范圍是5～50方/天，流量計渦輪為等重疊度渦輪(不同半徑位置葉柵重疊度相同)，具體參數如表1所示，這種流量計在大慶油田的生產(chǎn)測井中廣泛應用，其結構如圖2(a)所示。來(lái)流經(jīng)過(guò)一段導流葉片整流后進(jìn)入渦輪的環(huán)空空間，驅動(dòng)葉輪轉動(dòng)，輸出響應信號。在柴油計量表的線(xiàn)性響應區間中，處于穩定轉動(dòng)時(shí)渦輪葉片對流體的干擾較少，相對來(lái)流攻角較小，對流體的軸向速度分布基本沒(méi)有影響，僅會(huì )稍稍增加其周向轉速。因而為了測量柴油計量表入口速度分布，特別制作了各參數與實(shí)際渦輪相同但并無(wú)葉片的透明外殼輪轂模型，如圖2(b)所示，通過(guò)PIV手段，對管道中軸面上輪轂和管壁之間的區域的軸向速度分布進(jìn)行剖面測量。輪轂模型安裝在待測柴油計量表的上游，相距超過(guò)2m以保證二者之間無(wú)相互干擾。實(shí)驗流速范圍在5～25方/天，在管路下游采用時(shí)間-質(zhì)量法獲得真實(shí)流速，通過(guò)光學(xué)觀(guān)測獲得柴油計量表葉輪的真實(shí)轉動(dòng)頻率，同時(shí)采用PIV技術(shù)測量輪轂模型中的速度分布。


本文實(shí)驗中所使用的PIV系統為作者單位自行研制的PIV系統(見(jiàn)圖3(a))，激光器發(fā)出的激光依次通過(guò)凸透鏡聚焦，經(jīng)柱面鏡發(fā)散成片光，再通過(guò)平面反射鏡反射成豎直片光，進(jìn)入實(shí)驗觀(guān)察區。示蹤粒子跟隨流體流過(guò)實(shí)驗段，由高速攝影記錄實(shí)驗過(guò)程，通過(guò)相關(guān)計算處理得到速度分布結果。其中所用的激光器為可連續發(fā)射532mm激光(綠光)，發(fā)射*大輸出功率為2W的半導體激光器。實(shí)驗采用的相機為每秒可拍攝5000幅的高速攝影。示蹤粒子采用的是空心玻璃微球，粒徑為20～40μm，密度1．05g/cm3。圖像的互相關(guān)處理程序由作者所在單位自行在MATLAB軟件平臺中編寫(xiě)成。
本文進(jìn)行圖像采集的方法均為多幀單曝光，即相繼2次曝光的粒子圖像分別記錄在相繼的2幅照片上，因此采取互相關(guān)算法進(jìn)行圖像處理?；驹硎怯孟嗬^2幀粒子圖像I1(x珋)，I2(x珋)進(jìn)行相關(guān)計算：

Ｒc(sˉ)的理想空間分布如圖3(b)所示，僅有一個(gè)明顯的級大峰值，其中珋s為判讀小區內粒子的平均位移矢量。本文算法采用16×16的矩形像素作為判讀小區，對應的空間分辨率為0．35mm×0．35mm，時(shí)間分辨率達0．2ms。整個(gè)圖像在輪轂與管壁之間的速度剖面可取的32個(gè)流速點(diǎn)，從而可以較準確地得到其間的速度分布情況。

4實(shí)驗結果與分析
通過(guò)對PIV實(shí)驗中所拍攝的照片(見(jiàn)圖4(a))進(jìn)行后處理，可以得到各個(gè)流量點(diǎn)下輪轂與管壁之間軸截面流場(chǎng)軸向速度分布信息如圖4(b)所示。鑒于試驗模型的軸對稱(chēng)性，從原理上說(shuō)該軸截面的速度分布可以推廣到周向環(huán)形區域。


實(shí)驗中流量計的渦輪輪轂半徑為4mm，而管道內徑為10mm，因而速度分布都在這6mm的區間內。通過(guò)圖像處理可以獲得32個(gè)不同位置的速度，在此基礎上進(jìn)行插值即可獲得整個(gè)環(huán)空流場(chǎng)的軸向速度分布。圖5(a)反映了實(shí)驗所測得的幾個(gè)不同工況點(diǎn)的軸向速度分布，從圖中可以看出，流速在中間位置較高，由于邊界層的影響，在靠近輪轂和管壁附近流速逐漸趨近于零。隨著(zhù)流速的升高，整體速度分布向管壁方向偏移，速度*大值位置半徑增大，輪轂表面邊界層厚度增加，管壁表面邊界層厚度減少。與Xu采用LDA測量的結果相比，本文結果在半徑較大處速度較高，二者的不同結果也反映了不同設計的柴油計量表入口速度分布存在差異。相比LDA而言，PIV可以更加全面地獲得輪轂與管壁之間的流速分布信息。


將實(shí)驗中PIV測得的速度分布與同流量下的完全發(fā)展的環(huán)形通道速度分布以及此流量下的均勻分布進(jìn)行對比，如圖5(b)所示，從中可以看出，用PIV測得的速度分布與完全發(fā)展的環(huán)形通道速度分布有明顯不同。其中前者的峰值比較靠近管道內壁方向，而后者的峰值較靠近輪轂方向。另外，完全發(fā)展的環(huán)形通道速度分布比用PIV測得的速度分布更加平緩。由于不同位置的流體對渦輪葉片作用效果不同，實(shí)際流速中峰值在不同位置對渦輪產(chǎn)生的驅動(dòng)效果可能會(huì )有很大差異，如圖5(b)中所示的完全發(fā)展速度分布和均勻速度分布都很平緩，不能完全反映實(shí)際流動(dòng)中不同位置的流場(chǎng)信息，計算的結果中自然也就將這些差異對渦輪響應可能產(chǎn)生的特殊貢獻有所體現。
分別用3種速度分布作為渦輪入口速度分布求解流量計響應，與實(shí)際測得的響應進(jìn)行對比，如圖6(a)所示。從圖中可以看出，采用完全發(fā)展的環(huán)形速度分布和均勻速度分布計算的渦輪響應值明顯低于柴油計量表的真實(shí)響應。在渦輪正常響應時(shí)，在葉片中上部(即靠近管壁部分)的流體驅動(dòng)渦輪轉動(dòng)，而在葉片底部(即靠近輪轂部分)的流體阻礙渦輪轉動(dòng)，因而流體分布越靠近管壁，帶來(lái)的驅動(dòng)力矩越大，使得葉輪的轉速越快。從圖5(b)中來(lái)看，真實(shí)速度分布更靠進(jìn)管壁，應為理論求解結果轉速偏低的原因。
分別計算各個(gè)響應在不同流量點(diǎn)處與真實(shí)響應的相對誤差，結果如圖6(b)所示。從圖中可以明顯看出，用PIV獲得的速度分布計算的結果與實(shí)際響應的相對誤差*小，在3%以?xún)?用均勻入口速度分布計算的結果誤差*大;用完全發(fā)展的環(huán)形通道速度分布當雷諾數超過(guò)2000時(shí)，由于計算模型假設由層流的速度分布直接轉變?yōu)橥牧鞯乃俣确植寄Ｐ?，未能合理地反映?shí)際流動(dòng)中逐步轉變的過(guò)渡階段，導致理論結果與實(shí)際速度分布有較大差異，所以誤差較大。通過(guò)這些對比不難看出，獲得真實(shí)的速度分布能更為準確地計算柴油計量表的實(shí)際響應。面對復雜的上游來(lái)流條件，PIV結果更能反映柴油計量表內部流動(dòng)的主要特征，這也從另一個(gè)側面表明，對柴油計量表內部復雜流動(dòng)的精細測量和深入認識也將是完善相關(guān)理論和進(jìn)一步優(yōu)化其性能的重要途徑。

5結論
本文通過(guò)PIV技術(shù)觀(guān)測了柴油計量表入口軸向速度分布并代入TG模型進(jìn)行計算。結果表明，PIV技術(shù)可以作為柴油計量表的入口速度的觀(guān)測手段。PIV技術(shù)代入模型后計算所得的渦輪轉速與實(shí)際較為吻合，而采用均勻速度入口或是充分發(fā)展的環(huán)空空間速度分布均與實(shí)際存在些許差異，代入模型后所得誤差較大，也反映了不同的入口速度分布對流量計響應具有十分重要的影響。
由于柴油計量表入口速度分布受到多種因素的影響，難以完全依賴(lài)簡(jiǎn)單的理論計算，因而PIV技術(shù)可以有的放矢地用于柴油計量表的內流場(chǎng)觀(guān)察，獲得真實(shí)的流速分布信息，進(jìn)而改進(jìn)理論模型的計算和分析，在新一代柴油計量表的研制和完善相關(guān)理論中發(fā)揮重要的作用。