王方劍 王宏偉 李曉輝 董 磊 黃 湛 陳 蘭

(中國(guó)航天空氣動(dòng)力技術(shù)研究院,北京 100074)

引言

空空導(dǎo)彈作為現(xiàn)代空戰(zhàn)的主要攻擊手段,要求比目標(biāo)飛機(jī)更高的機(jī)動(dòng)性和敏捷性.新型空空導(dǎo)彈在面對(duì)新一代飛機(jī)時(shí)必須具備全方位攻擊能力,尤其對(duì)來(lái)自后方目標(biāo)的威脅,則需要更高轉(zhuǎn)彎率和更大機(jī)動(dòng)包絡(luò)線的航向反轉(zhuǎn)機(jī)動(dòng)等先進(jìn)高效機(jī)動(dòng)方法.為了保證高效機(jī)動(dòng)的順利完成,要求導(dǎo)彈在超大攻角(α=0°～ 180°)范圍內(nèi)具有飛行和機(jī)動(dòng)控制能力[1].以往對(duì)超大攻角流動(dòng)的觀測(cè)和研究大多集中在α=40°～ 60°范圍內(nèi),最大角度不超過(guò)90°[2-14].在超大攻角條件下,復(fù)雜的非定常分離流問(wèn)題廣泛存在,包括非定常渦的產(chǎn)生、脫落、干擾和非定常非對(duì)稱(chēng)等氣動(dòng)現(xiàn)象[2-10].因此,提高對(duì)細(xì)長(zhǎng)體復(fù)雜非定常氣動(dòng)特性的認(rèn)識(shí)尤為重要.

早在20 世紀(jì)50 年代[3],對(duì)細(xì)長(zhǎng)體(錐柱體)分離流動(dòng)特性的探索就已經(jīng)開(kāi)始.在小攻角處,細(xì)長(zhǎng)體背風(fēng)側(cè)形成對(duì)稱(chēng)的軸向集中渦[4].隨著攻角的增大(達(dá)到30°左右),流型由對(duì)稱(chēng)型轉(zhuǎn)變?yōu)榉菍?duì)稱(chēng)型.隨著旋渦沿軸向向后發(fā)展,旋渦一個(gè)接一個(gè)離開(kāi)物面,形成旋渦脫落流動(dòng)結(jié)構(gòu),并以一定角度向下游流動(dòng).此時(shí)如果觀察其橫截面,會(huì)發(fā)現(xiàn)形成類(lèi)卡門(mén)渦街的流動(dòng)形式[5].

在無(wú)側(cè)滑來(lái)流條件下的側(cè)向力是非對(duì)稱(chēng)渦的重要標(biāo)志[9-10].Lamont 和Hunt[11]的研究表明,在大部分攻角范圍內(nèi)(0°～ 90°),沿旋成體軸向的側(cè)向力呈振蕩分布.Dexter[12]發(fā)現(xiàn)側(cè)向力強(qiáng)烈依賴(lài)于旋成體滾轉(zhuǎn)角.在一個(gè)可獨(dú)立滾轉(zhuǎn)的頭尖部上進(jìn)行的風(fēng)洞試驗(yàn)表明,滾轉(zhuǎn)角所造成的對(duì)側(cè)向力的影響主要頭尖部帶來(lái).Luo 等[13]在一個(gè)橢圓形截面頭尖部的旋成體上進(jìn)行了壓力分布測(cè)量,結(jié)果表明,橢圓形截面頭尖部側(cè)向力隨滾轉(zhuǎn)角的變化規(guī)律比圓形更容易預(yù)測(cè),且兩種頭尖部的最大側(cè)向力基本相同.通過(guò)在頭尖部附近安裝一根可伸縮的鋼絲[14]或沙礫[15],可以通過(guò)改變頭尖部的周向角來(lái)有效地控制側(cè)向力的方向.大多數(shù)文獻(xiàn)采用的試驗(yàn)?zāi)Ｐ褪穷^部尖銳的旋成體外形,然而頭部形狀也會(huì)對(duì)流動(dòng)產(chǎn)生一定影響.Hsieh 等[16-17]對(duì)半球形頭部旋成體在攻角從0°～50°的工況下的分離流進(jìn)行了試驗(yàn)和數(shù)值研究.令人驚訝的是,分離流包含了許多復(fù)雜的流動(dòng)結(jié)構(gòu)特征,這是三維流所獨(dú)有的.結(jié)果表明,半球形頭部旋成體同時(shí)存在兩個(gè)分離區(qū)(頭部分離泡和橫流分離區(qū)).

20 世紀(jì)80 年代以前,由于缺乏強(qiáng)大的計(jì)算軟件和先進(jìn)的風(fēng)洞試驗(yàn)方法,學(xué)者們關(guān)注點(diǎn)主要集中在附著流和穩(wěn)定分離區(qū)上[18].隨著技術(shù)的發(fā)展,一些非定常流動(dòng)特性被揭示出來(lái).Zeiger 等[19]對(duì)細(xì)長(zhǎng)旋成體的流動(dòng)特性進(jìn)行流動(dòng)顯示水洞試驗(yàn),認(rèn)為在大攻角時(shí),細(xì)長(zhǎng)體旋成體上的流動(dòng)可主要分為3 種流動(dòng)形態(tài),第1 種流動(dòng)形態(tài)是頭尖部附近的產(chǎn)生的較為集中的渦旋,隨著流動(dòng)沿軸向向后發(fā)展,流動(dòng)結(jié)構(gòu)逐漸體現(xiàn)為旋渦脫落流動(dòng)形態(tài)(第2 種流動(dòng)形態(tài)),渦脫落方向與旋成體呈一定角度,當(dāng)流動(dòng)發(fā)展至接近底部時(shí)(橫截面位置X/D 從10 到32),脫落的旋渦渦軸與旋成體平行,形成第3 流動(dòng)形態(tài).管小榮和徐誠(chéng)[20]、楊云軍和周偉江[21]、張贏和劉超峰[22]、劉仙名和符松[23],采用數(shù)值模擬方法細(xì)致的研究了非對(duì)稱(chēng)旋渦的演化與發(fā)展,研究表明不同的Re、湍流模型均會(huì)對(duì)旋渦非對(duì)稱(chēng)流動(dòng)的模擬產(chǎn)生一定影響.Degani 和Zilliac[24]的研究表明,其非定常流動(dòng)主要包括類(lèi)卡門(mén)渦街的低頻流動(dòng)、剪切層失穩(wěn)引起的高頻流動(dòng),以及介于以上兩種頻率之間的旋渦干擾流動(dòng).Ayoub 和Karamcheti[25]同時(shí)測(cè)量了細(xì)長(zhǎng)旋成體表面壓力和尾流速度波動(dòng).結(jié)果表明,在頭部區(qū)域存在一定的旋渦脫落現(xiàn)象.在離旋成體頭部不同距離處,渦脫落產(chǎn)生了一個(gè)接近10.5 Hz 的頻率.

除此之外,馬赫數(shù)對(duì)旋成體分離流動(dòng)細(xì)節(jié)也有一定影響,Keneer[26]對(duì)旋成體前體在不同馬赫數(shù)及攻角細(xì)長(zhǎng)體流動(dòng)進(jìn)行研究,在風(fēng)洞試驗(yàn)紋影照片中可以觀察到,第一脫落渦在高亞聲速下(Ma=0.6)的軌跡比低速(Ma=0.25)時(shí)更接近表面,側(cè)力系數(shù)從Ma=0.25 時(shí)的1.9 降至Ma=0.6 時(shí)的0.8,同時(shí)也可以看到,在Ma=0.6 時(shí),攻角55°時(shí)后半部分出現(xiàn)不對(duì)稱(chēng),旋渦較為分散,而Ma=0.25 在此攻角下攻角還較為集中,這體現(xiàn)出旋成體高亞音速流動(dòng)與低速旋渦流動(dòng)相比,分離流動(dòng)的旋渦位置、旋渦強(qiáng)度、渦破裂位置、分離位置均有一定的不同.

在傳統(tǒng)的復(fù)雜分離流動(dòng)數(shù)值模擬工作中,基于雷諾平均Navier-Stokes (RANS)方程的流場(chǎng)預(yù)測(cè)方法對(duì)附著流以及小分離流動(dòng)有較好的模擬能力,計(jì)算成本較低,是工程運(yùn)用中最廣泛的一種方法,但由于RANS 方法平均了小尺度的旋渦脈動(dòng),難以精細(xì)的模擬較大范圍的分離與復(fù)雜的旋渦流動(dòng).大渦模擬方法(LES)能夠較好的模擬大范圍分離旋渦流動(dòng)現(xiàn)象,但對(duì)邊界層的網(wǎng)格密度有非常高的要求,帶來(lái)很大的計(jì)算成本.所以為了結(jié)合RANS 方法與LES方法的優(yōu)勢(shì),Spalart[27]提出了基于S-A 湍流模型的DES 方法,引入特征長(zhǎng)度參數(shù),根據(jù)特征長(zhǎng)度選擇不同的流動(dòng)模擬方法,在近壁面區(qū)域,采用RANS 方法,在遠(yuǎn)離物面的區(qū)域,選擇LES 方法進(jìn)行模擬,能夠兼顧數(shù)值計(jì)算精度與效率.深入研究發(fā)現(xiàn),DES 方法會(huì)發(fā)生模型應(yīng)力耗散現(xiàn)象(modeled stress depletion,MSD),導(dǎo)致在流動(dòng)中出現(xiàn)網(wǎng)格誘導(dǎo)分離問(wèn)題(grid induced separa-tion,GIS).2006 年,Spalart 等[28-37]提出了延遲脫體渦模擬方法(DDES)有效解決MSD與GIS 問(wèn)題,該方法被廣泛應(yīng)用于圓柱擾流、空腔流動(dòng)和起落架分離流動(dòng)能計(jì)算研究中.

從細(xì)長(zhǎng)體大攻角的研究情況來(lái)看,之前對(duì)細(xì)長(zhǎng)體大攻角流動(dòng)研究速域主要集中在低速,并且研究攻角范圍主要在α=0°～ 90°之間,本文針對(duì)細(xì)長(zhǎng)旋成體在高速(Ma=0.6)超大攻角下(α=0°～ 180°)存在的旋渦非定常性、旋渦脫落和旋渦非對(duì)稱(chēng)等復(fù)雜流動(dòng)現(xiàn)象,采用基于結(jié)構(gòu)網(wǎng)格的DDES 數(shù)值方法,開(kāi)展細(xì)長(zhǎng)旋成體復(fù)雜氣動(dòng)非定常特性研究,為新型機(jī)動(dòng)導(dǎo)彈超大攻角氣動(dòng)設(shè)計(jì)提供理論支撐.

1 研究方法

本文采用數(shù)值模擬與油流顯示風(fēng)洞試驗(yàn)相結(jié)合的方法來(lái)研究細(xì)長(zhǎng)體超大攻角流動(dòng)特性問(wèn)題.其中數(shù)值模擬采用Roe 方法求解Navier-Stokes 方程,細(xì)致分析細(xì)長(zhǎng)體超大攻角在不同攻角下的非定常流動(dòng)特性及壓力脈動(dòng)特性,油流顯示風(fēng)洞試驗(yàn)用于獲得細(xì)長(zhǎng)體模型表面的流動(dòng)拓?fù)?通過(guò)對(duì)比試驗(yàn)結(jié)果與數(shù)值結(jié)果,驗(yàn)證數(shù)值模擬結(jié)果的準(zhǔn)確性.

1.1 數(shù)值方法與算例驗(yàn)證

控制方程為三維可壓縮Navier-Stokes 方程,廣義坐標(biāo)表示為

空間離散采用基于MUSCL 方法插值方法的FDS-Roe 格式,黏性項(xiàng)采用二階中心差分離散,對(duì)流和壓力項(xiàng)使用3 階迎風(fēng)格式.時(shí)間推進(jìn)格式為隱式LUSGS 方法,保證了較高的時(shí)間計(jì)算精度.非定常計(jì)算采用雙時(shí)間步長(zhǎng)方法,同時(shí)采用多重網(wǎng)格算法加快子迭代的收斂速度.

在用湍流模型進(jìn)行流動(dòng)模擬時(shí),控制方程中應(yīng)加入雷諾應(yīng)力相關(guān)項(xiàng).本文采用的湍流模型是基于Menter 的兩方程剪切應(yīng)力輸運(yùn)湍流模型(SST),表達(dá)式為

DES 方法的擴(kuò)展形式是基于SST 湍流模型,其中湍流尺度為

DES 長(zhǎng)度尺度dDES表示為

式中,Δ 為當(dāng)前網(wǎng)格單元與周邊單元的最大距離,CDES是SST 湍流模型中附加的經(jīng)驗(yàn)?zāi)Ｐ统?shù),CDES表示為

為了解決MSD 與GIS 問(wèn)題問(wèn)題,Spalart[27]提出了延遲函數(shù)fd,將SST 湍流模型的湍流能量輸運(yùn)方程的耗散項(xiàng)重新定義為與網(wǎng)格尺度和湍流黏度有關(guān)的DES 長(zhǎng)度尺度,結(jié)合延遲函數(shù)fd的長(zhǎng)度尺度dDES被重新定義為

因此,dDES不僅與網(wǎng)格尺度有關(guān),而且與湍流黏度有關(guān),可以防止湍流模擬過(guò)早地切換到LES 模式.

算例驗(yàn)證采用高雷諾數(shù)圓柱擾流算例進(jìn)行計(jì)算,圓柱直徑為D=457 mm,圓柱長(zhǎng)度為5.66D,圓柱兩端采用對(duì)稱(chēng)邊界條件,其他邊界采用自由來(lái)流條件.圖1 為圓柱擾流計(jì)算域,為保證滿足自由來(lái)流邊界條件,在流向計(jì)算域?yàn)?0D,圖2 是圓柱橫截面網(wǎng)格.計(jì)算參數(shù)為:Ma=0.25,Re=3 × 106,物理時(shí)間步長(zhǎng)為Δt=4 × 10-4s,一共計(jì)算5000 個(gè)物理時(shí)間步,計(jì)算物理時(shí)間為2 s.

圖1 圓柱擾流計(jì)算域Fig.1 Domain of calculation

圖2 圓柱橫截面網(wǎng)格Fig.2 Zoomed-in-view of mesh around the cylinder

圖3 是圓柱擾流壓力分布試驗(yàn)值[38-40]與計(jì)算值壓力分布對(duì)比,其中縱軸是壓力系數(shù),橫軸是方位角θ,計(jì)算值是截取了圓柱中間截面(y=2.83D)的壓力分布時(shí)均值,通過(guò)計(jì)算獲得的圓柱升力斯特勞哈爾數(shù)St=0.32,并將取一個(gè)周期內(nèi)的60 個(gè)瞬時(shí)值進(jìn)行平均,得到時(shí)均值.圖中顯示,計(jì)算值與試驗(yàn)值吻合較好,說(shuō)明本文所采用的DDES 方法能夠較好的捕捉圓柱尾跡的旋渦結(jié)構(gòu),見(jiàn)圖4.

圖3 壓力分布試驗(yàn)值[38-40]與計(jì)算值壓力分布對(duì)比Fig.3 Comparisons of pressure coefficients distribution[38-40]

圖4 圓柱擾流尾跡流動(dòng)結(jié)構(gòu)Fig.4 Q-criterion iso-surface of wake

1.2 油流流動(dòng)顯示技術(shù)

為了便于油流的光學(xué)測(cè)量技術(shù)的應(yīng)用,試驗(yàn)選定在中國(guó)航天空氣動(dòng)力技術(shù)研究院FD-12 風(fēng)洞常規(guī)亞跨聲速試驗(yàn)段作為流場(chǎng)顯示試驗(yàn)的主要試驗(yàn)環(huán)境,該試驗(yàn)段長(zhǎng)度為3.8 m,試驗(yàn)段橫截面尺寸為1.2 m ×1.2 m,兩側(cè)沿來(lái)流方向?qū)ΨQ(chēng)布置4 個(gè)窗口,內(nèi)芯頂部也具有光學(xué)窗口,是目前能夠?qū)崿F(xiàn)常規(guī)攻角機(jī)構(gòu)試驗(yàn)光學(xué)條件最好的試驗(yàn)段.由于模型需要安裝在常規(guī)攻角機(jī)構(gòu)上,攻角調(diào)節(jié)范圍是-15°～ 25°,無(wú)法實(shí)現(xiàn)連續(xù)的大攻角變化,需要通過(guò)預(yù)置攻角模塊來(lái)實(shí)現(xiàn)不同階段的攻角調(diào)節(jié).

如圖5,模型為了實(shí)現(xiàn)180°可翻轉(zhuǎn)設(shè)計(jì)安裝,頭部和尾部均可與支撐段連接,因此,頭部需要設(shè)計(jì)可拆卸段,適合尾部的支撐連接.

圖5 模型的支撐設(shè)計(jì)Fig.5 Support design of experimental model

由前述風(fēng)洞試驗(yàn)條件可知,常規(guī)試驗(yàn)段和攻角機(jī)構(gòu)的攻角調(diào)節(jié)范圍是-15°～ 25°,流動(dòng)顯示試驗(yàn)?zāi)Ｐ蜑榱嗽陲L(fēng)洞中實(shí)現(xiàn)0°～ 180°的攻角變化,需要對(duì)模型的攻角進(jìn)行預(yù)置.在0°～ 30°范圍內(nèi),預(yù)置-15°的攻角,在60°～ 90°范圍內(nèi),預(yù)置-75°的攻角,配合模型頭尾更換,可以實(shí)現(xiàn)模型0°～ 180°的攻角變化.

傳統(tǒng)油流試驗(yàn),一般在風(fēng)洞運(yùn)行結(jié)束后,對(duì)表面油流圖像進(jìn)行拍攝,獲取表面流動(dòng)結(jié)構(gòu).在高速、大攻角條件下,較強(qiáng)的流動(dòng)分離、風(fēng)洞臺(tái)階波干擾、油流自身重力等時(shí)刻在影響著表面流動(dòng)圖譜,而暫沖式風(fēng)洞關(guān)車(chē)時(shí)氣流沖擊對(duì)表面流動(dòng)影響非常大,風(fēng)洞運(yùn)行結(jié)束后的表面流動(dòng)圖譜可能與風(fēng)洞運(yùn)行時(shí)的表面流動(dòng)圖譜純?cè)诿黠@差別.因此,本試驗(yàn)采用實(shí)時(shí)化表面油流圖像拍攝方法,在風(fēng)洞運(yùn)行時(shí),按照下圖布局進(jìn)行拍攝.同時(shí),選取特定工況進(jìn)行精細(xì)化的熒光油膜摩阻測(cè)量試驗(yàn),采用紫外光源的照明系統(tǒng)、混有熒光指示劑的油膜、帶有高通濾波片的圖像采集系統(tǒng),針對(duì)模型局部區(qū)域進(jìn)行試驗(yàn)拍攝,見(jiàn)圖6.

圖6 油流流動(dòng)顯示試驗(yàn)測(cè)量布局Fig.6 Oil-flow visualization test measurement layout

2 細(xì)長(zhǎng)旋成體外形及網(wǎng)格

現(xiàn)代高機(jī)動(dòng)導(dǎo)彈除了細(xì)長(zhǎng)彈身外,還具有鴨舵、尾翼等控制部件,Wang 等[41]的研究結(jié)果顯示,對(duì)于細(xì)長(zhǎng)導(dǎo)彈大攻角流動(dòng),其非定常特性、頻率特性主要由細(xì)長(zhǎng)彈身誘導(dǎo),彈翼貢獻(xiàn)較小,所以為了能夠更清晰的研究超大攻角下的流動(dòng),本文采用的細(xì)長(zhǎng)旋成體外形.如圖7,彈身直徑為60 mm,彈身全長(zhǎng)為707.7 mm,彈身頭部采用的是尖型頭部設(shè)計(jì),坐標(biāo)設(shè)置為X 軸向后,Z 軸向上.

圖7 細(xì)長(zhǎng)旋成體外形 (單位:mm)Fig.7 Slender revolutionary body (unit:mm)

計(jì)算網(wǎng)格采用結(jié)構(gòu)網(wǎng)格進(jìn)行求解,如圖8,采用O 型網(wǎng)格拓?fù)?外邊界距離的設(shè)置以能夠滿足自由來(lái)流條件為判斷準(zhǔn)則.貼近物面的第一層網(wǎng)格厚度保持Y+～ 1,以確保邊界層的準(zhǔn)確模擬.為了能夠精確的捕捉細(xì)長(zhǎng)旋成體背風(fēng)側(cè)的分離渦流動(dòng),在模型的法向進(jìn)行了整體加密,在模型背風(fēng)側(cè)展向網(wǎng)格進(jìn)行適當(dāng)加密,同時(shí)當(dāng)攻角超過(guò)90°時(shí),底部超前,這時(shí)底部引起的回流區(qū)模擬也非常重要,此時(shí)會(huì)在尾部附近進(jìn)行加密,總網(wǎng)格量為1500 萬(wàn).

圖8 細(xì)長(zhǎng)體網(wǎng)格Fig.8 Mesh of slender revolutionary body

由于本文針對(duì)大攻角下非定常流動(dòng)特性問(wèn)題,其對(duì)旋渦尺度的分辨以及旋渦頻率特性的捕捉十分重要,所以需要對(duì)網(wǎng)格與時(shí)間步長(zhǎng)無(wú)關(guān)性進(jìn)行檢查,即確定數(shù)值模擬結(jié)果不再隨網(wǎng)格規(guī)模與時(shí)間步長(zhǎng)大小而變化.因此本文對(duì)典型工況(Ma=0.6,α=60°)進(jìn)行模擬,并選取升力CL與阻力Cd的平均值,以及側(cè)向力的斯特勞哈爾數(shù)St 進(jìn)行分析.首先進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證,采用了3 套計(jì)算網(wǎng)格,3 套網(wǎng)格在流向、周向和法向均進(jìn)行變化,網(wǎng)格量分別為530 萬(wàn)、1500 萬(wàn)和4200 萬(wàn),圖9 展示了3 套網(wǎng)格的物面網(wǎng)格分布.表1 為3 套網(wǎng)格下的升力CL、阻力Cd、側(cè)向力St 計(jì)算結(jié)果.從表1 我們可以看到,網(wǎng)格規(guī)模對(duì)CL,Cd與St 均有一定的影響,隨著網(wǎng)格量的增加,當(dāng)網(wǎng)格量達(dá)到1500 萬(wàn)時(shí),其計(jì)算結(jié)果(CL,Cd,St)與網(wǎng)格量4200 萬(wàn)較為接近,證明網(wǎng)格量1500 萬(wàn)對(duì)于研究該問(wèn)題已滿足要求.

表1 網(wǎng)格無(wú)關(guān)性Table 1 Sensitivity of mesh resolution

圖9 3 套不同規(guī)模網(wǎng)格示意Fig.9 Three grids with different cell counts

在時(shí)間步長(zhǎng)無(wú)關(guān)性方面,分別采用物理時(shí)間步長(zhǎng)為dt=1.2 × 10-5s,2.3 × 10-5s,5 × 10-5s 進(jìn)行計(jì)算,網(wǎng)格量為1500 萬(wàn).表2 為不同時(shí)間步長(zhǎng)下的計(jì)算那結(jié)果,從表2 中我們可以看到,時(shí)間步長(zhǎng)的變化對(duì)升力CL、阻力Cd影響不大,但對(duì)St 影響相對(duì)較大,隨著時(shí)間步長(zhǎng)不斷減小,當(dāng)時(shí)間步長(zhǎng)達(dá)到dt=2.3 × 10-5s 時(shí),CL,Cd與St 基本不隨時(shí)間步長(zhǎng)而變化,證明時(shí)間步長(zhǎng)dt=2.3 × 10-5s 能夠滿足研究要求.

表2 時(shí)間步長(zhǎng)無(wú)關(guān)性Table 2 Sensitivity of physical time step

3 超大攻角瞬時(shí)流動(dòng)特性

本文采用CFD 數(shù)值模擬手段,計(jì)算了細(xì)長(zhǎng)旋成體在攻角α=0°～ 180°下的背風(fēng)側(cè)流動(dòng)形態(tài),側(cè)滑角β=0°,雷諾數(shù)Re=1.29 × 106m-1,參考長(zhǎng)度為直徑D=60 mm.由于細(xì)長(zhǎng)體在大攻角下具有較強(qiáng)的非定常效應(yīng),所以這里采用數(shù)值模擬得到的瞬時(shí)流場(chǎng)與試驗(yàn)進(jìn)行對(duì)比.

圖10 為攻角30°,CFD 模擬的表面極限流線分布,以及CFD 與油流試驗(yàn)圖像比較,從圖中我們可以看到,油流試驗(yàn)與數(shù)值模擬得到的主分離線位置較為一致.主分離線在模型頭尖部起始并較為靠近背風(fēng)側(cè),隨著流動(dòng)沿軸向向后發(fā)展,分離線周向位置逐漸向迎風(fēng)側(cè)移動(dòng),當(dāng)流動(dòng)發(fā)展至彈身圓柱段時(shí),分離線周向位置沿著軸向幾乎不變.

圖10 CFD 物面流線與油流試驗(yàn)對(duì)比(α=30°)Fig.10 Surface streamlines comparisons between CFD and oil-flow test (α=30°)

圖11 是攻角30°是旋成體背風(fēng)側(cè)流動(dòng)及物面流線,其中背風(fēng)側(cè)流動(dòng)通過(guò)采用不同截面下的渦量值來(lái)體現(xiàn).從圖中可以看到,旋成體背風(fēng)側(cè)主要呈現(xiàn)出了對(duì)稱(chēng)渦的流動(dòng)形態(tài),旋渦從頭部起始,一直沿軸向向后發(fā)展至尾部.并且在主旋渦下方靠近物面的位置還存在二次渦流動(dòng)結(jié)構(gòu).從物面流線分布也能夠看到明顯的主分離線以及二次分離線.圖12 是X/L=0.1,0.2,0.4,0.6,0.8 截面的壓力分布情況,圖中顯示5 個(gè)截面下壓力分布均為對(duì)稱(chēng)分布形式,體現(xiàn)出旋成體在該攻角的對(duì)稱(chēng)旋渦流動(dòng)結(jié)構(gòu).

圖11 物面流線及空間流動(dòng)(α=30°)Fig.11 Surface streamlines and flow (α=30°)

圖12 各截面壓力分布(α=30°)Fig.12 Pressure distribution of cross sections (α=30°)

圖13 是CFD 物面流線與油流試驗(yàn)對(duì)比圖,圖中有4 張圖,從左至右分別為油流試驗(yàn)拍攝的照片(圖13(a)左)、通過(guò)試驗(yàn)照片處理出的物面流線圖(圖13(a)右)、CFD 手段獲得物面流線圖(圖13(b)左)、CFD 獲得的橫截面流線圖(圖13(b)右).從圖中試驗(yàn)獲得的物面流線圖(圖13(a)左)與CFD 獲得的物面流線圖(圖13(b)左)進(jìn)行對(duì)比發(fā)現(xiàn),存在清晰的主分離線、二次分離線、二次再附線.試驗(yàn)與CFD 獲得分離位置與再附位置對(duì)比較好.圖13(b)也展示了橫截面流線,從圖中可以看到較為清晰的主分離渦流動(dòng)結(jié)構(gòu),同時(shí)在物面附近存在較小的二次渦流動(dòng)結(jié)構(gòu).

圖13 CFD 物面流線與油流試對(duì)比(α=60°)Fig.13 Surface streamlines comparisons between CFD and oil-flow test (α=60°)

圖14 是攻角在α=60°下的背風(fēng)側(cè)旋渦流動(dòng)形態(tài),圖中可以看到,在旋成體的前半段(約X/L=0.4 之前),其旋渦流動(dòng)還較為對(duì)稱(chēng),隨著旋渦流動(dòng)沿軸向向后發(fā)展,背風(fēng)側(cè)旋渦流動(dòng)逐漸轉(zhuǎn)化為非對(duì)稱(chēng)旋渦流動(dòng)形態(tài).當(dāng)達(dá)到X/L=0.6 截面左右時(shí),旋渦流動(dòng)變的不再集中,逐漸發(fā)展為旋渦脫落的流動(dòng)形態(tài),此時(shí)二次旋渦也變的不再穩(wěn)定,在物面上也觀察不到較為明顯的二次分離線.圖15 是各截面的壓力分布情況,從圖中也可以看到,在截面X/L=0.1～ 0.2,其壓力分布較為對(duì)稱(chēng).在X/L=0.4 位置,旋渦流動(dòng)開(kāi)始出現(xiàn)輕微的非對(duì)稱(chēng)形態(tài),隨著旋渦向后發(fā)展,非對(duì)稱(chēng)愈加明顯清晰.

圖14 物面流線及空間流動(dòng)(α=60°)Fig.14 Surface streamlines and flow (α=60°)

圖15 各截面壓力分布(α=60°)Fig.15 Pressure distribution of cross sections (α=60°)

圖16 是攻角α=90°工況下的流動(dòng)圖像.圖中我們可以看到,由于此時(shí)旋成體軸線與來(lái)流成垂直關(guān)系,在背風(fēng)側(cè)的軸向流動(dòng)速度很低,所以背風(fēng)側(cè)無(wú)法形成集中渦流動(dòng)結(jié)構(gòu),在不同截面下基本都體現(xiàn)出了旋渦脫落的流動(dòng)現(xiàn)象,并且具有較強(qiáng)的非對(duì)稱(chēng)性.從壓力分布中(圖17)我們也可以看到,由于背風(fēng)側(cè)旋渦流動(dòng)主要是渦脫落流動(dòng)形態(tài),旋渦在物面附近形成的吸力較弱,沒(méi)有出現(xiàn)吸力峰值.

圖16 物面流線及空間流動(dòng)(α=90°)Fig.16 Surface streamlines and flow (α=90°)

圖17 各截面壓力分布(α=90°)Fig.17 Pressure distribution of cross sections (α=90°)

圖18 是攻角α=120°下油流試驗(yàn)物面流線與CFD 物面流線的對(duì)比.從圖中可以看到明顯的主分離線,CFD 計(jì)算結(jié)果與油流試驗(yàn)結(jié)果吻合很好,沒(méi)有看到比較明顯的二次分離線在模型底部附近均出現(xiàn)了主分離線“彎曲”的現(xiàn)象,這主要是旋成體底部誘導(dǎo)的分離流動(dòng)的影響.圖18 還展示了旋成體在對(duì)稱(chēng)面的流動(dòng)特性,圖中我們看到,當(dāng)來(lái)流流過(guò)模型的底部時(shí)會(huì)背風(fēng)側(cè)產(chǎn)生一個(gè)回流區(qū),該回流區(qū)在軸向的影響范圍約為1.2D.

圖18 CFD 物面流線與油流試驗(yàn)對(duì)比(α=120°)Fig.18 Surface streamlines comparisons between CFD and oil-flow test (α=120°)

圖19 是α=120°攻角下空間流動(dòng)與物面流線.圖中我們可以看到,當(dāng)自由來(lái)流在底部附近產(chǎn)生回流區(qū)后,隨著流動(dòng)沿軸向向后發(fā)展,背風(fēng)側(cè)流動(dòng)逐漸發(fā)展非對(duì)稱(chēng)旋渦流動(dòng),并伴隨著一定的旋渦脫落現(xiàn)象.圖20 是不同截面下的壓力分布情況,可以看到在不同截面下壓力分布始終為非對(duì)稱(chēng).

圖19 物面流線及空間流動(dòng)(α=120°)Fig.19 Surface streamlines and flow (α=120°)

圖20 各截面壓力分布(α=120°)Fig.20 Pressure distribution of cross sections (α=120°)

圖21 是攻角150°下試驗(yàn)與數(shù)值模擬結(jié)果的對(duì)比情況.從圖中我們可以看到,隨著攻角進(jìn)一步增加,底部分離帶來(lái)的回流區(qū)在模型軸向的影響范圍將進(jìn)一步加大,影響范圍約為2.5D.隨著流動(dòng)沿軸向向后發(fā)展,模型背風(fēng)側(cè)逐漸發(fā)展為旋渦流動(dòng),油流試驗(yàn)中的分離線逐漸清晰起來(lái).圖22 是在攻角150°下的空間流動(dòng),從圖中我們可以看到,旋渦流動(dòng)在回流區(qū)后開(kāi)始逐漸發(fā)展起來(lái),在流動(dòng)剛流過(guò)回流區(qū)時(shí),旋渦渦位較低,并且仍然受底部回流區(qū)的影響.隨著流動(dòng)繼續(xù)向后發(fā)展,非對(duì)稱(chēng)旋渦流動(dòng)越來(lái)越清晰起來(lái),受底部回流區(qū)的影響逐漸減弱.

圖21 CFD 物面流線與油流試驗(yàn)對(duì)比(α=150°)Fig.21 Surface streamlines comparisons between CFD and oil-flow test (α=150°)

圖22 物面流線及空間流動(dòng)(α=150°)Fig.22 Surface streamlines and flow (α=150°)

圖23 是攻角180°工況下物面流線與油流試驗(yàn)的對(duì)比情況,圖中一共分為上下兩張附圖,上面附圖為CFD 得到物面流線以及對(duì)稱(chēng)面的渦量分布圖,下面附圖為油流試驗(yàn)照片,在此攻角下主要流動(dòng)特性為底部分離引起的復(fù)雜回流流動(dòng).從試驗(yàn)以及數(shù)值模擬結(jié)果我們都可以看到,在底部附近會(huì)出現(xiàn)一個(gè)十分明顯的分離線,并且回流區(qū)再附位置主要在距離底部約2.5D 的位置,試驗(yàn)與計(jì)算結(jié)果對(duì)比較好.同時(shí)在對(duì)稱(chēng)面渦量分布中我們可以看到,來(lái)流經(jīng)過(guò)底部時(shí)會(huì)產(chǎn)生較強(qiáng)的剪切層,隨著流動(dòng)向后發(fā)展逐漸失穩(wěn),與回流區(qū)中的小尺度旋渦流動(dòng)相互作用,形成較為復(fù)雜的具有較強(qiáng)非線性的局部流動(dòng).

圖23 CFD 物面流線與油流試驗(yàn)對(duì)比(α=180°)Fig.23 Surface streamlines comparisons between CFD and oil-flow test (α=180°)

4 超大攻角旋渦流動(dòng)時(shí)序演化

上一節(jié)對(duì)超大攻角瞬時(shí)流動(dòng)進(jìn)行分析,但細(xì)長(zhǎng)體在超大攻角下體現(xiàn)為復(fù)雜的非定常流動(dòng)現(xiàn)象,其背風(fēng)側(cè)旋渦的強(qiáng)度、位置、渦脫落均會(huì)隨著時(shí)間產(chǎn)生變化,所以本節(jié)對(duì)超大攻角旋渦流動(dòng)的時(shí)序演化進(jìn)行分析.本文每個(gè)計(jì)算物理時(shí)間步為Δt=2.3 ×10-5s,共計(jì)算2000 個(gè)物理時(shí)間步,總計(jì)算物理時(shí)間為0.046 s,以保證足夠的時(shí)長(zhǎng)使得流動(dòng)進(jìn)入穩(wěn)態(tài),并選取后1200 步進(jìn)行分析.

圖24 是攻角α=60°下在X=350 mm 截面背風(fēng)側(cè)渦量圖隨時(shí)間的變化.從圖中我們可以看到,在該截面下背風(fēng)側(cè),呈現(xiàn)出了兩個(gè)集中渦以及兩個(gè)二次渦,集中渦的渦量和渦位隨著時(shí)間增加較小,呈現(xiàn)出輕微的非對(duì)稱(chēng)流動(dòng)形態(tài).隨著流動(dòng)沿著軸向向后發(fā)展,當(dāng)流動(dòng)發(fā)展至X=550 mm 截面,如圖25,集中渦體現(xiàn)出明顯的非對(duì)稱(chēng)流動(dòng)特性,同時(shí)隨著時(shí)間的變化,體現(xiàn)出旋渦脫落特性,產(chǎn)生類(lèi)卡門(mén)渦街式的旋渦交替脫落現(xiàn)象.

圖24 截面渦量隨時(shí)間變化(α=60°,X=350 mm)Fig.24 Vortice of cross section (α=60°,X=350 mm)

圖25 截面渦量隨時(shí)間變化(α=60°,X=550 mm)Fig.25 Vortice of cross section (α=60°,X=350 mm)

圖26 是攻角α=90°下在X=350 mm 截面下的渦量圖隨時(shí)間變化.從圖中我們可以看到,此時(shí)彈身垂直于來(lái)流,來(lái)流流過(guò)彈身后會(huì)產(chǎn)生較長(zhǎng)的剪切層,隨著剪切層逐漸發(fā)展,慢慢卷起旋渦,旋渦強(qiáng)度相對(duì)較弱,并且旋渦距離物面較遠(yuǎn),旋渦對(duì)物面的誘導(dǎo)能力較低,沒(méi)有發(fā)現(xiàn)明顯的二次渦結(jié)構(gòu).彈身兩側(cè)剪切層誘導(dǎo)的旋渦呈現(xiàn)較為明顯的非對(duì)稱(chēng)性,并且隨著時(shí)間增加,旋渦逐漸脫離剪切層向后脫落,左右旋渦的渦脫落現(xiàn)象交替進(jìn)行.

圖26 截面渦量隨時(shí)間變化(α=90°,X=350 mm)Fig.26 Vortice of cross section (α=90°,X=350 mm)

圖27 是攻角α=120°下的X=350 mm 截面渦量圖.從圖中我們可以看到,此時(shí)旋渦距離物面較近,出現(xiàn)了明顯的非對(duì)稱(chēng)旋渦流動(dòng)結(jié)構(gòu),同時(shí)在物面附近會(huì)誘導(dǎo)出較小尺度的二次渦結(jié)構(gòu),隨著時(shí)間發(fā)展,左右旋渦的渦位會(huì)出現(xiàn)切換現(xiàn)象,渦位較高的旋渦逐漸脫落并飄向尾跡區(qū).圖28 是Y=0 mm 的渦量分布圖,圖中可以看到,在底部附近的回流區(qū)存在較多的小尺度旋渦結(jié)構(gòu),隨著時(shí)間增加,距離物面較近小尺度渦在物面附近形成回流區(qū),距離物面較遠(yuǎn)的小尺度渦慢慢脫落飄向尾跡區(qū).彈身(圓柱等直段)附近也存在較強(qiáng)的渦量分布,這是圓柱段集中渦渦量在Y=0 mm 平面的投影,反應(yīng)出了超大攻角下旋渦具有復(fù)雜的三維效應(yīng).

圖27 截面渦量隨時(shí)間變化(α=120°,X=350 mm)Fig.27 Vortice of cross section (α=120°,X=350 mm)

圖28 截面渦量隨時(shí)間變化(α=120°,Y=0 mm)Fig.28 Vortice of cross section (α=120°,Y=0 mm)

圖29 是攻角α=150°下的X=350 mm 截面渦量圖.圖中我們可以看到,旋成體背風(fēng)側(cè)產(chǎn)生了兩個(gè)集中渦,渦位較低,隨著時(shí)間增加,旋成體背風(fēng)側(cè)的渦位出現(xiàn)左右切換的現(xiàn)象,并且隨著時(shí)間的增加逐漸脫體向下游運(yùn)動(dòng).圖30 是Y=0 mm 截面渦量圖,來(lái)流流過(guò)底部后會(huì)出現(xiàn)較為明顯的剪切層,剪切層逐漸向后發(fā)展逐漸失穩(wěn)變成小尺度渦飄向下游,同時(shí)在底部回流區(qū)內(nèi)依然為很多小尺度旋渦流動(dòng).

圖29 截面渦量隨時(shí)間變化(α=150°,X=350 mm)Fig.29 Vortice of cross section (α=150°,X=350 mm)

圖30 截面渦量隨時(shí)間變化(α=150°,Y=0 mm)Fig.30 Vortice of cross section (α=150°,Y=0 mm)

圖31 是攻角α=180°時(shí)在Y=0 mm 截面的渦量分布圖,從圖中我們可以看到,來(lái)流會(huì)在底部處拖出明顯的剪切層,剪切層逐漸向后發(fā)展逐漸失穩(wěn),一部分渦量會(huì)變成小尺度旋渦向下游發(fā)展,另一部分會(huì)再附至物面,在頭部形成一個(gè)回流區(qū),回流區(qū)內(nèi)的小尺度旋渦相互干擾作用,形成復(fù)雜的非定?；亓髁鲃?dòng)結(jié)構(gòu).

圖31 截面渦量隨時(shí)間變化(α=180°,Y=0 mm)Fig.31 Vortice of cross section (α=180°,Y=0 mm)

5 超大攻角物面壓力脈動(dòng)特性

為了更好的對(duì)細(xì)長(zhǎng)體超大攻角下非定常流動(dòng)頻率特性進(jìn)行分析,如圖32 所示,取兩個(gè)主要截面,一個(gè)是X=350 mm 截面,位于模型中間位置,能夠較好的描述背風(fēng)側(cè)由圓柱引起的旋渦流動(dòng);另一個(gè)是X=550 mm 截面,能夠較好的描述模型底部附近流動(dòng).在這兩個(gè)平面內(nèi),分別取了一個(gè)檢測(cè)點(diǎn)(m1 點(diǎn)和m2 點(diǎn)),通過(guò)提取檢測(cè)點(diǎn)的壓力脈動(dòng),可以更深入的分析非定常流動(dòng)的頻域特性.

圖32 監(jiān)測(cè)點(diǎn)及截面位置示意圖Fig.32 Location of monitored point and cross section

圖33 是攻角α=60°下監(jiān)測(cè)點(diǎn)m2 的壓力脈動(dòng),在此截面流動(dòng)已經(jīng)體現(xiàn)出了明顯的非定常特性,從測(cè)壓點(diǎn)的壓力脈動(dòng)就可以看出來(lái),通過(guò)對(duì)m2 點(diǎn)壓力脈動(dòng)進(jìn)行頻譜分析,發(fā)現(xiàn)壓力脈動(dòng)具有主頻,主頻為f=543 Hz.

圖33 監(jiān)測(cè)點(diǎn)壓力脈動(dòng)(α=60°,X=550 mm)Fig.33 Fluctuating pressure of monitored point (α=60°,X=550 mm)

圖34 是兩個(gè)監(jiān)測(cè)點(diǎn)(m1 點(diǎn)和m2 點(diǎn))的壓力系數(shù)脈動(dòng)數(shù)據(jù).從圖中可以看到,壓力脈動(dòng)較為明顯,體現(xiàn)出了明顯的非定常性,通過(guò)壓力脈動(dòng)進(jìn)行頻譜分析,發(fā)現(xiàn)在攻角α=90°下沒(méi)有發(fā)現(xiàn)明顯的主頻現(xiàn)象,這主要是由于該攻角下,旋渦距離物面較遠(yuǎn),對(duì)物面誘導(dǎo)能力較弱的原因.

圖34 監(jiān)測(cè)點(diǎn)壓力脈動(dòng)(α=90°)Fig.34 Fluctuating pressure of monitored point (α=90°)

圖35 展示了監(jiān)測(cè)點(diǎn)m1 和m2 的壓力脈動(dòng)情況,m1 點(diǎn)的壓力脈動(dòng)體現(xiàn)出了圓柱段旋渦的切換情況,從頻譜分析中可以看到,該點(diǎn)的壓力脈動(dòng)有一定的主頻,主頻為f=717 Hz 左右,但是該主頻并沒(méi)有特別集中,這里面也體現(xiàn)出了圓柱段的旋渦結(jié)構(gòu)不僅僅存在較大尺度的集中渦結(jié)構(gòu),同時(shí)也存在較小尺度的旋渦,小尺度渦與大尺度渦相互干擾影響,造成一定程度上的頻率分散.m2 點(diǎn)體現(xiàn)出了底部分離回流區(qū)的壓力脈動(dòng)特性(如圖35(b)),可以看到在m2 點(diǎn)在較低頻率(1000 Hz 以內(nèi))沒(méi)有發(fā)現(xiàn)主頻,但是在高頻段出現(xiàn)了較為明顯的主頻,主頻為5391 Hz左右,這主要是在該截面的壓力脈動(dòng)主要由底部分離回流區(qū)中的小尺度旋渦帶來(lái),小尺度旋渦頻率較高,對(duì)應(yīng)壓力脈動(dòng)頻率也較高.

圖35 監(jiān)測(cè)點(diǎn)壓力脈動(dòng)(α=120°)Fig.35 Fluctuating pressure of monitored point (α=120°)

圖36 為α=150°監(jiān)測(cè)點(diǎn)m1 和m2 的壓力脈動(dòng)情況,圖中可以看到,集中渦所產(chǎn)生的脈動(dòng)存在主頻,主頻為543 Hz 左右.從監(jiān)測(cè)點(diǎn)m2 的頻譜分析來(lái)看,類(lèi)似α=120°所發(fā)現(xiàn)的高頻脈動(dòng),在α=150°下的高頻脈動(dòng)為5086 Hz 左右,但是此壓力脈動(dòng)能量較低.

圖36 監(jiān)測(cè)點(diǎn)壓力脈動(dòng)(α=150°)Fig.36 Fluctuating pressure of monitored point (α=150°)

圖37 是攻角α=180°下m2 點(diǎn)壓力脈動(dòng),用于監(jiān)測(cè)底部回流區(qū)的壓力脈動(dòng)情況.從圖中可以看到,底部回流區(qū)引起了較高頻的壓力脈動(dòng),脈動(dòng)頻率為f=5173 Hz.

圖37 監(jiān)測(cè)點(diǎn)壓力脈動(dòng)(α=180°)Fig.37 Fluctuating pressure of monitored point (α=180°)

從以上分析可以看到,不同攻角下,超大攻角非定常流動(dòng)體現(xiàn)出了不同的物面壓力脈動(dòng)特性,主要有兩類(lèi)主頻,一類(lèi)為圓柱段非對(duì)稱(chēng)旋渦誘導(dǎo)的較低頻率,另一類(lèi)為攻角超過(guò)90°時(shí),底部分離流動(dòng)誘導(dǎo)的較高的頻率.本文將特征頻率進(jìn)行無(wú)量綱化,提取出脈動(dòng)主頻的斯特勞哈爾數(shù)St

式中f 為流動(dòng)頻率,D 為直徑,V 為特征速度.對(duì)于非對(duì)稱(chēng)旋渦誘導(dǎo)的頻率,類(lèi)比圓柱擾流流動(dòng),特征速度取橫截面流動(dòng)速度V=V∞sinα,對(duì)于底部分離流動(dòng)誘導(dǎo)的流動(dòng)頻率,特征速度取V=V∞.圖38 是不同攻角下的流動(dòng)St,圖中我們可以看到,由非對(duì)稱(chēng)渦引起的物面壓力脈動(dòng)主要集中在攻角60°～ 150°,其中流動(dòng)主頻St 范圍為St=0.19～ 0.33,而由底部分離誘導(dǎo)的非定常流動(dòng)在120°～ 180°都有主頻,其脈動(dòng)頻率較高,無(wú)量綱主頻范圍為St=1.55～ 1.64.

圖38 不同攻角下的非定常流動(dòng)頻率St 數(shù)Fig.38 St of unsteady flow versus angle of attack

6 結(jié)論

本文結(jié)合數(shù)值模擬DDES 方法與油流顯示試驗(yàn)方法,研究了細(xì)長(zhǎng)旋成體在Ma=0.6 時(shí)超大攻角范圍內(nèi)(α=0°～ 180°)的流動(dòng)特性,得到以下結(jié)論.

(1)數(shù)值模擬與油流顯示試驗(yàn)獲得的物面流線吻合較好,在攻角范圍α=30°～ 60°和α=150°背風(fēng)側(cè)旋渦較為集中并且離物面較近時(shí),均能捕捉到明顯的分離線、再附線等物面流動(dòng)拓?fù)湫问?在攻角范圍α=90°～ 120°,背風(fēng)側(cè)旋渦距離物面較遠(yuǎn),物面流動(dòng)拓?fù)湟砸淮畏蛛x線為主,當(dāng)攻角α=180°時(shí),物面的分離線與再附線主要由底部分離回流區(qū)產(chǎn)生.

(2)在攻角范圍α=0°～ 90°,細(xì)長(zhǎng)體背風(fēng)側(cè)流動(dòng)特性主要體現(xiàn)為由圓柱段產(chǎn)生的分離渦流動(dòng)結(jié)構(gòu),當(dāng)攻角α=30°時(shí),分離渦位對(duì)稱(chēng)集中渦,攻角至α=60°時(shí),細(xì)長(zhǎng)體頭部附近依然為對(duì)稱(chēng)渦,在尾部附近非定常非對(duì)稱(chēng)渦結(jié)構(gòu),攻角至90°時(shí),細(xì)長(zhǎng)體背風(fēng)側(cè)均為非對(duì)稱(chēng)旋渦結(jié)構(gòu),并伴有多尺度旋渦干擾、渦脫落等現(xiàn)象.

(3)攻角α=120°～ 150°范圍內(nèi),細(xì)長(zhǎng)體底部朝前,底部產(chǎn)生較為明顯的回流區(qū),回流區(qū)內(nèi)存在較多小尺度旋渦相互作用干擾,隨著流動(dòng)逐漸沿軸向向后發(fā)展,流動(dòng)逐漸發(fā)展為非對(duì)稱(chēng)渦流動(dòng).當(dāng)攻角達(dá)到α=180°時(shí),來(lái)流經(jīng)過(guò)底部時(shí)會(huì)產(chǎn)生較強(qiáng)的剪切層,隨著流動(dòng)向后發(fā)展逐漸失穩(wěn),與回流區(qū)中的小尺度旋渦流動(dòng)相互作用.

(4)由非對(duì)稱(chēng)渦引起的物面壓力脈動(dòng)主要集中在攻角60°～ 150°,其中流動(dòng)主頻St 范圍為St=0.19～0.33,而由底部分離誘導(dǎo)的非定常流動(dòng)在120°～180°都有主頻,其脈動(dòng)頻率較高,無(wú)量綱主頻范圍為St=1.55～ 1.64.

本文此次研究旨在對(duì)超大攻角下細(xì)長(zhǎng)體的流動(dòng)特性有個(gè)較為清晰的認(rèn)識(shí),下一步將對(duì)非定常氣動(dòng)特性對(duì)機(jī)動(dòng)飛行的影響及其控制策略做進(jìn)一步研究.