波形板汽水分離器流場的計算模擬

2014-08-24 00:41:10，，，

應用能源技術 2014年6期

關鍵詞：模型

，，，

(哈爾濱汽輪機廠有限責任公司，哈爾濱 150046)

0 引言

目前，汽水分離再熱器(MSR)廣泛應用波形板進行汽水分離。波形板汽水分離器可以分離蒸汽中的水分，當汽流攜帶著微小液滴在波形板構成的通道內作曲線運動，由于離心力、慣性力以及附著力的作用，水滴不能隨汽流偏轉而撞擊波形板壁面，離開波形板。波形板汽水分離器對MSR的安全性和經濟性產生很大的影響，根據相關研究[1-2]，通過MSR汽水分離對蒸汽進行除濕及加熱，可以提高汽輪機機組熱耗2.2%～2.4%，波形板汽水分離器是MSR的關鍵部件，因此對波形板汽水分離器的速度和壓力分布進行數值模擬是十分必要的。

文中對不同波形板汽水分離器進行了整體的三維數值模擬，對波形板汽水分離器的速度分布、壓降分布等情況進行了分析。分析了流經波形板汽水分離器壓力損失較大的區域，為降低壓力損失進行結構優化，從而提高機組經濟性具有指導意義。

1 計算模型

為了對波形板分離器的計算進行簡化，假設條件如下：

(1)由于重力及浮升力垂直于流動方向，不考慮重力及浮升力。

(2)熱態計算中，汽水兩相為飽和狀態，蒸汽和水滴之間無換熱，蒸汽和水滴不發生相變。

(3)定常流動，汽相不可壓。

(4)采用離散相模型，蒸汽或空氣是連續相，水滴為離散相。

(5)忽略散熱損失。

FLUENT 提供了以下湍流模型[3]：Spalart-Allmaras 模型、k-ε模型(包括標準k-ε模型、Renormalization-group (RNG)k-ε模型、帶旋流修正的Realizablek-ε模型)、k-ω模型(包括標準k-ω模型、壓力修正k-ω模型)、雷諾應力模型和大漩渦模擬模型等。值得指出的是，上述模型中沒有一個湍流模型對于所有問題都表現出絕對明顯的適用性。選擇模型時主要考慮以下幾點：流體是否可壓、流場的復雜性、精度的要求、計算機的能力和時間的限制等。

標準k-ε模型自從被Launder and Spalding提出之后，由于其具有適用范圍廣、經濟、合理的精度，所以在工程流場計算中和熱交換模擬中被廣泛應用。但對于波形板汽水分離器，由于流場存在著強烈的流線彎曲，標準k-ε模型很難滿足要求。

標準k-ε模型是個半經驗的公式，它是從實驗現象中總結出來的。研究者在確定了k-ε模型適用的范圍后，通過對其進行改造，發展了RNGk-ε模型和帶旋流修正的k-ε模型(Realizablek-ε模型)。

RNGk-ε模型來源于嚴格的統計技術。它和標準k-ε模型很相似，但是有以下改進：

(1)RNG模型在ε方程中加了一個修正項，有效的改善了精度。

(2)考慮了湍流漩渦，提高了精度。

(3)RNG理論為湍流Prandtl數提供了一個解析公式，然而標準k-ε模型使用的是用戶提供的常數。

(4)標準k-ε模型是一種高雷諾數的模型，RNG理論提供了一個考慮低雷諾數流動粘性的解析公式。通過精確處理近壁條件，RNGk-ε模型比標準k-ε模型在更廣泛的流動中具有更高的可信度和精度。

2 近壁面處理

湍流流動受壁面的影響很大。與主流方向的湍流脈動不一樣，在離壁面很近的地方，粘性耗散將抑制流體沿壁面平行方向的速度脈動，同時湍流動能的傳輸由于受到壁面的阻礙也抑制了流體垂直于壁面方向的速度脈動。所以對于僅適用于湍流核心區域(一般都遠離壁面)的k-ε模型，應該考慮怎樣使它也適用于壁面邊界層處的流動。由于近壁區域的速度分布存在較大的梯度，而且還存在著旺盛的湍流動能、耗散等的傳輸現象，因此，準確的求解近壁區域的流場是成功模擬波形板內湍流流動的關鍵條件。

近壁區域大體可以分為三層，最靠近壁面的那一層，粘性起到主導作用，呈現為層流的流態，因此稱為粘性子層；最外的一層，湍流起到主導作用，稱為完全湍流層；而在這兩者之間，粘性和湍流共同起作用，稱為過渡層。

通常，近壁面區域建模可以采用兩種方法。第一種方法并不對受粘性力影響的區域(粘性子層及過渡層)劃分網格，即避開近壁面處的粘性底層和液膜相的復雜計算，把第一個節點布置到旺盛的湍流區內。這種方法通常采用被稱為“壁面函數”的半經驗公式。壁面函數的運用能夠很好地修正湍流模型，從而解決壁面的存在對流動的影響。另一種方法中，湍流模型被修正，從而使壁面處受粘性力影響的區域也需要用網格來劃分計算，這種方法被稱為“近壁面模型”法。

對于大多數高雷諾數的流動，壁面函數法能充分節省計算資源，這種方法經濟，實用而且很精確，很受歡迎。然而對于所求解的波形板汽水分離問題，需要精確求解出壁面附近包括層流底層的流動，這樣才能夠為后面研究氣場對液膜行為的影響以及近壁液滴的行為奠定堅實的基礎。

圖1 Fluent中的近壁面處理

增強壁面處理是一種近壁面模型方法，它將增強壁面函數和雙層模型相結合來進行近壁面處理。在FLUENT 的近壁面模型中，有粘性影響的近壁面區域總是完全地當作粘性底層來處理。雙層模型方法是增強壁面處理的一個不可缺少的部分，它被用來指定近壁面單元格的湍流粘性和耗散率。此時，整個流場區域被分為一個主要受粘性影響的層流區域和一個完全的湍流區域。在建立雙層模型后，要想將其應用范圍拓展到貫穿整個近壁面區域(即包括層流邊界子層、過渡層和完全的湍流層)，需要對整個壁面區域擬定統一的壁面規則。為此，FLUENT 采用了一個由Kader提出的增強壁面函數，將線性的層流邊界子層和對數的湍流壁面規則結合起來，從而為精確求解整個湍流邊界層提供了巧妙的解決辦法。

3 網格劃分

采用gambit軟件，首先根據波形板的尺寸建立幾何模型，取兩塊波形板之間的空間為計算域，可以初步劃分網格。考慮到流道結構的復雜性，采用了非結構化的四邊形單元格、Pave算法來生成網格。在Fluent 中非結構化網格特征的兩個重要優點是：

(1)跟結構化網格相比減少創建時間。

(2)不會帶來因插入新的網格點而使其他網格計算品質變壞，這個性質對于后面利用Fluent的自適應調整功能來優化網格是非常重要的。

采用初步生成的網格進行模擬計算并不能實現良好的計算精確度和網格無關性。由于時均流和湍流的強烈相互作用，使得湍流流動的數值結果比層流流動更依賴于網格的劃分。

雖然增強壁面處理能夠拓展到精確求解湍流邊界下層的粘性子層，但是也必須要在近壁面的區域，構造出足夠細的網格，以便能夠分辨出這一薄薄的粘性子層。因此，在采用增強壁面處理時，近壁面處的網格劃分必須滿足以下條件：

(1)緊鄰壁面的單元格滿足y+≈1。

(2)為了能夠求解出近壁區域內的時均速度和湍流度，在有粘性影響的近壁面區域(Rey<200)內至少有10 個單元格。

圖2 波形板局部在自適應調整前后網格劃分的變化

以上過程可以通過Fluent軟件針對y+的網格自適應功能來調整近壁面處網格的劃分。該功能使得能夠根據數值計算的結果來修改網格疏密布置或網格走向。網格的自適應性能夠充分地利用計算資源來優化網格，可以在現有的網格上直接根據計算結果來進行網格優化，而不需要每次都重新生成網格，這樣可以大大縮短計算周期。但這也并不意味著Fluent的網格自適應功能是萬能的，仍然需要保證最初生成的網格能夠充分體現流場的幾何特征。圖4給出了針對y+采用自適應調整后，波形板局部在近壁處網格劃分變化的例子。從圖中可以看出，經過自適應調整后，在壁面上的網格劃分變細了。

5 模擬結果分析

在模擬計算中[4]，采用k-epsilon Realizable模型，近壁面處理采用 Enhanced Wall Treatment。求解器中的壓力離散格式采用“PRESTO!”，其余動量項采用“Second Order Upwind”。在求解過程中，采用Yplus最大值為2的網格自適應策略細化近壁面的網格。

針對板型1的流道仿真計算結果如圖3所示，通過流道的蒸汽壓降為880 Pa，計算初始網格數49621。

針對板型2維流道仿真計算結果如圖4所示，通過流道的蒸汽壓降為420 Pa，計算初始網格數49621。

針對板型3的流道仿真計算結果如圖5所示，通過流道的蒸汽壓降為200 Pa，計算初始網格數60485。

圖6、7、8顯示沿著流動方向和高度方向三維流道的流場示意圖。