• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

      亞跨聲速邊界層增長因子輸運(yùn)模式研究進(jìn)展

      2022-12-06 09:57:14徐家寬王玉軒張揚(yáng)喬磊劉建新白俊強(qiáng)
      航空學(xué)報(bào) 2022年11期
      關(guān)鍵詞:橫流壓力梯度雷諾數(shù)

      徐家寬,王玉軒,張揚(yáng),喬磊,劉建新,白俊強(qiáng),*

      1. 西北工業(yè)大學(xué) 航空學(xué)院,西安 710072

      2. 西安交通大學(xué) 機(jī)械強(qiáng)度與振動(dòng)國家重點(diǎn)實(shí)驗(yàn)室,西安 710049

      3. 天津大學(xué) 機(jī)械工程學(xué)院 高速空氣動(dòng)力學(xué)研究室,天津 300072

      20世紀(jì)50年代,Smith[1]和Van Ingen[2]等提出了二維邊界層的eN方法,該方法被廣泛應(yīng)用于航空航天領(lǐng)域。后來,Malik[3]、Arnal[4]、Cebeci[5]、蘇彩虹[6]、趙磊[7]等對(duì)該方法進(jìn)行了發(fā)展,使得其能夠應(yīng)用于三維邊界層,完善了基于線性穩(wěn)定性理論的eN方法。由于eN方法基于線性化理論,因而可以預(yù)測(cè)小擾動(dòng)引起的轉(zhuǎn)捩。而航空航天飛行器的飛行環(huán)境大多是低湍流度環(huán)境,因此eN方法比較適合航空航天飛行器繞流的轉(zhuǎn)捩預(yù)測(cè)。

      Saric等[8]總結(jié)了在三維邊界層中常見的轉(zhuǎn)捩類型。對(duì)于大氣中飛行的航空航天飛行器,所處的環(huán)境湍流度較低,來流比較平穩(wěn),失穩(wěn)主要由流向Tollmien-Schlichting(T-S)波轉(zhuǎn)捩、層流分離泡轉(zhuǎn)捩和橫流不穩(wěn)定波轉(zhuǎn)捩來主導(dǎo)。根據(jù)Obremski[9]和Klebanoff[10]等的研究結(jié)果,T-S波的線性發(fā)展區(qū)域幾乎能占據(jù)大約80%的層流區(qū),雖然橫流(Crossflow,CF)的非線性階段飽和區(qū)較長,但是大部分區(qū)域?yàn)榫€性發(fā)展區(qū)。因此,線性穩(wěn)定性理論能夠較好地預(yù)測(cè)T-S波轉(zhuǎn)捩和橫流轉(zhuǎn)捩。線性穩(wěn)定性理論計(jì)算得到擾動(dòng)增長因子的分布以后,還需要給定轉(zhuǎn)捩臨界N值才能判定轉(zhuǎn)捩。通常臨界N值則依靠經(jīng)驗(yàn)(試驗(yàn)等)給定,具有一定的經(jīng)驗(yàn)性。因此,eN方法是一種半經(jīng)驗(yàn)方法。

      Bippes[11]、Saric[12]、Kloker[13]和符松[14]等通過后掠平板和后掠機(jī)翼的試驗(yàn)和穩(wěn)定性計(jì)算,總結(jié)出如下規(guī)律:對(duì)于二維邊界層流動(dòng),一般當(dāng)來流湍流度Tu<(0.5%~1%),邊界層流動(dòng)由T-S波主導(dǎo)轉(zhuǎn)捩。在中、高來流湍流度下,邊界層流動(dòng)由條帶(Streaks)主導(dǎo)轉(zhuǎn)捩。對(duì)于三維橫流邊界層流動(dòng),在低湍流度環(huán)境下由定常橫流渦主導(dǎo)轉(zhuǎn)捩,而在中高湍流度環(huán)境下由非定常橫流渦主導(dǎo)轉(zhuǎn)捩。對(duì)于三維橫流邊界層,來流湍流度和壁面粗糙度是主要的擾動(dòng)源,分別誘導(dǎo)產(chǎn)生非定常橫流渦和定常橫流渦。Radeztsky等[15]進(jìn)一步研究了壁面粗糙度形狀、大小和位置對(duì)橫流定常渦感受性和失穩(wěn)的影響。低湍流度環(huán)境被認(rèn)為與實(shí)際飛行器飛行環(huán)境相似。在此環(huán)境下,盡管非定常橫流渦的增長率更大,但是三維橫流邊界層流動(dòng)主要由定常橫流渦主導(dǎo)轉(zhuǎn)捩。

      為了能夠在工程問題中應(yīng)用eN方法,美國國家航空航天局(NASA)的Malik[16]和Chang[17]、俄羅斯的Boiko等[18]、比利時(shí)的Pinna[19]、天津大學(xué)的黃章峰等[20]、西北工業(yè)大學(xué)的宋文萍等[21]和密歇根大學(xué)的史亞云等[22]都發(fā)展了雷諾平均方程(Reynolds Averaged Navier-Stokes,RANS)求解器耦合eN方法的分析程序。法國宇航研究中心(ONERA)的Perraud等[23]和Bégou等[24]、德國宇航中心的Krumbein[25]和瑞典國防研究局(FOI)的Eliasson等[26]基于大量的穩(wěn)定性分析結(jié)果構(gòu)造了eN數(shù)據(jù)庫,并利用該數(shù)據(jù)庫與CFD(Computational Fluid Dynamics)求解器進(jìn)行耦合求解,最終獲得擾動(dòng)增長因子。根據(jù)國內(nèi)外的科研人員大量的試驗(yàn)分析驗(yàn)證,在飛行器表面出現(xiàn)的T-S波轉(zhuǎn)捩、層流分離泡轉(zhuǎn)捩和橫流不穩(wěn)定性轉(zhuǎn)捩,基于線性穩(wěn)定性方法理論的eN方法都可以對(duì)其進(jìn)行預(yù)測(cè),且預(yù)測(cè)精度能夠滿足工程需求。將eN方法求解模塊與現(xiàn)有的RANS框架下的CFD求解器進(jìn)行耦合使用,拓展了線性穩(wěn)定性分析的應(yīng)用,使得其能夠應(yīng)用到常規(guī)氣動(dòng)外形上,不再局限于簡單外形。

      另一方面,眾多研究者致力于建立N值包絡(luò)與層流基本流動(dòng)的特征參數(shù)之間的顯式函數(shù)映射關(guān)系。1987年,Drela和Giles[27]根據(jù)線性穩(wěn)定性分析方法,對(duì)不同壓力梯度下二維邊界層的層流相似性解,進(jìn)行了擾動(dòng)增長因子N的計(jì)算,獲得了T-S波的擾動(dòng)增長因子與邊界層形狀因子和動(dòng)量厚度雷諾數(shù)的顯式關(guān)系式,首次對(duì)線性穩(wěn)定性理論進(jìn)行了簡化,使得eN方法在轉(zhuǎn)捩預(yù)測(cè)中的應(yīng)用更加簡單方便。

      隨著近年來基于當(dāng)?shù)刈兞康霓D(zhuǎn)捩預(yù)測(cè)模式的發(fā)展,將eN方法變成CFD轉(zhuǎn)捩模式的做法越來越流行。從2013年開始,基于Drela和Giles[27]的研究成果,Coder和Maughmer[28-30]首次提出了流向擾動(dòng)增長因子NTS的輸運(yùn)模式。他們通過對(duì)不同壓力梯度下的層流邊界層相似性解進(jìn)行大量的變參數(shù)穩(wěn)定性分析,找出了擾動(dòng)增長因子包絡(luò)線與當(dāng)?shù)氐男螤钜蜃雍蛣?dòng)量厚度雷諾數(shù)以及邊界層特征尺度之間的函數(shù)關(guān)系式。對(duì)于出現(xiàn)在關(guān)系式中的非當(dāng)?shù)刈兞浚眠m當(dāng)?shù)闹甘疽蜃?,?gòu)造出新的關(guān)系式,將這些非當(dāng)?shù)刈兞浚饌€(gè)進(jìn)行當(dāng)?shù)鼗?jì)算,最終與湍流模式耦合使用,形成轉(zhuǎn)捩-湍流預(yù)測(cè)模式。該模式的應(yīng)用非常簡便,可以與CFD求解器在計(jì)算中一起使用,一次性將層流區(qū)和湍流區(qū)的分布以及氣動(dòng)力等流場(chǎng)數(shù)據(jù)求解出來。之后,Coder[31-32]又提出了滿足伽利略不變性的形狀因子參數(shù),使得該模式在非慣性參考系中的應(yīng)用有了理論支持。

      然而,Coder[31-32]提出的預(yù)測(cè)模式只能對(duì)流向T-S波轉(zhuǎn)捩和層流分離泡轉(zhuǎn)捩進(jìn)行預(yù)測(cè),不能對(duì)橫流轉(zhuǎn)捩進(jìn)行預(yù)測(cè)。為了完善這一功能,徐家寬等[33-34]提出了橫流駐波擾動(dòng)增長因子輸運(yùn)模式,以橫流強(qiáng)度因子和螺旋度雷諾數(shù)為根基發(fā)展了完全基于當(dāng)?shù)刈兞康腘CF輸運(yùn)模式,在較多低速算例上進(jìn)行了驗(yàn)證,證實(shí)了其合理性。最近,該模式又進(jìn)一步被添加了壁面溫度效應(yīng)[35-36],使其可以對(duì)跨聲速的部分變壁面溫度算例進(jìn)行預(yù)測(cè)且預(yù)測(cè)精度較高。至此,基于線性穩(wěn)定性理論構(gòu)造的NTS-NCF低速和亞聲速轉(zhuǎn)捩模式的框架基本成型。2020年,徐家寬等[37]對(duì)超聲速邊界層的Oblique T-S波[37]和可壓縮橫流駐波的擾動(dòng)增長因子[38]進(jìn)行了輸運(yùn)模式構(gòu)建,但結(jié)果當(dāng)前僅在平板和無限展長后掠翼等簡單構(gòu)型上精度較高,在真實(shí)復(fù)雜三維外形表面的預(yù)測(cè)還需進(jìn)一步改進(jìn)和完善。

      下面將對(duì)NTS-NCF模式的發(fā)展歷程、特點(diǎn)和存在的不足等方面進(jìn)行描述。

      1 流向擾動(dòng)增長因子輸運(yùn)模式

      1.1 發(fā)展歷程和特點(diǎn)分析

      NTS模式原始版本詳見Coder和Maughmer在2014年發(fā)表的論文[29],其輸運(yùn)方程形式為

      (1)

      式中:ρ為密度;ui為xi方向的速度;PTS為源項(xiàng);σTS為擴(kuò)散項(xiàng)系數(shù);μ和μt分別為層流和湍流黏性系數(shù)。源項(xiàng)PTS的構(gòu)造主要基于Drela和Giles在1987年發(fā)表的成果[27],包含T-S波開始失穩(wěn)點(diǎn)的判據(jù)、隨著動(dòng)量厚度雷諾數(shù)增加而增大的斜率以及動(dòng)量厚度雷諾數(shù)沿流向的變化導(dǎo)數(shù)等。沿流向的積分效應(yīng)可以通過輸運(yùn)方程的性質(zhì)去實(shí)現(xiàn),但是所有的計(jì)算公式都是形狀因子H12的函數(shù),而形狀因子是一個(gè)非當(dāng)?shù)氐姆e分變量。構(gòu)造RANS方程物理模式最關(guān)鍵的點(diǎn)在于所有變量可當(dāng)?shù)鼗蠼?,因此,如何讓形狀因子H12當(dāng)?shù)鼗蠼獬蔀榻5碾y點(diǎn)。

      Coder和Maughmer 2014年[29]提出的當(dāng)?shù)鼗瘔毫μ荻葏?shù)為

      (2)

      式中:S為剪切應(yīng)變率的模;d為距離壁面的最小距離;邊界層邊緣速度Ue由可壓縮伯努利方程獲得

      (3)

      其中:U∞為自由來流速度;p為當(dāng)?shù)貕毫?;p∞為自由來流壓力;ρ∞為自由來流的密度;γH=1.4。

      圖1給出了不可壓邊界層層流相似性解中當(dāng)?shù)貕毫μ荻葏?shù)HL的分布[39],圖中η為無量綱的壁面法向坐標(biāo),βH為Hartree壓力梯度因子??梢钥闯觯琀L在邊界層中部的最大值與形狀因子呈單調(diào)變化關(guān)系,因此,可以用其表征形狀因子的大小,計(jì)算公式為

      圖1 當(dāng)?shù)鼗瘔毫μ荻葏?shù)HL在Falkner-Skan相似速度型中的分布[39]

      (4)

      該模式的整體流程框架如圖2所示,得到T-S波的增長因子之后,需要人為給定轉(zhuǎn)捩臨界N值才能判斷轉(zhuǎn)捩,一般在航空航天小擾動(dòng)飛行環(huán)境里臨界N值一般使用Mack[40]提出的與來流湍流度之間的關(guān)系式:

      NTS,crit=-8.43-2.4ln(Tu/100)

      (5)

      γeff=max[γ,exp(2(NTS-NTS,crit))]

      (6)

      (7)

      (8)

      式中:θ為動(dòng)量損失厚度;U為平板流向速度;V為平板壁面法向速度;n為壁面法向單位矢量;u為速度矢量。然后用任意網(wǎng)格點(diǎn)距離壁面的最小距離d替換動(dòng)量損失厚度θ,得到

      (9)

      Coder[31-32]在此基礎(chǔ)上發(fā)展了類似定義的當(dāng)?shù)貕毫μ荻戎甘疽蜃雍蛯?duì)應(yīng)的形狀因子計(jì)算函數(shù):

      (10)

      H12=min[max(0.26HL+2.4,2.2),20.0]

      (11)

      該版本對(duì)應(yīng)的其他函數(shù)公式詳見圖2,圖中變量的定義詳見文獻(xiàn)[29,32],其中HL和H12在S&K平板算例中的分布如圖3所示[32],圖中z為壁面法向坐標(biāo),Rex為距離平板前緣距離為x位置的流向雷諾數(shù)。

      圖2 T-S波擾動(dòng)增長因子輸運(yùn)方程和計(jì)算框架

      圖3 滿足伽利略不變性的HL和H12在S&K平板算例中的分布[32]

      1.2 存在的問題和需要改進(jìn)的地方

      最新的NTS模式還存在3個(gè)方面的缺陷:

      1) 形狀因子H12的計(jì)算結(jié)果不夠物理。仔細(xì)分析可以發(fā)現(xiàn),該版本計(jì)算出的形狀因子H12在零壓力梯度平板的值并不是理論上的2.59,而在3附近。因此本文作者對(duì)此進(jìn)行了修正,其表達(dá)式為

      H12=min[max(0.093 5HL+2.382,2.2),8.0]

      (12)

      這樣才能使得形狀因子H12的計(jì)算結(jié)果符合基本的物理規(guī)律:例如零壓力梯度平板的值為2.59,臨界分離點(diǎn)的形狀因子為4.0等。

      2) 跨聲速高雷諾數(shù)工況計(jì)算所得N值過大。在NASA進(jìn)行的層流標(biāo)模計(jì)算對(duì)比驗(yàn)證過程中發(fā)現(xiàn)Langtry的經(jīng)驗(yàn)關(guān)系模式[45]和NTS模型[46-48]在基于平均氣動(dòng)弦長的雷諾數(shù)(ReMAC)超過 1 500萬時(shí),就已經(jīng)失效。NTS轉(zhuǎn)捩模式計(jì)算所得的N值遠(yuǎn)遠(yuǎn)大于標(biāo)準(zhǔn)LST分析的值,導(dǎo)致轉(zhuǎn)捩非常靠近前緣,如圖4所示[48],其中AOA表示迎角。因此,當(dāng)前的NTS模式需要做高雷諾數(shù)增長函數(shù)修正。而且Drela和Giles[27]用的是運(yùn)動(dòng)學(xué)形狀因子Hk,其對(duì)跨聲速邊界層的可壓縮性并不敏感,如果只參照不可壓的函數(shù)關(guān)系進(jìn)行形狀因子的求解,則會(huì)造成跨聲速邊界層內(nèi)形狀因子的過大估計(jì),導(dǎo)致順壓梯度的H12值過大,增長斜率向逆壓梯度的公式傾斜。這也是造成N值預(yù)測(cè)偏大的重要因素,所以還需要進(jìn)行可壓縮性修正。

      圖4 NASA跨聲速高雷諾數(shù)層流標(biāo)模試驗(yàn)和計(jì)算結(jié)果對(duì)比(ReMAC=1.5×107,AOA=1.44°, 1 in=2.54 cm)[48]

      3) 不具備橫流轉(zhuǎn)捩的預(yù)測(cè)能力。由于Coder只對(duì)T-S波進(jìn)行了建模,因此,該模式無法預(yù)測(cè)橫流轉(zhuǎn)捩。其在文章中還展示了該模式對(duì)于橫流現(xiàn)象主導(dǎo)的算例的預(yù)測(cè)結(jié)果,如圖5所示[32],其中Cf為摩擦力系數(shù),L為橢球長軸的長度。

      圖5 NTS模式在變后掠機(jī)翼和6∶1標(biāo)準(zhǔn)橢球大迎角算例中的預(yù)測(cè)結(jié)果與試驗(yàn)數(shù)據(jù)的對(duì)比[32]

      2 橫流擾動(dòng)增長因子輸運(yùn)模式

      2.1 發(fā)展歷程和特點(diǎn)分析

      由于第1節(jié)的流向擾動(dòng)增長因子輸運(yùn)方程不具備橫流轉(zhuǎn)捩的預(yù)測(cè)能力,而在亞跨聲速的小擾動(dòng)來流環(huán)境中,橫流駐波主導(dǎo)橫流失穩(wěn),因此需要對(duì)橫流駐波進(jìn)行NCF因子的輸運(yùn)模式構(gòu)造。徐家寬等[33]2016年初步構(gòu)造了適用于無限展長后掠翼的NCF因子輸運(yùn)方程,于2019年基于大量線性穩(wěn)定性分析結(jié)果,構(gòu)造了適用于復(fù)雜氣動(dòng)外形的NCF擾動(dòng)增長因子輸運(yùn)模式[34-35],方程形式為

      (13)

      式中:σCF表示擴(kuò)散項(xiàng)系數(shù)。源項(xiàng)PCF的建模也包含了橫流駐波開始失穩(wěn)點(diǎn)的判定、擾動(dòng)開始增長以后的增長斜率以及增長率修正函數(shù),具體的公式如圖6所示。

      圖6 橫流駐波擾動(dòng)增長因子輸運(yùn)方程和計(jì)算框架

      得到橫流駐波的N因子之后,與1.1節(jié)的流向擾動(dòng)增長因子就形成了NTS-NCF轉(zhuǎn)捩預(yù)測(cè)框架。有效間歇因子的表達(dá)式也變?yōu)?/p>

      γeff=max{γ,exp[2max(NTS-NTS,crit,

      NCF-NCF,crit)]}

      (14)

      哪種失穩(wěn)模態(tài)比較強(qiáng)就會(huì)率先觸發(fā)轉(zhuǎn)捩,從而實(shí)現(xiàn)多種失穩(wěn)模態(tài)同時(shí)預(yù)測(cè)的能力。

      如圖7[34]所示,橫流速度型中橫流雷諾數(shù)Recf和橫流形狀因子Hcf的定義為

      Recf=|wmax|δCF/νe

      (15)

      Hcf=ymax/δCF

      (16)

      式中:νe為邊界層邊緣的運(yùn)動(dòng)黏度。

      需要指出的是δCF的高度為最大橫流速度的10%位置(靠近邊界層邊緣的一側(cè))。其中橫流速度的當(dāng)?shù)鼗褂脵M流強(qiáng)度因子,橫流雷諾數(shù)的當(dāng)?shù)鼗褂寐菪壤字Z數(shù)進(jìn)行。

      本文螺旋度因子首先由Müller和Herbst[49]提出,后經(jīng)Langtry[45]改造成橫流強(qiáng)度因子,經(jīng)德國宇航中心(DLR)研究人員[50]改造成螺旋度雷諾數(shù)。該因子不滿足伽利略不變性,但是可以表征不同當(dāng)?shù)睾舐咏窍碌臋M流分量強(qiáng)度。

      本文作者團(tuán)隊(duì)對(duì)Langtry提出的橫流強(qiáng)度因子和DLR研究人員提出的螺旋度雷諾數(shù)做一點(diǎn)微調(diào)[35],使得最大橫流速度分量與橫流強(qiáng)度因子、橫流雷諾數(shù)與螺旋度雷諾數(shù)之間的函數(shù)映射可以規(guī)避掉后掠角的影響。微調(diào)后的螺旋度雷諾數(shù)和橫流強(qiáng)度因子定義為

      (17)

      (18)

      新定義的螺旋度雷諾數(shù)與橫流雷諾數(shù)之間的函數(shù)關(guān)系和新定義的橫流強(qiáng)度因子與最大橫流速度分量之間的函數(shù)關(guān)系如圖8[34]和圖9[34]所示,圖中βH為Hartree壓力梯度,F(xiàn)SC表示Falkner-Skan-Cooke三維層流邊界層相似性解,詳細(xì)計(jì)算公式見文獻(xiàn)[34]和圖6。其中可壓縮修正項(xiàng)包含了邊界層邊緣馬赫數(shù)Mae=Ue/ae,ae為邊界層邊緣的聲速。該參數(shù)的計(jì)算繼續(xù)使用可壓縮伯努利方程以及總溫守恒關(guān)系[51-52]

      圖8 新定義的螺旋度雷諾數(shù)與橫流雷諾數(shù)的函數(shù)關(guān)系[34]

      圖9 新定義的橫流強(qiáng)度因子與最大橫流速度分量的函數(shù)關(guān)系[34]

      (19)

      (20)

      最后,該模式還剩下Hartree壓力梯度βH的當(dāng)?shù)鼗?jì)算。其可以使用更為精細(xì)的不同后掠角下的三維邊界層壓力梯度因子計(jì)算公式,但是二維邊界層壓力梯度因子的計(jì)算精度已經(jīng)能夠滿足當(dāng)前的計(jì)算需求,其與Thwaites壓力梯度因子的關(guān)系式為

      (21)

      詳細(xì)計(jì)算公式見文獻(xiàn)[34-35]和圖6。

      圖10[34]給出了NLF(2)-0415無限展長后掠翼在弦長雷諾數(shù)為3.81×106時(shí)預(yù)測(cè)所得NCF云圖和與標(biāo)準(zhǔn)LST結(jié)果的對(duì)比,圖中βr表示橫流駐波的展向波數(shù)。在每個(gè)站位取最大N值組成的N值包絡(luò)線與標(biāo)準(zhǔn)LST分析結(jié)果吻合良好。該模式在鐮刀翼上的結(jié)果如圖11[34]所示,也在6∶1 橢球上進(jìn)行了詳細(xì)測(cè)試。針對(duì)相同雷諾數(shù),不同迎角狀態(tài)下預(yù)測(cè)所得摩擦力系數(shù)分布如圖12[35]所示,其中雷諾數(shù)ReL是基于橢球長軸長度L的特征雷諾數(shù)。而針對(duì)相同迎角,不同雷諾數(shù)下的轉(zhuǎn)捩預(yù)測(cè)結(jié)果如圖13[35]所示,其中Φ為周向角度。由圖可知,采用NTS-NCF耦合框架,可以很好地捕捉到橫流轉(zhuǎn)捩和T-S波轉(zhuǎn)捩。

      圖10 NLF(2)-0415算例中NCF模式預(yù)測(cè)結(jié)果與標(biāo)準(zhǔn)LST結(jié)果的對(duì)比[34]

      圖11 NCF模式在鐮刀翼上表面的預(yù)測(cè)結(jié)果[34]

      圖12 不同迎角下NCF模式在6∶1橢球上的預(yù)測(cè)結(jié)果(Ma=0.136, ReL=6.5×106) [35]

      圖13 NCF模式在6∶1橢球上不同雷諾數(shù)下的轉(zhuǎn)捩預(yù)測(cè)結(jié)果與試驗(yàn)數(shù)據(jù)和標(biāo)準(zhǔn)LST結(jié)果的對(duì)比(Ma=0.136, AOA=29.5°) [35]

      2020年作者團(tuán)隊(duì)對(duì)當(dāng)前NCF模式進(jìn)行了可壓縮修正[35],該模式在DLR-F4機(jī)翼上的預(yù)測(cè)結(jié)果如圖14所示,圖中CL表示升力系數(shù)。此外,該模式還應(yīng)用于ARA層流機(jī)翼[36],NCF模式預(yù)測(cè)所得的N值包絡(luò)與標(biāo)準(zhǔn)LST計(jì)算結(jié)果吻合良好,如圖15[36]所示。在臨界值給7.0的情況下與風(fēng)洞試驗(yàn)結(jié)果吻合,如圖16[36]所示。

      圖14 NCF模式在DLR-F4機(jī)翼上表面的預(yù)測(cè)結(jié)果與風(fēng)洞試驗(yàn)結(jié)果的對(duì)比[35]

      圖15 NCF模式在ARA構(gòu)型上的預(yù)測(cè)結(jié)果與標(biāo)準(zhǔn)LST結(jié)果的對(duì)比(AOA=-3.34°)[36]

      圖16 NCF模式與NTS模式預(yù)測(cè)結(jié)果的對(duì)比[36]

      2.2 存在的問題和需要改進(jìn)的地方

      NCF模式在亞跨聲速邊界層的表現(xiàn)還是比較出色,其進(jìn)行了可壓縮修正。但是當(dāng)前NCF模式還存在以下問題需要進(jìn)一步研究和改進(jìn):

      1) 壓力梯度參數(shù)、螺旋度雷諾數(shù)和橫流強(qiáng)度因子3個(gè)關(guān)鍵變量均不滿足伽利略不變性。壓力梯度參數(shù)可以利用Coder[31-32]提出的壓力梯度因子進(jìn)行替換,解決伽利略不變性的問題。橫流強(qiáng)度因子和螺旋度雷諾數(shù)則需要繼續(xù)尋找其他可替代的變量。突破口為“避免使用絕對(duì)速度,嘗試用流動(dòng)變量的法向梯度來替換”。

      2) 橫流駐波主導(dǎo)的轉(zhuǎn)捩臨界N值判定準(zhǔn)則還需要進(jìn)一步構(gòu)造,盡管Crouch和Ng[53]提出了初步的經(jīng)驗(yàn)關(guān)系,但還遠(yuǎn)遠(yuǎn)不夠普適。因此,橫流駐波的臨界N值選取準(zhǔn)則還需要進(jìn)一步研究和探索。

      3 總結(jié)和展望

      文章對(duì)擾動(dòng)增長因子NTS-NCF模式的發(fā)展歷程進(jìn)行了詳細(xì)的梳理和分析,總結(jié)了其克服當(dāng)?shù)鼗蠼夂唾だ圆蛔冃缘难芯克悸?,并指出了?dāng)前該預(yù)測(cè)框架存在的不足:

      1) 當(dāng)前流向NTS因子模式需要進(jìn)行可壓縮性修正和高雷諾數(shù)修正,從而避免在跨聲速邊界層出現(xiàn)過大估計(jì)的現(xiàn)象。

      2)NCF因子模式也需要進(jìn)行高雷諾數(shù)驗(yàn)證并且需要對(duì)表征橫流強(qiáng)度因子的當(dāng)?shù)鼗甘疽蜃舆M(jìn)行改進(jìn),使其滿足伽利略不變性。

      3) 超/高聲速邊界層中出現(xiàn)的三維T-S斜波失穩(wěn)、第二模態(tài)失穩(wěn)和可壓縮橫流模態(tài)失穩(wěn)都需要進(jìn)一步研究或改進(jìn)。

      4) 判定轉(zhuǎn)捩的臨界N值如何選取,也是一個(gè)開放性的問題,需要大量的試驗(yàn)數(shù)據(jù)和理論分析來構(gòu)建一個(gè)較為普適的選取準(zhǔn)則。

      猜你喜歡
      橫流壓力梯度雷諾數(shù)
      橫流熱源塔換熱性能研究
      煤氣與熱力(2021年3期)2021-06-09 06:16:20
      基于Transition SST模型的高雷諾數(shù)圓柱繞流數(shù)值研究
      基于橫流風(fēng)扇技術(shù)的直升機(jī)反扭驗(yàn)證
      壓力梯度在油田開發(fā)中的應(yīng)用探討
      失穩(wěn)初期的低雷諾數(shù)圓柱繞流POD-Galerkin 建模方法研究
      基于轉(zhuǎn)捩模型的低雷諾數(shù)翼型優(yōu)化設(shè)計(jì)研究
      民機(jī)高速風(fēng)洞試驗(yàn)的阻力雷諾數(shù)效應(yīng)修正
      脊下橫流對(duì)PEMFC性能影響的數(shù)值分析
      疊加原理不能求解含啟動(dòng)壓力梯度滲流方程
      致密砂巖啟動(dòng)壓力梯度數(shù)值的影響因素
      斷塊油氣田(2014年5期)2014-03-11 15:33:45
      凤山市| 上饶市| 运城市| 于都县| 锡林郭勒盟| 兴海县| 贵德县| 宜兰市| 成武县| 毕节市| 齐齐哈尔市| 乐业县| 江安县| 广昌县| 息烽县| 吴堡县| 奉节县| 上饶县| 宜君县| 五家渠市| 堆龙德庆县| 裕民县| 类乌齐县| 新野县| 南昌市| 新沂市| 牙克石市| 顺平县| 阿克陶县| 义马市| 大城县| 中江县| 元江| 万宁市| 呼玛县| 长武县| 宜黄县| 桂林市| 改则县| 库伦旗| 丰都县|