田 濱，李華星，孟宣市，羅時(shí)鈞，劉 鋒

（1.西北工業(yè)大學(xué) 翼型葉柵國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710072；2.美國(guó)加州大學(xué)爾灣分校 機(jī)械與宇航工程系，爾灣 92697－3975，美國(guó)）

0 引 言

具有高機(jī)動(dòng)性能的戰(zhàn)斗機(jī)、導(dǎo)彈等飛行器，頭部通常設(shè)計(jì)成細(xì)長(zhǎng)的尖頭旋成體。在大迎角下，飛行器的這種細(xì)長(zhǎng)前體會(huì)產(chǎn)生一對(duì)分離渦。當(dāng)迎角增加到一定程度時(shí)，分離渦會(huì)變得非對(duì)稱(chēng)，同時(shí)產(chǎn)生很大的側(cè)向力和力矩［1－2］。這種由于分離渦的非對(duì)稱(chēng)而產(chǎn)生的側(cè)向力和力矩，其大小和方向是不可預(yù)估的，對(duì)飛行器的機(jī)動(dòng)性和穩(wěn)定性有很大的影響。由于這種大迎角非對(duì)稱(chēng)分離渦對(duì)尖頭的擾動(dòng)是非常敏感的［3－4］。因此，通過(guò)在尖頭很小的區(qū)域注入少部分能量就能實(shí)現(xiàn)對(duì)飛行器的主動(dòng)控制，從而使得飛行器的機(jī)動(dòng)性和穩(wěn)定性得到很大的提高。

近年來(lái)，研究者們應(yīng)用各種新型流動(dòng)控制技術(shù)對(duì)細(xì)長(zhǎng)體大迎角非對(duì)稱(chēng)分離渦的流動(dòng)控制進(jìn)行了研究。E.Hanff等［5］在細(xì)長(zhǎng)尖拱體頭部安裝了兩個(gè)向前吹氣的氣管，實(shí)現(xiàn)了對(duì)細(xì)長(zhǎng)尖拱體側(cè)向力的比例控制。顧蘊(yùn)松等［6］通過(guò)在尖拱體頭部安裝小擾動(dòng)片的方式實(shí)現(xiàn)了側(cè)向力的比例控制。

近年來(lái)，劉鋒、孟宣市等［7－8］在風(fēng)洞中利用等離子流動(dòng)控制雙側(cè)占空循環(huán)技術(shù)實(shí)現(xiàn)了對(duì)細(xì)長(zhǎng)圓錐前體側(cè)向力和力矩的近似比例控制。本文在此基礎(chǔ)上，利用單側(cè)脈沖等離子流動(dòng)控制技術(shù)，改變不同的占空比，研究其對(duì)細(xì)長(zhǎng)圓錐前體分離渦流場(chǎng)的控制特性。

1 模型及試驗(yàn)裝置

試驗(yàn)?zāi)Ｐ陀蓤A錐段、圓弧過(guò)渡段和圓柱整流段3部分組裝而成（見(jiàn)圖1）。圓錐段半頂角為10°，長(zhǎng)度L＝463.8mm，底面直徑D＝163.6mm。在圓錐段布置有九個(gè)測(cè)壓截面，沿軸線(xiàn)等距分布。其中第八截面為動(dòng)態(tài)測(cè)壓截面，周向布置24個(gè)測(cè)壓孔。其他截面為靜態(tài)測(cè)壓截面，每個(gè)截面周向均勻布置36個(gè)測(cè)壓孔。另外，圓錐尖端頭部開(kāi)始的150mm長(zhǎng)度部分為絕緣塑料加工制成，用于粘貼等離子體激勵(lì)器。

在圓錐尖端頭部，對(duì)稱(chēng)安裝一對(duì)單電極介質(zhì)阻擋放電等離子激勵(lì)器，分別位于方位角θ＝±120°處，θ＝0°的位置是在迎風(fēng)面的正中。電極材料為銅箔制，厚度0.03mm。電極沿圓錐母線(xiàn)的長(zhǎng)度為20mm，前緣距圓錐頂點(diǎn)9mm，中間為厚度0.4mm的Kapton絕緣層，水平間距1.5mm。裸露電極和掩埋電極的寬度分別為1mm和2mm（見(jiàn)圖2）。激勵(lì)器所誘導(dǎo)的氣流從θ＝±120°起，與圓錐相切，方向是逆著來(lái)流的。根據(jù)Coanda效應(yīng)（附壁效應(yīng)），誘導(dǎo)氣流將附著于圓錐表面。

圖1 試驗(yàn)?zāi)Ｐ虵ig.1 Test model

圖2 激勵(lì)器的安裝Fig.2 Arrangement of the plasma actuators

南京蘇曼電子有限公司生產(chǎn)的CTP－2000K介質(zhì)阻擋放電等離子體發(fā)生器和1臺(tái)PC－07數(shù)字脈沖調(diào)制器為等離子體激勵(lì)器提供脈沖激勵(lì)和占空循環(huán)調(diào)制。試驗(yàn)使用的是單側(cè)脈沖控制技術(shù)，因此只有左舷等離子激勵(lì)器工作，右舷激勵(lì)器始終是關(guān)閉的，一直處于不工作狀態(tài)。

試驗(yàn)在西北工業(yè)大學(xué)NF－3直流式低速風(fēng)洞進(jìn)行。風(fēng)洞試驗(yàn)段尺寸長(zhǎng)×寬×高＝8m×3m×1.6m，最大風(fēng)速為120m／s，氣流湍流度ε≤0.045%。試驗(yàn)風(fēng)速為5m／s，基于圓錐段底面直徑的雷諾數(shù)Re＝5×104，迎角α＝35°。

試驗(yàn)所用的壓力掃描閥包括靜態(tài)掃描閥和動(dòng)態(tài)掃描閥。靜態(tài)壓力掃描閥選取的是PSI公司產(chǎn)的9816型，采樣頻率100Hz，測(cè)量精度為0.05%，采樣時(shí)間為15s。動(dòng)態(tài)壓力掃描閥選用Kulite型XCS－093，單通道采樣速率100K／s，測(cè)量精度為0.1%。試驗(yàn)中動(dòng)態(tài)掃描閥的采樣頻率為5000Hz，采樣時(shí)間均為30s。

2 試驗(yàn)結(jié)果

2.1 靜態(tài)壓力分布及側(cè)力

試驗(yàn)探討的是左舷（面向來(lái)流，模型左側(cè)）激勵(lì)器單側(cè)脈沖放電下，對(duì)流場(chǎng)進(jìn)行控制的情況。因此，右舷激勵(lì)器始終是關(guān)閉的，一直處于不工作狀態(tài)。加載在左舷等離子激勵(lì)器上的輸出電壓Vp－p＝14.5kV，頻率F＝11.8kHz，脈沖調(diào)制頻率為50Hz。τp＝0.01時(shí)，輸入功率最小，基本為0。τp＝0.99時(shí)，輸入功率達(dá)到最大值14.8W。

圖3分別給出了在不同占空比下第5截面方位角0°～360°和60°～300°之間的壓力分布情況。從圖中可以看出，隨著占空比的變化，吸力峰逐漸從左舷高、右舷低的狀態(tài)轉(zhuǎn)變?yōu)橛蚁细摺⒆笙系偷臓顟B(tài)，左右舷壓力分布發(fā)生了轉(zhuǎn)換。這說(shuō)明通過(guò)單側(cè)脈沖放電，改變占空比同樣可以實(shí)現(xiàn)對(duì)非對(duì)稱(chēng)渦的有效控制。τp＝0.01時(shí)的壓力分布和激勵(lì)器關(guān)閉時(shí)的壓力分布基本是重合的。這是因?yàn)楫?dāng)τp＝0.01時(shí)激勵(lì)器在一個(gè)脈沖周期內(nèi)工作的時(shí)間非常短，注入的能量很小，幾乎可以忽略。因此，對(duì)流場(chǎng)的影響微乎其微。τp＝0.01時(shí)的輸入功率也說(shuō)明了這一點(diǎn)，當(dāng)τp＝0.01時(shí)，輸入電流很小，基本為0，輸入的功率也近似等于0。當(dāng)τp＝0.6時(shí)，壓力分布曲線(xiàn)基本處于對(duì)稱(chēng)位置，這時(shí)的當(dāng)?shù)貍?cè)向力CYd也近似為0。τp＝0.01和τp＝0.99時(shí)的壓力分布實(shí)現(xiàn)了兩個(gè)雙穩(wěn)態(tài)之間的轉(zhuǎn)換，但沒(méi)有完全的鏡像對(duì)稱(chēng)，左舷吸力峰要略高于右舷。

圖3 第五截面周向壓力分布隨占空比的變化Fig.3 Circumference pressure distribution variations with different duty cycle ratios at Section 5th

圖4給出了奇數(shù)測(cè)量截面當(dāng)?shù)貍?cè)向力系數(shù)CYd隨占空比的變化情況。圖5描述的是當(dāng)?shù)貍?cè)向力系數(shù)增量ΔCYd隨占空比的變化情況。當(dāng)?shù)貍?cè)向力系數(shù)增量ΔCYd是通過(guò)不同占空比的當(dāng)?shù)貍?cè)向力系數(shù)減去激勵(lì)器關(guān)閉時(shí)的當(dāng)?shù)貍?cè)向力系數(shù)得來(lái)的。從圖4和圖5中可以看出，CYd的變化呈現(xiàn)良好的線(xiàn)性度，并且越是靠近圓錐前端的截面，線(xiàn)性度越好，到了圓錐的后端截面線(xiàn)性度有所減弱。但整體看來(lái)，這次的控制效果仍然是很好的，達(dá)到了良好的線(xiàn)性度。

圖4 奇數(shù)截面當(dāng)?shù)貍?cè)向力系數(shù)隨占空比的變化Fig.4 Local side force coefficient variation of odd sections with different duty cycle ratios

圖5 奇數(shù)截面當(dāng)?shù)貍?cè)向力系數(shù)增量隨占空比的變化Fig.5 The increment of local side force coefficient variation of odd sections with different duty cycle ratios

圖6給出了圓錐段總的側(cè)向力系數(shù)和偏航力矩系數(shù)隨占空比的變化情況。隨著占空比的變化，總的側(cè)向力和偏航力矩基本上是線(xiàn)性變化的。這說(shuō)明通過(guò)單側(cè)脈沖放電，改變不同的占空比同樣可以控制圓錐段的側(cè)向力和偏航力矩，并且可以達(dá)到雙穩(wěn)態(tài)之間的任意一個(gè)數(shù)值，包括零值。與文獻(xiàn)［7－8］研究結(jié)果相比較，這次的線(xiàn)性度有了很大的提高，這主要與脈沖頻率的提高有關(guān)。激勵(lì)器制作精度的提高也是一個(gè)重要方面。

圖6 圓錐段側(cè)向力系數(shù)和偏航力矩系數(shù)隨占空比的變化Fig.6 The variations of overall side force and yawing moment coefficients on cone with different duty cycle ratios

2.2 動(dòng)態(tài)結(jié)果分析

圖7描述的是第八截面動(dòng)態(tài)壓力分布30s全時(shí)間平均結(jié)果隨占空比的變化情況。從圖中可以看出，隨著占空比的變化，吸力峰從左舷高，右舷低的狀態(tài)逐漸變?yōu)橛蚁细?，左舷低的狀態(tài)。這與圖3給出的靜態(tài)壓力分布相一致。圖8給出了第八截面動(dòng)態(tài)壓力分布在τp＝0.8時(shí)不同時(shí)間下的全時(shí)間平均結(jié)果，分別取1s、10s、15s、20s、25s、30s進(jìn)行全時(shí)間平均。從圖中可以看出，10s以后壓力分布基本達(dá)到收斂。

圖7 第八截面周向壓力分布隨占空比的變化Fig.7 Circumference pressure distribution variations with different duty cycle ratios at Section 8th

圖8 第八截面周向壓力分布隨時(shí)間的變化（τp＝0.8）Fig.8 Circumference pressure distribution variations with different time at Section 8th（τp＝0.8）

本次試驗(yàn)中動(dòng)態(tài)傳感器的采集頻率是5000Hz，采集時(shí)間為30s，等離子脈沖調(diào)制頻率是50Hz。因此，每一個(gè)脈沖調(diào)制周期有100個(gè)相位角。在30s內(nèi)，每個(gè)相位角采集了1500個(gè)數(shù)據(jù)。圖9描述的是第八截面τp＝0.4時(shí)在給定的相位角下隨不同時(shí)間相位鎖定平均的壓力分布，選取的相位角為Ψ＝0°、144°、288°。從圖9可以看出，各個(gè)相位角下的壓力分布曲線(xiàn)在30s都達(dá)到了很好的收斂，只有Ψ＝144°時(shí)160°～240°方位角之間的測(cè)量值在30s內(nèi)沒(méi)有收斂。這是因?yàn)檫@個(gè)區(qū)域正處于兩個(gè)分離渦的中間，對(duì)擾動(dòng)比其他區(qū)域更加敏感。其結(jié)果體現(xiàn)了本次試驗(yàn)單側(cè)脈沖激勵(lì)的非定常特性。因此，從流場(chǎng)穩(wěn)定性角度看，這次的收斂效果是很好的。

圖9 相位鎖定平均壓力分布隨時(shí)間的收斂情況（τp＝0.4）Fig.9 Convergence of phase－locked－averaged pressure distribution with different times（τp＝0.4）

圖10描述的是第八截面τp＝0.4時(shí)在不同相位角下30s相位鎖定平均的壓力分布曲線(xiàn)與30s全時(shí)間平均的壓力分布曲線(xiàn)，選取的相位角為Ψ＝0°、72°、144°、216°、288°。從圖9的分析已經(jīng)得出，各個(gè)相位角在30s基本都已收斂。圖10中除θ＝160°～240°的區(qū)域外，相位鎖定的壓力分布曲線(xiàn)與全時(shí)間平均的壓力分布曲線(xiàn)都趨于重合。而θ＝160°～240°區(qū)域壓力分布的跳動(dòng)正體現(xiàn)了單側(cè)脈沖激勵(lì)的非定常特性。當(dāng)τp＝0.4時(shí)，在一個(gè)脈沖調(diào)制周期里，激勵(lì)器開(kāi)啟的時(shí)間是40%，其中有60%的時(shí)間是不工作的。因此，若流場(chǎng)的響應(yīng)是即時(shí)的，那么在激勵(lì)器不工作的時(shí)間里，壓力分布應(yīng)立即回到激勵(lì)器關(guān)閉的狀態(tài)，左右舷吸力峰應(yīng)該發(fā)生變化。比較不同相位角下的壓力分布，可以看到曲線(xiàn)并沒(méi)有發(fā)生明顯變化。這與文獻(xiàn)［9－10］的結(jié)果有所不同。這說(shuō)明，在脈沖調(diào)制頻率50Hz的情況下，流動(dòng)存在滯后現(xiàn)象［11］。第八截面處流場(chǎng)的響應(yīng)滯后于此次脈沖放電的頻率。

圖10 相位鎖定和全時(shí)間平均的壓力分布（τp＝0.4）Fig.10 Phase－locked and ensemble average pressure distribution（τp＝0.4）

3 結(jié) 論

試驗(yàn)對(duì)圓錐－圓柱組合體進(jìn)行了主動(dòng)流動(dòng)控制研究，利用單側(cè)脈沖等離子流動(dòng)控制技術(shù)實(shí)現(xiàn)了對(duì)圓錐－圓柱組合體的有效控制。試驗(yàn)證明了：

（1）當(dāng)非對(duì)稱(chēng)渦流場(chǎng)處于雙穩(wěn)態(tài)時(shí)，通過(guò)單側(cè)脈沖放電同樣可以實(shí)現(xiàn)對(duì)細(xì)長(zhǎng)圓錐前體非對(duì)稱(chēng)側(cè)向力和力矩的比例控制，且具有良好的線(xiàn)性度。

（2）第八截面動(dòng)態(tài)壓力數(shù)據(jù)的全時(shí)間平均和相位鎖定平均結(jié)果表明，在全時(shí)間平均和相位鎖定平均下，壓力分布都達(dá)到收斂。

（3）在脈沖調(diào)制頻率50Hz的情況下，流場(chǎng)的響應(yīng)滯后于此次脈沖放電的頻率，存在滯后現(xiàn)象。

致謝：本文得到了西北工業(yè)大學(xué)NF－3風(fēng)洞高永衛(wèi)教授、惠增宏高級(jí)工程師、肖春生工程師和課題組王健磊、趙帥、馬沖、秦浩的幫助，在此向他們表示感謝。

［1］ERICSSON L.Sources of high alpha vortex asymmetry at zero sideslip［J］.JournalofAircraft，1992，29（6）：1086－1090.

［2］LOWSON M，PONTON A.Symmetry breaking in vortex flows on conical bodies［J］.AIAAJ.，1992，30：1576－1583.

［3］ZILLIAC G G，DEGANI D，TOBAK M.Asymmetric vortices on a slender body of revolution［J］.AIAAJ.，1991，29（5）：667－675.

［4］陳學(xué)銳，鄧學(xué)鎣.旋成體頭部擾動(dòng)對(duì)非對(duì)稱(chēng)背渦特性的影響［A］.全國(guó)第九屆分離流、旋渦和流動(dòng)控制會(huì)議論文集［C］，北京，2002：253－262.

［5］HANFF E，LEE R，KIND R J.Investigations on a dynamic forebody flow control system［C］.Proceedings of the 18th International Congress on Insrumentation in Aerospace Simulat ion Facilities.1999，28：1－9.

［6］顧蘊(yùn)松，明曉.大迎角細(xì)長(zhǎng)體側(cè)向力的比例控制［J］.航空學(xué)報(bào)，2006，27（5）：746－750.

［7］LIU F，LUO S J，GAO C，et al.Flow control over a conical forebody using duty－cycled plasma actuators［J］.AIAAJ.，2008，46（11）：2969－2973.

［8］孟宣市，郭志鑫，劉鋒，等.細(xì)長(zhǎng)圓錐前體非對(duì)稱(chēng)渦流場(chǎng)的等離子體控制［J］.航空學(xué)報(bào)，2010，31（3）：500－505.

［9］MENG X S，GUO Z X，LIU F，et al.Ensemble and phase－locked averaged loads controlled by plasma duty cycles［R］.AIAA Flow Control Conference.Orlando，F(xiàn)lorida，010.AIAA 2010－878.

［10］LIU F，LUO S J，GAO C，et al.Mechanisms for conical forebody flow control using plasma actuators［R］.AIAA 2009－4284，2009.

［11］李應(yīng)紅，梁華，馬清源，等.脈沖等離子體氣動(dòng)激勵(lì)抑制翼型吸力面流動(dòng)分離的試驗(yàn)［J］.航空學(xué)報(bào)，2008，29（6）：1429－1435.