低液位非滿管流量校準(zhǔn)及其測量影響因素分析

劉 瑛， 陳 旭， 徐明剛， 陳銀銀

(北方工業(yè)大學(xué) 機械與材料工程學(xué)院，北京 100144)

非滿管流量測量在市政排污、農(nóng)業(yè)灌溉、廢水處理等領(lǐng)域中被廣泛應(yīng)用。非滿管指液體未充滿管道的狀態(tài)，管內(nèi)為氣液混合物。在實際情況中，液體流經(jīng)管道時，很難不留空氣的充滿整個管道，同時管道的管壁厚度、截面尺寸、電極之間間距等因素的變化，都很可能引起管內(nèi)液體的液位變化[1]。

非滿管電磁流量計最早于20世紀(jì)90年代被Fisher & Porter公司研制成功[2]。它的出現(xiàn)，使管道橫截面為圓形情況下的非滿管流量測量精度得到提高[3]。目前，非滿管流量測量主要通過德國科隆公司的TDALFLUX、ABB公司的Pati-MagII、上海大學(xué)的非滿管流量測量系統(tǒng)[4～7]進行流量測量。流量測量的主要參數(shù)是管道內(nèi)流體的液位高度和平均流速，文獻(xiàn)[4]～文獻(xiàn)[7]所述流量測量方法均是通過測量主要參數(shù)進行流量測量的，但其電極必須位于被測流體中，流體液位高度為較低充滿度時(管道直徑的10%)，電極無法全部浸入到被測流體中，難以實現(xiàn)對信號的采集[7-9]。

為解決上述問題，本文提出了一種低液位非滿管流量測量方法：基于連通管原理結(jié)合高精度激光測距儀測得流體液位高度值，管道上方通氣管結(jié)合流速儀精確測量管道內(nèi)流體平均流速，通過體積流量計算公式得到流量值。

影響流速的因素有流體黏度和管道粗糙度，文獻(xiàn)[4]～文獻(xiàn)[7]所述流量測量方法均未對流體黏度和管道粗糙度進行分析。晏飛等[10]基于Herschel-Bulkley模型，運用CFD對潤滑脂的流動進行數(shù)值模擬，得到黏度和速度的關(guān)系分布規(guī)律。Santos[11]通過玻爾茲曼方程和Bhatnagar-Gross-Krook動力學(xué)模型對簡單縱向流動中的黏度和速度進行了分析，得到黏度和速度的非線性函數(shù)。張桂欣等[12]、高楊等[13]分別應(yīng)用激光測距儀和三維超聲波流速儀，得到粗糙度和速度的分布規(guī)律。

本文運用FLUENT仿真軟件對管道粗糙度和流體黏度進行模擬，建立粗糙度-平均流速、黏度-平均流速模型，分析出粗糙度、黏度與流量的具體關(guān)系，通過多次實驗驗證模型的準(zhǔn)確性及測量方法的可行性。預(yù)計本文提出的測量方法適用于流體液位高度為較低充滿度時(管道直徑的10%)的流量測量，對低液位非滿管流量測量的發(fā)展與創(chuàng)新具有重要意義。

1 流量測量原理

管道中間部分平均流速的測量：將流速傳感器經(jīng)管道中間上方通氣管放入待測流體中，旋槳位于當(dāng)前液位高度的50%～60%，進行平均流速的測量。測量原理為

(1)

測量平均流速時，槳葉在水力推動作用下旋轉(zhuǎn)，內(nèi)置信號裝置產(chǎn)生轉(zhuǎn)數(shù)信號，槳葉水力螺距K和流速儀常數(shù)C均為常數(shù)，測流歷時T和該時段內(nèi)產(chǎn)生的信號數(shù)N測得后，即可算出測流時段內(nèi)的平均流速值。

管道內(nèi)液位高度的測量：將激光測距儀垂直置于U型連通管正上方，連通管中置一浮漂，進行液位高度的測量。

根據(jù)連通管原理可得

ρ1·g·h1=ρ2·g·h2

(2)

式中，所測流體相同、密度相等(ρ1=ρ2)、兩邊壓強一致時，管內(nèi)液位高度與U型管內(nèi)液位高度相等?；诖嗽?，可在不干擾管道內(nèi)部流體流動的情況下，通過測量U型管內(nèi)液位高度值，實現(xiàn)管道內(nèi)液位高度的測量。同時，外置型的安裝使激光測距儀的更換更加簡便，降低了維修成本。

2 仿真分析

2.1 流量計算

根據(jù)體積流量公式可得

(3)

(4)

式中，D為管道內(nèi)徑，單位為m；h為液位高度，單位為m。

2.2 建立模型

在SolidWorks中建立進水管直徑18 mm，管道直徑100 mm，長3 m的圓形截面直管道，前1.5 m為湍流發(fā)展前置段。在Workbench中將建立好的模型導(dǎo)入Geometry中進行優(yōu)化后，將模型導(dǎo)入Mesh中進行網(wǎng)格劃分，劃分方式為四面體網(wǎng)格，網(wǎng)格總數(shù)為1870715，管道中間部分網(wǎng)格劃分如圖1所示。對網(wǎng)格質(zhì)量進行檢查，最小體積為正值且最小正交質(zhì)量為0.23，大于0.2，符合仿真網(wǎng)格質(zhì)量要求，F(xiàn)LUENT中求解器選取分離式，采用SIMPLE算法，采用隱式算法。

圖1 局部網(wǎng)格劃分圖

由于水、植物油和空氣多相物質(zhì)同時存在，因此選擇多相流(Multiphase)下的VOF模型進行計算，流體介質(zhì)包括空氣、水、植物油，各介質(zhì)基本參數(shù)如表1所示。

表1 介質(zhì)基本參數(shù)取值表

設(shè)置進水口為速度進口，速度大小為0.4 m/s，設(shè)置湍流強度為5%，水力直徑取入水管直徑18 mm[15～16]。通氣管為空氣進口，U型管、水平管出口設(shè)置為壓力出口，管道主體為固體壁面，采用無滑移邊界條件，設(shè)置壁面粗糙度分別為0 mm，0.033 mm，0.049 mm，0.139 mm，0.150 mm，0.445 mm。在計算過程中，當(dāng)管道內(nèi)部流動形態(tài)變化較小時，認(rèn)為流動趨于穩(wěn)定，以此時計算結(jié)果作為最終流場進行分析。

2.3 仿真結(jié)果分析

對6種管道粗糙度及3種流體黏度進行仿真分析，管道1.4～1.6 m處的速度分布數(shù)據(jù)通過Origin進行擬合處理，得到粗糙度-平均流速模型和黏度-平均流速模型。粗糙度-平均流速模型見式(5)，擬合曲線如圖2所示。

(5)

式中，Ra為粗糙度，單位為mm；A1，A2，A3為常數(shù)項，具體數(shù)值如表2所示。

表2 粗糙度-平均流速模型常數(shù)取值表

對粗糙度-平均流速模型進行分析，可得到流量與粗糙度的具體關(guān)系，見式(6)。

(6)

黏度-平均流速模型見式(7)，擬合曲線如圖3所示。

(7)

式中，μ為流體黏度，單位為kg/(m·s)；B1，B2，B3為常數(shù)項，具體數(shù)值如表3所示。

圖2 不同流體的粗糙度-平均速度曲線圖

圖3 不同粗糙度的黏度-平均速度曲線圖

流體黏度μ與雷諾數(shù)Re的關(guān)系見式(8)。

Re=ρvd/μ

(8)

式中，Re為雷諾系數(shù)；ρ為流體密度，單位為kg/m3；v為流體特征速度，單位為m/s；d為特征長度，即流體液位高度h，單位為m。對式(7)、式(8)進行分析，可得到雷諾系數(shù)與平均流速的關(guān)系：

(9)

對黏度-平均流速模進行分析，可得到流量與黏度的具體關(guān)系：

(10)

表3 黏度-平均流速模型常數(shù)取值表

實際應(yīng)用中，需增加數(shù)據(jù)查表功能。數(shù)據(jù)表存儲不同介質(zhì)及對應(yīng)參數(shù)值，根據(jù)不同介質(zhì)，選擇對應(yīng)參數(shù)。

3 實驗與分析

3.1 實驗平臺

為驗證測量方法的可行性，搭建實驗平臺，設(shè)計方案示意圖如圖4所示，實物圖如圖5所示。管道外置U型連通管，激光測距儀沿垂直方向朝下發(fā)射激光，激光終點為浮漂上表面，以此實現(xiàn)液位高度的測量。管道中間正上方焊接一個通氣管，將流速儀垂直放入待測流體中進行測量，旋槳位于當(dāng)前液位高度的50%～60%，液位不得低于旋槳直徑的一半(6 mm)，測量信號經(jīng)流量計算模型處理后得到對應(yīng)流量值。

圖4 實驗平臺設(shè)計方案示意圖

圖5 實驗平臺實物圖

圖4中實驗器材如下。

① 自制0.5 m×0.5 m×0.5 m的不銹鋼連通水箱；

同樣是金枝玉葉的段譽，第一次來燕子塢吃的那些：“茭白蝦仁”“龍井茶葉雞丁”，看看就教人饞涎欲滴。段譽的當(dāng)時心理評判是這樣的：“魚蝦肉食之中混以花瓣鮮果，色彩既美，自別有天然清香?！?/p>

② 額定功率為250 W的水泵；

③ L20的彎頭、溢流閥；

④ QS智能型電動調(diào)節(jié)閥；

⑤ 精度等級為±0.2%FS的DN20電磁流量計；

⑥ 精度等級為±0.01%FS的激光測距傳感器；

⑦ 定制的U型連通管；

⑧ 長3 m、粗糙度0.15 mm的DN100標(biāo)準(zhǔn)鍍鋅管道以及配對法蘭盤；

⑨ 精度等級為±1.0%FS的便攜式流速傳感器；

⑩ ARM公司STM32F103系列的微控制器信號處理單元。

3.2 實驗測量

在μC/OS II操作系統(tǒng)下對信號處理單元進行編程，激光測距傳感器的采集頻率為5 Hz，每1 h采集一次流速傳感器的測量值，信號處理單元將采集到的液位測量值進行均值處理，連續(xù)測量1 h后，經(jīng)算法處理輸出該時段的累計流量值，同時根據(jù)Modbus標(biāo)準(zhǔn)協(xié)議由串口向DN20電磁流量計發(fā)送指令得到1 h實際流過的流量值。

最后通過串口將計算值和電磁流量計的標(biāo)準(zhǔn)流量值發(fā)送到PC端，實驗記錄結(jié)果如表4所示。

表4 實驗結(jié)果記錄表

對電磁調(diào)節(jié)閥進行調(diào)節(jié)，使DN20電磁流量計流速顯示0.4 m/s，與仿真進水口流速一致，進行測量實驗。通過實驗發(fā)現(xiàn)，調(diào)節(jié)閥開合度為28%時，DN20電磁流量計流速顯示0.4 m/s。

表5 仿真結(jié)果與實驗結(jié)果對比表

4 結(jié)束語

本文通過FLUENT仿真軟件對低液位非滿管流量測量影響因素進行分析，基于提出的測量方法進行實驗，實驗結(jié)果表明：

① 管道內(nèi)液位高度不大于0.1D時，流量測量誤差不大于0.85%，通過粗糙度-平均流速、黏度-平均流速模型計算得到的流量值，誤差不大于3.03%；

② 實驗驗證了該方法的可行性和粗糙度-平均流速、黏度-平均流速模型的可靠性，適用于低液位的非滿管流量計量。

本文創(chuàng)新點如下：

① 提出一種低液位非滿管流量測量方法，實現(xiàn)對低液位流量的精確測量；

② 通過建立粗糙度-平均流速、黏度-平均流速模型，進一步得到粗糙度、黏度與流量的具體關(guān)系。

論文有待進一步解決的問題如下：

① 受流速傳感器本身結(jié)構(gòu)的影響無法對含大顆粒的流體進行測量，在對結(jié)構(gòu)的優(yōu)化方面有待進一步解決；

② 受實驗室環(huán)境限制，無法對大口徑管道進行模擬和實驗，大口徑管道的非滿管流量測量有待進一步研究。