楊璐璐，盧朝陽

(南京航空航天大學(xué)民航學(xué)院，南京 211106)

出流孔傾角對層板冷卻結(jié)構(gòu)流動與換熱特性的影響

楊璐璐，盧朝陽

(南京航空航天大學(xué)民航學(xué)院，南京 211106)

設(shè)計了5種不同層板冷卻結(jié)構(gòu)出流孔傾角模型，應(yīng)用Fluent軟件對其進行流動與換熱的耦合計算，并對數(shù)值計算結(jié)果進行了對比分析。結(jié)果表明:當(dāng)冷流沿出流孔流入主流時，冷流與熱流摻混較強，湍動能沿程逐漸增大;隨后摻混逐漸趨于平穩(wěn)，湍動能逐漸減小;隨著出流孔傾角的逐漸增大;出流孔中心線兩側(cè)的流體不斷被卷吸進來，不斷擠壓出流孔中心線上的流體，冷流與熱流的混合增強，摻混核心區(qū)的面積逐漸增大。隨著出流孔傾角的逐漸增大，沿出流孔流入主流的冷流覆蓋在熱側(cè)壁面上的溫度逐漸升高，換熱逐漸減弱，熱側(cè)壁面溫度逐漸升高。

傾角;層板結(jié)構(gòu);湍動能;流動換熱

未來航空發(fā)動機朝著高推重比的方向發(fā)展，為增加發(fā)動機的推力和功率，需要不斷增加直接參與燃燒的空氣量，從而導(dǎo)致用于冷卻的空氣量逐漸減少，同時燃燒室出口溫度逐漸升高［1］，給高溫部件的冷卻帶來一定困難，這就要求人們探索更為高效的冷卻結(jié)構(gòu)。

在燃燒室中采用層板冷卻結(jié)構(gòu)是解決火焰筒冷卻問題的先進技術(shù)之一。層板冷卻結(jié)構(gòu)能夠在減少冷卻用氣量的情況下降低火焰筒壁面溫度梯度和結(jié)構(gòu)熱應(yīng)力水平，延長其使用壽命。因此，具有層板冷卻結(jié)構(gòu)的火焰筒是未來燃燒室的主要結(jié)構(gòu)形式，具有良好的發(fā)展前景［2-3］。

典型的層板冷卻結(jié)構(gòu)(如 Lamilloy層板［4-5］和Transply層板［6］)，是集沖擊冷卻、強化對流換熱和氣膜冷卻于一體的高效冷卻方式。國內(nèi)外眾多學(xué)者對層板冷卻開展了研究工作［7-11］。許全宏［12］研究表明:發(fā)散壁熱側(cè)換熱增強系數(shù)隨著吹風(fēng)比的增大而增大，長菱形排布優(yōu)于正菱形排布，換熱增強系數(shù)隨著孔排數(shù)的增加而增加。宋雙文［13］研究表明:層板冷卻效率隨吹風(fēng)比的增大而增加，當(dāng)吹風(fēng)比大于1時，吹風(fēng)比對層板冷卻效率沒有影響。袁星［14］研究了填充比和直徑比對層板結(jié)構(gòu)流阻和換熱特性的影響。曹?。?5］研究表明:擾流柱的排列方式對層板冷卻效果沒有影響。譚曉茗［16］研究了7種不同層板結(jié)構(gòu)的流動和換熱特性。研究表明:沖擊孔、氣膜孔和擾流柱的排布方式以及擾流柱直徑對層板結(jié)構(gòu)的流動與換熱有顯著的影響。

國內(nèi)對層板冷卻結(jié)構(gòu)的流動與換熱特性影響因素的研究大多限于主擾流柱直徑、通道高度、不同沖擊孔、氣膜孔、吹風(fēng)比以及排布方式等參數(shù)的變化，對出流孔傾角的研究不多［17-24］。而出流孔傾角變化是層板冷卻結(jié)構(gòu)的流動與換熱特性的重要影響因素，急需開展對應(yīng)的研究。為此，本文設(shè)計了5種不同出流孔傾角的層板冷卻結(jié)構(gòu)，以Fluent商用軟件為工具，研究了出流孔傾角變化對層板冷卻結(jié)構(gòu)的流動與換熱特性的影響。

1 計算模型

1.1 物理模型

如圖1所示，本文選取叉排的2列擾流柱的局部區(qū)域進行研究。由于燃燒室半徑較大，將以上角度很小的區(qū)域簡化為矩形通道，從而得到本文研究的層板冷卻結(jié)構(gòu)模型。在如圖2所示的層板冷卻結(jié)構(gòu)示意圖中:沖擊板長L=35 mm;沖擊板厚δ=1 mm;外壁厚δ0=1.5 mm;沖擊孔距離沖擊板末端L0=35 mm;沖擊孔直徑D0=1.5 mm;第1排擾流柱距沖擊板前緣的距離L1=14.75 mm;擾流柱直徑d=0.8 mm;擾流柱軸向間距S=1.75 mm;擾流柱展向間距P=1.0 mm;擾流柱高度和通道高度均為H=2.5 mm;出流孔中心距離沖擊板前緣L2=8 mm;出流孔直徑D=1.5 mm。由于出流孔傾角的變化對熱側(cè)壁面溫度分布有很大影響，因此在保證通過沖擊孔流量相等的條件下，本文設(shè)計了層板冷卻結(jié)構(gòu)中出流孔的5種傾角模型，孔傾角α分別為20°，30°，45°，60°和90°。

1.2 計算模型與計算網(wǎng)格

圖3給出了本文的計算域模型。根據(jù)航空發(fā)動機燃燒室的實際通道高度，確定本文的熱流通道高度為42 mm，熱流進口、冷流進口和以及混合流出口如圖3所示。

圖4給出了計算域周期面局部網(wǎng)格分布。計算模型全部采用結(jié)構(gòu)化網(wǎng)格，在流體近壁處均實施了網(wǎng)格局部加密，保證了計算結(jié)果的可靠性。圖5給出了冷側(cè)壁面網(wǎng)格分布。本文分別對網(wǎng)格數(shù)為85萬、130萬、160萬、210萬和240萬的模型進行網(wǎng)格獨立性驗證。當(dāng)網(wǎng)格數(shù)為160萬時，算例滿足網(wǎng)格獨立性要求。因此，本文所有的計算模型網(wǎng)格數(shù)量均為160萬左右。

圖1 擾流柱的排布方式

圖2 層板冷卻結(jié)構(gòu)示意圖

圖3 計算域模型示意圖

圖4 周期面局部網(wǎng)格分布

圖5 冷側(cè)壁面網(wǎng)格分布

2 邊界條件與湍流模型

表1給出了冷熱流體的進出口參數(shù)，其中冷流進口采用流量進口，流量m˙=0.125 g/s。冷流體沖擊到?jīng)_擊板后，沿放置有擾流柱的冷流通道流動，最后從出流孔流出與主流混合。冷、熱流體均設(shè)置為理想流體;計算域的端面均設(shè)置為絕熱面;前后面均設(shè)置為周期面;流體與固體壁面的交界面設(shè)置為流固耦合面。計算采用Fluent商業(yè)軟件;湍流模型采用Realizablek-ε模型;控制方程采用控制容積法進行離散，采用隱式分離求解;壓力速度耦合采用SIMPLE算法;對流項采用二階迎風(fēng)進行離散，解收斂的判斷標準是相對殘差小于1×10-6。

表1 冷熱流體的進出口參數(shù)

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

3.1 速度分布

冷流沖擊到?jīng)_擊板面后，沿有擾流柱的冷卻通道流動，并從出流孔流出與熱流發(fā)生摻混。圖6 (a)～(e)分別給出了模型1～5中x/D=7截面局部區(qū)域的速度分布云圖。以圖6(a)模型1為例，x/D=7截面上冷流出流中心線y=0上速度uZ為正值，且沿Z軸正向逐漸減小。這是因為冷流沿出流孔流入主流沿程動量逐漸減小，使得冷流穿透主流的速度逐漸減小;在中心線y=0處，由于沿出流孔流出的冷流的卷吸作用，使得中心線y= 0左右兩側(cè)的流體微團發(fā)生轉(zhuǎn)動，速度uZ的方向發(fā)生改變。不同模型的x/D=7截面局部速度分布規(guī)律相似。

比較不同模型可以看出:隨著出流孔傾角的增大，中心線y=0上逐漸出現(xiàn)2個uZ為正值的峰值區(qū)域，底部區(qū)域的速度uZ峰值逐漸增大，上方區(qū)域的速度uZ峰值逐漸減小。這是因為沿出流孔流出的冷流卷吸中心線y=0兩側(cè)的流體，使中心線y=0兩側(cè)的流體微團發(fā)生轉(zhuǎn)動，隨著出流孔傾角的逐漸增大，這種轉(zhuǎn)動逐漸加劇，并逐漸擠壓中心線y=0附近區(qū)域，將中心線y=0附近的區(qū)域分為上、下2個速度uZ峰值中心。

3.2 速度矢量與流線分布

圖7(a)～(e)分別給出了模型1～5中x/D= 7截面局部區(qū)域的速度矢量與流線分布。以圖7(a)模型1為例，冷流沿出流孔流入主流，中心線y=0兩側(cè)的流體被卷吸進來，在中心線y=0兩側(cè)形成旋轉(zhuǎn)方向相反的左右對稱分布的一對渦;冷流覆蓋在熱側(cè)壁面上，隨著冷流與熱流的摻混以及換熱，其溫度沿Z軸正向逐漸升高。

圖6 x/D=7截面uZ分布云圖(單位:m/s)

比較不同模型可以看出:同一Z軸位置處，隨著出流孔傾角的逐漸增大，覆蓋在熱側(cè)壁面上的氣膜溫度逐漸降低。這是因為出流孔傾角的增大，使得冷流穿透熱流的距離增大，冷流射入主流的距離增大;隨著出流孔傾角的增大，中心線y=0兩側(cè)的渦對中心間的距離逐漸減小，這是由于中心線y=0兩側(cè)的流體被卷吸進來擠壓中心線y= 0附近區(qū)域;隨著出流孔傾角的增大，渦對中心逐漸向Z軸正向偏移，這是由于中心線y=0兩側(cè)的流體被卷吸進來，使渦對中心被抬升;隨著出流孔傾角的增大，對渦影響的區(qū)域面積逐漸增大，發(fā)生轉(zhuǎn)動的流體微團的區(qū)域逐漸增大。

圖7 x/D=7截面速度矢量與流線分布

3.3 湍動能分布

湍動能的大小和分布在一定程度上能夠反映流體間的混合速率。圖8(a)～(e)分別給出了模型1～5中x/D=7、x/D=8、x/D=9、x/D=10和x/D=11截面上局部區(qū)域的湍動能分布。

以圖8(a)模型1為例，在x/D=7截面上，冷流沿出流孔流入主流后，靠近熱側(cè)壁面附近的區(qū)域湍動能較小，此時冷流與熱流的混合較弱;隨著與熱側(cè)壁面距離的增大，湍動能逐漸增大，冷流與熱流的摻混較強;此后，隨著與熱側(cè)壁面距離的增大，湍動能逐漸減小直至0，冷流與熱流的摻混逐漸減弱。這是因為冷流沿出流孔流入主流，在壁面附近脫離壁面流動，冷流與熱流摻混較弱。在冷流流入熱流并與熱流摻混的區(qū)域內(nèi)，摻混較強，此時為摻混核心區(qū)。由于沿出流孔流入主流的冷流沿程動量逐漸減小，使得冷流穿透主流的速度和穿透距離逐漸減小，摻混逐漸減弱，在冷流與熱流接觸的區(qū)域之外，此時冷流與熱流不再發(fā)生摻混。比較不同截面可以看出:隨著x/D的增大，壁面附近流體的湍動能逐漸增大，這是由于沿出流孔流入主流的冷流，在出流孔出口處脫離壁面流動，在與主流逐漸混合過程中，流動逐漸貼緊壁面。隨著x/D的增大，摻混核心區(qū)湍動能先增大后逐漸減小，這是由于沿出流孔流入主流的冷流不斷卷吸兩側(cè)流體并與之摻混，隨后由于流動逐漸趨于平穩(wěn)，摻混逐漸減弱。

比較不同模型可以看出:隨著出流孔傾角的增大，靠近熱側(cè)壁面附近區(qū)域的湍動能逐漸增大，摻混核心區(qū)的湍動能逐漸增大，摻混核心區(qū)域的面積逐漸增大。這是因為冷流沿出流孔流入主流，不斷卷吸兩側(cè)流體進入摻混核心區(qū)，加強了冷流與熱流的摻混，同時使冷熱流體摻混核心區(qū)域的面積不斷增大。

3.4 努塞爾數(shù)分布

努塞爾數(shù)是衡量對流換熱強烈程度的重要參數(shù)，其表達式為Nu=hl/λ，其中:h表示對流換熱系數(shù)(W/m2·K);λ表示空氣的導(dǎo)熱系數(shù)(W/m·K)。

圖9給出了模型1～5中熱側(cè)壁面上努塞爾數(shù)沿流向分布。以模型1(20°)為例，熱側(cè)壁面出流孔前緣，由于出流孔的存在，流過出流孔的冷流與壁面進行了對流換熱，強化了對流換熱，努塞爾數(shù)逐漸增大;冷流沿出流孔流入主流，并形成一層保護膜完全覆蓋在熱側(cè)壁面上，此時熱側(cè)壁面與流體對流換熱量較小，努塞爾數(shù)逐漸減小;隨后由于冷流體與熱流體不斷進行摻混，擾流柱和沖擊對流的存在強化了對流換熱，努塞爾數(shù)逐漸增大。

比較不同模型可以看出:隨著出流孔傾角的增大，努塞爾數(shù)逐漸減小。這是由于隨著出流孔傾角的增大，在出流孔出口處，沿出流孔流入主流的冷流與熱側(cè)壁面發(fā)生分離，換熱量逐漸減小。

圖8 截面湍動能分布(單位:k2/s2)

圖9 熱側(cè)壁面努塞爾數(shù)沿軸向分布

3.5 熱側(cè)壁面溫度分布

圖10(a)～(e)分別給出了模型1～5中熱側(cè)壁面上溫度分布云圖。

圖10 熱側(cè)壁面溫度分布

比較不同模型可以看出:隨著出流孔傾角的增大，出流孔附近溫度等值線逐漸稀疏，溫度梯度逐漸減小，溫度變化劇烈程度減弱;隨著出流孔傾角的增大，熱側(cè)壁面上溫度等值線逐漸向下游偏移，熱側(cè)壁面上高溫區(qū)域的面積逐漸增大，即隨著出流孔傾角的增大，沿軸向同一位置處的熱側(cè)壁面溫度逐漸升高。

綜上所述，冷流沿出流孔流入主流時，冷流與熱流摻混較強，湍動能沿程逐漸增大;隨后摻混逐漸趨于平穩(wěn)，湍動能逐漸減小。隨著出流孔傾角的逐漸增大，出流孔中心線兩側(cè)的流體不斷被卷吸進來，不斷擠壓出流孔中心線上的流體，冷流與熱流的混合增強，摻混核心區(qū)的面積逐漸增大。隨著出流孔傾角的逐漸增大，沿出流孔流入主流的冷流覆蓋在熱側(cè)壁面上的溫度逐漸升高，換熱逐漸減弱，熱側(cè)壁面溫度逐漸升高。

4 結(jié)論

本文針對層板冷卻結(jié)構(gòu)出流孔傾角的改變，設(shè)計了5種不同出流孔傾角的模型，在冷流流量一定的情況下，數(shù)值模擬研究了其流動和換熱特性。在本文研究條件下，有如下結(jié)論:

1)冷流沿出流孔流入主流時，冷流與熱流摻混較強，湍動能沿程逐漸增大;隨后摻混逐漸趨于平穩(wěn)，湍動能逐漸減小。

2)隨著出流孔傾角的逐漸增大，出流孔中心線兩側(cè)的流體不斷被卷吸進來，不斷擠壓出流孔中心線上的流體，冷流與熱流的混合增強，摻混核心區(qū)的面積逐漸增大。

3)隨著出流孔傾角的逐漸增大，沿出流孔流入主流的冷流覆蓋在熱側(cè)壁面上的溫度逐漸升高，換熱逐漸減弱，熱側(cè)壁面溫度逐漸升高。

［1］Han J C，Dutta S，Ekkad S.Gas turbine heat transfer and cooling technology［M］.Florida:CRC Press，2012.

［2］侯曉春，季鶴鳴.高性能航空燃氣輪機燃燒技術(shù)［M］.北京:國防工業(yè)出版社，2002.

［3］劉軍.航空發(fā)動機的先進燃燒室技術(shù)［J］.國際航空，1997(12):48-50.

［4］Nealy D B，Reider S.Evaluation of Laminated Porous Wall Materials for Combustor Liner Cooling［J］.Journal of Engineering for Power-Transactions of the ASME，1980，102(2):268-276.

［5］Eriksen V L.Film Cooling Effectiveness and Heat Transfer with Injection through Holes［Z］.USA:［s.n.］，1971:179.

［6］Rolls-Royce AeroSpace Group.噴氣發(fā)動機［M］.北京:羅爾斯-羅伊斯公司，1986.

［7］Funazaki K，Hachiya K.Systematic numerical studies on heat transfer and aerodynamic characteristics of impingement cooling devices combined with pins［C］//ASME Turbo Expo.USA:［s.n..］，2003:16-19.

［8］Favaretto C F F，F(xiàn)unazaki K.Application of genetic algorithms to design of an internal turbine cooling system［C］//ASME Turbo Expo 2003.USA:［s.n.］，2003:273 -278.

［9］郁新華，全棟梁，劉松齡，等.層板結(jié)構(gòu)冷卻有效性的研究［J］.航空動力學(xué)報，2003，18(5):615-622.

［10］陶智，魏豪杰，丁水汀，等.典型層板冷卻結(jié)構(gòu)中流體流阻與換熱特性的實驗［J］.航空動力學(xué)報，2007，22 (2):193-198.

［11］全棟梁，劉松齡，李江海，等.多孔層板冷卻有效性的研究［J］.航空動力學(xué)報，2004(4):520-524.

［12］許全宏，林宇震，劉高恩.沖擊/發(fā)散復(fù)合冷卻方式發(fā)散壁換熱系數(shù)研究［J］.航空動力學(xué)報，2004，19(2): 213-218.

［13］宋雙文，楊衛(wèi)華，胡好生等.沖擊+逆向?qū)α?氣膜冷卻傳熱特性的研究［J］.航空動力學(xué)報，2007，22(9): 1417-1422.

［14］袁星，呂東.擾流柱幾何參數(shù)對層板流阻和換熱特性的影響［J］.航空動力學(xué)報，2009，24(5):1016-1021.

［15］曹俊，熊純，楊衛(wèi)華，等.回流燃燒室彎曲段沖擊擾流柱+逆向?qū)α?氣膜冷卻效率［J］.南京航空航天大學(xué)學(xué)報，2011，43(2):178-183.

［16］譚曉茗，李業(yè)芳，張靖周.致密多孔層板冷卻結(jié)構(gòu)研究［J］.機械工程學(xué)報，2012，48(12):144-149.

［17］全棟梁，劉松齡，李江海，等.層板冷卻特性的實驗與數(shù)值模擬研究［J］.推進技術(shù)，2004(2):134-138.

［18］全棟梁，劉松齡，李江海.層板冷卻結(jié)構(gòu)強化換熱機理［J］.航空動力學(xué)報，2004(6):860-865.

［19］王鳴，盧元麗，吉洪湖.沖擊孔對層板冷卻葉片前緣傳熱影響的數(shù)值研究［J］.航空動力學(xué)報，2013(10): 2240-2247.

［20］全棟梁，劉松齡，李江海.層板冷卻結(jié)構(gòu)傳熱優(yōu)化數(shù)值模擬研究［J］.西北工業(yè)大學(xué)學(xué)報，2004(6):816 -820.

［21］聶建豪，朱惠人，任戰(zhàn)鵬.對層板冷卻效率隨主流流動方向變化的數(shù)值研究［J］.機械與電子，2013(9):71 -73.

［22］郁新華，劉松齡，董志銳.航空發(fā)動機層板冷卻結(jié)構(gòu)的內(nèi)流阻特性的試驗研究［J］.燃氣渦輪試驗與研究，2000(3):22-26.

［23］吳向宇，黎旭，時艷，等.典型層板冷卻結(jié)構(gòu)熱疲勞破壞特性研究［J］.航空動力學(xué)報，2014(5):1177 -1183.

［24］胡婭萍，盧元麗，吉洪湖，等.層板冷卻渦輪葉片前緣內(nèi)部流動與傳熱特性實驗［J］.航空動力學(xué)報，2014 (2):241-249.

(責(zé)任編輯陳 艷)

Influence of Flow Angle of Aole on Heat Transfer and Flow Characteristics of Lamilloy Configuration

YANG Lu-lu，LU Chao-yang
(College of Civil Aviation，Nanjing University of Aeronautics and Astronautics，Nanjing 211106，China)

Five different lamilloy configurations model in which the flow angle of hole were changed were designed to conjugate simulate the flow and heat transfer characteristics with fluent software，and the results were analyzed by comparison.The results showed that when the cool air flow into the mainstream along the hole of flow，cool air and hot air mix stronger and turbulence kinetic energy increases along the way.Then the mixing gradually stabilizes，the turbulent kinetic energy decreases gradually. With the gradually increase of the flow angle of hole，the fluid on the both sides of the center line of the hole is sucked in continuous，and the fluid on the center line of the hole is squeezed out in contin-uous，which enhances cool air and hot air mixing and increases the proportion of mixing core area. With the gradually increase of the flow angle of hole，the temperature of cool air which is covered on the hot side surface increases gradually，and heat transfer are weakened and the temperature of hot side surface are increased gradually.

angle;lamilloy;turbulence kinetic energy;flow and heat transfer

V231.1

A

1674-8425(2015)11-0066-07

10.3969/j.issn.1674-8425(z).2015.11.011

2015-07-10

楊璐璐(1988—)，女，河南漯河人，碩士研究生，主要從事空中交通運輸規(guī)劃與管理研究。

楊璐璐，盧朝陽.出流孔傾角對層板冷卻結(jié)構(gòu)流動與換熱特性的影響［J］.重慶理工大學(xué)學(xué)報:自然科學(xué)版，2015(11):66-72.

format:YANG Lu-lu，LU Chao-yang.Influence of Flow Angle of Aole on Heat Transfer and Flow Characteristics of Lamilloy Configuration［J］.Journal of Chongqing University of Technology:Natural Science，2015(11):66 -72.