席雷,高建民,徐亮,趙振,李云龍
(西安交通大學機械制造系統(tǒng)工程國家重點實驗室,710049,西安)
高效的葉片冷卻技術是提高燃氣輪機循環(huán)熱效率和輸出功率的重要手段之一。肋片擾流冷卻是燃氣輪機渦輪葉片中弦區(qū)域最常用的高效冷卻技術之一[1]。Han發(fā)展了不同結構參數(shù)下空氣冷卻帶肋通道的流動和傳熱關聯(lián)式[2]。Mohammadi等研究指出肋頂和肋后區(qū)域的當?shù)貍鳠嵯禂?shù)對操作和幾何不確定性更為敏感[3]。Dritselis采用大渦模擬方法對帶有方形、圓形和三角形組合肋片的冷卻通道內(nèi)的流動特性進行了詳細的研究,并分析了通道內(nèi)細微流動結構[4]。Kaewchoothong等研究了傾斜肋、V型肋和倒V型肋對靜止方形通道傳熱系數(shù)的影響,指出60°V型肋的傳熱性能最優(yōu),其平均努塞爾數(shù)比直肋提高了30%[5]。Wang等對帶波紋肋直通道內(nèi)冷氣的流動和傳熱特性進行了數(shù)值模擬研究,指出優(yōu)化后的波紋肋的傳熱性能優(yōu)于45°V型肋[6-7]。Jang等實驗研究指出凹窩或凸起增強了帶肋通道的傳熱性能[8]。Moon等對16種不同形狀肋片的傳熱性能進行了研究,指出靴子型肋片的傳熱性能最優(yōu)[9]。劉聰?shù)妊芯苛私孛嫘螤顬槿切蔚腣肋和反向V肋的傳熱性能,指出三角形反向V肋的傳熱效果較好[10]。但是,帶肋通道的結構越來越復雜,使得葉片加工困難,制造成本昂貴。
針對空氣冷卻結構復雜的問題,可以通過對簡單帶肋通道結構進行優(yōu)化來解決。目前,許多學者將神經(jīng)網(wǎng)絡和遺傳算法應用到冷卻通道流動和傳熱性能的建模、預測和優(yōu)化等研究中,指出該方法具有很高的預測效率和尋優(yōu)能力[11]。Seo等采用遺傳算法和近似模型對冷卻通道中的靴型肋進行了多目標優(yōu)化[12]。Kim等通過數(shù)值模擬方法分析了帶V型肋冷卻通道的流動和傳熱特性,之后基于數(shù)值結果采用多項式響應面模型和NSGA-II(帶精英策略的非支配排序遺傳算法)對V型肋片的結構參數(shù)進行了優(yōu)化[13]。Darvish等使用神經(jīng)網(wǎng)絡和遺傳算法對冷卻通道中的V型肋進行多目標優(yōu)化,獲得了良好的結果[14]。Xi等采用神經(jīng)網(wǎng)絡和遺傳算法相結合對蒸汽冷卻厚壁帶肋通道傳熱性能進行了預測和優(yōu)化[15]。Haasenritter等指出遺傳算法是一種整體性能良好的可用于渦輪葉片內(nèi)冷通道肋片結構優(yōu)化的好方法[16]。阮有鋼等采用Kriging模型和多目標遺傳算法開展了渦輪葉片U型通道內(nèi)部流動與傳熱的多目標優(yōu)化研究[17]。龔建英等基于代理模型和遺傳算法優(yōu)化得到了U型通道的最佳肋片角度、肋片方向和肋片距離[18]??梢?神經(jīng)網(wǎng)絡結合遺傳算法是一種可用于對冷卻通道傳熱性能的建模、預測和尋優(yōu)的合適的、可靠的方法。
綜上,盡管一些學者采用神經(jīng)網(wǎng)絡結合遺傳算法對帶肋通道開展了優(yōu)化研究,然而針對渦輪葉片厚壁帶肋通道的冷卻性能預測和結構參數(shù)優(yōu)化的研究還鮮見。因此本文基于已有實驗數(shù)據(jù)開展了雷諾數(shù)為10 000~60 000、寬高比為0.25~4和肋角度為30°~90°的渦輪葉片厚壁帶肋通道冷卻性能研究。利用Sobol方法分析了厚壁帶肋通道流動及傳熱性能對雷諾數(shù)、寬高比、肋角度的全局敏感性;采用BP神經(jīng)網(wǎng)絡對60組實驗數(shù)據(jù)進行了多輸入、多輸出的非線性擬合,得到了可以準確預測厚壁帶肋通道摩擦系數(shù)、平均努塞爾數(shù)和綜合熱力系數(shù)的神經(jīng)網(wǎng)絡模型,并結合遺傳算法對厚壁帶肋通道在研究范圍內(nèi)的不同工況下進行結構優(yōu)化。研究結果可為未來重型燃機葉片冷卻結構設計提供參考。
本文研究對象為?;阅硿u輪葉片內(nèi)冷通道的厚壁帶肋通道。針對?;暮癖趲Ю咄ǖ?本課題組曾開展了10種典型厚壁帶肋通道流動及傳熱性能的實驗測量工作[19],相關實驗平臺、測量過程及實驗數(shù)據(jù)分析參見文獻[20]。
圖1給出了厚壁帶肋通道的結構示意圖,通道的材質(zhì)為304不銹鋼(0Cr18Ni9),長度L為1 000 mm,壁面厚度為3 mm。在通道的兩個尺寸較寬的壁面上等間距布置了方形截面的擾流肋片。實驗測量了通道寬高比和肋角度對帶肋通道流動和傳熱性能的影響規(guī)律,通道寬高比W/H變化范圍為0.25~4.00,肋角度α變化范圍為30°~90°,肋高比e/D恒定為0.047,肋間距比P/e恒定為10。
圖1 厚壁帶肋通道
實驗測量在本課題組所在的渦輪葉片雙工質(zhì)冷卻實驗室內(nèi)的單元帶肋通道冷卻測試平臺完成,實驗數(shù)據(jù)的特性分析可參見文獻[20]。
表1 實驗測量得到的和G
雷諾數(shù)Re
Re=ρuD/μ
(1)
式中:ρ為空氣密度;u為進口空氣速度;D為通道當量直徑;μ為空氣動力黏度。
當?shù)嘏麪枖?shù)Nu
Nu=qD/[(Tw-Tf)λ]
(2)
式中:q為壁面熱流密度;Tw為壁面溫度;Tf為根據(jù)進出口空氣溫度線性插值得到的參考溫度;λ為空氣導熱系數(shù)。
摩擦系數(shù)f
f=ΔpD/(2ρLu2)
(3)
式中:Δp為通道進出口的壓差。
綜合熱力系數(shù)G
(4)
(a)努塞爾數(shù)
敏感性分析是研究輸入?yún)?shù)的變化對輸出參數(shù)的影響程度的方法。通過敏感性分析可以獲得各輸入?yún)?shù)對輸出參數(shù)的重要程度,敏感性系數(shù)越大,輸入?yún)?shù)越重要。本文采用目前文獻中應用廣泛的基于方差分解的Sobol全局敏感性分析方法開展厚壁帶肋通道流動及傳熱性能的敏感性分析。Sobol方法的原理是通過計算得到每個輸入?yún)?shù)對輸出參數(shù)方差的貢獻度,從而確定各輸入?yún)?shù)對輸出參數(shù)的相對影響程度。因此該方法可以穩(wěn)健、可靠地計算出高度非線性模型中各輸入?yún)?shù)之間相互作用所產(chǎn)生的敏感性[21]。Sobol方法可以計算得到各輸入?yún)?shù)的一階、二階和總敏感性系數(shù),本文關注厚壁帶肋通道各輸入?yún)?shù)的總體影響程度,因此只進行各輸入?yún)?shù)的總敏感性分析,其計算公式為[22]
(5)
在文獻[20]中擬合得到了厚壁帶肋通道壁面平均努塞爾數(shù)和摩擦系數(shù)有關于寬高比、肋角度和進口雷諾數(shù)的經(jīng)驗關聯(lián)式,具體如下
0.243(W/H)-0.61](0.007 3α2-0.82α-4.19)
(6)
f=0.012 5Re-0.012(W/H)0.41α0.21
(7)
基于式(6)和(7)開展了厚壁帶肋通道的流動與傳熱性能對通道寬高比、肋角度和進口雷諾數(shù)的敏感性分析,結果如圖3所示。
圖3 敏感性系數(shù)分布
從圖3中可以看出,厚壁帶肋通道的壁面平均努塞爾數(shù)對進口雷諾數(shù)變化的敏感性很高,而對通道寬高比和肋角度的變化的敏感性較低且較為接近。通道摩擦系數(shù)對通道寬高比變化的敏感性最高,其次為肋角度,而對進口雷諾數(shù)變化的敏感性非常低(為0.004),這是因為雷諾數(shù)增大時通道進出口的壓差和通道進口冷氣的速度平方同時增大,導致了摩擦系數(shù)隨雷諾數(shù)的變化很小。從圖3中還可以看出,通道壁面平均努塞爾數(shù)對雷諾數(shù)的敏感性遠高于摩擦系數(shù)對雷諾數(shù)的敏感性,說明雷諾數(shù)的變化對厚壁帶肋通道傳熱性能的影響遠大于對通道流動性能的影響。通道壁面平均努塞爾數(shù)對通道寬高比的敏感性明顯小于摩擦系數(shù)對通道寬高比的敏感性,說明寬高比的變化對厚壁帶肋通道傳熱性能的影響明顯小于對通道流動性能的影響。通道壁面平均努塞爾數(shù)和摩擦系數(shù)對肋角度的敏感性較為接近,說明肋角度的變化對厚壁帶肋通道傳熱和流動性能的影響大致相同。綜上,當希望強化厚壁帶肋通道的傳熱性能時,通道寬高比和肋角度這兩個參數(shù)同樣重要;當希望降低厚壁帶肋通道的流動損失時,應該重點關注通道的寬高比。
人工神經(jīng)網(wǎng)絡是一種模仿生物神經(jīng)網(wǎng)絡結構和功能的計算模型,通過調(diào)整內(nèi)部大量節(jié)點之間相互連接的關系來處理輸入輸出的對應關系,被廣泛應用于復雜非線性模型的擬合。
在眾多神經(jīng)網(wǎng)絡中,BP神經(jīng)網(wǎng)絡是一種按誤差反向傳播(BP)算法訓練的多層前饋網(wǎng)絡[23],具有很強的非線性擬合能力,且精度高,非常適合處理非線性回歸問題。因此,BP神經(jīng)網(wǎng)絡在各研究領域中得到了廣泛的應用。BP神經(jīng)網(wǎng)絡的非線性回歸過程如圖4所示。本文使用Python完成了神經(jīng)網(wǎng)絡的非線性擬合和預測。
圖4 神經(jīng)網(wǎng)絡擬合和預測流程圖
為了提高神經(jīng)網(wǎng)絡模型的訓練速度和擬合精度,對樣本數(shù)據(jù)進行了歸一化處理,采用歸一化處理之后的數(shù)據(jù)進行神經(jīng)網(wǎng)絡模型的訓練。由于本文數(shù)據(jù)量較小,因此選擇最大-最小標準化方法來實現(xiàn)數(shù)據(jù)的歸一化處理。最大-最小標準化的原理是對原始數(shù)據(jù)進行線性變換,并映射到[0,1]區(qū)間。這種方法魯棒性較差,最大值與最小值非常容易受異常點影響,因此只適合傳統(tǒng)精確的小數(shù)據(jù)場景。其計算式如下
x′=(x-minX)/(maxX-minX)
(8)
式中:x表示原始數(shù)據(jù);x′表示歸一化之后的數(shù)據(jù);X表示數(shù)據(jù)樣本集。
在訓練神經(jīng)網(wǎng)絡模型時,將實驗數(shù)據(jù)樣本點劃分為兩部分,即訓練集和測試集,一般情況下,訓練集占比75%~90%,測試集占比10%~25%。本文選取10個通道的60組實驗數(shù)據(jù),其中54組用于網(wǎng)絡訓練,6組用于預測。BP神經(jīng)網(wǎng)絡由3層神經(jīng)元組成,即輸入層、隱藏層和輸出層。本文中有3個輸入?yún)?shù),分別為進口Re、通道寬高比W/H和肋角度α,3個非線性函數(shù)的輸出參數(shù),即平均努塞爾數(shù)、摩擦系數(shù)和綜合熱力系數(shù)。因此,本文BP神經(jīng)網(wǎng)絡的結構為3-5-3,即輸入層3個節(jié)點,隱藏層5個節(jié)點,輸出層3個節(jié)點,具體的神經(jīng)網(wǎng)絡模型如圖5所示。本文選擇的隱含層節(jié)點轉移函數(shù)為Logsig函數(shù),輸出層節(jié)點轉移函數(shù)為Purelin函數(shù);訓練函數(shù)為Trainlm,學習函數(shù)為Learngdm,最大訓練次數(shù)設為10 000,訓練目標偏差設為0.02,學習速率設為0.1,為防止神經(jīng)網(wǎng)絡模型過擬合,使用L2正則化方法,其他參數(shù)均取默認值。
圖5 神經(jīng)網(wǎng)絡模型
圖6給出了BP神經(jīng)網(wǎng)絡的預測偏差分布??梢钥闯?BP神經(jīng)網(wǎng)絡對平均努塞爾數(shù)預測的最大偏差為3.9%,對摩擦系數(shù)f預測的最大偏差為2.8%,對于綜合熱力系數(shù)G預測的最大偏差為4.6%。BP神經(jīng)網(wǎng)絡對平均努塞爾數(shù)、摩擦系數(shù)和綜合熱力系數(shù)預測偏差的平均值均小于2%,表明預測偏差分布較為合理。因此,訓練后的BP神經(jīng)網(wǎng)絡模型可以較為準確地預測厚壁帶肋通道平均努塞爾數(shù)、摩擦系數(shù)和綜合熱力系數(shù),可以用來考察和研究厚壁帶肋通道的流動、傳熱以及綜合熱力性能。本文獲得的BP神經(jīng)網(wǎng)絡模型可為燃氣輪機渦輪葉片內(nèi)部帶肋通道的結構設計提供高效、便利的工具。
圖6 BP神經(jīng)網(wǎng)絡預測偏差
本節(jié)采用遺傳算法對厚壁帶肋通道的結構參數(shù)進行優(yōu)化研究,并以文獻[20]中傳熱性能和綜合熱力性能最優(yōu)的實驗通道7為參考通道進行優(yōu)化結果的對比分析,參考通道的結構參數(shù)為:W/H=2.00,α=60°。
遺傳算法是基于自然選擇和基因遺傳機理的全局搜索優(yōu)化方法,由于其具有很強的全局尋優(yōu)能力,被廣泛地應用于各學科、領域之中,近年來,也被不斷地應用于傳熱領域的優(yōu)化問題之中[24]。遺傳算法的原理是在每一次迭代中都按照一定的適應度函數(shù)并通過遺傳中的選擇、交叉和變異對個體進行篩選,將適應度高的個體保留并組成新群體,經(jīng)過多次迭代,新群體“演變”為最佳解決方案。由于是隨機搜索過程,因此遺傳算法尋優(yōu)時不易陷入局部最優(yōu)值,且收斂速度快,具有較強的魯棒性。遺傳算法優(yōu)化流程圖如圖7所示。
圖7 遺傳算法的優(yōu)化流程
遺傳算法主要包括物種初始化、計算適應度函數(shù)、選擇、交叉和變異等操作。根據(jù)本文數(shù)據(jù)樣本的特點,其參數(shù)設置為:種群個體數(shù)量為10,交叉概率為0.4,變異概率為0.2,最大迭代次數(shù)為1 000。
由文獻[20]可知,在不同雷諾數(shù)下,通道摩擦系數(shù)隨著寬高比和肋角度的減小均大致呈現(xiàn)出降低的趨勢,即W/H=0.25、α=30°時通道的摩擦系數(shù)最低,故不需要針對通道寬高比和肋角度進行摩擦系數(shù)最小化的優(yōu)化設計。
2.3.1 基于平均努塞爾數(shù)的優(yōu)化 根據(jù)文獻[20]中對帶肋通道傳熱性能的分析可知,在不同的雷諾數(shù)下,均存在著最優(yōu)的寬高比W/H和肋角度α,使得帶肋通道的傳熱性能達到最優(yōu)。因此,需要對W/H和α進行以平均努塞爾數(shù)最大化為目標函數(shù)的參數(shù)優(yōu)化。首先在不同雷諾數(shù)時基于平均努塞爾數(shù)最大化對寬高比和肋角度進行了尋優(yōu),建立了如下的單目標優(yōu)化問題。
最大化
約束條件
0.25≤W/H≤4.00
30°≤α≤90°
Re=10 000~60 000
圖8給出了雷諾數(shù)為30 000和60 000時基于平均努塞爾數(shù)最大化的遺傳算法尋優(yōu)的收斂曲線,可以看出,迭代步數(shù)為200步時已經(jīng)大致達到了收斂,200步以后平均努塞爾數(shù)的值不再變化。
圖8 基于平均努塞爾數(shù)優(yōu)化的遺傳算法收斂曲線
圖9給出了不同雷諾數(shù)時基于平均努塞爾數(shù)最大化得到的最優(yōu)寬高比和肋角度、相應的壁面平均努塞爾數(shù)以及相同雷諾數(shù)下優(yōu)化后通道的平均努塞爾數(shù)值較參考通道的平均努塞爾數(shù)值提高的百分比ηNu。從圖9可見,雷諾數(shù)為10 000~60 000時,基于平均努塞爾數(shù)最大化的最優(yōu)寬高比和肋角度分別在2.23~2.75和41.12°~55.16°。以實驗中傳熱性能最優(yōu)的參考通道為參考,不同雷諾數(shù)時基于平均努塞爾數(shù)最大化的優(yōu)化通道的傳熱性能較參考通道提高了5.67%~50.21%。
圖9 不同雷諾數(shù)時基于平均努塞爾數(shù)最大化的優(yōu)化結果
從圖9中還可以看出,雷諾數(shù)對基于平均努塞爾數(shù)最大化得出的帶肋通道最優(yōu)寬高比和肋角度有較大的影響。隨著雷諾數(shù)的增大,最優(yōu)寬高比呈指數(shù)形式增大,而最優(yōu)肋角度呈對數(shù)形式增大,相應的平均努塞爾數(shù)幾乎呈線性增大。這一結果的指導意義在于當以強化傳熱為設計目標時,在通道進口雷諾數(shù)較大的情況下,應選擇本文得出的最優(yōu)寬高比(2.23~2.75)和肋角度(41.12°~55.16°)范圍內(nèi)相對較大的寬高比和肋角度。此外,在燃氣輪機渦輪葉片內(nèi)部冷卻通道的設計中,當冷氣的流量或流量范圍確定時,可以通過本文訓練得到的BP神經(jīng)網(wǎng)絡模型,采用遺傳算法進行優(yōu)化設計,獲得傳熱性能最優(yōu)所對應的寬高比和肋角度。
一般來講,在燃氣輪機渦輪葉片冷卻結構的設計中,由于葉型已經(jīng)確定,而前緣需要設計復雜的沖擊、氣膜或其復合冷卻方式,尾緣需要設置尾緣槽等,因此留給中弦區(qū)的空間較小。對于特定的葉片,中弦區(qū)帶肋通道的寬高比已經(jīng)大致確定。這就需要得到特定寬高比時的最優(yōu)肋角度,研究清楚寬高比對最優(yōu)肋角度的影響規(guī)律。通過本文訓練得到的BP神經(jīng)網(wǎng)絡模型,采用遺傳算法對特定Re和特定寬高比時的肋角度進行尋優(yōu),即可實現(xiàn)此目標。本小節(jié)下文將以雷諾數(shù)60 000為例,基于平均努塞爾數(shù)最大化對不同寬高比時(W/H=0.25,0.5,1,2,4)的肋角度進行優(yōu)化。遺傳算法的設置與前文相同。優(yōu)化的目標函數(shù)同樣為平均努塞爾數(shù)最大化,約束條件改寫如下。
約束條件
30°≤α≤90°
W/H=0.25~4.00
Re=60 000
為更加直觀地分析,圖10給出了Re為60 000時基于平均努塞爾數(shù)最大化的帶肋通道最優(yōu)肋角度隨寬高比的變化規(guī)律。從圖10中可以看出,厚壁帶肋通道最優(yōu)的肋角度隨著寬高比的增大而增大。這意味著在以強化傳熱為目標設計燃氣輪機渦輪葉片內(nèi)部冷卻通道時,如果選擇的通道寬高比較小(在0.25~4范圍內(nèi)),則應選擇最優(yōu)肋角度范圍內(nèi)(48.82°~55.13°)相對較小的肋角度,而寬高比較大時,應選擇相對較大的肋角度。
圖10 不同寬高比時基于平均努塞爾數(shù)最大化的最優(yōu)肋角度
2.3.2 基于綜合熱力系數(shù)的優(yōu)化 根據(jù)文獻[20]中對帶肋通道綜合熱力性能的分析可知,在不同的雷諾數(shù)下,也存在著最優(yōu)的寬高比W/H和肋角度α,使得帶肋通道獲得最佳的綜合熱力性能。因此,還需要對W/H和α進行以綜合熱力系數(shù)最大化為目標函數(shù)的參數(shù)優(yōu)化。
首先在不同雷諾數(shù)時基于綜合熱力系數(shù)最大化對寬高比和肋角度進行尋優(yōu),約束條件與基于平均努塞爾數(shù)最大化的優(yōu)化相同。表2給出了不同雷諾數(shù)時基于綜合熱力系數(shù)最大化的最優(yōu)寬高比和肋角度、相應的最優(yōu)綜合熱力系數(shù),以及相同雷諾數(shù)下優(yōu)化后通道的綜合熱力系數(shù)值較參考通道的綜合熱力系數(shù)值提高的百分比ηG。從表2中可以發(fā)現(xiàn),雷諾數(shù)為10 000~60 000時,基于綜合熱力性能最大化得到的最優(yōu)寬高比和肋角度分別在1.75~2.11和56.68°~65.2°,且最優(yōu)寬高比和最優(yōu)肋角度隨雷諾數(shù)的變化趨勢與基于平均努塞爾數(shù)最大化的優(yōu)化結果相同。不同雷諾數(shù)時優(yōu)化后的通道綜合熱力性能較參考通道提高了12.46%~67.03%。
表2 不同雷諾數(shù)時基于綜合熱力系數(shù)最大化的優(yōu)化結果
同樣地,本小節(jié)也以雷諾數(shù)為60 000為例,基于綜合熱力系數(shù)最大化對不同寬高比時的肋角度進行優(yōu)化,約束條件與基于平均努塞爾數(shù)最大化的優(yōu)化相同。表3給出了Re為60 000時的優(yōu)化結果,與基于平均努塞爾數(shù)的優(yōu)化結果相類似,基于綜合熱力系數(shù)優(yōu)化得到的厚壁帶肋通道的最優(yōu)肋角度也隨著寬高比的增大而增大,最優(yōu)肋角度的范圍為59.14°~66.06°。
表3 不同寬高比時基于G最大化的最優(yōu)肋角度
本文基于Sobol方法、BP神經(jīng)網(wǎng)絡和遺傳算法進行了厚壁帶肋通道流動及傳熱性能的敏感性分析、預測模型構建以及結構參數(shù)優(yōu)化研究,得到如下主要結論。
(1)厚壁帶肋通道的傳熱性能對進口雷諾數(shù)變化的敏感性很高,而對通道寬高比和肋角度變化的敏感性較低且較為接近;流動性能對通道寬高比變化的敏感性非常高,對肋角度變化的敏感性較低,而對進口雷諾數(shù)變化不敏感。
(2)訓練得到的BP神經(jīng)網(wǎng)絡模型對平均努塞爾數(shù)預測的最大偏差為3.9%,對摩擦系數(shù)預測的最大偏差為2.8%,對于綜合熱力系數(shù)預測的最大偏差為4.6%,可以用來考察和研究厚壁帶肋通道的流動、傳熱以及綜合熱力性能。
(3)最優(yōu)寬高比和最優(yōu)肋角度都隨著雷諾數(shù)的增大而增大。不同雷諾數(shù)時,基于平均努塞爾數(shù)優(yōu)化的最優(yōu)寬高比和肋角度分別在2.23~2.75和41.12°~55.16°,優(yōu)化通道的傳熱性能較參考通道提高了5.67%~50.21%。
(4)基于綜合熱力系數(shù)優(yōu)化得到的最優(yōu)寬高比和肋角度分別在1.75~2.11和56.68°~65.2°,優(yōu)化后的通道綜合熱力性能較參考通道提高了12.46%~67.03%。