趙

(上海理工大學(xué) 機(jī)械工程學(xué)院，上海 200093)

基于Fluent的流量標(biāo)準(zhǔn)裝置穩(wěn)壓罐流場特性研究

(上海理工大學(xué) 機(jī)械工程學(xué)院，上海 200093)

流量標(biāo)準(zhǔn)裝置是檢測和標(biāo)定流量儀表的重要設(shè)備，其中穩(wěn)壓裝置的設(shè)計(jì)及流場分析關(guān)系到是否可以保證裝置的標(biāo)定精確度。文中基于流場數(shù)值模擬基礎(chǔ)，利用Fluent軟件對穩(wěn)壓罐的內(nèi)部流場進(jìn)行了仿真分析，經(jīng)過迭代計(jì)算，得到了穩(wěn)壓罐出口處的流量波動系數(shù)，然后分析了穩(wěn)壓罐罐內(nèi)隔板結(jié)構(gòu)對穩(wěn)壓性能的影響。結(jié)果表明，在兼顧功能和經(jīng)濟(jì)性的要求下，具有3塊橫隔板的穩(wěn)壓罐結(jié)構(gòu)比較合理，同時，豎隔板位于穩(wěn)壓罐前部1/3處穩(wěn)壓特性最佳。

Fluent；流量標(biāo)準(zhǔn)裝置；穩(wěn)壓罐；流場分析

流量儀表被廣泛應(yīng)用于農(nóng)業(yè)、能源和汽車領(lǐng)域，而新出廠的或者經(jīng)過維修的流量儀表如各類流量計(jì)，必須由流量標(biāo)準(zhǔn)裝置進(jìn)行標(biāo)定之后方可投入使用。在大口徑高性能的流量標(biāo)準(zhǔn)裝置中，通常使用穩(wěn)壓罐來穩(wěn)定管道及被檢表前后段的壓力波動，從而保證裝置標(biāo)定的精確度。本文介紹了針對此類穩(wěn)壓罐體的流場流場特性研究，及在此基礎(chǔ)上的穩(wěn)壓特性仿真對比。

1 流場的數(shù)值模擬基礎(chǔ)

1.1 計(jì)算流體動力學(xué)和Fluent軟件

計(jì)算流體動力學(xué)(Computational Fluid Dynamics，CFD)[1]是通過計(jì)算機(jī)數(shù)值計(jì)算和圖像顯示，對包含有流體流動和熱傳導(dǎo)等相關(guān)物理現(xiàn)象的系統(tǒng)所做的分析。CFD的基本思想可以歸結(jié)為[2-3]：把原來在時間域及空間域上連續(xù)的物理量的場，例如速度場和壓力場，用一系列有限個離散點(diǎn)上的變量值的集合來代替，通過一定的原則和方式建立起關(guān)于這些離散點(diǎn)上場變量之間關(guān)系的代數(shù)方程組，然后求解代數(shù)方程組獲得場變量的近似值。

本文研究主要基于Fluent軟件，此軟件是CFD軟件中相對成熟和應(yīng)用最為廣泛的，用于模擬和分析復(fù)雜幾何區(qū)域內(nèi)的流體流動與傳熱現(xiàn)象[4]，能夠提供靈活的網(wǎng)絡(luò)特性，可以支持多種網(wǎng)絡(luò)，且以用戶界面友好而著稱。

1.2 流體動力學(xué)控制方程

流體流動要受物理守恒定律的支配，基本的守恒定律包括：質(zhì)量守恒定律、動量守恒定律、能量守恒定律。對這些守恒定律的數(shù)學(xué)描述就是控制方程[5]。

由于本文所述研究應(yīng)用的數(shù)學(xué)模型包含應(yīng)力項(xiàng)，故簡要介紹動量守恒方程。動量守恒定律即微元體中流體的動量對時間的變化率等于外界作用在該微元體上的各種力之和。

(1)

式中，ρ為流體密度，t為時間，p為流體的靜壓強(qiáng)，u、v和w是速度矢量u在x、y和z方向上的分量，Su、SV和Sw是動量守恒方程的廣義源項(xiàng)。

而上述3種定律對應(yīng)的控制方程可以用如下通用形式表示

(2)

式中，u為速度矢量；φ是通用變量，可代表u、v、w和T等求解變量；Γ是廣義擴(kuò)散系數(shù)；S是廣義源項(xiàng)。3個守恒方程的相關(guān)符號和各個方程對應(yīng)關(guān)系如表1所示。

表1 通用控制方程中各符號的具體形式

其中，cp為比熱容；T為溫度；k為流體的傳熱系數(shù)；ST為流體的內(nèi)熱源及由于粘性作用使流體機(jī)械能轉(zhuǎn)換為熱能的部分，簡稱粘性耗散項(xiàng)。

1.3 湍流數(shù)值模擬方法

對于流體流動，定義Reynolds數(shù)(雷諾數(shù))

(3)

式中，u為液體速度；ν為運(yùn)動粘度；R為水力半徑。

自然界中的流體流動狀態(tài)主要有兩種形式，即層流和湍流。層流是指流體在流動過程中兩層之間沒有相互混摻，而湍流[6]是指流體不處于分層流動狀態(tài)。一般說來，湍流是普遍的，層流則屬于個別情況。

當(dāng)雷諾數(shù)小于2 300時，流體一定為層流；當(dāng)雷諾數(shù)大于8 000～12 000時，流體一定為湍流；當(dāng)雷諾數(shù)在2 300～8 000之間時，流動處于層流和湍流之間。

本文的研究基于大口徑水流量標(biāo)準(zhǔn)裝置，流體狀態(tài)為湍流。湍流的數(shù)值模擬方法可以分為直接數(shù)值模擬方法和非直接數(shù)值模擬方法[7]。直接數(shù)值模擬方法是指求解瞬時湍流控制方程。非直接數(shù)值模擬方法則是對湍流作某種程度的近似和簡化處理，非直接數(shù)值模擬方法分為大渦模擬、統(tǒng)計(jì)平均法、雷諾時均(RANS，Reynolds)法。目前，直接數(shù)值模擬法、大渦模擬法和RANS法是目前求解湍流問題的3種主要方法，其中RANS法是目前應(yīng)用最多、范圍最廣的數(shù)值模擬求解方法[8]。

RANS法的核心[9]是不直接求解瞬時的控制方程，而是想辦法求解時均化的雷諾方程，即將瞬態(tài)的動量通過某種模型在時均化的方程中體現(xiàn)出來。與控制方程相比，RANS方程多出了雷諾應(yīng)力項(xiàng)，雷諾應(yīng)力項(xiàng)可用來表示湍流效應(yīng)，為使方程組可解，須做出合適假設(shè)，并建立相關(guān)模型。

根據(jù)對雷諾應(yīng)力做出的假定或處理方式的不同，目前常用的湍流模型有兩大類[10]：Reynolds應(yīng)力模型和渦粘膜型。Reynolds應(yīng)力模型包括Reynolds應(yīng)力方程模型和代數(shù)應(yīng)力方程模型。渦粘膜型包括零方程模型、一方程模型和兩方程模型。其中，兩方程模型在工程中應(yīng)用最為廣泛，最基本的兩方程模型是標(biāo)準(zhǔn)k-ε模型，此外，還有RNGk-ε模型和Realizablek-ε模型。

RNG k-ε模型適合于求解非均勻湍流問題，如旋轉(zhuǎn)流動、較大速度梯度的流場等，它的應(yīng)用范圍更廣，計(jì)算精度更高。湍動能k的輸運(yùn)方程和耗散率ε如式(4)所示。

(4)

在Fluent中，有3種多相流模型，分別為VOF模型，混合模型和歐拉模型[11]。其中，VOF模型適合于分層或自由表面流，而混合模型和歐拉模型適合于流動中有相混合或分離，或者分散相的體積分?jǐn)?shù)超過10%的情形，此模型也是本文研究在軟件中應(yīng)用的模型。

2 穩(wěn)壓罐的流場仿真

2.1 計(jì)算模型和網(wǎng)格處理

數(shù)值模擬[12]的可行性、收斂性和計(jì)算精度與計(jì)算模型的結(jié)構(gòu)和網(wǎng)格劃分質(zhì)量有直接的關(guān)系。在建立計(jì)算域模型以及對模型進(jìn)行網(wǎng)格劃分時，要充分考慮可以接受的計(jì)算機(jī)承受能力和計(jì)算精度，選擇合適的建模方式。

考慮到計(jì)算機(jī)的性能，本研究對穩(wěn)壓罐橫隔板上的孔進(jìn)行簡化處理，即增大孔的直徑和減少孔的數(shù)量，同時保證一定的流通面積，使用Solidworks建立其計(jì)算域模型，根據(jù)尺寸建立穩(wěn)壓罐的計(jì)算域模型如圖1所示。建立好穩(wěn)壓罐的計(jì)算域模型之后，導(dǎo)入Gambit中劃分網(wǎng)格。如圖2所示，流體流動方向與網(wǎng)格線方向一致，為了提高數(shù)值模擬結(jié)果的穩(wěn)定性和加快迭代收斂，把穩(wěn)壓罐網(wǎng)格進(jìn)行分區(qū)域劃分，在封頭和橫隔板部分使用非結(jié)構(gòu)性網(wǎng)格，在其它部分采用結(jié)構(gòu)性網(wǎng)格。

2.2 確定邊界類型

在建立模型計(jì)算時，假設(shè)進(jìn)入穩(wěn)壓罐的水為純水，不考慮進(jìn)口水流混入少量氣泡的情況，將穩(wěn)壓罐左側(cè)進(jìn)水管的入口面設(shè)為質(zhì)量流量入口(mass-flow-inlet)。對于出口邊界條件，一般選在離機(jī)械結(jié)構(gòu)對流場造成擾動足夠遠(yuǎn)的地方。為了得到準(zhǔn)確的結(jié)果，出口邊界必須位于最后一個障礙物后10倍于障礙高度或更遠(yuǎn)的位置。在穩(wěn)壓罐中對流動起到阻礙作用的主要障

礙物是豎隔板和3個橫隔板，穩(wěn)壓罐與出水管相連接的位置是一個突然縮小的出口管，為了更真實(shí)的模擬實(shí)流情況，在出水管20倍直徑處選擇一個垂直于流動方向的面，施加流動出口邊界條件(outlet)，同時，此面也是穩(wěn)壓罐的出口流量監(jiān)測面。

圖1 穩(wěn)壓罐的計(jì)算模型

圖2 穩(wěn)壓罐的網(wǎng)格劃分

2.3 流體模型和邊界條件的設(shè)置

穩(wěn)壓罐入口最大流速為7 m/s，根據(jù)雷諾數(shù)方程計(jì)算得Re=48 935。雷諾數(shù)48 935>8 000，因此選擇湍流模型。通過比較發(fā)現(xiàn)，RNG k-ε模型的數(shù)值模擬結(jié)果更接近，而標(biāo)準(zhǔn)k-ε模型的偏差較大且更難收斂。所以，選擇RNG k-ε模型作為流場仿真的湍流模型，近壁區(qū)域采用壁面函數(shù)法進(jìn)行處理[13]。

穩(wěn)壓罐內(nèi)上部是壓縮空氣，下部是水，氣液分界面為自由界面，這種自由層流動的形式適合采用VOF模型。在運(yùn)算環(huán)境設(shè)置中，打開Implicit Body Force，部分平衡壓力梯度和動量方程中的體積力，提高解的收斂性。穩(wěn)壓容器利用空氣的壓縮性穩(wěn)定壓力，因此將空氣設(shè)置為理想氣體，默認(rèn)打開能量方程。為了提高解的穩(wěn)定性，將可壓縮的空氣設(shè)置為主相，不可壓縮的水設(shè)置為次相。

在操作環(huán)境設(shè)置中，把參考壓力點(diǎn)設(shè)置在壓縮空氣一側(cè)，提高解的收斂速度和穩(wěn)定性，水在穩(wěn)壓容器中的流動要考慮重力的影響，在y軸的負(fù)方向上設(shè)置重力加速度。在對邊界條件設(shè)置時，需要涉及一些流動參數(shù)，如壓力、湍流參數(shù)等。在流場數(shù)值計(jì)算中，壓力總是按相對值表示的，實(shí)際求解的壓力并不是絕對值，而是相對參考壓力場而言的，因此，進(jìn)出口的壓力可以通過操作壓力估算出來。

2.4 初始化和迭代參數(shù)的設(shè)置

使用patch命令對穩(wěn)壓罐進(jìn)行初始化，分別把穩(wěn)壓罐中的氣相部分和液相部分標(biāo)識出來，罐內(nèi)上部是氣體，下部是水，由于是非穩(wěn)態(tài)問題的計(jì)算，需要對迭代參數(shù)進(jìn)行設(shè)置，包括時間步長(Time Step Size)，時間步數(shù)(Number of Time Steps)和每步的最大迭代次數(shù)等。時間步長的選取存在一個恰當(dāng)值，過大或過小都會影響模擬結(jié)果的準(zhǔn)確性。一般要小于脈動周期的1/10，以網(wǎng)格的特征尺寸除以特征速度的值為基準(zhǔn)進(jìn)行調(diào)節(jié)。

2.5 計(jì)算結(jié)果

對于流量穩(wěn)定度可以用量化指標(biāo)流量波動系數(shù)IPi進(jìn)行描述[14]。

(5)

圖3 穩(wěn)壓罐縱剖面速度云圖

從圖3可以看出，進(jìn)口的水流速度較快，水流經(jīng)進(jìn)水管進(jìn)入穩(wěn)壓罐后撞擊到豎隔板上，流速迅速下降，其中一部分水流向下流動，并在穩(wěn)壓罐的左下方形成漩渦，大部分水流向上流動，經(jīng)過3層橫隔板后，水流速度降到最低，整個穩(wěn)壓罐內(nèi)看不出明顯的流動。在出口位置，水流速度增大，可以認(rèn)為水流是從罐內(nèi)溢流出來的。再監(jiān)測面得到的數(shù)據(jù)代入式(5)進(jìn)行計(jì)算，結(jié)果如表2所示。

表2 數(shù)值模擬結(jié)果

從表2可以看出，經(jīng)過穩(wěn)壓罐的穩(wěn)壓之后，流出穩(wěn)壓罐的水流波動符合要求。

3 流場仿真結(jié)果拓展

結(jié)合上述仿真基礎(chǔ)和參數(shù)設(shè)定，本研究亦針對不同的擋板設(shè)計(jì)[15]分別進(jìn)行了流場仿真。圖4和圖5是有無橫隔板時穩(wěn)壓罐內(nèi)的流體湍流粘度的仿真結(jié)果。由圖可知橫隔板對脈動水流有二次緩沖作用，當(dāng)穩(wěn)壓罐具有3塊橫隔板時，衰減作用變大，穩(wěn)壓罐出口處水流的湍流粘度離散成都變小，從而得到更好的穩(wěn)壓新能。

對于豎隔板位置的多種方案，由計(jì)算機(jī)數(shù)值模擬得出豎隔板處于D/3處時，波動系數(shù)最小，即穩(wěn)壓效果最優(yōu)。圖6和圖7是Fluent軟件對不同方案進(jìn)行仿真得到的水粒子流動跡線圖，由圖可知，當(dāng)豎隔板位置過近或過遠(yuǎn)時，其對水流的緩沖作用也相應(yīng)過強(qiáng)或者過弱，印證了數(shù)值模擬的結(jié)果。

圖4 3塊橫隔板湍流粘度模擬

圖5 無橫隔板湍流粘度模擬

圖6 板在距離進(jìn)水端面D/6處

圖7 隔板在距離進(jìn)水端面D/3處

4 結(jié)束語

場的數(shù)值模擬知識進(jìn)行了簡單的概述，針對某流量標(biāo)準(zhǔn)裝裝置的穩(wěn)壓罐，利用Fluent軟件，采用RNG k-ε湍流模型和VOF氣液兩相流模型對其脈動水流進(jìn)行了非穩(wěn)態(tài)數(shù)值模擬計(jì)算，為了得到準(zhǔn)確的模擬效果，對入口編寫了正弦流量脈動信號，經(jīng)過迭代計(jì)算，得到了穩(wěn)壓罐監(jiān)測面的流量波動系數(shù)，結(jié)果表明，所設(shè)計(jì)的穩(wěn)壓罐在性能上符合要求。然后分析了穩(wěn)壓罐罐內(nèi)隔板結(jié)構(gòu)對穩(wěn)壓性能的影響，結(jié)果表明，在兼顧功能和經(jīng)濟(jì)性的要求下，具有3塊橫隔板的穩(wěn)壓罐結(jié)構(gòu)比較合理，同時，豎隔板位于穩(wěn)壓罐前部1/3處最佳。此經(jīng)驗(yàn)同樣適用于此類大口徑流量標(biāo)準(zhǔn)裝置的穩(wěn)壓罐設(shè)計(jì)。

[1] Lim K W, Chung M K. Numerical investigation on the installation effects of electromagnetic flowmeter downstream of a 90 elbow-laminar flow[J].Flow Measurement and Instrumentations,1999,10(3):167-174.

[2] Jaiswal S K,Yadav S,Bandyopadhyay A K,et al.Global water flow measurement and calibration facilities: review of methods and instrumentations[J].Journal of Metrology Society of India,2012,27(2):63-76.

[3] 王福軍.計(jì)算流體動力學(xué)分析-CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社,2004.

[4] 王瑞金,張凱,王剛.Fluent技術(shù)基礎(chǔ)與應(yīng)用實(shí)例[M].北京:清華大學(xué)出版社,2007.

[5] 張功暉,黎志航,周志鴻.基于Fluent流場數(shù)值仿真的管路流量計(jì)算[J].液壓氣動與密封,2010.30(12):41-43.

[6] Chou P Y,Chou R L.Fifty years of turbulence research in china[J].Annual Review of Fluid Mechanics,1995,27(1):1-16.

[7] 韓占忠,王敬,蘭小平.Fluent流體工程仿真計(jì)算實(shí)例與應(yīng)用[M].北京:北京理工大學(xué)出版社,2004.

[8] 白俊春,李培咸,郝躍,等.GaN-MOCVD系統(tǒng)反應(yīng)室流場的數(shù)值仿真[J].電子科技,2009,22(5):62-65.

[9] Guo Lanlan.Numerical study of flow fluctuation attenuation performance of a surge tank[J].Journal of Hydrodynamics,2013,25(6):938-943.

[10] 任志安,郝點(diǎn),謝紅杰.幾種湍流模型及其在Fluent中的應(yīng)用[J].化工裝備技術(shù),2009,30(2):38-40.

[11] 樊亞明,翁國華,岳堅(jiān).基于Fluent的流場分析在穩(wěn)壓罐設(shè)計(jì)中的作用[J].電子科技,2015,28(11):100-103.

[12] 王川,陸偉剛,施衛(wèi)東,等.不銹鋼沖壓潛水井泵的數(shù)值計(jì)算與試驗(yàn)驗(yàn)證[J].江蘇大學(xué)學(xué)報,2012,33(2):176-180.

[13] 李崢.水流量標(biāo)準(zhǔn)裝置不確定度和流量穩(wěn)定性研究[D].天津:天津大學(xué):2009.

[14] 王建中,梁魯林,倪聰明.液體流量標(biāo)準(zhǔn)裝置穩(wěn)壓罐性能模擬[C].成都:全國流量測量學(xué)術(shù)交流會,2012.

[15] 郭蘭蘭.影響水流量標(biāo)準(zhǔn)裝置性能的關(guān)鍵問題研究[D].濟(jì)南:山東大學(xué),2014.

Research on Flow Field Features of Regulator Tank in Flow Standard Device Based on Fluent

ZHAO Yan

(School of Mechanical Engineering,University of Shanghai for Science and Technology,Shanghai 200093,China)

The flow standard device is an important equipment for detecting and calibrating the flow instruments. The design of the regulator tank and the flow field analysis are related to whether the calibration accuracy of the device can be ensured. Based on the numerical simulation of flow field, the Fluent software is used to simulate the internal flow field of the regulator tank. After the iteration calculation, the flow fluctuation coefficient at the outlet of the regulator tank is obtained. Then, effect of structure on regulating performance is analyzed. The results show that, with the balance of function and economy, the structure of the regulator tank with three transverse partitions is reasonable, and the vertical partition is the best in the front 1/3 of the regulator tank.

Fluent;flow standard device;regulator tank;flow field analysis

2017- 03- 23

上海市科委科研計(jì)劃基金(12DZ0512600)

TN06;TB937

A

1007-7820(2018)02-020-05