齊龍舟,馮和英*,彭葉輝,趙鯤
1 湖南科技大學(xué) 機(jī)械設(shè)備健康維護(hù)湖南省重點(diǎn)實(shí)驗(yàn)室,湖南 湘潭 411201
2 湖南科技大學(xué) 數(shù)學(xué)與計(jì)算科學(xué)學(xué)院,湖南 湘潭 411201
3 中國(guó)空氣動(dòng)力研究與發(fā)展中心 氣動(dòng)噪聲控制重點(diǎn)實(shí)驗(yàn)室,四川 綿陽(yáng) 621000
艦載機(jī)起飛時(shí),尾噴管?chē)姵龅母咚?、高溫射流?huì)嚴(yán)重影響其附近區(qū)域的人員和設(shè)備,因此,航母甲板上通常設(shè)有偏流板用來(lái)改變高溫氣流的方向,以達(dá)到保護(hù)人員及設(shè)備安全的目的。
由于尾噴管射流速度可達(dá)到音速甚至是超音速,致使偏流板沖擊噪聲最高可達(dá)140 dB以上[1],這不僅增加了周?chē)h(huán)境的噪聲污染,嚴(yán)重威脅長(zhǎng)時(shí)間、近距離工作人員(聲壓級(jí)在90 dB以上將影響人的聽(tīng)覺(jué))的身心健康,而且還會(huì)導(dǎo)致飛行器機(jī)體結(jié)構(gòu)的疲勞損傷,降低飛行器的安全性能和使用壽命。2010年,美國(guó)政府給聽(tīng)力受損的航母甲板上工作人員支付了高達(dá)14億美元的補(bǔ)貼,其中還不包括治療費(fèi)用[2]??梢?jiàn),偏流板的設(shè)計(jì)對(duì)艦載機(jī)周?chē)h(huán)境及飛行器的安全影響重大。
在艦載機(jī)偏流板設(shè)計(jì)研究領(lǐng)域,人們很早就開(kāi)始關(guān)注偏流板附近流場(chǎng)及溫度場(chǎng)對(duì)周?chē)h(huán)境與偏流板性能的影響。何慶林等[3]以國(guó)外某型艦載機(jī)和噴氣偏流板為研究對(duì)象,從理論上證明了該噴氣偏流板的設(shè)計(jì)方案能夠滿(mǎn)足艦載機(jī)起飛的需求。趙留平[4]分析了發(fā)動(dòng)機(jī)噴流沖擊偏流板的流場(chǎng)特性,得到了尾流場(chǎng)流動(dòng)特性隨偏流板角度的變化規(guī)律。郭濤[5]研究了偏流板對(duì)航母飛行甲板、艦載人員以及周?chē)h(huán)境的影響,并基于偏流板性能的計(jì)算結(jié)果分析了偏流板的安全性及設(shè)計(jì)要求。馬彩東等[6]探究了不同偏轉(zhuǎn)角對(duì)偏流板周?chē)鲌?chǎng)的影響,結(jié)果顯示當(dāng)偏轉(zhuǎn)角β=30°時(shí),尾噴管?chē)姵龅牟糠指邷厝細(xì)饨?jīng)偏流板導(dǎo)流向偏流板后方的流動(dòng)對(duì)板后設(shè)備產(chǎn)生了影響;當(dāng)β=60°時(shí),高溫燃?xì)饣亓鲗?duì)偏流板前的設(shè)備產(chǎn)生了影響。因此,現(xiàn)役航母艦載機(jī)起飛時(shí)偏流板的偏轉(zhuǎn)角皆在30°~60°之間。
上述研究主要是圍繞流動(dòng)與熱輻射問(wèn)題對(duì)偏流板進(jìn)行安全性設(shè)計(jì),近年來(lái),研究人員開(kāi)始逐漸關(guān)注艦載機(jī)的噪聲輻射問(wèn)題。謝軍龍等[7]采用大渦模擬(LES)結(jié)合FW-H聲類(lèi)比方法,研究了低速?zèng)_擊射流和非定常周期速度沖擊射流的氣動(dòng)聲場(chǎng),并利用實(shí)驗(yàn)數(shù)據(jù)檢測(cè)了低速射流沖擊斜板聲學(xué)模擬方法的適用性。黃河源等[8]從流場(chǎng)出發(fā),通過(guò)渦聲原理分析了流動(dòng)噪聲的產(chǎn)生機(jī)理。高速?zèng)_擊射流氣動(dòng)聲場(chǎng)研究的關(guān)注熱點(diǎn)主要集中在自由射流噪聲和沖擊大平板噪聲[9]上。Dauptain等[10]利用LES方法成功捕獲了超音速平板沖擊射流中的激波結(jié)構(gòu)和內(nèi)外剪切層中不同尺度的渦結(jié)構(gòu),模擬結(jié)果表明,環(huán)形激波的產(chǎn)生和消失是由馬赫盤(pán)、斜激波、內(nèi)剪切層渦和沖擊平板共同作用引起的。Martin等[11]圍繞高速射流沖擊大平板的噪聲指向性進(jìn)行了實(shí)驗(yàn)分析,實(shí)驗(yàn)結(jié)果表明,射流沖擊噪聲指向性具有方向性,當(dāng)監(jiān)測(cè)角為45°時(shí)噪聲最小,在噪聲源中捕捉到了較為典型的四極子聲源特性。
綜上所述,在艦載機(jī)起飛背景下,偏流板偏轉(zhuǎn)角度的設(shè)計(jì)關(guān)系到艦載機(jī)及其周?chē)h(huán)境的安全性。為此,本文將基于某型發(fā)動(dòng)機(jī)尾噴管的形狀尺寸,采用LES結(jié)合FW-H聲類(lèi)比方法建立超音速射流沖擊斜板聲流場(chǎng)計(jì)算模型,研究不同偏轉(zhuǎn)角(β=35°,45°,55°)對(duì)尾噴管出口流場(chǎng)、溫度場(chǎng)、聲場(chǎng)特性的影響,探索導(dǎo)流效果好、遠(yuǎn)場(chǎng)噪聲小、安全區(qū)域大的最佳偏轉(zhuǎn)角度。
本文所研究的噴管為收縮擴(kuò)張型圓口噴管,噴管入口直徑為120 mm,出口直徑D=100 mm,收縮段長(zhǎng)度為27 mm,擴(kuò)張段為38 mm。以噴管出口直徑D為長(zhǎng)度基準(zhǔn),使用Pro/E軟件建立物理模型,計(jì)算域如圖1所示。坐標(biāo)原點(diǎn)位于噴管入口圓心,計(jì)算域?qū)挾葂和高度y均為15D,長(zhǎng)度z為30D,沖擊距離X=8D,偏流板寬度為9D、高度為5D,偏轉(zhuǎn)角β為偏流板中心法線(xiàn)與地面的夾角,聲積分面寬度和高度均為L(zhǎng)=12D,長(zhǎng)度H=15D。
圖 1 物理模型及計(jì)算域Fig.1 Physical model and computational domain
采用ICEM軟件對(duì)模型進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分。因偏流板相對(duì)于計(jì)算域的尺寸來(lái)說(shuō)是一個(gè)微小量,故將偏流板簡(jiǎn)化為一無(wú)厚度的面。噴嘴內(nèi)外壁面附近、偏流板壁面附近、噴嘴到斜板間的區(qū)域以及聲源面內(nèi)部都需要采用局部網(wǎng)格進(jìn)行加密。在壁面附近,第1層網(wǎng)格厚度為0.1 mm,整體網(wǎng)格質(zhì)量在0.65以上。
采用LES結(jié)合FW-H聲類(lèi)比方法建立超音速射流沖擊斜板聲流場(chǎng)計(jì)算模型。數(shù)值計(jì)算采用Fluent軟件。由于噴射速度超過(guò)了音速,故其氣體狀態(tài)為可壓縮流體,密度選用理想氣體模型。將噴管進(jìn)口設(shè)為壓力進(jìn)口邊界條件,輸入總壓為260 000 Pa,總溫為2 030 K,靜壓為258 000 Pa。噴管內(nèi)外壁面、偏流板壁面和底面均設(shè)置為無(wú)滑移固壁邊界條件。計(jì)算域的對(duì)稱(chēng)面設(shè)置為對(duì)稱(chēng)面邊界條件,其余為壓力出口邊界條件,壓力和溫度均為環(huán)境壓力和溫度,分別為101 325 Pa和300 K。
穩(wěn)態(tài)計(jì)算采用基于密度的求解器,耦合求解矢量形式的控制方程組;計(jì)算格式采用AUSM的通量差分分裂格式,離散格式采用二階迎風(fēng)格式。穩(wěn)態(tài)計(jì)算中,湍流模型選用RNGk-?湍流模型,該模型為高雷諾數(shù)模型,認(rèn)為各項(xiàng)殘差小于10-4時(shí)計(jì)算收斂。將穩(wěn)態(tài)計(jì)算結(jié)果作為初場(chǎng),開(kāi)始非穩(wěn)態(tài)數(shù)值計(jì)算。非穩(wěn)態(tài)計(jì)算湍流模型選用LES,亞格子模型選用WALE模型,動(dòng)量方程的離散方式采用Second Order Upwind格式,為提高計(jì)算精度,時(shí)間步長(zhǎng)設(shè)置為 2.5×10-5s。
聲場(chǎng)計(jì)算采用FW-H聲類(lèi)比模型,通過(guò)快速傅里葉變換,可以獲得各個(gè)監(jiān)測(cè)點(diǎn)處聲壓級(jí)。遠(yuǎn)場(chǎng)噪聲監(jiān)測(cè)點(diǎn)的分布如圖2(圖中,Z為遠(yuǎn)場(chǎng)噪聲監(jiān)測(cè)平面高度,S為監(jiān)測(cè)距離)所示,每隔10°設(shè)置1個(gè)監(jiān)測(cè)點(diǎn),監(jiān)測(cè)點(diǎn)距離噴嘴出口軸心20 m,圍繞噴嘴出口軸心1/4圓共設(shè)置10個(gè)監(jiān)測(cè)點(diǎn)。聲場(chǎng)頻率為20 Hz~10 kHz,參考聲壓為2×10-5Pa。根據(jù)確定的邊界條件和計(jì)算參數(shù),對(duì)不同偏轉(zhuǎn)角條件下沖擊射流的流場(chǎng)與氣動(dòng)聲場(chǎng)進(jìn)行數(shù)值計(jì)算。
圖 2 噪聲監(jiān)測(cè)點(diǎn)分布Fig.2 Noise monitoring points distribution
射流速度是影響沖擊射流流場(chǎng)的一個(gè)重要參數(shù),且對(duì)網(wǎng)格數(shù)較為敏感,因此,選擇射流速度隨網(wǎng)格數(shù)的變化進(jìn)行網(wǎng)格無(wú)關(guān)性考核較為合適。圖3示出了測(cè)點(diǎn)C和測(cè)點(diǎn)D處流動(dòng)速度v隨網(wǎng)格數(shù)的變化關(guān)系,計(jì)算選取的網(wǎng)格規(guī)模分別為216萬(wàn)、253萬(wàn)、295萬(wàn)和349萬(wàn)。由圖3可知,網(wǎng)格數(shù)為349萬(wàn)和295萬(wàn)時(shí)的軸向和縱向速度基本一致,說(shuō)明自由射流段(軸向)和壁面射流段(縱向)的流場(chǎng)均趨于穩(wěn)定,295萬(wàn)網(wǎng)格數(shù)足以獲得較為準(zhǔn)確的計(jì)算結(jié)果。為安全起見(jiàn),最終網(wǎng)格總數(shù)定為318萬(wàn)。
圖 3 不同監(jiān)測(cè)點(diǎn)處流動(dòng)速度隨網(wǎng)格總數(shù)的變化Fig.3 Variation of flow velocity at different monitoring points with the total number of grids
為了驗(yàn)證所建立超聲速射流沖擊偏流板計(jì)算模型的有效性,將本文數(shù)值計(jì)算結(jié)果與文獻(xiàn)[12]和文獻(xiàn)[13]中的實(shí)驗(yàn)結(jié)果進(jìn)行了比較,計(jì)算中所采用物理模型及進(jìn)出口參數(shù)與文獻(xiàn)[12]及文獻(xiàn)[13]中一致。圖4示出了射流軸向速度時(shí)均值沿尾噴管軸線(xiàn)分布的數(shù)值計(jì)算結(jié)果及文獻(xiàn)[12]和文獻(xiàn)[13]中的實(shí)驗(yàn)結(jié)果。圖中,Y/D為噴口軸線(xiàn)坐標(biāo)與噴口直徑的比值。由圖可以看出,數(shù)值計(jì)算的軸向時(shí)均速度與實(shí)驗(yàn)結(jié)果變化趨勢(shì)一致,大部分區(qū)域吻合良好,最大誤差不超過(guò)10%,數(shù)值計(jì)算可信度較高,證明了模型的有效性與準(zhǔn)確性。
圖 4 軸向時(shí)均速度沿軸線(xiàn)的分布Fig.4 Axial time average velocity distribution along the axis
發(fā)動(dòng)機(jī)尾噴管?chē)姵龅纳淞鳛楦邷?、高速氣流,偏流板附近的溫度分布?huì)直接影響周?chē)O(shè)備及人員的安全。圖5所示為不同偏轉(zhuǎn)角β下噴管軸心至偏流板的溫度云圖(此處為側(cè)視圖)。由圖5可以看出,尾噴管?chē)姵龅囊徊糠指邷亓黧w經(jīng)過(guò)偏流板順流側(cè)導(dǎo)流后,向板后區(qū)域流動(dòng),另一部分高溫流體則撞擊偏流板逆流側(cè)后往尾噴管區(qū)域回流。通過(guò)對(duì)比不同偏轉(zhuǎn)角情況可以看出,β=35°時(shí)偏流板的導(dǎo)流效果較差,板后大部分區(qū)域的溫度都較高,這將對(duì)板后設(shè)備及人員安全造成威脅,在該偏轉(zhuǎn)角下回流區(qū)相對(duì)較小,高溫流體尚未流至噴嘴下方;β=55°時(shí)偏流板的導(dǎo)流效果較好,板后高溫區(qū)域較小,但噴嘴下方回流區(qū)域的溫度較高,回流的高溫流體有可能灼傷飛機(jī)后機(jī)輪;β=45°時(shí)偏流板導(dǎo)流效果較好,且回流區(qū)域的高溫流體尚未影響到噴嘴。
圖 5 不同偏轉(zhuǎn)角下偏流板附近的溫度云圖Fig.5 Temperature contours near the deflector under different deflection angles
表1給出了3種偏轉(zhuǎn)角條件下尾噴管下方區(qū)域以及偏流板后方區(qū)域測(cè)量點(diǎn)A,B的具體溫度值(測(cè)點(diǎn)布置詳見(jiàn)圖2(a))。從表1中可以看出,當(dāng)β=45°,55°時(shí),偏流板后方測(cè)點(diǎn)B的溫度已接近環(huán)境溫度,不會(huì)對(duì)板后人員和設(shè)備造成威脅;β=45°時(shí),尾噴管下方測(cè)點(diǎn)A的溫度為415.87 K,低于后輪機(jī)可承受的溫度423.15 K。因此,無(wú)論是考慮偏流板后方的安全性,還是飛機(jī)后輪機(jī)的耐高溫能力,β=45°都是最佳選擇。
表 1 不同偏轉(zhuǎn)角條件下流場(chǎng)內(nèi)測(cè)量點(diǎn)的靜溫值Table 1 Static temperature at measurement points in the flow field under different deflection angles
沖擊射流的湍流流場(chǎng)對(duì)溫度場(chǎng)和噪聲輻射均有較大影響。Q準(zhǔn)則判據(jù)是一種直接觀(guān)察湍流運(yùn)動(dòng)規(guī)律的有效手段。為此,圖6展示了3種偏轉(zhuǎn)角條件下從噴嘴到偏流板區(qū)域渦量的Q準(zhǔn)則分布。從圖6中可以看出,射流與偏流板的撞擊強(qiáng)度是隨偏轉(zhuǎn)角的增大而增強(qiáng)的,這使得回流區(qū)的面積和渦流強(qiáng)度也隨偏轉(zhuǎn)角的增大而加強(qiáng)。因此,測(cè)點(diǎn)A的溫度會(huì)隨偏轉(zhuǎn)角的增大而升高(表1)。同時(shí),當(dāng)射流撞擊35°偏轉(zhuǎn)角的偏流板時(shí),板后產(chǎn)生的渦結(jié)構(gòu)明顯多于另外2種偏轉(zhuǎn)角情況,因此,β=35°時(shí)板后面的溫度明顯高于β=45°和55°情況(表1)。
圖 6 射流漩渦分布圖Fig.6 Jet vortex distribution
另一方面,由射流沖擊引發(fā)的漩渦總量、剪切層區(qū)域的渦量、撞擊區(qū)域渦量的集中程度都將隨偏轉(zhuǎn)角的增大而增加,這就預(yù)示著向空間輻射的沖擊噪聲也將隨著偏流板的增加而增大[14]。
高馬赫數(shù)時(shí),四極子噪聲的影響不可忽略。為了研究高馬赫數(shù)條件下四極子聲源的噪聲特性及其貢獻(xiàn)大小,選取了2種源面形式進(jìn)行聲學(xué)計(jì)算:一種是同時(shí)選取聲積分面內(nèi)部壁面及聲積分面(包含四極子噪聲源貢獻(xiàn));另一種只選取聲積分面內(nèi)部壁面(不包含四極子噪聲源貢獻(xiàn))。
如文獻(xiàn)[15]所述,射流噪聲的測(cè)量可以分為近場(chǎng)聲測(cè)量和遠(yuǎn)場(chǎng)聲測(cè)量2種。Viswanathan[16]通過(guò)對(duì)比由不同距離測(cè)量得到的噪聲的頻譜譜形發(fā)現(xiàn),遠(yuǎn)場(chǎng)測(cè)量距離應(yīng)至少為35倍噴口直徑才能滿(mǎn)足射流頻譜譜形的相似性要求。因此,圖7和圖8分別描述了不同偏轉(zhuǎn)角條件下,與尾噴管距離40D監(jiān)測(cè)點(diǎn)處,包含四極子聲源和不包含四極子聲源的1/3倍頻程特性(圖中,SPL為聲壓級(jí))。通過(guò)對(duì)比圖7和圖8可以看出,四極子聲源主要表現(xiàn)為低頻段噪聲,因?yàn)閳D7和圖8中的中、高頻段噪聲幾乎沒(méi)有差別,而兩圖中的低頻段噪聲則相差較大。而且,低頻段內(nèi)噪聲是隨偏轉(zhuǎn)角度的增大而增加的,這是因?yàn)樯淞骷羟袑又饕a(chǎn)生帶有低頻特性的四極子噪聲,而射流剪切層區(qū)域渦量則是隨斜板角度的增大而增加的(圖6)。因此,在射流沖擊斜板的過(guò)程中,產(chǎn)生于剪切層的四極子噪聲不容忽視[13]。
圖 7 不同偏轉(zhuǎn)角下的1/3倍頻程(偶極子聲源+四極子聲源)Fig.7 The 1/3 octave of different deflection angle (dipole sound source+quadrupole sound source)
圖 8 不同偏轉(zhuǎn)角下的1/3倍頻程(偶極子聲源)Fig.8 The 1/3 octave of different deflection angle (dipole sound source)
圖9和圖10分別為3種偏轉(zhuǎn)角下,與尾噴管距離40D的監(jiān)測(cè)點(diǎn)在不同監(jiān)測(cè)角θ處,包含四極子聲源和不包含四極子聲源的總聲壓級(jí)(OASPL)指向性圖。從圖中可以看出,不管是純粹的偶極子噪聲,還是偶極子和四極子噪聲的總和,都是隨偏轉(zhuǎn)角的增大而增加的,這與圖6的預(yù)測(cè)相符[17]。從圖9和圖10中還可以看出,監(jiān)測(cè)角θ越大,越靠近沖擊板面,聲壓級(jí)也越大;在圖10中該現(xiàn)象表現(xiàn)得比圖9更明顯,說(shuō)明沖擊板面噪聲主要為偶極子聲源。圖9中結(jié)果表明,在不同偏轉(zhuǎn)角下,射流沖擊噪聲的指向性具有相同的方向性,當(dāng)監(jiān)測(cè)角θ=60°時(shí),噪聲最小。
研究表明,當(dāng)環(huán)境噪聲超過(guò)80 dB時(shí),人們會(huì)感覺(jué)心煩意亂,無(wú)法正常工作;而長(zhǎng)時(shí)間處于噪聲超過(guò)90 dB的環(huán)境中,就會(huì)損壞聽(tīng)力[18]。因此,有必要分析不同偏轉(zhuǎn)角下的安全工作區(qū)域范圍。
圖11描述了不同偏轉(zhuǎn)角情況下的80 dB等總聲壓線(xiàn),等值線(xiàn)右側(cè)區(qū)域?yàn)樵肼暤陀?0 dB的安全區(qū)域。由圖可知,β=35°時(shí)安全區(qū)域最大,距噴嘴中心3 m以外的區(qū)域噪聲輻射均小于80 dB;β=45°時(shí) 的 安 全 區(qū) 域 略 小 于β=35°情 況;而當(dāng)β=55°時(shí),安全區(qū)域急劇縮小,且不同監(jiān)測(cè)角下的不安全半徑變化較大。
圖 9 總聲壓級(jí)指向性(偶極子聲源+四極子聲源)Fig.9 Directivity of total sound pressure level (dipole sound source+quadrupole sound source)
圖 10 總聲壓級(jí)指向性(偶極子聲源)Fig.10 Directivity of total sound pressure level (dipole sound source)
圖 11 不同偏轉(zhuǎn)角下80 dB等總聲壓線(xiàn)Fig.11 80 dB isoline of OASPL under different deflection angels
綜合考慮不同偏轉(zhuǎn)角下的溫度及噪聲分布情況發(fā)現(xiàn),β=45°時(shí)偏流板附近人員及設(shè)備所受的高溫及噪聲威脅最小。因此,下面將進(jìn)一步分析45°偏轉(zhuǎn)角下的安全范圍。圖12所示為β=45°時(shí),不同監(jiān)測(cè)距離 (S=20D,30D,50D,80D)下總聲壓級(jí)的大小及指向性。從圖12中可以看出,當(dāng)β=45°時(shí),與尾噴管保持30D及以上的距離方能保證所有角度下的噪聲輻射都在90 dB以下。而為了保證人員在任何角度下都能正常工作,需與尾噴管保持80D以上的距離。
圖 12 不同監(jiān)測(cè)距離總聲壓級(jí)指向性Fig.12 Directivity of total sound pressure level at different monitoring distances
本文采用LES結(jié)合FW-H聲類(lèi)比方法研究了不同偏轉(zhuǎn)角下艦載機(jī)發(fā)動(dòng)機(jī)沖擊射流流場(chǎng)、溫度場(chǎng)及聲場(chǎng)的特性,分析了艦載機(jī)的安全區(qū)域范圍,主要得出以下結(jié)論:
1) 由不同偏轉(zhuǎn)角的溫度場(chǎng)對(duì)比分析可知,當(dāng)β=55°時(shí)高溫回流氣體會(huì)灼傷飛機(jī)后機(jī)輪,而β=35°時(shí)板后的高溫氣體會(huì)威脅到人員及設(shè)備的安全。因此,無(wú)論是考慮偏流板后方的安全性還是飛機(jī)后機(jī)輪的耐高溫能力,β=45°都是最佳選擇。
2) 通過(guò)采用不同的聲積分面來(lái)區(qū)分偶極子和四極子聲源發(fā)現(xiàn),射流沖擊斜板時(shí),產(chǎn)生于沖擊板面的偶極子聲源具有高頻特性;產(chǎn)生于剪切層的四極子聲源具有低頻特性,且對(duì)總聲壓級(jí)的貢獻(xiàn)較大,四極子噪聲不容忽視。
3) 從噪聲輻射的角度,無(wú)論是總聲壓級(jí)還是安全區(qū)域,β=35°下的控制效果都是最好的,此時(shí),同一監(jiān)測(cè)角度下的總聲壓級(jí)最小,安全區(qū)域最大。同時(shí),不同偏轉(zhuǎn)角下的噪聲指向性均具有相同的方向性,在與尾噴管距離40D監(jiān)測(cè)點(diǎn)處,60°監(jiān)測(cè)角下的噪聲最小。
4) 綜合考慮溫度場(chǎng)及噪聲輻射,發(fā)現(xiàn)β=45°時(shí)效果最佳。在該角度下,與尾噴管保持30D及以上距離方能保證所有角度下的噪聲輻射都在90 dB以下。而為了保證人員能正常工作,需與尾噴管保持80D以上的距離。