（東方汽輪機有限公司，四川德陽，618000）

三角形渦發(fā)生器增強U型通道換熱性能研究

梁崇治，陶志，宋立明

（東方汽輪機有限公司，四川德陽，618000）

將三角形渦發(fā)生器及導(dǎo)葉結(jié)構(gòu)應(yīng)用于U型通道中，以提高燃氣輪機高溫葉片內(nèi)部換熱性能。在驗證分析方法正確性的基礎(chǔ)上，數(shù)值研究了采用不同結(jié)構(gòu)時U型通道內(nèi)的流動及換熱特性。與典型的布置斜肋的U型通道不同，三角形渦發(fā)生器通過誘導(dǎo)產(chǎn)生下洗渦對，從而有效地強化U型通道的換熱能力。而將導(dǎo)葉布置于U型通道轉(zhuǎn)折處，可進一步提高三角形渦發(fā)生器在U型通道流出段的換熱性能，并降低通道的阻力損失。和典型的布置斜肋的U型通道相比，在不同雷諾數(shù)下帶導(dǎo)葉且布置三角形渦發(fā)生器的U型通道內(nèi)換熱的均勻性顯著提高，在通道換熱能力相當(dāng)?shù)那闆r下，阻力損失顯著降低，綜合換熱性能提高20.0%左右。

三角形渦發(fā)生器，U型通道，導(dǎo)葉，強化換熱

1 前言

為了提高燃氣輪機的熱經(jīng)濟性，透平進口溫度越來越高，已遠遠超過了葉片材料的熔點。因此，采用多種高效的冷卻技術(shù)對燃機葉片進行冷卻和熱防護是保證燃氣輪機安全可靠運行的關(guān)鍵。葉片內(nèi)部通道肋片擾流冷卻是目前常用的冷卻方式之一。最早使用在內(nèi)部通道中的肋片是正方形截面的90°直肋[1]，在氣流流過肋片時，在肋片下游產(chǎn)生流動分離和再附，增強了湍流混合，減薄了邊界層，從而提高了通道換熱性能。Han等[2-3]指出當(dāng)肋片與主流呈一定角度（如：斜肋、V型肋、W型肋）時，還能產(chǎn)生“肋片誘導(dǎo)渦”，能進一步增強主流流體和壁面低能流體的湍流混合，從而提高換熱。Han等[4]研究表明，相比于連續(xù)肋片，間斷肋片破壞了正常的二次流結(jié)構(gòu)，產(chǎn)生了額外的“誘導(dǎo)渦”，導(dǎo)致湍流混合的增強，從而提高換熱。Rallbandi等[5-6]實驗研究了波紋狀截面肋片的換熱和阻力特性，指出圓角導(dǎo)致肋片換熱效果的改變并不顯著，但是阻力損失降低明顯。

除了傳統(tǒng)的肋片擾流結(jié)構(gòu)，采用渦發(fā)生器強化通道換熱性能的研究也得到了廣泛關(guān)注。Han等[7]實驗研究了楔形和三角形渦發(fā)生器的換熱和壓力損失特性，結(jié)果表明，三角形肋片和楔形肋片相比，強化換熱效果更好，壓損更小。Liou等[8]實驗研究了在直通道內(nèi)分別布置12種渦發(fā)生器時的換熱和阻力特性，指出布置與通道等寬的三角形渦發(fā)生器具有最高的等泵功條件下的換熱性能和換熱均勻性。Henze等[9]實驗測量了在直通道內(nèi)布置單個三角形渦發(fā)生器時的換熱和流動特性，指出下洗渦對的產(chǎn)生是其強化換熱的本質(zhì)原因，并研究了渦發(fā)生器尺寸和雷諾數(shù)對換熱的影響。Henze等[10]實驗研究了一列和一排布置的三角形渦發(fā)生器對直通道內(nèi)的換熱特性的影響，指出渦系之間的相互作用是影響換熱的主要原因。近年來，在U型通道中布置導(dǎo)葉的研究也逐漸開展。Chen等[11]實驗和數(shù)值研究了導(dǎo)葉對斜肋U型通道換熱和阻力特性的影響，指出導(dǎo)葉對U型通道內(nèi)換熱的改善并不顯著，但是明顯降低了通道內(nèi)的阻力損失。Lee等[12]實驗研究了不同結(jié)構(gòu)的導(dǎo)葉對U型通道換熱和阻力特性的影響。

本文將三角形渦發(fā)生器和導(dǎo)葉應(yīng)用于U型通道中，并在驗證分析方法正確性的基礎(chǔ)上，數(shù)值研究了不同雷諾數(shù)下三角形渦發(fā)生器及導(dǎo)葉對U型通道內(nèi)的流動及換熱特性的影響，闡明了其強化換熱機理。

2 幾何模型

為了研究三角形渦發(fā)生器及導(dǎo)葉對U型通道內(nèi)流動和換熱特性的影響，構(gòu)建了3個計算模型。模型1采用文獻[13]給出的45°斜肋U型通道，一方面驗證本文的數(shù)值方法，另一方面作為參考與三角形渦發(fā)生器進行對比。另外兩個模型中U型通道尺寸均與第一個模型保持一致，U型通道的寬高比為1∶2，截面尺寸為12.7 mm×25.4 mm。模型1中，在U型通道流入段和流出段各布置9根45°斜肋，如圖1（a）所示，肋片截面為正方形，肋高e=1.59 mm，肋片間距p/e=10。模型2采用三角形渦發(fā)生器來強化U型通道換熱性能，在U型通道流入段和流出段各布置8個三角形渦發(fā)生器，三角形渦發(fā)生器距U型通道進出口距離為6.35 mm，三角形渦發(fā)生器間距為15.24 mm，如圖1（b）所示。三角形渦發(fā)生器的尺寸如圖1（c）所示。模型3在模型2的基礎(chǔ)上，在通道轉(zhuǎn)折處加裝了導(dǎo)葉。導(dǎo)葉的型線為圓弧[14]，尺寸如圖1（d）所示。

圖1 幾何模型

3 數(shù)值方法

3.1計算網(wǎng)格

本文采用ICEM軟件生成結(jié)構(gòu)化網(wǎng)格。為了保證網(wǎng)格質(zhì)量，在U型通道轉(zhuǎn)折處采用C型網(wǎng)格劃分，在三角形渦發(fā)生器區(qū)域采用Y型網(wǎng)格劃分，在導(dǎo)葉區(qū)域采用O型網(wǎng)格劃分，其他區(qū)域采用H型網(wǎng)格劃分，局部網(wǎng)格細節(jié)如圖2所示。所有壁面附近均進行網(wǎng)格加密，保證離開壁面第一層網(wǎng)格值小于1，以滿足湍流模型的要求。經(jīng)過網(wǎng)格無關(guān)性驗證，確定3個計算模型的網(wǎng)格數(shù)量分別為320萬、370萬和390萬。

圖2 計算網(wǎng)格

3.2邊界條件

本文以空氣作為工質(zhì)，采用商業(yè)軟件CFX求解三維穩(wěn)態(tài)RANS方程，湍流模型采用SST k-ω模型。本文所有計算模型的邊界條件設(shè)定與文獻[13]一致。圖3給出了帶導(dǎo)葉的三角形渦發(fā)生器U型通道的計算模型，進口給定速度Ui，靜溫298.15 K，湍流度5%，出口給定大氣壓，換熱面給定恒定壁溫338.15 K。由于在上、下壁面對稱布置擾流結(jié)構(gòu)，故計算模型只選取一半，采用對稱面邊界設(shè)定。其余壁面設(shè)定為絕熱無滑移邊界。為了保證數(shù)值計算的準(zhǔn)確性，在U型通道進、出口處均添加了長度為10倍水力直徑的延長段。其中進口速度Ui由當(dāng)量直徑Dh和雷諾數(shù)Re決定，即：

本文對雷諾數(shù)為10 000、25 000和40 000的3個工況進行了計算。

圖3 計算模型示意圖

3.3數(shù)值方法驗證

圖4給出了Re=25 000工況下的數(shù)值計算得到的斜肋通道區(qū)域平均Nu/Nu0與實驗值[13]的對比。其中，表征通道換熱性能的努賽爾數(shù)Nu的定義為：

式中：qw為壁面熱流密度；λ為流體的導(dǎo)熱系數(shù)；Tw和Tb分別為壁面溫度和冷卻氣體進、出口平均溫度；Dh為通道進口截面水力直徑。

以Nu/Nu0表示強化換熱結(jié)構(gòu)引起的換熱增強因子，其中，Nu0為相應(yīng)流動參數(shù)下光滑通道中充分發(fā)展流動的努賽爾數(shù)[13]：

從圖4可以看出，計算結(jié)果與實驗值在趨勢上吻合較好，在通道流入段計算結(jié)果最大偏大11.0%，在通道流出段計算結(jié)果最大偏大8.0%。

圖5給出了3種雷諾數(shù)下斜肋冷卻通道換熱面平均Nu/Nu0與實驗值的對比，計算值相對實驗值分別偏大8.5%、10.5%和10.8%。由于實驗[13]中換熱測量的不確定性為7%～21%，因此本文所采用的數(shù)值方法是可靠的。

圖4 Re=25 000時區(qū)域平均Nu/Nu0對比

圖5 不同雷諾數(shù)下通道平均Nu/Nu0對比

4 計算結(jié)果與分析

4.1不同結(jié)構(gòu)的換熱和阻力特性比較

圖6給出了Re=25 000時，布置不同結(jié)構(gòu)的U型通道Nu/Nu0分布云圖。相比于斜肋冷卻通道，三角形渦發(fā)生器冷卻通道中的換熱分布更加均勻。此外，布置導(dǎo)葉能顯著提高U型冷卻通道流出段的換熱性能。

圖6 不同結(jié)構(gòu)的Nu/Nu0分布云圖

圖7給出了不同雷諾數(shù)下，3種結(jié)構(gòu)的通道平均Nu/Nu0分布。可以看出，在3種雷諾數(shù)下，帶導(dǎo)葉的三角形渦發(fā)生器冷卻通道的換熱性能與典型的斜肋冷卻通道相當(dāng)，而不帶導(dǎo)葉的渦發(fā)生器通道換熱性能略有降低。

圖7 不同結(jié)構(gòu)的通道平均Nu/Nu0分布

圖8給出了不同雷諾數(shù)下，3種結(jié)構(gòu)的無量綱摩擦系數(shù)f/f0分布。其中，表征通道中的阻力特性的摩擦系數(shù)f定義為：

式中：Δp為通道進、出口壓差；u為進口速度；Dh為通道水力直徑；L為通道的流向總長度。

f0表示對應(yīng)光滑通道中充分發(fā)展流動的摩擦系數(shù)，其表達式[13]為：

圖8 不同結(jié)構(gòu)的f/f0分布

從圖8可以看出，相比于斜肋冷卻通道，采用三角形渦發(fā)生器能顯著降低U型通道的摩擦系數(shù)，而布置導(dǎo)葉能夠進一步降低U型通道的摩擦系數(shù)。

為了全面評估強化換熱結(jié)構(gòu)的換熱與阻力特性，定義綜合換熱性能[13]為：

圖9給出了不同雷諾數(shù)下，3種結(jié)構(gòu)的綜合換熱性能比較。可以看出，相比于斜肋冷卻通道，布置三角形渦發(fā)生器的兩種U型通道的綜合換熱性能均顯著提高。3種雷諾數(shù)下，不帶導(dǎo)葉的渦發(fā)生器冷卻通道相比斜肋冷卻通道分別提高9.8%、6.8%、9.1%。帶導(dǎo)葉的渦發(fā)生器冷卻通道相比斜肋冷卻通道，分別提高18.4%、17.1%、21.0%。表明三角形渦發(fā)生器和導(dǎo)葉的組合布置能顯著提高U型通道內(nèi)的綜合換熱性能。

圖9 不同結(jié)構(gòu)的綜合換熱性能對比

4.2渦發(fā)生器冷卻通道內(nèi)流動與換熱特性分析

U型通道內(nèi)復(fù)雜的三維流動結(jié)構(gòu)對其換熱有重要影響。圖10給出了布置三角形渦發(fā)生器的U型通道內(nèi)渦系結(jié)構(gòu)。渦發(fā)生器誘導(dǎo)產(chǎn)生了反向旋轉(zhuǎn)的下洗渦對，該渦對的兩個分支在向下游發(fā)展過程中逐漸向兩側(cè)分開，最終會向上抬升與上游的下洗渦對卷吸并匯合到一起，并繼續(xù)往下游運動。在渦發(fā)生器的側(cè)壁會形成角渦，角渦沿著側(cè)壁向下游運動，從渦發(fā)生器脫離以后，最終與通道側(cè)壁的角渦匯合，且兩者旋向一致。這些渦系結(jié)構(gòu)的產(chǎn)生和交互作用產(chǎn)生了強化換熱的效果。

圖10 渦核結(jié)構(gòu)示意圖

圖11給出了斜肋和三角形渦發(fā)生器通道流向截面的流線分布。在斜肋的下游會產(chǎn)生較大流動分離，這是其流動阻力大的主要原因。而三角形渦發(fā)生器產(chǎn)生的縱向渦系結(jié)構(gòu)，對主流的阻滯作用更小，在下游基本上不會產(chǎn)生流動分離區(qū)，因而三角形渦發(fā)生器具有很低的流動阻力。

圖11 流向截面流線

圖12給出了流入段4個橫截面上的二次流流線和流向渦量分布。下洗渦對形成于渦發(fā)生器的頂部，然后沿著頂面下洗運動并不斷發(fā)展壯大。經(jīng)過渦發(fā)生器以后，下洗渦對在向下游發(fā)展過程中，下洗渦對的兩個分支逐漸向兩側(cè)分開，并且渦對的強度逐漸減弱。從給出的4個截面上均可以看到，通道底面和兩側(cè)壁面的夾角處形成了壁角渦，其旋向與對應(yīng)側(cè)的下洗渦對分支的旋向相反。另外，在渦發(fā)生器左右兩側(cè)存在范圍較大且強度較弱的渦對，這是上游的渦系在下游的延續(xù)發(fā)展形成的，這會導(dǎo)致下游下洗渦對的強度不斷增加，因此使得其換熱性能不斷加強。

圖12 流入段橫截面處流線及渦量分布

圖13給出了截面3和截面4與通道換熱面的交線上的Nu/Nu0分布。交線1和交線2上Nu/Nu0最大的位置與圖12（d）和12（e）中箭頭所示位置是對應(yīng)的。在交線1上，換熱最強的位置出現(xiàn)在通道中心線附近，而隨著下洗渦對向下游兩側(cè)發(fā)展，換熱最強的位置出現(xiàn)在通道兩側(cè)下洗渦對分支的渦核附近，如圖12（e）所示。這是由于下洗渦對增強了主流流體和近壁面低能流體之間的湍流混合，顯著提高了強化換熱效果。由于下洗渦對往下游兩側(cè)的發(fā)展過程中，其強度減弱，故其增強換熱的效果有所降低，但是換熱性能的橫向分布更加均勻。另外，交線2上靠近通道兩側(cè)壁面位置出現(xiàn)Nu/Nu0升高，且高于交線1對應(yīng)位置，這是由于渦發(fā)生器側(cè)壁角渦與通道側(cè)壁角渦匯合并共同作用的結(jié)果。

圖13 不同位置橫向Nu/Nu0分布

U型通道轉(zhuǎn)折處流動復(fù)雜，導(dǎo)致流出段的流動與換熱特性不同于流入段[15]。圖14對比了Re= 25 000工況下不帶導(dǎo)葉和帶導(dǎo)葉渦發(fā)生器冷卻通道沿流動方向的區(qū)域Nu/Nu0分布?？梢钥闯觯瑢?dǎo)葉對流入段的流動與換熱特性影響甚微，與文獻[11]的結(jié)果一致。在流入段，兩種通道的Nu/Nu0均逐漸升高，這是由于下洗渦對在向下游運動過程中與上游渦系匯合并逐漸增強所導(dǎo)致的。導(dǎo)葉對流出段的換熱有顯著的影響，采用導(dǎo)葉使得U型通道轉(zhuǎn)折區(qū)域的Nu/Nu0提高約4%，流出段Nu/Nu0提高約17%。

圖14 區(qū)域平均Nu/Nu0分布

圖15給出了不帶導(dǎo)葉和帶導(dǎo)葉的渦發(fā)生器冷卻通道轉(zhuǎn)折處的壁面極限流線和壓力分布。從圖中可以看出，導(dǎo)葉抑制了轉(zhuǎn)折處流動分離的發(fā)生，從而降低了流動阻力。氣流經(jīng)過通道轉(zhuǎn)折區(qū)域時，發(fā)生沖擊、偏轉(zhuǎn)、流動分離，導(dǎo)致流線向流出段內(nèi)側(cè)偏移。而安裝了導(dǎo)葉以后，氣流方向與通道中心線平行，更加有利于下洗渦對的形成與發(fā)展。

圖15 轉(zhuǎn)折處極限流線和壓力

圖16給出了不帶導(dǎo)葉和帶導(dǎo)葉的渦發(fā)生器冷卻通道的二次流渦核結(jié)構(gòu)。采用導(dǎo)葉結(jié)構(gòu)后，削弱了氣流180°轉(zhuǎn)折導(dǎo)致的復(fù)雜二次流結(jié)構(gòu)，增強了流出段下洗渦對靠近通道外側(cè)的分支的強度，使得下洗渦對的強度和對稱性得到明顯提升，因而具有和流入段相似的分布。

圖16 U型通道二次流渦核結(jié)構(gòu)

5 結(jié)論

本文采用數(shù)值方法，對不同雷諾數(shù)下布置導(dǎo)葉及三角形渦發(fā)生器結(jié)構(gòu)的U型通道的換熱性能及阻力特性進行了研究，主要結(jié)論如下：

（1）在不同雷諾數(shù)下，不帶導(dǎo)葉和帶導(dǎo)葉的三角形渦發(fā)生器U型冷卻通道的綜合換熱性能均優(yōu)于斜肋U型通道。相比于典型的斜肋U型通道，帶導(dǎo)葉且布置三角形渦發(fā)生器的U型通道內(nèi)的換熱均勻性顯著提高，并且在換熱能力相當(dāng)?shù)那闆r下，阻力損失顯著降低，綜合換熱性能提高20%左右。

（2）三角形渦發(fā)生器通過誘導(dǎo)產(chǎn)生下洗渦對，增強了主流與近壁面低能流體之間的湍流混合，能有效強化U型通道內(nèi)的換熱性能。另外，不同于斜肋結(jié)構(gòu)，三角形渦發(fā)生器產(chǎn)生的縱向渦系結(jié)構(gòu)，對主流的阻滯作用更小，且不會導(dǎo)致下游分離區(qū)的產(chǎn)生，從而顯著降低了流動阻力。

（3）U型通道轉(zhuǎn)折處的流動非常復(fù)雜，在U型通道的轉(zhuǎn)折處布置導(dǎo)葉，可以明顯改善轉(zhuǎn)折處的流動，顯著提高布置三角形渦發(fā)生器的U型通道流出段的換熱性能，并進一步降低通道內(nèi)的流動阻力損失。

[1]Webb R L,Eckert E R G,Goldstein R J.Heat transfer and friction in tubes with repeated-rib roughness[J].Internation?al Journal of Heat and Mass Transfer,1971,14(4):601-617.

[2]Han J C,Zhang Y M,Lee C P.Augmented heat transfer in square channels with parallel,crossed,and V-shaped angled ribs[J].Journal of Heat Transfer,1991,113(3):590-596.

[3]Park J S,Han J C,Huang Y,et al.Heat transfer perfor?mance comparisons of five different rectangular channels with parallel angled ribs[J].International Journal of Heat and Mass Transfer,1992,35(11):2891-2903.

[4]Han J C,Zhang Y M.High performance heat transfer ducts with parallel broken and V-shaped broken ribs[J].Interna?tional Journal of Heat and Mass Transfer,1992,35(2):513-523.

[5]Rallabandi A P,Yang H,Han J C.Heat transfer and pres?sure drop correlations for square channels with 45 deg ribs at high reynolds numbers[J].Journal of Heat Transfer,2009, 131(7):071703.

[6]Rallabandi A P,Alkhamis N,Han J C.Heat transfer and pressure drop measurements for a square channel with 45 deg round-edged ribs at high reynolds numbers[J].Journal of Turbomachinery,2011,133(3):031019.

[7]Han J C,Huang J J,Lee C P.Augmented heat transfer in square channels with wedge-shaped and delta-shaped tur?bulence promoters[J].Journal of Enhanced Heat Transfer, 1993,1(1).

[8]Liou T M,Chen C C,Tsai T W.Heat transfer and fluid flow in a square duct with 12 different shaped vortex gen?erators[J].Journal of Heat Transfer,2000,122(2):327-335.

[9]Henze M,Dietz C,Neumann S O,et al.Heat transfer en?hancement from single vortex generators[C]//ASME Turbo Expo 2005:Power for Land,Sea,and Air.American Society of Mechanical Engineers,2005:121-129.

[10]Henze M,Dietz C,von Wolfersdorf J,et al.Heat transfer around longitudinal and parallel arranged wedge-shaped vortex generators[C]//ASME Turbo Expo 2007:Power for Land,Sea,and Air.American Society of Mechanical Engineers,2007:693-702.

[11]Chen W,Ren J,Jiang H.Effect of turning vane configura?tions on heat transfer and pressure drop in a ribbed internal cooling system[J].Journal of Turbomachinery,2011,133 (4):041012.

[12]Lee D M,Park J S,Lee D H,et al.Heat-transfer character?istics of a non-rotating two-pass rectangular duct with var?ious guide vanes in the tip turn region[J].Journal of Tur?bomachinery,2012,134(5):051039.

[13]Liu Y H,Wright L M,Fu W L,et al.Rib spacing effect on heat transfer and pressure loss in a rotating two-pass rectangular channel(AR=1:2)with 45-degree angled ribs [C]//ASME Turbo Expo 2006:Power for Land,Sea,and Air.American Society of Mechanical Engineers,2006: 363-373.

[14]Schüler M,Zehnder F,Weigand B,et al.The effect of turning vanes on pressure loss and heat transfer of a ribbed rectangular two-pass internal cooling channel[J].Journal of Turbomachinery,2011,133(2):021017.

[15]Luo J,Razinsky E H.Analysis of turbulent flow in 180 deg turning ducts with and without guide vanes[J].Journal of Turbomachinery,2009,131(2):021011.

Study on Heat Transfer Enhancement of U-shaped Channel Equipped with Delta-shaped Vortex Generator

Liang Chongzhi，Tao Zhi，Song Liming
（Dongfang Turbine Co.,Ltd.,Deyang Sichuan,618000）

In order to improve the heat transfer performance in gas turbine blade internal cooling system,the delta-shaped vortex gen?erators and turning vane are used in the U-shaped channel.Upon numerical validation,the flow and heat transfer characteristic of the U-shaped channel equipped with different structures are analyzed.Unlike the U-shaped channel equipped with typical angled ribs, the delta-shaped vortex generators can effectively enhance the heat transfer capability in the U-shaped channel by inducing downwash vortex pairs in the flow passage.Then,by adding a turning vane at the turning of the U-shaped channel equipped by vortex gener?ators,the heat transfer capability in outlet channel is further improved,and the friction loss in the passage is reduced.The uniformity of heat transfer in U-shaped channel equipped with turning vane and delta-shaped vortex generators is significantly improved at differ?ent Reynolds numbers,when compared with that of the U-shaped channel equipped with typical angled ribs.Provided that the heat transfer enhancement in flow passage of U-shaped channel is equal,the friction loss of the former design equipped with turning vane and delta-shaped vortex generators is significantly reduced.Consequently,the thermal performance of the former design is increased by 20%.

delta-shaped vortex generator，U-shaped channel，turning vane，heat transfer enhancement

V232

A

1674-9987（2016）04-0018-07

10.13808/j.cnki.issn1674-9987.2016.04.005

梁崇治（1977-），男，工學(xué)學(xué)士，工程師，畢業(yè)于西安交通大學(xué)熱動力專業(yè)，主要從事熱動力機械設(shè)計、質(zhì)量管理工作。