摘要:超臨界二氧化碳循環(huán)發(fā)電技術(shù)采用超臨界狀態(tài)下的二氧化碳代替?zhèn)鹘y(tǒng)水蒸氣工質(zhì),由于物性發(fā)生巨大變化,質(zhì)量流量測量與計算方法需重新構(gòu)建。本文開展了孔板流量計數(shù)值模擬研究,使用fluent軟件模擬計算了孔板直徑比、節(jié)流孔厚度、孔板厚度等結(jié)構(gòu)參數(shù)對流量系數(shù)的影響,結(jié)果表明現(xiàn)行標(biāo)準(zhǔn)對超臨界二氧化碳并不適用。本文提出了針對超臨界二氧化碳工質(zhì)的孔板流量計結(jié)構(gòu)參數(shù)推薦范圍,在該范圍內(nèi)絕大多數(shù)計算結(jié)果相對誤差小于2%,并針對入口邊緣鈍化提出了新修正系數(shù),修正后計算結(jié)果相對誤差為0.11%~1.85%,滿足測量精度要求。
伴隨著經(jīng)濟社會的不斷發(fā)展,我國發(fā)電機組在發(fā)電效率、能源結(jié)構(gòu)、環(huán)境效益等多方面都面臨著轉(zhuǎn)型升級的嚴(yán)峻挑戰(zhàn)。中國在2007年已經(jīng)成為全球溫室氣體第一大排放國家川。習(xí)近平總書記在聯(lián)合國-般性辯論會.上承諾[2],中國碳排放量力爭分別于2030年和2060年前實現(xiàn)達峰和中和,作為一種變革性火力發(fā)電技術(shù)超臨界二氧化碳循環(huán)是實現(xiàn)能源結(jié)構(gòu)轉(zhuǎn)型進程中的重要手段。
隨著“雙碳”目標(biāo)的逐步落實，超臨界二氧化碳(溫度高于303.98K、壓力高于7.38MPa)被廣泛應(yīng)用,其作為工質(zhì)的布雷頓循環(huán)具有極高的熱源適用性,可應(yīng)用于太陽能、核能、余熱等多種場景。且由于超臨界二氧化碳密度大、黏性小、壓縮性好、循環(huán)過程無相變，相比于傳統(tǒng)水工質(zhì),超臨界二氧化碳循環(huán)珂以實現(xiàn)更高的循環(huán)效率,Dostal等[3]指出在透平人口工質(zhì)溫度高于550℃條件下,超臨界二氧化碳循環(huán)發(fā)電系統(tǒng)性能顯著高于水循環(huán)系統(tǒng)。此外，配合間歇性、隨機性強的可再生能源供電以保障社會用電穩(wěn)定是未來火力發(fā)電重要任務(wù),超臨界二氧化碳系統(tǒng)靈活性高、能實現(xiàn)完全熱電解耦的特點也使其更能滿足未來火力發(fā)電的深度調(diào)峰需求。
系統(tǒng)內(nèi)部流動工質(zhì)流量的正確測量是其得以應(yīng)用的基礎(chǔ)。當(dāng)前超臨界二氧化碳主要用于藥物/化學(xué)試劑萃取、油田驅(qū)油等溫度、壓力相對較低的場景,而超臨界二氧化碳循環(huán)系統(tǒng)需要二氧化碳工質(zhì)達到極高的溫度與壓力,二氧化碳的密度、比熱、黏度等物性參數(shù)發(fā)生了顯著變化,對于該條件下二氧化碳流量測量,傳統(tǒng)流量測量方法將不再適用。孔板流量計是--種技術(shù)成熟且適合于高溫高壓流體流量測:量的方法,經(jīng)過多年發(fā)展孔板流量計已形成標(biāo)準(zhǔn)化形式,主要包括兩部分,分別是具有直角邊緣的-段節(jié)流孔,以及在節(jié)流孔后具有一斜角的錐形擴流段,其結(jié)構(gòu)如圖1所示。然而有關(guān)測量超臨界二氧化碳循環(huán)中工質(zhì)流量的孔板流量計設(shè)計方案,國內(nèi)外并無經(jīng)驗借鑒。因此需要針對超臨界二氧化碳工質(zhì)的全新特性,探究孔板結(jié)構(gòu)參數(shù)的變化對于流量系數(shù)的影響,同時驗證現(xiàn)有標(biāo)準(zhǔn)中的相關(guān)規(guī)定對于超臨界二氧化碳工質(zhì)是否適用。.
我國學(xué)者采用數(shù)值模擬為主,實驗驗證為輔的研究方式,以水或天然氣為研究對象,開展了管徑、孔徑厚度等結(jié)構(gòu)參數(shù)對孔板流量計的影響研究�？装逯睆奖�、厚度等參數(shù)會顯著影響孔板的節(jié)流特性,從而影響流量計的計量性能。當(dāng)直徑比小于0.3時,流量系數(shù)隨直徑比增加而快速下降,當(dāng)直徑比大于0.3時，流量系數(shù)逐漸遞增,但增速較緩;直徑比在0.2~0.8范圍內(nèi)時，流量系數(shù)隨β增大呈先減小后增大的趨勢,并以0.55為分界點,其中β在0.45~0.65之間時可控制誤差在3%以內(nèi)。與直徑比不同,流量系數(shù)隨孔板厚度的變化特性較--致。厚度e增加，流出系數(shù)直線上升;林棋等人[4-5]也認(rèn)為流出系數(shù)隨縮徑孔厚度增大而增大;在模型中考慮了引壓管的存在,結(jié)果顯示,e變化0.15mm時，流出系數(shù)變化1.56%;e變化1mm時，流出系數(shù)變化2.125%。
近年來的理論知識、不斷優(yōu)化的算法以及不斷更新擴充的實驗數(shù)據(jù)庫等都保證了數(shù)值模擬研究的正確率與精度,因而逐漸成為主流研究方法之一�？装辶髁坑嫻艿纼�(nèi)部介質(zhì)流動復(fù)雜,參數(shù)變化劇烈,采用數(shù)值模擬方法可以有效捕捉到管道內(nèi)部的細微變化,因此是孔板流量計研究的有力工具。部分學(xué)者利用數(shù)值模擬對孔板流量計結(jié)構(gòu)進行了優(yōu)化設(shè)計。利用Fluent模擬了一種半雙曲線型的新式孔板流量計,并同時利用牛頓流體和非牛頓流體進行驗證，發(fā)現(xiàn)這種流量計可使內(nèi)部介質(zhì)近似無剪切流動,大大消除了渦流和停滯區(qū)等流動結(jié)構(gòu);研究發(fā)現(xiàn)在孔板流量計下游插入-個環(huán)可以有效減少壓力損失,并利用數(shù)值模擬和遺傳算法優(yōu)化結(jié)構(gòu),可減少33.5%的壓力損失,極大的降低了能耗和成本。
因此,本文進行了孔板流量計結(jié)構(gòu)參數(shù)對于流量系數(shù)影響的模擬研究,包括直徑比、節(jié)流孔厚度、孔板厚度等結(jié)構(gòu)參數(shù),明確了在超臨界二氧化碳工質(zhì)典型工況下不同結(jié)構(gòu)參數(shù)對流量系數(shù)的影響,同時將通過現(xiàn)行孔板流量計國際標(biāo)準(zhǔn)文件中經(jīng)驗公式計算得到的結(jié)果與數(shù)值模擬結(jié)果進行比較,提出了針對超臨界二氧化碳工質(zhì)的孔板流量計結(jié)構(gòu)參數(shù)推薦范圍與推薦設(shè)計值,提升了其測量精度。除此之外，還探究了孔板人口直角邊緣鈍化對孔板流量計測量精度的影響,并據(jù)此提出了新的針對現(xiàn)行孔板流量計國際標(biāo)準(zhǔn)文件中經(jīng)驗公式計算得到的流量系數(shù)的修正系數(shù)。
1.計算模型與模擬方法
1.1模型建立與網(wǎng)格劃分
根據(jù)標(biāo)準(zhǔn)文件[1]規(guī)定的孔板流量計結(jié)構(gòu)設(shè)計與參數(shù)要求，本文分別建立了DN25和DN200兩種管徑的孔板流量計,結(jié)構(gòu)參數(shù)如表1所示,在后文進行相關(guān)研究時均以該表中的結(jié)構(gòu)參數(shù)為基礎(chǔ)參數(shù),依據(jù)該參數(shù)使用SolidWorks軟件建立了孔板流量計及其前后一定長度管道的幾何模型,如圖2所示。
本文采用非結(jié)構(gòu)化網(wǎng)格進行模擬計算,利用AnsysMeshing軟件將孔板流量計管道劃分為四面體網(wǎng)格和六面體網(wǎng)格相結(jié)合的混合形式。除此之外,為了準(zhǔn)確捕捉到流場內(nèi)的細微變化,在介質(zhì)與管道內(nèi)壁接觸處進行邊界層的網(wǎng)格劃分，采用平滑過渡法,第一層高度根據(jù)面網(wǎng)格和過渡比進行確定,最大層數(shù)為5層，增長率為1.2,這時邊界層總厚度是變化的,對于復(fù)雜流動更有效，結(jié)果如圖3所示。
為提高節(jié)流孔板內(nèi)部及其到前后取壓截面處的模擬精度,利用影響球?qū)装迩昂箝L度為D的流場范圍內(nèi)進行了局部網(wǎng)格加密,網(wǎng)格數(shù)量過少會導(dǎo)致計算精度不足,而過多的網(wǎng)格數(shù)量則會無謂地加大計算工作量,降低計算速度。本文對DN25和DN200兩種管徑不同的管道進行網(wǎng)格數(shù)量與計算結(jié)果無關(guān)性的驗證，綜合計算精度與計算速度考慮,對于DN25管道,選取網(wǎng)格數(shù)分別為956036和1190483時,在各點測出的壓力相差均小于0.01%，因此選擇劃分網(wǎng)格數(shù)為956036;對于DN200管道,選取網(wǎng)格數(shù)分別為2308874與4328293時,在各點測出的壓力相差均小于0.01%，因此選擇劃分網(wǎng)格數(shù)為2308874。
1.2模擬程序參數(shù)設(shè)置
1.2.1物性參數(shù)設(shè)置
Refprop軟件由NIST開發(fā)，該軟件含有豐富的數(shù)據(jù)庫以及適用于超臨界CO,的多個狀態(tài)方程。本文通過在Fluent軟件中激活NISTRealGas模型[川進行調(diào)用,計算超臨界二氧化碳流體的物性參數(shù)。其中物性參數(shù)采用FEK狀態(tài)方程模型計算，黏度采用VS1模型,導(dǎo)熱系數(shù)采用TC1模型,各模型的相關(guān)參數(shù)如表2所示。
1.2.2邊界條件設(shè)置
本文針對超臨界二氧化碳鍋爐人口處的循環(huán)工質(zhì)進行流量系數(shù)測量的數(shù)值模擬研究,設(shè)置了質(zhì)量流量人口與壓力出口,溫度、壓力等參數(shù)的選取為超臨界二氧化碳鍋爐入口處工質(zhì)典型參數(shù),即750K、21MPa。由于循環(huán)系統(tǒng)運行于高壓環(huán)境，管道的壓力損失相較而言很小，因此可認(rèn)為管道壓力為恒定.值,壓力出口參數(shù)設(shè)置與人口相同,其余參數(shù)保持默認(rèn)不變;由于超臨界二氧化碳鍋爐人口處管道一般采取嚴(yán)格保溫措施,因此可忽略壁面與工質(zhì)間的換熱,設(shè)置為絕熱邊界。
1.2.3數(shù)值模型設(shè)置
本文主要模擬超臨界二氧化碳工質(zhì)流經(jīng)孔板流量計的流動過程,基本方程包含質(zhì)量、動量和能量輸運方程,由雷諾數(shù)的定義公式
計算可知本文針對的超臨界二氧化碳工況下雷諾數(shù)均遠大于4000,因此管道內(nèi)的流動均處于湍流狀態(tài),在進行數(shù)值模擬時需進行湍流模型的設(shè)置,本文選擇SSTh-ɷ湍流模型。
1.2.4求解參數(shù)設(shè)置
FLUENT中的亞松弛因子主要控制計算過程中每次迭代的變化量,可以通過減少兩層次之間計算.結(jié)果的差值從而促進收斂。本文設(shè)置的亞松弛因子如表3所示。
1.3模型的驗證
基于上述設(shè)置，本文針對溫度為535.1~642.5K、壓力為19MPa、質(zhì)量流量為1.28kg/s的實驗工況進行了模擬研究,模擬的孔板結(jié)構(gòu)參數(shù)、溫度、壓力、流量等參數(shù)以及數(shù)據(jù)處理方法均與實驗保持--致,得到了超臨界二氧化碳工質(zhì)的流量系數(shù)。將模擬計算得到的流量系數(shù)與實流測量結(jié)果進行對比,結(jié)果如圖4所示。,通過數(shù)值模擬得到的流量系數(shù)與實驗數(shù)據(jù)總體趨勢相似,在數(shù)值上均高于實驗數(shù)據(jù),但相對于實驗數(shù)據(jù)的偏差較小,偏差為1.62%~2.69%。
造成偏差的原因可能有多種,如實驗選用測量儀表具有一定的不確定度、模擬參數(shù)的設(shè)置無法與真實情況完全對應(yīng)、收斂判據(jù)設(shè)置不嚴(yán)格等等。為了降低模擬結(jié)果與實驗數(shù)據(jù)的偏差,本文分別按各模擬結(jié)果相對同工況下實驗數(shù)據(jù)的偏差平均值進行修正。模擬得到流量系數(shù)相對實驗數(shù)據(jù)平均增大0.013,因此對模擬結(jié)果減去該修正值,修正后相對偏差為0.016%~0.674%。
上述結(jié)果說明數(shù)值模擬方法與實驗結(jié)果的一致性較好,因此本文建立的數(shù)值模擬方法可用于后續(xù)進一步的研究。
2孔板結(jié)構(gòu)參數(shù)對流量系數(shù)影響
2.1直徑比的影響
直徑比會顯著影響孔板對于介質(zhì)流過的節(jié)流效果,改變介質(zhì)流過的速度、壓力等參數(shù),是影響孔板流量計測量性能的首要因素。ISO國際標(biāo)準(zhǔn)中規(guī)定,孔板流量計的直徑比一般在0.1~0.75內(nèi)變化,本文分別選取直徑比在0.3~0.9之內(nèi)的7個工況進行了模擬計算,探究孔板流量計直徑比對流量系數(shù)的影響，得到的結(jié)果如圖5所示。
對數(shù)據(jù)進行分析可知:
(1)孔板流量系數(shù)隨直徑比的變化趨勢與管徑無關(guān)。隨著孔板直徑比增大,DN25和DN200管道內(nèi)孔板流量系數(shù)呈現(xiàn)近似相同變化趨勢;上升-平穩(wěn)-.上升,主要區(qū)別在于前者在β為0.4~0.8范圍內(nèi)較平穩(wěn)、而后者在0.5~0.8范圍內(nèi)較平穩(wěn);
(2)標(biāo)準(zhǔn)文件[10]中經(jīng)驗計算公式的結(jié)果隨直徑比增加而逐漸下降,其中直徑比在0.3~0.6范圍內(nèi)時下降趨勢較平緩,當(dāng)超過0.6時下降值逐漸增大;
(3)孔板流量系數(shù)在β為0.3~0.6時小于經(jīng)驗公式計算值,此范圍內(nèi)使用經(jīng)驗計算公式會使測量結(jié)果較真實值大2.45%~47.03%;β在0.7~0.9時孔板流量系數(shù)大于經(jīng)驗公式計算值，此時使用經(jīng)驗計算公式會使測量值比實際值小0.5%~60.19%;
(4)當(dāng)直徑比在標(biāo)準(zhǔn)文件[10]規(guī)定的0.1~0.75范圍內(nèi)時,孔板流量系數(shù)的模擬結(jié)果與經(jīng)驗公式計算結(jié)果的相對誤差波動較大,如對DN25管道而言,β為0.3時相對偏差達到47.03%,而β為0.7時相對偏差僅為0.5%。因此對于超臨界二氧化碳工質(zhì)而言，孔板直徑比的選擇范圍應(yīng)較標(biāo)準(zhǔn)規(guī)定范圍縮小;對于超臨界二氧化碳工質(zhì)而言,直徑比在0.6~0.7范圍內(nèi)時孔板流量系數(shù)的模擬結(jié)果與經(jīng)驗公式計算結(jié)果的相對誤差較小，其中DN25管道相對誤差為0.5%~2.45%,DN200管道相對誤差為2.27%~3.6%。
2.2節(jié)流孔厚度的影響
孔板節(jié)流孔厚度決定了超臨界二氧化碳工質(zhì)流過收縮管道的長度,是影響孔板節(jié)流能力的主要參數(shù)之一，會對工質(zhì)流過孔板的流速、壓力等參數(shù)產(chǎn)生影響,進而影響測量結(jié)果。標(biāo)準(zhǔn)文件[10]規(guī)定標(biāo)準(zhǔn)孔板節(jié)流孔厚度應(yīng)在0.005D~0.02D之間,對應(yīng)DN25管道的e應(yīng)在0.115~0.46mm,DN200管道的e應(yīng)在0.695~2.78,本文分別模擬了DN25管道e為0.1~0.7mm、DN200管道e為0.2~4.2mm時超臨界二氧化碳工質(zhì)流過節(jié)流孔板的流量系數(shù)變化,為便于對比,以e/D為橫坐標(biāo)將結(jié)果表示在圖6中。
對數(shù)據(jù)進行分析可知:
(1)孔板流量系數(shù)隨e/D的變化趨勢與管徑無關(guān)。隨著e/D逐漸增加,DN25和DN200管道內(nèi)孔板的流量系數(shù)均呈現(xiàn)先減小后增大的趨勢，分別在e/D為0.017和0.023時達到最小值，此后流量系數(shù)先急劇增大,隨后保持平緩增長;
(2)標(biāo)準(zhǔn)文件[I0]中經(jīng)驗計算公式的結(jié)果不隨節(jié)流孔厚度而發(fā)生改變,其中DN25管道的經(jīng)驗公式計算結(jié)果略大-一些,模擬得到的DN25和DN200管道的流量系數(shù)均小于經(jīng)驗公式計算結(jié)果,其中前者的相對誤差為0.18%~1.84%,后者的相對誤差為0.31%~2.05%;
(3)在標(biāo)準(zhǔn)文件[10]規(guī)定孔板節(jié)流孔厚度范圍內(nèi),孔板流量系數(shù)模擬結(jié)果與經(jīng)驗公式的相對誤差均在2%以下,因此標(biāo)準(zhǔn)中規(guī)定的孔板節(jié)流孔厚度范圍可以接受;同時還發(fā)現(xiàn)當(dāng)節(jié)流孔厚度超過規(guī)定范圍一定值后,相對誤差仍可接受,甚至相對誤差還可能減少，如DN25管道的e為0.6mm、0.7mm時,均超出了規(guī)定上限0.46mm,但相對誤差分別達到了0.3%和0.18%，因此標(biāo)準(zhǔn)中規(guī)定的節(jié)流孔厚度范圍在針對超臨界二氧化碳工質(zhì)時可以適當(dāng)擴大,推薦DN25管道孔板節(jié)流孔厚度可在0.004D~0.03D內(nèi)變化，DN200管道在0.005D~0.03D范圍內(nèi);
(4)基于模擬結(jié)果給出相對誤差更小時對應(yīng)孔板節(jié)流孔厚度的推薦值,其中DN25管道孔板在e/D為0.004~0.008及0.02~0.03之間,即e為0.1~0.2mm.0.5~0.7mm時，相對誤差小于1.5%;DN200管道孔板在e/D為0.005~0.012及0.027~0.03時,對應(yīng)e為0.7~1.7mm及3.7~4.2mm時，相對誤差小于等于1.5%。
2.3孔板厚度的影響
由圖1可知，標(biāo)準(zhǔn)孔板在節(jié)流孔之后還設(shè)置一定長度的錐形擴流段,與節(jié)流孔段共同組成孔板的節(jié)流段,該擴流段長度也會對孔板的節(jié)流能力產(chǎn)生影響,從而改變工質(zhì)流過孔板后的速度、壓力等參數(shù),對測量精度產(chǎn)生影響。標(biāo)準(zhǔn)文件[10]規(guī)定孔板厚度E應(yīng)在e~0.05D之間,對應(yīng)DN25管道的E應(yīng)不大于1.15mm,DN200管道的E不超過6.95mm。
本文在保持節(jié)流孔厚度不變的情況下，分別設(shè)置了不同的孔板厚度用以探究流量系數(shù)的變化,其中DN25管道孔板厚度E為0.5~1.4mm,DN200管道孔板厚度E為3~8mm,模擬結(jié)果如圖7所示。
對數(shù)據(jù)進行分析可知:
(1)孔板流量系數(shù)隨E/D的變化趨勢與管徑無關(guān)。隨著E/D逐漸增加,DN25管道和DN200管道內(nèi)孔板流量系數(shù)呈現(xiàn)近似相同的變化趨勢:即下降上升-平穩(wěn)-下降,主要區(qū)別在于DN200管道內(nèi)孔板流量系數(shù)下降和上升的趨勢更加明顯;
(2)流量系數(shù)經(jīng)驗計算公式的結(jié)果不隨孔板厚.度而發(fā)生變化，其中DN25管道的經(jīng)驗公式計算結(jié)果偏大--些,DN25和DN200管道的流量系數(shù)均小于經(jīng)驗公式計算值,因此當(dāng)使用經(jīng)驗公式進行工質(zhì)流量計算時會造成計算結(jié)果偏大;
(3)在標(biāo)準(zhǔn)文件[10]規(guī)定孔板厚度范圍內(nèi),DN25和DN200管道內(nèi)孔板流量系數(shù)與經(jīng)驗計算公式的相對誤差均在2%以下，因此標(biāo)準(zhǔn)中的規(guī)定范圍可以接受,但該規(guī)定范圍對于超臨界二氧化碳工質(zhì)可適當(dāng)擴充,如模擬結(jié)果所示，當(dāng)DN25和DN200兩種管徑的孔板厚度E達到0.06D左右時，雖然已經(jīng)超出了規(guī)定的0.05D這一限值,但相對誤差仍小于2%，處于可接受的范圍,但依據(jù)變化趨勢可以合理預(yù)測，當(dāng)孔板厚度繼續(xù)增加時，相對誤差將大于2%，因此建議對于超臨界二氧化碳工質(zhì)而言,孔板.厚度可設(shè)置在0.02D~0.06D之間。
2.4孔板流量計結(jié)構(gòu)參數(shù)設(shè)計建議
通過對孔板流量計各結(jié)構(gòu)參數(shù)的模擬研究,明確了在進行超臨界二氧化碳工質(zhì)質(zhì)量流量測量時,孔板流量系數(shù)隨各結(jié)構(gòu)參數(shù)的變化趨勢與相對誤差,本節(jié)主要對以上模擬結(jié)果進行總結(jié)分析,參考《用能單位能源計量器具配備和管理通則》中的精度規(guī)定，，給出了針對超臨界二氧化碳工質(zhì)的孔板流量計結(jié)構(gòu)參數(shù)推薦設(shè)計范圍,在該范圍內(nèi)經(jīng)驗計算公式的計算結(jié)果可滿足2.5精度等級要求,還進一步提出了該范圍內(nèi)精度相對更高的結(jié)構(gòu)參數(shù)推薦值，將以上結(jié)果與現(xiàn)行國際標(biāo)準(zhǔn)ISO5167-2:2003中標(biāo)準(zhǔn)孔板流量計各結(jié)構(gòu)參數(shù)的規(guī)定范圍進行對比,如表4所示。
可以看出,對于超臨界二氧化碳工質(zhì)而言,標(biāo)準(zhǔn)文件[10]規(guī)定的孔板流量計各結(jié)構(gòu)參數(shù)的設(shè)計范圍并不完全適用,其中直徑比的規(guī)定范圍過大,對應(yīng)的流量系數(shù)的相對誤差也波動較大,從0.5%到47.03%不等,而當(dāng)直徑比在0.6~0.7范圍內(nèi)時,可將相對誤差有效降低至0.5%~3.6%;在標(biāo)準(zhǔn)文件[10]規(guī)定的節(jié)流孔厚度、孔板厚度等參數(shù)范圍內(nèi),絕大多數(shù)流量系數(shù)的相對誤差可控制在2%以下，因此其規(guī)定范圍可以繼續(xù)使用,同時本文的數(shù)值模擬結(jié)果顯示，當(dāng)孔板的以，上幾個結(jié)構(gòu)參數(shù)的數(shù)值超.出其規(guī)定范圍時,最大相對誤差也僅為2%左右，因此對于超臨界二氧化碳工質(zhì)而言，孔板的節(jié)流孔厚度、孔板厚度等參數(shù)均可一定程度上超出標(biāo)準(zhǔn)中的規(guī)定范圍，相對誤差也可接受。
3入口直角邊緣尖銳度及其修正系數(shù)的模擬研究
一般而言,孔板人口邊緣應(yīng)該是尖銳的,其與超臨界二氧化碳工質(zhì)首先直接接觸,如果其尖銳度不夠的話則無法保證對于工質(zhì)的節(jié)流作用達到預(yù)期,因而會對測量精度產(chǎn)生影響。在孔板實際工作過程:中,可能存在加工精度不足、工質(zhì)磨損、腐蝕等問題的存在,造成直角邊緣變鈍,故標(biāo)準(zhǔn)文件[10]規(guī)定,孔板人口邊緣的圓弧半徑應(yīng)小于等于0.0004D,在此限值之內(nèi)的誤差是可以接受的,若超過這一-限值,則無法保證測量精度,應(yīng)進行相應(yīng)的維修、更換或修正等。本文模擬了孔板人口邊緣圓弧半徑為0~0.015D時孔板的流量系數(shù)變化趨勢,結(jié)果如圖8所示。
分析結(jié)果可以得出:
(1)孔板人口邊緣尖銳度對于孔板流量系數(shù)的影響趨勢與管徑無關(guān)。隨著孔板人口邊緣逐漸變鈍,DN25和DN200兩種管道的孔板流量系數(shù)呈現(xiàn)近似相同的變化趨勢，均隨著人口圓弧半徑的增大而先增大后減小,分別在r達到0.01D和0.008D時流量系數(shù)達到最大;
(2)孔板人口邊緣開始鈍化時,流量系數(shù)顯著增加,遠大于經(jīng)驗計算公式結(jié)果,因此造成使用經(jīng)驗計算公式時得到的工質(zhì)質(zhì)量流量相對真實值很小，其中DN25管道內(nèi)相對誤差為5.16%~12.61%,DN200管道的相對誤差為5.13%~11.96%;
(3)對于超臨界二氧化碳工質(zhì)而言,當(dāng)孔板人口直角邊緣變鈍后,應(yīng)立即進行相應(yīng)的處理或流量系數(shù)的修正,否則誤差將會變得很大。
當(dāng)使用標(biāo)準(zhǔn)文件[10]給出的不同邊緣尖銳度對應(yīng)的修正系數(shù)b進行孔板流量系數(shù)的修正時,可以--定程度上減少流量系數(shù)的相對誤差。本文使用表5中的修正系數(shù)b對經(jīng)驗公式計算結(jié)果進行修正,結(jié)果如圖9所示。
結(jié)果顯示，對于超臨界二氧化碳工質(zhì)而言,當(dāng)使用修正系數(shù)b進行孔板流量系數(shù)的修正后,僅可使部分邊緣圓弧半徑對應(yīng)的孔板流量系數(shù)相對誤差降低到可以接受的程度,而大部分情況下相對誤差仍比較大,如DN25管道多數(shù)情況下的流量系數(shù)相對誤差在4.41%~6.94%之間,DN200管道的相對誤差多數(shù)在3.74%~7.11%,因此該修正系數(shù)b對于超臨界二氧化碳工質(zhì)并不適用。
對DN25和DN200管道的模擬流量系數(shù)及經(jīng)驗公式計算結(jié)果求平均,得到針對超臨界二氧化碳工質(zhì)的不同孔板邊緣尖銳度對應(yīng)的修正系數(shù)b,如表6所示。將更正的修正系數(shù)應(yīng)用于模擬數(shù)據(jù)，結(jié)果如圖10所示。
可以看出,當(dāng)使用更正后的修正系數(shù)b進行孔板的流量系數(shù)經(jīng)驗公式計算結(jié)果的修正后,得到的流量系數(shù)與模擬結(jié)果擬合較好,其中DN25管道相對誤差為0.11%~1.21%,DN200管道相對誤差為0.11%~1.85%。
4結(jié)論
本文開展了孔板流量計的數(shù)值模擬研究，探究了孔板的各結(jié)構(gòu)參數(shù)對超臨:界二氧化磯工質(zhì)流量系數(shù)的影響,基于此給出了孔板流量計結(jié)構(gòu)參數(shù)的設(shè)計建議,并且探究了人口直角邊緣尖銳度對流量系數(shù)的影響,得到的主要結(jié)論如下:
(1)現(xiàn)行標(biāo)準(zhǔn)文件中的孔板流量計結(jié)構(gòu)參數(shù)的規(guī)定范圍測量相對誤差在0.5%~47%的較大范圍內(nèi)波動,對于超臨界二氧化碳工質(zhì)并不適用。
(2)本文針對超臨界二氧化碳工質(zhì)提出了孔板.流量計結(jié)構(gòu)參數(shù)推薦設(shè)計范圍,其中直徑比應(yīng)為0.6~0.7,節(jié)流孔厚度應(yīng)為0.004~0.03倍的管道內(nèi)徑,孔板厚度應(yīng)為0.02~0.06倍的管道內(nèi)徑,在該范圍內(nèi)絕大多數(shù)工況下流量系數(shù)的相對誤差可控制在2%以下;
(3)孔板人口邊緣鈍化會使流量系數(shù)顯著增加,且修正系數(shù)b并不能使相對誤差降低至可以接受的范圍,修正后相對誤差仍有約3.74%-7.11%,本文針對不同工況提出不同修正參數(shù),修正后經(jīng)驗公式的相對誤差降低為0.11%~1.85%。
本文來源于網(wǎng)絡(luò),如有侵權(quán)聯(lián)系即刪除！
