摘 要：以50 mm口徑旋進(jìn)旋渦流量計(jì)為研究對象，利用計(jì)算流體動(dòng)力學(xué)軟件FLUENT,采用歐拉- 歐拉模型，結(jié)合RNG k-ε湍流模型與SIMPLEC算法對流量計(jì)內(nèi)部進(jìn)行氣固兩相數(shù)值模擬和實(shí)驗(yàn)研究，對旋進(jìn)旋渦流量計(jì)氣固兩相流動(dòng)中固相的體積分布規(guī)律、速度場分布、壓力場分布進(jìn)行了探討，并做了相應(yīng)的實(shí)驗(yàn)研究。結(jié)果表明：顆粒濃度越大,監(jiān)測點(diǎn)的壓力頻率越小，流量計(jì)讀數(shù)越小，流體流動(dòng)越復(fù)雜，在流道壁面上的固體顆粒分布越多。

  多相流廣泛地存在各大工業(yè)領(lǐng)域與自然界，而氣固兩相流動(dòng)是其中重要的組成部分，有很多的中外學(xué)者對氣固兩相流現(xiàn)象進(jìn)行了大量的研究[9]。林金賢等[10]利用數(shù)值模擬的方法，分析了彎管內(nèi)的氣固兩相流現(xiàn)象，指出了顆粒大小與顆粒濃度對固體顆粒在彎管內(nèi)分布及單位管長壓力損失的影響。楊鈺辰等[11]利用CFD方法對電除塵器進(jìn)行數(shù)值模擬研究，比較全面地分析了氣固兩相流動(dòng)對電除塵器效率的影響。蔣夢婷等[12]利用數(shù)值模擬計(jì)算與實(shí)驗(yàn)驗(yàn)證的方法對旋風(fēng)分離器內(nèi)部流動(dòng)及顆粒運(yùn)動(dòng)特性進(jìn)行了研究，并分析了顆粒的受力情況對旋風(fēng)分離器性能的影響。

1、數(shù)學(xué)模型與計(jì)算方法：
1.1、物理模型：
  本次研究的對象是口徑為50 mm的旋進(jìn)旋渦流量計(jì)，流道的二維幾何模型以及消旋器的二維圖如圖1所示；起旋器、消旋器以及全流域幾何三維模型如圖2所示。
1.2、網(wǎng)格劃分：
  本文利用ICEM軟件對流量計(jì)模型進(jìn)行網(wǎng)格劃分，入口延伸段、收縮段、喉部及出口延伸段采用結(jié)構(gòu)網(wǎng)格劃分，由于起旋器部分結(jié)構(gòu)比較復(fù)雜，故采用非結(jié)構(gòu)網(wǎng)格劃分，如圖3所示。各部分網(wǎng)格數(shù)與節(jié)點(diǎn)數(shù)分別為：入口延伸段45 623、128 600，起旋器部分297 244、74 898，喉部74 377、186 428，出口延伸段384 903、1 029 028。

1.3、計(jì)算方法：
  進(jìn)口邊界條件選擇速度進(jìn)口（velocity），出口邊界條件設(shè)為壓力出口，其值設(shè)為1個(gè)標(biāo)準(zhǔn)大氣壓；壁面邊界條件采用無滑移邊界條件，在流量計(jì)管道近壁面處采用標(biāo)準(zhǔn)壁面函數(shù)。本次計(jì)算采用的是歐拉-歐拉方法，選取RNG k-ε湍流模型與SIMPLEC算法相結(jié)合的計(jì)算方法。計(jì)算中的湍流模型采用默認(rèn)的混合湍流模型，該模型中的常數(shù)項(xiàng)均采用默認(rèn)的數(shù)值。在計(jì)算方法中各動(dòng)量參數(shù)與體積分?jǐn)?shù)計(jì)算均采用一階迎風(fēng)格式。

圖 3 網(wǎng)格劃分

  流量計(jì)對流量Q=10 m3/h、Q=22.5 m3/h、Q=37.5 m3/h、Q=60 m3/h、Q=105 m3/h、Q=150 m3/h六個(gè)工況點(diǎn)進(jìn)行了標(biāo)定實(shí)驗(yàn)，實(shí)驗(yàn)結(jié)果與數(shù)值模擬結(jié)果如圖4所示。由圖4可以看出，在誤差允許范圍內(nèi)，模擬結(jié)果與實(shí)驗(yàn)結(jié)果是吻合的，說明本次模擬采用的模型是正確的。

圖4 標(biāo)定結(jié)果

2.2、氣體單相的數(shù)值模擬：
  本次模擬對工況點(diǎn)Q=60 m3/h、Q=45 m3/h、Q=37.5 m3/h進(jìn)行了數(shù)值模擬，模擬采用RNG k-ε湍流模型與SIPPLEC算法相結(jié)合的數(shù)值計(jì)算方法，監(jiān)測點(diǎn)P點(diǎn)的壓力脈動(dòng)如圖5所示。

圖 5 P 點(diǎn)壓力脈動(dòng)
  在Q=60 m3/h時(shí)，流量計(jì)橫截面與監(jiān)測點(diǎn)截面一個(gè)周期的壓力云圖如圖6、圖7所示。由圖5可以看出隨著工況點(diǎn)的增大流量計(jì)產(chǎn)生的壓力脈動(dòng)頻率增大，流量越大，氣體的速度就越大，在監(jiān)測點(diǎn)的壓力值就越大。從流體壓力云圖6、圖7可以看出，流量計(jì)渦核的壓力是***小的，沿著壁面逐漸增加。橫截面在t=0.05 s時(shí)的速度云圖如圖8所示。

圖8 一個(gè)周期內(nèi)速度云圖（m/s）圖9 一個(gè)周期內(nèi)流線圖

2.3、氣固兩相數(shù)值模擬：
  本次模擬計(jì)算的固相為玻璃粉，平均粒徑d=0.104mm，密度ρ2=2600 kg/m3，固氣質(zhì)量比濃度為C=0.21、C=0.42、C=0.63、C=0.84四種。當(dāng)Q=60 m3/h時(shí)，監(jiān)測點(diǎn)的壓力脈動(dòng)如圖10所示。

圖10 Q=60 m3/h時(shí)P點(diǎn)氣體壓力脈動(dòng)

圖11 各工況點(diǎn)頻率變化圖圖 12 橫截面壓力云圖（Pa）
  當(dāng)Q=60 m3/h時(shí)，監(jiān)測點(diǎn)P點(diǎn)在t=0.035 64 s時(shí)刻截面氣體單相模擬與氣固兩相模擬的氣相壓力云圖如圖14所示，由于固體顆粒對流量計(jì)管道壁有碰撞作用，氣固兩相流對監(jiān)測點(diǎn)的壓力要大一些。由云圖可知，在旋進(jìn)旋渦流量計(jì)的喉部與擴(kuò)張段的管道中心氣固兩相氣流的壓力要比氣體單相的壓力整體偏大，氣流的低壓區(qū)要比氣體單項(xiàng)時(shí)更多，氣流的壓力梯度更為復(fù)雜。當(dāng)Q=60 m3/h，固氣體積比C=0.01時(shí)在一個(gè)周期內(nèi)氣體在旋進(jìn)旋渦流量計(jì)流道內(nèi)的體積分布與監(jiān)測點(diǎn)橫截面的體積分布情況如圖15、圖16所示。由圖15、圖16可知，氣體的體積分布跟氣體的流動(dòng)規(guī)律一樣，在流量計(jì)收縮段與喉部的旋渦流旋渦中心分布

圖13 P點(diǎn)截面壓力云圖圖14 t=0.035 64 s時(shí)刻P點(diǎn)截面的壓力云圖（Pa）圖15 橫截面氣相體積分布圖圖16 P點(diǎn)截面氣相體積分布云圖
***高。固體顆粒由于與流量計(jì)流道壁面以及起旋器壁面之間相互力的作用，有很大一部分固體顆粒分布在流量計(jì)流道壁面與起旋器壁面。

3、實(shí)驗(yàn)研究：
  為了與數(shù)值模擬計(jì)算進(jìn)行對比分析，本文對流量計(jì)進(jìn)行了氣固兩相實(shí)驗(yàn)，本次是在一個(gè)氣力輸送的管網(wǎng)系統(tǒng)中進(jìn)行的，固體顆粒為玻璃粉。數(shù)值模擬與實(shí)驗(yàn)結(jié)果對比分析結(jié)果如圖17所示。由圖17可知，隨著顆粒濃度的增大，顆粒對實(shí)驗(yàn)結(jié)果的影響越大，數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果之間的差值有增大的趨勢，但在誤差允許范圍內(nèi)，數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果是相符合的。

4、結(jié)論：
1）隨著顆粒濃度的增加，旋進(jìn)旋渦流量計(jì)各工況點(diǎn)在監(jiān)測點(diǎn)P點(diǎn)的壓力脈動(dòng)頻率呈現(xiàn)線性下降的趨勢。當(dāng)固氣質(zhì)量比濃度C=0時(shí)（氣體單相時(shí)）氣流的脈動(dòng)頻率值***大。
2）當(dāng)氣流流入旋進(jìn)旋渦流量計(jì)的收縮段時(shí)，氣流流速會(huì)增大，并形成旋渦流，經(jīng)過擴(kuò)展段時(shí)氣流流速會(huì)急劇減小，旋渦流渦核繞著流量計(jì)軸線做順時(shí)針運(yùn)動(dòng)，此時(shí)旋渦流的流動(dòng)***為復(fù)雜，當(dāng)流過消旋器后，旋渦流的流動(dòng)逐漸平穩(wěn)，旋渦逐漸消失，***后氣流平穩(wěn)地流過流量計(jì)出口。
3）流量計(jì)的氣固兩相流動(dòng)中，氣體的體積分布與氣流的流動(dòng)規(guī)律一致，在收縮段與喉部主要分布在旋渦中心，有很大一部分固體顆粒分布在流量計(jì)流道壁面以及起旋器壁面。