[導(dǎo)讀] 為深入研究金屬管浮子流量計的工作機理，采用計算流體力學(xué)和流體力學(xué)中湍流模式理論，對金屬管流量計流場進行了數(shù)值仿真研究，獲得了低黏度流體介質(zhì)中浮子受力及浮子在受力平衡下的流量。數(shù)值模擬與物理實驗標定數(shù)據(jù)的對比表明，模擬計算的流量最大滿度誤差為5.4695%，平均滿度誤差為2.4731%，說明數(shù)值仿真模型能滿足金屬管浮子流量計設(shè)計的需要。
金屬管浮子流量計是一種傳統(tǒng)的變截面流量計，具有結(jié)構(gòu)簡單、工作可靠、壓力損失小且穩(wěn)定、可測低速流體介質(zhì)等諸多優(yōu)點，廣泛應(yīng)用于測量高溫、高壓及腐蝕性流體介質(zhì)。目前，金屬管浮子流量計設(shè)計采用經(jīng)典浮子流量公式計算，其中流量系數(shù)α受到浮子形狀和來流雷諾數(shù)等多種因素的影響，需要實驗標定來確定。該設(shè)計方法需進行大量的實驗，費用昂貴，周期很長，且獲取的數(shù)據(jù)有限。為深入研究金屬管浮子流量計的工作機理，筆者利用計算流體力學(xué)方法對金屬管浮子流量計進行數(shù)值仿真，通過構(gòu)造仿真模型為研究金屬管浮子流量計的機理建立數(shù)值實驗平臺，優(yōu)化傳感器的結(jié)構(gòu)；同時為金屬管浮子流量計提供了低成本、短周期的設(shè)計方法。
1 基本原理
1.1 金屬管浮子流量計的工作原理
圖1為金屬管浮子流量計的工作原理示意。在垂直的錐形管中，放置一阻力件，即浮子。當流體自下而上穿過錐管時，受到浮子迎流體積阻擋而產(chǎn)生一個壓差，從而對浮子形成向上的作用力，同時由于流體本身的黏性，對浮子產(chǎn)生黏性力，當這兩個力的合力大于浮子本身的重力時，浮子就會向上升，同時浮子與錐形管問的環(huán)形流通面積增大，流速減低，此時浮子對流體阻力作用減小.當浮子受到的力達到平衡時，浮子就會停留在某一高度。傳統(tǒng)的金屬管浮子流量計設(shè)計采用經(jīng)典浮子流量公式計算，即
（1）
式中：qv為浮子流量計的體積流量；α為流量系數(shù)；D0為浮子最大迎流面的直徑；h為浮子在錐管中的垂直位置；φ為錐形管夾角；Vf為浮子體積；ρf為浮子材料密度；ρ為流體密度；Af為浮子垂直于流向的最大截面面積。
1.2 計算流體力學(xué)及其控制方程
計算流體力學(xué)（computational fluid dynamics，CFD）是一門利用計算機求解描述流體流動規(guī)律的控制方程組技術(shù)，涉及到流體力學(xué)、計算方法和計算機圖形處理等技術(shù)。為了簡便，以不可壓縮湍流流動為例寫出筆者使用的k-ε模型通用形式的流體控制方程。在直角坐標系中，流動可由連續(xù)性方程和雷諾時均N-S方程描述，即
連續(xù)方程
（2）
雷諾方程
（3）
式中：ui為平均速度；P為平均壓力；ν和νt分別為分子黏性系數(shù)和渦黏性系數(shù)。
對高數(shù)湍流，渦黏性系數(shù)的計算公式為
νT=Cμk2/ ε （4）
式中：k、ε分別為湍動能和耗散率
Cμ為無量綱常數(shù)。
k和ε由各自的輸運方程得到，對高Re數(shù)問題有K方程
（5）
ε方程
（6）
式中：Gk為湍動能產(chǎn)生項
Dÿ為平均應(yīng)變率張量。
筆者的仿真模型建立在已有仿真實驗的基礎(chǔ)上，并針對以往模型中存在的入口速度剖面為等值面這一不合理設(shè)置，給出合理的圓管入口速度剖面；同時按照實際流量計的構(gòu)造，在流量計的入口與出口處加入導(dǎo)流架，從而提高了仿真結(jié)果的準確性。
2 模型建立和數(shù)值計算
針對25mm口徑（DN25）浮子流量計內(nèi)部機械結(jié)構(gòu)進行流場計算和分析，金屬管浮子流量計結(jié)構(gòu)剖面如圖2（a）所示，金屬管浮子流量計采用1Cr18NI9Ti材質(zhì)。流量計設(shè)計流量范圍為0.4~4.0 m3，量程比為10：1，為了檢測錐管中浮子的位置，在浮子中安裝有永久磁鐵，同時利用網(wǎng)格生成器建立了浮子流量計內(nèi)部結(jié)構(gòu)仿真模型，如圖2（b）所示，實驗所用的低黏度流體介質(zhì)是在常溫下密度為998.0 kg/m3、黏度為0.001002 kg/m.s）的水。金屬管浮子流量計內(nèi)部流場是高雷諾數(shù)完全發(fā)展湍流流動，故采用湍流模式理論中標準一模型來計算。
2.1 入口速度剖面設(shè)定
入口速度由經(jīng)驗性光滑圓管湍流速度分布指數(shù)公式表示，即
（7）
式中：ur為入口截面半徑為r處的速度；R為圓管入口處半徑；m為常數(shù)，由雷諾數(shù)的大小決定，見表1；umax為圓管入口中心處的最大速度。
（8）
式中：q為入口流量；A為金屬管入口的計算面積，對本研究，A=2πR2m2/[（2m+1）（m+1）]。
表1 常數(shù)m與雷諾數(shù)的關(guān)系
Re
m
4.0×103
6.0
2.3×104
6.6
1.1×105
7.0
1.1×106
8.8
2.0×106
10.0
3.2×106
10.0
為使構(gòu)建模型能滿足不同口徑金屬管浮子流量計的仿真結(jié)果，引入流量修正因子Cq，通過修正入口流量來適應(yīng)多種口徑金屬管浮子流量的仿真模型。
流量的修正公式如下：
qc=Cqq （9）
式中：q為實際流量；qc為修正流量；Cq為流量修正因子，對于25mm口徑的金屬管浮子流量計，Cq=1.13。
根據(jù)式（8）~式（11）編寫速度分布函數(shù)。圖3是q=4.8 m3/h時金屬管浮子流量計入口速度剖面的仿真。圖中色標顏色由冷色調(diào)到暖色調(diào)表示速度由小到大，從圖中可以清楚地看到從邊壁到中心的速度是由小到大的非線性分布。
2.2 其他邊界條件的設(shè)定
1）流量計壁面條件的設(shè)置
由于浮子流量計內(nèi)表面采用不銹鋼材料，根據(jù)壁面特征設(shè)定：粗糙高度Ks=0.04；粗糙常數(shù)Cs=1。
2）出入口的湍流參數(shù)
湍動能
（10）
耗散率
（11）
式中：uav為入口平均速度；Cμ為常數(shù)；L為金屬管浮子流量計內(nèi)部特征尺度。
3 仿真結(jié)果及分析
圖4給出了不同入口流量下，仿真計算得到的金屬管浮子流量計內(nèi)部截面上的速度偽色圖?？梢郧宄乜吹浇饘俟芨∽恿髁坑嬛辛黧w在浮子周圍以及出入口的速度分布。隨著流量的增加，浮子在管中的位置升，浮子與管道之間環(huán)隙變大，流體在管中的速度分布也隨之發(fā)生明顯變化。通過觀察偽色圖中的速度分布，并根據(jù)流體力學(xué)基本原理，可以初步判斷出計算所得結(jié)果是合理的。
根據(jù)物理實驗結(jié)果，仿真計算了浮子在錐管中某一高度處流體對它表面的壓力、黏性力和浮子重力Gf 3項的合力，仿真結(jié)果中浮子所受合力結(jié)果見表2，合力向上為正方向。
表2 浮子所受合力Ff的仿真數(shù)據(jù)以及數(shù)值計算與物理實驗所得流量的比較
浮子高度h/mm
計算合力Ff/N
實驗流量qp/(m3·h-1)
計算流量qs/(m3·h-1)
滿度誤差δs/%
4.86
0.1986
0.4
0.3853
-0.3660
9.36
0.3955
O.8
0.7746
-0.6342
14.41
-0.0790
1.2
1.2298
0.7447
19.28
-0.2751
1.6
1.7246
3.1159
23.97
-0.4258
2.0
2.2188
5.4695
28.65
-O.3175
2.4
2.5670
4.1739
32.77
-0.354 7
2.8
3.0029
5.0727
36.89
0.3755
3.2
3.3898
4.7447
40.44
-0.3435
3.6
3.8056
5.1409
44.00
-0.2853
4.0
4.2134
5.3355
47.19
-O.O518
4.4
4.4605
1.5125
50.00
0.329 2
4.8
4.6147
-4.6317
由表2可知計算所得浮子在各位置的合力接近于零，趨于平衡，為驗證仿真模型所能達到的仿真精度，根據(jù)浮子只有在所受合力為零時才能平衡，又計算了浮子在各個高度保持平衡，即合力為零時，所對應(yīng)的入口流量。表2還列出了數(shù)值計算與物理實驗所得流量的比較，并計算出仿真實驗的滿度誤差占，即
（12）
式中：qs為數(shù)值模擬流量，m3/h；qp為物理實驗流量。
由表2可知，數(shù)值計算與物理實驗所得流量的最大滿度誤差為5.4695%，平均滿度誤差為2.4731%。利用本模型對50mm、80mm口徑金屬管浮子流量計仿真數(shù)據(jù)與物理實驗數(shù)據(jù)的比較結(jié)果也相當接近。表明筆者構(gòu)建的模型得到令人滿意的結(jié)果。
4 結(jié)語
通過計算流體力學(xué)方法對金屬管浮子流量計進行仿真，并對流場進行了定量分析，計算出浮子在金屬管中不同垂直位置的受力大小及受力平衡時所對應(yīng)流量。仿真實驗結(jié)果與物理實驗數(shù)據(jù)相比，最大滿度誤差為5.4695%，平均誤差為2.4731%，表明筆者構(gòu)建的金屬管浮子流量計數(shù)值模型能夠滿足低黏度介質(zhì)金屬管浮子流量計設(shè)計要求，為金屬管浮子流量計傳感器結(jié)構(gòu)進一步優(yōu)化提供了可靠的數(shù)值仿真平臺。