時差法超聲波流量計流速修正系數(shù)的數(shù)值模擬.docx

1、時差法超聲波流量計流速修正系數(shù)的數(shù)值模擬吳志敏，蘇滿紅，葉瑋淵（深圳職業(yè)技術(shù)學(xué)院，廣東深圳518055）摘要文章利用了計算機數(shù)值模擬的方法得出了超聲波流量計流速測雖的修正系數(shù)，這種數(shù)值模擬的方法是基于計算流體力學(xué)（CFD）軟件Fluent±實現(xiàn)的。通過計算在不同的超聲波波束直徑和不同常諾系數(shù)下流速測最的修正系數(shù),與利用數(shù)學(xué)分析計算方法得出的結(jié)果進行了相應(yīng)的比較，結(jié)果表明，兩種方法得出的修正系數(shù)具有很好的一致性。因此，對于在更為復(fù)雜的測故環(huán)境中，流速的修正系數(shù)用數(shù)值模擬的方法具有很好的輔助作用。關(guān)鍵詞時差法超聲波流域；數(shù)值模擬；CFD；計算分析中圖分類號TB551文獻標識碼B文章編號1

2、002-1183（2008）04-0001-04NumericalsimulationofflowvelocitycorrectionfactorforultrasonicflowmeterwithtransittimemethodWUZhimin,SUMan-hong,YEWeiyuan(ShenzhenPolytechnicColle取，Shenzhen518055,China)Abstract:Correctionfactorofflowvelocitymeasurementinultrasonicflowmeterwasobtainedbyusingcomputernumericals

3、imulationapproach.Theproposedtechniqueisbasedoncomputationalfluiddynamics(CFD)softwarefluent.FlowvelocitycorrectionfactorswithdifferentultrasonicbeamdimensionsandReynoldsnumberswerecalculatedbyanalyticalcalculationmethod.TheywerecomparedwithcorrectionfactorsobtainedbyFluentsimulation.Resultsshowthat

4、correctionfactorsobtainedbytwokindsofmethodhaveagoodconsistency.So,incomplicatedflowmetermeasurementenvironment,flowvelocitycorrectionfactorobtainedbynumericalsimulationmethodhasagoodeffectKeywords：transittimeultrasonicflowmeter;numericalSimulation；CFD；analyticalcalculation收稿日期2008-02-26-作者簡介吳志敏（196

5、7-）,男，河北石家莊人，博上研究生，副教授，畢業(yè)于西安交通大學(xué)，研究方向為智能儀表與機電-體化的研究與開發(fā)。工業(yè)計量2008年第18卷第4期】時差法超聲波流量計測量流量是目前廣泛應(yīng)用的一種方法，它具有安裝簡單、使用方便、測量范圍寬等特點。近幾年來，隨著硬件和數(shù)字信號處理技術(shù)的改進使流量的測量精度得到了不斷提高。在實際應(yīng)用中，時差法超聲波流量計應(yīng)用于大口徑管道水流量的測量中具有明顯的優(yōu)越性。但是，超聲波流量計的應(yīng)用對管道的條件有一定的要求。一般，換能器應(yīng)安裝在上游直管段為10。（為管徑）、下游直管段為50的位置。因此，選擇良好的安裝位省和安裝條件的管道才能獲得滿意的應(yīng)用效果。對于同一個超聲波流

6、量計在應(yīng)用于不同的安裝條件時，流量的測量結(jié)果應(yīng)使用不同的修正系數(shù)加以修正。通常，修正系數(shù)的確定是在實驗室中測試得出的，但這只能針對有限的安裝條件，而且實驗室的建設(shè)和測試的過程既費時又費錢，效率很低。本文針對上面的問題，通過建立超聲波流量計的理論模型和利用Fluent軟件，對不同雷諾系數(shù)的流體和在不同超聲波波束育徑的測最進行了數(shù)值模擬，計算出流量計的流速修正系數(shù)，并通過分析計算的方法計算出修正系數(shù)，兩者進行比較，比較結(jié)果共有很好的一致性。1基本原理1.1基本物理模型對于所有的流體流動，計算機數(shù)值模擬都是解質(zhì)量和動量守恒方程。質(zhì)量守恒方程又稱為連續(xù)性方程：齊況(p.%"(1)該方程是質(zhì)量

7、守恒方程的一般形式，它適用于可壓流動和不可壓流動。源項Sg是從分散的二級相中加入到連續(xù)相的質(zhì)量，源項也可以是如何自定義源項。動量守恒方程在慣性(非加速)坐標系中i方向上表示為：%如)+散媽)=嘿+襄+使"(2)式中：P是靜壓，為是應(yīng)力張，昭和F,分別為，方向上的重力體積力和外部體積力(如離散相相互作用產(chǎn)生的升力)，?，包含了其它的模型相關(guān)源項，如多孔介質(zhì)和自定義源項。1.2Fluent中的紊流模型Fluent中提供了三種紊流模型：Spalart-AHmaras模型、I模型和模型。沒有一個紊流模型對于所有的問題是通用的。選擇模型時主要依靠以下幾點：流體是否可壓、建立特殊的可行的問題、精

8、度的要求、計算機的能力、時間的限制等等。在這里采用的是標準的J模型，它適用范圍廣、經(jīng)濟，有合理的精度，是個半經(jīng)驗公式。標準的k-E模型主要是基于紊流動能和擴散率，其中A方程是個精確方程，£方程是個由經(jīng)驗公式導(dǎo)出的方程°紊流動能方程婦斜幻嚎(她)噫幣+勻)衰+G*+Gb-Ym+S*(3)擴散方程£：.8+安(*)噫S針菁+g4"(&+g&)-Cm，+S.(4)C&是方程中表示由層流速度梯度而產(chǎn)生的紊流動能，&是由浮力產(chǎn)生的紊流動能，匕是由于在可壓縮紊流中過渡的擴散產(chǎn)生的波動，G、和C3是常量，g和J是方程和方程的紊流Pran

9、dtl數(shù)，S&和S,是用戶定義的。1.3修定系數(shù)由于流體的流速沿管道直徑的不均勻分布，致使時差法超聲波流量:計測得的流速并不等于實際的流速，因此，測置出的結(jié)果需乘上一個修正系數(shù)艮IJ：=k(5)式中：虹為超聲波流量計測得的平均流速，Ms為實際的流體流速，其中2k=ku(x9y9z),v(x,y,z),w(x,y,z)(6)這里的k系數(shù)是基于CFD數(shù)值模擬的預(yù)測模型模擬出來的，因此需根據(jù)超聲波流量計的測量原理建立合適的模型，測的原理圖如圖1所示。圖1超聲波流速測量原理圖從圖1可以看出，結(jié)合時差法超聲波流髯計的流速計算方法，可以得出測得的平均流速為：&=虹心,"),心,”

10、)出也勺(7)式中：為超聲波傳播的距離，Q為超聲波傳播方向與管道軸線的夾角，勺和分別是超聲波逆流和順流方向的傳播時間。從圖1中分析，由于流體流速的兩個分量"和"對超聲波傳播方向E都有作用，僅垂直于傳播方向的流速分量沒有作用。因此，要得到A(u,V)則需了解管道內(nèi)流體的流動狀態(tài)，這個可以利用CFD軟件根據(jù)流體在管道入口的初始狀態(tài)和管道的安裝狀況模擬得出。2分析計算和CFD建模2.1分析計算當流體在紊流狀態(tài)時，流速曲線方程為：=卜(1一普)(8)式中：r為截面中的任意點到中心軸的距離，兒皿為中心軸上的流體流速，即截面上的最大流速，p是由雷諾系數(shù)決定的參數(shù)。p=0.25-0.02

11、3lg/?e(9)當流體在層流狀態(tài)時，流速曲線方程為：")="蘋(1-東)(10)從圖1中可以得知，由超聲波測得的流體流速與管徑d、夾角a和超聲波傳播波束的直徑6有關(guān)，因此利用上述的流速曲線方程和超聲波傳播的路徑可以得出流體的平均流速："京b,d)=1Vgd)JJJu(r)dV(H)為了計算出流速的修正系數(shù)，還要得出流體的實際流速。由紊流理論可知，流體截面上各點的流速不(12)(13)Rev/d.圖2管道截面網(wǎng)格劃分樣式圖是均勻分布的，從Prandtl定律中的描述可以知道,流體的實際流速可以由雷諾系數(shù)得出，即：Rev式中：S為管道內(nèi)流體的平均流速，d為管道內(nèi)徑,&

12、quot;為運動粘度(p=8xl0-7m2/s,水溫為30大)。從上面幾式中可以得出人系數(shù)：k(b,d)=%通過改變超聲波傳播范圍的半徑b與管道半徑的比例，可以得出在不同規(guī)格超聲波換能器下測量的流速修正系數(shù)A。這種分析計算模型可以適用于層流和紊流兩種狀態(tài)下。2.2CFD模擬這里利用計算流體力學(xué)軟件Fluent進行管道流體模擬運算，F(xiàn)luent軟件是個工程運用的CFD軟件，針對每一種流動的物理問題的特點，采用適合于它的數(shù)值解法在計算速度、穩(wěn)定性和精度方面達到最佳，可以計算流場、傳熱和化學(xué)反應(yīng)。其思想實際上就是做很多模塊，這樣只要判斷是哪一種流場和邊界就可以拿已有的模型來計算。Fluent軟件能推

13、出多種優(yōu)化的物理模型，如定常和非定常流動、層流(包括各種非牛頓流模型)、紊流(包括最先進的紊流模型)、不可壓縮和可壓縮流動、傳熱、化學(xué)反應(yīng)等等。對每一種物理問題的流動特點，有適合它的數(shù)值解法，用戶可對顯式或隱式差分格式進行選擇，以期在計算速度、穩(wěn)定性和精度等方面達到最佳。Fluent將不同領(lǐng)域的計算軟件組合起來，成為CFD計算機軟件群，軟件之間可以方便地進行數(shù)值交換，并采用統(tǒng)一的前、后處理工具，這就省卻了科研工作者在計算方法、編程、前后處理等方面投入的重復(fù)、低效的勞動，而可以將主要精力和智慧用于物理問題本身的探索上。在Fluent軟件中，采用Gambit的專用前處理軟件，使網(wǎng)格可以有多種形狀。

14、對二維流動，可以生成三角形和矩形網(wǎng)格；對于三維流動，則可生成四面體、六面體、三角柱和金字塔等網(wǎng)格；結(jié)合具體計算，還可生成混合網(wǎng)格，其自適應(yīng)功能，能對網(wǎng)格進行細分或粗化，或生成不連續(xù)網(wǎng)格、可變網(wǎng)格和滑動網(wǎng)格。網(wǎng)格劃分的好壞將直接影響模擬的精度和收斂性，本文管道截面的網(wǎng)格劃分采用的是矩形網(wǎng)格，外層采用的是同心圓的排列方式，里層采用的是交叉排列方式，式樣如圖2所示，這種混合網(wǎng)格排列的方式有利于管道邊界層的模擬效果。對于管道流速修正系數(shù)A3,v)可以通過Fluent軟件模擬的流速數(shù)據(jù)計算出來，根據(jù)網(wǎng)格劃分的體積權(quán)重來計算超聲波傳播路徑上的平均流速&“:布=；(14)式中：耳和與可以從Fluen

15、t輸出的數(shù)據(jù)中得到，*為體積，如圖1中坐標系所標的u和v兩個速度分量。根據(jù)Fluent的數(shù)值模擬可以得出流速修正系數(shù)v):UCOSGC,、k(uyv)=(15)Utim式中：偵是流體的實際流速，&g是流體在超聲波傳播路徑上的平均流速，a是超聲波傳播路徑和管道軸線的夾角。3試驗結(jié)果分析設(shè)定的邊界條件為：管徑為100mm,長度L為50D,以保證在數(shù)值模擬中流體能夠得到充分流動，入口流體流速設(shè)定為lm/so采用三種方式通過分析計算和數(shù)值模擬來計算流速修正系數(shù)K值：(1) 當流體雷諾系數(shù)Rc=100,流體流動狀態(tài)在層流的情況下，通過改變超聲波傳播波束的直徑與管徑的比率b/d,計算出流速修正系數(shù)

16、左值，結(jié)果如圖3所示。(2) 當流體雷諾系數(shù)Re=10000,流體流動狀態(tài)在紊流的情況下，同樣通過改變6/d值，計算出流速修正系數(shù)&值，結(jié)果如圖4所示。(3) 在超聲波傳播波束的直徑與管徑的比率b/d=0.5的情況下，通過改變流體雷諾系數(shù)Re的大小來計算流速修正系數(shù)&值，結(jié)果如圖5所示。1.000.980.960.940.920.900.20.40.6bld0.81.0圖4Re=10000時流速修正系數(shù)曲線圖1.000.980.960.940.92。分析計算。數(shù)值模擬O-O-o0.9011!1,.20000400006000080000100000雷諾系數(shù)：Re圖5不同雷諾系數(shù)

17、下的流速修正系數(shù)圖從圖3可以看出，當超聲波波束直徑較小時，修正系數(shù)約為0.78,而當超聲波波束直徑較大時，約為0.95,兩種方法計算出的結(jié)果相差小于1%。從圖4中可以看出，修正系數(shù)最小約為0.93,最大達到約0.98,兩種方法計算出的結(jié)果具有很好的一致性。從圖3、圖4中分析，修正系數(shù)接近1,由于超聲波傳播路徑與管軸存在著一定的夾角關(guān)系，因此修正系數(shù)不會等于1。從圖3、圖4中還可以看出，當6/d值變小時,分析計算和數(shù)值模擬兩種方法得出的修正系數(shù)相差也逐漸加大，這是因為在數(shù)值模擬中的網(wǎng)格劃分導(dǎo)致的，數(shù)值模擬區(qū)域的實際體積與它包含的網(wǎng)格體積存在若一定的偏差，當b/d值較小時，這個偏差也較大，因此出現(xiàn)

18、這種情況。從圖5中可以看出，當雷諾系數(shù)從20000到100000變化時，分析計算和數(shù)值模擬得出流速修正系數(shù)結(jié)果都很接近，而且受雷諾系數(shù)變化的影響很小。4結(jié)論Fluent軟件在我國已經(jīng)獲得了廣泛的應(yīng)用。由于其減少了研究者在計算方法、編程、前后處理等方面投入的重復(fù)、低效的勞動，將更多的精力和時間投入到考慮問題的物理本質(zhì)、優(yōu)化算法選用、參數(shù)的設(shè)定，因而提高了工作效率，其必然會獲得越來越多的應(yīng)用，應(yīng)用效果也必將越來越好。從試驗結(jié)果中可以看出，分析計算和數(shù)值模擬兩種方法得出的修正系數(shù)具有很好的一致性，采用數(shù)值模擬的方法可以很好地模擬預(yù)測出時差法超聲波流量計在不同的安裝環(huán)境和條件下流速的修正系數(shù)，對于在非標準而復(fù)雜的測量環(huán)境中，數(shù)值模擬