王思維, 王掩剛, 劉漢儒, 陳延俊, 劉淑麗

1.西北工業(yè)大學(xué) 動力與能源學(xué)院, 陜西 西安 710072; 2.中國航空工業(yè)空氣動力研究院, 黑龍江 哈爾濱 150001

隨著航空電氣化及民用航空動力的不斷發(fā)展,未來城市交通系統(tǒng)必將由地面延伸至廣闊的空域以緩解當(dāng)前的城市交通壓力。這一應(yīng)用前景對未來城市空中交通工具提出了嚴(yán)苛的要求,既要具備垂直起降功能又要限制噪音水平[1]。與開式螺旋槳相比,涵道風(fēng)扇在低速領(lǐng)域中優(yōu)秀的氣動性能[2-3]以及更低的噪音水平得到了國內(nèi)外學(xué)者的廣泛關(guān)注[4-7]。當(dāng)飛行器處于懸停或垂直起降狀態(tài)時(shí),涵道風(fēng)扇將為飛行器提供升力從而被稱為升力風(fēng)扇[8-9]。涵道唇口作為升力風(fēng)扇的關(guān)鍵氣動部件,理想情況下可以為升力風(fēng)扇提供50%的升力[10-12],因此對涵道唇口的深入研究有利于進(jìn)一步推進(jìn)未來航空動力的發(fā)展[13-14]。

涵道唇口的增升效應(yīng),最早是由Stipa[15]在1932年通過試驗(yàn)研究得到的,研究發(fā)現(xiàn)給螺旋槳外加裝環(huán)形翼型的涵道之后可以提高靜止?fàn)顟B(tài)下的升力。后來Platt[6]在靜態(tài)試驗(yàn)中測得,安裝環(huán)形翼型的涵道后,風(fēng)扇的升力提高了一倍。Zhang等[16]通過理論分析指出,等直徑的升力風(fēng)扇,在提供相同升力時(shí),其功率僅為同直徑螺旋槳的70%,在消耗相同功率時(shí),其升力為同直徑螺旋槳的1.26倍,這一結(jié)論體現(xiàn)了升力風(fēng)扇的性能優(yōu)勢。

Taylor[17]以涵道螺旋槳為研究對象,用試驗(yàn)方法研究了圓形唇口半徑對升力的影響,得出圓形唇口半徑是影響唇口升力的主要參數(shù)之一。當(dāng)唇口半徑大于螺旋槳直徑的6%時(shí),唇口升力保持不變,小于6%時(shí),唇口升力開始降低。Graf等[18]以橢圓、圓等幾何形狀作為唇口截面曲線,設(shè)計(jì)了5種唇口并分別研究了它們的力學(xué)特性,結(jié)果表明大半徑的唇口在懸停狀態(tài)下有更好的升力性能,小半徑的唇口在橫風(fēng)的狀態(tài)下有更好的氣動穩(wěn)定性。Deng等[19]通過試驗(yàn)方法測得升力風(fēng)扇雙紐線唇口的壓力分布,指出轉(zhuǎn)子在風(fēng)扇唇口誘導(dǎo)的低壓區(qū)使唇口產(chǎn)生了升力。

雖然唇口升力與其前緣的半徑有關(guān),但是在工程應(yīng)用中,唇口的厚度和高度直接影響著升力風(fēng)扇的徑向尺寸及軸向長度[20]。因此如何在有限唇口厚度和高度約束下實(shí)現(xiàn)升力的最大化,成為了研究人員必須面對的問題之一。Wang等[10]應(yīng)用面元方法對升力風(fēng)扇唇口進(jìn)行了優(yōu)化設(shè)計(jì)。葉坤等[21]采用動量源方法進(jìn)行了涵道氣動力的計(jì)算,分別用響應(yīng)面模型和神經(jīng)網(wǎng)絡(luò)模型對涵道進(jìn)行優(yōu)化設(shè)計(jì)。但是他們的優(yōu)化方法均是通過不斷調(diào)整唇口截面曲線控制點(diǎn)位置,并基于優(yōu)化算法完成唇口的優(yōu)化工作。這種方法雖然可以完成唇口的優(yōu)化,但是其優(yōu)化結(jié)果都是基于大量樣本及優(yōu)化策略。為了在進(jìn)行唇口的優(yōu)化設(shè)計(jì)之前可以確定性能更好的基礎(chǔ)唇口,以提高唇口優(yōu)化效率,需要建立升力風(fēng)扇唇口特征幾何與其氣動性能的關(guān)聯(lián)性。

本文從構(gòu)成涵道唇口的特征幾何出發(fā),設(shè)計(jì)了4種唇口,在試驗(yàn)校核的基礎(chǔ)上,采用數(shù)值計(jì)算方法對比研究了4種唇口在懸停狀態(tài)下的性能差異并進(jìn)行了流動分析。唇口升力由氣流在其表面各微元面積處氣動力軸向分量的合力產(chǎn)生,斜率反映了各微元軸向分力占?xì)鈩恿Φ谋戎?。因?以斜率對4種唇口截面曲線進(jìn)行參數(shù)化表示,研究了曲線斜率對性能參數(shù)的影響規(guī)律,并對這一規(guī)律進(jìn)行了驗(yàn)證。以此規(guī)律設(shè)計(jì)了新的唇口曲線,并研究了新唇口對升力風(fēng)扇的總升力貢獻(xiàn),從而確認(rèn)了研究成果的工程應(yīng)用價(jià)值。

1 研究方法與對象

1.1 計(jì)算方法校驗(yàn)

升力風(fēng)扇動力裝置作為一個(gè)整體,各部件在共同工作中存在著相互耦合關(guān)系,唇口為轉(zhuǎn)子提供良好進(jìn)氣條件的同時(shí),轉(zhuǎn)子驅(qū)動氣流在唇口表面繞流使唇口產(chǎn)生升力。因此,本文以某升力風(fēng)扇作為計(jì)算方法的校驗(yàn)對象,采用內(nèi)外流同時(shí)求解的方式對其進(jìn)行數(shù)值計(jì)算,以獲取更加貼合真實(shí)流動的仿真結(jié)果。升力風(fēng)扇周圍流場沿周向均勻分布,因此為了降低網(wǎng)格數(shù)量以提高計(jì)算效率,采用圖1a)所示的計(jì)算域,僅對1/2的靜止域進(jìn)行數(shù)值計(jì)算,并且對轉(zhuǎn)風(fēng)扇葉片所在旋轉(zhuǎn)域采用單通道的計(jì)算域。

圖1 計(jì)算方法校驗(yàn)

為了使升力風(fēng)扇周圍的流場充分發(fā)展,上游總壓進(jìn)口邊界、側(cè)面的進(jìn)出口邊界和下游靜壓出口邊界與升力風(fēng)扇的距離分別為8倍、5倍和12倍的涵道直徑。其中對轉(zhuǎn)風(fēng)扇葉片所在空間為旋轉(zhuǎn)計(jì)算域,外流場計(jì)算域?yàn)殪o止計(jì)算域,采用混合平面模型處理動靜交界面和動動交界面的數(shù)值傳遞問題。具體參數(shù)設(shè)置上,環(huán)境總溫設(shè)置為300 K,參考壓力為101 325 Pa,計(jì)算采用二階中心格式結(jié)合SSTk-ω湍流模型。使用結(jié)構(gòu)化網(wǎng)格分別對對轉(zhuǎn)風(fēng)扇及遠(yuǎn)場計(jì)算域進(jìn)行空間離散,第一層網(wǎng)格劃分保證了y+值小于2。轉(zhuǎn)子域網(wǎng)格劃分時(shí)考慮 0.5 mm 轉(zhuǎn)子葉頂間隙的影響。

調(diào)整風(fēng)扇葉片與唇口表面的網(wǎng)格尺寸,分別劃分了350萬、400萬、500萬以及600萬的網(wǎng)格數(shù)量進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證。設(shè)置R1與R2計(jì)算域的旋轉(zhuǎn)速度分別為10 000 r/min與-10 000 r/min。計(jì)算結(jié)果如圖1b)所示,隨著網(wǎng)格數(shù)量的增多,唇口升力及涵道流量均逐漸增大,定量來看,當(dāng)網(wǎng)格數(shù)量由400萬增加到500萬時(shí),唇口升力與涵道流量分別增大了1.2%與0.7%。網(wǎng)格數(shù)量增多會伴隨著計(jì)算資源的大量消耗,為了平衡計(jì)算精度與計(jì)算資源之間的矛盾,選擇500萬網(wǎng)格的配置策略作為本文后續(xù)的網(wǎng)格劃分依據(jù)。

上文所述的校驗(yàn)方法雖然保證了網(wǎng)格劃分策略的準(zhǔn)確性,但是計(jì)算方法的正確性仍需要進(jìn)一步驗(yàn)證。采用圖2所示的試驗(yàn)裝置完成升力風(fēng)扇的升力測量。唇口與涵道殼體之間留有2 mm的安裝間距,采用柔性薄膜將兩部分結(jié)構(gòu)進(jìn)行連接,從而既保證了流道的完整性,又留有涵道與唇口間的變形余量,使兩部分結(jié)構(gòu)因受力產(chǎn)生軸向變形時(shí)不會進(jìn)行力的傳遞。使用ATI Mini40和ATI Mini45高精度力學(xué)天平分別采集唇口及涵道主體的升力值。試驗(yàn)時(shí)通過PID調(diào)節(jié)使R1與R2保持同一轉(zhuǎn)速,因此下文僅用一個(gè)轉(zhuǎn)速描述對轉(zhuǎn)風(fēng)扇的旋轉(zhuǎn)速度。

圖2 唇口升力試驗(yàn)

數(shù)值計(jì)算與試驗(yàn)總計(jì)對比了8 000,9 000,10 000 r/min這3組轉(zhuǎn)速下的唇口升力。結(jié)果對比如圖3所示,3組轉(zhuǎn)速的數(shù)值計(jì)算與試驗(yàn)結(jié)果的差異均在3%以內(nèi)。由此進(jìn)一步驗(yàn)證了前文所述的網(wǎng)格劃分策略及數(shù)值計(jì)算方法的正確性。

圖3 唇口升力結(jié)果對比

1.2 研究對象

從以往的研究中可發(fā)現(xiàn),研究人員常常選取3種類型的幾何形狀作為唇口截面曲線。Stipa[15]和Platt[6]選擇翼型作為唇口的截面曲線,Wang[10]和Deng[19]選擇雙紐線作為唇口的截面曲線,Taylor[17]和Graf[18]選擇橢圓或圓作為唇口的截面曲線。本文則同時(shí)選取以上3類幾何形狀作為唇口截面曲線。以往的研究中各位學(xué)者選擇的翼型不盡相同,本文使用的翼型選擇為NACA0012及NACA65。

各唇口截面曲線的原始及最終幾何形狀如圖4所示,為了使NACA0012及NACA65翼型滿足厚度及高度的約束條件,本文僅選取翼型前緣處滿足厚度與高度比值的一段作為唇口的截面曲線形狀。選取四分之一雙紐線、二分之一橢圓作為唇口截面的部分曲線,其余部分用曲線光滑過渡至直線段。為了明確各構(gòu)成唇口截面曲線的原始形狀,下文均采用NACA0012、NACA65、雙紐線和橢圓指代所設(shè)計(jì)的4種唇口。

為了避免唇口厚度及高度差異對結(jié)果的影響,設(shè)計(jì)時(shí)保證4種曲線的厚度及高度分別為20 mm和31 mm。本文的研究對象,除了唇口截面曲線存在差異外,進(jìn)口導(dǎo)流錐、轉(zhuǎn)子以及涵道殼體的幾何參數(shù)均保持一致。詳細(xì)參數(shù)如圖4和表1所示。

表1 涵道唇口幾何參數(shù)

2 結(jié)果與討論

針對所設(shè)計(jì)的4種唇口,應(yīng)用經(jīng)過試驗(yàn)校驗(yàn)的數(shù)值分析方法,開展了無來流條件下唇口性能參數(shù)及流動研究。使用無量綱升力LND(non-dimensional lift)及無量綱流量QND(non-dimensional mass flow)對唇口性能進(jìn)行定量表示。

式中:Ln及Lref,n分別是轉(zhuǎn)速為n時(shí)唇口的實(shí)際升力與參考升力,升力方向沿涵道軸向;Qn及Qref,n分別是轉(zhuǎn)速為n時(shí)通過唇口的實(shí)際流量與參考流量。

參考升力及參考流量的定義為

式中,pn為唇口出口截面處的平均相對靜壓力。

2.1 結(jié)果分析

在數(shù)值計(jì)算方法被充分驗(yàn)證的基礎(chǔ)上,得到圖5所示的無量綱流量及無量綱升力特性曲線。無量綱流量與無量綱升力結(jié)果均表明,雙紐線唇口無量綱參數(shù)最高,其次是橢圓唇口,最后是NACA0012和NACA65唇口。在10 000 r/min的轉(zhuǎn)速下,橢圓唇口的無量綱升力與無量綱流量比NACA65唇口分別提高33%與9.5%。

圖5 性能參數(shù)

翼型之所以不適合作為升力風(fēng)扇的唇口曲線,是因?yàn)榇娇诹鲃咏Y(jié)構(gòu)與二維機(jī)翼的流動結(jié)構(gòu)完全不同。二維機(jī)翼的流場為氣流在機(jī)翼前緣處分開,沿機(jī)翼上下表面流動并在尾緣處匯集,是一種由前緣向尾緣的流動。而氣流在升力風(fēng)扇唇口表面的流動為氣流由唇口外表面出發(fā),流過唇口前緣并向內(nèi)表面流動,是一種繞前緣的流動。除上述原因之外,也不排除因?yàn)楸疚膬H選取了翼型的前緣區(qū)域用于研究從而使翼型唇口性能參數(shù)降低的可能性。但這一問題并不是本文研究重點(diǎn),后續(xù)不再進(jìn)行深入討論。

對于升力風(fēng)扇唇口,氣流沿唇口周向均勻分布,因此僅需分析通過中軸線截面上的流動信息即可。10 000 r/min的轉(zhuǎn)速下,各唇口截面速度云圖如圖6所示。

圖6 截面速度云圖

氣流在雙紐線及橢圓唇口表面形成了明顯的加速區(qū),兩者加速區(qū)的差異在于橢圓唇口表面的高速區(qū)位于唇口前緣附近,而雙紐線唇口表面的高速區(qū)更靠近唇口內(nèi)部。相比于雙紐線和橢圓唇口,氣流在2種翼型唇口表面的加速作用明顯減弱,氣流加速區(qū)范圍明顯減小。正是氣流在唇口表面加速作用的差異引起了無量綱升力的差異。

從圖6中同樣可以看出,4種唇口內(nèi)表面均形成了流動分離區(qū)。其中NACA65和NACA0012唇口表面分離區(qū)最強(qiáng),橢圓和雙紐線唇口最弱,這一趨勢與唇口性能參數(shù)趨勢相一致。說明唇口內(nèi)表面的流動分離影響著性能參數(shù)的變化,但是流動分離的產(chǎn)生原因仍需進(jìn)一步研究。

流動分離一般由逆壓梯度產(chǎn)生,圖7顯示了各唇口表面的相對靜壓力分布,4種唇口外表面的壓力分布幾乎一致,說明了唇口外表面曲線差異不會對其流動產(chǎn)生影響。

橢圓、NACA0012以及NACA65唇口的壓力峰值均位于唇口前緣頂點(diǎn)處,雙紐線唇口的壓力峰值則更靠近唇口內(nèi)部。因此,氣流在橢圓、NACA0012以及NACA65唇口內(nèi)表面的逆壓流動自前緣頂點(diǎn)開始。而雙紐線唇口表面的逆壓流動更靠近唇口內(nèi)部,這一流動現(xiàn)象延緩了雙紐線唇口內(nèi)表面的流動分離,由此解釋了雙紐線唇口氣流加速區(qū)更靠近唇口內(nèi)部的原因。

雙紐線唇口的氣流加速區(qū)及分離起始位置更靠近唇口內(nèi)部,低壓區(qū)覆蓋更大的面積,使其升力高于其他3種唇口。對于均從前緣頂點(diǎn)開始具有逆壓流動的NACA0012、NACA65及橢圓唇口,壓力峰值越大說明了在唇口可提供更大的升力,因此壓力峰值的大小也表示了升力的大小,圖7中壓力峰值與升力由大到小分別為橢圓、NACA0012和NACA65。

唇口可為升力風(fēng)扇提供兩方面的作用,一方面是提供升力,另一方面則是對進(jìn)口氣流進(jìn)行整理,改善轉(zhuǎn)子的進(jìn)氣品質(zhì)。進(jìn)口氣流的均勻性越好意味著進(jìn)氣品質(zhì)越高,為了對比各唇口對進(jìn)口氣流均勻性的影響,從圖8所示的唇口出口處軸向速度沿徑向的分布可以看出,氣流在唇口出口截面處的速度差異僅存在于唇口內(nèi)表面附近,由于黏性的作用,唇口內(nèi)表面處的氣流速度將降至0。結(jié)合圖6中的云圖可以得到,流動分離使靠近各唇口內(nèi)表面附近的速度分布產(chǎn)生了差異,并且分離區(qū)的存在減小了涵道的有效流通面積,降低了氣流流通能力,使各唇口的無量綱流量產(chǎn)生差異。

圖8 出口截面軸向速度沿徑向分布

完成對4種唇口的性能參數(shù)及流動現(xiàn)象的分析,得到了影響唇口性能參數(shù)的流場信息。但是,無法得出唇口特征幾何與其氣動性能的影響關(guān)系。

結(jié)合圖4～5可以發(fā)現(xiàn)唇口性能參數(shù)與幾何形狀之間存在著關(guān)聯(lián)性,可以定性看出,在唇口頂部,唇口截面內(nèi)表面曲線越高,其無量綱流量與無量綱升力越大。而唇口截面內(nèi)表面曲線越高,意味著自唇口頂部開始,其截面曲線切線與X軸夾角的變化率越小。在數(shù)學(xué)中,斜率表示了曲線的切線關(guān)于橫坐標(biāo)的傾斜程度,因此,本文提出以斜率對唇口截面曲線進(jìn)行參數(shù)化表示,研究斜率與性能參數(shù)之間的關(guān)聯(lián)性。

2.2 斜率對唇口性能參數(shù)的影響分析

唇口升力為唇口表面的各微元面積處產(chǎn)生的力在涵道軸線方向分力的合力,而各微元面積處沿軸向的力則是由唇口表面壓力及曲線切角決定,表達(dá)式為

Llip=∑Δp·Δs·cosθ

(6)

式中:Δp為微元上的相對靜壓力;Δs為微元面積;θ為微元處的切線與涵道徑向的夾角。

唇口曲線斜率

(7)

式中,dZ與dX分別為唇口曲線在Z與X方向上的微分。

對于唇口曲線,θ滿足-900≤θ≤900,同時(shí)滿足

(8)

因此可用斜率k即tanθ的相對大小對唇口曲線進(jìn)行參數(shù)化表示,進(jìn)而研究其對唇口性能參數(shù)的影響規(guī)律。

由2.1節(jié)的研究發(fā)現(xiàn),當(dāng)Z小于等于24 mm時(shí),4種唇口均已產(chǎn)生流動分離現(xiàn)象,所以分離的起始位置均位于Z大于24 mm處。根據(jù)公式(7)計(jì)算得到唇口曲線斜率,如圖9所示。唇口外表面曲線不會對其流動產(chǎn)生影響,因此下文僅關(guān)心唇口內(nèi)表面曲線斜率(大于等于0部分)與其性能參數(shù)之間的影響關(guān)系。

圖9 唇口曲線斜率

從圖9中可以看出,唇口內(nèi)表面斜率由大到小分別為NACA65唇口、NACA0012唇口、橢圓唇口和雙紐線唇口,這一結(jié)果與唇口性能參數(shù)形成了一定的關(guān)聯(lián)性,即唇口內(nèi)表面曲線斜率越小,其性能越好。對于關(guān)系式(6)而言,唇口斜率越小意味著cosθ越大,同樣也會使唇口升力增大。

2.3 斜率對唇口性能參數(shù)影響規(guī)律驗(yàn)證

由于上文的結(jié)論是在4組有限數(shù)據(jù)的基礎(chǔ)上得到的,結(jié)果不具有普適性。為了驗(yàn)證上述猜想,以雙紐線為基礎(chǔ),設(shè)計(jì)多組不同斜率的唇口截面曲線對其進(jìn)行對比研究。新曲線的設(shè)計(jì)依據(jù)為,保證唇口外表面曲線不變,即圖10中C區(qū)域(89

圖10 雙紐線唇口區(qū)域劃分

使用公式(9)中的n次多項(xiàng)式對雙紐線B區(qū)域的斜率進(jìn)行表示

f′=b0+b1x+b2x2+…+bnxn

(9)

新設(shè)計(jì)的唇口曲線斜率由系數(shù)λ確定,即新曲線在B區(qū)域的斜率表達(dá)式為

(10)

新設(shè)計(jì)的唇口曲線在B區(qū)域的表達(dá)式為

(11)

其滿足邊界條件fb(89)=31。

新設(shè)計(jì)的曲線在A區(qū)域的表達(dá)式由公式(12)中的3次多項(xiàng)式進(jìn)行參數(shù)化

fa=α+a0x+a1x2+a2x3

(12)

其滿足邊界條件fa(77)=0,fa(78)=fb(78)。

采用上述方法,令系數(shù)λ分別等于4/7,6/7,8/7,10/7,12/7,14/7,16/7,得到新的曲線斜率及曲線幾何如圖11所示。當(dāng)λ=16/7時(shí),新設(shè)計(jì)的曲線斜率大于橢圓唇口斜率,當(dāng)8/7≤λ≤14/7時(shí),曲線斜率位于橢圓唇口及雙紐線唇口斜率中間,當(dāng)λ≤6/7時(shí),曲線斜率小于雙紐線唇口斜率。

圖11 新唇口設(shè)計(jì)

新設(shè)計(jì)唇口的性能參數(shù)對比結(jié)果如圖12所示,可以看出,當(dāng)唇口內(nèi)表面斜率小于雙紐線斜率時(shí),唇口性能參數(shù)均差于原始的雙紐線唇口。當(dāng)唇口內(nèi)表面斜率大于雙紐線斜率時(shí),唇口性能參數(shù)均好于原始的雙紐線唇口,且隨著唇口斜率的增大,其性能參數(shù)逐漸變好。但是,當(dāng)系數(shù)λ≥14/7時(shí),新設(shè)計(jì)的唇口性能參數(shù)不再有明顯提升,甚至當(dāng)λ=16/7時(shí),其無量綱升力低于λ=14/7時(shí)的值。

圖12 新設(shè)計(jì)唇口性能參數(shù)對比

結(jié)合2.2節(jié)內(nèi)容可以得到,對于本文的研究對象,雙紐線及橢圓唇口的曲線斜率之間存在最優(yōu)的唇口曲線斜率,這一唇口曲線斜率為2倍雙紐線唇口斜率。當(dāng)唇口內(nèi)表面斜率大于該斜率時(shí),隨唇口內(nèi)表面斜率的減小唇口性能參數(shù)逐漸變好,當(dāng)小于該斜率時(shí),隨唇口內(nèi)表面斜率的減小唇口性能參數(shù)逐漸變差。

上述研究仍留有一個(gè)問題需要解決,即當(dāng)λ=14/7和16/7時(shí),橢圓唇口曲線斜率位于新設(shè)計(jì)的唇口曲線斜率之間,但是其性能參數(shù)明顯差于新設(shè)計(jì)的唇口曲線,那么上文中提到的斜率與性能間的影響規(guī)律是否成立需要進(jìn)一步分析。為了分析該現(xiàn)象的產(chǎn)生原因,獲取λ=4/7,8/7,14/7,16/7唇口與橢圓唇口表面相對靜壓力分布,進(jìn)行如圖13所示的對比?？梢钥吹?λ=14/7和16/7時(shí)的唇口內(nèi)表面壓力分布與橢圓唇口幾乎重合。由此可以得到,唇口斜率直接影響了內(nèi)表面壓力分布,進(jìn)而影響著性能參數(shù)的變化。

圖13 唇口表面壓力分布對比

雖然λ=14/7,16/7與橢圓唇口曲線的斜率和壓力分布幾乎一致,但是從圖11b)中可以看出,λ=14/7,16/7唇口前緣頂點(diǎn)位于橢圓唇口前緣頂點(diǎn)右側(cè),說明λ=14/7,16/7唇口內(nèi)表面占據(jù)更大的X方向上的長度,從而使其Δs比橢圓唇口更大,這也就是λ=14/7,16/7唇口性能參數(shù)好于橢圓唇口的原因。從這一結(jié)果中可以得到啟示,在進(jìn)行升力風(fēng)扇唇口設(shè)計(jì)時(shí),可以適當(dāng)壓縮唇口外表面所占的厚度(X方向上的長度),即盡可能將唇口前緣頂點(diǎn)向唇口外側(cè)移動以增大內(nèi)表面面積從而提高其性能參數(shù)。

對比分析λ=4/7,8/7,14/7這3種唇口的壓力分布可以看出,唇口斜率對其內(nèi)表面壓力分布的影響規(guī)律為:隨著斜率的減小,壓力峰值會逐漸降低,同時(shí)伴隨著壓力峰值逐漸向唇口內(nèi)部移動。這一結(jié)果將有助于從曲線斜率出發(fā)對唇口表面壓力梯度進(jìn)行控制,從而控制流動分離現(xiàn)象。

2.4 唇口對升力風(fēng)扇的升力貢獻(xiàn)

雖然上文發(fā)現(xiàn)并驗(yàn)證了斜率對唇口性能參數(shù)影響規(guī)律,但是升力風(fēng)扇作為一個(gè)整體的動力單元,當(dāng)唇口的升力貢獻(xiàn)發(fā)生改變時(shí),其對升力風(fēng)扇的升力又會產(chǎn)生何種影響需要進(jìn)行進(jìn)一步的驗(yàn)證。

圖14為采用2種唇口升力風(fēng)扇各部件的升力貢獻(xiàn),可以看出,當(dāng)升力風(fēng)扇采用基于2倍雙紐線斜率設(shè)計(jì)的唇口時(shí),其總升力和唇口升力均高于采用雙紐線唇口升力風(fēng)扇的總升力和唇口升力,且在10 000 r/min的轉(zhuǎn)速下,其總升力比雙紐線唇口升力風(fēng)扇的總升力提高5.6%。

圖14 升力風(fēng)扇各部件升力貢獻(xiàn)

3 結(jié) 論

本文對構(gòu)成升力風(fēng)扇唇口的截面曲線進(jìn)行研究,分析了4種唇口曲線的性能參數(shù)差異,使用斜率對唇口進(jìn)行參數(shù)化表示,研究了斜率與唇口性能參數(shù)的關(guān)聯(lián)關(guān)系,得到了以下結(jié)論:

1) 對雙紐線、橢圓、NACA0012和NACA65作為截面曲線的唇口性能參數(shù)進(jìn)行了對比研究,結(jié)果表明,在10 000 r/min的轉(zhuǎn)速下,雙紐線唇口的無量綱升力與無量綱流量分別比NACA65唇口提高33%與9.5%。

2) 唇口內(nèi)表面曲線斜率與性能參數(shù)的關(guān)聯(lián)關(guān)系為:存在2倍雙紐線斜率的最優(yōu)唇口斜率,當(dāng)唇口內(nèi)表面斜率大于或小于該斜率時(shí),唇口的性能參數(shù)均會降低。

3) 與基于雙紐線設(shè)計(jì)的唇口相比,基于2倍雙紐線斜率設(shè)計(jì)的唇口,可使升力風(fēng)扇在10 000 r/min轉(zhuǎn)速下的總升力提高5.6%。