摘要:渦輪流量計(jì)在脈動(dòng)流量下動(dòng)態(tài)特性，利用FLUENT軟件對(duì)渦輪流量計(jì)內(nèi)脈動(dòng)流場(chǎng)進(jìn)行仿真計(jì)算。研究中獲得了流量計(jì)在不同脈動(dòng)幅值和頻率下的瞬時(shí)輸出流量，通過(guò)正弦函數(shù)擬合獲得各工況輸出流量的平均值、脈動(dòng)幅值和初始相位，進(jìn)而對(duì)渦輪流量計(jì)幅頻特性和相頻特性進(jìn)行了分析，幅頻特性隨脈動(dòng)頻率成線性降低的趨勢(shì)，而相頻特性隨頻率增加而增大而后趨于穩(wěn)定。
0引言
　　渦輪流量計(jì)是典型的速度式流量計(jì)「1-2)，通過(guò)測(cè)量葉輪的轉(zhuǎn)速來(lái)計(jì)算管道中流體體積流量。穩(wěn)定流動(dòng)下渦輪流量計(jì)具有重復(fù)性好、量程范圍寬、適應(yīng)性強(qiáng)、精度高、體積小等特點(diǎn)。但在工業(yè)領(lǐng)域的流量測(cè)量中，非穩(wěn)態(tài)流動(dòng)工況廣泛存在，如脈動(dòng)流量。脈動(dòng)流量可由旋轉(zhuǎn)式、往復(fù)式或其他可運(yùn)動(dòng)傳送設(shè)備而產(chǎn)生，脈動(dòng)流動(dòng)--旦形成就會(huì)在流體中傳播，將會(huì)對(duì)渦輪流量計(jì)產(chǎn)生較大的影響田。其脈動(dòng)幅度和脈動(dòng)頻率的變化均會(huì)對(duì)渦輪流量計(jì)動(dòng)態(tài)特性產(chǎn)生一定的影響。頻率低、幅度小的脈動(dòng)流，一般情況下對(duì)流量測(cè)量影響不大，但若脈動(dòng)頻率較高或幅度較大時(shí)則將對(duì)流量計(jì)的幅頻特性和相頻特性產(chǎn)生很大影響。
　　利用CFD軟件對(duì)渦輪流量計(jì)在正弦脈動(dòng)流動(dòng)下的流場(chǎng)進(jìn)行仿真計(jì)算，獲得不同脈動(dòng)頻率和幅度下渦輪流量計(jì)角速度曲線，進(jìn)而計(jì)算瞬時(shí)流量，對(duì)渦輪流量計(jì)動(dòng)態(tài)特性進(jìn)行分析研究。
1渦輪流量計(jì)建模
1.1幾何模型
　　研究中針對(duì)DN32口徑的液體渦輪流量計(jì)展開(kāi)，其機(jī)芯內(nèi)部結(jié)構(gòu)如圖1所示，主要包括前后導(dǎo)向件和葉輪三部分。對(duì)渦輪流量計(jì)實(shí)物的外形尺寸和位置關(guān)系等關(guān)鍵尺寸進(jìn)行測(cè)繪(表1為渦輪流量計(jì)主要尺寸參數(shù))，忽略流量計(jì)進(jìn)出口連接形式等次要因素的影響，繪制三維結(jié)構(gòu)圖如圖2所示。
 
 
1.2流體仿真模型
　　利用Gambit軟件進(jìn)行幾何建模和網(wǎng)格劃分，為使流動(dòng)在流量計(jì)入口處達(dá)到充分發(fā)展的流動(dòng)狀態(tài)，分別在流量計(jì)的.上游設(shè)置10倍長(zhǎng)直管段，下游設(shè)置5倍直管段。將葉輪所在區(qū)域定義為旋轉(zhuǎn)區(qū)域，其余部分定.義為靜止區(qū)域，通過(guò)interface面將動(dòng)靜區(qū)域進(jìn)行連接。旋轉(zhuǎn)區(qū)域及結(jié)構(gòu)較復(fù)雜的區(qū)域采用非結(jié)構(gòu)化網(wǎng)格;直管段等結(jié)構(gòu)較簡(jiǎn)單的區(qū)域采用結(jié)構(gòu)化網(wǎng)格。
仿真中選用ReynoldsStress(S-BLS)湍流模型，該模型是最符合物理現(xiàn)象的模型，各向異性，輸運(yùn)中的雷諾應(yīng)力可直接計(jì)算出來(lái)。模型計(jì)算時(shí)間較長(zhǎng)，適合大彎曲流線、漩渦以及三維轉(zhuǎn)動(dòng)流動(dòng)。選用RP-3航空煤油作為流體介質(zhì)，以實(shí)測(cè)介質(zhì)密度和粘度并設(shè)置仿真流體參數(shù)。設(shè)置下游直管段出口為自由出流(out-flow)條件，直管段及流量計(jì)殼體為靜止壁面(wall);上游直管段入口為速度入口(velocity-inlet)。
1.3仿真方法
　　渦輪流量計(jì)葉輪受到流體沖擊時(shí)，對(duì)葉輪產(chǎn)生驅(qū)動(dòng)力矩，同時(shí)受到各種阻力矩的影響圖3是葉輪所受力矩示意圖。葉輪轉(zhuǎn)動(dòng)過(guò)程中受到流體對(duì)葉片的驅(qū)動(dòng)力矩Td，輪轂側(cè)面受到流體粘性阻力矩Th，輪轂端面受到流體粘性阻力矩Tw，葉片頂端受到流體粘性阻力矩Tt，同時(shí)葉片輪軸與軸承之間存在機(jī)械阻力矩Tb，磁電信號(hào)檢出器產(chǎn)生磁電阻力矩Tm。
 
　　式中:J為葉輪轉(zhuǎn)動(dòng)慣量;ɷ為葉輪轉(zhuǎn)動(dòng)角速度;t為運(yùn)行時(shí)間，Trf為流體阻力矩，Trf=Th+Tw+Tt。
　　采用6DOF模型實(shí)現(xiàn)對(duì)葉輪6個(gè)自由度的控制，包括X，Y，Z方向的移動(dòng)自由度和圍繞X，Y，Z軸的旋轉(zhuǎn)自由度。通過(guò)DEFINE_SDOF_PROPERTIES宏文件約束葉輪在X，Y，Z方向上的移動(dòng)以及圍繞X，Y軸的旋轉(zhuǎn)，只能圍繞Z軸(流動(dòng)方向)進(jìn)行轉(zhuǎn)動(dòng)。葉輪三維建模中選用的材質(zhì)與實(shí)際材質(zhì)相同，即可獲得葉輪的轉(zhuǎn)動(dòng)慣量。流體仿真中系統(tǒng)自行計(jì)算流體驅(qū)動(dòng)力矩、流體阻力矩，忽略軸承摩擦阻力矩和磁電阻力矩。
　　仿真中通過(guò)UDF函數(shù)對(duì)上游直管段入口流速進(jìn)行設(shè)置，流速通過(guò)式(2)和式(3)計(jì)算。
 
　　式中:qt為入口瞬時(shí)體積流量;q0為體積流量平均值，Q0=16m3/h;A為脈動(dòng)幅度;ƒ為脈動(dòng)頻率;t0為脈動(dòng)流量起始時(shí)刻，t0=0.0132s;Vt為瞬時(shí)入口速度平均值;Ain為上游直管段入口截面積。
　　仿真初始葉輪轉(zhuǎn)速為零，根據(jù)葉輪運(yùn)動(dòng)方程自動(dòng)計(jì)算和調(diào)整旋轉(zhuǎn)角速度，直到葉輪趨于穩(wěn)定狀態(tài)，仿真中各工況在t,前葉輪轉(zhuǎn)速已經(jīng)達(dá)到穩(wěn)定狀態(tài)。
2仿真結(jié)果與試驗(yàn)驗(yàn)證
　　仿真計(jì)算中調(diào)整脈動(dòng)幅度(A=5，10，15,20L/min)和脈動(dòng)頻率(ƒ=5，10，15，20，30，40，50Hz)兩參數(shù)設(shè)置，計(jì)算28個(gè)不同工況下流場(chǎng)，獲得葉輪轉(zhuǎn)速隨時(shí)間變化曲線，圖4時(shí)脈動(dòng)頻率為20Hz和50Hz時(shí)的葉輪轉(zhuǎn)速曲線。
　　待流量計(jì)葉輪旋轉(zhuǎn)平穩(wěn)后，取to時(shí)刻前的葉輪轉(zhuǎn)速通過(guò)式(4)計(jì)算儀表系數(shù)。
 
　　式中:K為儀表系數(shù)，L-1;ɷ∞為穩(wěn)定葉輪轉(zhuǎn)速，rad/s;N為渦輪葉片數(shù)量，N=6;q0為入口平均流量，m3/h。
　　經(jīng)計(jì)算DN32渦輪流量計(jì)在16m3/h流量下，儀表系數(shù)K=160.6L-1。依據(jù)JJG1037-2008《渦輪流量計(jì)檢定規(guī)程》，利用流量標(biāo)準(zhǔn)裝置對(duì)渦輪流量計(jì)進(jìn)行校準(zhǔn)試驗(yàn)，獲得16m3/h流量下儀表系數(shù)K=161.7L-1，與仿真結(jié)果的儀表系數(shù)相差僅0.7%，仿真與校準(zhǔn)試驗(yàn)一致性非常好。
 
3動(dòng)態(tài)特性分析
　　渦輪流量計(jì)瞬時(shí)輸出流量可由式(5)計(jì)算。
 
　　式中:qm為渦輪流量計(jì)瞬時(shí)流量，m3/h;ɷ為葉輪瞬時(shí)轉(zhuǎn)速，rad/s.
　　由于流量計(jì)前后導(dǎo)流件的作用，導(dǎo)致葉輪轉(zhuǎn)速和流量計(jì)瞬時(shí)流量存在小幅波動(dòng)，圖5是ƒ=40Hz渦輪流量計(jì)瞬時(shí)流量曲線，通過(guò)正弦函數(shù)對(duì)流量計(jì)瞬時(shí)流量進(jìn)行擬合，擬合函數(shù)為
 
　　式中:qm0為流量計(jì)輸出瞬時(shí)流量平均值，m3/h;ƒm為流量計(jì)輸出脈動(dòng)頻率，Hz;Am為流量計(jì)輸出脈動(dòng)幅度;tm0為流量計(jì)輸出脈動(dòng)流量起始時(shí)刻，:S。
　　仿真中28組工況擬合結(jié)果決定系數(shù)都大于0.995，擬合效果非常好，擬合獲得了各工況渦輪流量計(jì)輸出流量的平均值qm0、脈動(dòng)頻率ƒm、脈動(dòng)幅度Am和脈動(dòng)起始時(shí)刻tm0四個(gè)參數(shù)。流量計(jì)輸出幅頻特性和相頻特性分別通過(guò)式(7)和式(8)計(jì)算。
 
 
　　將各工況正弦擬合曲線作為動(dòng)態(tài)輸出，與入口瞬時(shí)流量進(jìn)行比較，圖6是脈動(dòng)頻率分別為5，50Hz時(shí)渦輪流量計(jì)動(dòng)態(tài)信號(hào)曲線圖，可見(jiàn)脈動(dòng)幅值對(duì)相位差基本無(wú)影響。各工況流量計(jì)輸出流量的平均值比較恒定，變化范圍是16.079~16.094m3/h，比輸入流量平均值偏大不足0.6%。圖7是渦輪流量計(jì)幅頻特性與相頻特性曲線圖。A*隨脈動(dòng)頻率ƒ基本成線性降低的趨勢(shì)，低頻情況下A*≈1，且略大于1;隨脈動(dòng)頻率增加，A*逐漸減小，在ƒ=50Hz時(shí)A*≈0.8。對(duì)于相位差而言，在ƒ<40Hz的范圍內(nèi)，相位差隨脈動(dòng)頻率增加而增大，相位差由3.7°~4.8°增大至20.9°~24.2°;在ƒ=40Hz之后，相位差趨于平穩(wěn)。葉輪是渦輪流量計(jì)內(nèi)的唯---可動(dòng)部件，葉輪在流體驅(qū)動(dòng)力矩和阻力矩作用下產(chǎn)生加速度，由于葉輪自身轉(zhuǎn)動(dòng)慣量的影響，導(dǎo)致葉輪角速度變化滯后于流量脈動(dòng)流動(dòng)變化，產(chǎn)生流量示值與真實(shí)之間的相位差，而相位差受到流量計(jì)自身時(shí)間常數(shù)和流量脈動(dòng)輸入雙方面的綜合影響，渦輪流量計(jì)一般認(rèn)為是一階非線性系統(tǒng)[1o]，其時(shí)間常數(shù)是隨流量輸入而改變的。
 
 
4結(jié)論
　　研究中利用FLUENT軟件中的6DOF模型對(duì)DN32口徑渦輪流量計(jì)進(jìn)行流體仿真，仿真過(guò)程中改變?nèi)肟诹髁棵}動(dòng)幅值和脈動(dòng)頻率，共獲得28組工況脈動(dòng)流動(dòng)下的流場(chǎng)。分析獲得葉輪轉(zhuǎn)速隨時(shí)間變化曲線，并利用正弦函數(shù)對(duì)流量計(jì)輸出進(jìn)行擬合，獲得流量計(jì)瞬時(shí)流量的平均值、脈動(dòng)幅值和初始相位等參數(shù)，進(jìn)一步獲得了流量計(jì)幅頻特性和相頻特性。各工況流量計(jì)輸出流量平均值可認(rèn)為是定值，幅頻特性隨脈動(dòng)頻率成線性降低的趨勢(shì)，而相頻特性受到渦輪流量計(jì)時(shí)間常數(shù)和脈動(dòng)輸入綜合影響，在ƒ<40Hz的范圍內(nèi)隨頻率增加而增加，在ƒ=40Hz之后，相頻特性趨于平穩(wěn)。
　　本仿真研究中，渦輪流量計(jì)規(guī)格單一且工況較少，未來(lái)還需要結(jié)合理論分析、實(shí)流試驗(yàn)、流體仿真等多種手段對(duì)渦輪流量計(jì)在脈動(dòng)流動(dòng)下動(dòng)態(tài)特性開(kāi)展研究，進(jìn)一步探究渦輪流量計(jì)動(dòng)態(tài)響應(yīng)機(jī)理和修正方法，提高渦輪流量計(jì)在脈動(dòng)流量測(cè)量中的精度。

本文來(lái)源于網(wǎng)絡(luò),如有侵權(quán)聯(lián)系即刪除！