環(huán)保在線APP
環(huán)保在線手機站
環(huán)保在線小程序
微信公眾號
淮安金鼎自動化儀表有限公司作者
金屬管浮子流量計是一種傳統(tǒng)的變截面流量計，具有結(jié)構(gòu)簡單、工作可靠、壓力損失小且穩(wěn)定、可測低速流體介質(zhì)等諸多優(yōu)點，廣泛應用于測量高溫、高壓及腐蝕性流體介質(zhì)。目前，金屬管浮子流量計設計采用經(jīng)典浮子流量公式計算，其中流量系數(shù)α受到浮子形狀和來流雷諾數(shù)等多種因素的影響，需要實驗標定來確定。該設計方法需進行大量的實驗，費用昂貴，周期很長，且獲取的數(shù)據(jù)有限。為深入研究金屬管浮子流量計的工作機理，筆者利用計算流體力學方法對金屬管浮子流量計進行數(shù)值仿真，通過構(gòu)造仿真模型為研究金屬管浮子流量計的機理建立數(shù)值實驗平臺，優(yōu)化傳感器的結(jié)構(gòu)；同時為金屬管浮子流量計提供了低成本、短周期的設計方法。
1 基本原理
1.1 金屬管浮子流量計的工作原理
圖1為金屬管浮子流量計的工作原理示意。在垂直的錐形管中，放置一阻力件，即浮子。當流體自下而上穿過錐管時，受到浮子迎流體積阻擋而產(chǎn)生一個壓差，從而對浮子形成向上的作用力，同時由于流體本身的黏性，對浮子產(chǎn)生黏性力，當這兩個力的合力大于浮子本身的重力時，浮子就會向上升，同時浮子與錐形管問的環(huán)形流通面積增大，流速減低，此時浮子對流體阻力作用減小.當浮子受到的力達到平衡時，浮子就會停留在某一高度。傳統(tǒng)的金屬管浮子流量計設計采用經(jīng)典浮子流量公式計算，即
（1）
式中：qv為浮子流量計的體積流量；α為流量系數(shù)；D0為浮子zui大迎流面的直徑；h為浮子在錐管中的垂直位置；φ為錐形管夾角；Vf為浮子體積；ρf為浮子材料密度；ρ為流體密度；Af為浮子垂直于流向的zui大截面面積。
1.2 計算流體力學及其控制方程
計算流體力學（computational fluid dynamics，CFD）是一門利用計算機求解描述流體流動規(guī)律的控制方程組技術，涉及到流體力學、計算方法和計算機圖形處理等技術。為了簡便，以不可壓縮湍流流動為例寫出筆者使用的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ÿ為平均應變率張量。
筆者的仿真模型建立在已有仿真實驗的基礎上，并針對以往模型中存在的入口速度剖面為等值面這一不合理設置，給出合理的圓管入口速度剖面；同時按照實際流量計的構(gòu)造，在流量計的入口與出口處加入導流架，從而提高了仿真結(jié)果的準確性。
2 模型建立和數(shù)值計算
針對25mm口徑（DN25）浮子流量計內(nèi)部機械結(jié)構(gòu)進行流場計算和分析，金屬管浮子流量計結(jié)構(gòu)剖面如圖2（a）所示，金屬管浮子流量計采用1Cr18Ni9Ti材質(zhì)。流量計設計流量范圍為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 入口速度剖面設定
入口速度由經(jīng)驗性光滑圓管湍流速度分布指數(shù)公式表示，即
（7）
式中：ur為入口截面半徑為r處的速度；R為圓管入口處半徑；m為常數(shù)，由雷諾數(shù)的大小決定，見表1；umax為圓管入口中心處的zui大速度。
（8）
式中：q為入口流量；A為金屬管入口的計算面積，對本研究，A=2πR2m2/[（2m+1）（m+1）]。
表1 常數(shù)m與雷諾數(shù)的關系
| 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，通過修正入口流量來適應多種口徑金屬管浮子流量的仿真模型。
流量的修正公式如下：
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 其他邊界條件的設定
1）流量計壁面條件的設置
由于浮子流量計內(nèi)表面采用不銹鋼材料，根據(jù)壁面特征設定：粗糙高度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ù)流體力學基本原理，可以初步判斷出計算所得結(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ù)浮子只有在所受合力為零時才能平衡，又計算了浮子在各個高度保持平衡，即合力為零時，所對應的入口流量。表2還列出了數(shù)值計算與物理實驗所得流量的比較，并計算出仿真實驗的滿度誤差占，即
（12）
式中：qs為數(shù)值模擬流量，m3/h；qp為物理實驗流量。
由表2可知，數(shù)值計算與物理實驗所得流量的zui大滿度誤差為5.4695%，平均滿度誤差為2.4731%。利用本模型對50mm、80mm口徑金屬管浮子流量計仿真數(shù)據(jù)與物理實驗數(shù)據(jù)的比較結(jié)果也相當接近。表明筆者構(gòu)建的模型得到令人滿意的結(jié)果。
4 結(jié)語
通過計算流體力學方法對金屬管浮子流量計進行仿真，并對流場進行了定量分析，計算出浮子在金屬管中不同垂直位置的受力大小及受力平衡時所對應流量。仿真實驗結(jié)果與物理實驗數(shù)據(jù)相比，zui大滿度誤差為5.4695%，平均誤差為2.4731%，表明筆者構(gòu)建的金屬管浮子流量計數(shù)值模型能夠滿足低黏度介質(zhì)金屬管浮子流量計設計要求，為金屬管浮子流量計傳感器結(jié)構(gòu)進一步優(yōu)化提供了可靠的數(shù)值仿真平臺。
全年征稿/資訊合作
聯(lián)系郵箱：hbzhan@vip.qq.com
- 版權(quán)與免責聲明
- 1、凡本網(wǎng)注明"來源：環(huán)保在線"的所有作品，版權(quán)均屬于環(huán)保在線，轉(zhuǎn)載請必須注明環(huán)保在線，http://www.hg1112.cn。違反者本網(wǎng)將追究相關法律責任。
- 2、企業(yè)發(fā)布的公司新聞、技術文章、資料下載等內(nèi)容，如涉及侵權(quán)、違規(guī)遭投訴的，一律由發(fā)布企業(yè)自行承擔責任，本網(wǎng)有權(quán)刪除內(nèi)容并追溯責任。
- 3、本網(wǎng)轉(zhuǎn)載并注明自其它來源的作品，目的在于傳遞更多信息，并不代表本網(wǎng)贊同其觀點或證實其內(nèi)容的真實性，不承擔此類作品侵權(quán)行為的直接責任及連帶責任。其他媒體、網(wǎng)站或個人從本網(wǎng)轉(zhuǎn)載時，必須保留本網(wǎng)注明的作品來源，并自負版權(quán)等法律責任。
- 4、如涉及作品內(nèi)容、版權(quán)等問題，請在作品發(fā)表之日起一周內(nèi)與本網(wǎng)聯(lián)系，否則視為放棄相關權(quán)利。
環(huán)保在線官微
環(huán)保在線APP
環(huán)保在線小程序