• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于Gauss羽流模型低階預(yù)估旋流燃燒室中守恒標(biāo)量的空間分布*

    2023-10-20 01:11:10吳子恒張世紅王柏森
    關(guān)鍵詞:模型

    吳子恒, 張 弛, 張世紅, 王柏森

    (北京航空航天大學(xué) 航空發(fā)動(dòng)機(jī)研究院 航空發(fā)動(dòng)機(jī)氣動(dòng)熱力國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室, 北京 102206)

    0 引 言

    隨著先進(jìn)燃?xì)廨啓C(jī)性能提升,其對(duì)燃燒室高效穩(wěn)定燃燒和高出口溫度場(chǎng)品質(zhì)的需求也在增加[1],對(duì)渦輪葉片可靠性與壽命也具有重要意義.燃燒室內(nèi)復(fù)雜的旋流流動(dòng)和油氣混合是影響燃燒高溫區(qū)和出口溫度分布的關(guān)鍵因素,空氣通過(guò)頭部旋流器進(jìn)入到火焰筒中,與燃料混合燃燒后產(chǎn)生局部高溫區(qū),其在旋流作用下在下游發(fā)生遷移,影響燃燒室出口溫度分布的不均勻性.為了改善燃燒室出口溫度場(chǎng)品質(zhì),通過(guò)頭部旋流流動(dòng)和燃料-空氣混合方法來(lái)控制旋流燃燒局部高溫區(qū)的生成具有重要意義.混合分?jǐn)?shù)是表征燃料-空氣混合效果的守恒標(biāo)量[2],其三維空間分布對(duì)于燃?xì)廨啓C(jī)燃燒室混合設(shè)計(jì)具有重要的指導(dǎo)意義.混合分?jǐn)?shù)也是燃燒室湍流燃燒建模的關(guān)鍵參考標(biāo)量,如概率密度函數(shù)輸運(yùn)(PDF)[3]、條件矩(CMC)[4]、層流小火焰等模型[5]通過(guò)混合分?jǐn)?shù)構(gòu)建小尺度混合模型或基于溫度、組分與混合分?jǐn)?shù)的關(guān)聯(lián)函數(shù)封閉湍流-火焰的相互作用.國(guó)內(nèi)外學(xué)者利用非接觸光學(xué)測(cè)量了簡(jiǎn)單典型燃燒器的混合分?jǐn)?shù)[6-8],然而目前還無(wú)法準(zhǔn)確測(cè)量復(fù)雜非典型燃燒器的混合分?jǐn)?shù)場(chǎng),且通過(guò)實(shí)驗(yàn)指導(dǎo)燃燒室設(shè)計(jì)成本高、效率低,測(cè)試方案數(shù)量有限,參數(shù)顆粒度跨度大,對(duì)方案的精細(xì)化尋優(yōu)難度極大[9].燃?xì)廨啓C(jī)燃燒室等復(fù)雜燃燒器的混合分?jǐn)?shù)場(chǎng)空間分布往往通過(guò)三維數(shù)值仿真獲得.Kim等[10]采用Reynolds時(shí)均(RANS)數(shù)值模擬方法探究了不同入口速度和當(dāng)量比條件下直葉片和扭曲葉片旋流對(duì)燃料-空氣混合和燃燒特性的影響.Lv等[11]通過(guò)大渦模擬(LES)方法研究了燃料-空氣混合和分級(jí)旋流火焰燃燒穩(wěn)定性.Steinhausen等[12]消耗了1 800萬(wàn)核時(shí)(計(jì)算時(shí)處理器核心數(shù)量乘以小時(shí)數(shù)),通過(guò)直接數(shù)值模擬(DNS)的方法,采用詳細(xì)化學(xué)反應(yīng)機(jī)理,研究了湍流火焰-壁面的相互作用.目前三維數(shù)值模擬方法廣泛應(yīng)用于燃燒室混合、燃燒高溫區(qū)、出口溫度分布評(píng)估和氣動(dòng)設(shè)計(jì)的過(guò)程中,但對(duì)于復(fù)雜燃燒器,存在算不動(dòng)、算不起的問(wèn)題,這對(duì)燃燒室方案設(shè)計(jì)的迭代效率產(chǎn)生很大的影響.低階預(yù)估模型可以規(guī)避上述劣勢(shì),目前鮮有對(duì)于燃燒室混合分?jǐn)?shù)場(chǎng)低階預(yù)估的研究.發(fā)展混合分?jǐn)?shù)場(chǎng)的低階預(yù)估模型,以加速燃料-空氣混合策略的評(píng)估和燃燒室參數(shù)化設(shè)計(jì)過(guò)程,具有很高的工程應(yīng)用價(jià)值.

    Gauss羽流模型可以用來(lái)描述點(diǎn)源在主流來(lái)流中的擴(kuò)散過(guò)程[13].Sánchez-Sosa等[14]將Gauss羽流模型應(yīng)用于移動(dòng)機(jī)器人室內(nèi)氣體源定位.李萬(wàn)莉[15]將Gauss羽流模型應(yīng)用到了復(fù)雜地形下天然氣泄露這一場(chǎng)景,對(duì)模型進(jìn)行了校驗(yàn).目前,鮮有研究將Gauss羽流模型應(yīng)用于燃燒室中湍流燃燒混合過(guò)程的預(yù)測(cè).傳統(tǒng)Gauss羽流模型沒有考慮徑向?qū)α鲗?duì)混合過(guò)程的影響,且不能應(yīng)用于存在旋流來(lái)流、多點(diǎn)源和限制域壁面的場(chǎng)景.

    本文從混合分?jǐn)?shù)控制方程出發(fā),同時(shí)考慮了標(biāo)量的對(duì)流與擴(kuò)散,推導(dǎo)出考慮對(duì)流影響的Gauss煙團(tuán)方程,從而得到了考慮對(duì)流影響的Gauss羽流模型.進(jìn)一步考慮了旋流來(lái)流和多點(diǎn)源場(chǎng)景,發(fā)展了鏡像反射模型來(lái)模擬壁面-羽流霧相互作用,并引入相關(guān)修正來(lái)確保質(zhì)量守恒,將新推導(dǎo)的Gauss羽流模型應(yīng)用于甲烷旋流燃燒室混合分?jǐn)?shù)場(chǎng)的低階預(yù)測(cè).對(duì)甲烷旋流燃燒室開展三維數(shù)值仿真計(jì)算,獲得了數(shù)值收斂的混合分?jǐn)?shù)場(chǎng)數(shù)據(jù)庫(kù).采用最小二乘法對(duì)模型參數(shù)進(jìn)行優(yōu)化,在寬范圍條件下驗(yàn)證了模型的的預(yù)測(cè)精度.

    1 理 論 模 型

    1.1 模型推導(dǎo)

    圖1所示為點(diǎn)源釋放氣體在直流空氣來(lái)流速度U=(u,v,w)條件下、開放空間中的發(fā)展示意圖,針對(duì)該過(guò)程開展模型推導(dǎo),對(duì)Gauss羽流模型進(jìn)行二次發(fā)展.點(diǎn)源來(lái)流在空氣來(lái)流輸運(yùn)過(guò)程中,守恒標(biāo)量混合分?jǐn)?shù)ξ的控制方程可以寫為

    圖1 點(diǎn)源釋放氣體在直流空氣來(lái)流中的發(fā)展示意圖Fig. 1 Schematic diagram of the development of a point source releasing gas in a straight air stream

    (1)

    其中ρ為密度,D為擴(kuò)散系數(shù)[16],在Gauss羽流模型中兩者均假設(shè)為定值,根據(jù)來(lái)流工況確定.將ρξ記為C(x,y,z,t),引入三個(gè)方向速度分量u,v,w,同樣滿足Gauss羽流模型假設(shè),u,v,w為定值.式(1)可以寫為

    (2)

    (3)

    其中F為Fourier變換,i為虛數(shù)單位.

    進(jìn)一步解析式(3)得[18]

    (4)

    通過(guò)Fourier逆變換獲得

    C(x,y,z,t)=F-1[F(λ1,λ2,λ3,t)]=

    (5)

    基于式(4),定義函數(shù)G(λ1,λ2,λ3,t):

    (6)

    設(shè)函數(shù)g(x,y,z,t)的Fourier變換為G(λ1,λ2,λ3,t),即

    F[g(x,y,z,t)]=G(λ1,λ2,λ3,t).

    (7)

    在t=0 s時(shí)刻,上述定義的Fourier變換函數(shù)為

    F(λ1,λ2,λ3,0)=F[C(x,y,z,0)].

    (8)

    根據(jù)Fourier卷積定理[19],對(duì)于任意函數(shù)f和h有

    f*h=F-1[F(f)·F(h)].

    (9)

    其中*表示卷積運(yùn)算.

    根據(jù)式(5)—(9),可以推導(dǎo)獲得

    C(x,y,z,t)=C(x,y,z,0)*g(x,y,z,t).

    (10)

    為獲得C(x,y,z,t)的解析式,則需要求解函數(shù)g(x,y,z,t):

    g(x,y,z,t)=F-1[G(λ1,λ2,λ3,t)]=

    (11)

    由式(10)和(11)可得

    (12)

    (13)

    式(13)便是考慮對(duì)流影響的Gauss煙團(tuán)方程(Gaussian puff equation).

    需要說(shuō)明的是:式(13)表征點(diǎn)源瞬態(tài)釋放的氣體的空間分布.若要得到式(13)的等效三維穩(wěn)態(tài)方程,即考慮徑向?qū)α饔绊懙腉auss羽流模型,可將點(diǎn)源上游來(lái)流空間離散為厚度為dx的薄板,并對(duì)每個(gè)薄板采用Lagrange觀點(diǎn)進(jìn)行分析[21].圖2為將來(lái)流離散為薄板掃略點(diǎn)源的示意圖.

    圖2 將來(lái)流離散為薄板掃略點(diǎn)源的示意圖Fig. 2 Schematic diagram of discretizing the flow stream into thin sheets sweeping the point source

    參考式(2),離散薄板上的二維輸運(yùn)方程為

    (14)

    參考式(3)—(13)的推導(dǎo)過(guò)程可得到

    (15)

    (16)

    式(16)可表征不同軸向位置薄板上的濃度分布.每個(gè)薄板厚度為dx,將式(16)除以dx,即可得到三維穩(wěn)態(tài)空間分布,即考慮徑向?qū)α鞯腉auss羽流模型

    (17)

    Gauss羽流模型的有效范圍為x軸下游x≥x0.基于C=ρξ,同時(shí)考慮點(diǎn)源處混合分?jǐn)?shù)ξ上限為1,將式(17)進(jìn)一步變換為

    (18)

    其中a,b為待定模型參數(shù).

    如圖3所示,當(dāng)考慮旋流來(lái)流的影響時(shí),該點(diǎn)源位置y和z方向上的速度分量v和w分別記為

    圖3 點(diǎn)源釋放氣體在旋流空氣來(lái)流中的發(fā)展簡(jiǎn)圖Fig. 3 Schematic diagram of the development of a point source releasing gas in a swirling air stream

    v=vc,

    (19)

    w=wc.

    (20)

    由于燃燒室實(shí)際的燃燒室燃料噴射方式存在一個(gè)噴嘴具有多個(gè)燃料噴射點(diǎn)的情況,因此需要發(fā)展多點(diǎn)源模型,以圖4所示三點(diǎn)源為例.通過(guò)“像源法”[15]建模,引入鏡像羽流(mirror-image plume)[21]對(duì)應(yīng)的虛擬等效“點(diǎn)源”,對(duì)壁面-羽流相互作用建模.

    圖4 多點(diǎn)源釋放氣體在旋流空氣來(lái)流中的發(fā)展和鏡像反射模型示意圖Fig. 4 Schematic diagram of the development of multi-point sources releasing gas in a swirling air stream and the mirror reflection model

    若引入N個(gè)點(diǎn)源,則點(diǎn)源k處y和z方向上的速度分量分別為

    vk=vc,k,

    (21)

    wk=wc,k.

    (22)

    由此可獲得由多點(diǎn)源表征的混合分?jǐn)?shù)場(chǎng)預(yù)估模型

    (23)

    其中(x0,y0+pk,z0+qk)為點(diǎn)源的位置.引入多點(diǎn)源的同時(shí)會(huì)產(chǎn)生質(zhì)量守恒的問(wèn)題[21],為了模型在預(yù)估混合分?jǐn)?shù)場(chǎng)的同時(shí)保證質(zhì)量守恒,并將混合分?jǐn)?shù)限定于0≤ξ(x,y,z)≤1,可將式(23)修正為

    (24)

    m為待定參數(shù).由于燃燒室限制域的影響,最小混合分?jǐn)?shù)為非0定值,進(jìn)一步引入燃燒室下游混合分?jǐn)?shù)邊界ξfinal,得到

    (25)

    α,β為待定參數(shù),x′為基于燃燒室長(zhǎng)度的無(wú)量綱軸向距離.

    2 三維數(shù)值仿真數(shù)據(jù)

    2.1 模型燃燒室

    圖5所示為甲烷旋流燃燒室?guī)缀谓Y(jié)構(gòu)圖,其基于GT2500燃?xì)廨啓C(jī)燃燒室頭部旋流方案構(gòu)建.火焰筒長(zhǎng)為345 mm,直徑為77 mm,旋流器分內(nèi)外兩級(jí),內(nèi)旋流葉片角度和數(shù)目分別為35°和12,旋流數(shù)為0.8.外旋流葉片角度和數(shù)目分別為30°和18,旋流數(shù)為0.6.噴嘴前端有噴孔24個(gè),每個(gè)噴孔直徑為2 mm.表1所示為燃燒室工況條件,進(jìn)口壓力為1 554 800 Pa,進(jìn)口空氣流量為0.625 kg/s,通過(guò)改變?nèi)剂狭髁縼?lái)改變?nèi)紵耶?dāng)量比.分別在燃燒室當(dāng)量比為0.8,0.9,1.0,1.1的工況下(class 1)開展Reynolds平均計(jì)算,獲得三維數(shù)值仿真數(shù)據(jù)庫(kù),校驗(yàn)?zāi)P蛥?shù)并驗(yàn)證模型準(zhǔn)確性;在當(dāng)量比為0.7,0.75,1.2三個(gè)外推工況條件下(class 2)驗(yàn)證模型的寬范圍適用性.

    表1 工況條件

    圖5 甲烷旋流燃燒室?guī)缀谓Y(jié)構(gòu)Fig. 5 The configuration of the methane swirl combustor

    在三維數(shù)值仿真中[22],采用realizablek-ε模型來(lái)描述帶旋流特征的湍流流動(dòng),壁面處湍流通過(guò)可伸縮壁面函數(shù)實(shí)現(xiàn),同時(shí)采用壁面絕熱、無(wú)滑移假設(shè).燃燒室進(jìn)口空氣速度脈動(dòng)和噴嘴出口速度脈動(dòng)設(shè)置為平均速度的5%,密度通過(guò)理想氣體狀態(tài)方程求解,黏度值計(jì)算遵循Sutherland定律,分別定義Schmidt數(shù)和Prandtl數(shù)將濃度和溫度擴(kuò)散系數(shù)與氣體黏度關(guān)聯(lián).湍流燃燒模型采用FGM(flamelet generated manifold)模型,化學(xué)反應(yīng)機(jī)理為GRI3.0,壓力-速度耦合算法采用SIMPLE算法.采用有限體積法求解,離散精度為二階.

    2.2 網(wǎng)格無(wú)關(guān)性驗(yàn)證

    針對(duì)上述燃燒室流體域,分別生成了384萬(wàn)、520萬(wàn)、890萬(wàn)和1 153萬(wàn)四套離散網(wǎng)格.在當(dāng)量比φ為1.1的工況下開展計(jì)算,取燃燒室中軸線上的CH4質(zhì)量分?jǐn)?shù)進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證,如圖6所示.從圖中可以看出,當(dāng)網(wǎng)格數(shù)量大于890萬(wàn)時(shí),燃燒室軸線上CH4質(zhì)量分?jǐn)?shù)分布不再發(fā)生改變,滿足網(wǎng)格無(wú)關(guān)性的基本要求.因此選用890萬(wàn)的網(wǎng)格生成三維數(shù)值仿真數(shù)據(jù),作為低階模型的驗(yàn)證依據(jù).

    圖6 中軸線上甲烷質(zhì)量分?jǐn)?shù)的分布Fig. 6 Mass fractions of CH4 on the centerline of the axial direction

    3 模 型 驗(yàn) 證

    3.1 點(diǎn)源數(shù)量選取

    從低階預(yù)估模型的推導(dǎo)過(guò)程看,點(diǎn)源數(shù)量與噴嘴的噴射點(diǎn)數(shù)量有一定關(guān)聯(lián),此外考慮壁面反射的影響,通過(guò)鏡面反射模型引入虛擬等效“點(diǎn)源”.為了保證混合分?jǐn)?shù)場(chǎng)的預(yù)估精度,點(diǎn)源數(shù)量可能需要大于等于噴嘴實(shí)際噴射點(diǎn)的數(shù)量.為進(jìn)一步測(cè)試點(diǎn)源數(shù)量對(duì)預(yù)估結(jié)果的敏感性,以當(dāng)量比1.1的工況為例,分別設(shè)置點(diǎn)源數(shù)量為N=12,24,36開展預(yù)估模型的測(cè)試,并選取軸向距離x=30 mm和x=60 mm的低階預(yù)估與三維數(shù)值仿真云圖對(duì)比,如圖7和圖8所示.在x=30 mm處,N=12,24,36的預(yù)測(cè)結(jié)果差異不大.隨著軸向距離增大,如在x=60 mm處,N=12導(dǎo)致模型預(yù)測(cè)的混合分?jǐn)?shù)場(chǎng)出現(xiàn)不平滑、“失真”的現(xiàn)象,而N=24和36的結(jié)果基本一致,可以看到混合分?jǐn)?shù)較高的區(qū)域均為光滑的同心圓,與CFD結(jié)果相符.

    (a) 混合分?jǐn)?shù)模型預(yù)估結(jié)果(N=12) (b) 混合分?jǐn)?shù)模型預(yù)估結(jié)果(N=24)(a) Model results of the mixture fraction(N=12)(b) Model results of the mixture fraction(N=24)

    (a) 混合分?jǐn)?shù)模型預(yù)估結(jié)果(N=12)(b) 混合分?jǐn)?shù)模型預(yù)估結(jié)果(N=24)(a) Model results of the mixture fraction(N=12)(b) Model results of the mixture fraction(N=24)

    如圖9所示,取x=60 mm處的徑向分布曲線,對(duì)不同N的預(yù)估結(jié)果進(jìn)行定量比較,N=24和36的預(yù)估結(jié)果幾乎一致,與N=12的結(jié)果略有差異.綜上所述,將點(diǎn)源數(shù)量N=24應(yīng)用于該案例的模型驗(yàn)證是合理的.

    圖9 不同點(diǎn)源數(shù)量的低階模型預(yù)估結(jié)果與三維數(shù)值模擬結(jié)果對(duì)比 (x=60 mm,徑向分布, N=12, N=24, N=36)Fig. 9 Comparison of low-order model prediction results of different numbers of point sources and the 3D numerical simulation results (x=60 mm, radial distribution,N=12, N=24, N=36)

    3.2 參數(shù)尋優(yōu)以及確定

    將式(25)中的N取為24,根據(jù)當(dāng)量比0.8,0.9,1.0,1.1四個(gè)工況的混合分?jǐn)?shù)場(chǎng),使用最小二乘法對(duì)模型的參數(shù)進(jìn)行尋優(yōu),確定的參數(shù)α,β,a,b,m,vc,k,wc,k,pk,qk的取值見附錄.

    3.3 預(yù)估結(jié)果與三維數(shù)值模擬結(jié)果對(duì)比分析

    3.3.1 參考工況對(duì)比分析

    取中截面云圖,將模型的預(yù)估結(jié)果與對(duì)應(yīng)的三維數(shù)值模擬結(jié)果進(jìn)行對(duì)比,如圖10所示.模型預(yù)估的云圖與CFD結(jié)果具有一致性,混合分?jǐn)?shù)從噴嘴開始向下游沿著頭部壁面不斷減小直至消失,整體呈“V”形區(qū)域,下游回流區(qū)同一軸向位置混合分?jǐn)?shù)最高點(diǎn)不在中軸線上.本文發(fā)展的低階模型考慮旋流來(lái)流和徑向?qū)α鞯挠绊?可以較好地預(yù)測(cè)上述特征.

    (a) 中截面混合分?jǐn)?shù)低階模型預(yù)估結(jié)果(φ=0.8)(b) 中截面混合分?jǐn)?shù)CFD計(jì)算結(jié)果(φ=0.8) (a) Model results of the mixture fraction (b) CFD results of the mixture fraction at the central plane (φ=0.8) at the central plane (φ=0.8)

    將模型預(yù)估與三維數(shù)值模擬結(jié)果進(jìn)行定量對(duì)比,如圖11所示.從結(jié)果中可以看出,低階模型可以準(zhǔn)確預(yù)估混合分?jǐn)?shù)沿中軸線的分布,混合分?jǐn)?shù)在x<0.05 m的區(qū)域迅速降低,之后基本保持很小的定值不變,且隨頭部當(dāng)量比的升高而升高.

    (a) 模型與CFD結(jié)果對(duì)比(φ=0.8)(b) 模型與CFD結(jié)果對(duì)比(φ=0.9)(a) Comparison of the results from the model and (b) Comparison of the results from the model and the CFD(φ=0.8) the CFD(φ=0.9)

    圖12所示為x=0.05 m處的混合分?jǐn)?shù)徑向分布,模型的計(jì)算結(jié)果和CFD的計(jì)算結(jié)果均是中部較低,兩側(cè)較高,且隨著頭部當(dāng)量比的減小,峰值逐漸增大,徑向邊界值逐漸降低.

    (a) 模型與CFD結(jié)果對(duì)比(φ=0.8)(b) 模型與CFD結(jié)果對(duì)比(φ=0.9)(a) Comparison of the results from the model and (b) Comparison of the results from the model and the CFD(φ=0.8) the CFD(φ=0.9)

    同樣由于低階模型考慮了旋流來(lái)流對(duì)標(biāo)量遷移的影響,模型可以準(zhǔn)確地預(yù)估徑向混合分?jǐn)?shù)“雙峰”和“低谷”值及其所在位置.

    3.3.2 外推工況對(duì)比分析

    在參考工況下將預(yù)估模型的參數(shù)進(jìn)行校驗(yàn),并與三維數(shù)值模擬結(jié)果對(duì)比分析后,為驗(yàn)證寬范圍下模型的預(yù)測(cè)精度,對(duì)當(dāng)量比為0.7,0.75,1.2的數(shù)據(jù)開展驗(yàn)證分析.如圖13所示,預(yù)估模型依然可以獲得混合分?jǐn)?shù)場(chǎng)的“V”形分布特征.圖14和圖15軸向和徑向混合分?jǐn)?shù)變化趨勢(shì)可以被預(yù)估模型準(zhǔn)確捕捉.與參考工況相比,模型預(yù)估誤差和最大誤差的位置相似,并對(duì)定量對(duì)比數(shù)據(jù)進(jìn)行分析,模型的預(yù)估結(jié)果和三維數(shù)值仿真結(jié)果的平均誤差為11.19%.

    (a) 中截面混合分?jǐn)?shù)模型計(jì)算結(jié)果(φ=0.7) (b) 中截面混合分?jǐn)?shù)CFD計(jì)算結(jié)果(φ=0.7)(a) Model results of the mixture fraction (b) CFD results of the mixture fraction at the central plane (φ=0.7) at the central plane (φ=0.7)

    (a) 模型與CFD結(jié)果對(duì)比(φ=0.7)(a) Comparison of the results from the model and the CFD (φ=0.7)

    3.4 計(jì)算效率

    圖16所示為低階預(yù)估模型與三維數(shù)值仿真模型在參數(shù)化設(shè)計(jì)過(guò)程中的耗時(shí)對(duì)比.本研究需要四個(gè)工況下的三維數(shù)值仿真數(shù)據(jù)對(duì)所發(fā)展的低階預(yù)估模型將進(jìn)行模型參數(shù)校驗(yàn),三維數(shù)值仿真每個(gè)算例開銷為120核并行迭代2 800步,3 h之后結(jié)果收斂,之后低階模型可以應(yīng)用于外推工況的混合分?jǐn)?shù)場(chǎng)預(yù)估,每算例僅耗時(shí)1 s,故低階模型的開銷基本由模型校驗(yàn)所需的算例數(shù)量決定.若使用三維數(shù)值仿真開展參數(shù)化研究,計(jì)算開銷隨算例數(shù)量呈線性增長(zhǎng).本研究發(fā)展的低階模型可以極大加快混合分?jǐn)?shù)的計(jì)算效率,加快燃燒室混合方案評(píng)估和參數(shù)化研究.

    圖16 兩種方法耗時(shí)隨算例個(gè)數(shù)的對(duì)比Fig. 16 Comparison of the computational time costs with the number of cases

    3.5 低階預(yù)估模型誤差

    為了更好地評(píng)估模型誤差,我們將模型外推應(yīng)用于DLR旋流燃燒室混合分?jǐn)?shù)場(chǎng)預(yù)估[23].如圖17所示,與高精度實(shí)驗(yàn)數(shù)據(jù)(experiment)進(jìn)行定量比對(duì),低階預(yù)估模型(model)的平均誤差為10.02%,且低階模型預(yù)估與三維數(shù)值仿真(CFD)預(yù)測(cè)[23]效果相當(dāng).綜上所述,該低階預(yù)估模型在追求效率的同時(shí),仍然具有較好的預(yù)測(cè)精度,符合工程應(yīng)用的需求.

    4 結(jié)論與展望

    1) 本研究對(duì)傳統(tǒng)的Gauss羽流模型進(jìn)行了二次開發(fā),考慮了徑向?qū)α?、旋流?lái)流和壁面的影響,建立了能夠應(yīng)用于甲烷旋流燃燒室的混合分?jǐn)?shù)場(chǎng)低階預(yù)估模型.

    2) 基于模型燃燒室當(dāng)量比0.8,0.9,1.0,1.1四個(gè)工況的三維數(shù)值模擬數(shù)據(jù)庫(kù),通過(guò)最小二乘法對(duì)低階預(yù)估模型開展參數(shù)尋優(yōu),驗(yàn)證了混合分?jǐn)?shù)場(chǎng)的預(yù)估精度,在外推工況條件下驗(yàn)證了模型的適用性.

    3) 由于最小二乘法對(duì)于初值的敏感性, 且容易陷入局部最優(yōu)解, 會(huì)影響混合分?jǐn)?shù)場(chǎng)局部細(xì)節(jié)的預(yù)測(cè), 接下來(lái)可以發(fā)展更精確的模型參數(shù)尋優(yōu)方法[24],進(jìn)一步增強(qiáng)模型的寬適用性及其在燃燒室設(shè)計(jì)應(yīng)用時(shí)的魯棒性.

    附 錄

    參數(shù)α,β,a,b,m的值如表A1所示.

    表A1 參數(shù)α,β,a,b,m的值

    參數(shù)vc,k,wc,k,pk,qk的值如表A2—A5所示.

    表A2 參數(shù)vc,k的值

    表A3 參數(shù)wc,k的值

    表A4 參數(shù)pk的值

    表A5 參數(shù)qk的值

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    国产一区二区在线观看日韩| 久久精品久久精品一区二区三区| 韩国av在线不卡| 黑人巨大精品欧美一区二区蜜桃 | 免费不卡的大黄色大毛片视频在线观看| 日韩成人av中文字幕在线观看| 97精品久久久久久久久久精品| 在线观看人妻少妇| 色网站视频免费| 又大又黄又爽视频免费| 看免费av毛片| a级毛色黄片| 天天影视国产精品| 波野结衣二区三区在线| 少妇人妻 视频| 亚洲 欧美一区二区三区| 99热全是精品| 人妻一区二区av| 草草在线视频免费看| 久久青草综合色| 国产精品一二三区在线看| 日韩欧美一区视频在线观看| 纵有疾风起免费观看全集完整版| 欧美国产精品一级二级三级| 日韩一区二区三区影片| 午夜福利,免费看| 亚洲精品乱码久久久久久按摩| xxx大片免费视频| 国产在视频线精品| 国产亚洲av片在线观看秒播厂| 国产精品一区二区在线观看99| 国产一区二区在线观看日韩| 亚洲国产最新在线播放| 在线观看一区二区三区激情| 寂寞人妻少妇视频99o| 亚洲av.av天堂| 美女福利国产在线| 新久久久久国产一级毛片| 成人亚洲精品一区在线观看| 亚洲国产av影院在线观看| 一区二区三区四区激情视频| 国产成人欧美| 黄色一级大片看看| 国产一区二区在线观看av| 日日摸夜夜添夜夜爱| 久热久热在线精品观看| 波野结衣二区三区在线| 久久人人爽人人爽人人片va| 51国产日韩欧美| 亚洲丝袜综合中文字幕| 免费人妻精品一区二区三区视频| 9色porny在线观看| 香蕉精品网在线| 久久久久精品人妻al黑| 日韩中字成人| 日韩av免费高清视频| 久久 成人 亚洲| 久久99热6这里只有精品| 国产精品麻豆人妻色哟哟久久| 乱码一卡2卡4卡精品| 久热这里只有精品99| 男女午夜视频在线观看 | 日韩av免费高清视频| 国产亚洲最大av| 国产成人精品婷婷| 中文字幕av电影在线播放| 欧美日韩av久久| 涩涩av久久男人的天堂| 久久久久久久久久成人| 日本黄大片高清| 精品国产乱码久久久久久小说| 激情五月婷婷亚洲| 亚洲,一卡二卡三卡| 免费观看性生交大片5| 亚洲高清免费不卡视频| 2021少妇久久久久久久久久久| 日韩一区二区三区影片| av线在线观看网站| 精品人妻在线不人妻| 国产精品.久久久| 永久免费av网站大全| 日韩伦理黄色片| 国产成人精品婷婷| 久久青草综合色| 中文乱码字字幕精品一区二区三区| 国产乱来视频区| 久久精品国产鲁丝片午夜精品| 国产精品一区二区在线观看99| 亚洲av男天堂| 亚洲欧美清纯卡通| 国产av一区二区精品久久| 精品人妻偷拍中文字幕| av国产精品久久久久影院| 丰满饥渴人妻一区二区三| 国产日韩欧美亚洲二区| 国产精品麻豆人妻色哟哟久久| 巨乳人妻的诱惑在线观看| 另类亚洲欧美激情| 欧美日本中文国产一区发布| 国产有黄有色有爽视频| 91午夜精品亚洲一区二区三区| 日韩不卡一区二区三区视频在线| 高清av免费在线| 国产免费视频播放在线视频| a 毛片基地| 18禁观看日本| 少妇的逼水好多| 国产有黄有色有爽视频| 在线观看国产h片| 999精品在线视频| 久久久久国产网址| 菩萨蛮人人尽说江南好唐韦庄| 日韩欧美精品免费久久| 热99久久久久精品小说推荐| 国产一区二区三区av在线| av国产精品久久久久影院| av线在线观看网站| 男女边摸边吃奶| 久久久久人妻精品一区果冻| 日韩制服丝袜自拍偷拍| 亚洲av欧美aⅴ国产| 国产片内射在线| 国产亚洲一区二区精品| 国产男女超爽视频在线观看| www.熟女人妻精品国产 | 人人妻人人澡人人爽人人夜夜| 亚洲国产日韩一区二区| 美国免费a级毛片| 啦啦啦啦在线视频资源| 国产精品人妻久久久久久| a级毛片在线看网站| 免费播放大片免费观看视频在线观看| 国产激情久久老熟女| 一本大道久久a久久精品| 黄色视频在线播放观看不卡| 国产亚洲一区二区精品| 久久精品国产a三级三级三级| 青春草视频在线免费观看| 亚洲,欧美精品.| 大陆偷拍与自拍| 我的女老师完整版在线观看| 91国产中文字幕| 久久久亚洲精品成人影院| 日韩av在线免费看完整版不卡| 国产一区二区三区av在线| 亚洲av日韩在线播放| 汤姆久久久久久久影院中文字幕| 久久久久精品久久久久真实原创| 亚洲精品久久久久久婷婷小说| 国产永久视频网站| 香蕉国产在线看| 99国产综合亚洲精品| 18+在线观看网站| 在线观看人妻少妇| 国产午夜精品一二区理论片| 99国产精品免费福利视频| 永久网站在线| 精品一区在线观看国产| 韩国av在线不卡| 国产老妇伦熟女老妇高清| 国产国拍精品亚洲av在线观看| 国产 一区精品| 欧美xxⅹ黑人| 亚洲色图 男人天堂 中文字幕 | 日韩三级伦理在线观看| 七月丁香在线播放| 国产精品一区二区在线不卡| 丰满乱子伦码专区| 性色avwww在线观看| 春色校园在线视频观看| 五月天丁香电影| www日本在线高清视频| 欧美xxxx性猛交bbbb| 黑人巨大精品欧美一区二区蜜桃 | 99热网站在线观看| 久久久久久伊人网av| 亚洲国产欧美在线一区| 成人国产麻豆网| 男女国产视频网站| 丰满少妇做爰视频| 51国产日韩欧美| a 毛片基地| 亚洲精品一区蜜桃| 蜜桃在线观看..| 久久久精品免费免费高清| 男人爽女人下面视频在线观看| 精品99又大又爽又粗少妇毛片| 午夜免费观看性视频| 看十八女毛片水多多多| 免费观看a级毛片全部| 亚洲精品国产色婷婷电影| 久久午夜福利片| 国产爽快片一区二区三区| 人人妻人人澡人人爽人人夜夜| 高清av免费在线| 国产成人欧美| 午夜福利网站1000一区二区三区| 久久久久久久精品精品| 日韩欧美一区视频在线观看| 日日撸夜夜添| 欧美日韩亚洲高清精品| 亚洲精品日韩在线中文字幕| 国内精品宾馆在线| 亚洲经典国产精华液单| 青春草亚洲视频在线观看| 亚洲少妇的诱惑av| 欧美人与性动交α欧美精品济南到 | 欧美亚洲 丝袜 人妻 在线| www.色视频.com| 亚洲图色成人| 亚洲精华国产精华液的使用体验| av播播在线观看一区| 一区二区三区乱码不卡18| 国产成人午夜福利电影在线观看| 最近最新中文字幕免费大全7| 国产精品一二三区在线看| 国产av码专区亚洲av| 丰满迷人的少妇在线观看| 男女午夜视频在线观看 | 亚洲欧美一区二区三区国产| 欧美日韩综合久久久久久| 亚洲欧洲精品一区二区精品久久久 | 成人国语在线视频| 欧美3d第一页| 亚洲精品乱码久久久久久按摩| 久久精品国产综合久久久 | 国产精品久久久久久av不卡| 国产午夜精品一二区理论片| 一级片'在线观看视频| 欧美 日韩 精品 国产| 国产福利在线免费观看视频| 成人无遮挡网站| 久久av网站| 日韩三级伦理在线观看| 亚洲精品一区蜜桃| 99热国产这里只有精品6| 最近2019中文字幕mv第一页| 成人亚洲欧美一区二区av| 青春草亚洲视频在线观看| 国产精品麻豆人妻色哟哟久久| 免费观看a级毛片全部| 国产精品久久久久久精品电影小说| 不卡视频在线观看欧美| 国产视频首页在线观看| 久久狼人影院| 日韩三级伦理在线观看| 十八禁高潮呻吟视频| 免费人妻精品一区二区三区视频| 日韩 亚洲 欧美在线| 99国产综合亚洲精品| 午夜免费鲁丝| 亚洲,欧美精品.| 亚洲国产欧美在线一区| 久久久久久久国产电影| 国产熟女午夜一区二区三区| 国产日韩欧美视频二区| 天天操日日干夜夜撸| 美女内射精品一级片tv| 亚洲国产日韩一区二区| 欧美性感艳星| 97超碰精品成人国产| 亚洲国产欧美在线一区| 国产高清三级在线| 国产欧美日韩综合在线一区二区| 国产成人欧美| 国产黄色免费在线视频| 女性被躁到高潮视频| 亚洲国产av新网站| 国产成人欧美| 国产一区二区三区av在线| 欧美日韩av久久| 精品国产露脸久久av麻豆| 一区二区三区四区激情视频| 十八禁网站网址无遮挡| 久久久精品94久久精品| 亚洲三级黄色毛片| 在线观看三级黄色| 九九爱精品视频在线观看| 中文欧美无线码| 亚洲欧美中文字幕日韩二区| 国产成人精品无人区| 精品卡一卡二卡四卡免费| 在线精品无人区一区二区三| 国产精品女同一区二区软件| 久久亚洲国产成人精品v| 午夜激情久久久久久久| 国产xxxxx性猛交| 亚洲国产av影院在线观看| 激情视频va一区二区三区| 中文天堂在线官网| 久久久久网色| 寂寞人妻少妇视频99o| 好男人视频免费观看在线| 国产探花极品一区二区| 天堂8中文在线网| h视频一区二区三区| 最近最新中文字幕免费大全7| 中文字幕免费在线视频6| 日本与韩国留学比较| 午夜激情av网站| 寂寞人妻少妇视频99o| 卡戴珊不雅视频在线播放| 丝袜喷水一区| 亚洲五月色婷婷综合| 青春草国产在线视频| 亚洲少妇的诱惑av| 欧美日韩亚洲高清精品| 欧美精品高潮呻吟av久久| 99re6热这里在线精品视频| 人体艺术视频欧美日本| 汤姆久久久久久久影院中文字幕| 一区二区av电影网| 成人18禁高潮啪啪吃奶动态图| 亚洲精品国产色婷婷电影| 亚洲欧美中文字幕日韩二区| 国产精品不卡视频一区二区| 日本午夜av视频| 亚洲成人手机| 一级爰片在线观看| 久久久亚洲精品成人影院| 啦啦啦中文免费视频观看日本| 免费黄网站久久成人精品| 国产精品久久久久久久电影| 国产国拍精品亚洲av在线观看| 国产乱来视频区| 看十八女毛片水多多多| 各种免费的搞黄视频| 午夜激情av网站| 国产探花极品一区二区| 亚洲av.av天堂| 亚洲人与动物交配视频| 人妻一区二区av| 国产片内射在线| 日韩视频在线欧美| 精品午夜福利在线看| 极品少妇高潮喷水抽搐| 精品久久久久久电影网| 精品99又大又爽又粗少妇毛片| 如日韩欧美国产精品一区二区三区| 国产黄色视频一区二区在线观看| 久久av网站| 女人精品久久久久毛片| 日韩欧美精品免费久久| 亚洲欧美色中文字幕在线| 黄色一级大片看看| 国产成人精品福利久久| 国产精品久久久久久久电影| 久久99蜜桃精品久久| 菩萨蛮人人尽说江南好唐韦庄| 一本色道久久久久久精品综合| 久久精品久久精品一区二区三区| 在现免费观看毛片| 成人午夜精彩视频在线观看| 午夜激情久久久久久久| 国产高清三级在线| 在线观看三级黄色| 青春草视频在线免费观看| 色吧在线观看| 亚洲国产毛片av蜜桃av| 国产男人的电影天堂91| 男女无遮挡免费网站观看| 日本wwww免费看| 国产黄色视频一区二区在线观看| 国产免费视频播放在线视频| 久久精品久久久久久久性| 亚洲丝袜综合中文字幕| 亚洲色图 男人天堂 中文字幕 | 久久婷婷青草| 色哟哟·www| 国产日韩欧美在线精品| 亚洲综合色网址| 亚洲国产毛片av蜜桃av| 国产欧美另类精品又又久久亚洲欧美| 午夜影院在线不卡| 欧美精品国产亚洲| 青春草视频在线免费观看| 久久人人爽人人爽人人片va| 在线观看免费视频网站a站| 国产麻豆69| 建设人人有责人人尽责人人享有的| 国产午夜精品一二区理论片| 男人舔女人的私密视频| 宅男免费午夜| 成人亚洲欧美一区二区av| 亚洲国产欧美在线一区| 亚洲成人一二三区av| 热99久久久久精品小说推荐| 制服丝袜香蕉在线| xxx大片免费视频| 99九九在线精品视频| 国产亚洲一区二区精品| 午夜福利,免费看| 免费人妻精品一区二区三区视频| 在线观看美女被高潮喷水网站| 丝袜喷水一区| 亚洲精品乱久久久久久| 日韩精品免费视频一区二区三区 | 亚洲欧美日韩另类电影网站| 制服诱惑二区| 黑人巨大精品欧美一区二区蜜桃 | 国产精品国产av在线观看| 在线观看国产h片| 久久久久久久久久成人| 热re99久久国产66热| 99国产综合亚洲精品| 欧美 亚洲 国产 日韩一| 男女下面插进去视频免费观看 | 看非洲黑人一级黄片| 午夜精品国产一区二区电影| 精品一区二区免费观看| 精品一区二区三卡| 久久这里只有精品19| 999精品在线视频| 大话2 男鬼变身卡| 啦啦啦中文免费视频观看日本| h视频一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 国产爽快片一区二区三区| 欧美日韩精品成人综合77777| 欧美精品亚洲一区二区| 久久久久久久久久成人| 少妇的逼好多水| 美女脱内裤让男人舔精品视频| 亚洲国产av影院在线观看| 欧美成人午夜精品| 中文欧美无线码| 纯流量卡能插随身wifi吗| 最黄视频免费看| 国产爽快片一区二区三区| 桃花免费在线播放| 秋霞伦理黄片| 国产成人av激情在线播放| 性高湖久久久久久久久免费观看| 春色校园在线视频观看| 亚洲激情五月婷婷啪啪| 欧美日韩亚洲高清精品| 国产精品国产三级国产av玫瑰| 久久精品熟女亚洲av麻豆精品| 久久ye,这里只有精品| 亚洲,欧美精品.| 亚洲欧美日韩卡通动漫| 亚洲经典国产精华液单| 日本欧美国产在线视频| 日韩人妻精品一区2区三区| 麻豆乱淫一区二区| 日韩不卡一区二区三区视频在线| 国产又爽黄色视频| 久久久久人妻精品一区果冻| 女性生殖器流出的白浆| 99久国产av精品国产电影| 男的添女的下面高潮视频| 国产亚洲一区二区精品| 极品人妻少妇av视频| 日本欧美视频一区| 欧美性感艳星| 夫妻性生交免费视频一级片| 九草在线视频观看| 你懂的网址亚洲精品在线观看| 日本欧美国产在线视频| 免费高清在线观看视频在线观看| 熟妇人妻不卡中文字幕| 一级片'在线观看视频| 蜜桃在线观看..| 亚洲欧美日韩卡通动漫| 制服丝袜香蕉在线| 国产亚洲一区二区精品| 搡女人真爽免费视频火全软件| 国产一区有黄有色的免费视频| 多毛熟女@视频| 亚洲在久久综合| 日韩成人av中文字幕在线观看| 国产av一区二区精品久久| 精品一区二区三卡| 欧美日韩精品成人综合77777| 亚洲精品第二区| 久久久久久久久久人人人人人人| 久久精品久久精品一区二区三区| 国产又色又爽无遮挡免| 少妇的丰满在线观看| 99久久人妻综合| 国产极品粉嫩免费观看在线| 下体分泌物呈黄色| 性高湖久久久久久久久免费观看| 久久热在线av| 久久精品熟女亚洲av麻豆精品| 少妇的丰满在线观看| 男女午夜视频在线观看 | 一本久久精品| 亚洲av电影在线进入| 久久青草综合色| 99精国产麻豆久久婷婷| 九色亚洲精品在线播放| 99香蕉大伊视频| 日本91视频免费播放| 狂野欧美激情性xxxx在线观看| 久久久久国产网址| 久久国内精品自在自线图片| 国产精品成人在线| 国产亚洲精品第一综合不卡 | 国产一区二区激情短视频 | 男人添女人高潮全过程视频| 黄色怎么调成土黄色| 精品国产露脸久久av麻豆| 侵犯人妻中文字幕一二三四区| 最新中文字幕久久久久| 久久狼人影院| 国产av国产精品国产| 午夜福利在线观看免费完整高清在| 国产精品国产三级国产专区5o| 精品一区二区三区四区五区乱码 | 最近最新中文字幕免费大全7| 久久久欧美国产精品| 国产精品久久久久久av不卡| 日韩一本色道免费dvd| 在线亚洲精品国产二区图片欧美| 国产国拍精品亚洲av在线观看| 不卡视频在线观看欧美| 色婷婷久久久亚洲欧美| 国产 一区精品| tube8黄色片| 青春草亚洲视频在线观看| 久久久久久久久久人人人人人人| 秋霞在线观看毛片| 亚洲av男天堂| 国产成人精品一,二区| 少妇的丰满在线观看| 少妇熟女欧美另类| 中文字幕免费在线视频6| 香蕉国产在线看| 欧美3d第一页| 五月伊人婷婷丁香| 欧美成人精品欧美一级黄| 街头女战士在线观看网站| 国产乱来视频区| 如日韩欧美国产精品一区二区三区| 亚洲精品456在线播放app| 少妇被粗大猛烈的视频| 日韩制服丝袜自拍偷拍| 亚洲av男天堂| 日本欧美视频一区| 国产精品三级大全| 欧美老熟妇乱子伦牲交| 天堂8中文在线网| 国产有黄有色有爽视频| 国产无遮挡羞羞视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 免费看光身美女| av女优亚洲男人天堂| 26uuu在线亚洲综合色| 亚洲国产看品久久| 久久97久久精品| 成人毛片60女人毛片免费| 如何舔出高潮| 男女边吃奶边做爰视频| 夫妻性生交免费视频一级片| 建设人人有责人人尽责人人享有的| 成人国语在线视频| 日本黄色日本黄色录像| 国产老妇伦熟女老妇高清| 日本午夜av视频| 午夜福利乱码中文字幕| 9191精品国产免费久久| 街头女战士在线观看网站| 免费少妇av软件| 啦啦啦视频在线资源免费观看| 一级片免费观看大全| 国产在线一区二区三区精| 青春草亚洲视频在线观看| 免费观看av网站的网址| 亚洲国产成人一精品久久久| 免费av中文字幕在线| 欧美日韩综合久久久久久| 热re99久久精品国产66热6| 精品午夜福利在线看| 精品酒店卫生间| 久久人人爽人人片av| 国产成人免费观看mmmm| 人人妻人人爽人人添夜夜欢视频| 在线 av 中文字幕| 国产极品粉嫩免费观看在线| kizo精华| 热99国产精品久久久久久7| 69精品国产乱码久久久| 少妇被粗大的猛进出69影院 | 免费黄频网站在线观看国产| 观看美女的网站| 久久久久久久久久久久大奶| 最新的欧美精品一区二区| 日韩伦理黄色片| 高清毛片免费看| 菩萨蛮人人尽说江南好唐韦庄| 丝瓜视频免费看黄片| 一级片'在线观看视频| 午夜免费观看性视频| 亚洲欧洲精品一区二区精品久久久 | 高清在线视频一区二区三区| 国产一区亚洲一区在线观看| 久久毛片免费看一区二区三区| 久久青草综合色| 你懂的网址亚洲精品在线观看| 欧美xxxx性猛交bbbb| 国产熟女午夜一区二区三区| 丁香六月天网| 国产精品欧美亚洲77777| 深夜精品福利| 一本大道久久a久久精品| 日本午夜av视频| 久久久欧美国产精品| videosex国产| 欧美丝袜亚洲另类| √禁漫天堂资源中文www| 国产熟女欧美一区二区| 午夜精品国产一区二区电影| 国产日韩一区二区三区精品不卡| 国产女主播在线喷水免费视频网站|