動力渦輪葉輪流場分析

祁博武，王 強，黃秋萍（.中北大學，太原03005；.東南汽車，福州3509）

動力渦輪葉輪流場分析

祁博武1，王 強1，黃秋萍2
（1.中北大學，太原030051；2.東南汽車，福州350119）

為了研究動力渦輪關鍵部件的流場，借助GT-Power一維模型獲取的概括性數據作為建模的邊界條件，建立了動力渦輪關鍵部位的模型，劃分網格，利用CFX分析得出子午面，葉片不同葉高，氣動載荷，熵分布分析的云圖與分析結果，分析了流場內氣流分離和尾流現象的產生機理，有助于改進動力渦輪整體氣動造型和局部氣動造型，為動力渦輪的改進和優化提供了理論依據。

動力渦輪葉輪流場分析

1　　引言

動力渦輪作為渦輪復合柴油機的關鍵部件，對整個聯合循環系統的性能有重要影響。同時，動力渦輪作為能量回收裝置，將廢氣的能量轉換為機械能，而廢氣能量的利用率與動力渦輪內的氣體流動有關，因而研究動力渦輪關鍵部件內的氣體流動具有重要的意義。

2　　建模

根據動力渦輪的相關幾何參數，運用BLADE-GEN軟件繪制出動力渦輪葉輪的三維幾何結構模型，如圖1（a）所示。由于渦輪葉輪具有周期對稱性，簡化模型，減少網格的計算量，節省計算時間，本文任意選取渦輪的一個葉片建立單通道流場計算域，即兩側由周期邊界界定的一個圍繞葉片的流場，適當地延伸進出口的流體域，如圖1 （b）所示。

3　　網格劃分

將建立的葉片單通道流體域模型導入到ICEM軟件中劃分網格，選用六面體結構化網格，結構化網格具有比較高的網格質量，對復雜邊界能夠保持較好的貼體性，如圖2所示。

4　　邊界條件

將網格模型導入CFX軟件中，使用軟件自帶的Turbomachinery旋轉機械模式，選擇機械類型為徑流渦輪，設定相應的邊界條件。在動力渦輪葉片流場計算中，邊界條件有進出口邊界條件、周期性邊界條件、壁面邊界條件等。據GT-Power一維模型確定的動力渦輪性能參數模擬計算邊界條件見表1。

圖1　動力渦輪模型

圖2　　網格模型和局部加密

表1　　動力渦輪性能參數

5　計算模型

本文將對動力渦輪葉輪內氣體的紊流黏性流動進行定常計算。該流動過程遵循三大守恒方程，通過增加一個合適的湍流模型使方程組封閉，求解雷諾時均Navier-Stokes控制方程。計算模型見表2。求解控制設置，殘差收斂準則設為1×10-6，計算時間步長0.02 s。當全流場計算殘差曲線到趨于1×10-6，且平直，認為計算收斂。

表2　　計算模型及求解器設置

6　　計算結果

6.1子午面流場分析

圖3至圖5分別渦輪葉輪流道子午面的溫度、壓力和馬赫數的分布云圖。從圖3中可以看到，氣流由進口到出口溫度變化趨勢大致沿著氣流方向溫度由高到低的分布。流道進口處溫度最高，與設定的邊界值條件吻合，流場進口前端至葉片前緣部分的溫度變化小，此后溫度逐漸降低，呈梯度性變化。由圖4中可以看出，壓力的變化趨勢大致是沿著氣流方向壓力由高到低的分布，與溫度分布趨勢相類似，壓力最大值出現在流場的進口前端、葉片前緣以及葉片根部，而后壓力逐漸降低，呈梯度性變化。從圖5中可以看到，馬赫數分布沒有呈現出從大到小的階梯變化，相反馬赫數的變化趨勢有些復雜，隨著流場方向整體呈現先增大再減小的趨勢。

圖3　　子午面溫度分布云圖

6.2葉高截面流場分析

圖6至圖8列出了20%、50%和80%葉高截面處溫度、壓力和馬赫數的分布云圖。

圖4　子午面壓力分布云圖

圖5　子午面馬赫數分布云圖

圖6　不同葉高截面溫度分布

圖7　不同葉高截面壓力分布

圖8　不同葉高截面馬赫數分布

從圖6中可看出從進口到出口溫度的大致變化趨勢是先逐漸降低至趨于平穩。在葉片尾緣處出現溫度極大值，在靠近尾緣的吸力面上，則出現局部壓力過低。從圖7可看到壓力與溫度近似成正比例關系，即壓力分布與溫度分布的大體趨勢相似，壓力從葉片前緣至尾緣逐漸減小，在靠近尾緣的吸力面處壓力出現局部極小值。從圖8中可看出馬赫數分布的大體趨勢是先增大后減小。在流場進口處馬赫數最低，當高溫氣體從噴嘴出口進入渦輪葉片中，在渦輪中膨脹加速，在葉片的壓力面與吸力面上氣體流速均增加，溫度壓力下降，但是在葉片的壓力面與吸力面上的具體分布則不同，在吸力面馬赫數緩慢升高緩慢下降；而在壓力面上，馬赫數上升和下降的趨勢較快。

6.3氣動載荷與熵分布分析

從圖9的葉片氣動載荷分布可以看出，葉片兩側的氣動載荷變化不相同，一側壓力載荷沿著氣流方向從高到低逐漸減小，另一側壓力載荷波動性大，先減小后增大再較小。由圖10的熵增分布看到，從進口到出口，氣體的熵一直在增加。因為，理想值是基于渦輪內的氣體運動是理想狀態設定，不考慮葉片通道中流體摩擦損失、湍流影響、與外界介質熱交換、熱輻射等。在三維仿真的過程中，模擬的是氣體流動的實際過程，考慮到氣體進入葉輪時的流入損失、湍流影響、氣體在通道內的轉彎損失、氣體在通道內的摩擦損失等。

圖9　葉片氣動載荷分布圖

7　　結論

（1）動力渦輪葉片前緣和尾緣處能量損失較大，而等熵效率值比較低，可能是葉片造型局部設計不合理，如入口攻角等不合適，需進一步優化改進，才能提高動力渦輪的熱效率，進一步提升動力渦輪的輸出功率。

（2）動力渦輪內部的流場數值模擬結果表明了子午面，不同葉高截面溫度、壓力、馬赫數的分布云圖，氣動載荷分布，熵增分布的大體走勢，對流場中產生的現象有了定性的分析，對產生機理有了概括性的分析，有助于動力渦輪下一步的開發和改進。

圖10　熵增分布圖

［1］朱大鑫.渦輪增壓與渦輪增壓器［M］.北京：機械工業出版社，1992：117-130.

［2］柯玄齡，趙士杭，孟昭勇等.絕熱渦輪復合發動機動力渦輪的設計與試驗［J］.兵工學報（坦克裝甲車與發動機分冊），1991，44（4）：14-19.

［3］趙士杭，柯玄齡，孟昭勇等.絕熱渦輪復合發動機用動力渦輪的研制［J］.內燃機工程，1991，12（4）：34-39.

［4］鄔平，李雷，王錫斌.動力渦輪復合增壓柴油機性能的數值模擬研究［C］.2011年中國汽車工程學會年會論文集，2011：384-388.

［5］趙士杭，柯玄齡，孟昭勇等.絕熱渦輪復合發動機用動力渦輪的結構設計和改進［J］.熱能動力工程，1991，6（5）：245-250.

［6］王井山，葛煒，王振彪等.動力渦輪增壓系統試驗研究［J］.拖拉機與農用運輸車，2013，40（2）：29-35.

［7］王強.基于GT-POWER的渦輪-活塞式柴油機熱力分析及仿真研究［D］.中北大學，2012：48-67.

Power Turbine and Impeller Flow Field Analysis

Qi Bowu1，Wang Qiang1，Huang Qiuping2
（1.The North University of China，Taiyuan 030051，China；2.South East（Fujian）Motor Corp.，Fuzhou 350119，China）

Study the flow field of the Power turbine key parts，application GT-Power，one-dimensional model obtain a general data as boundary conditions，establishment the model of the Power turbine key parts and divide mesh，using CFX analysis obtain meridian plane，different vane is high，the aerodynamic load，the entropy distribution nephogram and analysis results，analysis the airflow field flow separation and formation mechanism of the wake flow phenomenon is helpful to improving Power turbine overall aerodynamic shape and locality aerodynamic shape，for improving and optimization Power turbine provide theory basis.

power turbine，impeller，flow field analysis

10.3969/j.issn.1671-0614.2016.01.004

來稿日期：2015-08-06

祁博武（1988-），男，碩士研究生，研究方向為增壓器性能仿真分析。