賈 華， 郭 朔， 王月明

(內(nèi)蒙古科技大學(xué) 信息工程學(xué)院,內(nèi)蒙古 包頭 014010)

?

研究與探討

超聲法測量氣固兩相流傳感器選取研究*

賈 華， 郭 朔， 王月明

(內(nèi)蒙古科技大學(xué) 信息工程學(xué)院,內(nèi)蒙古 包頭 014010)

超聲法因其穿透性好、精度高等優(yōu)點(diǎn)成為氣固兩相流參數(shù)檢測的一個(gè)新趨勢。由于超聲波在空氣中衰減較大，因此，選擇合適的超聲波傳感器對(duì)于氣固兩相流參數(shù)檢測至關(guān)重要。通過運(yùn)用有限元軟件COMSOL構(gòu)建超聲法測量氣固兩相流的仿真模型，研究接收裝置接收信號(hào)的大小與氣固兩相流固相顆粒的粒徑和發(fā)射頻率的關(guān)系，為實(shí)際測量選取超聲波傳感器提供了理論依據(jù)和指導(dǎo)建議。

氣固兩相流； 有限元； 粒徑； 頻率； 超聲波傳感器

0 引 言

在電力、化工、冶金等工業(yè)生產(chǎn)過程中，應(yīng)用管道氣力輸送技術(shù)輸送煤粉、水泥、礦石等均屬于典型的氣固兩相流動(dòng)[1,2]。由于氣固兩相流涉及領(lǐng)域的廣泛性以及過程參數(shù)的重要性使得對(duì)于氣固兩相流的研究工作得到了迅速發(fā)展。

目前，國內(nèi)外學(xué)者對(duì)于氣固兩相流相關(guān)參數(shù)的測量所提出的方法可分為兩類，接觸式測量和非接觸式測量。接觸式測量方法主要有：傳熱法、摩擦電法、能量法、差壓法等[3]。由于接觸式測量會(huì)對(duì)管內(nèi)流體產(chǎn)生干擾和增加原有阻力，在測量時(shí)會(huì)降低測量的精度。非接觸式測量方法主要有：光學(xué)法[4]、電容法[5,6]、微波法[7]、超聲法等[8]。超聲法作為一種非接觸式測量方法，具有穿透性好、結(jié)構(gòu)簡單、無污染、精度高、應(yīng)用范圍廣等較多優(yōu)勢，在近年來成為研究熱點(diǎn)并得到一定的發(fā)展[9,10]。上海理工大學(xué)的蘇明旭和蔡小舒在超聲法測量氣固兩相流方面做了大量研究工作并取得較大進(jìn)展[11～13]。

本文運(yùn)用有限元軟件COMSOL構(gòu)建了超聲法測量氣固兩相流的仿真模型，通過仿真，研究接收裝置接收信號(hào)的大小與氣固兩相流固相顆粒的粒徑和發(fā)射頻率的關(guān)系，為實(shí)際測量選取超聲波傳感器提供理論依據(jù)和指導(dǎo)建議。

1 仿真建模

在實(shí)際檢測過程中，固相顆粒隨機(jī)分布在氣固兩相流體中。為了獲取固相顆粒粒徑d和超聲波發(fā)射頻率f對(duì)超聲波接收信號(hào)的影響，仿真過程中，在測量管道中心放置不同粒徑的固相顆粒，研究不同粒徑和不同發(fā)射頻率對(duì)接收裝置接收信號(hào)的影響情況。

本文選取610 mm直徑的工業(yè)輸粉管道作為研究對(duì)象，為簡化模型，取管道長度1 000 mm，并對(duì)其縱剖面進(jìn)行仿真建模，基于超聲法測量管道中存在固相顆粒的氣固兩相流流體仿真模型如圖1所示。超聲波發(fā)射和接收傳感器剖面半徑為10 mm，厚度為20 mm，傳感器由復(fù)合壓電材料PZT—5H構(gòu)成，在仿真建模時(shí)將發(fā)射和接收傳感器進(jìn)行良好對(duì)中。輸粉管道中連續(xù)相介質(zhì)采用軟件內(nèi)置材料空氣，離散相選用電廠煤粉顆粒。常溫下，煤粉—空氣兩相流的其余相關(guān)參數(shù)如表1所示。在進(jìn)行有限元分析時(shí)，網(wǎng)格的質(zhì)量將直接影響整個(gè)過程的計(jì)算精度，因此，在網(wǎng)格劃分時(shí)，對(duì)輸粉管道、超聲波發(fā)射和接收傳感器以及固相顆粒分別進(jìn)行細(xì)化、較細(xì)化和特別細(xì)化處理。圖2為網(wǎng)格劃分后的仿真模型，本文采用自由剖分三角形方式。

圖1 管道中存在固相顆粒的仿真模型

參數(shù)空氣煤粉密度/(kg/m3)1.231702聲速/(m/s)3431800比熱容/(J/(kg·K))1004850導(dǎo)熱系數(shù)0.0250.89膨脹系數(shù)3.47×10-31.0×10-5剪切粘度1.8×10-59.04×10-4

圖2 網(wǎng)格劃分后的仿真模型

1.1 多物理場耦合方程

在仿真過程中，由激勵(lì)裝置產(chǎn)生激勵(lì)電壓驅(qū)動(dòng)超聲波發(fā)射傳感器產(chǎn)生超聲波，超聲波經(jīng)管道傳播到達(dá)接收傳感器，接收傳感器將衰減后的聲波信號(hào)轉(zhuǎn)換為電信號(hào)以供采集。因此,有限元分析主要涉及電—結(jié)構(gòu)—聲三物理場耦合，其中聲場波動(dòng)方程可表示為

(1)

結(jié)構(gòu)力學(xué)方程可表示為

(2)

式中 ρ為壓電材料密度； u為位移；s為應(yīng)力；FV為體積力。

電場的麥克斯韋方程可表示為

(3)

式中 D為電位移; ρV為體積電荷密度。

在進(jìn)行有限元分析時(shí)，設(shè)定發(fā)射傳感器內(nèi)側(cè)邊界為發(fā)射端，其電邊界條件為：幅值為300V的瞬時(shí)激勵(lì)電壓。設(shè)定接收傳感器內(nèi)側(cè)邊界為接收端，并設(shè)定為終端。發(fā)射端與接收端的聲邊界條件為聲—結(jié)構(gòu)耦合，外側(cè)邊界的電邊界條件為接地，結(jié)構(gòu)邊界條件為輥支撐。

2 仿真過程與分析

2.1 仿真過程

在測量管道中心放置不同粒徑的固相顆粒，研究不同粒徑和不同發(fā)射頻率對(duì)接收裝置接收信號(hào)的影響情況。為了研究顆粒粒徑d對(duì)接收傳感器接收電壓的影響，仿真過程中，設(shè)定d由小到大分別為0，0.02，0.04，0.06，0.08，0.1，0.15，0.2mm，并對(duì)接收傳感器接收電壓的數(shù)據(jù)進(jìn)行采集。同時(shí)為了研究發(fā)射頻率f對(duì)接收傳感器接收電壓的影響，設(shè)定f由小到大分別為20，40，60，80，100，150，200kHz，并對(duì)接收傳感器接收電壓的數(shù)據(jù)進(jìn)行記錄。限于篇幅，本文只給出部分仿真結(jié)果。

圖3(a)為d=0mm,f=40kHz時(shí)接收電壓的變化情況。圖3(b)為d=0.1mm,f=40kHz時(shí)接收電壓的變化情況。從圖3可以看出，存在顆粒時(shí)對(duì)接收電壓產(chǎn)生一定影響，但無法作出具體判斷。

圖3 有無顆粒存在時(shí)的接收電壓

2.2 分析與討論

為了進(jìn)一步研究氣固兩相流固相顆粒粒徑和發(fā)射頻率對(duì)接收電壓的影響，需要對(duì)仿真數(shù)據(jù)進(jìn)行處理，圖4和圖5為進(jìn)行數(shù)據(jù)處理后得到的結(jié)果。由圖4可以看出：在同一頻率下，接收電壓隨著顆粒粒徑的增大而減小，但衰減較為緩慢；由圖5可以看出：在相同的顆粒粒徑下，接收電壓隨著發(fā)射頻率的增大而迅速減小，說明接收電壓對(duì)于頻率更為敏感。仿真發(fā)現(xiàn)：當(dāng)超聲波的頻率在一定范圍內(nèi)(20～100 kHz)，接收傳感器的接收電壓相對(duì)較大。在實(shí)際測量中，由于不可避免地會(huì)引入噪聲，而大多數(shù)噪聲信號(hào)頻率在50 kHz以下，因此，在測量中應(yīng)當(dāng)選用頻率在50～100 kHz的超聲波傳感器。

圖4 粒徑與接收電壓的關(guān)系

圖5 發(fā)射頻率與接收電壓的關(guān)系

3 結(jié) 論

運(yùn)用有限元軟件COMSOL構(gòu)建了超聲法測量氣固兩相流的仿真模型，研究接收傳感器的接收電壓與氣固兩相流固相顆粒的粒徑和發(fā)射頻率的關(guān)系。仿真結(jié)果表明：選用頻率在50～100 kHz的超聲波傳感器不僅能夠得到較高幅值的接收信號(hào)還能有效避免噪聲信號(hào)的引入，為實(shí)際測量中超聲波傳感器的選取提供了理論依據(jù)和指導(dǎo)建議。

[1] 連桂森.多相流動(dòng)基礎(chǔ)[M].杭州:浙江大學(xué)出版社,1989:1-2.

[2] 李海青.兩相流參數(shù)檢測及應(yīng)用[M].杭州:浙江大學(xué)出版社,1991:1-2.

[3] 李光輝.電廠鍋爐一次風(fēng)煤粉濃度測量系統(tǒng)研究[D].南京:南京理工大學(xué),2007:4-6.

[4] 高 陽,楊 斌,何 淵,等.電廠煤粉細(xì)度在線監(jiān)測系統(tǒng)研制與應(yīng)用[J].中國粉體技術(shù),2015,21(2):28-31.

[5] 李文濤,葉 俊,李忠虎,等.氣/固兩相流固體顆粒流速測量方法研究[J].儀表技術(shù)與傳感器,2015(3):88-91.

[6] 李 鵬,楊道業(yè),史華瑩,等.電容式氣固兩相流的參數(shù)測量系統(tǒng)[J].儀表技術(shù)與傳感器,2015(5):80-83.

[7] 張春堂,王軍政,王國棟.高爐噴煤在線監(jiān)測系統(tǒng)的設(shè)計(jì)[J].化工自動(dòng)化及儀表,2011(6):695-698.

[8] 金秀章,陳小剛,吳元林.基于氣固兩相流體理論求解風(fēng)粉濃度[J].計(jì)算機(jī)仿真,2013,30(12):121.

[9] 田 昌,蘇明旭,蔡小舒.基于超聲法測量氣固兩相流濃度實(shí)驗(yàn)研究[J].工程熱物理學(xué)報(bào),2013,34(8):1487-1490.

[10] 李永明,蘇明旭,袁安利,等.基于超聲的管道內(nèi)粉體體積分?jǐn)?shù)的測量[J].中國科學(xué)院大學(xué)學(xué)報(bào),2016,33(2):277-282.

[11] 薛明華,蘇明旭,蔡小舒.超聲波多次回波反射法測量兩相流密度實(shí)驗(yàn)研究[J].工程熱物理學(xué)報(bào),2008,29(8):1343-1346.

[12] 李永明,蘇明旭,周健明,等.超聲波測量電廠煤粉管道內(nèi)煤粉體積分?jǐn)?shù)的研究[J].動(dòng)力工程學(xué)報(bào),2016,36(3):201-206.

[13] 田 昌,蘇明旭,顧建飛.氣固兩相流超聲過程層析成像系統(tǒng)研究[J].儀器儀表學(xué)報(bào),2016,37(3):586-592.

Research on selection of gas-solid two phase flow transducer with ultrasonic measurement*

JIA Hua， GUO Shuo， WANG Yue-ming

(School of Information Engineering，Inner Mongolia University of Science and Technology，Baotou 014010，China)

The ultrasonic method has become a new trend of gas-solid two phase flow parameter measurement because of its excellent penetration and high precision.Due to the large attenuation of ultrasonic wave in air,it is of great importance to select the suitable ultrasonic transducer for the parameter measurement of gas-solid two phase flow.The simulation model for gas-solid two phase flow measurement is established by the finite element software COMSOL,the relationship between signal magnitude recieved by receiving device and the solid phase particle size of gas-solid two phase flow and transmitting frequency is studied,which provides theoretical basis and guidance for choosing the ultrasonic transducer in practical measurement.

gas-solid two phase flow； finite element； particle size； frequency； ultrasonic transducer

10.13873/J.1000—9787(2017)07—0005—03

2016—09—05

國家自然科學(xué)基金資助項(xiàng)目(61463042)； 內(nèi)蒙古自然科學(xué)基金資助項(xiàng)目(2014MS0609)

TP 212

A

1000—9787(2017)07—0005—03

賈 華(1961-)，男，副教授，從事控制理論與控制工程研究工作。

郭 朔，通訊作者，碩士研究生，研究方向?yàn)闅夤虄上嗔鳈z測,E—mail：cavicent@163.com。