插入式電磁流量計的理論研究
插入式電磁流量計的理論研究
摘.. 要.. 從電磁流量計理論出發, 建立插入式電磁流量計的物理模型, 求出相應的電勢、磁勢和權函數的分布, 并進行了定量計算, 分析不同流場下插入式電磁流量計的輸出電勢差。計算研究表明, 如果把電極放置于平均流速點位置, 這種流量計可以測得通過管道的流量, 電極位置偏差所產生的測量相對誤差約為1%~ 2%。
關鍵詞.. 插入式電磁流量計; 權函數; 湍流系數; 輸出電勢差
0 .. 引言
電磁流量計是一種重要的測量導電性液體體積流量的儀表, 在城市用水、工業廢水、漿液測量及食品等多方面得到廣泛應用。但是高精度的電磁流量計價格昂貴, 特別是大管徑的, 不僅加工困難, 而且給安裝、維修帶來很多不便。因此, 在大管徑管道的流量測量方面常使用插入式電磁流量計代替傳統的管道式流量計。本文從流量計理論出發研究該插入式流量計的特性與可行性。
1 .. 電磁流量計測量理論
描述智能電磁流量計的積分式由Bevir 在1970 年給出: U2 - U1 = ..W.. .. V.. dA ( 1) 式中: U2- U1 是兩電極之間的電勢差; A 表示對所有的空間積分; W.. 稱為矢量權函數, 是一個只有電磁流量計本身結構決定的量, 其表達式為: W.. = j.. v .. B.. ( 2) 而且: j.. v= - ..G .. .. B.. = - ..F ( 3) 式中: G 和F 分別是j.. v 和B.. 的標量勢, 它們滿足Laplace 方程: ..2G = 0 .. .. ..2F = 0 ( 4) 由以上分析可知, 電勢差的測量不受流體的溫度、壓力、密度、電導率( 高于某閾值) 變化的影響, 具有很大的*性。
2 .. 插入式電磁流量計的理論計算
典型的插入式流量計結構如圖1 所示, 將電極插入管道內, 磁極留在管道外, 在電極周圍產生一個局部磁場。圖1 .. 插入式流量計結構示意圖圖2 .. 簡化的物理模型建立物理模型如圖2 所示: e1、e2 為插入管道.. 理論與實驗.的兩個電極, 電極位置由插入深度b 以及電極開角..0 決定, B.. 是由外部磁極產生的磁場。基于此模型, 計算G 、F 、W 的分布。
2..1 .. 虛電勢G 的計算
由于管道內有插入的電極, 所以不能直接使用式( 4) 的Laplace 方程求解虛電勢。我們可將該模型的虛電勢分布認為是分別只有電極和邊界產生的虛電勢的疊加, 即G = G 0+ G r 。
2..1..1 .. 只有電極的虛電勢分布假設邊界無窮遠, 根據虛電流的定義有: ..i .. d.. = e ( 5) 普通電場中的高斯定律: ..E .. d.. Q ( 6) 根據對稱原則, 虛電勢應具有和電勢相似的形式, 即: G 0= - kQ/ r ( 7) 根據圖2 的幾何關系, 不難求出G 0 的解析表達式: G 0( r , ..) = k r 2+ k 2/ cos2 ..- 2 rb / cos..sin( ..+ ..0) - k r 2+ k 2/ cos2 ..- 2 rb / cos..sin( ..- ..0)( 8)
2..1..2 .. 只有邊界的虛電勢分布由于測量管壁絕緣, 因此有邊界條件..G ..r r = R = 0, 即: ..Gr ..r r = R = - ..G 0 ..r r = R ( 9) 式( 9) 即為Gr 的邊界條件。由于此時虛電勢僅由邊界決定, 所以有方程: ..2G r = 0 ( 10) 這是一個定解條件的Laplace 方程, 使用分離變量及傅立葉系數公式可進行求解。由于很難求得邊界條件的解析解, 我們在徑向使用差分方法求得Gr 的邊界條件來求得Gr 的數值解。
3..2 .. 磁勢F 的計算
由于電極的插入深度一般僅為管道直徑的10% ~ 12 ..5% , 因此假設在電極附近的磁感強度是均勻的, 即: B.. = - y.. ( y..是y 軸方向的單位矢量) ( 11)
3..3 .. 權函數W 的計算由梯度的定義可求得: j..= - ..G = - ..G ..x x..+ ..G ..y y.. ( 12) 由于磁場是均勻的, 不難得到: W.. = B.. .. j..= ..G ..x z.. ( 13)
3..4 .. 輸出電勢差的計算
假設管道中的流動為充分發展的湍流, 我們選用經典湍流模型, 其流場分布為: v ( r ) = v max..( 1+ 2n) ( 1+ n) ( 2n2) 1- r R 1/ n ( 14) 與求得的W 在二維圓面內做數值積分即可求得輸出電勢差U。
3 .. 編程計算
綜合上述討論可以看出, 問題的關鍵在于虛電勢函數G 的計算, 考慮到精度要求以及資源消耗, 使用離散方法計算G 。具體實現步驟如下:
1) 將感興趣的區域在二維直角坐標上劃分網格, 使用式( 8) 求出每一微元上的G 0 值;
2) 使用差分方法計算式( 9) 中邊界處網格的G0 法向方向偏導值, 作為計算G r 的邊界條件;
3) 通過分離變量、利用傅立葉系數公式, 以及離散的Simphson 積分法計算式( 10) 得到Gr 的半解析表達式, 計算每一網格的G r 值, 并合成G ;
4) 按照式( 13) 計算G 在x 方向的差分, 求得每一網格的W 值;
5) 結合式( 14) 的流場模型, 計算輸出電壓。編寫程序計算不同流場, 不同電極位置的輸出電壓, 并繪制G 、W 的等勢分布圖。
4 .. 結果與分析
4..1 .. 虛電勢G 分布( 取電極間距為0..1R ) 取b= 0..9R ( R 為管道半徑) , ..= 0..0555rad, 繪制G 分布并放大電極附近區域如圖3 所示。圖3 中的黑點為電極, 可以明顯的看出G 主要分布在電極周圍并且在邊界處分布發生顯著的變化。
4..2 .. 權函數W 分布( 取電極間距為0..1R ) 取b= 0..9 R , ..= 0..0555rad, 繪制W 分布如圖.. 理論與實驗.. 4 所示。圖3.. b= 0.. 9R 時的G 分布及局部放大圖圖4.. b= 0.. 9R 時的W 分布及局部放大圖.. .. 從圖4 中可以看出W 主要分布在電極附近, 并且成對稱分布。
4..3 .. 輸出電勢差通過計算可以發現, 權函數W 主要分布在電極附近。選擇b= 0..752R , 對W.. .. V.. 進行全空間積分, 求得輸出電勢差U= 0..1475V( 為規一起見, 假定v max= 1m/ s, R = 1m, 電極處B = 1T) ; 對距離電極所在圓周0..05R 的環狀區域進行積分, 求得輸出電勢差U= 0..1231。因此, 對最終輸出電勢差起作用的主要是電極附近的流場。說明我們假設的磁場模型是可用的。選擇模擬計算中常用的湍流模型v ( r ) = v max.. ( 1+ 2n ) ( 1+ n ) ( 2n 2) 1- r R 1/ n 進行計算, 取v max = 1, 在不同的插入深度對于不同的湍流系數n 進行求解, 得到結果如表1 所示。表1 不同電極位置和不同湍流系數下的輸出電勢差湍流系數n 不同深度對應的輸出電勢差U( V) 0.. 6R 0.. 7R 0..752R 0..8 R 0.. 9R 6 0..157982 0.. 150387 0..147322 0.. 139948 0.. 123205 6.. 6 0..157027 0.. 150134 0..147501 0.. 140576 0.. 125018 7 0..156468 0.. 149975 0..147590 0.. 140923 0.. 126058 8.. 8 0..154517 0.. 149364 0..147817 0.. 142032 0.. 129582 10 0..153570 0.. 149038 0..147882 0.. 142521 0.. 131233 .. 理論與實驗.. 繪制湍流系數- 輸出電勢差曲線如圖5 所示。圖5 .. 湍流系數- 輸出電勢差擬合曲線對各組數據做最小二乘擬合, 計算斜率及線性度如表2 所示。表2 不同電極位置的電勢差擬合直線斜率及線性度插入深度0. 6R 0..7 R 0..752 R 0..8 R 0.. 9R 擬合斜率- 0.. 00109- 0.. 00034 0.. 00013 0..00063 0.. 00197 擬合線性度- 0.. 99394- 0.. 99761 0.. 96565 0..98883 0.. 99235 由圖5 可以看出, 取v max= 1, 即同量下, 不同的湍流系數n 對應了不同的輸出電壓。但當b= 0..752R , 也就是常說的平均流速點位置, 輸出的電勢差U 值基本不變。因此, 只要將電極插至該位置, 即可用來測量流量。為了研究插入深度偏離平均流速點所產生的測量誤差, 假設平均流速點位置的輸出電勢差為標準值, 計算得到: 插入深度與平均流速點偏差在0..1 R 范圍內, 輸出電勢與該標準值的相對誤差約為1% ~ 2% 。
5 .. 結論
本文完成了以下工作:
1) 建立了插入式電磁流量計的物理模型, 并編寫程序計算出虛電勢、權函數的數值解, 用于指導插入式電磁流量計的實際生產與運用;
2) 引入經典湍流模型, 對不同湍流系數, 不同電極位置的輸出電壓進行模擬計算, 給出關系曲線, 從理論上給出電極工作位置。希望在進一步的工作中能加工制作出插入式流量計的實物, 通過流量標定實驗來驗證理論分析結果。參考文獻[ 1] 金艷, 周泉, 傅新.. 插入式電磁流量計的研制.. 工業儀表與自動化裝置, 2006( 1) : 63- 64 [ 2] Shercliff J A.. Relat ion betw een th e velocity profile and the sensit iv-i t y of elect romagnet ic f lowmet ers.. J.. Phys.. D: Appl.. Phys, 1954 ( 25) : 817- 818 [ 3] Hem p J , Verst eeg H K.. Predict ion of elect romagnet ic f lowmet er charact erist ics.. J.. Phys.. D: ..Appl.. Phys, 1986 ( 19) : 1459 - 1476 [ 4] Zhang X Z.. A method f or solving Laplace.. s equat ion w ith mixed boundary condit ion in electromagnet ic flowmet ers.. J.. Phys.. D: Appl.. Phys, 1989 ( 22) 573- 576 [ 5] 王永成.. 數學物理方程.. 北京師范大學出版社, 2000: 72- 86--擴展閱讀：開封中儀流量儀表有限公司專業生產電磁流量計、孔板流量計、渦街流量計、文丘里流量計、v錐流量計、v型錐流量計、噴嘴流量計、插入式電磁流量計、智能電磁流量計、分體式電磁流量計、一體式電磁流量計、標準孔板流量計、標準孔板、一體化孔板流量計、標準噴嘴流量計、長徑噴嘴流量計、標準噴嘴、長徑噴嘴、插入式渦街流量計、智能渦街流量計、錐型流量計、v錐型流量計、節流裝置、節流孔板、限流孔板等流量產品，更多有關電磁流量計、孔板流量計、渦街流量計的信息請訪問開封中儀網站：
