王雪東 王銀燕 楊傳雷 王賀春

摘要：渦輪的流量特性和效率特性，對渦輪與柴油機的匹配有重要影響。文章中保持渦輪邊界條件相同，改變軸流渦輪噴嘴的沖角，并采用三維數(shù)值模擬，對某型軸流渦輪噴嘴沖角特性進行了數(shù)值計算和流場分析。結(jié)果表明，發(fā)動機兩個工況下，在-40°-50°沖角下，渦輪內(nèi)氣體流動穩(wěn)定，流量收斂，在-40°-40°范圍，渦輪流量值相差不到0.4%;兩個工況下渦輪最高效率沖角約為-30°，增加或降低該沖角，渦輪效率逐漸下降，工況越低，效率變化越明顯。

Abstract： The flow rate and efficiency characteristics of the turbine have an important influence on? turbine and the diesel engine matching. This paper keeps the same turbine boundary conditions， changes the incidence of a axial turbine nozzle， the numerical calculation and flow field analysis on incidence of a turbine nozzle were carried out by the use of 3D numerical simulation. The results show that under the two working conditions of the engine， at theincidence of -40°-50°， the flow in the turbine is stable and the flow rate converges. Within the range of -40°-40°， the turbine flow rate difference is less than 0.4%; Under the two working conditions， the incidence of best efficiency is about -30°， if the incidence is increased or decreased， the turbine efficiency will gradually decrease. The lower the engine working condition is， the more obvious the efficiency change will be.

關(guān)鍵詞：噴嘴沖角;渦輪特性;流場;數(shù)值模擬

Key words： incidence of nozzle;characteristics of turbine;flow field;numerical simulation

0? 引言

發(fā)動機是機械系統(tǒng)的核心部件，進氣形式有自然吸氣、機械傳動式渦輪增壓和廢氣渦輪增壓三種[1][2]。采用機械傳動渦輪增壓會損耗發(fā)動機功率，增加油耗。為減小功率損失，采用廢氣渦輪增壓來驅(qū)動壓氣機是比較可行的做法。隨著發(fā)動機朝著高功率密度、高效率和低油耗發(fā)展，除了保證增壓系統(tǒng)與發(fā)動機實現(xiàn)良好地匹配，還應(yīng)使增壓器渦輪向高效率方向發(fā)展。

此前的研究大多集中于渦輪內(nèi)部流動情況、渦輪葉片的分析，并且主要研究對象為徑流式渦輪。軸流式渦輪相比徑流式向心渦輪的主要優(yōu)點是流通能力強、高效率的流動范圍寬。

柴油機是船舶中使用最廣的動力裝置，現(xiàn)代柴油機大多采用廢氣渦輪增壓。廢氣渦輪增壓是應(yīng)用普遍的增壓技術(shù)，可以提升柴油機功率，降低油耗。廢氣渦輪增壓器與柴油機匹配效果的好壞主要與增壓器中渦輪、壓氣機的流量特性和效率特性有關(guān)[3]。

如果渦輪與壓氣機或柴油機匹配運行時的絕熱效率偏低，則說明渦輪內(nèi)部流動情況惡化，可能會造成渦輪做功量無法滿足壓氣機需求，增壓壓力降低，進氣量減少，柴油機氣缸燃燒惡化，排溫過高等問題[4]。

若提高渦輪與柴油機匹配時的絕熱效率，可以使渦輪功率增加，提升柴油機的進氣量，使柴油機進氣與噴油量同步，改善燃燒，進而提升柴油機的性能[5]。

本文應(yīng)用三維流場數(shù)值計算模擬方法，利用流體仿真計算軟件NUMECA對渦輪內(nèi)部流動進行了分析。NUMECA是計算流體力學(xué)軟件的后起之秀，軟件核心算法使用Runge-kutta法，求解雷諾平均N-S方程。并采用殘差光順和多重網(wǎng)格加快收斂速度?；诓裼蜋C不同運行工況，對某型渦輪噴嘴環(huán)不同沖角下的流動狀態(tài)進行了對比，研究了渦輪效率與渦輪噴嘴進氣沖角的關(guān)系。

1? 數(shù)值模擬方法

1.1 渦輪葉片實體模型的建立

渦輪三維模型的建立有很多種方法。由于本文渦輪運行狀態(tài)是基于某型四缸直列低速柴油機，因此，渦輪設(shè)計參數(shù)的選取對與低速柴油機的良好匹配有至關(guān)重要的意義。本文渦輪葉片模型的設(shè)計是基于葉輪機械造型軟件NREC。NREC是一款新型的與流體機械密切相關(guān)的設(shè)計軟件。內(nèi)含Rital（徑流渦輪一維設(shè)計）、Axial（軸流渦輪一維設(shè)計）、Compal（離心式壓氣機一維設(shè)計）、Axcent（一維到三維轉(zhuǎn)化）等模塊，本文采用Axial及Axcent模塊，完成軸流渦輪葉片的三維模型設(shè)計。

渦輪參數(shù)的選取自某型船用低速機標(biāo)定工況下渦輪進口及出口的參數(shù)。具體參數(shù)見表1。

NREC的優(yōu)點在于，用戶無需進行復(fù)雜且耗時的葉片型線的設(shè)計。只需將渦輪進出口條件，及影響渦輪性能的關(guān)鍵參數(shù)，如反動度、流量系數(shù)、負(fù)荷系數(shù)等輸入，便會得到想要的渦輪幾何數(shù)據(jù)。圖1為渦輪葉片子午視圖。圖2為渦輪三維模型示意圖。表2為渦輪基本參數(shù)。

1.2 網(wǎng)格劃分

網(wǎng)格劃分采用NUMECA軟件包中的Autogrid模塊，Autogrid是一款專為葉輪機械設(shè)計的網(wǎng)格劃分軟件，生成速度快，網(wǎng)格質(zhì)量好，入口段為噴嘴前25mm處，出口段為動葉都25mm處。渦輪網(wǎng)格劃分才用HCH網(wǎng)格，將網(wǎng)格通道分成5個網(wǎng)格區(qū)域，繞葉片區(qū)域單獨采用C型結(jié)構(gòu)，來滿足流體繞葉片流動特征的捕捉。噴嘴網(wǎng)格數(shù)843087，動葉網(wǎng)格數(shù)1312505，最小正交性33.24，最大長寬比202.53，最大延展比1.8579。網(wǎng)格圖3為渦輪噴嘴網(wǎng)格圖。圖4為噴嘴尾緣處網(wǎng)格示意圖。

1.3 邊界條件

本次計算通過給定進口總壓、進口總溫、轉(zhuǎn)速、出口靜壓，來對渦輪內(nèi)部流動進行數(shù)值模擬。改變渦輪噴嘴入口沖角，來研究沖角對渦輪效率的影響。

1.4 湍流模型及收斂條件

本次計算采用的是Spalart-Allmaras湍流模型，來求解定常粘性Navier-Stokes方程，該模型是Reynolds平均法中介于零方程模型和兩方程模型間的一方程模型，能在保證計算精度的同時加快運算速度[6]。

計算中的收斂條件：全局殘差下降三個量級以上;進出口流量相對誤差不超過0.5%，且不再發(fā)生變化。

2? 計算結(jié)果及分析

2.1 渦輪工作過程焓熵圖

廢氣流經(jīng)渦輪做功時，其壓力、焓值、熵值的變化可以用焓熵圖清楚的表示，如圖5所示。

渦輪效率的定義為：實際過程氣體膨脹到出口壓力對渦輪所做的功與理想的氣體定熵膨脹到出口壓力氣體做功的比值，渦輪效率計算公式為：

2.2 渦輪流量隨沖角的變化

分別在柴油機100%負(fù)荷及50%負(fù)荷工況，保持進出口邊界條件一致。沖角范圍從-70°-60°，每隔10°計算一個渦輪流動工況點，觀察渦輪流量的變化。

經(jīng)計算發(fā)現(xiàn)，在渦輪噴嘴大沖角（-70°，-60°）下，迭代600步之后，通過渦輪的流量仍未收斂，并有不規(guī)律的波動情況，圖6為100%工況下渦輪-30°-60°、-70°沖角渦輪? ?流量收斂史追蹤圖。圖7顯示了迭代結(jié)束時流量變化。

圖6顯示了渦輪大沖角及小沖角情況流量隨迭代次數(shù)的變化情況。從圖中明顯看出，噴嘴沖角為-60°時，進口流量和出口流量從迭代開始至迭代終止都在波動，沖角為-70°時，出口流量隨著迭代次數(shù)增加趨于穩(wěn)定，而進口流量自始至終都在無規(guī)則變化。噴嘴沖角為30°時，流量隨迭代收斂于某一定值，且不在隨迭代改變。

經(jīng)計算發(fā)現(xiàn)，除了大沖角情況，其他情況下迭代結(jié)束時流量都已收斂，在柴油機100%負(fù)荷及50%負(fù)荷下流量收斂的沖角范圍均為-40°-50°，而大沖角情況，無法給出準(zhǔn)確流量值。

圖8說明，在流量收斂的沖角范圍內(nèi)，通過渦輪的流量值相差不大。在-40°-40°沖角下流量值相差不超過0.4%。其他沖角下的渦輪流量無法收斂，無法在圖8中給出。

2.3 渦輪效率隨沖角的變化

渦輪流量隨迭代次數(shù)的變化表明了渦輪內(nèi)氣體流動的穩(wěn)定性。渦輪絕熱效率則反映了氣體流經(jīng)渦輪的做功能力，效率高則渦輪能量的利用程度高，渦輪損失低。

由2.2可知，渦輪流量收斂的沖角范圍為-40°-50°，其他沖角下，流動狀態(tài)不穩(wěn)定，因此本文只對穩(wěn)定流動下渦輪效率進行分析。渦輪效率的計算方式見式（1）。柴油機不同負(fù)荷下渦輪效率隨沖角的變化由圖9和圖10給出。

圖9和圖10表明，在柴油機不同負(fù)荷下，渦輪效率最高時的沖角不是0°，而是在-30°，并且隨著沖角逐漸增加或降低，渦輪效率呈下降的趨勢，且在50%負(fù)荷時，效率下降更明顯。

2.4 仿真云圖

在發(fā)動機100%工況下，在流量收斂的沖角范圍內(nèi)，選取-30°、0°、30°、50°沖角，分析渦輪在四個噴嘴沖角下葉高中部截面壓力、熵的分布。

2.4.1 壓力分布圖

從圖11-圖14可看出，渦輪壓力從進口到出口呈總體下降的趨勢，在經(jīng)過工作輪之后壓力的降低尤為明顯。

在噴嘴沖角為30°和50°時，在靠近噴嘴前緣及吸力面處，出現(xiàn)了局部小壓力區(qū)域，且越靠近前緣和吸力面，壓力越小，說明在沖角由0°增大時，由于來流氣體與軸向成一定角度進入噴嘴環(huán)，在噴嘴環(huán)前緣及吸力面?zhèn)劝l(fā)生撞擊，造成氣體壓力急劇下降，流動損失增加，渦輪效率降低。在噴嘴沖角為-30°及0°時，未出現(xiàn)明顯的局部小壓力區(qū)域，氣體流動相比于前者較好。

渦輪動葉中的壓力分布較復(fù)雜。四個沖角下，氣體經(jīng)過噴嘴的膨脹加速，撞擊在動葉前緣及動葉壓力面中部，使局部壓力升高（如上圖動葉紅色區(qū)域），并形成靠近動葉前緣面與靠近動葉尾緣面兩個局部低壓區(qū)域（如上圖綠色區(qū)域），造成壓力損失，使渦輪效率有所降低。

綜上所述，噴嘴沖角由-30°變化到50°，渦輪內(nèi)部流動情況逐漸變差，因沖角造成的損失主要發(fā)生于噴嘴前緣及吸力面，不同沖角下動葉壓力分布基本相同。這符合了圖9中渦輪效率的變化曲線，即沖角由-30°到50°，渦輪效率逐漸降低[7]。

2.4.2 熵分布圖

熵值是物體無規(guī)則運動的體現(xiàn)，在渦輪中反映了氣體流經(jīng)渦輪時的損失，熵越大，表明氣體無規(guī)則運動越明顯，渦輪內(nèi)損失越大，良好情況下渦輪內(nèi)流動熵的分布應(yīng)是從進口到出口熵平緩增加。

從圖15-圖18可以看出，渦輪熵從進口到出口呈總體增加的趨勢，只在動葉壓力面處出現(xiàn)局部小熵區(qū)域。

噴嘴流道內(nèi)，50°沖角時，在噴嘴吸力面處，出現(xiàn)局部熵增區(qū)域，該區(qū)域從前緣吸力面開始，一直蔓延到尾緣出口，使局部損失增大，并將影響氣體流經(jīng)動葉時的流動。其他沖角下，噴嘴吸力面未出現(xiàn)明顯的熵增區(qū)域。

動葉流道內(nèi)，各個沖角下，在靠近前緣的壓力面處出現(xiàn)局部熵增區(qū)域（如動葉紅黃色區(qū)），且隨沖角從-30°變化到50°，局部熵增變化越明顯。由2.4.1中壓力分布圖可知，該區(qū)域的流體極可能出現(xiàn)渦流。渦流產(chǎn)生的原因是噴嘴來流與動葉接觸面主要在動葉壓力面中部，還有一部分在動葉前緣吸力面處，在這兩區(qū)域中間部分，由于來流氣體不直接作用在其表面，故而產(chǎn)生渦流[8]。

在動葉壓力面的熵值明顯小于吸力面?zhèn)?。?30°、0°、30°沖角的小熵值區(qū)大于50°沖角時的小熵值區(qū)。各個沖角下動葉壓力面和吸力面熵值分布的不均，從動葉上游區(qū)一直延續(xù)到了動葉下游區(qū)。

由上可知，渦輪因沖角變化引起的損失主要來源于噴嘴吸力面?zhèn)鹊臍怏w撞擊影響、動葉壓力面的渦流及動葉壓力面、吸力面流動的不均勻。沖角由-30°變化到50°，渦輪內(nèi)的損失逐漸增加，與2.3中圖9的效率曲線相符。

3? 結(jié)論

①通過對柴油機100%負(fù)荷和50%負(fù)荷下單級軸流渦輪的研究得出，在大沖角下，渦輪內(nèi)流動不穩(wěn)定，無法收斂。使流量收斂的沖角范圍為-40°-50°，且在-40°到40°沖角下，流量值相差不到0.4%。

②渦輪最高效率的噴嘴沖角在-30°附近，隨著沖角逐漸增大或減小，渦輪效率逐漸下降。發(fā)動機工況越低，效率的變化越明顯。沖角變化造成的損失主要來源于噴嘴吸力面氣體撞擊、動葉壓力面渦流及動葉壓力面、吸力面不均勻的流動。

參考文獻：

[1]Shuo Liu，Yi Cui，Yi Wang，Kangyao Deng，Sheng Liu. An evaluation method for transient response performance of turbocharged diesel engines[J]. Energy，2019，182.

[2]Sunil Kumar Pandey，Suryanarayana Vandana，S.R. Sarma Akella，R.V. Ravikrishna. Potential of Early Direct Injection （EDI） for simultaneous NO x and soot emission reduction in a heavy duty turbocharged diesel engine[J]. Applied Thermal Engineering，2019，158.

[3]楊偉.慢速船上柴油機和渦輪增壓器配合影響因素與特性研究[J].科技資訊，2019，17（13）：76-77.

[4]張海磊.發(fā)動機不同工況下渦輪增壓器軸向氣動作用力形成機制研究[D].中北大學(xué)，2019.

[5]邢世凱，馬朝臣，呂偉，陳山.葉輪沖角變化對葉輪效率的影響[J].科技導(dǎo)報，2012，30（33）：27-30.

[6]X.F. Yuan，W.R. An，Y.W. Ju，S. Antonov，Z.N. Bi，W. Li，J.T. Wu. Evaluation of microstructural degradation and its corresponding creep property in integral cast turbine rotor made of K424 alloy[J]. Materials Characterization，2019，158.

[7]蘆成英.可變截面渦輪增壓器（VGT）流場分析[D].大連交通大學(xué)，2012.

[8]Spence S W T， Artt D W. An experimental assessment of incidencelosses in a radial inflow turbine rotor [J]. Journal of Power and Energy，1998， 212（A1）： 43-53.