王方，王煜棟，姜勝利，陳軍，唐軍，徐華勝，李象遠，邢競文，高東碩，金捷

1.北京航空航天大學(xué) 能源與動力工程學(xué)院，北京 100191 2.北京應(yīng)用物理與計算數(shù)學(xué)研究所，北京 100088 3.中物院高性能數(shù)值模擬軟件中心，北京 100088 4.中國燃?xì)鉁u輪研究院，成都 610599 5.四川大學(xué) 化學(xué)工程學(xué)院，成都 610065

現(xiàn)代航空發(fā)動機研究需要應(yīng)用高精度數(shù)值模擬技術(shù)。燃燒室?guī)缀谓Y(jié)構(gòu)復(fù)雜，其中的湍流燃燒現(xiàn)象存在強烈非線性關(guān)系，實驗研究難度大、耗時長、風(fēng)險高，迫切需要高保真、高精度的數(shù)值模擬技術(shù)。擁有先進發(fā)動機制造技術(shù)的國家在燃燒室數(shù)值模擬軟件的研發(fā)中都投入了大量資金與科研力量，開發(fā)了可用于燃燒室設(shè)計的數(shù)值模擬軟件，如美國的國家燃燒仿真代碼(NCC)等。近年來，中國越來越重視發(fā)展航空發(fā)動機數(shù)值模擬技術(shù)和軟件?；谥袊売嬎銠C算力的優(yōu)勢，針對燃燒室復(fù)雜結(jié)構(gòu)及其內(nèi)部湍流流動、化學(xué)反應(yīng)及燃料液滴蒸發(fā)相互非線性耦合關(guān)系，開發(fā)大規(guī)模并行、適用于航空發(fā)動機燃燒室等復(fù)雜結(jié)構(gòu)的兩相湍流燃燒仿真軟件具有重要意義。

北京航空航天大學(xué)航空發(fā)動機數(shù)值仿真研究中心(北航仿真中心)與英國帝國理工William Jones教授合作在Boffin程序基礎(chǔ)上開發(fā)出基于LES-TPDF(Large-Eddy Simulation-Transported Probability Density Function)方法的兩相湍流燃燒數(shù)值模擬軟件AECSC 2.0(Aero Engine Combustor Simulation Code 2.0)，已被應(yīng)用于多種場合的兩相湍流燃燒計算并驗證了其準(zhǔn)確性[1-3]。北航仿真中心與北京應(yīng)用物理與計算數(shù)學(xué)研究所、中國工程物理研究院高性能數(shù)值模擬軟件中心合作研發(fā)的AECSC-JASMIN，將大渦模擬(LES)與概率密度函數(shù)輸運方程湍流燃燒模擬方法(TPDF)和JASMIN(并行自適應(yīng)結(jié)構(gòu)網(wǎng)格應(yīng)用支撐軟件)大規(guī)模并行計算框架、浸沒邊界方法及正交均勻結(jié)構(gòu)網(wǎng)格耦合，形成幾何高保真湍流燃燒高分辨率的仿真軟件。其串行部分由AECSC 2.0軟件發(fā)展而來，二者在兩相湍流燃燒模擬中算法一致；網(wǎng)格生成與并行框架采用北京應(yīng)用物理與計算數(shù)學(xué)研究所研制的面向結(jié)構(gòu)網(wǎng)格應(yīng)用的JASMIN框架。JASMIN框架面向現(xiàn)代高性能計算機體系結(jié)構(gòu)高效并行計算，由底至上采用并行自適應(yīng)支撐層、數(shù)值共性層、應(yīng)用接口層的分層結(jié)構(gòu)，已成功應(yīng)用于武器物理、多介質(zhì)輻射流體力學(xué)[4]、激光聚變[5]、計算電磁學(xué)和位錯動力學(xué)[6-7]等多種領(lǐng)域。

高分辨率、高精度的燃燒室數(shù)值模擬是燃燒仿真的熱點之一。莫妲等[8]采用ICEMCFD(The Integrated Computer Engineering and Manufacturing Code for Computational Fluid Dynamics)劃分423萬網(wǎng)格對三旋流燃燒室進行模擬，出口溫度分布與實驗值符合較好。丁勇能等[9]采用ICEMCFD繪制1.7億網(wǎng)格對重型燃?xì)廨啓C的單筒燃燒室進行數(shù)值模擬，得到分辨率較高、準(zhǔn)確性較好的燃燒場模擬結(jié)果。對于燃燒室復(fù)雜幾何結(jié)構(gòu)，采用ICEMCFD等軟件劃分大規(guī)模網(wǎng)格難度較高，耗費大量時間，且網(wǎng)格質(zhì)量難以保證。大規(guī)模高質(zhì)量網(wǎng)格的劃分是燃燒室高分辨率仿真的難點之一。采用大規(guī)模網(wǎng)格使計算量增大，需要性能極高的計算設(shè)備以及大量的計算時間[10]，提高非定常數(shù)值模擬的可擴展性和計算效率十分必要[11]。在燃燒數(shù)值模擬中，湍流流動計算方法呈現(xiàn)以大渦模擬為主的趨勢。Zhou等[12]采用動態(tài)亞網(wǎng)格模型和大渦模擬對兩級反向旋流燃燒室單頭部構(gòu)型進行數(shù)值模擬，得出高分辨率、高精度的計算結(jié)果。

本文首先介紹AECSC-JASMIN軟件研發(fā)工作，其次采用AECSC-JASMIN軟件模擬射流火焰算例、支板火焰穩(wěn)定器算例和單頭部燃燒室算例，用實驗數(shù)據(jù)檢驗軟件。在射流火焰算例中模擬3種不同的工況，對比多個截面的溫度脈動均方根以及時均溫度的實驗數(shù)據(jù)和模擬結(jié)果檢驗軟件湍流燃燒模擬是否準(zhǔn)確。在兩相支板火焰穩(wěn)定器算例中，對比新研發(fā)軟件AECSC-JASMIN、商用軟件Fluent、舊版軟件AECSC 2.0模擬得到的出口時均溫度分布數(shù)據(jù)及實驗數(shù)據(jù)，檢驗新軟件在較復(fù)雜模型兩相湍流燃燒模擬中是否可行。之后，基于對真實燃燒室單頭部模型的模擬檢驗AECSC-JASMIN軟件在航空發(fā)動機燃燒室數(shù)值模擬應(yīng)用的可行性，探索燃燒室模擬技術(shù)的技術(shù)途徑。

1 計算模型與方法

1.1 湍流和湍流燃燒模型

大渦模擬通過濾波函數(shù)對流動變量過濾，尺度較大的變量直接求解，尺度較小的變量使用亞網(wǎng)格模型模擬。模擬精度高于雷諾平均方法而計算量小于直接數(shù)值模擬(DNS)，在湍流燃燒模擬的工程應(yīng)用中較為實用。軟件采用動態(tài)Smagorinsky-Lilly亞網(wǎng)格應(yīng)力模型進行湍流場計算、隨機場概率密度函數(shù)輸運方程方法計算湍流燃燒相互作用，可以耦合詳細(xì)化學(xué)反應(yīng)機理。

對于湍流和化學(xué)反應(yīng)都重要的燃燒過程，概率密度函數(shù)輸運方程方法能夠比較準(zhǔn)確模擬湍流與化學(xué)反應(yīng)強烈相互作用，結(jié)合詳細(xì)化學(xué)反應(yīng)機理，在不假設(shè)組分空間低維流形時封閉平均化學(xué)反應(yīng)速率，對化學(xué)反應(yīng)機理的剛性要求較其他湍流燃燒模型更低。

Jaberi等[13]在大渦模擬中應(yīng)用等擴散系數(shù)假設(shè)，將Favre濾波后的隨機場輸運方程表示為

(1)

(2)

最后經(jīng)大渦濾波后得到

(3)

AECSC軟件采用了二階中心差分的空間離散格式和Crank-Nicholson的時間離散格式。軟件中還采用了隱式低馬赫數(shù)方程、壓力平滑算法、半隱式壓力方程SIMPLE算法等。

1.2 網(wǎng)格生成及并行計算

基于JASMIN并行框架重構(gòu)AECSC程序，采用分層架構(gòu)設(shè)計，將數(shù)值計算與并行方法分離。采用構(gòu)件模型，定制數(shù)值構(gòu)件并組裝計算流程以規(guī)范數(shù)值計算子程序。重新設(shè)計的AECSC-JASMIN架構(gòu)如圖1所示。

圖1 AECSC-JASMIN軟件架構(gòu)圖

AECSC-JASMIN加入了正交結(jié)構(gòu)網(wǎng)格的自動生成算法，通過指定計算域并導(dǎo)入幾何模型，自動構(gòu)建正交結(jié)構(gòu)網(wǎng)格。數(shù)據(jù)結(jié)構(gòu)方面，替換原AECSC中module型數(shù)組以適應(yīng)高性能計算體系結(jié)構(gòu)，提高并行效率。改造原軟件中數(shù)值計算模塊，將內(nèi)存調(diào)度、輸入輸出、數(shù)據(jù)通信功能從數(shù)值計算子程序中分離。

浸沒邊界方法(IBM)由Peskin[14-16]最早提出，用于計算流固耦合中物理邊界與網(wǎng)格邊界不一致時的流固耦合問題。采用正交結(jié)構(gòu)網(wǎng)格，極大程度上減少了網(wǎng)格劃分和加密的工作量。浸沒邊界方法在流體網(wǎng)格的邊界處添加δ函數(shù)作為強迫外力源項實現(xiàn)無滑移邊界條件。AECSC-JASMIN程序中設(shè)計了針對半隱式速度壓強耦合求解的浸沒邊界處理方法并在源程序中增加邊界條件處理模塊，支撐程序使用框架自動生成的正交結(jié)構(gòu)網(wǎng)格，該網(wǎng)格包含固體與流體，整體呈立方體形，通過網(wǎng)格點的變量“tflag”的值區(qū)分流場區(qū)域和固壁區(qū)域。以射流火焰算例為例，其網(wǎng)格外觀如圖2所示，其他算例的網(wǎng)格與之相似。

圖2 射流火焰算例網(wǎng)格

AECSC-JASMIN中還重構(gòu)了輸入輸出模塊，使輸入文件格式清晰。輸出模塊中采用了JASMIN框架的JaVis數(shù)據(jù)輸出器，提高海量數(shù)據(jù)實時輸出性能。

2 算例模擬結(jié)果及討論

為測試采用第1節(jié)算法的自研軟件AECSC-JASMIN的準(zhǔn)確性及實際應(yīng)用的可靠性，按照氣相、兩相湍流燃燒測試的順序，分別選取射流火焰Flame-D、Flame-E、Flame-F算例，兩相支板穩(wěn)焰算例進行測試，并將模擬結(jié)果分別與前人實驗結(jié)果或已驗證正確的模擬結(jié)果進行比對。

2.1 射流火焰算例

為驗證AECSC-JASMIN版本軟件對于湍流燃燒計算的準(zhǔn)確性，選取Sandia實驗室的湍流射流火焰實驗[17]的Flame-D、Flame-E、Flame-F作為檢驗算例，并將計算溫度場與Barlow和Frank[18]的實驗溫度測量結(jié)果以及舊版軟件AECSC 2.0的模擬結(jié)果進行對比。

按照Sandia實驗室的實驗數(shù)據(jù)，射流火焰噴射燃燃料的內(nèi)噴口直徑為7.2 mm，噴射高溫導(dǎo)流氣體的外噴口直徑為18.2 mm，內(nèi)外噴管壁厚均為0.25 mm。

圖2展示了射流火焰算例所使用的全場尺寸均勻的630萬網(wǎng)格。考慮到本算例中網(wǎng)格尺寸為1.33 mm，大于實驗所用噴管壁厚，因此算例中直接在進口的管壁處設(shè)置一層速度為0的網(wǎng)格，不單獨設(shè)置管壁的厚度。本算例噴嘴結(jié)構(gòu)如圖3所示，內(nèi)噴口噴射溫度為293 K、比例為1∶3的甲烷燃料與空氣的混合氣，外噴口噴射溫度為1 880 K(Flame-D、Flame-E)或1 860 K(Flame-F)的空氣，高溫空氣將直接引燃內(nèi)噴口噴射的燃料。Flame-D、Flame-E、Flame-F的雷諾數(shù)分別為22 400、33 600、44 800，射流速度分布如圖4所示。

圖3 射流火焰實驗裝置簡化結(jié)構(gòu)圖

圖4 Flame-D、Flame-E和Flame-F的速度分布

在射流火焰算例的湍流燃燒模擬中，3個算例計算所用網(wǎng)格皆為630萬正方體網(wǎng)格。在基本保證計算精度的條件下，為了盡量加快計算速度，采用2個隨機場計算燃燒反應(yīng)。每個算例使用192個核心計算至2萬時間步，在每個時間步中CFL數(shù)(收斂條件判斷數(shù))取0.2～0.3計算得到當(dāng)前時間步長。射流火焰算例中，平均每個算例需消耗132 225核小時(并行計算總核數(shù)與以小時為單位的計算時間之積)計算得到，每千步約需6 611核小時。

圖5、圖6展示了3種射流火焰算例計算至最終2萬時間步時的瞬態(tài)速度云圖和溫度云圖。

圖5 Flame-D、Flame-E和Flame-F的速度云圖

圖6 Flame-D、Flame-E和Flame-F的溫度云圖

在計算達到穩(wěn)定燃燒狀態(tài)后，開始統(tǒng)計溫度數(shù)據(jù)，直至計算結(jié)束，得到時均溫度以及溫度脈動數(shù)據(jù)。根據(jù)實驗數(shù)據(jù)分別選取3種火焰中部的4個截面，對AECSC-JASMIN模擬結(jié)果的平均溫度場，作平均溫度分布曲線，并與實驗值及AECSC 2.0的計算值對比。取樣的4個截面Z/D分別為7.5、15.0、30.0、45.0，其中Z為射流方向坐標(biāo)，D為內(nèi)噴口直徑。取以上4個截面并對比模擬結(jié)果中的時均溫度分布以及溫度脈動分布。

從圖7～圖9中可以看出，AECSC-JASMIN計算得到的時均溫度場分布與實驗值及AECSC 2.0計算值基本相符。下面繼續(xù)對模擬中與實驗中射流火焰算例Flame-D、Flame-E、Flame-F的溫度脈動分布進行對比。

圖7 AECSC-JASMIN、AECSC 2.0模擬Flame-D算例的時均溫度分布與實驗對比

圖8 AECSC-JASMIN、AECSC 2.0模擬Flame-E算例的時均溫度分布與實驗對比

圖9 AECSC-JASMIN、AECSC 2.0模擬Flame-F算例的時均溫度分布與實驗對比

由圖10～圖12可知，AECSC-JASMIN軟件模擬結(jié)果中的溫度脈動與AECSC 2.0計算值、實驗值基本相符。表明AECSC-JASMIN對湍流燃燒時均溫度場和溫度脈動場的模擬具有較好的準(zhǔn)確性和較高的精度，能夠準(zhǔn)確地模擬湍流流動與化學(xué)反應(yīng)的相互作用。

圖10 AECSC-JASMIN、AECSC 2.0模擬Flame-D算例的溫度脈動分布與實驗對比

圖11 AECSC-JASMIN、AECSC 2.0模擬Flame-E算例的溫度脈動分布與實驗對比

圖12 AECSC-JASMIN、AECSC 2.0模擬Flame-F算例的溫度脈動分布與實驗對比

經(jīng)驗證，AECSC-JASMIN軟件對湍流燃燒溫度場的模擬較為準(zhǔn)確，下文將模擬支板火焰穩(wěn)定器算例以進一步驗證該軟件對于較復(fù)雜模型兩相湍流燃燒模擬的精度。

2.2 支板火焰穩(wěn)定器算例

支板火焰穩(wěn)定器是加力燃燒室中常用的穩(wěn)焰結(jié)構(gòu)。采用圖13所示的凹腔支板火焰穩(wěn)定器作為加力燃燒室的簡化模型，并在該模型的基礎(chǔ)上進行兩相湍流燃燒的模擬。

圖13 凹腔支板火焰穩(wěn)定器結(jié)構(gòu)

凹腔支板火焰穩(wěn)定器采用流線型翼面設(shè)計，由3段圓弧構(gòu)成外緣輪廓，整體長度108 mm，最大寬度30 mm，尾緣寬度16 mm。在支板弧度最大位置開深度為7 mm、寬度為12 mm的凹腔。凹腔后端面與中軸線夾角30°。支板兩相湍流燃燒算例使用RP3航空煤油燃料，計算工況為設(shè)計工況，油氣比為0.007?？諝膺M口(計算域坐標(biāo)X=0處)為馬赫數(shù)Ma=0.2、溫度1 173 K、壓強101 325 Pa的高溫高速空氣，初始全場壓強為101 325 Pa，通過高溫空氣進行熱射流點火。實驗時，噴油流道尺寸為700 mm×150 mm×150 mm，由3個供油點同時供油，實驗照片如圖14所示。

圖14 支板火焰穩(wěn)定器噴霧燃燒實驗照片

數(shù)值模擬時，將支板展向流道尺寸縮小為原尺寸的1/3，相應(yīng)地采用單個供油點，以減少計算量。計算時，支板展向邊界采用對稱邊界條件。支板上下兩側(cè)為固壁，前部邊界為入口，后部邊界為出口。

本算例網(wǎng)格總數(shù)270萬，全場網(wǎng)格尺寸為1.25 mm。網(wǎng)格由軟件自動生成，大規(guī)模并行計算的負(fù)載均衡由JASMIN框架實現(xiàn)。計算至4萬時間步共需8 320核小時，每1 000步需208核小時。

圖15為在AECSC-JASMIN軟件中采用四步機理模擬支板火焰穩(wěn)定器某穩(wěn)定燃燒工況得到的瞬態(tài)溫度分布云圖。四步機理模擬凹腔支板火焰穩(wěn)定器兩相湍流燃燒時，蒸發(fā)后的煤油在支板前緣有燃燒放熱現(xiàn)象。支板的凹腔部分有明顯的高溫區(qū)，存在駐定火焰，尾緣部分火焰溫度最高。這證明支板火焰穩(wěn)定器的凹腔和尾緣都對湍流火焰有穩(wěn)定作用，在本文工況下，煤油燃燒放熱過程主要在支板尾緣進行。

圖15 AECSC-JASMIN計算得到的支板火焰穩(wěn)定器噴霧燃燒瞬態(tài)溫度云圖

將AECSC-JASMIN采用烷烴四步機理計算得到的出口溫度分別與AECSC 2.0烷烴四步機理的計算結(jié)果、Fluent擴散火焰面的計算結(jié)果及實驗數(shù)據(jù)進行對比，結(jié)果如圖16所示。其中，AECSC 2.0采用貼體結(jié)構(gòu)網(wǎng)格及粒子反彈的壁面條件，其他參數(shù)及反應(yīng)機理與AECSC-JASMIN一致。對比可以發(fā)現(xiàn)，AECSC-JASMIN算例結(jié)果與AECSC 2.0計算結(jié)果基本一致。

在出口截面沿計算域Y軸方向(圖13～圖15所示平面中垂直于流向的方向)統(tǒng)計時均溫度沿Y坐標(biāo)的分布。

結(jié)合圖16、表1、表2，AECSC-JASMIN軟件模擬得到的出口溫度分布計算結(jié)果與采用貼體網(wǎng)格的舊版軟件模擬結(jié)果一致并且與實驗值相符。

表1 不同軟件計算的出口溫度分布及實驗數(shù)據(jù)

圖16 不同軟件計算得到的出口溫度分布曲線及實驗數(shù)據(jù)

表2中，AECSC-JASMIN平均相對誤差為5.514%，最大相對誤差為6.38%。AECSC 2.0平均相對誤差為2.964%，最大相對誤差為7.76%。AECSC-JASMIN平均誤差比AECSC 2.0更大，AECSC-JASMIN軟件所計算的算例網(wǎng)格全場尺寸一致，無局部加密算法，網(wǎng)格尺寸為1.25 mm，而AECSC 2.0算例的網(wǎng)格進行過局部加密，最小尺寸為0.2 mm，網(wǎng)格分辨率更高，在支板模型附近的計算精度高于AECSC-JASMIN。

表2 各軟件計算的出口溫度分布與實驗值的相對誤差

在表2中對比Fluent計算結(jié)果以及AECSC-JASMIN計算結(jié)果。Fluent火焰面模型模擬得到的出口溫度分布與實驗值相比，平均相對誤差為9.454%，最大相對誤差為13.01%。Fluent計算結(jié)果的平均相對誤差與最大相對誤差皆遠大于AECSC-JASMIN計算結(jié)果的平均相對誤差與最大相對誤差，由于Fluent采用了大渦模擬和火焰面模型，而AECSC-JASMIN采用了LES-TPDF方法，誤差對比結(jié)果證明了該模型的優(yōu)越性。AECSC-JASMIN軟件中采用自動構(gòu)建正交結(jié)構(gòu)網(wǎng)格，并在其中應(yīng)用浸沒邊界方法以解決復(fù)雜結(jié)構(gòu)網(wǎng)格生成問題，實際測試表明其計算精度仍較高。由于JASMIN框架自動生成的網(wǎng)格全場尺寸一致，在保證網(wǎng)格分辨率與AECSC 2.0對應(yīng)算例相同時，所需總網(wǎng)格量與AECSC 2.0相比更大，計算量更高。但由于軟件數(shù)據(jù)結(jié)構(gòu)重構(gòu)以及JASMIN并行框架的應(yīng)用使得AECSC-JASMIN軟件支持大規(guī)模并行計算，因此該軟件可以進行更加復(fù)雜的工程探索。

2.3 真實主燃燒室流場模擬

針對燃?xì)廨啓C燃燒室等復(fù)雜結(jié)構(gòu)中的兩相湍流燃燒仿真，AECSC-JASMIN軟件加入正交結(jié)構(gòu)網(wǎng)格自動生成算法、浸沒邊界處理方法以及高效率大規(guī)模并行框架。應(yīng)用該軟件采用1.6億網(wǎng)格對圖17(對部分圖的旋流器等結(jié)構(gòu)進行了高斯模糊處理)所示的某型航空發(fā)動機環(huán)形燃燒室的1/18構(gòu)型進行冷態(tài)流動模擬。

圖17 環(huán)形燃燒室1/18構(gòu)型圖

本算例中，進口位于擴壓器前部，即計算域中X=0平面，按無引氣時的試驗工況進行模擬，經(jīng)壓氣機壓縮過的空氣壓強P3為136 610 Pa，溫度T3為290.22 K，以V3=74.4 m/s的速度流入擴壓器，馬赫數(shù)Ma3為0.218，進口流量W3為2.07 kg/s。

本算例采用的網(wǎng)格與圖2中射流火焰算例網(wǎng)格圖相似，為全場尺寸均勻的立方體網(wǎng)格，網(wǎng)格總數(shù)為164 123 870，即1.6億。分別將壁面網(wǎng)格、流體區(qū)域網(wǎng)格標(biāo)記后可得燃燒室?guī)缀文Ｐ途W(wǎng)格如圖18、圖19所示。

圖18 燃燒室固體壁面網(wǎng)格

由圖18、圖19可知，軟件自動構(gòu)建的大規(guī)模(1.6億)網(wǎng)格可精確反應(yīng)燃燒室內(nèi)部細(xì)節(jié)結(jié)構(gòu)。

圖19 燃燒室流體區(qū)域網(wǎng)格(中央截面)

本算例中，軟件自動生成1.6億網(wǎng)格后采用480個核心計算至3.8萬步。關(guān)閉噴霧模塊與PDF模塊，僅計算流動場時，共需519 840核小時，每千步約需13 680核小時。

圖20展示了仿真得到的燃燒室內(nèi)部流線。單頭部燃燒室中流動復(fù)雜，旋流器附近的流場中含大量渦結(jié)構(gòu)。

圖20、圖21展示了采用AECSC-JASMIN模擬真實環(huán)形燃燒室1/18構(gòu)型內(nèi)部流動得到的流線圖和瞬態(tài)速度場，可知模擬結(jié)果流場宏觀特性與實驗相符。

圖20 燃燒室算例的流線圖

圖21 燃燒室算例的瞬態(tài)速度分布

表3展示了主燃燒室單頭部模型流場部分參數(shù)的實驗值與AECSC-JASMIN軟件的模擬值，由表3可知實驗測量的總壓恢復(fù)系數(shù)σ34等關(guān)鍵流場數(shù)值和軟件模擬值一致。

表3 實驗值與模擬值的流場參數(shù)對比

3 結(jié) 論

介紹了基于LES-TPDF方法及浸沒邊界法的自研兩相湍流燃燒軟件AECSC-JASMIN的湍流燃燒相關(guān)算法以及軟件架構(gòu)。之后應(yīng)用該軟件計算射流、支板、燃燒室算例，將計算結(jié)果與其他軟件模擬結(jié)果以及實驗數(shù)據(jù)對比，檢驗了AECSC-JASMIN軟件在兩相湍流燃燒數(shù)值模擬中的準(zhǔn)確性及工程應(yīng)用可行性。

數(shù)值計算算法方面，AECSC-JASMIN軟件采用LES-TPDF方法計算湍流燃燒，采用正交結(jié)構(gòu)網(wǎng)格和浸沒邊界方法簡化網(wǎng)格生成并處理邊界條件。軟件架構(gòu)方面，采用分層架構(gòu)設(shè)計，將數(shù)值計算與并行方法分離，采用構(gòu)件模型，定制數(shù)值構(gòu)件并組裝計算流程以規(guī)范數(shù)值計算子程序，采用JASMIN并行框架，滿足高分辨率高效大規(guī)模并行數(shù)值模擬的工程需要。

在軟件檢驗部分，應(yīng)用該軟件計算氣相射流火焰算例和兩相支板火焰穩(wěn)定器算例并將計算結(jié)果與其他仿真軟件的模擬結(jié)果以及實驗數(shù)據(jù)對比，發(fā)現(xiàn)該軟件在射流火焰模擬中多截面溫度脈動和平均溫度分布與AECSC 2.0模擬結(jié)果及實驗值相符，并且在兩相支板算例中該軟件模擬得到的出口溫度與其他數(shù)值模擬結(jié)果及實驗數(shù)據(jù)貼合較好，表明該軟件在湍流燃燒模擬精度較高，驗證了該軟件的正確性和準(zhǔn)確性。最后采用AECSC-JASMIN實現(xiàn)了復(fù)雜結(jié)構(gòu)環(huán)形燃燒室1/18構(gòu)型的1.6億網(wǎng)格大規(guī)模并行流動模擬，得到了燃燒室精細(xì)流場結(jié)構(gòu)，燃燒室流動宏觀特性與實驗結(jié)構(gòu)基本相符，流場中總壓損失系數(shù)等重要參數(shù)與實驗值一致。算例檢驗表明AECSC-JASMIN軟件對航空發(fā)動機燃燒室等復(fù)雜結(jié)構(gòu)內(nèi)的兩相湍流燃燒數(shù)值仿真具有較高的工程應(yīng)用價值，且隨著計算機算力的提高，在大規(guī)模并行湍流燃燒數(shù)值模擬中具有良好的發(fā)展前景，可以進一步探索。