濕式汽輪機葉柵熱噴射最佳位置的綜合研究

王慶生

（哈爾濱汽輪機廠有限責任公司，黑龍江哈爾濱 150046）

0 引言

在許多工業設備中會發生凝結現象，如超音速噴管、熱壓機、汽輪機葉片、噴射器和冷凝器等。在低壓汽輪機的末級膨脹過程中，過熱蒸汽通過飽和線進入兩相流區，并產生較大的流動，形成許多細小的液滴。渦輪內液相的存在，一方面會造成葉片的熱力學損失、氣動損失、葉片侵蝕、侵蝕損傷和熱效率降低，另一方面濕度增加會降低濕式汽輪機效率。

為更好地理解蒸汽凝結流動，研究了近幾年的大量文獻以及各種式蒸汽流動數值模型的求解，總結實驗數據驗證解析方法。例如，Sharifi 和Mazzelli 團隊基于非平衡冷凝的影響因素來研究噴射器以及熱壓機的流動參數，Dykas 團隊對單流體模型與雙流體模型進行研究、比較，Wroblewski 采用矩量法計算雙流體液滴半徑模型，山本研究組基于進口位置濕度來研究渦輪葉片內濕蒸汽流動對壓縮機的影響。

研究人員的重點是如何減少工程系統中的不可逆性，以及如何將傳熱和粘性應力降低到最小狀態、獲得最佳性能。因此，一致認為熵產是衡量的標準工程系統設計中的不可逆性的重要指標。顯然，有一些方法可以減少熱力學損失和熵產，最近進行了大量的研究。Gerber 和Kermani 團隊提出濕蒸汽流動的熱力學對氣動損失的影響，Lakzian 團隊在聚結過程對熵的影響方向進行了深入研究，也對發散噴嘴中冷凝濕蒸汽流的產生規律進行了模擬，并研究了濕蒸汽中汽液相之間滑移的影響，以及一維超音速噴管內的流動。Mahpeykar 研究了體積冷卻和加熱對兩相參數的影響流動，冷凝激波定位和熵產，利用一維收斂—發散噴嘴解析。Ahmadpour 等人調查了入口過熱和體積冷卻的影響作為兩種實際的控制選擇濕蒸汽流中的凝結。他們證明這兩種方法都能有效減少濕損失和提高熱壓機和汽輪機的性能。Lakzian 團隊研究表明，可以通過加熱渦輪容積來減少葉片的濕量，這是因為容積加熱可以減少熵的生成。Han 團隊設計了定子葉柵端壁柵欄，并安裝在通道中，探索了不同加熱強度和不同放置位置對級聯性能的影響。Boroomand 和Mirhoseini 注入熱將蒸汽在會聚—發散型噴嘴中轉化為濕蒸汽流動，研究其對兩相流動參數的影響。通過遺傳算法，他們提出了最優的HSI（Latent Heat of Vaporization，熱蒸汽注入）量，結果表明，在收斂—發散型噴嘴中，濕度和液滴半徑都會影響恒定生成指數。

如何控制冷凝現象是影響汽輪機葉片性能的關鍵問題。HSI進入葉片通道以及降低出口液質量分數，可以控制凝結現象。上述文獻提供了研究基礎，但并沒有能詳細介紹氣體噴射的最佳位置。在已知壓強分布、馬赫數的條件下，靠近吸力側的液體質量分數與靠近壓力側的液體質量分數需要考慮到不同位置的吸力和壓力，以選擇最佳位置。為此，采用EEHIC（Energy Efficient Home Improvement Credit）方法注射位置與未注射情況下進行比較，采用湍流模型模擬渦輪葉片內的粘性流動，來選擇得出最佳位置。

1 數值方法

1.1 控制方程

在本研究中，采用單流體模型求解，建立了液滴的質量、動量和能量守恒方程、數值方程和濕度方程。該模型同時計算了液相和氣相求解了液滴和蒸汽流體相結合的控制方程。由于液滴半徑小或阻力小，兩相之間滑移速度被忽略。

1.2 邊界條件

總壓力和溫度、液滴數量和液體質量分數，由于超音速流動的所有熱物理性質為葉片亞音速入口邊界和噴射槽指定。

1.3 侵蝕率

液滴導致的侵蝕被認為是蒸汽輪機設計中最嚴重的問題之一，因為它會造成機械損壞。Lee 等人提出了一個模型計算侵蝕速率，該速率基于液滴流的碰撞速率、碰撞速度液滴、液滴尺寸和材料硬度：在水滴成核和生長過程中，熱力學和動力學蒸汽的性質也會改變，導致蒸汽的工作能力降低。這些變化稱為濕度損失，在設計蒸汽輪機過程中對濕度損失的計算非常重要。當成核現象發生時，蒸汽液滴表面的分子凝結、液滴生長，導致冷凝損失。冷凝過程是一個不可逆過程。

1.4 數值模擬

控制方程采用基于密度的有限體積法求解，采用Roe 方法計算對流通量，最后采用二階精度的迎風方法進行空間離散。Bakhtar 等人目標的實驗數據用于驗證濕空氣的數值解汽輪機葉片中的蒸汽流量，列出了幾何規格和該葉片的邊界條件。

（1）入口的滯止壓力和滯止溫度。葉片的壓力分別為172 kPa和380.66 kPa，葉片入口的流量為亞音速、出口處為超音速。

（2）3 組四邊形計算網格：精細（30 840 個單元），中等（14 440 個細胞）和粗（7432 個細胞）通過網格收斂進行計算和分析索引以獲得最佳計算網格。

（3）關于濕蒸汽冷凝流量的規范，取質量流率為參數。

（4）噴嘴A 和B 的幾何形狀，具有指定的入口邊界條件：噴嘴A 的停滯壓力為25 kPa，停滯溫度為354.6 K，噴嘴B 的停滯溫度為357.6 K。

為了獲得最佳計算網格，使用GCI（Grid Convergence Index，網格收斂指數）方法確定適當的網格分辨率。噴嘴A 采用具有9200 個單元的四邊形網格、噴嘴B 使用15 000 個單元的四邊形網格，沿噴嘴中心線，從模擬和實驗摩爾噴嘴A 和B 獲得數據。

2 結果和討論

2.1 熱蒸汽注入（HSI）

本研究的主要目的是將熱蒸汽注入適當位置以消除濕度和液滴半徑的影響，并減少產生的熵。自從bakhtar 的渦輪葉片具有實驗半徑尺寸，該葉片用于研究HSI 在吸入側和壓力側的不同位置執行。對于每個注射位置完成中心線上的壓力分布后，通過選擇最佳網格檢查注射部位的影響。

2.2 評估HSI 位置的影響

當濕蒸汽在渦輪葉片上流動時，當馬赫數上升時橫截面、壓力和溫度在葉片收斂部分減小。之后，盡管壓力和溫度馬赫數和流動橫截面增加，但仍繼續減小。此外，過冷的程度是非平衡的標準蒸汽增加，因此形成液相以恢復平衡。

首先，增加液滴，然后是蒸汽溫度。根據瑞利流特性，向超音速流提供熱量會增加壓力并降低馬赫數，之間液相的影響是每個渦輪級的效率降低約1%（在1%濕度下）。在本研究中，由于侵蝕和濕度的存在，使用HSI 來降低濕度，考慮到壓力、溫度和HSI 速率的限制，選擇適當的HSI 位置對EEHIC 有很大影響。

本項目在不同位置評估HSI 位置對吸力和壓力側，噴射角為90°，噴射槽寬度為0.09 E#c、注入，注入蒸汽溫度為500 K，壓力側噴射壓力為170 kPa，壓力側為160 kPa。在吸入側，考慮了5 個注射位置，以檢查影響注射位置。此外，1 號和2 號的注射位置位于葉片處收斂部分（亞音速部分），3 號和4 號位于葉片發散部分（超音速部分5 號是在吸入側發生第一次空氣動力沖擊之后靠近葉片出口，在壓力側6 號～10 號的注射位置位于葉片收斂部分（亞音速）。熱量和質量被添加到HSI 的流動中。盡管使用了簡單的假設流動，二維數值解用于評估HSI 的影響，因為兩相流條件和渦輪葉片幾何形狀的復雜性。噴射位置對吸入側和壓力側的影響葉片中心線的壓力和溫度比，并與沒有的情況進行比較。注射葉片收斂部分吸入側和壓力側的HSI（注入1 和在吸入側注入2 次，在壓力側注入6～10 次）噴射位置處葉片中心線的溫度。此外，超音速部分（3 號、4 號和5 號噴射）顯示了不同噴射中葉片中心線的馬赫數位置。葉片收斂部分的HSI在領先時降低馬赫數葉片發散部分的馬赫數突然變化。氣流在喉部不能變成音速，馬赫數在喉部之后和葉片的發散部分增加。原因是HSI 為氣流創造了一個理想的喉部，導致噴射中心線的馬赫數增加（3 號和4 號噴嘴）。在幾乎位于葉片末端的噴射5中，顯示出氣流具有類似于葉片收斂部分無噴射的情況，導致擴散部分的壓力和溫度增加，在注射位置附近馬赫數減少。

比較不同位置的HSI 對成核圖、液體質量分數和葉片中心線處的平均液滴半徑。HSI 略微延遲葉片中心線處的成核開始，在成核過程中在吸入側，注射1、注射2 和注射3 的最大成核值較大，而4 號注入的最大成核值小于不注射的情況。在注射5中，沒有顯著變化與不注射的情況相比。在成核圖的壓力側，HSI 略微增加了最大成核量。然而，吸入側的HSI 減少了葉片中心線處的液體質量分數，注射4 中的液體質量分數為0。

3 結論

在驗證了渦輪葉片的數值解之后，本項目研究了HSI 對兩相流參數的影響，以選擇最佳注入位置。多次注射選擇吸入側和壓力側的位置，是因為在吸入側和壓力側附近流動。此外，考慮采用EEHIC 方法選擇最佳注射位置。然后，研究了HSI 對每次注射的5 個參數的影響，未來將研究壓力、溫度、槽角的影響。結果表明，HSI 降低了葉片中的入口質量流率壓力、溫度和收斂部分的馬赫數，以及突然上升壓力、溫度和發散部分的馬赫數。建議熱蒸汽應注入沖擊前的一個位置且溫度最低流。在本研究中，4 號注入位置減少了冷凝損失和侵蝕與未注射的情況相比，分別提高了81%和99%以上。