雷玉昌，張登成，張艷華，蘇光旭

（1.空軍工程大學(xué) 研究生院，西安 710051； 2.空軍工程大學(xué) 航空工程學(xué)院，西安 710038）

飛行器上的襟翼、副翼等活動(dòng)舵面一般被用來改變機(jī)翼的彎度和周圍的流場特性，進(jìn)而改變機(jī)翼上的氣動(dòng)力和力矩。活動(dòng)舵面對(duì)提高飛行器的升力、改善升阻比起到了重要作用，但是復(fù)雜的舵面機(jī)械系統(tǒng)導(dǎo)致飛行器結(jié)構(gòu)復(fù)雜、質(zhì)量增加，隱身性能變差。為解決復(fù)雜的高升力舵面機(jī)械系統(tǒng)存在的弊端，國內(nèi)外陸續(xù)使用環(huán)量控制技術(shù)來解決這個(gè)矛盾［1］。環(huán)量控制技術(shù)的發(fā)展是備受關(guān)注的主動(dòng)流動(dòng)控制方式之一。這種流動(dòng)控制是指在機(jī)翼鈍后緣表面開縫，在后緣上噴射切向射流，利用Coanda效應(yīng)（壓力與離心力的平衡使射流沿著后緣曲面流動(dòng)），射流與外流混合后繞曲面外形流動(dòng)，使流線偏折，并增大繞流速度，增加氣流繞翼型的環(huán)量。由庫塔儒科夫升力定理可知翼型升力會(huì)增加［2-3］。風(fēng)洞試驗(yàn)和數(shù)值仿真的結(jié)果表明，環(huán)量控制技術(shù)對(duì)圓形Coanda后緣穩(wěn)定切向吹氣具有明顯的氣動(dòng)效益，對(duì)縮短飛行器的起飛著陸距離、提高機(jī)動(dòng)性等提供了可能性［3-4］。盡管有諸多好處，但是這項(xiàng)技術(shù)并沒有應(yīng)用于實(shí)際飛行器上，因?yàn)楹缶壣淞餍枰獜陌l(fā)動(dòng)機(jī)引氣或設(shè)置額外供氣裝置，由此帶來的發(fā)動(dòng)機(jī)推力損失及經(jīng)濟(jì)效益制約了該技術(shù)的實(shí)際應(yīng)用。而脈沖射流能夠大幅度降低環(huán)量控制射流所需的質(zhì)量流量，并減少引氣量［5］?？紤]到節(jié)省1%的商用噴氣燃料就相當(dāng)于節(jié)省100萬美元的運(yùn)營成本，脈沖射流能夠帶來可觀的經(jīng)濟(jì)效益［6］。為此，脈沖射流對(duì)環(huán)量控制翼型氣動(dòng)性能的影響規(guī)律，在著重分析其大迎角下的實(shí)際應(yīng)用價(jià)值中有重要的研究價(jià)值。

國內(nèi)外對(duì)環(huán)量控制射流都進(jìn)行了一定程度的研究。Yaros等［7］證明了環(huán)量控制可使機(jī)翼最大升力系數(shù)增大4倍。Jones等［8］也證明了采用環(huán)量控制的機(jī)翼可以獲得與傳統(tǒng)高升力系統(tǒng)相當(dāng)甚至更高的升力，但隨著升力增大，所需要的射流質(zhì)量流量也迅速增大。Kanistras等［9］通過實(shí)驗(yàn)的方法研究了不同后緣曲率對(duì)氣動(dòng)特性的影響，結(jié)果表明后緣曲率對(duì)升力系數(shù)的影響較大，對(duì)阻力系數(shù)的影響并不明顯。Jones等［10］提出了升力增長存在分離區(qū)和超環(huán)量控制區(qū)2個(gè)控制階段，這取決于后緣形狀和射流特性。AVT-239任務(wù)組［11］設(shè)計(jì)制作了利用環(huán)量控制技術(shù)替代傳統(tǒng)舵面的無人飛行器，驗(yàn)證了該技術(shù)可用于改善升阻特性。朱自強(qiáng)和吳宗成［12］總結(jié)了之前環(huán)量控制技術(shù)研究的最新進(jìn)展，并指出了國內(nèi)對(duì)環(huán)量控制技術(shù)研究經(jīng)驗(yàn)和成果的缺乏。張艷華等［13-14］研究了環(huán)量控制的作用機(jī)理，特別是對(duì)等離子體環(huán)量控制進(jìn)行了深入研究，相較于普通射流環(huán)量的控制效費(fèi)比較大。南京航空航天大學(xué)、西北工業(yè)大學(xué)等也進(jìn)行了相關(guān)研究［15-17］。上述研究都集中指出了定常射流只能在較低迎角下帶來巨大的升力效益，高動(dòng)量系數(shù)導(dǎo)致失速迎角提前，并在大迎角下氣動(dòng)效益變差；同時(shí)，高動(dòng)量系數(shù)需要消耗大量的質(zhì)量流量，影響飛行器的推力效益。這些問題都嚴(yán)重制約了環(huán)量控制技術(shù)的進(jìn)一步應(yīng)用。

Jones和Englar［3，18］用 脈 沖 射 流 替 代 定 常 射流用以環(huán)量控制吹氣，證明了脈沖射流能夠在取得與定常射流相當(dāng)?shù)纳η闆r下，降低約50%的質(zhì)量流量，初步證明了脈沖射流的部分優(yōu)越性，但是并未解釋射流作用機(jī)理，也未解決脈沖射流帶來的氣動(dòng)力脈動(dòng)問題。隨后，Shah等［5］進(jìn)一步分析了脈沖射流對(duì)后緣不同位置的速度影響，以及后緣射流分離點(diǎn)和邊界層厚度的變化。王萬波和Warsop等［19-20］也進(jìn)行了相關(guān)研究，分析了0°迎角下不同頻率時(shí)脈沖吹氣對(duì)翼型氣動(dòng)性能的影響，同時(shí)指出仍需進(jìn)行動(dòng)量系數(shù)、占空比等因素的影響研究。已有的研究多關(guān)注脈沖射流作用下的時(shí)均氣動(dòng)力效應(yīng)，對(duì)相應(yīng)的氣動(dòng)力脈動(dòng)現(xiàn)象和大迎角下的流動(dòng)情況研究較少，難以指導(dǎo)實(shí)際應(yīng)用。

為此，本文采用雷諾平均Navier-Stokes方法對(duì)二維環(huán)量控制翼型流場進(jìn)行非定常數(shù)值模擬。對(duì)比研究脈沖射流和定常射流參數(shù)對(duì)翼型氣動(dòng)特性的影響，對(duì)大迎角下的流動(dòng)情況進(jìn)行分析，探究減弱氣動(dòng)力脈動(dòng)現(xiàn)象的方法，為脈沖射流的具體應(yīng)用提供一定的參考。

1 數(shù)值計(jì)算方法與驗(yàn)證

1.1 計(jì)算模型和網(wǎng)格

超臨界翼型有鈍前緣、大厚弦比的特點(diǎn)，最大升力系數(shù)較大，內(nèi)部能夠容納環(huán)量控制技術(shù)所需要的供氣機(jī)構(gòu)或管道等，可以更充分地發(fā)揮環(huán)量控制增大升力的優(yōu)勢，因而成為環(huán)量控制技術(shù)研究的良好翼型。本文采用相對(duì)厚度為17%的超臨界翼型NASA LS（1）-0417作為基準(zhǔn)翼型。

翼型相對(duì)厚度為17%，弦長c=240 mm。采用文獻(xiàn)［18］中的修形方法，將距后緣25%弦長的部分修形處理，按照后緣半徑r/c=2%繪制柯恩達(dá)后緣曲面，射流口偏離中軸線15°，射流口高度h=0.001c。圖1為修形過的超臨界翼型示意圖，研究僅涉及上射流口吹氣。

圖1 翼型相關(guān)參數(shù)示意圖Fig.1 Sechematic map of airfoil shape parameters

計(jì)算區(qū)域選取翼型弦長的30倍，網(wǎng)格剖分采用結(jié)構(gòu)網(wǎng)格，生成O型網(wǎng)格拓?fù)浣Y(jié)構(gòu)，邊界層內(nèi)第一層的網(wǎng)格高度約為1×10-5m，射流口附近第一層網(wǎng)格高度為2×10-6m，保證第一層高度的壁面率y+均小于1，以滿足黏性底層的計(jì)算要求。網(wǎng)格總數(shù)約為58萬，圖2為翼型計(jì)算網(wǎng)格劃分情況。

圖2 翼型計(jì)算網(wǎng)格劃分Fig.2 Airfoil’s computational grid generation

1.2 計(jì)算方法

數(shù)值模擬方法采用二維雷諾平均Navier-Stokes方程，積分表達(dá)式為

式中：Q為流動(dòng)變量矩陣；Ω為控制體；S為控制體面積；FC和FV分別為無黏性通量和黏性通量；n為控制體面的外法線向量。采用k-w SST湍流模型，該模型對(duì)有較大逆壓梯度的邊界層流動(dòng)、分離預(yù)測性能較好。Swanson等［21-22］也曾指出該模型能較好地預(yù)測環(huán)量控制流動(dòng)分離。借助商業(yè)軟件FLUENT進(jìn)行相關(guān)計(jì)算，采用有限體積法離散控制方程，無黏性通量采用三階MUSCL格式，黏性通量采用二階迎風(fēng)格式離散。遠(yuǎn)場邊界為壓力遠(yuǎn)場，壁面邊界采用無滑移壁面條件。

早期研究表明，方波脈沖射流比傳統(tǒng)的正弦脈沖射流效率更高［23］。因此，本節(jié)采用方波脈沖射流進(jìn)行相關(guān)研究，其中脈沖射流在一個(gè)周期內(nèi)對(duì)應(yīng)的射流速度為

式中：T為脈沖周期；topen為一個(gè)脈沖周期中射流口產(chǎn)生射流的持續(xù)時(shí)間；DC為占空比；Vjet.max為在一個(gè)周期中射流產(chǎn)生的最大射流速度。射流動(dòng)量系數(shù)為衡量射流強(qiáng)度、表征氣流動(dòng)量的重要無量綱參數(shù)，定義為

非定常計(jì)算中，物理時(shí)間步長對(duì)數(shù)值計(jì)算起著重要的作用。特別是對(duì)脈沖射流而言，出口速度隨時(shí)間變化呈現(xiàn)出周期性變化，時(shí)間步長必須足夠小以解析脈沖射流的周期性變化。雖然時(shí)間步長過小會(huì)嚴(yán)重影響計(jì)算效率，但是對(duì)結(jié)果沒有太大影響。圖3為不同時(shí)間步長Δt下時(shí)均升力系數(shù)CL和阻力系數(shù)CD的計(jì)算結(jié)果。當(dāng)Δt≤5×10-5s時(shí)，時(shí)均升力系數(shù)CL的最大變化為0.9%，時(shí)均阻力系數(shù)CD的最大變化為2.1%，對(duì)高升力系統(tǒng)而言，存有一定計(jì)算誤差是可以被接受的。因此，選定計(jì)算時(shí)間步長Δt=5×10-5s，足以解析每個(gè)脈沖射流的周期性變化。

圖3 不同時(shí)間步長下的時(shí)均升阻力系數(shù)Fig.3 Time-averaged lift and drag coefficient at different time step sizes

1.3 計(jì)算方法驗(yàn)證

為了驗(yàn)證本文中數(shù)值計(jì)算方法的準(zhǔn)確性，對(duì)0°迎角下不同吹氣動(dòng)量系數(shù)的超臨界翼型進(jìn)行了數(shù)值模擬。計(jì)算條件為：來流速度為30 m/s，來流溫度為293.15 K，雷諾數(shù)為5×105，后續(xù)研究均采用此外流條件。設(shè)置動(dòng)量系數(shù)在0～0.06之間變化，將仿真結(jié)果與文獻(xiàn)［18］中的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，仿真采用穩(wěn)態(tài)方法，結(jié)果如圖4所示，仿真結(jié)果與實(shí)驗(yàn)數(shù)據(jù)基本吻合。

圖4 升力系數(shù)仿真與實(shí)驗(yàn)對(duì)比Fig.4 Comparison of lift coefficient between simulation and experiment

已有的環(huán)量控制脈沖射流實(shí)驗(yàn)中無高頻率內(nèi)容，因此選擇文獻(xiàn)［3］中35 Hz脈沖射流頻率下不同占空比的情況進(jìn)行對(duì)比。圖5給出了35 Hz脈沖射流頻率下的時(shí)均壓力系數(shù)Cp分布。從圖中可知，數(shù)值模擬得到的翼型上下表面的時(shí)均壓力系數(shù)Cp與實(shí)驗(yàn)結(jié)果基本吻合。

圖5 35 Hz脈沖射流下時(shí)均壓力系數(shù)分布Fig.5 Time-averaged pressure coefficient distribution under 35 Hz pulsed jet

2 結(jié)果分析

2.1 占空比的影響

占空比是脈沖射流對(duì)翼型氣動(dòng)特性的重要影響參數(shù)，通過調(diào)節(jié)占空比，可以獲得相應(yīng)質(zhì)量流量的變化。為了研究占空比對(duì)氣動(dòng)特性的影響，本節(jié)對(duì)35 Hz脈沖射流頻率下不同占空比情況的翼型氣動(dòng)參數(shù)進(jìn)行了數(shù)值模擬。

圖6為脈沖射流和定常射流下的時(shí)均升力系數(shù)隨質(zhì)量流量變化的對(duì)比曲線。其中，脈沖射流曲線對(duì)應(yīng)的質(zhì)量流量是在定常射流動(dòng)量系數(shù)下通過改變占空比實(shí)現(xiàn)的；定常射流動(dòng)量系數(shù)Cμ=0.02，不同占空比DC對(duì)應(yīng)的射流速度如表1所示。定常射流曲線對(duì)應(yīng)的質(zhì)量流量是通過改變射流速度實(shí)現(xiàn)的，動(dòng)量系數(shù)Cμ從0變化至0.02。

圖6 不同質(zhì)量流量下的時(shí)均升力系數(shù)對(duì)比Fig.6 Comparison of lift coefficient at different mass flow rates

表1 不同占空比對(duì)應(yīng)的射流速度Table 1 Jet velocity corresponding to different duty cycle

由圖6可以看出，在產(chǎn)生相同時(shí)均升力系數(shù)的基礎(chǔ)上，脈沖射流能夠大幅度降低射流所需要的質(zhì)量流量。在占空比變化階段，脈沖射流作用下的氣動(dòng)性能差別較大，較低占空比作用下，脈沖射流降低質(zhì)量流量的優(yōu)勢更加明顯。在給定時(shí)均升力系數(shù)為1的情況下，脈沖射流對(duì)應(yīng)的質(zhì)量流量減少了68.4%，同時(shí)，在占空比為40%、質(zhì)量流量相等的情況下，脈沖射流對(duì)應(yīng)的時(shí)均升力系數(shù)相較定常射流提高了65.3%。隨著占空比的增加，脈沖射流下的質(zhì)量流量效益逐漸降低。

圖7為不同占空比作用下的時(shí)均升力系數(shù)時(shí)域變化曲線。定義Δh為升力系數(shù)周期性變化過程中2個(gè)峰值之間的差值，即升力系數(shù)脈動(dòng)幅值。占空比DC越小，升力系數(shù)脈動(dòng)幅值Δh越大，這是因?yàn)榈驼伎毡葘?duì)應(yīng)的射流速度較大，是一種持續(xù)時(shí)間短、射流動(dòng)量大的穿透狀態(tài)，能在短時(shí)間內(nèi)在后緣移動(dòng)較大距離，而射流一旦消失，升力系數(shù)出現(xiàn)斷崖式下降的現(xiàn)象，射流消失后的升力響應(yīng)速度要遠(yuǎn)大于射流產(chǎn)生后的響應(yīng)速度，同時(shí)，由于占空比較低，射流消失后逐漸逼近無射流狀態(tài)下的升力系數(shù)。

圖7 35 Hz脈沖射流下時(shí)均升力系數(shù)時(shí)域曲線對(duì)比Fig.7 Comparison of lift coefficient time-domain curves at 35 Hz pulsed jet

綜上所述，在產(chǎn)生相同時(shí)均升力系數(shù)的前提下，脈沖射流需要的質(zhì)量流量更小，且占空比越小，這種優(yōu)勢越明顯。但是，低占空比下導(dǎo)致的射流穿透效應(yīng)將會(huì)直接導(dǎo)致升力系數(shù)脈動(dòng)幅值增大。當(dāng)DC=30%時(shí)，Δh可達(dá)1.3，這在實(shí)際飛行中將會(huì)導(dǎo)致劇烈的抖振現(xiàn)象，嚴(yán)重影響飛行安全，在實(shí)際應(yīng)用中是不可取的。因此，實(shí)際應(yīng)用中應(yīng)當(dāng)選用較高占空比的脈沖射流，選取DC=70%進(jìn)行后續(xù)研究，并進(jìn)一步探索脈沖射流的實(shí)際應(yīng)用前景及降低Δh的方法。

2.2 頻率的影響

圖8給出了在不同迎角α下頻率的改變對(duì)升力系數(shù)增量ΔCL的影響。低迎角（0°）下，隨頻率增加，升力系數(shù)增量呈現(xiàn)先增加后降低的趨勢，但整體變化幅度不大，在高頻率（200 Hz）下，脈沖射流甚至低于定常射流作用下的升力系數(shù)增量；高迎角（15°）下，隨頻率增加，升力系數(shù)增量呈現(xiàn)持續(xù)增加的趨勢，在高頻率（200 Hz）下，脈沖射流遠(yuǎn)遠(yuǎn)高于定常射流作用下的升力系數(shù)增量。在不同迎角下，隨頻率增加，脈沖幅值Δh均逐漸減小，約在100 Hz逐漸穩(wěn)定。綜上所述，高頻脈沖射流在較大迎角下有突出的應(yīng)用優(yōu)勢，具備替代定常射流進(jìn)行流動(dòng)控制的潛力。

圖8 不同頻率下的升力系數(shù)增量和脈動(dòng)幅值對(duì)比Fig.8 Comparison of lift coefficient and pulsation amplitude at different frequencies

為了分析脈沖射流作用下的流場作用機(jī)理，圖9給出了迎角分別為0°和15°時(shí)200 Hz脈沖射流與對(duì)應(yīng)定常射流作用下的翼型整體渦量圖和后緣局部流線放大圖。其中脈沖射流作用下各時(shí)刻對(duì)應(yīng)的時(shí)均升力系數(shù)如圖10所示，低頻率下射流消失后時(shí)均升力系數(shù)不再呈現(xiàn)出斷崖式下降的現(xiàn)象。

圖9 翼型渦量圖和流線圖Fig.9 Airfoil vorticity and streamlines

以15°迎角為例進(jìn)行分析，從定常射流作用下的翼型流線圖中可以看出，翼型上表面已經(jīng)出現(xiàn)大規(guī)模分離流動(dòng)，渦流區(qū)較大，繼續(xù)增大迎角會(huì)導(dǎo)致升力下降和阻力增加，翼型進(jìn)入失速分離階段。而脈沖射流作用下的翼型上表面分離渦與尾渦在一個(gè)脈沖周期內(nèi)逐漸相互耦合，共同向后移動(dòng)形成尾跡，上表面分離點(diǎn)明顯后移，渦流區(qū)明顯小于定常射流。而0°迎角下的翼型流線圖，定常射流作用下的流動(dòng)始終附體，并未出現(xiàn)分離流動(dòng)，脈沖射流作用下的上表面渦量分布與定常射流作用下的渦量分布差別也并不明顯。

從脈沖射流作用下的翼型流線圖中可以看出，在T/6時(shí)刻，翼型上表面存在分離流動(dòng)，并存在一個(gè)順時(shí)針的分離渦，由于后緣射流的產(chǎn)生，后緣上方存在一個(gè)逆時(shí)針尾渦，同時(shí)，后緣射流下存在一組較小的附著渦，翼型處于臨界失速狀態(tài)。當(dāng)進(jìn)入2T/6時(shí)刻，隨著后緣持續(xù)射流，后緣壓力逐漸減小，后緣上方的逆時(shí)針尾渦先向后移動(dòng)，此時(shí)上表面分離渦變化不大。當(dāng)進(jìn)入3T/6時(shí)刻，后緣壓力繼續(xù)減小，此時(shí)翼型上表面的分離渦也逐漸后移，分離渦幾乎已經(jīng)移動(dòng)至翼型尾端，翼型已經(jīng)遠(yuǎn)離失速狀態(tài)，后緣渦進(jìn)一步后移會(huì)導(dǎo)致升力系數(shù)下降，如圖10中2T/4時(shí)刻所示。

圖10 200 Hz脈沖射流下時(shí)均升力系數(shù)時(shí)域曲線對(duì)比Fig.10 Comparison of lift coefficient time-domain curves at 200 Hz pulsed jet

當(dāng)進(jìn)入4T/6時(shí)刻，原尾渦向后移動(dòng)至消散，上表面分離渦逐漸呈現(xiàn)脫離狀態(tài)，逐漸耦合形成新的尾渦。當(dāng)進(jìn)入5T/6時(shí)刻，射流已消失，此時(shí)分離渦已完全脫離，與原本位于射流下方的附著渦相互耦合，向后移動(dòng)。尾渦逐漸被外流吹散，上表面流動(dòng)逐漸恢復(fù)為附體。此時(shí)由于射流存在時(shí)滯效應(yīng)，后緣壓力并未迅速提高。當(dāng)進(jìn)入6T/6時(shí)刻，后緣無射流，此時(shí)流動(dòng)狀態(tài)趨向于15°迎角下的無射流定常狀態(tài)，鈍后緣形成新的一組附著渦，翼型上表面壓力較低，附著渦逐漸上移形成上表面的大分離區(qū)。隨后，后緣射流產(chǎn)生，射流上方產(chǎn)生新的尾渦，下方產(chǎn)生新的附著渦，流動(dòng)開始向下一個(gè)周期發(fā)展。

總體來看，翼型上表面的分離渦和后緣尾渦在一個(gè)周期的演化過程中，經(jīng)歷了形成、后移、脫落、消失、再形成的過程。而高頻脈沖射流在大迎角下的氣動(dòng)效益如此突出，是因?yàn)槠洳粩啻龠M(jìn)翼型上表面的分離渦脫落。

與15°迎角相比，0°迎角定常射流作用的翼型上表面并未分離，流動(dòng)始終附體。從整體流場來看，15°迎角上表面的分離渦不斷被吹除，渦流區(qū)明顯減小，較定常射流而言，升力增大；而0°迎角上表面不存在分離渦，尾跡受射流速度的影響跨度明顯較長，較定常射流而言升力略有降低。

從后緣流場來看，在整個(gè)周期并未形成如同15°迎角下明顯的分離渦與尾渦相互耦合向后移動(dòng)的狀態(tài)，更多表現(xiàn)為流線的上下偏折。這是因?yàn)?°迎角上表面渦強(qiáng)較弱，后緣難以耦合形成有效的渦旋。0°迎角下被不斷吹除的尾渦強(qiáng)度較低，連續(xù)性也較弱，因此，0°迎角下脈沖射流受頻率影響較低。

因此，對(duì)于低迎角而言，脈沖頻率更多表現(xiàn)為對(duì)升力系數(shù)的直接影響，與定常射流的作用機(jī)理是一致的，其效果類似于間歇性施加射流；而對(duì)于高迎角而言，脈沖頻率更多與分離渦移動(dòng)、尾緣渦脫落頻率等相關(guān)，這種流動(dòng)狀態(tài)已經(jīng)與定常射流完全不同。特別是在整個(gè)流動(dòng)周期中，分離渦的移動(dòng)都遠(yuǎn)離失速狀態(tài)，這一特性表明借助脈沖射流有望突破靜態(tài)失速迎角這一禁區(qū)。

2.3 動(dòng)量系數(shù)的影響

圖11為不同動(dòng)量系數(shù)下脈沖射流、定常射流、無射流及基準(zhǔn)翼型NASA LS（1）-0417的升力系數(shù)對(duì)比曲線?；谏鲜鰧?duì)占空比的研究，選取脈沖射流占空比為70%，動(dòng)量系數(shù)分別為0.02和0.04。從圖11（a）可以看出，在動(dòng)量系數(shù)為0.02的情況下，當(dāng)迎角小于10°時(shí)，定常射流能夠較大幅度地提高翼型升力系數(shù)，當(dāng)迎角大于10°時(shí)，提升幅度逐漸降低，而且失速迎角提前。為了衡量動(dòng)量系數(shù)帶來的影響，引用效費(fèi)比ΔCL/Cμ參數(shù)。0°迎角下，翼型效費(fèi)比為42.6，此時(shí)翼型升力系數(shù)提高144.5%；而15°迎角下，翼型效費(fèi)比降為15.2，此時(shí)翼型升力系數(shù)僅提高18.8%，同時(shí)，失速迎角由15°提前到10°。當(dāng)迎角小于5°時(shí)，定常射流升力系數(shù)大于脈沖射流；當(dāng)迎角大于5°時(shí)，定常射流升力系數(shù)小于脈沖射流。相較定常射流，脈沖射流能在更大的迎角范圍內(nèi)大幅度提高翼型升力系數(shù)，更加適合在大迎角下工作。在15°迎角時(shí)，ΔCL/Cμ=38.3，此時(shí)翼型升力系數(shù)仍能提高46.3%，同時(shí)，失速迎角大幅度延遲。

圖11 升力系數(shù)對(duì)比Fig.11 Comparison of lift coefficient

從圖11（b）可以看出，提高動(dòng)量系數(shù)至0.04后，對(duì)于定常射流而言，在0°迎角下，ΔCL/Cμ=36.9，而在15°迎角下，翼型效費(fèi)比降為11.9。隨著動(dòng)量系數(shù)的增加，定常射流的升力效益降低。定常射流對(duì)應(yīng)的失速迎角也進(jìn)一步提前，由15°提前至7.5°，同時(shí)，脈沖射流升力系數(shù)大于定常射流的轉(zhuǎn)折迎角提前至2.5°。此時(shí)相較定常射流，脈沖射流在大迎角下有更突出的升力表現(xiàn)，在15°迎角下，ΔCL/Cμ=36.7。隨著動(dòng)量系數(shù)的增加，定常射流下的失速迎角會(huì)進(jìn)一步提前，嚴(yán)重制約定常射流的升力效益。而脈沖射流大大擴(kuò)寬了飛行器的可用迎角，甚至超過了無射流狀態(tài)下的失速迎角，隨動(dòng)量系數(shù)增大，這種優(yōu)勢變得更加突出，完全彌補(bǔ)了定常射流在大迎角下存在的缺陷和不足。

從無射流狀態(tài)和基準(zhǔn)翼型NASA LS（1）-0417的升力系數(shù)對(duì)比中可以看出，修形后的環(huán)量控制翼型在無射流狀態(tài)下升力系數(shù)比NASA LS（1）-0417高，添加射流后升力系數(shù)進(jìn)一步提高。

為了更加全面地衡量脈沖射流、定常射流、無射流和基準(zhǔn)翼型4種狀態(tài)下的氣動(dòng)效益，并合理考量射流產(chǎn)生所需的能量消耗，參考文獻(xiàn)［24］，假設(shè)射流是通過安裝在機(jī)翼內(nèi)部的氣壓泵實(shí)現(xiàn)的，氣壓泵吸入外界自由流動(dòng)的氣體后，經(jīng)過壓縮送入噴射槽形成射流。因此，能量消耗功率可由質(zhì)量流量和總焓變化表示：

式中：Pc，jet為射流產(chǎn)生的能量消耗功率系數(shù)；Pt2和Tt2分別為自由來流的總壓和總溫；Pt1為射流出口處的總壓；η為氣壓泵的工作效率，本節(jié)取0.85。脈沖射流采用振蕩吹氣執(zhí)行器，即通過交替打開和關(guān)閉氣壓泵與噴射槽之間的連接，或打開和關(guān)閉氣壓泵與流體吸入通道之間的連接，通過調(diào)節(jié)閥門開度能夠?qū)崿F(xiàn)穩(wěn)定或振蕩射流。

為綜合分析氣動(dòng)特性和能量消耗特性對(duì)射流控制效果的影響，采用Seifert［6］提出的氣動(dòng)品質(zhì)因子AFM1進(jìn)行衡量。該指標(biāo)以升阻比為衡量基準(zhǔn)，將射流功率轉(zhuǎn)換為等效阻力，該值越大，則表明為執(zhí)行器提供功率比為動(dòng)力裝置提供功率更有效，在增升減阻方面更有利，并表征射流控制效率越高。氣動(dòng)品質(zhì)因子AFM1定義為

式中：（L/D）baseline為環(huán)量控制翼型無射流狀態(tài)下的升阻比。

圖12為不同情況下的AFM1對(duì)比曲線。環(huán)量控制翼型較基準(zhǔn)翼型NASA LS（1）-0417而言，升阻比較低，這表明修形后的環(huán)量控制翼型盡管升力系數(shù)較大，但對(duì)應(yīng)的阻力系數(shù)增長更為明顯，導(dǎo)致升阻比降低，特別是在大迎角下降低明顯。

圖12 氣動(dòng)品質(zhì)因子對(duì)比曲線Fig.12 Aerodynamic quality factor comparison

對(duì)定常射流而言，隨迎角增大，控制效率逐漸降低，然后基本保持不變。對(duì)脈沖射流而言，隨迎角增大，控制效率先降低，然后迅速提高，逐漸超越定常射流。這是因?yàn)樵谥械陀窍拢ǔＩ淞魈嵘ψ饔妹黠@，而耗能低于脈沖射流，控制效率較高；隨迎角增大，定常射流對(duì)升力的影響逐漸減弱，同時(shí)，翼型提前進(jìn)入失速狀態(tài)，控制效率降低。而脈沖射流在大迎角下升力系數(shù)提升明顯，失速迎角推遲，控制效率逐漸增高。在考慮了射流消耗的情況下，脈沖射流仍能夠部分改善環(huán)量控制翼型在大迎角下的氣動(dòng)效益，這與上文的研究是一致的。

2.4 疊加效應(yīng)的影響

上述研究已經(jīng)證明了脈沖射流在大迎角下卓越的性能優(yōu)勢，解決了定常射流在高動(dòng)量系數(shù)、大迎角下升力增益較差的問題。但是由脈沖射流所帶來的升力系數(shù)脈動(dòng)現(xiàn)象仍未完全解決，本節(jié)提出采用脈沖射流與定常射流相疊加的方式來降低脈沖射流的升力脈動(dòng)現(xiàn)象。維持其他參數(shù)不變，其中脈沖射流頻率為200 Hz，DC=70%，Cμ=0.04。

定義定常射流部分與脈沖射流部分的比例為占比rsp，即rsp=Vjet.min/Vjet.max。疊加效應(yīng)即研究占比rsp對(duì)氣動(dòng)特性的影響。圖13為在15°迎角下，定常射流部分與脈沖射流部分不同占比rsp對(duì)應(yīng)下的時(shí)均升力系數(shù)和脈動(dòng)幅值變化趨勢。隨rsp的增加，時(shí)均升力系數(shù)逐漸下降，并且呈現(xiàn)先慢后快的下降速度。這表明適當(dāng)添加定常射流并不會(huì)大幅度降低時(shí)均升力系數(shù)，同時(shí)能夠有效降低單一脈沖射流時(shí)的升力系數(shù)脈動(dòng)幅值。當(dāng)rsp為0.84時(shí)，時(shí)均升力系數(shù)僅降低約12%，升力系數(shù)脈動(dòng)幅值降低約80%，此時(shí)升力系數(shù)脈動(dòng)幅值約占時(shí)均升力系數(shù)的3.7%，已經(jīng)遠(yuǎn)遠(yuǎn)低于正常狀態(tài)下大迎角下升力系數(shù)的脈動(dòng)幅值［25］，能夠滿足飛行器在大迎角下的氣動(dòng)力要求。當(dāng)rsp約為0.8時(shí)，有效解決了脈沖射流帶來的氣動(dòng)力脈動(dòng)問題，擴(kuò)寬了環(huán)量控制技術(shù)的可用迎角范圍。

圖13 不同r sp下的時(shí)均升力系數(shù)與脈動(dòng)幅值對(duì)比曲線Fig.13 Comparison of lift coefficient and pulsation amplitude at different rsp

圖14給出了rsp分別為0.566和0.84時(shí)的后緣流線變化圖。從整體上來看，后緣流場中的脫體渦都經(jīng)歷了生成、后移、脫落的過程。當(dāng)rsp=0.566時(shí)，定常射流占比較低，定常射流速度較小，無法突破后緣外側(cè)流動(dòng)到后緣下方，當(dāng)產(chǎn)生高速的脈沖射流后，射流需要從翼型后緣上方移動(dòng)到后緣下方，導(dǎo)致升力系數(shù)脈動(dòng)程度劇烈；當(dāng)rsp=0.84時(shí)，從整體上來看，后緣流場中的脫體渦移動(dòng)并不明顯，此時(shí)定常射流占比較高，定常射流速度較大，已經(jīng)突破后緣外側(cè)流動(dòng)到了后緣下方，此時(shí)在整個(gè)射流階段，射流一直在翼型后緣下方進(jìn)行移動(dòng)，升力系數(shù)脈動(dòng)減弱，同時(shí)，能夠較大程度地保留脈沖射流不斷促進(jìn)翼型上表面的分離渦脫落特性，因此升力系數(shù)的下降并不明顯。

3 結(jié) 論

本文采用非定常數(shù)值模擬方法，對(duì)脈沖射流作用下的環(huán)量控制翼型進(jìn)行了相應(yīng)研究，特別是分析了脈沖射流的實(shí)際應(yīng)用前景，得到以下結(jié)論：

1）維持動(dòng)量系數(shù)不變的情況下，較定常射流而言，占空比越低，脈沖射流減小質(zhì)量流量的優(yōu)勢越明顯，但同時(shí)升力系數(shù)脈動(dòng)的劣勢也越明顯。

2）在占空比不變的情況下，小迎角下升力系數(shù)隨脈沖頻率增加而呈現(xiàn)先增加后降低的趨勢，但整體變化幅度不大；大迎角下升力系數(shù)隨脈沖頻率增加而持續(xù)增加，高頻脈沖射流在大迎角下有突出的應(yīng)用優(yōu)勢。

3）脈沖射流作用下，翼型上表面的分離渦和后緣尾渦通過形成、后移、脫落、消失、再形成的過程，延遲了翼型失速，改善了翼型在大迎角下的氣動(dòng)特性，并且隨動(dòng)量系數(shù)的增大，這種優(yōu)勢更加明顯。

4）采用脈沖射流和定常射流相疊加的方式能有效減弱升力系數(shù)脈動(dòng)幅值過大的問題，當(dāng)rsp約為0.8時(shí)，能有效擴(kuò)寬環(huán)量控制技術(shù)的使用迎角。