氣體渦輪流量計旋轉部件內流場模擬 發布時間:2018-02-01
0引言 氣體渦輪流量計是計量天然氣、液化氣、煤氣等介質的速度式儀表[1-2]。為了改善氣體渦輪流量計的性能,為設計提供指導和方向,近年來一些學者利用CFD技術對其內部流場進行了研究。LavanteEV等[3]利用FLUENT對氣體渦輪流量計內部流場進行數值模擬,并根據仿真結果解釋實驗過程中的現象。對前導流器引起的流量計壓力損失進行數值計算和實驗測量,從流動機理上解釋了結構和壓損之間的關系。LIZhifei等[6]利用數值模擬得到了導流器內部的速度場和壓力場,并以減小壓力損失為目標優化了導流器的結構。通過對氣體渦輪流量計進行CFD仿真,研究不同流量下的壓損值,并通過實驗證明了數值模擬的有效性。對渦輪傳感器內部的速度場和壓力場進行了數值仿真,提出對前后導流器、葉輪葉片形狀和頁頂間隙的改進。 上述研究中未涉及針對不同螺旋升角渦輪內流場的數值模擬,以及渦輪葉片螺旋升角的改變對儀表性能影響的研究。本文對安裝35°和45°葉片螺旋升角渦輪的DN150型氣體渦輪流量計的內流場進行數值模擬,通過模擬結果預測儀表的始動流量和壓力損失,并利用預測的正確性,為渦輪葉片螺旋升角的進一步提供數值方法。 1數學模型及邊界條件 利用FLUENT軟件對渦輪內流場進行數值模擬時,忽略天然氣的密度變化,在0~1200m3/h內,介質流動速度遠遠小于聲速(即馬赫數遠小于0.3),認為流體不可壓縮,且假設流動中無熱量交換,不考慮能量守恒方程。 1.1微分控制方程 氣體渦輪流量計內部流動為湍流黏性流動,滿足連續性方程和黏性流體運動方程。 基本微分方程[9]: 連續性方程: 1.2 湍流模型選擇由于雷諾應力項的加入使時均N-S方程不封閉,為了求解引入k-ε兩方程湍流模型。兩方程湍流模型有標準k-ε模型,Renormalization-group(RNG)k-ε模型,和可實現的k-ε模型。其中,RNGk-ε模型主要應用于旋轉機械的流動問題,在大范圍的湍流模擬中有較高的精度。該模型能夠比較準確地模擬各種復雜流動,其中湍流黏度由下式確定: 1.3網格劃分與定解條件 根據流量計的實際工況分別在介質入口和出口處添加10倍管徑的直管段,并把整個模型剖分為3個區域:入口管道,旋轉區,出口管道。旋轉區域又細分為渦輪轉子和支架定子兩個區域,定子和轉子之間的耦合采用多參考MRF(MultipleReferenceFrame)模型。利用GAMBIT前處理模塊對進、出口直管段采用結構化網格,而對旋轉區采用非結構化網格進行劃分以滿足對葉輪內部復雜區域的網格描述,各塊網格通過塊之間的交界面拼接在一起。網格總數為30多萬個四面體非結構化網格和100多萬個六面體結構化網格,旋轉區網格如圖1所示。 定解條件包括介質入口、出口和固壁邊界的設置。入口處給定相應流量(1200m3/h)下的主流速度值;出口采用壓力出口邊界條件,出口壓力相對大氣壓為0;進、出口管道內壁,支架均取無滑移固壁邊界條件。葉輪部分采用旋轉坐標系,給定相應流量下的葉輪轉速,將葉片的吸力面和壓力面以及輪轂定義為旋轉壁面條件,在旋轉壁面條件的定義中,按照MRF的要求,將旋轉壁面的旋轉速度定義為相對速度,并且相對周圍流體速度為0。 2計算結果分析 2.1壓力場分析 流量計全壓定義為入口全壓與出口全壓之差,通過全壓分析能夠直接反映儀表壓損的大小。全壓越大表明流體經過流量計后產生的壓損越大,壓損過大會導致流量計不能正常使用。進口全壓一定時,出口全壓Pout越大,則流量計的全壓△P越小,壓力損失越小。如圖2(a)和圖3(a)所示,35°渦輪出口全壓要明顯小于45°渦輪出口處的全壓,這說明相同的工況下45°渦輪所產生的壓損較小。 渦輪葉片動壓的分布和大小直接影響渦輪驅動力矩的大小,35°葉片所受動壓明顯小于45°葉片所受動壓,說明在相同工況下45°螺旋升角渦輪能獲得較大的驅動力矩,如圖2(b)和圖3(b)所示,與35°渦輪相比,較小的流量就可推動渦輪穩定旋轉,從而使儀表進入線性工作區。由此可預測安裝45°螺旋角渦輪的儀表能獲得較小的始動流量。 2.2速度場分析 當氣體介質以充分發展的湍流經過渦輪時,35°渦輪的速度矢量方向變化較大且向壁面集中,使得與葉片直接作用產生推動力矩的速度矢量減少,如圖(4a)所示,且在出口處速度衰減較大,間接說明介質流經渦輪后壓損的增加,如圖4(b)所示。而45°渦輪內部的速度矢量分布比較均勻,過流性較好,與葉片直接作用的速度矢量較多,產生較大的驅動力矩,如圖(5a)所示,且在出口處速度衰減較小,如圖(5b)所示。 3實驗對比 氣體渦輪流量計的檢定采用負壓檢測方法,如圖6所示,由標準吸風裝置產生負壓使標準羅茨流量計和被檢定的氣體渦輪流量計同時測量,安裝在被測儀表兩端取壓口處的U型管可以測量流量計進、出口處的壓力,從而得到儀表的壓力損失。 利用黃金分割法選取0~1200m3/h范圍8個流量點,在每一個流量點隨機采集3組不同時刻的數據,包括標準羅茨流量計和被檢定流量計的累積流量及其輸出脈沖數,對每組數據進行算術平均得到流量點處的平均儀表系數。通過采集U型管壓差裝置的指示值記錄每個流量點處的壓力損失,檢定結果如表2所示。 利用多項式插值對表2中的數據進行密化,得到20組插值數據,通過3次B樣條擬合得到儀表系數曲線和壓力損失曲線。 始動流量以儀表系數進入線性區的最小流量來確定,在小流量區內安裝45°螺旋升角渦輪的流量計在流量20m3/h左右即進入線性工作區;而安裝35°螺旋升角渦輪的流量計則在流量150m3/h左右時才進入線性工作區,而且在線性工作區內也存在著明顯的波動,如圖7所示. 35°渦輪流量計在各工況點處的壓損明顯大于45°渦輪流量計,最大壓損達到3500Pa以上,如圖8所示。上述分析表明安裝45°螺旋升角渦輪的流量計與安裝35°螺旋升角渦輪的流量計相比具有較小的始動流量,較小的壓力損失,而且儀表計量的線性度較好。 4結論 對螺旋升角為35°和45°的氣體渦輪流量計旋轉部件內流場進行數值模擬,分析描述其內部流動的壓力場和速度場,安裝45°螺旋升角渦輪的流量計比安裝35°螺旋升角渦輪的流量計具有較小的始動流量和壓力損失。 利用黃金分割法選取儀表流量范圍內的檢定點,通過儀表負壓檢定平臺獲得了儀表系數曲線和壓力損失曲線,與數值仿真中的預測相吻合,表明數值模在流量計性能預測中的有效性。 渦輪葉片的螺旋升角是影響儀表性能的關鍵參數,合理選擇渦輪的葉片螺旋升角,可進一步改善儀表的性能。
以上內容來源于網絡,如有侵權請聯系即刪除!