定子預(yù)旋角對泵噴推進器敞水脈動特性影響

李福正,黃橋高,潘 光,秦登輝

(1.西北工業(yè)大學(xué) 航海學(xué)院,西安 710072；2.無人水下運載技術(shù)工業(yè)和信息化部重點實驗室(西北工業(yè)大學(xué)),西安 710072)

泵噴推進器(Pumpjet,簡稱泵噴)是包含定子、轉(zhuǎn)子和導(dǎo)管的一種組合式推進器,一般以定子相對轉(zhuǎn)子前后順序分為前置定子式和后置定子式[1]。其中前置泵噴雖然相對后置雖然效率低,但由于定子工作在相對壓力波動較低的區(qū)域,避免了轉(zhuǎn)子高速尾流的沖擊,可極大減小推進器噪聲[2],因此目前被廣泛應(yīng)用于核潛艇,如以英國、法國、美國、俄羅斯分別為代表的最新研制的“無畏級”、“哥倫比亞級”、“梭魚級”和“北風(fēng)之神級”核潛艇均采用前置泵噴。

目前前置泵噴水動力性能獲取主要途徑為數(shù)值模擬和實驗方法,由于數(shù)值方法成本相對低,并能最大程度保留流場信息,是泵噴設(shè)計階段和校核階段普遍采用的方法。已有文獻如Huyer等[3]將定子安裝角度沿周向按照正弦函數(shù)排列,由于單個轉(zhuǎn)子進流條件不同,最終產(chǎn)生了側(cè)向力,可作為新型矢量推進器。Park等[4]采用低雷諾數(shù)k-ε方程研究了帶前置定子的導(dǎo)管槳,主要分析了槳葉表面流線和不同半徑處的壓力系數(shù)。Ahn等[5]以單方程S-A模型分析了轉(zhuǎn)子葉梢加葉冠后的水動力性能,結(jié)果表明雖然葉冠可有效控制葉梢泄露渦,但卻造成推力損失,而效率保持不變。趙威等[6]研究了前置定子彎曲對槳葉前三階非定常推力的性能,結(jié)果表明前彎和后彎均可有效降低推力脈動,并分析了不同軸向位置處定子尾流沿周向分布,文中指出可能是由于定子尾跡相位差引起。Wang等[7]針對葉梢渦模型,并結(jié)合面元法分析了水動力性能。彭云龍等[8]從敞水特性、空化特性和脈動特性3方面比較了前置和后置泵噴的性能,結(jié)果表明前置泵噴效率較低,空化性能較差,但導(dǎo)管進出口壓力脈動遠低于后置泵噴。Qin等[9]同樣分析了前置和后置泵噴區(qū)別,除非定常脈動性能外,重點分析了泵噴尾渦動力特性。結(jié)果表明前置泵噴轂渦相對穩(wěn)定,梢渦同樣受導(dǎo)管和葉梢相互作用而增強。由于前置泵噴特殊應(yīng)用,公開實驗研究較少,如Yu等[10]通過空泡水筒的實驗驗證了前置泵噴計算方法的準確性。

雖然前置泵噴在民用方面已有普及,如美國通用電氣公司(General electric Company,GE)和ABB分別研發(fā)的吊艙推進器[11],但目前仍以潛艇應(yīng)用居多,故噪聲指標極為重要。泵噴噪聲與流場低頻脈動息息相關(guān),改善泵噴噪聲的重點在于如何調(diào)整推進器的結(jié)構(gòu),以促使流場向更有利于聲隱性的方向發(fā)展。而前置定子結(jié)構(gòu)是泵噴水動力性能和流場特性的決定性因素之一,如饒志強等[12]采用RNGk-ε模型研究了定子葉片數(shù)、翼型厚度及安裝角對泵噴定常性能的影響,表明安裝角對推力、扭矩和效率影響最大。Yu等[13]采用SSTk-ω分析了定子弦長、安裝角、轉(zhuǎn)定距離對泵噴的非定常性能的影響,結(jié)果同樣表明安裝角對槳葉表面脈動特性和推力脈動影響最大?；谏鲜鲅芯?本文將以某泵噴為原型,采用Ansys/Fluent重點研究前置定子角度對包括槳葉在內(nèi)的各部件及轉(zhuǎn)定之間流場點的低頻脈動性能影響。不同于其他學(xué)者,本文將定子角度分為進流角、出流角和安裝角,并對其分別進行調(diào)整研究。

1 泵噴模型

1.1 幾何參數(shù)

本文采用的泵噴,其原始幾何模型(origin model)如圖1所示,定子和轉(zhuǎn)子葉片隨邊為鈍邊,葉根部分別有2.0,2.5 mm半徑的倒圓,為后續(xù)幾何模型處理及離散網(wǎng)格的方便,本文研究的模型取消了葉根倒圓,并且將葉片隨邊處理為圓角。改型后的泵噴(modified model)如圖1所示,其中定子和轉(zhuǎn)子葉片數(shù)目分別為8和6,轉(zhuǎn)子的直徑(Dr)為166.4 mm,間隙為1 mm,大約為0.6%Dr,轂徑比為0.3,投影盤面比Ap/A0=0.8。轉(zhuǎn)子各半徑處螺距比(P/D)見表1。圖2為定子和導(dǎo)管的主要幾何參數(shù)示意圖,具體參數(shù)見表2。由于泵噴進流段與潛艇尾端相匹配,故定子槳轂為圓錐面,導(dǎo)致前后定子直徑不同。圖2中Ds1和Ds2分別為定子進口和出口直徑,H為定子葉高,表示葉片軸向長度。Dd1、Dd2分別為導(dǎo)管進口和出口直徑,Ld為導(dǎo)管軸向總長,坐標原點位于轉(zhuǎn)子中間正對下方r=0的位置。

圖1 前置泵噴外觀

表1 轉(zhuǎn)子各半徑處螺距比

圖2 泵噴模型幾何參數(shù)

表2 定子和導(dǎo)管參數(shù)

1.2 定子角度調(diào)節(jié)

本文泵噴定子采用圓弧設(shè)計法,即骨線(中弧線)為圓弧段,通過控制沿進流和出流方向切線角度以確定圓弧形狀,再將其根據(jù)翼型沿弦長的厚度分布進行加厚。圖3(a)所示為定子在徑向為0.7展面時剖面翼型示意圖。其中進流角和出流角的幾何關(guān)系為αin+αout=90°+β。

在調(diào)節(jié)定子出流角和進流角時,控制另一個角度不變,通過幾何關(guān)系先求得圓弧半徑R,得到骨線。再根據(jù)原始翼型厚度分布沿半徑方向加厚,加厚方法為：以骨線某一點為圓心,以對應(yīng)位置厚度為直徑作圓,與半徑方向相交的點即為壓力面和吸力面的坐標點,由此得到新調(diào)整的翼型。在調(diào)節(jié)安裝角αsta時,則直接將翼型繞前緣點旋轉(zhuǎn)對應(yīng)的角度。定子角度調(diào)節(jié)為

(1)

此外角度的調(diào)整將引起定子投影面軸向長度不一致,為消除定子尾緣到轉(zhuǎn)子距離對壓力脈動的影響,將定子不同剖面處翼型等比例擴展到與原模型軸向長度一致。圖3(b)為采用不同預(yù)旋角調(diào)整方式后的翼型示意圖。

圖3 定子剖面翼型幾何示意

2 數(shù)值分析方法與驗證

2.1 控制方程

采用基于雷諾平均的N-S控制方程：

(2)

(3)

2.2 計算域及邊界條件設(shè)定

將求解區(qū)域劃分為內(nèi)流域和外流域,其中內(nèi)流域又可分為定子域和轉(zhuǎn)子域,考慮到計算成本及后處理時高保真度,計算域全部采用結(jié)構(gòu)化網(wǎng)格。如圖4、5分別展示了外流場和內(nèi)流程的網(wǎng)格,為避免進口及出口區(qū)域?qū)α鲌龅母蓴_,外流域的前后距離分別為5Dr和10Dr,在徑向處距離同樣設(shè)為5Dr,使流場充分延伸,具體邊界條件如圖4所示。在網(wǎng)格劃分過程中充分考慮到近壁面影響,通過布拉休斯公式,即

(4)

計算得到y(tǒng)+為1時,壁面第1層網(wǎng)格y應(yīng)為2.6×10-6m,其中Vref為特征速度,取槳葉在0.75半徑處剖面翼型的進流速度,ν為流體運動黏性系數(shù),即μ/ρ。Lref為特征長度,取轉(zhuǎn)子0.75半徑處翼型弦長,即128 mm。為得到精確的結(jié)果,本文采用y=2.0×10-6m,保證壁面處有足夠密的附面層網(wǎng)格。對于定常工況,基于多重參考系(MRF)模型迭代不小于1 000步,保證結(jié)果收斂,并提取最后一步作為泵噴敞水特性結(jié)果;而非定常工況則基于Sliding Mesh技術(shù)在定常流場的基礎(chǔ)上求解,并以槳葉轉(zhuǎn)過1°的時間為時間步長,共計算6個周期,以第4個周期以后結(jié)果作為有效數(shù)據(jù)。

圖4 計算域及邊界條件

圖5 內(nèi)流域結(jié)構(gòu)化網(wǎng)格

2.3 網(wǎng)格驗證

泵噴完全沉入深水中,且不受自由液面的影響,在邊界條件設(shè)置中出口處相對壓力設(shè)為0,在分析流場形態(tài)時可加上水壓頭。固定螺旋槳轉(zhuǎn)速n=20 r/s,改變進速VA來模擬在水槽中的實驗狀態(tài)。進速與z軸正方向同向,轉(zhuǎn)速根據(jù)右手定則為z軸負方向。分析泵噴的性能時,其各部件包括轉(zhuǎn)子、定子和導(dǎo)管的推力系數(shù)和扭矩系數(shù)及敞水效率可分別為

(5)

式中：Tr、Ts、Td、Qs、Qs分別為槳葉推力、定子推力、導(dǎo)管推力、槳葉扭矩、定子扭矩,總推力為T,η0為敞水效率,n為螺旋槳轉(zhuǎn)速,ρ為水在常溫25°時的密度,取997 kg/m3。

為保證計算結(jié)果的準確度,對比了不同數(shù)目網(wǎng)格的計算結(jié)果,表3所示為3種不同內(nèi)流場網(wǎng)格單個計算通道數(shù)目的對比。而外域網(wǎng)格一般只進行微調(diào),除保證導(dǎo)管近壁面足夠多的網(wǎng)格外,針對與之相匹配的轉(zhuǎn)子域和定子域的網(wǎng)格大小,網(wǎng)格要保證較好地連續(xù)過渡。最終外域網(wǎng)格數(shù)目約為259萬。

表3 不同數(shù)目網(wǎng)格

圖6所示為泵噴在n=20 r/s水動力性能曲線,圖6(a)為不同網(wǎng)格數(shù)目的槳葉推力系數(shù)和扭矩系數(shù)同實驗數(shù)據(jù)[14]的對比,可以看出,網(wǎng)格的細化對結(jié)果并不會造成太大的影響,對比中等網(wǎng)格與實驗結(jié)果,發(fā)現(xiàn)在進速系數(shù)J=0.1～0.8時,推力系數(shù)誤差為5%～7%,扭矩系數(shù)為2%～4%,而誤差值隨著進速系數(shù)的提高逐步呈非線性增加,在進速系數(shù)

圖6 泵噴水動力系數(shù)曲線

為1.2時誤差最大,推力系數(shù)和扭矩系數(shù)的誤差分別達到10%和8%。由于是非設(shè)計工況,故不予考慮?？紤]到計算成本,本文后續(xù)將采用中等網(wǎng)格做進一步研究。圖6(b)所示為泵噴各部件及總體性能,注意其中定子推力系數(shù)為負,其扭矩系數(shù)方向與槳葉相反,為方便后續(xù)分析,取其絕對值。導(dǎo)管在不同進速下對泵噴總推力所體現(xiàn)的作用不同,可以看到在進速系數(shù)約為0.7時導(dǎo)管推力由正轉(zhuǎn)負,并隨進速的增加逐步增加。對比不同進速系數(shù)下由總推力系數(shù)和槳葉扭矩系數(shù)計算得到的敞水效率η0,可以看出相比傳統(tǒng)螺旋槳,η0在較廣的進速系數(shù)范圍內(nèi)變化并不大,這也是泵噴優(yōu)于傳統(tǒng)螺旋槳的原因之一。本文泵噴最高效率大約為J=0.9,考慮到一般艦艇設(shè)計工況低于最高效率點,故本文后續(xù)將以J=0.8為主要分析工況。

圖7所示為泵噴在J為0.8時近壁面y+分布云圖,其中槳葉和定子從左向右分別為壓力面和吸力面,可以看到槳葉導(dǎo)邊靠近葉梢處y+值最大,但不超過1,而其他部位都在0.5左右,對于定子和導(dǎo)管,由于其特征速度較低,采用與槳葉相同的第一層附面層網(wǎng)格導(dǎo)致y+值遠小于1?？傮w來看,滿足后續(xù)分析的精度要求。

圖7 壁面y+分布云圖

3 結(jié)果分析

3.1 泵噴敞水特性

圖8所示為泵噴在不同定子角度時敞水特性曲線,可以看出,隨著泵噴的推力、扭矩和最高效率都隨著定子角度的增大而增大。其中αsta的影響最大,且最高效率點隨αsta的增大而向較高進速系數(shù)偏移,αout影響次之,αin影響最小。綜合分析認為：定子角度的增大將使流體預(yù)旋更強烈,即增加了槳葉進流位置周向速度,使槳葉各半徑翼型剖面攻角增大,與螺旋槳中增大槳葉螺距角作用相同。由于αout改變的是出流位置,更接近于槳葉進流位置,而調(diào)整αin雖然增大了定子首部預(yù)旋角,但流體卻在流過定子狹長的通道后導(dǎo)致出流位置的預(yù)旋損失,因此αout影響比αin更顯著;由于調(diào)整αsta是同時改變αin和αout,αsta每增大1.0°,αin和αout都要增大0.5°,因此αsta的效果更明顯。此外,雖然角度的增大會導(dǎo)致定子阻力增大,但由于定子負載相對轉(zhuǎn)子小的多,其絕對增大值遠小于轉(zhuǎn)子,故使得總推力、扭矩和效率均呈現(xiàn)增大趨勢。

圖8 不同定子角度泵噴敞水曲線

3.2 槳葉推力脈動

槳葉非定常力可以從一定程度上反應(yīng)葉片表面壓力脈動,并常用來分析葉片激振力是否滿足設(shè)計要求。圖9所示為J=0.8時,不同定子角度下轉(zhuǎn)子單葉片非定常推力在頻域的結(jié)果,圖9(a)為推力的頻譜曲線,縱坐標表示推力在對應(yīng)頻率下的脈動幅度?？梢园l(fā)現(xiàn),受定子干擾,轉(zhuǎn)子單葉片的推力脈動峰值對應(yīng)頻率主要為BPFs(定子數(shù)和轉(zhuǎn)速的乘積)的倍頻程。此外還存在大量離散低幅軸頻fa(即轉(zhuǎn)速)的倍頻程信息,并在4fa時有相對較大的脈動力。圖9(b)為進一步對比包括4fa和1BPFs～6BPFs特征頻率處推力脈動具體幅值大小,結(jié)果表明：1BPFs脈動幅度相對其他頻率最大,以αφ=0°為例,2BPFs下降了80%,其他頻率下降超過90%。在1BPFs處,推力脈動幅度最低值在αin和αsta為4°時分別達到最小;而至于αout,脈動幅度與之呈單調(diào)性減小趨勢,并且該趨勢隨αout的增大逐漸變緩;對比同一αφ下不同αφ調(diào)整方式的脈動幅度,可以發(fā)現(xiàn)當(dāng)αφ為正時,αsta影響最大,αin次之,αout最小;當(dāng)αφ為負時,影響規(guī)律不是特別明顯。在2BPFs處,脈動幅度分別在αin=8°、αout=0°和αsta=4°達到最小。在3BPFs處,最小值分別為αin=-8°、αout=-4°和αsta=-8°。其他頻率處脈動幅度最小值對應(yīng)的角度也均存在差異。

圖9 不同角度時單葉片推力頻域結(jié)果

圖10所示為J=0.8時,不同定子角度下槳葉總推力的非定常脈動幅度在頻域的結(jié)果,圖10(a)為對應(yīng)的頻譜曲線,可以發(fā)現(xiàn),受轉(zhuǎn)定子數(shù)目干擾,最大脈動幅度對應(yīng)頻率為4BPFr(轉(zhuǎn)子數(shù)目與轉(zhuǎn)速的乘積),同時也是3BPFs。根據(jù)文獻[15]的推導(dǎo),在對轉(zhuǎn)螺旋槳中,一般非定常推力和扭矩的變化頻率為：f=fa(mBF+kBA),其中BF和BA分別為前槳和后槳葉片數(shù),m和k為常數(shù),并滿足mBF=kBA。由于泵噴中定子為靜止?fàn)顟B(tài),相對于轉(zhuǎn)子以n的轉(zhuǎn)速反向旋轉(zhuǎn)。故可看作轉(zhuǎn)速為n/2對轉(zhuǎn)推進器,其特征頻率為：f=fa(mBF+kBA)/2,對應(yīng)本文的480 Hz。此外圖中波峰位置除數(shù)倍的fa外,在諸如1BPFr～5BPFr處均存在不同幅度的脈動峰值。進一步比較數(shù)倍于BPFr的特征頻率下脈動幅值,可以看出4BPFr處推力幅值最大,并遠大于其他頻率。以αφ=0°為例,1BPFr～5BPFr相對其下降幅度超過95%。當(dāng)角度為正時,尤其是αφ=8°時,αsta對脈動幅度的影響最大,αin次之,αout最小;而角度為負時,3種角度的調(diào)整方式影響不大。在4BPFr處,脈動幅度隨αin和αsta的增大而增大,而αout=-4°時脈動幅度最小。其他頻率的脈動幅度隨角度變化規(guī)律并不一致,最小幅值對應(yīng)的角度也各不相同。

圖10 不同角度時槳葉總推力頻域結(jié)果

3.3 葉片表面壓力脈動

葉片表面的壓力脈動是推進器重要的物理量之一,不僅反映水動力噪聲強度,而且是非定常激振力的成因之一。本文主要分析定子和轉(zhuǎn)子葉片表面的壓力脈動,并以式(5)中的壓力標準差來衡量。

(6)

圖11為J=0.8時槳葉表面壓力標準差云圖,其中圖11(a)為槳葉壓力面,圖11(b)為槳葉吸力面,從左到右分別對應(yīng)αφ為-8°,-4°,0°,4°和8°的云圖,從上到下分別為不同αin、αout、αsta的影響,流向從左向右。

圖11 不同角度時槳葉壓力標準差

可以看到,槳葉壓力面和吸力面均存在較強的壓力脈動。葉片前緣因與流體首先相接觸而受到?jīng)_擊,其壓力脈動最強,并且最大值超過1 000 Pa,隨著角度增大,前緣局部較大的壓力標準差范圍逐漸增大。葉梢附近由于受葉梢泄露渦(Tip leakage vortex，TLV)、葉梢泄露流(Tip leakage flow，TLF)[16]影響,在壓力面和吸力面均存在較大脈動壓力,定子角度的增大均使得葉梢壓力脈動增強;然而3種角度的調(diào)整方式對葉梢脈動影響規(guī)律在αφ為正、負時不同,在同一角度下,當(dāng)αφ為負時,脈動強度規(guī)律為：αin>αout>αsta,而當(dāng)αφ為正時,αin<αout<αsta。此外,在經(jīng)過葉片前緣較高的脈動強度區(qū)域后,靠近前緣部分范圍內(nèi)及葉根處均存在局部較高的脈動幅值;而葉片徑向中部存在一定壓力波動平穩(wěn)區(qū)域。整體來看,隨著定子角度增大,該區(qū)域范圍逐漸減小,尤其當(dāng)αsta增大到8°時,脈動較平穩(wěn)區(qū)域大幅減小,并且在αφ=8°時,該區(qū)域范圍大小規(guī)律為αsta>αin>αout,同時也印證了槳葉推力脈動中定子角度調(diào)整方式對單槳葉推力脈動幅度的影響規(guī)律。而對于其他αφ,3種角度調(diào)整方式對該區(qū)域范圍大小的影響規(guī)律并不完全一致。

圖12所示為J=0.8時不同定子葉片表面的壓力標準差云圖,該圖從左到右及從上到下布置對應(yīng)的角度大小及調(diào)整方式同圖11不同的是,每個角度對應(yīng)兩個云圖,左右分別為定子壓力面和吸力面。流向從左向右。

圖12 不同角度時定子壓力標準差

由于定子圍流場變化不是特別劇烈,其負載相對轉(zhuǎn)子小得多,故脈動幅度也較低。而定子葉片在較高徑向位置及后緣區(qū)域存在較大分層顯示的壓力標準差,并且越靠近后緣,其脈動幅度越大;但壓力面高脈動區(qū)范圍小于吸力面。分析原因可能為葉片吸力面邊界層出現(xiàn)分離,導(dǎo)致定子尾緣存在周期性脫落渦,進而使定子隨邊壓力脈動較大。

此外,當(dāng)αφ為正時,隨著角度增大,吸力面前緣位置壓力波動幅度及范圍均劇烈增加,所不同的是,調(diào)整αsta和αin均比αout影響大,尤其當(dāng)αsta=8°時,吸力面高壓力波動范圍區(qū)域大幅增加。而對于同一角度下,3種角度調(diào)整方式對葉表脈動強度的影響規(guī)律有所不同,如在壓力面,當(dāng)αφ為負時,葉片尾緣局部脈動強度規(guī)律為αin>αout>αsta,當(dāng)αφ為正時,αin<αout<αsta;而吸力面規(guī)律在αφ為負時并不明顯。

3.4 泵噴內(nèi)部流場點壓力脈動

為進一步分析定子角度對泵噴流場壓力脈動的影響,在非定常計算中引入多個監(jiān)控點,監(jiān)測位置如圖13所示,在包括定子前緣、轉(zhuǎn)定子之間及噴尾部不同軸向位置(A1～A5)和槳轂到槳罩之間的徑向位置(S1～S5)分別布置壓力監(jiān)控點。其中軸向坐標分別為z=-0.09,-0.06,-0.03,0.03,0.06;徑向截面位置分別為Span=0.1,0.2,0.3,0.4,0.5。監(jiān)測點以軸向和徑向位置命名,如圖13中的P點為：A1S3,其壓力頻譜曲線如圖14所示?？梢钥闯?壓力脈動最大幅值位于1BPFr處,2BPFr～4BPFr對應(yīng)的脈動幅值已經(jīng)下降了95%～98%。因此本文將主要分析各監(jiān)測點在1BPFr處的脈動幅值。

圖13 壓力監(jiān)控點布置

圖14 監(jiān)測點P(A1S3)的頻譜曲線

圖15所示為原模型,即αφ=0°時各監(jiān)測點經(jīng)傅里葉變換后在120 Hz的壓力信息,圖15(a)為所有監(jiān)測點的壓力幅度云圖,可以看出,在徑向方向,壓力幅度由槳轂到槳罩呈先增大后減小的規(guī)律,這是由于槳葉在位于中間及較高半徑處的葉元體為主要做功位置所致;在軸向方向,呈現(xiàn)從上游到下游逐步增大的趨勢,并且最大位置出現(xiàn)在泵噴尾部,這是由于泵噴尾渦結(jié)構(gòu)不穩(wěn)定作用所引起。圖15(b)為不同軸向位置處,展向S2～S5截面監(jiān)測點相對于S1處的相位差。由于葉表及附近壓力脈動與非定常流動直接相關(guān),而各展向截面處的流動表達式為

m(r,t)=m(r)+m′[r,ωt+φm(r)]

(6)

式中m為流場中物理量,可表示為流量、速度、壓力等,其中m為平均值,m′為脈動量,并且在不同半徑r處以同頻角速度ω,不同角速度φm呈正弦脈動。因此各監(jiān)測點相位信息也是泵噴低頻脈動的關(guān)鍵信息?？梢钥闯鲈诎?.5及以下的徑向位置處,各點相較于S1處相位均均為延遲;而大于0.5展向后均有延遲或超前,尤其軸向A2和A5位置,在0.9展向處相位變化劇烈,這是由于在葉頂端壁黏性渦漩及泵噴尾流復(fù)雜的渦流所致。

圖16所示為J=0.8時,所有監(jiān)測點在120 Hz的壓力幅度云圖,云圖色階及范圍同圖15(a)。比較不同定子角度調(diào)整方式對最大位置及附近壓力脈動幅度的影響,可以發(fā)現(xiàn)其與定子角度呈正相關(guān);當(dāng)角度為負時,角度調(diào)整方式對其影響規(guī)律為：αin>αout>αsta;而角度為正時,該處的壓力脈動幅度影響規(guī)律為：αin<αout<αsta。

圖15 各監(jiān)控點在1BPFr處壓力信息

圖16 監(jiān)測點在1BPFr壓力幅值云圖

槳葉入流條件是決定葉片表面非定常脈動特性的因素之一,因此需重點分析槳葉前緣的流動。圖17給出了A3位置不同徑向處120 Hz的相位關(guān)系,其中相位差計算時以0.1展向截面處即S1位置為基準?？梢钥闯?不同定子角度相位差在調(diào)整αin時變化相對較小,并且相位與原模型一致,均落后于基準點。而αout和αsta均以αφ的正、負為區(qū)別表現(xiàn)出兩極化的趨勢,其中αsta表現(xiàn)更為明顯,即當(dāng)αφ為正時,相位均落后基準點;當(dāng)αφ為負時,相位提前于基準點;并且角度絕對值越大,其相位絕對值也越大。由于較高半徑處(S5位置)監(jiān)測點靠近導(dǎo)管黏性壁面,而該區(qū)域以壁面黏渦和槳葉附近復(fù)雜二次流動主導(dǎo),其相位關(guān)系與其他監(jiān)測點相比差異較大。

圖17 A3位置監(jiān)測點相位差

4 結(jié) 論

1)數(shù)值分析了3種定子角度調(diào)整方式對泵噴性能影響,發(fā)現(xiàn)增大進流角,出流角和安裝角均能增大泵噴推力、扭矩及最高效率,但安裝角影響最大,出流角次之,進流角最小。

2)分析J為0.8時,不同定子角度下推力、葉片表面、各監(jiān)測點脈動幅度。發(fā)現(xiàn):①單槳葉受定子葉片數(shù)影響,在1倍定子葉頻處達到最大,并且在調(diào)整角度為正時,對幅度影響規(guī)律為：安裝角>進流角>出流角,影響規(guī)律不明顯。而總推力脈動幅度受定子和轉(zhuǎn)子葉片數(shù)影響,在4倍轉(zhuǎn)子葉頻達到最大,其受定子角度調(diào)整方式影響規(guī)律與單槳葉相同。②槳葉前緣和葉梢附近脈動強度最大,而定子則在尾緣最強,并且兩者都隨定子角度增大而增強,但角度調(diào)整方式對其影響規(guī)律在調(diào)整角分別為正、負時相反。③泵噴較強壓力脈動分布在徑向中部及定子尾部,最大位置在尾部中高半徑處,并隨角度增加而增加;角度調(diào)整方式對其影響規(guī)律同樣在角度為正和負時相反。槳葉入流位置不同監(jiān)測點相位關(guān)系在調(diào)整進流角時差異較小,但出流角和安裝角的調(diào)整使其按定子角度正、負呈兩極化分布。