鄭誠(chéng)毅, 杜旭之, 東喬天, 楊志剛 熊 兵, 徐 毅, 吳凌昊, 金哲巖,

(1. 同濟(jì)大學(xué) 航空航天與力學(xué)學(xué)院,上海200092; 2. 中國(guó)商用飛機(jī)上海飛機(jī)設(shè)計(jì)研究院,上海 201203; 3. 上海市地面交通工具空氣動(dòng)力與熱環(huán)境模擬重點(diǎn)實(shí)驗(yàn)室, 上海 201804;4. 中國(guó)航發(fā)四川燃?xì)鉁u輪研究院,四川 綿陽(yáng),621000)

結(jié)冰會(huì)改變飛機(jī)的氣動(dòng)外形,影響飛機(jī)的飛行安全.當(dāng)飛機(jī)飛行在一個(gè)低溫且濕度大的環(huán)境中,未凍結(jié)的過(guò)冷水滴會(huì)撞擊飛行器表面,通過(guò)質(zhì)量和熱量交換形成冰的吸積.飛機(jī)的大部分部件都有可能發(fā)生結(jié)冰,比如機(jī)翼、尾翼、發(fā)動(dòng)機(jī)進(jìn)氣口、天線罩等,嚴(yán)重影響飛機(jī)的氣動(dòng)性能,危害飛機(jī)的飛行安全[1-4].近年來(lái),過(guò)冷大水滴 (Supercooled Large Droplet, SLD) 產(chǎn)生的結(jié)冰問題越來(lái)越受到人們的關(guān)注.由于SLD的飛濺、反彈等特殊現(xiàn)象,液滴會(huì)越過(guò)防除冰系統(tǒng)的極限,形成更容易造成流動(dòng)分離的冰脊,并導(dǎo)致飛機(jī)過(guò)早失速[5].粗糙冰是在過(guò)冷水滴積累過(guò)程中附著在翼型前緣表面,隨著時(shí)間變化而不斷增長(zhǎng)的一種特殊冰形.粗糙冰的粗糙度會(huì)直接影響到翼型表面的局部對(duì)流換熱系數(shù)和摩擦系數(shù),它們反過(guò)來(lái)又會(huì)影響到最終冰形.除此之外,粗糙冰會(huì)顯著改變翼型輪廓和表面的無(wú)黏流場(chǎng),破壞翼型表面的附著流并影響翼型整體的氣動(dòng)性能.研究人員已對(duì)不同形狀的冰形(角冰、流向冰和展向冰)及其對(duì)翼型繞流流場(chǎng)的影響做了大量研究,卻缺乏對(duì)粗糙冰存在時(shí)翼型繞流流場(chǎng)的實(shí)驗(yàn)研究.因此使用SLD結(jié)冰條件下產(chǎn)生的粗糙冰的粗糙度,并且以實(shí)驗(yàn)手段來(lái)探究粗糙冰對(duì)翼型繞流流場(chǎng)結(jié)構(gòu)的影響,對(duì)機(jī)翼的防除冰設(shè)計(jì)來(lái)說(shuō)具有重要的研究?jī)r(jià)值及實(shí)際意義.

文獻(xiàn)[6-7]中研究了冰的粗糙度特征并將帶冰機(jī)翼前緣劃分為3個(gè)主要區(qū)域: 光滑區(qū)、粗糙區(qū)和羽毛區(qū).研究發(fā)現(xiàn)即使在結(jié)冰產(chǎn)生的早期階段,粗糙冰的特征高度也大于局部邊界層的厚度,對(duì)流動(dòng)產(chǎn)生破壞.文獻(xiàn)[8-10]中詳細(xì)研究了獨(dú)立分布的冰粗糙元對(duì)邊界層轉(zhuǎn)變的影響,結(jié)果表明邊界層會(huì)在冰粗糙元后緩慢過(guò)渡成湍流邊界層.Smith等[11]用具有相同氣動(dòng)效應(yīng)的均勻沙粒代替冰的粗糙度,發(fā)現(xiàn)粗糙元的密度對(duì)翼型的氣動(dòng)性能影響不明顯.Jackson[12]也研究了冰粗糙元大小和密度對(duì)翼型氣動(dòng)力的影響.黃冉冉等[13]對(duì)附有流向冰的翼型氣動(dòng)力進(jìn)行了數(shù)值模擬后發(fā)現(xiàn),在流向冰的冰形表面添加粗糙度會(huì)明顯降低翼型的氣動(dòng)性能.霍西恒等[14]研究發(fā)現(xiàn)結(jié)冰表面粗糙度很大程度增強(qiáng)了表面的對(duì)流換熱效果,影響最終結(jié)冰冰形的幾何特征.常士楠等[15]通過(guò)受力分析,確定了水膜、水珠和細(xì)流特征粗糙度的大小及其之間的界限劃分,發(fā)展和完善了結(jié)冰表面的粗糙度模型.Mcclain等[16]對(duì)過(guò)冷大水滴結(jié)冰條件下,翼型前緣在短時(shí)間內(nèi)積累的粗糙冰特征高度進(jìn)行了測(cè)量和記錄.在各類翼型和葉型關(guān)于粗糙度的研究中,有很多學(xué)者選擇砂紙來(lái)代替表面粗糙度進(jìn)行實(shí)驗(yàn).Milsch[17]采用葉片表面貼砂紙的方法研究了不同粗糙度影響下的二維葉柵損失.Boer等[18]在研究粗糙度的特征高度和在翼型表面分布的位置對(duì)翼型最大升力的影響時(shí),使用了砂紙作為研究對(duì)象.Dees等[19]使用不同目數(shù)的砂紙和規(guī)則排布的粗糙度進(jìn)行了對(duì)比研究,結(jié)果表明兩者在燃?xì)廨啓C(jī)葉片吸力側(cè)的傳熱和表面摩擦系數(shù)方面區(qū)別很小.高磊[20]選取了4個(gè)型號(hào)的砂紙研究了粗糙度量級(jí)對(duì)壓力機(jī)葉片總壓損失攻角特性和氣流折轉(zhuǎn)角的影響.孫爽等[21]使用3種硅碳化物砂紙來(lái)實(shí)現(xiàn)3種表面粗糙度,研究了表面粗糙度對(duì)低壓渦輪葉型損失及附面層特性的影響.

雷諾數(shù)對(duì)翼型氣動(dòng)性能的影響通常隨著冰的粗糙度大小而改變.Hoerner等[22]探索了雷諾數(shù)對(duì)具有前緣粗糙度的NACA0012翼型最大升力的影響,結(jié)果表明雷諾數(shù)的降低會(huì)導(dǎo)致翼型最大升力略微降低.Broeren等[23]研究表明,光滑翼型的最大升力系數(shù)對(duì)雷諾數(shù)和馬赫數(shù)有很強(qiáng)的依賴性,對(duì)于冰翼結(jié)構(gòu),雷諾數(shù)和馬赫數(shù)對(duì)最大升力系數(shù)的影響很小.此外,Rafael等[24]在實(shí)驗(yàn)條件限制下選取1.5×105到2.8×105范圍的雷諾數(shù)研究了不同類型的角冰對(duì)翼型性能的影響.

現(xiàn)有研究對(duì)冰的粗糙度形成機(jī)理[6-7]、幾何特征[6-7, 14]及其對(duì)邊界層的影響[8-13]做了解釋,對(duì)于在實(shí)驗(yàn)中如何表征冰的粗糙度也開發(fā)了相應(yīng)的等效粗糙度模型[15].此外,也有學(xué)者記錄了過(guò)冷大水滴結(jié)冰條件下冰的粗糙度數(shù)值[16].但關(guān)于結(jié)冰初期粗糙冰對(duì)于翼型氣動(dòng)特性影響的研究比較缺乏,已有研究主要以數(shù)值模擬的方式來(lái)進(jìn)行,通過(guò)風(fēng)洞實(shí)驗(yàn)獲得的翼型繞流流場(chǎng)結(jié)構(gòu)數(shù)據(jù)較少.本文使用粒子圖像測(cè)速(PIV)技術(shù)在風(fēng)洞中獲取翼型繞流流場(chǎng)的信息,并以翼型前緣的粗糙冰的粗糙度和攻角作為控制變量,來(lái)探索翼型在過(guò)冷大水滴結(jié)冰條件下,粗糙冰的粗糙度對(duì)不同攻角下翼型繞流流場(chǎng)結(jié)構(gòu)產(chǎn)生的影響.同時(shí)還研究了在0° 攻角和特定粗糙冰條件下,雷諾數(shù)的改變對(duì)翼型繞流流場(chǎng)結(jié)構(gòu)的影響.

1 實(shí)驗(yàn)方法

本項(xiàng)研究所開展的流場(chǎng)測(cè)量實(shí)驗(yàn)是在一個(gè)低速直流式風(fēng)洞內(nèi)進(jìn)行的.該風(fēng)洞具有一個(gè)450 mm×450 mm截面的試驗(yàn)段,試驗(yàn)段的四側(cè)為可替換的光學(xué)玻璃.該風(fēng)洞試驗(yàn)段的上游有一個(gè)收縮段和相應(yīng)的整流段,用來(lái)提供均勻、穩(wěn)定的氣流.實(shí)驗(yàn)所用翼型模型如圖1所示,為不銹鋼材質(zhì)的NACA0012翼型,弦長(zhǎng)150 mm,展長(zhǎng)250 mm,表面經(jīng)過(guò)拋光處理,并通過(guò)一系列支撐結(jié)構(gòu)和攻角結(jié)構(gòu)將翼型模型布置在風(fēng)洞試驗(yàn)段的中心位置.翼型模型前緣黃色部分為安裝砂紙的區(qū)域,砂紙覆蓋機(jī)翼前緣直至33.77 mm弦長(zhǎng)處.本項(xiàng)研究所采用的實(shí)驗(yàn)裝置如圖2所示.同步器連接照相機(jī)、激光器和計(jì)算機(jī),計(jì)算機(jī)通過(guò)控制同步器來(lái)協(xié)調(diào)并驅(qū)動(dòng)激光和照相機(jī)完成拍攝.導(dǎo)光臂能夠?qū)⒓す馄靼l(fā)出的點(diǎn)光源轉(zhuǎn)換成片光源.

圖1 翼型模型Fig.1 Airfoil model

圖2 PIV實(shí)驗(yàn)裝置Fig.2 Experimental facility of PIV

砂紙的表面粗糙度具有弱隨機(jī)性,在某種程度上與冰的粗糙度達(dá)成一致.在現(xiàn)有的等效粗糙度理論支持下,通過(guò)多次測(cè)量求平均的手段得到的砂紙粗糙度可以表征冰的粗糙度.這在實(shí)驗(yàn)中是一種即經(jīng)濟(jì)又有效的方式[17-21].Mcclain等[16]的研究表明,在過(guò)冷大水滴結(jié)冰條件下,弦長(zhǎng)為533 mm的NACA0012翼型模型上短時(shí)(94 s內(nèi))結(jié)冰的粗糙元的峰谷差集中在500 μm以下,且分布在距翼型前緣120 mm的區(qū)域內(nèi).使用表面粗糙度儀器對(duì)不同目數(shù)的砂紙的垂直方向上的粗糙度(Rz)進(jìn)行了測(cè)量,多次測(cè)量求平均值后如表1所示.由于本項(xiàng)研究的翼型弦長(zhǎng)為150 mm,通過(guò)等比縮放并考慮安裝的可行性,最終選擇180目、220目、280目和360目砂紙代替粗糙冰進(jìn)行實(shí)驗(yàn),同時(shí)設(shè)置沒有砂紙的光滑翼型作為對(duì)照組.實(shí)驗(yàn)所用雷諾數(shù)包括8.64×104、1.08×105和1.3×105,實(shí)驗(yàn)攻角包括0°、2°、5°、7°、9° 和12°.經(jīng)測(cè)量,在實(shí)驗(yàn)所使用的的雷諾數(shù)范圍內(nèi),風(fēng)洞來(lái)流湍流度小于0.5%,滿足實(shí)驗(yàn)所需精度.下文中以砂紙表征冰的粗糙度均以砂紙的型號(hào)表述.

表1 砂紙目數(shù)與Rz的關(guān)系Tab.1 Relation between mesh of sand paper and Rz

2 雷諾數(shù)對(duì)R1型粗糙冰作用下翼型流場(chǎng)的影響

取雷諾數(shù)Re=8.64×104,1.08×105,1.3×105,使用180目砂紙來(lái)代替R1型粗糙冰,在翼型的攻角為0° 的條件下進(jìn)行實(shí)驗(yàn).探索了雷諾數(shù)的不同對(duì)翼型繞流流場(chǎng)的流線圖、展向渦量和標(biāo)準(zhǔn)化雷諾應(yīng)力分布的影響.

圖3為R1型粗糙冰作用下不同雷諾數(shù)的流線圖.圖中:X為實(shí)際流場(chǎng)的橫坐標(biāo);Y為實(shí)際流場(chǎng)的縱坐標(biāo);C為翼型弦長(zhǎng).從圖中可以看出各雷諾數(shù)條件下,氣流在翼型前方的流線大體沿來(lái)流方向.氣流流過(guò)翼型上表面的流動(dòng)以附著流為主,未產(chǎn)生分離.說(shuō)明實(shí)驗(yàn)采用的雷諾數(shù)未對(duì)氣流的流動(dòng)分離產(chǎn)生明顯影響.

圖3 R1型粗糙冰作用下不同雷諾數(shù)的流線圖Fig.3 Streamline of R1 rough ice at different Reynolds numbers

圖4為R1型粗糙冰作用時(shí)不同雷諾數(shù)條件下的展向渦量圖.從圖中可以看出渦量在翼型的上表面和前緣并沒有產(chǎn)生.正值渦量產(chǎn)生于翼型尾部的上方,負(fù)值渦量產(chǎn)生于翼型尾部的下方,形成兩條渦量條帶.隨著雷諾數(shù)的增加,渦量條帶在X軸方向上的長(zhǎng)度顯著增加,且愈發(fā)清晰可見.流場(chǎng)中其他區(qū)域的渦量無(wú)明顯區(qū)別.

圖4 R1型粗糙冰作用下不同雷諾數(shù)的渦量圖Fig.4 Spanwise vorticity of R1 rough ice at different Reynolds numbers

圖5所示為不同雷諾數(shù)條件下X/C×100=-101 處尾流截面的渦量值對(duì)比.從圖中可以看出,以翼型弦線為中心,上方的渦量呈現(xiàn)正值,下方的渦量呈現(xiàn)負(fù)值,與渦量云圖相符.不同雷諾數(shù)條件下的渦量值均在Y/C×100=0附近發(fā)生了劇烈波動(dòng),且Re=1.3×105時(shí)渦量的極值明顯大于其他兩個(gè)雷諾數(shù).說(shuō)明雷諾數(shù)不改變?cè)摻孛嫣帨u量在Y軸方向上的變化趨勢(shì),但在一定程度上,雷諾數(shù)的增加會(huì)使尾流的渦量值增加.

圖5 不同雷諾數(shù)條件下X/C×100=-101截面的渦量對(duì)比Fig.5 Spanwise vorticity comparison of section X/C×100=-101 at different Reynolds numbers

圖6所示為Re=1.3×105時(shí)R1型粗糙冰作用下的標(biāo)準(zhǔn)化雷諾應(yīng)力圖.從圖中可以看出,在翼型的周圍無(wú)明顯的雷諾應(yīng)力的產(chǎn)生.在實(shí)驗(yàn)條件下,不同雷諾數(shù)時(shí),流場(chǎng)的標(biāo)準(zhǔn)化雷諾應(yīng)力分布沒有明顯的區(qū)別.

圖6 Re=1.3×105時(shí)R1型粗糙冰作用下的標(biāo)準(zhǔn)化雷諾應(yīng)力圖Fig.6 Normalized Reynolds stress of R1 rough ice at Re=1.3×105

圖7為不同雷諾數(shù)條件下X/C×100=-101 處尾流截面的雷諾應(yīng)力數(shù)值對(duì)比.從圖中可以看到Re=8.64×104時(shí)雷諾應(yīng)力的波動(dòng)最為劇烈,其最小值為-0.026,最大值為0.045,均大于其余兩個(gè)雷諾數(shù)條件下的雷諾應(yīng)力,Re=1.3×105時(shí)波動(dòng)最小.說(shuō)明雷諾數(shù)在翼型前緣覆粗糙冰情況下,影響了該截面處尾流的標(biāo)準(zhǔn)化雷諾應(yīng)力分布,且隨實(shí)驗(yàn)雷諾數(shù)增大,標(biāo)準(zhǔn)化雷諾應(yīng)力減小.

圖7 不同雷諾數(shù)條件下X/C×100=-101截面的雷諾應(yīng)力對(duì)比Fig.7 Normalized Reynolds stress comparison of section X/C×100=-101 at different Reynolds number

3 不同粗糙冰情況下攻角對(duì)翼型流場(chǎng)的影響

使用180目(R1型,Rz=158 μm)、220目(R2型,Rz=131.4 μm)、280目(R3型,Rz=119.4 μm)和360目(R4型,Rz= 94.6 μm)的砂紙代替不同的粗糙冰,選取0°、2°、5°、7°、9° 和12° 共6個(gè)攻角,在Re=1.3×105的條件下進(jìn)行實(shí)驗(yàn),同時(shí)設(shè)置沒有砂紙(Rz=0 μm)的光滑翼型作為對(duì)照組,探索不同粗糙冰情況下攻角對(duì)速度場(chǎng)、展向渦量分布和標(biāo)準(zhǔn)化雷諾應(yīng)力分布的影響.

圖8為0° 攻角時(shí)不同粗糙冰條件下的速度場(chǎng)(v).從圖中可以看出,速度場(chǎng)以翼型為軸對(duì)稱分布.翼型前緣處和后緣尾流的速度最小,上翼面表面的速度最大,并向外逐漸減小.不同粗糙冰條件下的速度場(chǎng)分布非常相近,無(wú)粗糙冰情況下的尾流相對(duì)更加平滑、穩(wěn)定.

圖8 不同粗糙冰條件下的速度場(chǎng)Fig.8 Velocity of different rough ice

表2提取了各粗糙冰作用下翼型上表面氣流速度的最大值和該最大值發(fā)生處的X軸坐標(biāo).從中可以看出,除R3型粗糙冰數(shù)據(jù)以外,粗糙度的存在降低了流場(chǎng)中氣流速度的最大值,且最大值發(fā)生的位置也向翼型尾端移動(dòng).R1型粗糙冰作用下流場(chǎng)中氣流速度的最大值僅為光滑翼型時(shí)的90.2%.流場(chǎng)中氣流速度的最大值越小意味著流場(chǎng)的平均速度越小,因此翼型上表面壓力值也越小.由此可見,粗糙冰的存在會(huì)影響翼型的升力,且粗糙度越大,對(duì)升力的破壞越嚴(yán)重.

通過(guò)對(duì)比各粗糙冰作用下翼型繞流流場(chǎng)的展向渦量分布和標(biāo)準(zhǔn)化雷諾應(yīng)力分布也可以發(fā)現(xiàn)粗糙冰使會(huì)使尾流流動(dòng)變得紊亂且具有旋度的流動(dòng)區(qū)域增大.但在所有粗糙冰中,R4型粗糙冰作用下翼型繞流流場(chǎng)分布與光滑翼型最為相近.因此,引入攻角對(duì)R4型粗糙冰和光滑翼型情況下的流場(chǎng)進(jìn)行下一步分析.

圖9為光滑翼型和R4型粗糙冰作用下不同攻角的流線圖.由圖可見,在0° 和9° 攻角下,氣流在翼型前方的流線大體沿來(lái)流方向.氣流流過(guò)翼型上表面的流動(dòng)以附著流為主,并未產(chǎn)生分離.在12°攻角下,氣流從翼型前緣往后緣運(yùn)動(dòng)的過(guò)程中,不再附著在壁面上,邊界層發(fā)生分離,出現(xiàn)回流并形成了分離泡.圖9很好地反映了翼型隨攻角增大而產(chǎn)生的流動(dòng)分離過(guò)程.表3所示為分離泡特征表.由表可見,帶粗糙冰的翼型氣流分離發(fā)生在X/C×100=-18處,而光滑翼型發(fā)生在X/C×100=-19處,說(shuō)明粗糙冰會(huì)誘導(dǎo)氣流提前分離,且提前量為1%的翼型弦長(zhǎng).同時(shí)可以發(fā)現(xiàn),粗糙冰的存在降低了分離泡的寬度,帶有粗糙冰的分離泡左極限在X/C×100=-121 處,而光滑翼型情況下分離泡左極限在X/C×100=-126 處.結(jié)合氣流的分離點(diǎn)可以發(fā)現(xiàn),粗糙冰存在時(shí)氣流會(huì)提前產(chǎn)生分離,但整體分離泡的寬度小于光滑翼型情況下.對(duì)比兩種情況下流線圖中分離泡的結(jié)構(gòu)可以發(fā)現(xiàn),光滑翼型的分離泡內(nèi)流動(dòng)較規(guī)則,而粗糙冰存在時(shí)分離泡中心區(qū)域的流線有一定程度的波動(dòng).

表3 分離泡特征表Tab.3 Feature of separated bubble

圖9 光滑翼型和R4型粗糙冰作用下不同攻角的流線圖Fig.9 Streamline of clean airfoil and R4 rough ice at different angles of attack

圖10為光滑翼型和R4型粗糙冰作用下12° 攻角時(shí)的渦量云圖.從圖10中可以看出,無(wú)論有無(wú)粗糙冰,正值渦量始于翼型前緣,負(fù)值渦量始于翼型尾緣.通過(guò)提取數(shù)據(jù)發(fā)現(xiàn),光滑翼型情況下,渦量的最小值為-5.33 s-1,發(fā)生在X/C×100=-103.25截面處;渦量的最大值為8.165 s-1,發(fā)生在X/C×100=-4.38 截面處.R4型粗糙冰情況下,渦量最小值為-3.94 s-1,發(fā)生在X/C×100=-102.27 截面處;渦量最大值為7.645 s-1,發(fā)生在X/C×100=-4.38 截面處.由此可見,粗糙冰的存在對(duì)渦量最大值發(fā)生的位置沒有影響,但略微影響了渦量最小值發(fā)生的位置.從渦量最大值和最小值的數(shù)值來(lái)看,粗糙冰存在時(shí),無(wú)論是正值渦量還是負(fù)值渦量均小于光滑翼型情況下的具體數(shù)值.

圖10 光滑翼型和R4型粗糙冰作用下12° 攻角時(shí)的渦量Fig.10 Spanwise vorticity of clean airfoil and R4 rough ice at a 12° angle of attack

圖11給出了不同流場(chǎng)截面處,光滑翼型和R4型粗糙冰的渦量值對(duì)比.在圖11(a)所述截面下,光滑翼型尾流的渦量絕對(duì)值明顯大于粗糙冰存在時(shí).由圖11(b)可見,在X/C×100=-20 截面處兩組數(shù)據(jù)的渦量分布差別不大.由此可見在氣流發(fā)生大分離的情況下,粗糙冰不改變翼型前緣產(chǎn)生的正值渦量在Y軸方向上的變化趨勢(shì):隨著Y軸坐標(biāo)增加,渦量值先略微減小再快速增加,最后快速降到零值.粗糙冰也不會(huì)改變翼型尾緣產(chǎn)生的負(fù)值渦量在Y軸方向上的變化趨勢(shì),但與光滑翼型相比,粗糙冰會(huì)使渦量的絕對(duì)值減小.

圖11 不同粗糙冰流場(chǎng)截面的渦量對(duì)比Fig.11 Spanwise vorticity comparison of different rough ices

圖12為光滑翼型和帶粗糙冰翼型在4個(gè)攻角下的標(biāo)準(zhǔn)化雷諾應(yīng)力分布.從圖中可以看出,在攻角為0°、2° 及5° 條件下,粗糙冰對(duì)標(biāo)準(zhǔn)化雷諾應(yīng)力的分布無(wú)明顯影響.從7° 攻角開始,由于砂紙具有一定厚度且和翼型表面不是平滑過(guò)渡,率先產(chǎn)生了應(yīng)力.在攻角為12° 時(shí),氣流產(chǎn)生了大分離,砂紙厚度帶來(lái)的誤差消失,標(biāo)準(zhǔn)化雷諾應(yīng)力的分布也歸于近似.由此可知砂紙的厚度不影響小攻角(0°,2°,5°)和氣流產(chǎn)生分離的大攻角(12°)時(shí)標(biāo)準(zhǔn)化雷諾應(yīng)力的分布.

圖12 光滑翼型和R4型粗糙冰作用下6個(gè)攻角下的標(biāo)準(zhǔn)化雷諾應(yīng)力Fig.12 Normalized Reynolds stress of clean airfoil and R4 rough ice at six angles of attack

圖13為X/C×100=-120處截面的標(biāo)準(zhǔn)化雷諾應(yīng)力分布對(duì)比,可見該截面處標(biāo)準(zhǔn)化雷諾應(yīng)力的負(fù)值區(qū)域分布類似,正值區(qū)域內(nèi)在Y/C×100=0附近差別較大,粗糙冰的存在降低了正值區(qū)域的雷諾應(yīng)力.從圖9(d)中可見,標(biāo)準(zhǔn)化雷諾應(yīng)力分布有差異的區(qū)域在分離泡與外界層流的交界處.粗糙冰存在導(dǎo)致該處標(biāo)準(zhǔn)化雷諾應(yīng)力減小,氣流流動(dòng)相較光滑翼型更為平穩(wěn),從另一個(gè)角度解釋了粗糙冰在12° 攻角下會(huì)使分離泡的左極限相較光滑翼型更靠近翼型前緣.

圖13 X/C×100=-120截面的標(biāo)準(zhǔn)化雷諾應(yīng)力分布對(duì)比Fig.13 Normalized Reynolds stress comparison of section X/C×100=-120

4 結(jié)論

探索粗糙冰對(duì)翼型氣動(dòng)力的影響對(duì)飛機(jī)的防/除冰設(shè)計(jì)來(lái)說(shuō)具有非常重要的價(jià)值和意義.目前,對(duì)粗糙冰的研究以數(shù)值模擬為主,本文使用PIV技術(shù)在風(fēng)洞內(nèi)獲取翼型繞流流場(chǎng)信息,使用砂紙表征粗糙冰的粗糙度.通過(guò)改變雷諾數(shù)、粗糙冰的粗糙度和翼型攻角,分析了速度場(chǎng)、流線圖、渦量分布、標(biāo)準(zhǔn)化雷諾應(yīng)力分布及相關(guān)信息的進(jìn)一步提取.主要研究結(jié)論如下:

(1) 實(shí)驗(yàn)使用的雷諾數(shù)未對(duì)帶R1型粗糙冰的翼型繞流流場(chǎng)的流動(dòng)分離產(chǎn)生明顯影響.雷諾數(shù)的增加不影響翼型尾流的渦量在Y軸方向上的變化趨勢(shì),但會(huì)增大尾流的渦量數(shù)值.雷諾數(shù)影響了翼型尾流處的標(biāo)準(zhǔn)化雷諾應(yīng)力分布,雷諾數(shù)越大則分布越均勻,尾流的流動(dòng)越平穩(wěn).

(2) 帶有不同粗糙冰的翼型流場(chǎng)在相同雷諾數(shù)和攻角條件下,氣流流動(dòng)的速度最大值隨著粗糙度增大而減小,最大值的位置也更靠近翼型尾端,翼型的升力受到的影響就越大.從渦量云圖和雷諾應(yīng)力云圖來(lái)看,粗糙冰的存在使尾流流動(dòng)變得紊亂且具有旋度的流動(dòng)區(qū)域增大.

(3) 與光滑翼型相比,在12° 攻角下帶R4型粗糙冰的翼型繞流流場(chǎng)中氣流提前分離,且提前量為1%的翼型弦長(zhǎng).粗糙冰作用下分離泡的寬度更小,中心位置的流動(dòng)更紊亂.從渦量云圖和標(biāo)準(zhǔn)化雷諾應(yīng)力云圖來(lái)看,粗糙冰不影響12° 攻角下渦量和雷諾應(yīng)力在Y軸方向上的變化趨勢(shì),但降低了相應(yīng)的數(shù)值,解釋了粗糙冰作用下分離泡寬度更小的原因.