韓中合， 張美鳳， 錢江波

（華北電力大學(xué) 電站設(shè)備狀態(tài)監(jiān)測與控制教育部重點實驗室，保定071003）

火電廠大功率凝汽式汽輪機(jī)低壓缸的末幾級和水冷堆核電汽輪機(jī)的全部級都工作在濕蒸汽區(qū).蒸汽帶水一方面會使機(jī)組效率降低，另一方面會對透平低壓級組動葉片產(chǎn)生強(qiáng)烈的侵蝕與沖擊，導(dǎo)致葉片表面出現(xiàn)凹坑甚至扭曲斷裂，嚴(yán)重威脅汽輪機(jī)的安全運行［1－4］.長期以來，各國都在致力于研究準(zhǔn)確測量蒸汽濕度的方法并設(shè)法減少濕蒸汽對汽輪機(jī)所造成的危害.

基于微波諧振腔的介質(zhì)微擾，諧振腔的諧振頻率隨腔內(nèi)電介質(zhì)的介電常數(shù)變化將發(fā)生偏移.在一定溫度（或壓力）下，濕蒸汽的介電常數(shù)只與其濕度有關(guān)，因此當(dāng)濕蒸汽流過微波諧振腔時，通過測量諧振腔諧振頻率的偏移，可以間接測量濕蒸汽的介電常數(shù)，進(jìn)而可以確定蒸汽的濕度［5－7］.

為使測量結(jié)果能準(zhǔn)確反映汽輪機(jī)末級后的排汽濕度，要求諧振腔傳感器取樣盡可能反映真實的流動情況.考慮流場的氣動性能，針對微波諧振腔傳感器兩種不同的取樣前端，分析了不同流速和不同水滴粒徑下由于非等動能取樣而產(chǎn)生的蒸汽濕度測量誤差.

1 數(shù)學(xué)模型的建立

微波諧振腔傳感器是一個兩端開口的圓筒形腔體，在兩端開口處設(shè)置等距離等厚度圓環(huán)組成的環(huán)形柵格分隔器.這種型式的分隔器結(jié)構(gòu)既可以保證濕蒸汽兩相流自由通過諧振腔，又可保證諧振腔具有較高的品質(zhì)因數(shù).取樣進(jìn)口附近是復(fù)雜的三維湍流流動，由于諧振腔和取樣前端是軸對稱結(jié)構(gòu)，而且主流方向與腔體軸線方向一致，故可簡化為縱向剖面的二維問題.帶楔形取樣前端的諧振腔縱向剖面圖見圖1.

圖1 楔形入口的諧振腔縱向剖面圖Fig.1 The longitudinal diagram of microwave cavity with wedge－shaped entrance

1.1 非等動能取樣誤差公式推導(dǎo)

汽輪機(jī)內(nèi)的濕蒸汽是由飽和蒸汽和小水滴組成的混合物.其中，大部分水滴是通過自發(fā)凝結(jié)、生長過程形成的一次水滴（直徑為0.01～2μm），這部分水滴的數(shù)量巨大，可達(dá)107個／cm3，占濕蒸汽中液相質(zhì)量的90%以上，其余為水膜在氣流作用下撕裂形成的二次水滴［4，8－9］（直徑為10～200μm）.

定義排汽的體積濕度φ與質(zhì)量濕度Y之間的關(guān)系為

式中：ρf為飽和水的密度；ρv為干飽和蒸汽的密度.

設(shè)圓柱形諧振腔的半徑為R，不存在取樣誤差時，離散相的體積流量為

式中：m為離散相的質(zhì)量流量.體積濕度為

式中：Vv為蒸汽的體積流量.

存在非等動能取樣誤差時，離散相的體積流量為

式中：r表示諧振腔半徑與取樣汽流半徑之差.體積濕度為

定義取樣偏差為

式中：Y′表示存在非等動能取樣誤差時的蒸汽質(zhì)量濕度.

1.2 諧振腔兩相流模擬的控制方程

1.2.1 汽相控制方程

汽相控制方程［10］的通用形式如下

式中：ρ、φ、Γ、v和s分別為流體的密度、通用變量、擴(kuò)散系數(shù)、速度矢量和源項.

在三維圓柱軸對稱坐標(biāo)系中，汽流在各個方向的分量為u、v和w，則式（7）展開得：連續(xù)性方程

動量方程

湍流模型k方程

湍流模型ε方程

式中：Gk為單位體積干飽和蒸汽由于層流速度梯度產(chǎn)生的湍流動能；系數(shù)c1、c2分別取為1.44和1.92［11］.

能量方程

式中：h、Γh分別為比焓和比焓交換系數(shù).

采用Yong［12］提出的蒸汽狀態(tài)方程：

采用基于控制體的有限體積法將上述微分方程組離散化轉(zhuǎn)換為代數(shù)方程組.該方法在每一個控制體內(nèi)對控制方程積分，從而產(chǎn)生基于控制體的每個變量都守恒的離散方程.以標(biāo)量ψ的定常狀態(tài)守恒方程說明控制方程的離散，對于控制體V的積分方程如下式中和Sψ分別為曲面面積矢量、ψ的擴(kuò)散系數(shù)、ψ的梯度和每一單位體積ψ的源項.

1.2.2 液相控制方程

采用拉格朗日方法計算水滴的運動軌跡，為了簡化問題，進(jìn)行如下假設(shè)［10］：

（1）水滴與蒸汽間無傳熱和傳質(zhì)過程，且不發(fā)生相變；

（2）水滴為均一尺寸的圓球，均勻分布；

（3）水滴與諧振腔壁面碰撞即被捕獲，不存在反彈；

（4）水滴之間不存在碰撞、聚合和撕裂.

水滴的運動軌跡與水滴的受力情況有關(guān)，蒸汽對水滴表面有曳力.在汽輪機(jī)末級排汽中大量存在的是亞觀尺寸（0.01～2μm）的水滴，適合采用Stokes曳力公式［13］.水滴的作用力平衡方程在直角坐標(biāo)系下（x方向）為：

對于亞觀粒子，其他作用力Fx中的布朗力分量不可忽視.

對式（16）積分就得到了水滴軌道上每一個位置的顆粒速度，顆粒軌道通過下式得到：

沿著每個坐標(biāo)方向求解式（17）就得到了離散相的軌跡.

2 數(shù)值模擬方法

2.1 取樣入口的網(wǎng)格結(jié)構(gòu)

諧振腔楔形取樣端入口角度為26.5°，其物理結(jié)構(gòu)和網(wǎng)格劃分見圖2.由于計算區(qū)域具有圓周對稱性，計算模型采用諧振腔縱向剖面的上半部分，諧振腔軸線為對稱邊界，計算區(qū)域采用四邊形結(jié)構(gòu)化網(wǎng)格，網(wǎng)格數(shù)量為248 871.

圖2 楔形入口網(wǎng)格結(jié)構(gòu)Fig.2 The grid structure of wedge－shaped entrance

根據(jù)楔形取樣端的水滴流線擬合曲線函數(shù)，確定諧振腔曲線取樣端結(jié)構(gòu)，以減少取樣端對流場的擾動，其曲線函數(shù)由式（18）給出，其物理結(jié)構(gòu)和網(wǎng)格劃分見圖3，計算區(qū)域采用四邊形結(jié)構(gòu)化網(wǎng)格，網(wǎng)格數(shù)量為246 046.

圖3 曲線入口網(wǎng)格結(jié)構(gòu)Fig.3 The grid structure of curved entrance

2.2 邊界條件的處理

采用壓力遠(yuǎn)場邊界，主蒸汽壓力為5 000Pa，主蒸汽溫度為306.3K；諧振腔軸心設(shè)為對稱性邊界.

3 流場的數(shù)值分析

采用計算流體力學(xué)通用軟件Fluent進(jìn)行數(shù)值模擬，計算工況的馬赫數(shù) Ma分別取0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8、0.9、1.0、1.1、1.2和1.3，濕蒸汽水滴半徑d分別取0.2μm、0.4μm、0.6μm、0.8μm和1.0μm.

3.1 楔形取樣入口計算結(jié)果及分析

當(dāng)Ma＝0.4、d＝0.4μm時，諧振腔楔形取樣前端的模擬結(jié)果見圖4和圖5；根據(jù)式（1）和式（6）計算不同工況下的取樣誤差，結(jié)果如圖6所示.

圖4 楔形入口水滴軌跡Fig.4 The droplets＇trajectory at wedge－shaped entrance

圖5 楔形入口的壓力分布Fig.5 The pressure distribution at wedge－shaped entrance

圖6 楔形入口誤差Fig.6 The measurement error at wedge－shaped entrance

3.2 曲線取樣入口計算結(jié)果及分析

當(dāng)Ma＝0.4、d＝0.4μm時，自定義函數(shù)曲線確定的取樣前端的模擬結(jié)果見圖7和圖8.根據(jù)式（1）和式（6）計算不同工況下的取樣誤差，結(jié)果如圖9所示.

Moore等關(guān)于非等動能取樣誤差［14］的計算示意圖見圖10.由于對流場的處理方法不同，選擇與取樣進(jìn)口的距離大于最小不干擾長度截面［15］（在該截面處氣流基本達(dá)到穩(wěn)定）處的數(shù)據(jù)，根據(jù)式（1）和式（6）計算取樣誤差.

圖7 曲線入口水滴軌跡Fig.7 The droplets＇trajectory at curved entrance

圖8 曲線入口的壓力分布Fig.8 The pressure distribution at curved entrance

圖9 曲線入口誤差Fig.9 The measurement error at curved entrance

圖10 非等動能取樣誤差分析Fig.10 Analysis of the non－isokinetic sampling deviation

圖4和圖7表明：由于受到諧振腔壁面的阻擋，含小尺寸水滴的蒸汽繞過諧振腔取樣前端，使得測量值偏小.圖5和圖8表明：不同的取樣前端對流場的擾動程度不同，流場的壓力分布也不同，楔形取樣前端對流場的擾動大于由自定義函數(shù)確定的曲線取樣前端對流場的擾動.圖6和圖9表明：取樣誤差與馬赫數(shù)及水滴尺寸有關(guān)；馬赫數(shù)較小時，水滴尺寸對取樣誤差的影響很小，取樣誤差主要受馬赫數(shù)的影響；馬赫數(shù)較大時，水滴尺寸對取樣誤差的影響較大，水滴尺寸較小時，測量值偏大，水滴尺寸較大時，測量值偏小.對于楔形取樣前端，馬赫數(shù)為1.0時，取樣誤差最大.對于自定義函數(shù)確定的曲線取樣前端，當(dāng)馬赫數(shù)為0.9時，含大尺寸水滴的蒸汽取樣誤差較大；當(dāng)馬赫數(shù)為1.1時，含小尺寸水滴的蒸汽取樣誤差較大.當(dāng)馬赫數(shù)較大時，蒸汽流速較高，流動阻力相應(yīng)增大，諧振腔前后的壓差增大，諧振腔對流體的抽吸能力增加，小尺寸的水滴會被吸入諧振腔內(nèi)部，由于受連續(xù)性方程的限制，進(jìn)入諧振腔取樣前端的大尺寸水滴數(shù)量將會減少，測量值偏小.此外，當(dāng)馬赫數(shù)較大時，由于流動的不穩(wěn)定性，諧振腔取樣誤差較大.

實際電廠排汽的馬赫數(shù)一般小于0.5，排汽濕度一般為6%～8%，采用帶楔形取樣前端的微波諧振腔測量蒸汽濕度時，取樣誤差小于0.24%，誤差較小，精確度較高；采用自定義函數(shù)確定的曲線取樣前端時，微波諧振腔測量蒸汽濕度的取樣誤差小于0.16%，取樣誤差很小，精確度高，能滿足工程實際的要求.

由于諧振腔壁面的阻擋，流線在進(jìn)口附近發(fā)生彎曲，改進(jìn)傳感器取樣前端的流型能有效減小取樣偏差，使取樣比較具有代表性，測量值能比較準(zhǔn)確地反映汽輪機(jī)末級后的排汽濕度.

4 結(jié) 論

（1）諧振腔微擾法測量蒸汽濕度的關(guān)鍵部件之一是諧振腔的取樣裝置，為使測量結(jié)果能準(zhǔn)確地反映汽輪機(jī)末級的排汽濕度，要求取樣盡可能反映真實的流動情況，必須充分考慮流場的氣動性能，使諧振腔對流場產(chǎn)生最小的擾動.

（2）不同的取樣前端對流場氣動性能的影響不同，非等動能取樣誤差不同.

（3）對于不同工況，取樣誤差的影響因素不同.當(dāng)馬赫數(shù)較小時，取樣誤差受水滴尺寸的影響較小，誤差較?。划?dāng)馬赫數(shù)較大時，取樣誤差受水滴尺寸影響明顯，取樣誤差較大；當(dāng)水滴尺寸較小時，測量值偏大，水滴尺寸較大時，測量值偏小.

（4）通過優(yōu)化設(shè)計諧振腔取樣前端能有效地減小取樣誤差，使諧振腔達(dá)到較高的取樣精度.

［1］韓中合，楊昆.汽輪機(jī)中蒸汽濕度測量方法的研究現(xiàn)狀［J］.華北電力大學(xué)學(xué)報，2002，29（4）：44－47.HAN Zhonghe，YANG Kun.Reviews on wetness measurement methods of wet steam in turbine［J］.Journal of North China Electric Power University，2002，29（4）：44－47.

［2］韓中合，張淑娥，田松峰，等.汽輪機(jī)排汽濕度諧振腔微擾測量法的研究［J］.中國電機(jī)工程學(xué)報，2004，23（12）：199－202.HAN Zhonghe，ZHANG Shue，TIAN Songfeng，et al.Study on a method of wetness measurement based on resonant cavity perturbation for steam turbine exhaust［J］.Proceedings of the CSEE，2004，23（12）：199－202.

［3］王升龍，楊善讓，王建國，等.汽輪機(jī)排汽濕度的在線監(jiān)測方法及工業(yè)試驗研究［J］.中國電機(jī)工程學(xué)報，2005，25（17）：83－87.WANG Shenglong，YANG Shanrang，WANG Jianguo，et al.Study on a method of wetness measurement on line and industrial test for steam turbine exhaust［J］.Proceedings of the CSEE，2005，25（17）：83－87.

［4］黃雪峰，盛德仁，陳堅紅，等.濕蒸汽兩相流濕度測量方法研究進(jìn)展［J］.電站系統(tǒng)工程，2006，22（5）：1－4.HUANG Xuefeng，SHENG Deren，CHEN Jianhong，et al.Development study on measuring moisture of the wet steam two phase flow［J］.Power System Engineering，2006，22（5）：1－4.

［5］韓中合，錢江波，楊昆，等.諧振腔微擾法測量汽輪機(jī)排汽濕度技術(shù)的關(guān)鍵問題［J］.動力工程，2005，25（3）：387－391.HAN Zhonghe，QIAN Jiangbo，YANG Kun，et al.Key problems of steam turbine exhaust wetness measurements by the resonant cavity perturbation method［J］.Journal of Power Engineering，2005，25（3）：387－391.

［6］韓中合，楊昆，田松峰，等.在線監(jiān)測透平油中微水含量的諧振腔微擾測量法［J］.汽輪機(jī)技術(shù)，2003，45（5）：273－275.HAN Zhonghe，YANG Kun，TIAN Songfeng，et al.Study on the method of on－line measurement of water content in the turbine oil based on resonant cavity perturbation［J］.Turbine Technology，2003，45（5）：273－275.

［7］田松峰，韓中合，楊昆.流動濕蒸汽濕度諧振腔微擾法測量的實驗研究［J］.動力工程，2005，25（2）：254－257.TIAN Songfeng，HAN Zhonghe，YANG Kun.Experimental research on wetness measurement of flowing steam by resonant cavity perturbation［J］.Journal of Power Engineering，2005，25（2）：254－257.

［8］汪麗莉，蔡小舒.汽輪機(jī)濕蒸汽兩相流中水滴尺寸研究進(jìn)展［J］.上海理工大學(xué)學(xué)報，2003，25（4）：307－312.WANG Lili，CAI Xiaoshu.Progress in research on water droplet size in two phase wet steam flow in steam turbine［J］.Journal of University of Shanghai for Science and Technology，2003，25（4）：307－312.

［9］蔡頤年，王乃寧.濕蒸汽兩相流［M］.西安：西安交通大學(xué)出版社，1983：41－43.

［10］錢江波.蒸汽濕度測量方法研究及實驗系統(tǒng)設(shè)計［D］.保定：華北電力大學(xué)能源動力與機(jī)械工程學(xué)院，2004.

［11］LAALI A R，COURANT J J，KLEITZ A.Calculation and measurement of fog droplet size：comparison between nuclear and fossil fired turbine［M］.Design，Repair and Refurbishment of Steam Turbine，1991，6（13）：133－138.

［12］YOUNG J B.An equation of status for steam for turbomachinery and other flow calculations［J］.Journal of Engineering for Gas Turbines and Power，1988，110（1）：1－7.

［13］OUNIS H，AHAMADI G，MCLAUGHLIN J B.Brownian diffusion of submicrometer particles in the viscous sublayer［J］.Journal of Colloid and Interface Science，1991，143（1）：266－277.

［14］英國中央電業(yè)研究實驗室，比利時馮.卡門流體力學(xué)研究所.透平和分離器中的雙相流［M］.蔡頤年，姚福生，譯.北京：機(jī)械工業(yè)出版社，1983.

［15］韓中合.汽輪機(jī)排汽濕度微波諧振腔測量技術(shù)的研究［D］.保定：華北電力大學(xué)能源動力與機(jī)械工程學(xué)院，2005.