■ 李維 潘尚能 孟博 肖為 / 中國航發(fā)動研所

中國航發(fā)動研所作為國內(nèi)航空渦軸發(fā)動機產(chǎn)品研發(fā)基地，不斷根據(jù)自身特點，深化仿真手段在渦軸發(fā)動機研發(fā)階段的應(yīng)用，在現(xiàn)有技術(shù)條件下，最大限度地發(fā)揮氣動仿真的作用，提升了渦軸發(fā)動機性能與壽命，降低了研發(fā)周期與費用。

近年來，隨著計算機技術(shù)的不斷發(fā)展，針對航空發(fā)動機性能仿真的工具也不斷成熟，在縮短發(fā)動機研制周期、降低研發(fā)成本方面的成效不斷顯現(xiàn)[1]。中國航發(fā)動研所作為國內(nèi)主要的渦軸發(fā)動機研制單位，先后承擔(dān)過多種型號的渦軸發(fā)動機研制工作，在渦軸發(fā)動機性能仿真方面有著扎實的理論研究基礎(chǔ)和豐富的工程實踐經(jīng)驗。目前，動研所在型號研制中不斷強化仿真的工具作用，已建成高性能計算平臺，具備億級網(wǎng)格的氣動三維數(shù)值仿真能力，針對渦軸發(fā)動機氣動仿真特點開展了葉輪機流場/性能仿真、聯(lián)合/耦合氣動及換熱仿真、起動過程全尺寸燃燒室三維點火仿真、多相流仿真等研究工作，如圖1所示，對產(chǎn)品研發(fā)的全面加速起到了關(guān)鍵作用。

圖1 渦軸發(fā)動機性能仿真主要內(nèi)容

渦軸發(fā)動機氣動仿真的特點

與大型渦扇、渦噴發(fā)動機相比，渦軸發(fā)動機的主要特點是：裝機對象主要為直升機，通過旋翼產(chǎn)生強大氣流的反作用力來拉起/拉動起降/飛行，機體周圍旋翼下洗流動顯著，給發(fā)動機進排氣帶來影響；軸功率一般不超過5000kW，多數(shù)都在2000kW及以下，因此，幾何尺寸及空氣流量相對較小，結(jié)構(gòu)緊湊、轉(zhuǎn)速高，“小尺寸效應(yīng)”及“小流量效應(yīng)”顯著；構(gòu)型多樣，包括多級軸流、單級離心、“軸流+離心”和雙級離心壓氣機，直流、斜流、回流、折流燃燒室，軸流、向心渦輪，軸向/徑向進氣與軸向/偏斜式排氣，有的還帶有整體式粒子分離器；另外，直升機起降靈活，需要很強的適應(yīng)復(fù)雜惡劣環(huán)境的能力，包括山區(qū)、沙漠、冰雪地面以及海面，尤其在起降過程吸入沙石、樹葉、雜草等外物，帶來葉片磨蝕，氣路、油路堵塞，軸承磨損等不利影響。

同時，新材料、新工藝的發(fā)展，以及適航性、安全性要求等均對渦軸發(fā)動機的性能仿真提出了新挑戰(zhàn)，例如，內(nèi)部流動的尺度效應(yīng)明顯，邊界層三維效應(yīng)強，小尺寸流動湍流邊界層厚度甚至接近流動尺度，黏性力影響大，壁面摩擦和熱交換現(xiàn)象均更劇烈；氣流折轉(zhuǎn)多，內(nèi)部流場的畸變和損失控制難度大；渦軸發(fā)動機最常見的裝機位置是在直升機旋翼下方，旋翼的下洗氣流會誘發(fā)發(fā)動機喘振和左右發(fā)功率不平衡等問題[2]。

此外，渦軸發(fā)動機仿真還面臨著逆壓梯度大、流動三維效應(yīng)和轉(zhuǎn)靜子干涉明顯、流動耦合互相關(guān)效應(yīng)強、內(nèi)流氣體物理特性變化大、存在氣/液/固三項流耦合情況等發(fā)動機仿真的共性問題[3]，這些問題均給渦軸發(fā)動機的流動模擬精度的提升帶來了挑戰(zhàn)，獲得高精度性能仿真結(jié)果的難度極大。因此，如何準(zhǔn)確模擬渦軸發(fā)動機典型部件和附屬系統(tǒng)內(nèi)部的流動現(xiàn)象，以提高性能仿真的精度和可靠性，一直以來都是渦軸發(fā)動機仿真研究的重中之重。

葉輪機流場/性能仿真

與一般航空發(fā)動機類似，渦軸發(fā)動機葉輪機氣動仿真主要包括壓氣機和渦輪的氣動仿真，是渦軸發(fā)動機仿真工作的重點。動研所通過相關(guān)研究工作，使得仿真的精度和可信度已滿足工程研制需要，為加速渦軸發(fā)動機研制提供了強有力的支撐。

采用S2流面計算分析與規(guī)范化的雷諾平均方程（RAN-S）方法，對壓氣機、渦輪開展常規(guī)的設(shè)計分析與迭代優(yōu)化工作，可以在較短時間內(nèi)以較低成本提供葉輪力學(xué)性能數(shù)據(jù)及對應(yīng)的詳細流場細節(jié)，從而能在數(shù)天時間內(nèi)完成方案設(shè)計與性能分析工作，加速渦軸發(fā)動機設(shè)計過程，如圖2所示。這對于不斷提升渦軸發(fā)動機整機性能、減少設(shè)計反復(fù)具有重大意義。

圖2 軸流-離心組合壓氣機氣動仿真

同時，基于非線性諧波法和全通道非定常RAN-S求解法的非定常仿真方法，是研究進氣畸變影響、葉輪機械端區(qū)流動、轉(zhuǎn)靜子干涉時序效應(yīng)等復(fù)雜流動問題的有效工具，可在試驗前獲得葉輪機械內(nèi)部流場的時域與頻域信息，從而指導(dǎo)試驗與測試方案，降低了型號研制對于試驗的依賴度和試驗需求，特別是減少了畸變試驗和喘振邊界摸底試驗等高風(fēng)險試驗的次數(shù)，提升了試驗的成功率。

另外，基于RAN-S的流熱耦合仿真廣泛應(yīng)用于渦輪葉片冷卻方案設(shè)計，并與葉片冷卻效果試驗相結(jié)合，解決小尺寸冷卻試驗測試布置困難、測量區(qū)域和精度有限的難題，獲得小尺寸葉片在真實環(huán)境下較為準(zhǔn)確的溫度分布，確保渦輪設(shè)計與試驗的一次成功。

近年來，以離散大渦模擬（DES）和大渦模擬（LES）的仿真方法在低雷諾數(shù)渦輪葉片、高壓比單級壓氣機流場分析和端區(qū)復(fù)雜流動分析和損失控制中也得到了應(yīng)用，解決小尺寸葉輪機流場測量困難的問題，在一定程度上可以替代葉柵試驗，優(yōu)化葉輪機設(shè)計，提高葉輪機全包線性能。

在仿真工作的支持下，動研所目前已經(jīng)具備壓比25∶1以上的組合壓氣機設(shè)計能力，以及單級膨脹比4.5∶1跨聲渦輪和進口總溫1900K量級、總效率92%的小尺寸冷卻渦輪設(shè)計能力，未來將繼續(xù)向負荷更大、效率更高、結(jié)構(gòu)更緊湊的葉輪機械高精度性能仿真方向發(fā)展。

聯(lián)合/耦合氣動及換熱仿真

渦軸發(fā)動機由于自身結(jié)構(gòu)緊湊和裝機位置特殊等原因，需要對氣動及換熱問題進行綜合考量，才能較為準(zhǔn)確地評估出設(shè)計方案的有效性及合理性，為發(fā)動機的安全運轉(zhuǎn)及性能優(yōu)化提供定量數(shù)據(jù)支撐。

空氣系統(tǒng)

渦軸發(fā)動機空氣系統(tǒng)與發(fā)動機軸向力、強度壽命以及變形等關(guān)系密切，而這些又與發(fā)動機的安全運行緊密相關(guān)。因此，只有準(zhǔn)確地仿真出轉(zhuǎn)子周圍的流動及換熱情況，對轉(zhuǎn)子件強度、壽命、變形進行準(zhǔn)確分析，才能為發(fā)動機的安全性設(shè)計提供保障。

由于渦軸發(fā)動機空氣流量小、溫度梯度大，冷氣沿程溫升很高，流動與傳熱熱效應(yīng)非常顯著。為滿足工程分析精度要求，必須進行空氣流動與零件傳熱耦合仿真分析。以一個渦軸發(fā)動機燃氣渦輪二級導(dǎo)向葉片及內(nèi)封嚴(yán)結(jié)構(gòu)的流熱耦合結(jié)果為例，其中考慮了冷氣沿程的沖擊換熱、肋片強化冷卻、空氣預(yù)旋、盤腔風(fēng)阻、篦齒射流等一系列復(fù)雜流動與換熱現(xiàn)象的相互耦合影響。空氣系統(tǒng)仿真在準(zhǔn)確提供轉(zhuǎn)子盤腔內(nèi)部氣流流動信息的同時，還有助于更精確地進行零部件冷卻效果、溫度場分析。這對提升渦軸發(fā)動機性能、壽命，保障渦軸發(fā)動機安全運轉(zhuǎn)提供了有力支撐。

整機聯(lián)合氣動仿真

渦軸發(fā)動機性能與進氣道性能和外界來流情況直接關(guān)聯(lián)，嚴(yán)重時發(fā)動機的壓比、效率和裕度均會大幅下降，甚至危及穩(wěn)定運行，這就需要通過各部件聯(lián)合氣動仿真和內(nèi)外流耦合仿真分別對不同因素進行評估[4]。為此，動研所開展了發(fā)動機進氣道—壓氣機—燃燒室—渦輪聯(lián)合氣動仿真研究工作，如圖3所示，可有效評估各部件對整機性能的影響程度，對各部件的耦合匹配意義重大。

圖3 渦軸發(fā)動機進氣道—壓氣機—燃燒室—渦輪聯(lián)合氣動仿真結(jié)果

直升機/發(fā)動機內(nèi)外流耦合仿真

為進一步探尋機體和直升機旋翼流動對發(fā)動機性能的影響規(guī)律以及發(fā)動機布局和進排氣相互干擾的影響，尤其是三發(fā)布局的直升機，需要進行內(nèi)外流耦合氣動仿真研究工作。目前，動研所通過采用多參考系模型，已開展了內(nèi)外流耦合流場仿真研究，如圖4所示，用以解決直升機的發(fā)動機進氣畸變、尾氣吸入等問題，避免飛發(fā)匹配與布局的重大調(diào)整，加速研制進程。

圖4 內(nèi)外流耦合流場仿真結(jié)果

起動過程全尺寸燃燒室三維點火仿真

動研所在基于RAN-S燃燒仿真[5]研究的基礎(chǔ)上，實現(xiàn)了燃燒室內(nèi)的速度、溫度和濃度等物理量的高精度預(yù)測，滿足了燃燒室多輪次、快速迭代工程設(shè)計的需要。同時，動研所采用的“低維模型+RAN-S冷態(tài)流場”的點火概率分析方法，可以低成本實現(xiàn)燃燒室點熄火等非穩(wěn)態(tài)工況的高精度仿真，進而實現(xiàn)對點火邊界的預(yù)測，已為多款發(fā)動機燃燒室的順利研發(fā)提供了技術(shù)保障，避免設(shè)計與試驗驗證的不斷反復(fù)，大大提升了燃燒室設(shè)計精準(zhǔn)度，減少了試驗次數(shù)。同時，針對下一代高性能燃燒室的研發(fā)需求，動研所還開展了LES燃燒仿真研究和基于可實現(xiàn)k-ε（RKE）湍流模型、火焰面生成流形（FGM）模型的低成本數(shù)值仿真研究，可分別實現(xiàn)精細化燃燒仿真和高精度點火概率仿真，如圖5所示，進一步為燃燒室設(shè)計的一次成功提供了保障。

圖5 RAN-S與LES仿真結(jié)果對比

多相流仿真

渦軸發(fā)動機在復(fù)雜環(huán)境使用中，其真實流動還包括氣/固、氣/液、氣/液/固等多相流動。典型的包括粒子分離器沙石分離、發(fā)動機進氣結(jié)冰、吞水以及滑油潤滑系統(tǒng)內(nèi)部的空氣/滑油等多相流。動研所的研究經(jīng)驗表明：高保真或高精度的多相流仿真工作在加速設(shè)計方案迭代、改進和優(yōu)化結(jié)構(gòu)、確保試驗驗證的成功等方面作用巨大。

粒子分離器氣/固兩相流仿真

粒子分離器的主要功能是避免沙石被渦軸發(fā)動機吸入，以減少對發(fā)動機葉片的撞擊、磨蝕和對渦輪冷卻通道的堵塞，避免發(fā)動機性能的衰減和故障的發(fā)生，在渦軸發(fā)動機上被廣泛應(yīng)用。目前針對粒子分離器的穩(wěn)態(tài)三維流動問題，動研所通過采用基于拉格朗日粒子追蹤原理的仿真方法，開展了大量的氣固兩相耦合仿真，模擬了粒子分離器內(nèi)的分叉流動現(xiàn)象（如圖6所示），并結(jié)合試驗進行了沙石模型的修正，提高了仿真精度，其分離效率的仿真預(yù)測精度已在5%以內(nèi)，基本滿足了工程需要。在此基礎(chǔ)上，動研所已實現(xiàn)了粗沙分離效率不低于90%、細沙分離效率不低于80%、總壓損失小于2.5%的粒子分離器的研發(fā)與驗證，為渦軸發(fā)動機的安全運轉(zhuǎn)提供了切實保障。

圖6 粒子分離器內(nèi)部沙粒運動軌跡仿真

防/結(jié)冰設(shè)計的氣/液及液固相變仿真

結(jié)冰作為航空器上的常見現(xiàn)象，給飛行器的安全帶來了極大隱患。目前，防/結(jié)冰設(shè)計分析涉及兩相流、相變、傳熱傳質(zhì)過程，技術(shù)難度大；冰風(fēng)洞試驗成本高、周期長，給防冰系統(tǒng)設(shè)計及驗證造成了較大困難。動研所基于FLUENT軟件二次開發(fā)的防冰/結(jié)冰仿真分析軟件，能夠?qū)^冷水滴撞擊特性、防冰熱平衡以及靜止部件的結(jié)冰冰形進行較為準(zhǔn)確的模擬，如圖7所示，極大地提升了航空發(fā)動機防冰系統(tǒng)的設(shè)計能力。通過防冰/結(jié)冰仿真分析，不僅提高了防冰系統(tǒng)設(shè)計成功率，還解決了民用渦軸發(fā)動機研發(fā)中遇到的結(jié)冰關(guān)鍵點選取等技術(shù)瓶頸，促進了民用渦軸發(fā)動機適航取證工作的順利開展，大大降低了防/結(jié)冰試驗周期和成本，為渦軸發(fā)動機的安全性設(shè)計提供了堅實的保障。

圖7 進氣支板結(jié)冰冰形仿真與試驗對比

滑油系統(tǒng)多相流動仿真

渦軸發(fā)動機滑油系統(tǒng)的性能對整機的可靠性、壽命、安全性都意義重大，其內(nèi)部涉及油氣兩相、液體撞擊固體壁面、旋轉(zhuǎn)湍流等復(fù)雜流動狀態(tài)，是典型的多相流動問題。目前，動研所采用體積函數(shù)法（VOF）兩相計算模型，開展了滑油潤滑系統(tǒng)的流動仿真分析，包括典型旋轉(zhuǎn)流動下的環(huán)下潤滑結(jié)構(gòu)內(nèi)部油/氣兩相流的流動狀態(tài)仿真，如圖8所示，仿真對收油效率的預(yù)測和試驗結(jié)果偏差基本達到工程需要。并且相較試驗而言，數(shù)值仿真周期短，比試驗研究更具靈活性，可以對試驗難以測量的內(nèi)部復(fù)雜流場做出預(yù)測，獲取內(nèi)部的油氣兩相流流動狀態(tài)和特性，為滑油系統(tǒng)管路、通道結(jié)構(gòu)優(yōu)化改進以及滑油消耗量精準(zhǔn)預(yù)估提供了依據(jù)，對提高滑油系統(tǒng)收油效率、保證發(fā)動機的安全運轉(zhuǎn)均有重要意義。

圖8 滑油系統(tǒng)空氣場幾何模型—仿真—試驗結(jié)果圖

結(jié)束語

綜合分析現(xiàn)有技術(shù)的研究現(xiàn)狀和發(fā)展趨勢可以預(yù)見：以機器學(xué)習(xí)為代表的低維度實時仿真技術(shù)，可以大大加快發(fā)動機氣動設(shè)計迭代速度；以LES大渦模擬為代表的高保真三維仿真技術(shù)，可大幅提高仿真精度并統(tǒng)一發(fā)動機氣動湍流模型。上述方法均對發(fā)動機詳細設(shè)計意義重大，如應(yīng)用得當(dāng)，則有望在計算量適中且精度較高的湍流仿真模型與方法的建立、飛發(fā)一體化耦合仿真、高精度流動分離預(yù)測及氣固熱耦合仿真等渦軸發(fā)動機性能仿真未來研究

中發(fā)揮巨大作用。