【摘要】通過ICEMCFD軟件將局部截面變?yōu)闄E圓形的異徑導流筒進行了三維建模，使用Fluent對不同入口速度下的流線場與速度分布進行仿真計算，建立了不同結構的導流筒所適用的速度范圍.結果表明，速度的大小和橢圓截面離心率對流場產(chǎn)生的影響較大.當速度減小或離心率變大時，導流筒尾部漸擴管容易發(fā)生回流，致使流場紊亂.本研究能為橢圓形管道電磁流量計的結構設計提供參考方案，為設計合理的導流筒提供理論依據(jù).
　　電磁流量計是工業(yè)過程中用于計量導電性流體體積流量的儀表［1］，當前國內使用大多電磁流量計為圓形截面導流筒.然而，電磁流量計對被測管道內的流場有一定的要求，流場的不穩(wěn)定會使得流量計示值不穩(wěn)定，致使測量誤差加大［2－3］.為了解決這些問題，本文提出橢圓形截面管道設計方案.
　　針對橫截面為不同離心率橢圓形的導流筒，對在不同入口速度下流場的流動性與速度分布進行Fluent仿真，欲為合理的導流筒提供理論依據(jù).
1異徑管結構的電磁理論分析
　　電磁流量計是基于法拉第電磁感應定律而開發(fā)的計量儀表［8］.通電后的勵磁線圈在導流筒垂直方位產(chǎn)生磁感應強度為Ｂ的工作磁場，待導電流體穿過時，在液體兩側產(chǎn)生感應電動勢Ｅ，通過對相應的電動勢進行信號處理而實現(xiàn)體積流量的準確測量.感應電動勢大小為
Ｅ＝ＢＶＤ.（1）
　　式（1）中：Ｂ為工作磁場中的磁感應強度；Ｖ為導電液體流速；Ｄ為測量導管內徑.
　　導電流體的速度Ｖ與工作磁場內的磁感應強度Ｂ都是有方向性的矢量，但各質點的速度為非均勻分布，當流體的流速很小時，會產(chǎn)生很小感應電動勢，與噪音混合后使得測量誤差增大，從而影響到設備的穩(wěn)定性和可靠性.其中Ｅ的數(shù)值由電極測量，單位時間內管道流量計算公式為

　　在電磁流量計的勵磁線圈中，電流為Ｉ，匝數(shù)為Ｎ，穿過工作區(qū)域的磁路長度均值為Ｌ，可得磁阻Ｒm與磁通勢Ｆ為

　　式中Ｓ為磁路的平均面積，μ為介質磁導率.由磁場歐姆定律［9］可得磁通量f

　　由（6）式可知，磁感應強度Ｂ與磁路長度平均值Ｌ成反比，與通過勵磁線圈的電流Ｉ成正比.相比起均勻的圓形管道，橢圓導流筒內的工作磁場縮小了Ｌ值，在產(chǎn)生同等磁感應強度Ｂ的條件下，勵磁線圈中的電流將小于前者，從而可降低電磁流量計的功耗.
2Fluent模型建立與參數(shù)設置
　　使用ICEMCFD建立橢圓截面導流筒的模型.導流筒的中間部分為橢圓管，兩側均為橢圓形漸變?yōu)閳A形的漸擴管.導流筒半長軸與Ｘ軸平行，長度35mm，半短軸與Ｙ軸平行，長度28mm，短長半軸之比為4/5，橢圓離心率為0.60，長88mm.兩端漸擴管最外側圓形的半徑為50mm，各長81mm.導流筒總長250mm.該模型的對象為在中間直管段具有均勻磁場分布的橢圓截面管道的電磁流量計將導流筒兩端分別定義為出口與出口.流體在入口邊界以固定速度垂直與入口邊界流入，在出口邊界自由流出，忽略重力.定義其他區(qū)域為壁面，最后以四面體結構對模型進行網(wǎng)格劃分，如圖1所示.單元格數(shù)量為204萬，網(wǎng)格質量評價系數(shù)為：0.65～0.70（2.5%）；0.70～0.90（8.6%）；0.90～1.0（86.2%）.該三維模型網(wǎng)格質量能夠滿足精度和收斂要求.文中其它結構的三維模型網(wǎng)格，其類型與上述一致，網(wǎng)格質量基本相同.

　　設置模型為ｋ－ｅｐｓｉｌｏｎ湍流模型［10］，模擬對象為液體水，仿真將以入口流速分別為小流速0.1m/ｓ、0.3m/ｓ與大流速5.0m/ｓ的條件下進行.
3速度場仿真結果分析
　　以不同進口速度對該結構導流筒進行流場仿真，求解后使用軟件提取數(shù)據(jù).由于磁場方向平行于Ｙ軸，故圖2至圖15是在選取了與Ｙ軸垂直的ＸＯＺ坐標平面，并觀察速度云與流線分布圖，計算結果如下.
3.1小流速下的仿真分析
　　取流入速度為0.1m/ｓ、0.3m/ｓ，設置仿真計算的迭代步數(shù)為300，過程中分別在第211步、第186步時計算結果收斂，流量計流道區(qū)域內可視為穩(wěn)態(tài)的定常流動.管內速度云圖如圖2、圖3，流線圖如圖4、圖5.


　　由圖2、圖3可知，在進口速度為0.1m/ｓ與0.3m/ｓ條件下，速度云圖無明顯差別，平面直管段的速度分布的上下對稱性較高，靠管壁速度小，中間大，出口流體向兩側流動，中間區(qū)流速小.

　　如圖4、圖5，當入口速度為0.1m/ｓ時，末端發(fā)生回流現(xiàn)象，但中間直管段流場平穩(wěn)，沒有受到尾部回流影響.當初始速度增加為0.3m/ｓ時尾部的回流減弱.
3.2大流速下的仿真分析
　　設置進口速度為5.0m/ｓ，設置仿真計算的迭代步數(shù)為300，過程中在第96步計算結果受斂，可視為定常流動.速度云圖如圖6.

　　中間直管段內靠管壁處速度小，中間大，速度分布的上下對稱性較高.在圖7中，當流速增加為5.0m/ｓ時，中間直管段與尾部漸擴管的流場非常平穩(wěn)，無回流現(xiàn)象.

　　綜合圖4、圖5、圖7可見，隨著流體速度增加，回流減弱.綜合3.1與3.2，流道域內均為穩(wěn)態(tài)的定常流動，且流場平穩(wěn)，速度分布對稱性較高，故該結構的電磁流量計在大小流速條件下的使用均是可行的.
4離心率對流場的影響
4.1離心率為0.8
　　中間橢圓截面直管段短長半軸之比為3/5，離心率0.8.分別定義入口速度為0.1m/ｓ、5.0m/ｓ，在此條件下使用Fluent進行模擬計算，過程中分別在第263步、192步時計算結果收斂，可視為定常流動.結果如圖8至圖10.


　　入口速度為0.1m/ｓ時（圖8、圖9），中間直管段內靠近但不接觸管壁的位置流速大，中間小.速度分布的上下對稱性較高，流道域尾部出現(xiàn)回流現(xiàn)象，但中間直管端的流場依然平穩(wěn).當入口速度增加至5.0m/ｓ時（圖10、圖11），中間直管段內速度分布基本均勻，尾部回流消失，流場整體平穩(wěn).

　　縮徑為0.8離心率的橢圓截面電磁流量計在初始流速為0.1m/ｓ與5.0m/ｓ條件下均為穩(wěn)態(tài)流動，速度分布對稱，直管內流場平穩(wěn)，那么該結構導流筒的電磁流量計在大小流速條件下的使用均是可行的.
4.2離心率為0.916
　　半長軸長35mm，半短軸長14mm，短長半軸之比2/5，離心率0.916.分別設置入口速度在0.1m/ｓ、5.0m/ｓ的條件下通過Fluent進行模擬仿真，設置計算迭代步數(shù)為1000，過程中各點的速度值隨時間產(chǎn)生無規(guī)律變化，無法收斂.圖11至圖14為步數(shù)等于1000時瞬時結果的抓取.

　　由圖12、圖13可知，當入口速度為0.1m/ｓ時，導流筒內速度分布無明顯規(guī)律，存在較大的流場畸變.因為導流筒兩側產(chǎn)生的感應電動勢與流速成正比，且流量計是根據(jù)流速值計算出一定時間內通過管道的體積流量，所以在非穩(wěn)態(tài)流場條件下流量計檢測到的是大小搖擺不定的感應電動勢，為體積流量的計算造成許多不確定因素，還降低了計量精度.
　　設置流入速度為5.0m/ｓ，計算過程中第117步收斂，流場可視為達到穩(wěn)定狀態(tài).如圖14、圖15所示，流場分布平穩(wěn)，中間直管段內速度場分布基本均勻，與其它結構導流筒在該速度下的分布無明顯區(qū)別.綜合圖12至圖15可知，截面離心率變?yōu)?.916時的導流筒在入口速度增大到一定值后，流場穩(wěn)定.

5不同結構導流筒所適應的速度區(qū)間
　　在完成不同結構導流筒在小流速與大流速情況下的仿真之后，對入口流速分別為0.03m/ｓ、0.5m/ｓ、0.8m/ｓ、1m/ｓ、3m/ｓ的條件下進行模擬計算.以流場速度分布為判據(jù)，得出了不同結構橢圓管所適應的速度區(qū)間.由表1可知：截面離心率為0.600和0.800的橢圓形導流筒的速度均適用于大流速與小流速，而截面離心率為0.916的導流筒卻不適用于小流量的條件，當該結構導流筒的入口流速達到0.8m/ｓ及以上時，內部流場分布才被接受.雖橢圓變扁，磁路長度平均值Ｌ減小，致使所需勵磁電流Ｉ減小，降低了設備功耗，但縮徑量過大會犧牲測速量程，導致量程下限升高，小流量的狀態(tài)下不再適用.

6結論
本文針對局部變?yōu)闄E圓形截面的異徑導流筒進行了模擬仿真計算.得出結論如下：
1）減小磁路長度平均值Ｌ，在產(chǎn)生同等磁感應強度Ｂ的條件下，可減小勵磁線圈的電流Ｉ，從而提升流量計的靈敏度，降低功耗.
2）當橢圓離心率增大到一定值時，尾部漸擴管便會出現(xiàn)明顯的回流現(xiàn)象，致使流量計量程下限升高，不再適用于低速計量.
3）入口速度對管內速度場的影響頗為重要，大流速在導流筒各部位的流場較平穩(wěn)，小流速則容易發(fā)生回流現(xiàn)象，隨著入口速度降低，回流更顯著.
4）離心率為0.8的橢圓截面導流筒可最大條件下滿足縮徑和流場要求，該尺寸適合在流量計中使用.

以上內容來源于網(wǎng)絡，如有侵權請聯(lián)系即刪除!