葉柵通道內(nèi)濕蒸汽非平衡凝結(jié)流動(dòng)的數(shù)值模擬

韓中合，韓　旭，李恒凡

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

?

葉柵通道內(nèi)濕蒸汽非平衡凝結(jié)流動(dòng)的數(shù)值模擬

韓中合，韓旭，李恒凡

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

過(guò)冷度是各種凝結(jié)現(xiàn)象產(chǎn)生、發(fā)展的直接驅(qū)動(dòng)力，葉柵通道內(nèi)濕蒸汽成核過(guò)程通常集中在喉部下游很窄的區(qū)域內(nèi)，水滴數(shù)目和水滴半徑分布則受到邊界層和尾跡影響。針對(duì)葉柵通道內(nèi)跨音速非平衡凝結(jié)流動(dòng)參數(shù)分布陡峭、變化敏感的特點(diǎn)，采用具有較好激波捕獲效果的高分辨率二階TVD格式進(jìn)行離散。利用時(shí)間推進(jìn)法對(duì)控制方程進(jìn)行求解，建立了凝結(jié)流動(dòng)的數(shù)值解法，模擬與實(shí)驗(yàn)結(jié)果相吻合，驗(yàn)證了模型的準(zhǔn)確性。研究了葉柵通道內(nèi)非平衡凝結(jié)流動(dòng)的基本物理現(xiàn)象，討論了進(jìn)口過(guò)冷度對(duì)凝結(jié)特性的影響，歸納了葉柵通道內(nèi)壓力、成核率、水滴數(shù)、水滴半徑、蒸汽濕度的變化規(guī)律。研究表明：進(jìn)口過(guò)冷度對(duì)非平衡凝結(jié)流動(dòng)特性有重要影響。

濕蒸汽；葉柵通道；跨音速；兩相流；凝結(jié)；熱力學(xué)

DOI：10.11949/j.issn.0438-1157.20151614

目前，許多研究人員通過(guò)數(shù)值模擬方法對(duì)微通道內(nèi)多相流動(dòng)進(jìn)行研究［3-11］。White等［7-8］運(yùn)用二維非黏性時(shí)間推進(jìn)法對(duì)葉柵非平衡凝結(jié)流動(dòng)進(jìn)行了計(jì)算，得到了精準(zhǔn)的計(jì)算結(jié)果，通過(guò)求解汽輪機(jī)內(nèi)過(guò)冷度、膨脹率、濕度等參數(shù)，得到了Wilson 點(diǎn)與過(guò)冷度之間的關(guān)系。Bakhtar 等［9］運(yùn)用經(jīng)典成核理論，采用有限體積法對(duì)跨音速二維葉柵凝結(jié)流動(dòng)進(jìn)行了計(jì)算。李亮等［10］采用二階TVD離散方程，以濕度和水滴數(shù)目作為控制參數(shù)，對(duì)非平衡流動(dòng)進(jìn)行了數(shù)值模擬，提出兩相耦合求解方法。林智榮等［11］采用LU-SGS-GE的隱式時(shí)間推算法，對(duì)高分辨率TVD差分格式進(jìn)行改良，較為準(zhǔn)確地得到了Laval噴管中凝結(jié)參數(shù)分布情況。

過(guò)熱蒸汽在葉柵通道內(nèi)不斷膨脹，表現(xiàn)出強(qiáng)烈的多維效應(yīng)和邊界效應(yīng)，在喉部下游很窄的截面厚度內(nèi)瞬間出現(xiàn)大量凝結(jié)核。濕蒸汽兩相流遠(yuǎn)比單相流復(fù)雜得多，實(shí)驗(yàn)測(cè)量很難得到精確數(shù)據(jù)且測(cè)量難度較大，目前文獻(xiàn)中實(shí)驗(yàn)數(shù)據(jù)多為壓力和濕度，還未見(jiàn)水滴數(shù)目和成核率的報(bào)道［7-9］。本文在凝結(jié)理論和傳熱傳質(zhì)機(jī)理的基礎(chǔ)上，針對(duì)葉柵通道內(nèi)跨音速非平衡凝結(jié)流動(dòng)參數(shù)分布陡峭、變化敏感的特點(diǎn)，采用具有較好激波捕獲效果的高分辨率二階TVD格式進(jìn)行離散，采用時(shí)間推進(jìn)法對(duì)控制方程進(jìn)行求解。建立了凝結(jié)流動(dòng)的數(shù)值模型并驗(yàn)證了模型的準(zhǔn)確性。討論了進(jìn)口過(guò)冷度對(duì)凝結(jié)特性的影響，歸納了葉柵通道內(nèi)壓力、成核率、水滴數(shù)、水滴半徑、蒸汽濕度的變化規(guī)律。

1　物理模型及數(shù)值算法

1.1物理模型

在涉及濕蒸汽非平衡凝結(jié)流動(dòng)實(shí)驗(yàn)的文獻(xiàn)中，以White等［7-8］和Bakhtar等［9］的實(shí)驗(yàn)最具代表性，White等實(shí)驗(yàn)給出了各種工況下葉片表面詳細(xì)的壓比分布和紋影圖。某660 MW汽輪機(jī)低壓缸共有6級(jí)，取第5級(jí)靜葉柵作為本文研究對(duì)象。該葉柵設(shè)計(jì)出口Mach數(shù)為1.2，設(shè)計(jì)出口角為71°，葉片間距為87.59 mm，葉片弦長(zhǎng)137.51 mm，葉片安裝角45.32°。葉片表面采用絕熱固壁邊界條件u=0，v=0，?t/?n=0，葉柵進(jìn)口和出口為壓力邊界條件，葉柵通道為周期性邊界條件，見(jiàn)圖1。

圖1　葉柵邊界條件設(shè)置Fig.1　Cascade boundary conditions

本文研究的葉柵通道采用非結(jié)構(gòu)化四邊形網(wǎng)格，對(duì)葉片表面和尾緣進(jìn)行加密處理。非結(jié)構(gòu)化四邊形網(wǎng)格對(duì)形狀不規(guī)則的葉柵通道模型適應(yīng)能力較強(qiáng)，絕大部分網(wǎng)格的排列方向與流動(dòng)方向相同，在一定程度上提高了收斂速度和計(jì)算精度。

1.2凝結(jié)流動(dòng)控制方程及定解條件

由于單位質(zhì)量濕蒸汽內(nèi)的水滴數(shù)目很多，約為1024/kg，因此可將液相視為連續(xù)介質(zhì)。鑒于自發(fā)凝結(jié)產(chǎn)生的水滴尺寸較小，可認(rèn)為水滴和周圍蒸汽處于慣性平衡狀態(tài)，忽略兩相間的速度滑移，即水滴和周圍蒸汽的速度相同。將濕蒸汽成核過(guò)程和水滴生長(zhǎng)過(guò)程視為孤立的黑箱，并將氣液兩相間的傳熱傳質(zhì)過(guò)程當(dāng)作開(kāi)口系統(tǒng)進(jìn)行質(zhì)量、動(dòng)量、能量交換過(guò)程，建立自發(fā)凝結(jié)數(shù)值模型，如圖2。

圖2　自發(fā)凝結(jié)流動(dòng)數(shù)值模型Fig.2 Spontaneous condensing flow numerical model

由圖2可知，該數(shù)值模型由氣相模塊、液相模塊和成核及水滴生長(zhǎng)模塊構(gòu)成。氣相模塊將n時(shí)刻的流動(dòng)參數(shù)ρg、U、p、T提供給成核及水滴生長(zhǎng)模塊；根據(jù)無(wú)滑移假設(shè)，液相速度分布與氣相相同；同時(shí)，液相模塊將濕蒸汽參數(shù)N、r、Y傳遞給成核及水滴生長(zhǎng)模塊。成核及水滴生長(zhǎng)模塊通過(guò)計(jì)算得到成核率J、水滴生長(zhǎng)速率、水滴密度等參數(shù)，由這些參數(shù)計(jì)算出質(zhì)量凝結(jié)速率m.、凝結(jié)潛熱hfg。并將m.和hfg反饋到氣相模塊，通過(guò)氣相控制方程計(jì)算n時(shí)刻兩相間傳熱傳質(zhì)對(duì)氣相流動(dòng)的影響；將m.和J傳遞給液相模塊，液相分布方程通過(guò)成核率計(jì)算出水滴數(shù)目，由質(zhì)量凝結(jié)速率計(jì)算出濕度分布。

假設(shè)葉柵通道與外界絕熱，濕蒸汽成核和水滴生長(zhǎng)過(guò)程釋放的凝結(jié)潛熱全部加到氣相中，汽相流動(dòng)控制方程為

其中

本文水滴尺寸分布函數(shù)由零階、一階、二階矩函數(shù)組成。3個(gè)函數(shù)分別為

平均水滴半徑

假設(shè)水滴生長(zhǎng)速率與成核和生長(zhǎng)過(guò)程無(wú)關(guān)，則液相控制參數(shù)為

由式（5）得出液相分布方程為

其中

濕度計(jì)算式

氣相方程組求解時(shí)，給定葉柵通道進(jìn)口總壓、總溫。對(duì)超音速流動(dòng)，出口邊界參數(shù)由上游節(jié)點(diǎn)計(jì)算值外推得到；對(duì)亞音速流動(dòng)，應(yīng)給定背壓，其他邊界參數(shù)由上游節(jié)點(diǎn)計(jì)算值外推得到。由于本文計(jì)算進(jìn)口濕度為0，液相方程組求解時(shí)，進(jìn)、出口邊界參數(shù)全部由內(nèi)點(diǎn)計(jì)算值外推得到。當(dāng)前，求解濕蒸汽非平衡凝結(jié)流動(dòng)的高分辨率差分格式主要有：TVD、NND、失通量分裂格式等［10-12］。本文采用具有較好激波捕獲效果的高分辨率二階TVD格式對(duì)凝結(jié)流動(dòng)控制方程進(jìn)行離散。

1.3凝結(jié)過(guò)程傳質(zhì)模型

單位質(zhì)量中的經(jīng)典成核率

經(jīng)典成核理論中認(rèn)為所有水滴溫度均與氣相溫度相同，并沒(méi)有考慮相變產(chǎn)生的凝結(jié)潛熱，對(duì)于濕空氣流動(dòng)，經(jīng)典成核理論的等溫假設(shè)是合理的。但對(duì)于濕蒸汽流動(dòng)，氣液兩相溫差不可忽略，非等溫效應(yīng)對(duì)成核率有較大影響。本文采用非等溫效應(yīng)修正成核率

非等溫修正因子［13］

水分子與凝結(jié)核的撞擊頻率與水滴半徑和氣體分子平均自由程有關(guān)，引入量綱1參數(shù)Knudsen數(shù)（Kn）衡量撞擊頻率［14-15］。水滴生長(zhǎng)過(guò)程用Young低壓修正模型進(jìn)行描述。由于葉柵通道內(nèi)的兩相流為非平衡狀態(tài)，蒸發(fā)系數(shù)與凝結(jié)系數(shù)并不相等，水滴與氣相間存在非等溫效應(yīng)，提出半經(jīng)驗(yàn)修正系數(shù)υ對(duì)水滴生長(zhǎng)速率進(jìn)行修正。

低壓修正水滴生長(zhǎng)模型［16］

氣體狀態(tài)方程采用Young提出的維里型方程

B、C、D為2～4階維里系數(shù)。對(duì)于濕蒸汽計(jì)算，2階維里方程計(jì)算精度即可滿足要求。

熱導(dǎo)率計(jì)算式［17］

2　模型驗(yàn)證

為了驗(yàn)證本文數(shù)值模型的準(zhǔn)確性，對(duì)低過(guò)熱蒸汽凝結(jié)流動(dòng)進(jìn)行計(jì)算，流動(dòng)進(jìn)、出口條件在表1中給出。低過(guò)熱蒸汽進(jìn)口來(lái)流方向與軸向夾角為0°，葉片表面壓力分布實(shí)驗(yàn)數(shù)據(jù)取自文獻(xiàn)［7］。

表1　R葉柵通道進(jìn)出口參數(shù)Table 1　Boundary conditions of cascade passage

圖3　葉片表面壓比分布與實(shí)驗(yàn)值比較Fig.3　Blade surfaces pressure ratio compared with experimental value

圖3給出了工況1和工況2葉片表面壓比計(jì)算值和實(shí)驗(yàn)值。兩種工況下的凝結(jié)沖波位置、波形計(jì)算結(jié)果與實(shí)驗(yàn)值吻合較好，說(shuō)明本文模型較為準(zhǔn)確。由于兩種工況進(jìn)口壓力值僅相差0.6 kPa，進(jìn)口溫度相同，因此葉片表面壓力分布規(guī)律較為相似。兩種工況壓力面壓比分布幾乎完全相同，可見(jiàn)非平衡凝結(jié)流動(dòng)對(duì)葉柵壓力面影響不大，但吸力面凝結(jié)沖波位置及形狀并不相同。在60%軸向弦長(zhǎng)附近，兩種工況吸力面壓力出現(xiàn)跳躍，此處即為Wilson點(diǎn)。該點(diǎn)位于葉柵通道喉部下游，此處曲率半徑變化較大，汽流劇烈膨脹并產(chǎn)生大量凝結(jié)核，伴隨釋放大量凝結(jié)潛熱，引發(fā)壓力突增，兩種工況計(jì)算結(jié)果較為準(zhǔn)確地預(yù)測(cè)了Wilson點(diǎn)位置。

3　數(shù)值計(jì)算結(jié)果和分析

3.1凝結(jié)流動(dòng)基本物理現(xiàn)象

根據(jù)實(shí)驗(yàn)和計(jì)算結(jié)果可以看出，葉柵通道內(nèi)濕蒸汽非平衡凝結(jié)流動(dòng)的本質(zhì)是凝結(jié)成核和水滴生長(zhǎng)過(guò)程，伴隨著熱力不平衡和動(dòng)力不平衡過(guò)程，是一種既出現(xiàn)膨脹同時(shí)又發(fā)生凝結(jié)的多相流動(dòng)［18-19］。液相為多元彌散相，其中存在大量不同尺度的微小水滴，具有連續(xù)性介質(zhì)和分散顆粒的雙重特性［20-21］。為了得到影響非平衡凝結(jié)位置和凝結(jié)沖波強(qiáng)度的熱力學(xué)參數(shù)，對(duì)葉柵通道內(nèi)的濕蒸汽凝結(jié)流動(dòng)進(jìn)行數(shù)值計(jì)算，本文共設(shè)計(jì)了3種工況，見(jiàn)表1。圖4～圖9給出了葉柵通道內(nèi)3種工況的壓力、過(guò)冷度、Mach數(shù)、成核率、水滴半徑、濕度分布。本節(jié)以工況2為例，分析非平衡凝結(jié)流動(dòng)的基本物理現(xiàn)象。

由圖4（b）和圖6（b）可知，葉片吸力面中部60%軸向弦長(zhǎng)附近，壓力和Mach數(shù)變化較為劇烈。凝結(jié)沖波主要為葉片尾緣處的斜激波，數(shù)值模擬得到的激波位置和尾跡渦形狀也與文獻(xiàn)［7］紋影圖中的相同，驗(yàn)證了模型的準(zhǔn)確性。圖5（b）給出了葉柵通道過(guò)冷度分布，隨著通道內(nèi)蒸汽不斷膨脹，過(guò)冷度不斷增加，在通道喉部達(dá)到最大值46 K，隨后自發(fā)凝結(jié)迅速出現(xiàn)，蒸汽迅速由非平衡狀態(tài)向平衡狀態(tài)過(guò)渡，過(guò)冷度迅速下降至2 K左右。圖7（b）為成核率分布，成核現(xiàn)象受通道膨脹率影響較大，成核主要集中在氣流迅速擴(kuò)張區(qū)域，峰值出現(xiàn)在壓力面尾緣，可達(dá)1023kg-1·s-1。成核區(qū)域剛好與壓比跳躍區(qū)域重合，說(shuō)明成核是導(dǎo)致凝結(jié)沖波的主要原因。

圖4　壓力分布Fig. 4　Pressure distribution/Pa

圖5　過(guò)冷度分布Fig. 5　Subcooled degree distribution/K

圖6　Mach數(shù)分布Fig. 6　Mach number distribution

圖7　成核率分布Fig. 7　Nucleation rate distribution［lg（nucleation rate/kg-1·s-1）］

圖8（b）和圖9（b）為水滴平均半徑和濕度分布。成核區(qū)域的水滴半徑很小，基本在13 nm以下，沿著葉柵通道，成核區(qū)后的水滴存在明顯生長(zhǎng)現(xiàn)象，濕度迅速增加至4%左右，水滴半徑增加到23 nm左右。雖然葉片尾緣處有大量水滴，但由于尾跡過(guò)冷度較低，水滴生長(zhǎng)速度很小，因此葉片尾跡區(qū)域的水滴半徑小于2 nm，濕度小于1%，與文獻(xiàn)［19］中的實(shí)驗(yàn)值相近。工況2葉片尾跡區(qū)域的水滴數(shù)目約為2.5×1018kg-1，主流區(qū)的水滴數(shù)目約為7.9× 1016kg-1。式（7）給出了濕度與水滴半徑、水滴數(shù)目的關(guān)系，雖然尾跡區(qū)域的水滴數(shù)目是主流區(qū)的30多倍，但主流區(qū)水滴半徑是尾跡區(qū)域的十幾倍，因此主流區(qū)濕度比尾跡大得多。

圖8　水滴半徑分布Fig. 8　Droplet radius distribution/nm

圖9　濕度分布Fig. 9　Wetness distribution

3.2凝結(jié)流動(dòng)變壓特性

工況1和工況2的進(jìn)口溫度均為354 K，進(jìn)口壓力值相差0.6 kPa，以這兩種工況為例，研究進(jìn)口壓力對(duì)凝結(jié)流動(dòng)的影響。

由于入口壓力對(duì)溫度場(chǎng)分布影響不大，在凝結(jié)發(fā)生前，蒸汽流動(dòng)狀態(tài)與等熵流動(dòng)相似。對(duì)比圖4（a）和圖4（b），工況1的壓力變化程度更大，凝結(jié)發(fā)生得更快，工況1的成核區(qū)域略小于工況2，但工況1的凝結(jié)沖波比工況2更劇烈。工況1的成核率峰值接近1.7×1024kg-1·s-1，工況2約為5.36×1023kg-1·s-1，說(shuō)明當(dāng)進(jìn)口溫度相同時(shí)，進(jìn)口壓力增加，凝結(jié)區(qū)域面積增加，成核率峰值降低。對(duì)比圖6（a）和圖6（b），兩種工況在葉柵通道喉部上游的膨脹過(guò)程基本相同，但喉部下游區(qū)域的Mach數(shù)分布差異明顯，工況1的最大Mach數(shù)接近1.357，工況2的最大Mach數(shù)為1.287。由于兩種工況的進(jìn)口過(guò)冷度差別不大，因此成核率和過(guò)冷度分布規(guī)律相似。成核率增大勢(shì)必影響水滴數(shù)目、水滴半徑和濕度分布，工況1的主流區(qū)水滴半徑約為14.6 nm，工況2為23 nm；工況1的最大濕度為4.76%，工況2為4.14%；工況1主流區(qū)水滴數(shù)目約為4×1017kg-1，工況2主流區(qū)的水滴數(shù)目約為7.9×1016kg-1。工況2主流區(qū)水滴半徑約為工況1的1.6倍，而工況1的水滴數(shù)目是工況2的5.06倍。

研究表明：隨著葉柵通道進(jìn)口壓力的增加，凝結(jié)區(qū)域面積增加，成核率降低，由此導(dǎo)致喉部下游水滴數(shù)量減少，水滴平均半徑增加。研究結(jié)果與文獻(xiàn)［20］中Laval噴管凝結(jié)流動(dòng)變壓特性相同。

3.3凝結(jié)流動(dòng)進(jìn)口過(guò)冷度特性

工況2和工況3的進(jìn)口壓力相同，進(jìn)口溫度分別為354 K和357.5 K，對(duì)應(yīng)的進(jìn)口過(guò)冷度分別為-4.46 K和-7.96 K，以這兩種工況為例，研究進(jìn)口過(guò)冷度對(duì)凝結(jié)流動(dòng)的影響。

在凝結(jié)沖波出現(xiàn)之前，兩種工況的壓力分布基本相同；凝結(jié)發(fā)生后，工況2的壓力回升程度比工況3高；工況3的過(guò)冷度變化梯度較小，但變化區(qū)域較大。喉部過(guò)冷度是影響成核的主要因素，因此工況3的成核區(qū)域面積大于工況2。工況2的成核率峰值為5.36×1023kg-1·s-1，工況3的成核率峰值為7.94×1023kg-1·s-1，因此工況3釋放的凝結(jié)潛熱更多。對(duì)比圖6（b）和圖6（c），兩種工況的Mach數(shù)分布基本相同，可見(jiàn)進(jìn)口過(guò)冷度對(duì)Mach數(shù)分布影響不大。入口過(guò)冷度降低后，雖然最大水滴半徑由25.9 nm增加到50.6 nm，但主流區(qū)平均水滴半徑減小，對(duì)比圖9（b）和圖9（c）可知，工況3的出口濕度比工況2降低0.3%。

研究表明：進(jìn)口過(guò)冷度是影響成核率的主要因素，不同工況下的水滴數(shù)目、水滴半徑、濕度差異均由成核率導(dǎo)致。進(jìn)口過(guò)冷度越小，葉柵通道中的成核率越大，凝結(jié)過(guò)程更為劇烈，水滴數(shù)目越多。

4　結(jié)　論

吸力面成核率對(duì)主流區(qū)液相參數(shù)分布影響較大，優(yōu)化吸力面膨脹率設(shè)計(jì)是控制葉柵通道中液相參數(shù)分布的最佳方法。針對(duì)葉柵通道內(nèi)跨音速非平衡凝結(jié)流動(dòng)參數(shù)分布陡峭、變化敏感的特點(diǎn)，采用具有較好激波捕獲效果的高分辨率二階TVD格式進(jìn)行離散，采用時(shí)間推進(jìn)法對(duì)控制方程進(jìn)行求解。討論了進(jìn)口過(guò)冷度對(duì)凝結(jié)特性的影響，歸納了葉柵通道內(nèi)壓力、成核率、水滴數(shù)、水滴半徑、蒸汽濕度的變化規(guī)律。研究表明：

（1）成核率受通道膨脹率影響，其峰值出現(xiàn)在壓力面尾緣，是導(dǎo)致凝結(jié)沖波的主要原因；

（2）隨著葉柵通道進(jìn)口壓力增加，凝結(jié)區(qū)域面積增加，但成核率降低，水滴數(shù)量減少，水滴平均半徑增加；

（3）進(jìn)口過(guò)冷度是影響成核率的主要因素，進(jìn)口過(guò)冷度越小，通道內(nèi)成核率峰值越大，凝結(jié)過(guò)程更為劇烈，水滴數(shù)目越多，但水滴平均半徑越小。

符號(hào)說(shuō)明

E ——能量密度，kJ·m-3

hfg——凝結(jié)潛熱，kJ·kg-1

ht——蒸汽焓值，kJ·kg-1

J ——單位體積成核率， （m3·s）-1

Kn ——Knudsen數(shù)

M ——Mach數(shù)

m ——單分子的質(zhì)量，kg

m. ——質(zhì)量凝結(jié)速率，kg·s-1

N ——單位質(zhì)量蒸汽水滴數(shù)

Prg——Prandtl數(shù)

p ——蒸汽實(shí)際壓力，Pa

qc——凝結(jié)系數(shù)

R——理想氣體常數(shù)，J·（mol·K）-1

r——水滴半徑，nm

S——過(guò)飽和度

T——溫度，K

ΔT——過(guò)冷度，K

u——軸向速度，m·s-1

v——徑向速度，m·s-1

Y——濕度

γ——蒸汽比熱容比

θ——量綱1表面張力

λg——蒸汽的熱導(dǎo)率，W·（m·K）-1

ρ——密度，kg·m-3ssσ——水滴表面張力，N·m-1

τ——微元控制體

υ ——水滴生長(zhǎng)半經(jīng)驗(yàn)修正系數(shù)

下角標(biāo)

c——水滴臨界半徑

g——?dú)庀鄥?shù)

i——第i個(gè)水滴

l——液相參數(shù)

s——飽和狀態(tài)

References

［1］ MOORE M J， SIEVERDING C H. Two-phase Steam Flow in Turbines and Separators［M］. London： Hemisphere Publishing Corporation， 1976.

［2］ BAUMANN K. Some recent developments in large steam turbine practice ［J］. J. Inst. Elec. Engrs.， 1921， 59（302）： 565-623.

［3］ 沈勝?gòu)?qiáng)， 楊勇， 張琨， 等. 水蒸氣超音速非平衡流動(dòng)的變壓凝結(jié)特性 ［J］. 化工學(xué)報(bào)， 2010， 61（4）： 820-824. SHEN S Q， YANG Y， ZHANG K， et al. Condensation characteristics of supersonic non-equilibrium steam flow under different pressures［J］. CIESC Journal， 2010， 61（4）： 820-824.

［4］ 趙家權(quán)， 劉培啟， 趙文靜， 等. 激波管中非定常凝結(jié)現(xiàn)象的數(shù)值分析 ［J］. 化工學(xué)報(bào)， 2012， 63（4）： 1050-1055. DOI：10.3969/j.issn. 0438-1157. 2012.04.009 ZHAO J Q， LIU P Q， ZHAO W J， et al. Numerical analysis of unsteady condensation during expansion in shock tube ［J］. CIESC Journal， 2012， 63（4）： 1050-1055. DOI：10.3969/j.issn.0438-1157. 2012.04.009.

［5］ 楊勇， 沈勝?gòu)?qiáng)， 董國(guó)海， 等. 水蒸氣超音速流動(dòng)中非平衡相變的溫度特性 ［J］. 化工學(xué)報(bào)， 2012， 63（2）： 401-407. DOI：10.3969/ j.issn.0438-1157.2012.02.010. YANG Y， SHEN S Q， DONG G H， et al. Effect of temperature on non-equilibrium phase change in transonic steam flow ［J］. CIESC Journal， 2012， 63（2）： 401-407. DOI：10.3969/j.issn.0438-1157. 2012.02.010.

［6］ 宿吉強(qiáng)， 孫中寧， 高力.中等過(guò)冷度下含不凝性氣體蒸汽冷凝傳熱特性 ［J］. 化工學(xué)報(bào)， 2014， 65（10）： 3884-3890. DOI： 10.3969/j.issn. 0438-1157.2014.10.019. SU J Q， SUN Z N， GAO L. Analysis of experiments for steamcondensation in presence of non-condensable gases with moderate wall subcooling ［J］. CIESC Journal， 2014， 65（10）： 3884-3890. DOI：10.3969/j.issn.0438-1157.2014.10.019.

［7］ WHITE A J， YOUNG J B， WALTERS P T. Experimental validation of condensing flow theory for a stationary cascade of steam turbine blades ［J］. Philos. Trans. R. Soc. Lond. ， 1996， 354（1704）： 59-88.

［8］ HUANG L， YOUNG J B. An analytical solution for the Wilson point in homogeneously nucleating flows ［J］. Proc. R. Soc.， 1996， A452：1459-1473.

［9］ BAKHTAR F， EBRAHIMI M， WEBB R A. On the performance of a cascade of turbine rotor tip setion blading in nucleating steam （Ⅰ）：Surface pressure distributions ［J］. Proc. Inst. Mech. Eng. Part C： J. Mech. Eng. Sci.， 1995， 209： 115-124.

［10］ 李瑜， 李亮， 鐘剛云， 等. 末級(jí)透平對(duì)低壓缸氣動(dòng)和凝結(jié)特性的影響 ［J］. 西安交通大學(xué)學(xué)報(bào)， 2012， 46（7）： 16-20. LI Y，LI L， ZHONG G Y ， et al. Influence of last stage turbine on aerodynamic performance and condensation characteristic of low pressure cylinder ［J］. Journal of Xi'an Jiaotong University， 2012，46（7）： 16-20.

［11］ 林智榮， 袁新. 自發(fā)凝結(jié)流動(dòng)數(shù)值模擬方法及其在Laval噴管中的應(yīng)用 ［J］. 工程熱物理學(xué)報(bào)， 2006， 27（1）： 42-44. LIN Z R， YUAN X. A numerical method for spontaneous condensing flow and its application on Laval nozzle ［J］. Journal of Engineering Thermophysics， 2006， 27（1）： 42-44.

［12］ GERBER A G， SIGG R， V?LKER L， et al. Predictions of non-equilibrium phase transition in a model low-pressure steam turbine ［J］. Proc. IMech E， Part A： J. Power and Energy， 2007，221（6）： 825-835.

［13］ KANTROWITZ A. Nucleation in very rapid vapor expansions ［J］. J. Chem. Phys.， 1951， 19： 1097-1100.

［14］ 安連鎖， 王智， 韓中合. 汽輪機(jī)葉柵內(nèi)濕蒸汽兩相凝結(jié)流動(dòng)的數(shù)值研究 ［J］. 中國(guó)電機(jī)工程學(xué)報(bào)， 2009， 29（11）： 70-74. DOI：10.13334/ j.0258-8013.pcsee.2009.11.013. AN L S， WANG Z， HAN Z H. Numerical simulation of wet Steam two-phase flow in turbine cascade ［J］. Proceedings of the CSEE， 2009，29（11）： 70-74. DOI：10.13334/j.0258-8013.pcsee.2009.11.013.

［15］ 韓中合， 陳柏旺， 劉剛， 等.濕蒸汽兩相凝結(jié)流動(dòng)中水滴生長(zhǎng)模型研究 ［J］.中國(guó)電機(jī)工程學(xué)報(bào)， 2011， 31（29）： 79-84. DOI：10.13334/ j.0258-8013.pcsee.2011.29.001. HAN Z H， CHEN B W， LIU G， et al. Droplets growth model in wet steam two-phase condensation flow ［J］. Proceedings of the CSEE，2011， 31（29）： 79-84. DOI：10.13334/j.0258-8013.pcsee.2011.29.001.

［16］ YOUNG J B. Two-dimensional non-equilibrium wet-steam calculations for nozzles and turbine cascades ［J］. Trans. Am. Soc. Mech. Engrs.： J. Turbomachinery， 1992， 114： 569-579.

［17］ 于新峰. 基于雙流體模型的濕蒸汽兩相流動(dòng)數(shù)值模擬［D］.哈爾濱：哈爾濱工業(yè)大學(xué)， 2011. YU X F. A numerical simulation of wet steam flow based on dual fluid model［D］. Harbin ：Harbin Institute of Technology， 2011.

［18］ PATEL Y， PATEL G， TURUNEN-SAARESTI T. The effect of turbulence and real gas models on the two phase spontaneously condensing flows in nozzle［C］//Proc. ASME Turbo Expo，GT2013-94778. San Antonio， USA， 2013： 1-8.

［19］ BOLN D， KERPICCI H， REN J， et al. Homogeneous and heterogeneous nucleation in a nozzle guide vane of a LP-steam turbine［C］//4th European Conference on Turbomachinery. Firenze，Italy， 2001.

［20］ 韓中合， 韓旭， 李鵬. 濕蒸汽非平衡凝結(jié)流動(dòng)的熱力學(xué)特性 ［J］.化工學(xué)報(bào)， 2015， 66（11）： 4312-4319. DOI：10.11949/j.issn.0438-1157.20150727. HAN Z H， HAN X， LI P. Effect of thermodynamic properties on wet steam non-equilibrium condensing flow ［J］. CIESC Journal， 2015，66（11）： 4312-4319. DOI：10.11949/j.issn.0438- 1157.20150727.

［21］ 張冬陽(yáng)， 劉建軍， 蔣洪德. 三維濕蒸汽流動(dòng)快速準(zhǔn)確數(shù)值模擬方法及應(yīng)用 ［J］. 工程熱物理學(xué)報(bào)， 2003， 24（2）： 262-264. ZHANG D Y， LIU J J， JIANG H D. The approach and application of fast and accurate numerical simulation on 3D wet steam flow ［J］. Journal of Engineering Thermophysics， 2003， 24（2）： 262-264.

Numerical simulation of wet steam non-equilibrium condensing flow in cascade passage

HAN Zhonghe， HAN Xu， LI Hengfan
（Key Laboratory of Condition Monitoring and Control for Power Plant Equipment， North China Electric Power University，Ministry of Education， Baoding 071003， Hebei， China）

Subcooled degree is the direct drive potential of condensing phenomena happening and developing. The nucleation process of cascade passage sharply occurs in a narrow flow section. Water droplets distribution is affected by boundary layers and wakes. Due to the condensing parameters distribution steep and sensitive， high resolution second-order TVD scheme and time-marching technique were adopted in numerical algorithm. Using numerical algorithm， the condensing parameter distributions was calculated and compared with measured data to verify the accuracy of the numerical algorithm. Basic physical phenomena of non-equilibrium condensing flow and the influence of inlet subcooled degree on the homogenous nucleation were examined. Varying pattern of pressure， nucleation rate， droplets number， droplets radius and wetness were investigated. It was found that the inlet subcooled degree had a significant influence on non-equilibrium phase change in transonic condensing flow. Key words：wet steam； cascade passage； transonic； two-phase flow； condensation； thermodynamics

引　言

濕蒸汽非平衡凝結(jié)流動(dòng)廣泛存在于航空、化工、發(fā)電等工程領(lǐng)域，在動(dòng)力推進(jìn)、氣液分離、動(dòng)力機(jī)械節(jié)能等方面有著較高的應(yīng)用價(jià)值。電站汽輪機(jī)通流部分在長(zhǎng)期水蝕作用下，出現(xiàn)流動(dòng)效率降低，

date： 2015-10-26.

HAN Zhonghe， han_zhonghe@163.com

supported by the National Natural Science Foundation of China （51576066） and the Fundamental Research Funds for the Central Universities （2015MS107）.機(jī)組經(jīng)濟(jì)性變差等問(wèn)題，葉片長(zhǎng)期工作在濕蒸汽工況下，還可能出現(xiàn)斷裂［1］。Baumann［2］指出：汽輪機(jī)低壓級(jí)1%的濕度增加，將導(dǎo)致1%的效率下降。因此，研究濕蒸汽兩相凝結(jié)流動(dòng)的傳熱傳質(zhì)機(jī)理，成為提高發(fā)電效率的重要課題。

O 354

A

0438—1157（2016）05—1806—08

2015-10-26收到初稿，2016-01-18收到修改稿。

聯(lián)系人及第一作者：韓中合（1964—），男，教授。

國(guó)家自然科學(xué)基金項(xiàng)目（51576066）；中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目（2015MS107）。