金屬管浮子流量計是一種傳統的變截面流量計，具有結構簡單、工作可靠、壓力損失小且穩定、可測低速流體介質等諸多優點，廣泛應用于測量高溫、高壓及腐蝕性流體介質。目前，金屬管浮子流量計設計采用經典浮子流量公式計算，其中流量系數α受到浮子形狀和來流雷諾數等多種因素的影響，需要實驗標定來確定。該設計方法需進行大量的實驗，費用昂貴，周期很長，且獲取的數據有限。為深入研究金屬管浮子流量計的工作機理，筆者利用計算流體力學方法對金屬管浮子流量計進行數值仿真，通過構造仿真模型為研究金屬管浮子流量計的機理建立數值實驗平臺，優化傳感器的結構；同時為金屬管浮子流量計提供了低成本、短周期的設計方法。
1 基本原理
1.1 金屬管浮子流量計的工作原理
圖1為金屬管浮子流量計的工作原理示意。在垂直的錐形管中，放置一阻力件，即浮子。當流體自下而上穿過錐管時，受到浮子迎流體積阻擋而產生一個壓差，從而對浮子形成向上的作用力，同時由于流體本身的黏性，對浮子產生黏性力，當這兩個力的合力大于浮子本身的重力時，浮子就會向上升，同時浮子與錐形管問的環形流通面積增大，流速減低，此時浮子對流體阻力作用減小.當浮子受到的力達到平衡時，浮子就會停留在某一高度。傳統的金屬管浮子流量計設計采用經典浮子流量公式計算，即
（1）
式中：qv為浮子流量計的體積流量；α為流量系數；D0為浮子zui大迎流面的直徑；h為浮子在錐管中的垂直位置；φ為錐形管夾角；Vf為浮子體積；ρf為浮子材料密度；ρ為流體密度；Af為浮子垂直于流向的zui大截面面積。
1.2 計算流體力學及其控制方程
計算流體力學（computational fluid dynamics，CFD）是一門利用計算機求解描述流體流動規律的控制方程組技術，涉及到流體力學、計算方法和計算機圖形處理等技術。為了簡便，以不可壓縮湍流流動為例寫出筆者使用的k-ε模型通用形式的流體控制方程。在直角坐標系中，流動可由連續性方程和雷諾時均N-S方程描述，即
連續方程
（2）
雷諾方程
（3）
式中：ui為平均速度；P為平均壓力；ν和νt分別為分子黏性系數和渦黏性系數。
對高數湍流，渦黏性系數的計算公式為
νT=Cμk2/ ε （4）
式中：k、ε分別為湍動能和耗散率
Cμ為無量綱常數。
k和ε由各自的輸運方程得到，對高Re數問題有K方程
（5）
ε方程
（6）
式中：Gk為湍動能產生項
Dÿ為平均應變率張量。
筆者的仿真模型建立在已有仿真實驗的基礎上，并針對以往模型中存在的入口速度剖面為等值面這一不合理設置，給出合理的圓管入口速度剖面；同時按照實際流量計的構造，在流量計的入口與出口處加入導流架，從而提高了仿真結果的準確性。
2 模型建立和數值計算
針對25mm口徑（DN25）浮子流量計內部機械結構進行流場計算和分析，金屬管浮子流量計結構剖面如圖2（a）所示，金屬管浮子流量計采用1Cr18Ni9Ti材質。流量計設計流量范圍為0.4~4.0 m3，量程比為10：1，為了檢測錐管中浮子的位置，在浮子中安裝有*磁鐵，同時利用網格生成器建立了浮子流量計內部結構仿真模型，如圖2（b）所示，實驗所用的低黏度流體介質是在常溫下密度為998.0 kg/m3、黏度為0.001002 kg/m.s）的水。金屬管浮子流量計內部流場是高雷諾數*發展湍流流動，故采用湍流模式理論中標準一模型來計算。
2.1 入口速度剖面設定
入口速度由經驗性光滑圓管湍流速度分布指數公式表示，即
（7）
式中：ur為入口截面半徑為r處的速度；R為圓管入口處半徑；m為常數，由雷諾數的大小決定，見表1；umax為圓管入口中心處的zui大速度。
（8）
式中：q為入口流量；A為金屬管入口的計算面積，對本研究，A=2πR2m2/[（2m+1）（m+1）]。
表1 常數m與雷諾數的關系
| 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 |
為使構建模型能滿足不同口徑金屬管浮子流量計的仿真結果，引入流量修正因子Cq，通過修正入口流量來適應多種口徑金屬管浮子流量的仿真模型。
流量的修正公式如下：
qc=Cqq （9）
式中：q為實際流量；qc為修正流量；Cq為流量修正因子，對于25mm口徑的金屬管浮子流量計，Cq=1.13。
根據式（8）~式（11）編寫速度分布函數。圖3是q=4.8 m3/h時金屬管浮子流量計入口速度剖面的仿真。圖中色標顏色由冷色調到暖色調表示速度由小到大，從圖中可以清楚地看到從邊壁到中心的速度是由小到大的非線性分布。
2.2 其他邊界條件的設定
1）流量計壁面條件的設置
由于浮子流量計內表面采用不銹鋼材料，根據壁面特征設定：粗糙高度Ks=0.04；粗糙常數Cs=1。
2）出入口的湍流參數
湍動能
（10）
耗散率
（11）
式中：uav為入口平均速度；Cμ為常數；L為金屬管浮子流量計內部特征尺度。
3 仿真結果及分析
圖4給出了不同入口流量下，仿真計算得到的金屬管浮子流量計內部截面上的速度偽色圖。可以清楚地看到金屬管浮子流量計中流體在浮子周圍以及出入口的速度分布。隨著流量的增加，浮子在管中的位置升，浮子與管道之間環隙變大，流體在管中的速度分布也隨之發生明顯變化。通過觀察偽色圖中的速度分布，并根據流體力學基本原理，可以初步判斷出計算所得結果是合理的。
根據物理實驗結果，仿真計算了浮子在錐管中某一高度處流體對它表面的壓力、黏性力和浮子重力Gf 3項的合力，仿真結果中浮子所受合力結果見表2，合力向上為正方向。
表2 浮子所受合力Ff的仿真數據以及數值計算與物理實驗所得流量的比較
| 浮子高度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可知計算所得浮子在各位置的合力接近于零，趨于平衡，為驗證仿真模型所能達到的仿真精度，根據浮子只有在所受合力為零時才能平衡，又計算了浮子在各個高度保持平衡，即合力為零時，所對應的入口流量。表2還列出了數值計算與物理實驗所得流量的比較，并計算出仿真實驗的滿度誤差占，即
（12）
式中：qs為數值模擬流量，m3/h；qp為物理實驗流量。
由表2可知，數值計算與物理實驗所得流量的zui大滿度誤差為5.4695%，平均滿度誤差為2.4731%。利用本模型對50mm、80mm口徑金屬管浮子流量計仿真數據與物理實驗數據的比較結果也相當接近。表明筆者構建的模型得到令人滿意的結果。
4 結語
通過計算流體力學方法對金屬管浮子流量計進行仿真，并對流場進行了定量分析，計算出浮子在金屬管中不同垂直位置的受力大小及受力平衡時所對應流量。仿真實驗結果與物理實驗數據相比，zui大滿度誤差為5.4695%，平均誤差為2.4731%，表明筆者構建的金屬管浮子流量計數值模型能夠滿足低黏度介質金屬管浮子流量計設計要求，為金屬管浮子流量計傳感器結構進一步優化提供了可靠的數值仿真平臺。
