高超聲速邊界層轉(zhuǎn)捩的油膜干涉測量技術(shù)研究

董昊, 劉是成, 耿璽, 程克明

(南京航空航天大學(xué) 航空宇航學(xué)院, 江蘇 南京 210016)

邊界層轉(zhuǎn)捩一直是流體力學(xué)研究中的熱點(diǎn)問題，也是人們長期關(guān)注但尚未解決的研究領(lǐng)域之一。對高超聲速飛行器來說，邊界層轉(zhuǎn)捩直接關(guān)系到表面摩擦阻力的大小，摩擦阻力對表面熱流及飛行器的熱防護(hù)系統(tǒng)有重要影響。因此，高超聲速邊界層轉(zhuǎn)捩的預(yù)測方法顯得十分重要。到目前為止，已經(jīng)發(fā)展改進(jìn)了各種高超聲速的轉(zhuǎn)捩預(yù)測方法，數(shù)值模擬方面包括直接數(shù)值模擬[1]和大渦模擬[2]等。而這些方法受到計(jì)算機(jī)條件和算法的制約，難以廣泛應(yīng)用于飛行器的工程設(shè)計(jì)中[3]。風(fēng)洞實(shí)驗(yàn)是研究高超聲速邊界層轉(zhuǎn)捩的重要手段，但是受到了實(shí)驗(yàn)方法的制約。傳統(tǒng)的紋影難以測量到三維的精確結(jié)果，新興的粒子圖像測速(PIV)[4]、平面激光誘導(dǎo)熒光(PLIF)[5]、基于納米技術(shù)的平面激光散射(NPLS)[6]方法雖然測量精度更高，但均存在校準(zhǔn)方法和實(shí)驗(yàn)方法復(fù)雜等問題。

油膜干涉法是一種適用于高超聲速邊界層轉(zhuǎn)捩區(qū)域測量的實(shí)驗(yàn)方法。該方法具有測量面積大、測量精度高和實(shí)驗(yàn)成本低的特點(diǎn)。目前該方法在低速、亞聲速及超聲速風(fēng)洞實(shí)驗(yàn)中已經(jīng)有了較多的應(yīng)用，并且取得很好的效果[7-9]。但是，在高超聲速中的發(fā)展卻十分緩慢。國內(nèi)劉志勇、代成果等[10-11]將油膜干涉法用于高超聲速平板的摩阻因數(shù)測量中，取得了一定的效果，實(shí)驗(yàn)中選用了黏度較大的硅油，導(dǎo)致吹風(fēng)時(shí)間較長，實(shí)驗(yàn)?zāi)Ｐ捅砻鏈囟茸兓^大，硅油黏性隨溫度發(fā)生變化，會導(dǎo)致計(jì)算的摩阻因數(shù)結(jié)果產(chǎn)生一定的誤差。因此，將現(xiàn)有的油膜干涉測量技術(shù)運(yùn)用到高超聲速風(fēng)洞實(shí)驗(yàn)?zāi)ψ枰驍?shù)精確測量，尤其是邊界層轉(zhuǎn)捩區(qū)域位置的判斷中還存在著一些問題，需要對其進(jìn)行進(jìn)一步改進(jìn)，保證既能夠提高摩阻因數(shù)測量精度，又能快速準(zhǔn)確地預(yù)測高超聲速邊界層轉(zhuǎn)捩的位置，為高超聲速飛行器設(shè)計(jì)提供技術(shù)支撐。

本文在南京航空航天大學(xué)高超聲速風(fēng)洞(NHW)中進(jìn)行了油膜干涉法測量摩擦阻力因數(shù)的實(shí)驗(yàn)研究，并且首次將其應(yīng)用于三維粗糙元誘導(dǎo)高超聲速平板轉(zhuǎn)捩的研究中。首先對硅油黏度特性隨溫度的變化進(jìn)行了標(biāo)定，并利用低黏度的硅油代替先前研究中使用的高黏度硅油，大幅度減少了實(shí)驗(yàn)時(shí)間，解決了模型表面溫度變化較大對測量結(jié)果精確度的影響。結(jié)果表明油膜干涉法可以很好地應(yīng)用于高超聲速邊界層轉(zhuǎn)捩區(qū)域測量的研究中。

1 實(shí)驗(yàn)設(shè)備和方法

1.1 風(fēng)洞及測量設(shè)備

實(shí)驗(yàn)是在南京航空航天大學(xué)高超聲速風(fēng)洞(NHW)中進(jìn)行的，NHW是一座噴管出口直徑Φ500 mm的高壓下吹、真空抽吸暫沖式高超聲速風(fēng)洞，其風(fēng)洞示意圖和照片見圖1。該風(fēng)洞主要包括：氣源系統(tǒng)及高壓閥門、金屬板蓄熱式加熱器、熱閥、穩(wěn)定段和噴管、試驗(yàn)段、擴(kuò)壓器、真空系統(tǒng)、電氣控制系統(tǒng)和數(shù)據(jù)采集處理系統(tǒng)、電動攻角支撐機(jī)構(gòu)、Φ300 mm彩色紋影系統(tǒng)和計(jì)算機(jī)視頻數(shù)字?jǐn)z錄采集系統(tǒng)等。NHW具有5副軸對稱超聲速型面噴管，實(shí)驗(yàn)名義馬赫數(shù)4,5,6,7,8，真空球容積為650 m3，每車次實(shí)驗(yàn)時(shí)間7～10 s。試驗(yàn)段及真空駐室上方有一個(gè)350 mm×350 mm的觀察窗，通過觀察窗可以使用CCD相機(jī)或高速攝影對模型和流場進(jìn)行相關(guān)光學(xué)測量。

圖1 NHW風(fēng)洞示意圖及照片

本實(shí)驗(yàn)主要用到的其他裝置及測量設(shè)備包括單色光源,高速相機(jī)以及高純度硅油。其中,單色光源采用單色LED鈉燈,其頻率為50 Hz,波長為589 nm。相機(jī)采用PCO公司所生產(chǎn)的S4高速相機(jī),本次實(shí)驗(yàn)采用的相機(jī)幀率為100幀/s,遠(yuǎn)低于其極限拍攝速率(1 000幀/s),圖2為NHW風(fēng)洞駐室上方安裝的PCO-S4相機(jī)。硅油則依據(jù)實(shí)驗(yàn)條件選用阿拉丁公司所生產(chǎn)的黏度為50 cSt的高純度分析級硅油。

圖2 實(shí)驗(yàn)用PCO-S4高速相機(jī)

1.2 實(shí)驗(yàn)?zāi)Ｐ?/h3>

根據(jù)光學(xué)原理可知,如果試驗(yàn)所用光為單色鈉燈發(fā)出的黃光,那么模型表面折射率為2.0時(shí)產(chǎn)生干涉條紋的效果最佳[11]。對多種表面材料干涉效果的研究結(jié)果如表1所示。

表1 模型表面干涉效果比較

以上材料都有各自的優(yōu)缺點(diǎn),SF11玻璃折射率最高,但是不易于制作試驗(yàn)?zāi)Ｐ?Mylar膜干涉效果好,但容易受到溫度的影響,不適用于高超聲速風(fēng)洞的高溫環(huán)境;鍍鎳的效果好,但是有污染不環(huán)保;不銹鋼和鋁的干涉效果較差,一般不建議采用。本實(shí)驗(yàn)采用在拋光后的不銹鋼模型表面噴涂黑色高溫反光漆,實(shí)驗(yàn)結(jié)果表明,黑色高溫反光漆可以產(chǎn)生清晰的干涉條紋,獲得準(zhǔn)確的結(jié)果,且工藝較為簡單,可作為高超聲速油膜干涉法測量的模型表面處理方法。

由于硅油黏度和溫度相關(guān),需對模型表面溫度進(jìn)行實(shí)時(shí)測量。測溫方式采用在模型表面布置多個(gè)熱電偶來獲取模型表面的溫度。熱電偶是一種常用的溫度檢測元件,具有測量精度高、測量范圍廣、構(gòu)造簡單、使用方便等特點(diǎn),易于布置在模型的任何位置。本實(shí)驗(yàn)選用廉金屬熱電偶中適合本實(shí)驗(yàn)測量溫度范圍以及精度等級較高的T型熱電偶,其精度為±0.5℃。

試驗(yàn)?zāi)Ｐ蜑槠桨迥Ｐ?模型尺寸為長190 mm,寬150 mm,平板前緣為尖楔形,倒角為20°。三維粗糙元為方形粗糙元,距離模型前緣50 mm,中心位于平板模型沿流向的對稱面上,粗糙元邊長為3 mm,高度為1.3 mm。將測溫用的T型熱電偶布置在測溫孔內(nèi),測溫端與平板上表面齊平并保持光滑,測溫孔直徑為1 mm。帶三維方形粗糙元的模型與測溫孔分布示意圖以及模型在風(fēng)洞中的安裝如圖3所示。

圖3 帶三維方形粗糙元的模型與測溫孔分布示意圖及模型在風(fēng)洞中安裝圖

1.3 硅油黏度標(biāo)定

前期研究表明,模型表面溫度變化帶來的油膜黏性因數(shù)的誤差是油膜干涉法誤差的主要來源[12]。在低速和亞跨聲速風(fēng)洞中,模型表面溫度接近常溫且變化范圍較小,此時(shí)溫度變化帶來的硅油黏性的變化不是很明顯,可按照常溫時(shí)硅油的黏度來處理。但常規(guī)高超聲速風(fēng)洞需要對來流進(jìn)行加熱以避免氣流過度膨脹產(chǎn)生的冷凝現(xiàn)象,因此來流的總溫會很高。本次風(fēng)洞實(shí)驗(yàn)前首先對使用硅油進(jìn)行了不同溫度下的黏度標(biāo)定,采用高精度旋轉(zhuǎn)黏度計(jì)對硅油黏度進(jìn)行了不同溫度下的測量。通過加熱改變硅油的溫度,獲得不同溫度下相應(yīng)的黏度值,圖4是對實(shí)驗(yàn)選用的黏度為50 cSt阿拉丁分析級硅油的標(biāo)定結(jié)果。

圖4 硅油黏度標(biāo)定結(jié)果

1.4 油膜干涉方法及光路設(shè)計(jì)

油膜干涉法的基本原理是利用硅油在流場中受到流體的推動而發(fā)生平移,從而形成楔形油膜。在流場中,受摩擦力影響,流場各處有各自的流動速度,從而導(dǎo)致油膜的厚度和斜率也都不相同,它們之間存在一定的算術(shù)關(guān)系[13]。于是通過干涉法,就可以得到油膜的厚度變化率,再通過一定的計(jì)算就能得到物體的表面摩擦力。

油膜厚度與表面摩擦應(yīng)力之間的關(guān)系如下:

(1)

式中坐標(biāo)系定義如下:x為模型表面上一點(diǎn)沿來流方向的坐標(biāo),y方向?yàn)樵擖c(diǎn)處垂直表面的法向方向,z方向滿足右手定則。因此,τw,x對應(yīng)x方向上的表面摩擦應(yīng)力分量,τw,z代表z方向的分量。h為油膜厚度,為實(shí)驗(yàn)實(shí)際測量所得。μ為硅油的黏性因數(shù),油膜厚度與干涉條紋間距之間關(guān)系如下所示:

(2)

式中,λ為實(shí)驗(yàn)所用的單色光波長,φ代表?xiàng)l紋級數(shù)。n0為硅油折射率。θi為光線入射角,可以通過測量得到。因而,只要將(2)式帶入(1)式,并去掉z方向的摩擦應(yīng)力分量即可得到所求的表面摩阻因數(shù)大小,結(jié)果如下:

(3)

式中，θ為光線折射角,可以由(2)式中的入射角經(jīng)折射定律得出。N為所測量的條紋級數(shù),Δx為N級條紋的總寬度,由實(shí)驗(yàn)實(shí)際測量得到。q∞為來流動壓。油膜干涉原理圖如圖5所示。

圖5 油膜干涉法原理示意圖[14]

圖6為NHW油膜干涉實(shí)驗(yàn)光路設(shè)計(jì)圖。平板模型安裝在流場中心位置,由于高速相機(jī)不能暴露在風(fēng)洞駐室的真空環(huán)境中,故其通過駐室上方的觀察窗對油膜干涉條紋進(jìn)行拍攝。為了獲得模型表面更大的拍攝范圍,單色鈉燈光源并沒有直接照射在模型上,而是通過風(fēng)洞駐室上方大面積的反光板進(jìn)行反射后,入射模型表面。這樣做的另一個(gè)優(yōu)點(diǎn)是減弱了入射光的強(qiáng)度,柔和的單色光能使干涉條紋更加清晰。

圖6 NHW油膜干涉實(shí)驗(yàn)光路設(shè)計(jì)圖

2 實(shí)驗(yàn)結(jié)果與分析

本次實(shí)驗(yàn)名義馬赫數(shù)為5,來流馬赫數(shù)Ma=4.97,來流總壓P0=3.909×105Pa,總溫T0=516 K,單位雷諾數(shù)為Re=4.7×106m-1。在進(jìn)行模型表面尺寸標(biāo)定后進(jìn)行風(fēng)洞實(shí)驗(yàn),測量范圍為距離平板前緣約50 mm到130 mm的局部區(qū)域。取風(fēng)洞流場建立且總壓穩(wěn)定時(shí)作為油膜條紋讀取的起始點(diǎn),每次實(shí)驗(yàn)車次時(shí)間設(shè)定為7 s。為了比對驗(yàn)證實(shí)驗(yàn)結(jié)果,采用文獻(xiàn)[15]中的基于有限體積的高階緊致格式(CFVM)并行算法,采用隱式大渦模擬的方法對相同來流條件下的光滑平板模型及帶三維方形粗糙元的平板模型進(jìn)行數(shù)值模擬,湍流模型采用kω-sst模型,根據(jù)實(shí)驗(yàn)條件,定義湍流黏性比為0.01。

圖7為光滑平板模型在風(fēng)洞流場建立3 s和4 s時(shí)的油膜干涉圖像,可以看出,硅油在剪切力的作用下向下游延展,產(chǎn)生硅油薄膜,經(jīng)單色鈉光照射形成清晰的干涉條紋,隨著時(shí)間的發(fā)展,干涉條紋逐步變寬,但相同時(shí)刻的條紋間距沿流向逐漸變窄,說明光滑平板模型表面的摩阻因數(shù)沿流向逐漸減小。

圖7 平板模型油膜干涉圖像

圖8為風(fēng)洞流場建立3～6 s,熱電偶測量的平板模型各測點(diǎn)處的表面溫度,可以看出,相同時(shí)刻沿流向模型表面溫度變化不大,實(shí)驗(yàn)中各點(diǎn)處相隔1 s前后的溫度變化也不明顯,3 s和4 s時(shí)的各點(diǎn)溫度相差約0.1℃,但是隨著時(shí)間的推進(jìn),各點(diǎn)溫度是逐漸升高的,因此,在能分辨出條紋寬度的前提下,盡量選取間隔較短的初始、最終時(shí)刻來減小溫度升高造成硅油黏度變化帶來的誤差。本文選取3 s和4 s時(shí)的平均溫度來對硅油黏度進(jìn)行修正,即平板表面溫度為48.3℃,由圖4可知,此時(shí)常溫下50 cSt的硅油黏度已經(jīng)變?yōu)?3.2 cSt。因此,在用油膜干涉法進(jìn)行高超聲速摩阻測量時(shí),對表面溫度進(jìn)行實(shí)時(shí)測量是有效和必要的。

圖8 平板模型表面溫度

圖9為由干涉條紋和溫度修正過的硅油黏度計(jì)算出的光滑平板模型沿流向中心線上相應(yīng)位置處的表面摩擦阻力因數(shù)和數(shù)值模擬結(jié)果的比較。

圖9 光滑平板數(shù)值模擬與風(fēng)洞實(shí)驗(yàn)結(jié)果比較

可以看出,光滑平板模型表面中心線處的摩阻因數(shù)沿流向逐漸減小,這與圖7結(jié)果一致,并且風(fēng)洞實(shí)驗(yàn)結(jié)果與層流數(shù)值模擬結(jié)果相吻合,Cf數(shù)值約為0.001,表明當(dāng)前平板邊界層流動處于層流流態(tài),在可視范圍內(nèi)沒有發(fā)生自然轉(zhuǎn)捩,可進(jìn)一步開展三維粗糙元誘導(dǎo)轉(zhuǎn)捩的實(shí)驗(yàn)和計(jì)算。

圖10為帶三維方形粗糙元的平板模型在Ma=5來流條件下的油膜干涉條紋,對比圖7可知,粗糙元后干涉條紋較光滑平板模型寬的多,這是由于三維粗糙元在其后方流場中形成尾渦,從而引起平板邊界層中強(qiáng)烈的擾動和剪切作用,使層流邊界層轉(zhuǎn)捩為湍流。圖11為帶三維方形粗糙元的平板模型油膜干涉法計(jì)算出的沿流向中心線上表面摩擦阻力因數(shù)和數(shù)值模擬結(jié)果的比較。

圖10 帶方形粗糙元平板模型油膜干涉圖像

圖11 帶方形粗糙元的平板數(shù)值模擬與風(fēng)洞實(shí)驗(yàn)結(jié)果比較

由圖可知,風(fēng)洞實(shí)驗(yàn)結(jié)果與數(shù)值模擬結(jié)果相吻合,三維粗糙元引起了平板表面摩阻因數(shù)的增加,在x=137 mm時(shí)Cf達(dá)到了峰值,約為0.002 75,之后沿流向趨于穩(wěn)定,這也說明三維粗糙元誘導(dǎo)了平板邊界層的轉(zhuǎn)捩,轉(zhuǎn)捩區(qū)直至x=137 mm處,之后流動完成發(fā)展為湍流。

3 結(jié) 論

本文首次在高超聲速風(fēng)洞中開展了油膜干涉法對邊界層轉(zhuǎn)捩區(qū)域位置預(yù)測的實(shí)驗(yàn)研究，通過與數(shù)值模擬結(jié)果的比對，得出以下結(jié)論：

1) 采用熱電偶可以準(zhǔn)確有效地測量模型表面的溫度，從而修正因?yàn)闇囟鹊淖兓瘜?dǎo)致實(shí)驗(yàn)結(jié)果的誤差；

2) 對于模型表面的處理，可采用在拋光后的不銹鋼模型表面做黑色高溫反光漆的方式，既能夠清晰地得到干涉條紋，又簡單方便；

3) 使用低黏度的硅油大大縮短了風(fēng)洞吹風(fēng)時(shí)間，使模型表面溫度變化較小，從而減小了溫度對硅油黏度的影響，提高了測量精度；

4) 對三維粗糙元誘導(dǎo)的平板模型邊界層轉(zhuǎn)捩實(shí)驗(yàn)研究表明，油膜干涉測量技術(shù)可應(yīng)用于高超聲速邊界層轉(zhuǎn)捩的研究中，通過表面摩阻因數(shù)的變化可對高超聲速邊界層轉(zhuǎn)捩區(qū)域位置進(jìn)行準(zhǔn)確的預(yù)測。