旋進旋渦流量計結(jié)構(gòu)參數(shù)優(yōu)化研究

借助計算流體力學(xué)（CFD）方法對旋進旋渦流量計進行了研究，獲得了流量計內(nèi)部流場信息，驗證了在流量計工作范圍內(nèi)旋渦進動頻率與流量之間良好的線性關(guān)系。重點對幾種旋進旋渦流量計結(jié)構(gòu)改進優(yōu)化方案進行了仿真分析，這些結(jié)構(gòu)優(yōu)化方案包括在旋進旋渦流量計起旋器入口加裝導(dǎo)流葉片、改變擴張段的擴張角和改變收縮比。結(jié)果發(fā)現(xiàn)，加裝導(dǎo)流葉片可以顯著減小流量計壓力損失，并可增加渦核的旋轉(zhuǎn)強度，使得檢測信號的強度得到增強，從而擴展了測量下限；收縮比和擴張角過大和過小對于流量計性能都不利，存在一個*值。

    旋進旋渦流量計是根據(jù)旋渦進動現(xiàn)象設(shè)計的一種流體振蕩式流量計，具有流量范圍寬、無可動部件、不易腐蝕、可靠性高、安裝使用方便、直管段要求短等優(yōu)點，適用于石油、蒸汽、天然氣、水等多種介質(zhì)的流量測量[1]。20世紀70年代，Dijsbergen[2]對旋進旋渦流量計性能進行了比較全面的實驗研究，驗證了該流量計線性輸出特性，并且發(fā)現(xiàn)該流量計不易受介質(zhì)黏度和密度影響，指出了旋進旋渦流量計在高壓氣體測量方面商業(yè)化應(yīng)用前景。Furio和Gianfranco[3]對旋進旋渦流量計做了實際工況下的儀表特征測試，探索該流量計在計量領(lǐng)域應(yīng)用的可行性。旋進旋渦流量計工作時振蕩頻率信號容易被外界振動與流體脈動噪聲所干擾[4]，對此宋開臣和傅新[5，6]提出安裝對稱兩個壓力信號探頭，相位差180°，用信號差分處理提高旋進旋渦流量計抗*力，取得了良好效果。

    對于旋進旋渦流量計內(nèi)部流動特性及流量計改進方面，科研人員也進行了一定探索。彭杰綱等人[7-9]對旋進旋渦流量計內(nèi)部流場進行了數(shù)值模擬分析，研究了旋渦進動效應(yīng)流場的演變情況，分析了流場干擾對旋進旋渦流量計流場進動效應(yīng)的影響。張濤等人[10]采用數(shù)值方針的方法對旋進旋渦流量計的結(jié)構(gòu)進行了優(yōu)化，使得流量計的壓力損失有了一定減小。何馨雨等[11]對旋進旋渦內(nèi)部流場進行了數(shù)值模擬分析，獲得了比較全面的流場信息，對這種流量計的內(nèi)部流動特性有了更加深入的理解。

    目前，旋進旋渦流量計應(yīng)用中還存在壓損較大、下限流量（始動流量）偏大、系列設(shè)計缺乏理論依據(jù)等問題。本文借助流體力學(xué)數(shù)值仿真的方法對旋進旋渦流量計的結(jié)構(gòu)參數(shù)做了系統(tǒng)研究，分析了旋進旋渦流量計壓力損失情況和輸出信號質(zhì)量。文中采用導(dǎo)流葉片來降低壓損，提高流量計性能，還著重考慮針對起旋器葉片數(shù)量、擴張角度等參數(shù)進行的優(yōu)化研究。希望通過此項研究工作能夠為旋進旋渦流量計的設(shè)計開發(fā)提供理論上的支持。

    1 旋進旋渦流量計工作原理與結(jié)構(gòu)優(yōu)化研究方案

    1.1 工作原理

    旋進旋渦流量計主要由起旋器、文丘里管、消旋器和檢測傳感器組成，其結(jié)構(gòu)原理如圖1。

    旋進旋渦流量計是基于旋渦進動現(xiàn)象工作的[12]。流體流入旋進旋渦流量計后，首先通過一組由固定螺旋形葉片組成的起旋器后被強制旋轉(zhuǎn)，使流體形成旋渦流。旋渦中心為“渦核”是流體旋轉(zhuǎn)運動速度很高的區(qū)域，其外圍是環(huán)流。流體流經(jīng)收縮段時旋渦加速，沿流動方向渦核直徑逐漸縮小，而強度逐漸加強。此時渦核與流量計的軸線相一致。當進入擴大段后，旋渦急劇減速，壓力上升，中心區(qū)域的壓力比周圍的壓力低，于是產(chǎn)生了局部回流。在回流作用下，渦核偏離中心軸像剛體一樣在擴張段壁面做螺旋進動，并且是圍繞中心軸進行的。進動頻率與流體的流速成正比。因此，測得旋進旋渦的頻率即能反映流速和體積流量的大小。

    1.2 結(jié)構(gòu)優(yōu)化方案

    針對現(xiàn)有旋進旋渦流量計壓損大、小流量信號弱等問題，提出了結(jié)構(gòu)改進和參數(shù)優(yōu)化研究方案。改進對象為一個DN50氣體旋進旋渦流量計，其內(nèi)部結(jié)構(gòu)和尺寸見圖2。流量計長度為232mm，入口管徑50mm，收縮段長度94.2mm，發(fā)展段管徑為36mm、長度為35.8mm，擴張段長為12mm，擴張角度為60°。具體改進和參數(shù)優(yōu)化研究如下：

    1）起旋器入口加裝導(dǎo)流葉片段。原有起旋器葉片在入口段沒有導(dǎo)流部分（圖1），葉片與來流之間夾角為60°，來流不是切向進入，會造成嚴重的流動分離。流動分離使得流動擾動和流動阻力增大。因此，改進方案中考慮把原來起旋器葉片延長并作彎曲（見圖3a），使得入口與來流夾角為0°，即流動切向進入，以期改善流動狀態(tài)。

    2）起旋器葉片數(shù)量增加。起旋器原來的葉片數(shù)量是6片，考慮增加1片（圖3b），同時為了保證流動面積不減小，將葉片厚度由2.5mm變?yōu)?.5mm（當前加工能力能夠做到）。螺旋角度保持不變，仍為30°。

    3）關(guān)于流量計殼體參數(shù)方面，對擴張角和收縮比進行研究。如前所述，旋渦進入擴張段后急劇減速，產(chǎn)生了局部回流，渦核偏離中心軸打在壁面上，產(chǎn)生壓力脈動。因此，擴張角大小對渦核運動會產(chǎn)生很大影響。本文將對30°、60°（原型參數(shù)）和90°（階梯型）三種情況進行模擬計算。收縮比也是比較重要的一個參數(shù)，收縮比越大，旋渦旋轉(zhuǎn)速度增加越多，但是同時壓損將變大，而且可能造成擴張段渦核擾動對沖擊壁面的影響減低。本文考察了50∶36（目前參數(shù)）和50∶32兩種收縮比的情況。需要說明的是，對擴張角和收縮比參數(shù)研究時，針對的起旋器模型是加裝了導(dǎo)流葉片（6個葉片）的情況。

    2 計算模型

    2.1 流體力學(xué)控制方程和湍流模型

    旋進旋渦流量計的流體動力特性可以用流體力學(xué)基本方程描述如下。

    連續(xù)性方程與動量方程[13，14]：

    式（1）（2）中：p—靜壓；ui—流動速度；f—質(zhì)量力，τij—應(yīng)力張量。

    流量計內(nèi)部為湍流流動，需引入湍流模型。標準的k2ε湍流模型用于強旋流或帶有彎曲壁面的流動時，會出現(xiàn)一定失真，因此本文選用RNGk2ε湍流模型。湍流模型和相關(guān)方程在文獻[11]中有詳細說明，感興趣的讀者可參考這篇文獻。

    2.2 仿真模型的建立與計算過程

    圖4為流量計仿真物理模型，計算區(qū)域總長度為447mm，網(wǎng)格劃分采用的是非結(jié)構(gòu)化混合網(wǎng)格，在局部流動復(fù)雜部分對網(wǎng)格進行加密。網(wǎng)格數(shù)大約為80萬。在旋進旋渦流量計的流通管道上設(shè)置五個壓力檢測截面1，2，3，4，5。截面1位于擴張段的中部，截面2位于擴張段上游5mm處，截面3位于擴張段上游10mm處，截面4擴張段上游15mm處，截面5位于擴張段下游5mm處。分別在每個截面距壁面6mm處設(shè)置對稱兩點作為壓力脈動信號采集點。

    本項工作采用了Fluent流體力學(xué)計算軟件。計算流程和條件參數(shù)設(shè)置如下：

    選用3d解算器；求解模型Solver設(shè)置為非耦合隱身算法，并設(shè)置為非定常流動。湍流模型是RNGk2ε湍流模型。設(shè)置流體物理屬性：選擇空氣密度1。225kg/m3。

    邊界條件：入口設(shè)為速度入口；出口邊界條件設(shè)定為壓力出口；壓力值固定為1個大氣壓（101325Pa）；管壁及起旋器表面設(shè)為無滑移邊界條件。

    流量從10m3/h開始，然后取20m3/h，40m3/h，再每間隔40m3/h取一次（80m3/h，120m3/h，160m3/h，200m3/h和240m3/h），計算中對非定常迭代時間步長的選取需注意要比脈動周期小得多，小于1/10脈動周期。

    3 數(shù)值模擬結(jié)果與分析

    何馨雨等計算發(fā)現(xiàn)在旋進旋渦流量計的流通管道截面1，2，3，4，5中，擴張段截面1及其上游5mm處（即截面2）處的旋渦脈動信號zui強。并且如果采用相位差為180°的兩個軸對稱點進行差分處理可以使脈動信號增強一倍，所以本研究選取截面1作為信號檢測面。

    首先考察了收縮比的影響。收縮比為50∶36時，在整個流量范圍內(nèi)，有正常數(shù)據(jù)獲得；而當收縮比為50∶32時（保持擴張角度60°不變），發(fā)現(xiàn)在流量為40m3/h以下和200m3/h以上時，壁面處無法檢測到壓力周期信號。也就是說，延長收縮段使得收縮比變大時，流量檢測范圍變小，達不到改善流量計測量特性的目的，因此，以下討論不再包含收縮比的內(nèi)容。

    3.1 流場特性和旋渦進動頻率

    圖5是針對原型流量計內(nèi)部流場的模擬結(jié)果（一個旋渦進動周期），流量為80m3/h。可以看到，在靠近壁面處速度明顯高于其他部位，而管道中心處由于回流的影響呈負速度分布，并且壁面高速流呈逆時針做圓周狀運動，即管道流體在擴張段做螺旋進動。

    表1列出了各種結(jié)構(gòu)參數(shù)情況下的旋渦進動頻率數(shù)據(jù)，這些數(shù)據(jù)繪制成旋渦進動頻率2流量關(guān)系圖（圖6）。圖中可以看出，對于研究涉及到的所有參數(shù)組合，旋渦進動頻率與流量之間都保持良好的線性關(guān)系。

    3.2 流量與壓損的關(guān)系

    在研究流量計壓損變化時，我們在計算模型的入口處與出口處分別設(shè)置監(jiān)測截面，入口壓力值減去出口壓力值即作為流量計壓力損失。不同流量，各種結(jié)構(gòu)參數(shù)情況下的壓損列于表2。

    圖7為壓力損失與流量的關(guān)系。圖中可以看出，原型流量計的起旋器前端加裝螺旋導(dǎo)流葉片后，壓損有了明顯的減小。流體進入流量計前為是充分發(fā)展管流，進入流量計后首先通過的是相當于一個阻流件的起旋器，流體受到起旋器前端葉片的阻擋，當未加導(dǎo)流葉片時，會發(fā)生嚴重的流動分離，從而帶來很大的能量損失；而前端增加螺旋導(dǎo)流葉片，流體以平行于流場的角度引導(dǎo)流體順利并呈螺旋狀流入起旋器起旋葉片，避免發(fā)生流動分離，壓力損失得以減小。當增加葉片數(shù)量至7片后，壓力損失比6個葉片時進一步減小，這是因為在增加葉片數(shù)量的同時減小了葉片厚度，流通面積面反而增大。從圖7還可以看出，將擴張段角度縮小至30°和增大至90°時，壓損曲線與60°（即原型）的情況相比基本沒有變化，說明了在流量計壓損方面，擴張角大小影響不大，起旋器的影響是主導(dǎo)因素。

    3.3 檢測探頭處的壓力信號

    對相位相差180°的兩路壓力信號做差分后能在一定程度上抑制管道同模干擾信號并且能使信號強度增大一倍，故而對小流量測量時更敏感，能有效提高測量精度。不同流量下原型流量計與參數(shù)改進后模型的差分壓力信號幅值列于表3。
圖7壓力損失與流量之間的關(guān)系

    圖8為各模型壓力信號幅值與流量的關(guān)系圖，可看出在起旋器前端加裝螺旋導(dǎo)流葉片后，壓力信號有了明顯的增強，平均增幅達到49%。葉片數(shù)增加至7片后壓力信號幅值比6片時有較大第2期周凱，等：旋進旋渦流量計結(jié)構(gòu)參數(shù)優(yōu)化研究111減弱，不可取。擴張角縮小至30°后信號較60°時變?nèi)?，擴張角增大至90°時信號zui弱，可見擴張角不是越大越好，也不是越小越好，就本文數(shù)據(jù)來說，60°是的一個選擇。

    4 結(jié)語

    通過數(shù)值計算的方法對旋進旋渦流量計內(nèi)部流場和輸出特性等進行了研究，針對幾種結(jié)構(gòu)改進方案進行了數(shù)值仿真，得到以下結(jié)果：

    1）對于各種結(jié)構(gòu)參數(shù)組合的模型，流量與測量信號頻率都能保持線性關(guān)系。

    2）在旋進旋渦流量計減小壓損方面，起旋器占主導(dǎo)因素?？s小擴張角至30°及增大擴張角至90°對改變壓損狀況不明顯，而通過前端增加螺旋導(dǎo)流葉片，流量計壓損有明顯降低。

    3）在測量信號強度方面，加裝導(dǎo)流葉片后*，達到近49%的增幅。特別在小流量時信號增強明顯，能有效降低測量下限。擴張角選取有*值，本研究中，擴張角為60°時測量信號zui強。