螺旋槳梢渦不穩(wěn)定性機(jī)理與演化模型研究1)

王戀舟 吳鐵成 郭春雨

* (西南交通大學(xué)力學(xué)與工程學(xué)院,成都 611756)

? (中山大學(xué)海洋工程與技術(shù)學(xué)院,廣東珠海 519082)

** (哈爾濱工程大學(xué)船舶工程學(xué)院,哈爾濱 150001)

引言

目前,水面艦船和水下航行器主要采用螺旋槳推進(jìn)方式.螺旋槳的水動(dòng)力性能、空化以及噪聲性能等一直是船舶研究和設(shè)計(jì)建造人員非常關(guān)心的問題[1-3].螺旋槳尾流場(chǎng)的測(cè)量和精準(zhǔn)模擬不僅關(guān)系到螺旋槳性能的評(píng)估,同時(shí)也是螺旋槳型線優(yōu)化、槳前后節(jié)能裝置設(shè)計(jì)時(shí)需要考量的關(guān)鍵因素[4-5].截止到目前,船舶螺旋槳尾流場(chǎng)這一高度復(fù)雜的湍流流動(dòng)結(jié)構(gòu)中依然包含著許多沒有被完全揭示和理解的演化機(jī)理,如橢圓不穩(wěn)定性(短波不穩(wěn)定性)、梢渦的融合、蛙跳、蜿蜒以及崩潰等復(fù)雜現(xiàn)象背后的物理機(jī)制.更好地揭示它們的流動(dòng)行為和特性往往能夠幫助人們更準(zhǔn)確地理解或預(yù)報(bào)螺旋槳的相關(guān)特性,如梢渦融合以及短波不穩(wěn)定性和長(zhǎng)波不穩(wěn)定性引起的渦系統(tǒng)崩潰,都會(huì)導(dǎo)致流噪聲的增加.螺旋槳尾流的機(jī)理研究在工程應(yīng)用中扮演著重要的角色,并引起了足夠的重視,究其原因是這一基礎(chǔ)問題與螺旋槳水動(dòng)力性能、振動(dòng)、噪聲以及船舶結(jié)構(gòu)問題等相關(guān)[6-7].另外,從物理的觀點(diǎn)來看,螺旋槳尾流的特征也是造成能量損失的關(guān)鍵因素,更好地理解螺旋槳尾流的物理本質(zhì),對(duì)螺旋槳的優(yōu)化設(shè)計(jì)以及新概念螺旋槳的發(fā)展十分重要[8].

螺旋槳尾流場(chǎng)計(jì)算的黏性方法仍主要采用計(jì)算流體力學(xué)方法進(jìn)行[9-10],但采用計(jì)算流體力學(xué)方法進(jìn)行螺旋槳尾流非定常特性的研究時(shí),也面臨著不小的挑戰(zhàn).當(dāng)螺旋槳尾流從近場(chǎng)向遠(yuǎn)場(chǎng)演化過程中的不穩(wěn)定性問題成為關(guān)注的焦點(diǎn)時(shí),這一物理過程中精細(xì)渦結(jié)構(gòu)的捕捉往往對(duì)數(shù)值模擬的時(shí)空分辨率要求較高[11-12].

國(guó)外在螺旋槳尾流場(chǎng)方面的研究較為深入和廣泛,對(duì)不同工況下螺旋槳的尾流場(chǎng)演化機(jī)理進(jìn)行了揭示.強(qiáng)光林等[13]研究了湍流模型效應(yīng)對(duì)螺旋槳尾流模擬的影響,分別討論了雷諾平均方法和分離渦模擬方法在模擬尾渦中的限制和優(yōu)勢(shì).王巍等[14]基于分離渦模擬方法研究了斜流中的螺旋槳尾渦動(dòng)力學(xué)特性,對(duì)螺旋槳尾渦的不穩(wěn)定性、湍動(dòng)能進(jìn)行了詳細(xì)分析.Dubbioso 等[15]采用雷諾平均方法研究了不同載況條件、不同斜流角情況下的螺旋槳水動(dòng)力性能.明尼蘇達(dá)大學(xué)的Mahesh 團(tuán)隊(duì)[16-18]在螺旋槳尾流場(chǎng)數(shù)值計(jì)算發(fā)面也開展了富有成效的研究.Jang 和Mahesh[16]基于大渦模擬方法對(duì)反轉(zhuǎn)螺旋槳的尾流特征進(jìn)行了詳細(xì)分析,通過研究發(fā)現(xiàn),螺旋槳重載荷主要由于流動(dòng)分離而產(chǎn)生.Kumar 和Mahesh[17]利用大渦模擬方法對(duì)五葉螺旋槳設(shè)計(jì)工況下的尾流場(chǎng)穩(wěn)定性進(jìn)行了研究,模擬結(jié)果揭示了梢渦相互干擾的不穩(wěn)定性機(jī)理.Verma 等[18]基于大渦模擬方法對(duì)兩個(gè)進(jìn)速系數(shù)條件下,螺旋槳上游存在和不存在船體時(shí)的反轉(zhuǎn)螺旋槳尾流場(chǎng)進(jìn)行了數(shù)值分析,重點(diǎn)考察了不同工況下反轉(zhuǎn)螺旋槳側(cè)向力產(chǎn)生的機(jī)理.愛荷華大學(xué)的Chase 和Carrica[19]探討了分別采用雷諾平均、分離渦模擬、延遲分離渦模擬以及層流模型時(shí),E1619 螺旋槳尾流場(chǎng)的計(jì)算結(jié)果,結(jié)果表明網(wǎng)格加密對(duì)推力和扭矩的計(jì)算結(jié)果影響不大,但對(duì)尾流細(xì)節(jié)有著非常大的影響,且采用雷諾平均方法時(shí)耗散較為嚴(yán)重.

國(guó)內(nèi)在螺旋槳尾流動(dòng)力學(xué)特性方面的研究不多.Zhu[20]基于雷諾平均方法對(duì)空化流中螺旋槳尾流場(chǎng)與螺旋槳側(cè)斜角特征的關(guān)系進(jìn)行了詳細(xì)分析.朱志峰等[21]利用多相流計(jì)算流體力學(xué)理論對(duì)螺旋槳模型的空化流特征進(jìn)行了數(shù)值仿真,其中滑移網(wǎng)格技術(shù)被用來實(shí)現(xiàn)螺旋槳的旋轉(zhuǎn)運(yùn)動(dòng),給出了不同空化模型和湍流模型的模擬結(jié)果.胡健等[22]基于雷諾平均方法計(jì)算并分析了螺旋槳尾流中速度分量的變化規(guī)律以及不同載荷工況下的流場(chǎng)特征.

本文針對(duì)螺旋槳尾流不穩(wěn)定性機(jī)理問題,基于粒子圖像測(cè)速試驗(yàn)以及不同湍流模擬方法(大渦模擬、延遲分離渦模擬、無湍流模型模擬)為手段,進(jìn)行了不同進(jìn)速系數(shù)條件下螺旋槳梢渦的不穩(wěn)定性觸發(fā)機(jī)制研究.根據(jù)探索得到的螺旋槳尾流不穩(wěn)定性觸發(fā)機(jī)理提出了螺旋槳梢渦演化模型,對(duì)梢渦的遷徙過程進(jìn)行模擬,并預(yù)測(cè)梢渦合并的時(shí)間和位置,以期為基于流動(dòng)控制技術(shù)的螺旋槳減振降噪設(shè)計(jì)提供基礎(chǔ)性支撐.

1 幾何模型

本文的研究對(duì)象為七葉螺旋槳INSEAN E1658[23],E1658 螺旋槳的主要幾何參數(shù)如表1 所示.圖1 給出了E1658 螺旋槳的幾何模型示意圖和參考系統(tǒng),其中x軸為旋轉(zhuǎn)軸,指向下游,z軸指向上,y軸由右手定則確定.

表1 E1658 螺旋槳的幾何參數(shù)Table 1 Main parameters of E1658 propeller

圖1 E1658 螺旋槳Fig.1 E1658 propeller

2 試驗(yàn)設(shè)置

螺旋槳的尾流場(chǎng)測(cè)量試驗(yàn)在意大利國(guó)家研究委員會(huì)海洋工程研究所(CNR-INM)的空泡水筒中進(jìn)行,其測(cè)試段尺寸為0.6 m × 0.6 m × 2.6 m (高 × 寬 ×長(zhǎng)),總長(zhǎng)度為13.9 m,總高度為9.3 m.該設(shè)備中來流的湍流度在2%以內(nèi),平均軸向速度和垂向速度的不均勻度分別在1%以內(nèi)和3%以內(nèi).

本文中螺旋槳的尾流場(chǎng)測(cè)量采用2D-PIV 系統(tǒng)進(jìn)行測(cè)量[24-25],該粒子圖像測(cè)速測(cè)量系統(tǒng)的照明由雙腔Nd-Yag 激光器(Quantel Twins,每個(gè)脈沖200 mJ,重復(fù)頻率12.5 Hz,波長(zhǎng)532 nm)以及圓柱發(fā)散透鏡(?15 mm 焦距)和球面透鏡組合提供(5000 mm 焦距),激光器的片光通過安裝在測(cè)試截面處底部的窗戶向上照射,散射光是由多相機(jī)記錄系統(tǒng)采集的,該系統(tǒng)由LaVision 的4 臺(tái)Imager sCMOS 相機(jī)組成,該相機(jī)能達(dá)到每秒50 幀的幀速率,其分辨率為2560 ×2160 像素.每個(gè)相機(jī)都配備了尼康的Nikkor AF-S 85 mm f/1.4G 鏡頭,其光圈設(shè)置為f/8.照相機(jī)在光軸垂直于激光片的縱向上并排布置,該試驗(yàn)設(shè)置方案已經(jīng)在類似的實(shí)驗(yàn)研究中采用[26-27].

粒子圖像測(cè)速采集系統(tǒng)會(huì)在螺旋槳被選定的參考葉片通過固定角位置時(shí)被觸發(fā),該鎖相測(cè)量是通過TTL 信號(hào)將粒子圖像測(cè)速相機(jī)系統(tǒng)和激光脈沖同步到螺旋槳每轉(zhuǎn)一圈工作一次的頻率來實(shí)現(xiàn)的,該信號(hào)由安裝在螺旋槳測(cè)力計(jì)上的每秒3600 脈沖的旋轉(zhuǎn)增量編碼器提供,這種同步是通過可編程計(jì)時(shí)單元(LaVision PTU)實(shí)現(xiàn)的.粒子圖像測(cè)速測(cè)量系統(tǒng)如圖2 所示.

圖2 粒子圖像測(cè)速測(cè)量的試驗(yàn)設(shè)置草圖Fig.2 Sketch of the experimental setup for PIV measurements

3 計(jì)算流體力學(xué)模擬

本文數(shù)值計(jì)算的工作主要基于REX 求解器開展,REX[28]是一款通用的應(yīng)用于船舶水動(dòng)力學(xué)領(lǐng)域的自編程求解器,它在愛荷華大學(xué)得到持續(xù)開發(fā)和維護(hù).該求解器中模型的運(yùn)動(dòng)通過動(dòng)態(tài)重疊網(wǎng)格技術(shù)實(shí)現(xiàn)[29-30].REX 基于多塊結(jié)構(gòu)化貼體網(wǎng)格進(jìn)行空間離散,并支持局部網(wǎng)格的加密算法.動(dòng)量方程中速度與壓力采用PISO 算法進(jìn)行解耦,控制方程中的對(duì)流項(xiàng)采用四階迎風(fēng)格式進(jìn)行離散,而擴(kuò)散項(xiàng)的離散則采用二階中心格式處理,時(shí)間項(xiàng)應(yīng)用的是二階隱式歐拉向后差分格式進(jìn)行離散化,延遲分離渦模擬模型中遠(yuǎn)場(chǎng)大渦模擬區(qū)域的正交網(wǎng)格區(qū)域采用的是二階中心對(duì)流格式.本文中采用的延遲分離渦模擬模型和大渦模擬模型的詳細(xì)信息請(qǐng)參考文獻(xiàn)[19]和文獻(xiàn)[31].

計(jì)算域采用多塊重疊/結(jié)構(gòu)網(wǎng)格進(jìn)行離散.螺旋槳敞水條件下的計(jì)算并沒有包含粒子圖像測(cè)速試驗(yàn)中空泡水筒的壁面,其計(jì)算域包括一個(gè)延伸到下游3.5 倍螺旋槳直徑的環(huán)形加密區(qū).這套網(wǎng)格以及笛卡爾尾流加密塊均是正交的,因此可以采用低耗散的二階中心對(duì)流算法.多塊結(jié)構(gòu)化網(wǎng)格劃分以及計(jì)算核數(shù)分配細(xì)節(jié)如表2 所示,表中“Ref1”為一個(gè)長(zhǎng)方體的背景網(wǎng)格加密區(qū),而“Ref2”為一個(gè)圓環(huán)狀的柱體,用來加密下游的梢渦演化區(qū).圖3 給出了E1658螺旋槳梢渦計(jì)算網(wǎng)格分布,網(wǎng)格總數(shù)為1.927 × 108.本文計(jì)算網(wǎng)格均滿足壁面到第一層網(wǎng)格的無量綱距離y+<1 這一條件.數(shù)值模擬的無量綱時(shí)間步長(zhǎng)Δt=4.478 9×10?5,它相當(dāng)于螺旋槳每旋轉(zhuǎn)0.4675°所需要的時(shí)間.數(shù)值計(jì)算一共采用458 核進(jìn)行并行計(jì)算.

表2 螺旋槳網(wǎng)格系統(tǒng)Table 2 Details of the grid system for the propeller

圖3 敞水條件下E1658 螺旋槳的計(jì)算網(wǎng)格Fig.3 Grids for propeller E1658 in open water condition

4 結(jié)果與分析

4.1 螺旋槳尾流特性分析

本文中相平均流場(chǎng)采用下式進(jìn)行計(jì)算

式中 〈f〉代表任一物理量f的相平均結(jié)果,N為所有參考槳葉在同一角位置時(shí)的瞬態(tài)流場(chǎng)總數(shù),t0為流動(dòng)充分發(fā)展之后開始采樣的時(shí)刻,n為螺旋槳轉(zhuǎn)速.

圖4 給出了基于計(jì)算流體力學(xué)計(jì)算和粒子圖像測(cè)速測(cè)量得到的不同工況條件下的相平均渦量場(chǎng),相平均結(jié)果是基于螺旋槳每旋轉(zhuǎn)一周的鎖相結(jié)果進(jìn)行統(tǒng)計(jì)的,目的是保證梢渦配對(duì)有一個(gè)螺旋槳旋轉(zhuǎn)周期的觀察時(shí)間.其中,計(jì)算流體力學(xué)計(jì)算的相平均結(jié)果是在螺旋槳旋轉(zhuǎn)20 圈、螺旋槳尾流達(dá)到穩(wěn)定狀態(tài)之后開始統(tǒng)計(jì)的,有效避免了前期非穩(wěn)定的結(jié)果對(duì)相平均統(tǒng)計(jì)的影響.另外,圖4 也包括了粒子圖像測(cè)速試驗(yàn)測(cè)量得到的瞬態(tài)渦量場(chǎng).從圖中可以看出,計(jì)算流體力學(xué)和工程流體力學(xué)結(jié)果中都能清楚的觀察到梢渦、槳葉隨邊泄出渦以及槳轂轂渦.隨邊泄出渦存在于葉梢和根部的范圍,并很快在下游尾流中耗散了.總體而言,在同樣的網(wǎng)格和物理?xiàng)l件下,相比較于延遲分離渦模擬,大渦模擬和無湍流模型模擬能求解出更多的湍流細(xì)節(jié),并且呈現(xiàn)出更早的梢渦不穩(wěn)定性和隨邊渦的崩潰.相比較于粒子圖像測(cè)速試驗(yàn)結(jié)果,大渦模擬和無湍流模型模擬似乎過早地失去了梢渦結(jié)構(gòu)的連貫性,且在超過x/D=1 范圍的渦結(jié)構(gòu)并沒有出現(xiàn).在x/D>1 的范圍內(nèi),大渦模擬和無湍流模型模擬計(jì)算中的梢渦分解成了一系列小的渦結(jié)構(gòu).對(duì)于J= 0.56 工況的試驗(yàn)結(jié)果,見圖4(a) 梢渦配對(duì)過程大約開始于x/D= 1 的位置,在x/D= 1.5 位置處結(jié)束.然而在相同的進(jìn)速系數(shù)條件下,延遲分離渦模擬預(yù)測(cè)的梢渦配對(duì)發(fā)生在下游更遠(yuǎn)的位置,較為滯后.產(chǎn)生這種現(xiàn)象的原因可能有兩個(gè),一個(gè)是試驗(yàn)中螺旋槳進(jìn)流中存在一定的湍流度,從而觸發(fā)了更早的尾流不穩(wěn)定性;另一個(gè)是試驗(yàn)中空泡水筒存在一定的阻塞效應(yīng).值得注意的是,計(jì)算流體力學(xué)中較差的統(tǒng)計(jì)(這是由于計(jì)算量和計(jì)算資源的限制,導(dǎo)致計(jì)算流體力學(xué)模擬中可以用于相平均統(tǒng)計(jì)的樣本遠(yuǎn)遠(yuǎn)沒有粒子圖像測(cè)速試驗(yàn)的多)導(dǎo)致梢渦配對(duì)發(fā)生時(shí)有噪音簇,但是依然能觀察到由于梢渦配對(duì)產(chǎn)生的較大渦處的渦量塊.粒子圖像測(cè)速試驗(yàn)的瞬態(tài)流場(chǎng)結(jié)果表明在x/D>1.4 的范圍內(nèi)存在豐富的小尺度渦結(jié)構(gòu),但是也存在一定有組織的大尺度渦結(jié)構(gòu),這與延遲分離渦模擬的計(jì)算結(jié)果相似.

圖4 不同工況下的相平均渦量場(chǎng)Fig.4 Phase-averaged out-of-plane vorticity at different loading conditions

槳葉隨邊泄出的渦結(jié)構(gòu)呈現(xiàn)出小的局部渦結(jié)構(gòu),它們的連貫性很快就隨著繼續(xù)演化后快速消失了,這表明梢渦和槳葉隨邊泄出的這些小的渦結(jié)構(gòu)的相互干擾效應(yīng)并沒有文獻(xiàn)[17]描述的那么強(qiáng),這可能是螺旋槳幾何形狀以及工況不同的原因造成的.但是這種干擾對(duì)梢渦不穩(wěn)定性的影響是不容忽視的.計(jì)算流體力學(xué)和工程流體力學(xué)結(jié)果在J= 0.74時(shí)表現(xiàn)出了與高載荷情況時(shí)相同的趨勢(shì),見圖4(b).在這種高進(jìn)速系數(shù)工況下,梢渦不穩(wěn)定性的觸發(fā)點(diǎn)更加延后,大概在螺旋槳盤面下游x/D= 1.4 的位置,并且可以在x/D>1.75 的位置觀察到梢渦對(duì).

另外可以看到,在J= 0.56 的工程流體力學(xué)結(jié)果中,轂渦在下游很長(zhǎng)的距離都能保持相對(duì)穩(wěn)定,然而延遲分離渦模擬的結(jié)果中,轂渦較早就出現(xiàn)了崩潰的現(xiàn)象.試驗(yàn)得到的瞬態(tài)渦量場(chǎng)中,轂渦在x/D<2.5 的范圍內(nèi)都能保持連貫性.值得注意的是,大渦模擬和無湍流模型模擬方法能夠捕捉到比延遲分離渦模擬方法連貫性更強(qiáng)的轂渦.在J= 0.74 時(shí),工程流體力學(xué)結(jié)果的轂渦耗散較為嚴(yán)重,這與延遲分離渦模擬模擬較為一致.

利用粒子圖像測(cè)速試驗(yàn)測(cè)量得到的所有槳葉相(36 個(gè))的相平均結(jié)果,進(jìn)行螺旋槳尾流3D 漩渦結(jié)構(gòu)的重構(gòu).如圖5 所示,在圓柱坐標(biāo)系統(tǒng)中導(dǎo)入流場(chǎng)試驗(yàn)得到的所有葉相的相平均流場(chǎng),并按照下式將網(wǎng)格插值到圓柱網(wǎng)格上

圖5 試驗(yàn)結(jié)果三維梢渦結(jié)構(gòu)的插值方法Fig.5 Interpolation methods for calculating 3D tip vortex structures of experiments

式中x,y,z為插值得到的新坐標(biāo),r=|y0(i,j,k)|,x0,y0為原始相平均流場(chǎng)中的坐標(biāo),i和j分別為兩個(gè)方向的索引,k為不同的槳葉相對(duì)應(yīng)的流場(chǎng)索引, φ 為不同葉相對(duì)應(yīng)的角度.流場(chǎng)按照下式插值到圓柱網(wǎng)格上

式中u,v,w和u0,v0分別為插值后的3 個(gè)速度分量和原始相平均流場(chǎng)(2D-2C)中的兩個(gè)速度分量.

圖6 給出了不同載荷條件下,采用上述方法對(duì)工程流體力學(xué)結(jié)果重構(gòu)的三維相平均漩渦結(jié)構(gòu),其中渦結(jié)構(gòu)采用無量綱軸向速度進(jìn)行著色.三維渦結(jié)構(gòu)采用渦量等值面進(jìn)行顯示,即 ωy=100的等值面.為了顯示清晰的梢渦結(jié)構(gòu),圖6 中內(nèi)半徑處的渦結(jié)構(gòu)即轂渦部分被移除了.圖6 也給出了基于延遲分離渦模擬方法模擬得到的相平均三維渦結(jié)構(gòu),漩渦分別定義為J= 0.56 的Q= 15000 等值面以及J=0.74 的Q= 20000 等值面.試驗(yàn)結(jié)果顯示梢渦的最初合并是靠近螺旋槳盤面成對(duì)的,然后是各對(duì)之間的合并.由于漩渦的蜿蜒會(huì)引起相位平均渦度分辨率的損失,使得很難識(shí)別到成對(duì)的、三重和四重漩渦.相比于兩種載況條件下的計(jì)算流體力學(xué)結(jié)果,試驗(yàn)結(jié)果的初始漩渦合并更靠近螺旋槳盤面,梢渦的不穩(wěn)定性更早地被觸發(fā).

圖6 不同工況下相平均三維渦結(jié)構(gòu)Fig.6 Three-dimensional phase-averaged vortical structures at different loading conditions

圖7 為基于延遲分離渦模擬、大渦模擬和無湍流模型模擬方法計(jì)算得到的瞬態(tài)螺旋槳尾渦結(jié)構(gòu),對(duì)于J= 0.56 工況,漩渦結(jié)構(gòu)由Q= 200000 等值面顯示,對(duì)于J= 0.74,漩渦結(jié)構(gòu)由Q= 100000 等值面顯示.這些瞬態(tài)視圖更好地展現(xiàn)了相位平均結(jié)果的梢渦配對(duì)和不穩(wěn)定性.從不同工況以及不同數(shù)值方法得到的結(jié)果可以看到構(gòu)成漩渦系統(tǒng)的梢渦、隨邊渦以及轂渦.與延遲分離渦模擬方法計(jì)算得到的結(jié)果相比,大渦模擬和無湍流模型模擬方法預(yù)測(cè)了更早的漩渦不穩(wěn)定性和崩潰,因此難以觀測(cè)到試驗(yàn)中的梢渦配對(duì)現(xiàn)象.在近場(chǎng)中可以觀察到橢圓不穩(wěn)定性,隨著梢渦彼此之間的距離越來越近,橢圓不穩(wěn)定性會(huì)不斷增強(qiáng).漩渦合并會(huì)產(chǎn)生顯著的順流方向渦量,并且在梢渦崩潰之前可以觀察到更大范圍的橢圓不穩(wěn)定性現(xiàn)象.在螺旋槳下游更遠(yuǎn)處,梢渦趨于分解為逐漸變小的渦結(jié)構(gòu),從而完全抑制了螺旋渦的合并.在相同時(shí)間和網(wǎng)格分辨率下,相比于大渦模擬模型,無湍流模型模擬模型似乎能求解更多的湍流,因此其預(yù)報(bào)的梢渦不穩(wěn)定性發(fā)生的更早.

延遲分離渦模擬方法得到的結(jié)果顯示了隨著合并過程的發(fā)展梢渦是如何相互纏繞的,但是梢渦的完全合并在瞬態(tài)計(jì)算結(jié)果中并未發(fā)生.在J= 0.56時(shí),梢渦配對(duì)會(huì)更早的發(fā)生,并且合并的漩渦在精細(xì)的計(jì)算域內(nèi)會(huì)變得不穩(wěn)定,從而展現(xiàn)了漩渦系統(tǒng)崩潰的過程,這在圖6 的相平均漩渦結(jié)構(gòu)中表現(xiàn)為多個(gè)梢渦合并時(shí)的嘈雜尾流.在J= 0.74 時(shí),梢渦的合并進(jìn)一步向下游延后,并且在計(jì)算流體力學(xué)結(jié)果無法精細(xì)求解漩渦結(jié)構(gòu)的梢渦加密區(qū)末端變得不穩(wěn)定.圖7 中延遲分離渦模擬結(jié)果在遠(yuǎn)場(chǎng)的漩渦結(jié)構(gòu)與文獻(xiàn)[26]中E779A 螺旋槳的試驗(yàn)結(jié)果相似,試驗(yàn)中E779A 螺旋槳的梢渦分組但并未合并,另外本文計(jì)算結(jié)果中的梢渦蛙跳現(xiàn)象和其他現(xiàn)象也與該試驗(yàn)中捕捉到的梢渦行為相近.

螺旋槳尾流結(jié)構(gòu)被漩渦合并所支配.圖8 給出了計(jì)算流體力學(xué)計(jì)算得到的梢渦合并過程,Q準(zhǔn)則等值面顯示的渦結(jié)構(gòu)突出了梢渦的合并過程.為了便于分析,圖8 對(duì)從不同槳葉泄出的梢渦進(jìn)行了編號(hào).從圖8 可以看到,4 號(hào)梢渦將3 號(hào)梢渦推進(jìn)內(nèi)部軸向速度較高的尾流區(qū),本質(zhì)上是包裹住梢渦以使其靠近并最終完成合并.由于初始合并成對(duì)發(fā)生且當(dāng)前螺旋槳有7 個(gè)葉片,因此合并對(duì)的角間隔應(yīng)該在大約103°處發(fā)生,這與1 號(hào)梢渦和2 號(hào)梢渦對(duì)相對(duì)于3 號(hào)梢渦和4 號(hào)梢渦對(duì)之間的距離接近.同樣,如果螺旋槳葉片的數(shù)量是偶數(shù),則所有成對(duì)合并的梢渦都應(yīng)該在螺旋槳盤面下游大約相同距離處產(chǎn)生互感.渦對(duì)1 + 2 和3 + 4 確實(shí)在圖8 中的相同軸向位置合并.在具有奇數(shù)個(gè)葉片的螺旋槳中,必定有一個(gè)梢渦與已經(jīng)合并的一對(duì)梢渦進(jìn)行合并,這是因?yàn)楹喜⒑蟮南噜彍u對(duì)會(huì)產(chǎn)生更強(qiáng)的感應(yīng).在圖8 中梢渦5 會(huì)與6 + 7 渦對(duì)合并,形成新的組合渦5 + 6 + 7.梢渦6 和梢渦7 在梢渦1 + 2 和梢渦3 + 4 合并的下游更遠(yuǎn)處合并,這很可能是由于奇數(shù)葉片配置的螺旋槳的不對(duì)稱效應(yīng)引起的.最后,漩渦1 + 2,3 +4 和5 + 6 + 7 會(huì)產(chǎn)生進(jìn)一步的互感效應(yīng),導(dǎo)致漩渦1 + 2 和漩渦3 + 4 合并為漩渦1 + 2 + 3 + 4.隨著漩渦強(qiáng)度由于擴(kuò)散而衰減(或漩渦由于不穩(wěn)定性而崩潰),并且由于合并而導(dǎo)致漩渦之間的距離增加,因此在尾流中出現(xiàn)單個(gè)(理想狀態(tài)下最終會(huì)合并成單個(gè)渦)渦之前,該過程可能會(huì)停止.另外,圖8 也顯示了漩渦合并之前和合并之后更加強(qiáng)烈的橢圓不穩(wěn)定性的存在.

圖9 給出了兩個(gè)工況條件下槳后的瞬態(tài)渦量云圖,圖中可以清晰的看到梢渦不穩(wěn)定性的觸發(fā)過程.梢渦和隨邊渦分別從槳葉的葉梢和隨邊泄出,在向下游演化的過程中,梢渦會(huì)與其上游相鄰梢渦的隨邊渦產(chǎn)生干擾(如圖9 中的E 區(qū)和F 區(qū)所示),并伴有各個(gè)梢渦之間的互感效應(yīng).在繼續(xù)向下游演化的過程中,梢渦的不穩(wěn)定性被觸發(fā),如圖9 所示的A 區(qū)、B 區(qū)、C 區(qū)和D 區(qū)所示,具體表現(xiàn)為這4 個(gè)區(qū)域中的下游梢渦將上游梢渦推入內(nèi)半徑處高速剪切層流中,這與圖7 中的分析一致,產(chǎn)生這一現(xiàn)象的原因是梢渦在向下游演化的過程中,其周向速度會(huì)逐漸出現(xiàn)擴(kuò)散,最終它們會(huì)有足夠的周向速度來抵達(dá)相鄰梢渦,這會(huì)導(dǎo)致梢渦互感效應(yīng),其中渦齡越大的梢渦(從槳葉泄出時(shí)間長(zhǎng),即下游梢渦)周向速度擴(kuò)散更為嚴(yán)重,直到梢渦演化到流場(chǎng)中的某一點(diǎn)處時(shí),下游梢渦有足夠的周向速度作用在相鄰的上游梢渦上,從而在高速剪切層流和隨邊渦的雙重影響下,觸發(fā)梢渦的不穩(wěn)定性.

圖9 不同工況下梢渦不穩(wěn)定性觸發(fā)過程Fig.9 Instability inception process of tip vortices under different conditions

4.2 梢渦演化模型

二維同向旋轉(zhuǎn)漩渦的不穩(wěn)定性和配對(duì)問題是湍流動(dòng)力學(xué)中的重要問題,受到了大量學(xué)者的關(guān)注[32-33].螺旋槳尾流中就包含著這一基本問題,因?yàn)槁菪龢覝u也是一系列天然的同向旋轉(zhuǎn)漩渦,但它們又在很多方面都比線性渦對(duì)復(fù)雜得多.它們除了呈現(xiàn)出螺旋狀的幾何結(jié)構(gòu),也同時(shí)沉浸在由螺旋槳推力和尾流中的高速區(qū)造成的剪切層中.另外,螺旋槳產(chǎn)生了一個(gè)尾流渦旋,它的旋轉(zhuǎn)速度取決于螺旋槳載荷,大概是螺旋槳轉(zhuǎn)速的1%～5%.更重要的是,槳葉升力產(chǎn)生的隨邊渦可以與螺旋槳梢渦相互作用,這都使得螺旋槳尾渦系統(tǒng)異常的復(fù)雜.螺旋槳尾渦結(jié)構(gòu)如圖10 所示,其中p為梢渦螺距.

圖10 螺旋槳尾流結(jié)構(gòu)Fig.10 Structure of propeller wake

本文4.1 節(jié)對(duì)螺旋槳尾流動(dòng)力學(xué)特性包括不穩(wěn)定性觸發(fā)和梢渦配對(duì)進(jìn)行了詳細(xì)的計(jì)算流體力學(xué)和工程流體力學(xué)分析,結(jié)論表明螺旋槳梢渦的互感效應(yīng)以及梢渦?葉尾流(隨邊渦)的相互干擾是螺旋槳尾流不穩(wěn)定性觸發(fā)的主要原因(計(jì)算流體力學(xué)計(jì)算中不考慮來流中的背景湍流度),其中梢渦周向速度隨時(shí)間的耗散在梢渦互感效應(yīng)中扮演者十分重要的作用,相關(guān)定性的機(jī)理分析能夠較好地對(duì)螺旋槳梢渦的不穩(wěn)定性行為進(jìn)行解釋,但有關(guān)螺旋槳梢渦演化規(guī)律的定量描述由于漩渦合并機(jī)制以及螺旋槳尾渦系統(tǒng)的復(fù)雜性至今仍沒有公開發(fā)表的文獻(xiàn)涉及,大多數(shù)文獻(xiàn)目前都只對(duì)這一問題進(jìn)行了多工況下的多角度定性分析.基于探索得到的均勻來流中螺旋槳梢渦的演化機(jī)理,本文提出了螺旋槳梢渦的演化模型,該模型基于渦?渦之間周向誘導(dǎo)速度相互作用的方式,首次定量地描述了梢渦之間的互感效應(yīng)和隨邊渦與梢渦之間的相互作用以及梢渦不穩(wěn)定性的觸發(fā).該模型包含了本文目前已探明的螺旋槳尾流不穩(wěn)定性機(jī)理,能夠較為準(zhǔn)確地模擬螺旋槳梢渦的演化規(guī)律,預(yù)測(cè)螺旋槳梢渦融合的時(shí)間和位置,相關(guān)研究為進(jìn)一步理解螺旋槳尾流的動(dòng)力學(xué)特征提供了指導(dǎo)作用.

通過4.1 節(jié)的分析可以知道,對(duì)于某給定位置的梢渦,其會(huì)分別受到一側(cè)更年輕(渦齡越小,即從槳葉泄出的時(shí)間越短)的梢渦和另一側(cè)更老的梢渦周向誘導(dǎo)速度的影響,且更老梢渦的影響更加強(qiáng)烈.螺旋槳下游不同渦齡梢渦示意圖如圖11 所示.

圖11 螺旋槳下游不同渦齡梢渦Fig.11 Downstream tip vortices with different ages

考慮二維梢渦周向速度由以下極坐標(biāo)系下的無量綱方程控制

式中vθ為周向速度,其初始速度剖面從計(jì)算流體力學(xué)結(jié)果中螺旋槳槳葉后的第一個(gè)參考梢渦中提取.邊界條件的選取也十分關(guān)鍵,為了符合本模型的物理機(jī)理,左右邊界的邊界條件均為vθ=0,需要注意的是,如果右邊界不為零,得到的速度剖面數(shù)值解為L(zhǎng)amb?Oseen 渦的速度剖面形式,顯然該形式不滿足本模型的假設(shè).雷諾數(shù)計(jì)算的參考速度和參考長(zhǎng)度與計(jì)算流體力學(xué)計(jì)算中保持一致,r為任意點(diǎn)到渦核中心的距離.

螺旋槳尾渦與傳統(tǒng)研究中同向旋轉(zhuǎn)漩渦不穩(wěn)定性問題的另一個(gè)區(qū)別是,螺旋槳尾流中的梢渦存在軸向遷徙速度,即螺旋槳尾渦系統(tǒng)中的渦結(jié)構(gòu)在相互作用的同時(shí)也存在軸向遷徙.風(fēng)機(jī)需要從流場(chǎng)中抽取能量,因此其尾流是向外徑擴(kuò)散的,而螺旋槳工作的原理是加速尾流產(chǎn)生推力,因此其尾流是收縮的.螺旋槳后收縮尾流中的渦結(jié)構(gòu)演化也給其不穩(wěn)定性的研究帶來了復(fù)雜性.為了得到更加真實(shí)的渦系演化模擬結(jié)果,本文提出的模型也考慮了中縱剖面內(nèi)梢渦沿軸向和垂向的基礎(chǔ)遷徙速度,這兩個(gè)速度從時(shí)間平均流場(chǎng)中雙線性插值而來,且取決于梢渦和隨邊渦在流場(chǎng)中的具體位置.為了便于計(jì)算渦結(jié)構(gòu)在每個(gè)時(shí)間步的位移,渦結(jié)構(gòu)上的誘導(dǎo)速度被分解為沿軸向和垂向的速度成分,速度分解的原則如圖12 所示.

圖12 兩個(gè)渦之間的誘導(dǎo)速度分解Fig.12 Induced velocity decomposition between two vortices

本模型中初始時(shí)刻漩渦系統(tǒng)配置的描述如圖13所示,圖13 中尺寸較大的圓代表槳后的梢渦,而尺寸較小的一系列圓代表隨邊渦,用來模擬槳葉尾流.

圖13 中,以大圓的圓心為垂向位置的零點(diǎn),小圓則在垂向方向上等間距的排列,與尾流中隨邊渦的特征尺度一致.每過一個(gè)螺旋槳槳葉通道的時(shí)間,就會(huì)初始化一套圖13 中的梢渦?隨邊渦系統(tǒng),這與螺旋槳在運(yùn)行時(shí)生成槳后尾渦系統(tǒng)的過程一致.

圖13 初始漩渦系統(tǒng)配置草圖Fig.13 Sketch of the initial configuration of vortex system

渦結(jié)構(gòu)之間的相互作用方式是本模型的核心問題之一,為了簡(jiǎn)化渦結(jié)構(gòu)之間的相互作用,多個(gè)漩渦之間的作用采用線性疊加的方式進(jìn)行評(píng)估.假設(shè)t時(shí)刻,流場(chǎng)中存在M個(gè)梢渦,N個(gè)隨邊渦,則基于線型疊加得到的流場(chǎng)中第j個(gè)梢渦的誘導(dǎo)速度可以表示為

式中uVi是第i個(gè)梢渦誘導(dǎo)的速度,uTi是第i個(gè)隨邊渦誘導(dǎo)的速度,xVi代表第i個(gè)梢渦的位置,xTi代表第i個(gè)隨邊渦的位置.um是從螺旋槳時(shí)間平均尾流場(chǎng)中獲得的速度.在每個(gè)時(shí)間步渦的速度計(jì)算完之后,渦的位移可已通過速度對(duì)時(shí)間的積分計(jì)算,例如流場(chǎng)中第j個(gè)梢渦的位移可通過下式計(jì)算

圖14 給出了演化模型運(yùn)行程序的流程圖,具體計(jì)算流程如下:

圖14 演化模型的計(jì)算策略Fig.14 Computation strategy of the evolution model

(1)準(zhǔn)備渦系統(tǒng).讀取初始條件,求解誘導(dǎo)速度偏微分方程,并給定初始漩渦的位置,計(jì)算各漩渦之間的距離和分解角度,讀取各自的周向誘導(dǎo)速度,開始迭代.

(2)更新誘導(dǎo)速度.主要分為三部分.第一,梢渦在梢渦上產(chǎn)生的誘導(dǎo)速度;第二,梢渦在隨邊渦上產(chǎn)生的誘導(dǎo)速度;第三,隨邊渦在梢渦上產(chǎn)生的誘導(dǎo)速度.并插值獲得基本遷徙速度,計(jì)算得到軸向和垂向速度分量.

(3)更新漩渦位置.根據(jù)上一步中計(jì)算得到的速度分量獲得軸向和垂向位移,根據(jù)此位移計(jì)算得到漩渦新的位置.

(4)更新渦距和角度.根據(jù)漩渦在流場(chǎng)中新的位置,計(jì)算各個(gè)漩渦之間的渦距并計(jì)算它們的速度分解角度.

(5)判定計(jì)算是否滿足停止條件,即有兩個(gè)梢渦在軸向上位置重合,如是則計(jì)算停止,如否則判定時(shí)間是否等于一個(gè)槳葉通道的時(shí)間,如否繼續(xù)從第2 步開始迭代,如是則向流場(chǎng)中添加新的梢渦和隨邊渦初始系統(tǒng).

(6)更新新進(jìn)梢渦和隨邊渦與流場(chǎng)中其他漩渦的渦距和速度分解角.

用來驗(yàn)證此模型的工況為計(jì)算流體力學(xué)計(jì)算中的J= 0.56 和J= 0.74 工況,梢渦的周向速度剖面隨時(shí)間的變化通過式(4)求解得到,隨邊渦的初始周向速度在計(jì)算流體力學(xué)計(jì)算中從對(duì)應(yīng)的垂向位置提取.漩渦的基本遷徙速度采用雙線性插值的方法在梢渦和隨邊渦對(duì)應(yīng)位置的時(shí)間平均流場(chǎng)中提取.圖15給出了采用梢渦演化模型預(yù)報(bào)得到的J= 0.56 和J=0.74 工況漩渦融合時(shí)對(duì)應(yīng)的瞬態(tài)漩渦系統(tǒng).從圖15中可以看出梢渦融合的位置和時(shí)間均與圖9 中的計(jì)算流體力學(xué)計(jì)算結(jié)果吻合較好,表明本文挖掘的梢渦不穩(wěn)定性觸發(fā)物理機(jī)制較為準(zhǔn)確,提出的梢渦演化模型和算法能夠有效地模擬梢渦的演化和梢渦的不穩(wěn)定性觸發(fā)過程.

圖15 采用模型預(yù)報(bào)得到的J = 0.56 和J = 0.74 時(shí)的瞬態(tài)漩渦系統(tǒng)Fig.15 Instantaneous vortex system predicted by proposed model at J = 0.56 and J = 0.74

5 結(jié)論

長(zhǎng)期以來,螺旋槳尾流的不穩(wěn)定性研究無論是從工程技術(shù)的角度,還是作為一個(gè)基礎(chǔ)的流體力學(xué)問題,都一直深受廣大學(xué)者的關(guān)注.但就目前公開發(fā)表的文獻(xiàn)來看,不管是粒子圖像測(cè)速試驗(yàn)還是高精度的計(jì)算流體力學(xué)模擬,大多數(shù)有關(guān)螺旋槳尾流不穩(wěn)定性的研究都是定性的,往往都是通過觀察到的流場(chǎng)信息和渦結(jié)構(gòu)行為來定性的描述其不穩(wěn)定性的機(jī)理,這主要是由槳后渦系和演化過程本身的復(fù)雜性及其所處的剪切層流環(huán)境造成的.計(jì)算流體力學(xué)模擬基于延遲分離渦模擬、大渦模擬和無湍流模型模擬模型開展,結(jié)果表明與延遲分離渦模擬方法相比,大渦模擬和無湍流模型模擬方法預(yù)測(cè)了更早的漩渦不穩(wěn)定性和崩潰,因此難以觀測(cè)到試驗(yàn)中的梢渦配對(duì)現(xiàn)象.延遲分離渦模擬模型能較好地捕捉到試驗(yàn)中的梢渦配對(duì)行為,但預(yù)報(bào)的梢渦合并過程與粒子圖像測(cè)速試驗(yàn)相比有所延后.本文基于探索得到的均勻來流中螺旋槳梢渦的演化機(jī)理,提出了螺旋槳梢渦的演化模型,該模型基于渦?渦之間周向誘導(dǎo)速度相互作用的方式,首次定量地描述了梢渦之間的互感效應(yīng)和隨邊渦與梢渦之間的相互作用以及梢渦不穩(wěn)定性的觸發(fā),并能預(yù)測(cè)螺旋槳梢渦融合的時(shí)間和位置.相關(guān)研究為進(jìn)一步理解螺旋槳尾流的動(dòng)力學(xué)特征、螺旋槳基于流動(dòng)控制方法的降噪措施實(shí)施、以及新型螺旋槳的優(yōu)化設(shè)計(jì)有著重要意義.