自由面對潛艇尾流場流動(dòng)特性影響研究1)

郭春雨 徐鵬 韓陽,2) 王超 郭欣雨 鄶云飛

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

?(水下測控技術(shù)重點(diǎn)實(shí)驗(yàn)室,大連 116013)

引言

潛艇的隱身性能不僅是保障完成其作戰(zhàn)使命的基本性能,更有助于潛艇在海戰(zhàn)中更好地發(fā)揮其強(qiáng)大的突襲能力.隨著如今作戰(zhàn)任務(wù)多樣性的提升,潛艇不可避免地會(huì)涉及到近自由面航行,艇體的分離渦與自由面相互作用,致使渦結(jié)構(gòu)會(huì)在原有結(jié)構(gòu)基礎(chǔ)上產(chǎn)生位置偏移,進(jìn)一步加劇流場的復(fù)雜程度,這對潛艇的隱身性能來說是一項(xiàng)更大的挑戰(zhàn).由于戰(zhàn)略戰(zhàn)術(shù)的需求,有關(guān)潛艇近水面航行的研究逐漸受到重視.近年來隨著水動(dòng)力學(xué)研究的深入,越來越多的科研人員通過流場流動(dòng)細(xì)節(jié)來尋找影響宏觀力及噪聲的本質(zhì)原因.對于精細(xì)流場的獲取主要有兩種途徑:試驗(yàn)測試與數(shù)值模擬.其中數(shù)值模擬由于成本較低且能獲取全場流動(dòng)信息而受到廣泛應(yīng)用,但其可信程度稍低,往往需要試驗(yàn)數(shù)據(jù)加以驗(yàn)證.隨著流場測試技術(shù)的不斷發(fā)展,粒子圖像測速(particle image velocimetry,PIV)因其無接觸、全局、瞬態(tài)測量等優(yōu)點(diǎn),迅速成為研究者的利器[1-4].

國內(nèi)外公開發(fā)表的運(yùn)用PIV 技術(shù)測量潛艇繞流場的文章很少.Ashok 等[5-6]對SUBOFF 標(biāo)模進(jìn)行了系列試驗(yàn)并進(jìn)行了分析.Fureby 等[7]對直航和10?斜航角狀態(tài)下的潛艇流場進(jìn)行測量,并與數(shù)值模擬結(jié)果進(jìn)行對比,兩者吻合良好.Lee 等[8]通過SPIV 技術(shù)測量了潛艇10?斜航角狀態(tài)下指揮臺圍殼分離渦情況.Wang 等[9]研究了自由面的存在對螺旋槳尾流的影響,并與CFD 計(jì)算結(jié)果相互驗(yàn)證.在國內(nèi),張軍等[10-11]利用隨車式PIV 系統(tǒng)測量了不同主體與附體連接形式的水下航行體尾部流場.張國平等[12]應(yīng)用隨車式2D-3C 水下PIV 測量系統(tǒng)進(jìn)行了SUBOFF潛艇流場測試,然而并未公開相關(guān)測試結(jié)果.薛慶雨等[13-15]在空泡水筒中測量了SUBOFF 模型直航、10?攻角、尾附體10?舵角等三種工況下槳盤面流場湍流強(qiáng)度和湍流積分長度.海軍工程大學(xué)張文照等[16]對潛艇模型的尾流速度場進(jìn)行測量,得到潛艇模型的尾流渦量場分布,為研究潛艇水動(dòng)力噪聲與尾流渦結(jié)構(gòu)之間的關(guān)系提供了參考.

由于流場測量困難,鮮有考慮自由面對潛艇尾流場流動(dòng)特性影響的試驗(yàn)研究,大都采用數(shù)值模擬的方法探究艇體與自由面相互作用的影響.Dawson[17]是詳細(xì)開展?jié)撏Ы杂擅婧叫袛?shù)值模擬較早的研究者之一.Shariati 等[18-19]應(yīng)用RANS 方法對SUBOFF 模型多種工況、航速、潛深做了詳細(xì)研究.近期,Amiri等[20-22]對自由面如何影響潛艇水動(dòng)力性能進(jìn)行了系列研究.2019 年Carrica 等[23]考慮靜水環(huán)境及波浪環(huán)境對Joubert BB2 潛艇近自由面自航的水動(dòng)力性能進(jìn)行了研究.Tian 等[24]研究了不同波浪作用以及不同潛深對水下航行器水動(dòng)力的影響.在國內(nèi),張楠等[25-28]開展了大量關(guān)于潛艇近水面航行數(shù)值模擬研究,奠定了國內(nèi)自由面對潛艇尾流影響數(shù)值研究的基礎(chǔ).隨后,逐漸受到各科研高校的關(guān)注,涌現(xiàn)出大批有關(guān)自由面對潛艇尾流影響的論文[29-34].

采用精細(xì)流場測量技術(shù)對潛艇繞流場的試驗(yàn)研究很少,有關(guān)自由面對潛艇尾流場流動(dòng)特性影響的研究更是少之又少.本文依托哈爾濱工程大學(xué)船模水池實(shí)驗(yàn)室,采用先進(jìn)的水下三維粒子圖像測速(stereoscopic particle image velocimetry,SPIV) 系統(tǒng)測量潛艇尾流場微觀結(jié)構(gòu),輔以數(shù)值模擬技術(shù)觀測近自由面航行時(shí)表面波系,以此來研究SUBOFF 潛艇無限潛深與近水面航行時(shí)尾流場特性,通過試驗(yàn)與數(shù)值模擬結(jié)果相互驗(yàn)證、補(bǔ)充,得到潛艇尾流場流動(dòng)特征以及隨航速與潛深變化的規(guī)律,以期為未來潛艇設(shè)計(jì)提供依據(jù).

1 試驗(yàn)方案

1.1 試驗(yàn)設(shè)備與試驗(yàn)?zāi)Ｐ瓦x取

1.1.1 試驗(yàn)設(shè)備

本文依托哈爾濱工程大學(xué)船模拖曳水池開展SUBOFF 模型尾流場SPIV 試驗(yàn),基本實(shí)驗(yàn)設(shè)備如下.

(1)拖曳水池:長度108 m;寬度7 m;水深3.5 m.

(2) 拖車:車速V≤6.5 m/s;穩(wěn)速范圍:0.1～6.5 m/s;精度:0.1%;通過自主設(shè)計(jì)并制作的固定裝置將試驗(yàn)?zāi)Ｐ凸潭ㄔ谕宪嚿?固定方式采用艇尾固定.

(3)水下粒子圖像測速系統(tǒng):2D-3C SPIV 系統(tǒng).

1.1.2 模型選取

本試驗(yàn)選用的模型是美國泰勒水池(David Taylor Model Basin,DTMB) 設(shè)計(jì)的SUBOFF 模型,該模型具有大量的試驗(yàn)結(jié)果進(jìn)行驗(yàn)證,因此比較適合做相關(guān)研究.由于水池深度及PIV 測試能力的限制,模型長度選為1 m,模型如圖1 所示,其主要參數(shù)如表1所示.

圖1 SUBOFF 模型圖Fig.1 SUBOFF model

試驗(yàn)所用模型由聚乳酸(PLA) 材料經(jīng)3D 打印技術(shù)制作而成.聚乳酸密度為1.20～1.30 kg/L,彈性模量為3000～4000 MPa,拉伸強(qiáng)度為40～60 MPa.由于打印機(jī)工作區(qū)域限制,整個(gè)艇體分三段打印,最后拼接而成.首先將艇體表面打磨光順并噴涂防水涂層,然后對艇體表面進(jìn)行亞光黑處理,最后對測量截面附近進(jìn)行噴涂吸光材料,避免試驗(yàn)過程中激光強(qiáng)反光影響.由于本次試驗(yàn)采用的潛艇模型縮尺比較小,故在艇體首部L/20 處加裝激流裝置以加快艇體表面的層流轉(zhuǎn)淚.經(jīng)3D 打印的實(shí)物模型以及噴涂效果如圖2 所示.

圖2 SUBOFF 實(shí)物模型Fig.2 SUBOFF solid model

1.2 試驗(yàn)工況的確定

為分析自由面對潛艇尾流場的影響,進(jìn)行三個(gè)潛深、多個(gè)航速下的潛艇尾流場測量.由文獻(xiàn)[17]可知,潛深在3.3 倍最大直徑D以上時(shí)可以忽略自由面的影響;同時(shí)由另一文獻(xiàn)[35]可知,除了自由面影響外,池底壁面也存在影響,要完全避免限制水域的干擾需要水下航行體距池底的距離不宜小于1/2 倍艇長.因此本次試驗(yàn)選取4D潛深作為無限潛深工況,此時(shí)可以忽略自由面和池底的影響.同時(shí),在潛艇靠近自由面時(shí)存在自由液面反光干擾,為減小反光現(xiàn)象,選取1.5D潛深作為近自由面工況,選取2.5D潛深作為中間工況.試驗(yàn)航速選取為Fr=0.3～0.7,傅汝德數(shù)Fr定義如下

式中,U為拖曳航速,L為潛艇的特征長度.

測量截面取為槳盤面處,模型參數(shù)及試驗(yàn)工況如表1 所示.

表1 模型參數(shù)及試驗(yàn)工況Table 1 Model parameters and test conditions

1.3 數(shù)值模型與網(wǎng)格劃分

采用基于有限體積法的STAR-CCM+商業(yè)軟件進(jìn)行潛艇流場數(shù)值模擬,應(yīng)用隱式非定常時(shí)間推進(jìn)格式進(jìn)行非穩(wěn)態(tài)求解,湍流模型選取分離渦模型(DES),時(shí)間步長的選取則充分考慮CFL～1 的條件,綜合網(wǎng)格空間尺度及時(shí)間尺度,最終選為0.001 s.由于涉及到近自由面航行,自由表面采用兩相體積VOF 方法建模.控制方程采用基于壓力的耦合求解,其中對流項(xiàng)采用二階迎風(fēng)格式進(jìn)行空間離散,耗散項(xiàng)采用二階中心差分格式進(jìn)行離散.計(jì)算時(shí)長為定為2 s.

計(jì)算域如圖3 所示,具體設(shè)定如下:速度入口邊界距艇首2 倍艇長L,壓力出口距艇尾4 倍艇長L,其余邊界設(shè)為壁面邊界條件,距艇中心10 倍艇徑D.

圖3 計(jì)算域示意圖Fig.3 The calculation domain schematic

網(wǎng)格劃分采用軟件中自帶的切割體網(wǎng)格形式,在艇體周圍進(jìn)行三層局部加密,在主附體交接部分及分離渦發(fā)展區(qū)域(指揮臺圍殼后端及十字舵后端)進(jìn)行局部加密,對潛艇近水面航行時(shí)的自由表面和開爾文波形等物理現(xiàn)象也進(jìn)行加密.邊界層網(wǎng)格劃分則采用棱柱層網(wǎng)格形式,并采用Ally+壁面處理,同時(shí)根據(jù)不同航速進(jìn)行棱柱層厚度調(diào)整,保證y+值處于1 左右范圍內(nèi).根據(jù)不同航速下的興波程度有針對性地加密自由表面及開爾文波形,最終將網(wǎng)格數(shù)控制在1000 萬左右,艇體周圍各加密區(qū)細(xì)節(jié)如圖4 所示.

2 潛艇尾流場測量試驗(yàn)

2.1 水下SPIV 測量系統(tǒng)

圖4 艇體周圍加密區(qū)示意圖Fig.4 The mesh refinement area

本試驗(yàn)所采用的流場測量設(shè)備為DANTEC 公司為哈爾濱工程大學(xué)船模拖曳水池定制的隨車式水下SPIV 系統(tǒng),布置于拖車右側(cè),SPIV 系統(tǒng)規(guī)格參數(shù)如下:CCD 相機(jī)分辨率為2048×2048 像素;雙脈沖Nd-YAG 激光器,最大脈沖頻率為15 Hz;最大脈沖激光能量為1200 mJ;常規(guī)工作能量為200 mJ;激光束持續(xù)時(shí)長為4 ns;激光波長為532 nm;脈沖激光片厚度為0.6 mm;最大測量范圍為400 mm×400 mm;示蹤粒子為聚酰胺顆粒.

該套系統(tǒng)包括一個(gè)浸沒水中的雷體、兩個(gè)翼型支架、脈沖式激光器、同步器、圖像采集裝置、粒子播撒裝置等.雷體包含相機(jī)及眾多反射鏡和柱面鏡結(jié)構(gòu),其內(nèi)部詳細(xì)光學(xué)結(jié)構(gòu)布置如圖5 所示.雷體內(nèi)的兩臺相機(jī)為Dantec 公司的FlowSense 4M MKII 12位精度相機(jī)組,其使用高性能的逐行掃描行間CCD芯片,相機(jī)的分辨率為2048×2048 像素并且具有相同數(shù)目的存儲(chǔ)單元.兩臺相機(jī)的鏡頭參數(shù)各不相同,靠近前端的Camera1 鏡頭為50 mm f/1.8 的佳能鏡頭,靠近后端的Camera2 鏡頭為85 mm f/1.8 的佳能鏡頭.支柱1 中包含激光導(dǎo)光臂,支柱2 中包含相機(jī)線纜支柱結(jié)構(gòu).

圖5 SPIV 系統(tǒng)中各元器件詳細(xì)布置圖[36]Fig.5 Detailed arrangement of the probes in the underwater SPIV system[36]

圖6 SPIV 系統(tǒng)實(shí)物圖Fig.6 Stereoscopic particle image velocimetry system

2.2 流場數(shù)據(jù)采集、測量與分析

采用平均直徑為50μm 聚酰胺顆粒(PSP-50)作為示蹤粒子,密度為1.03 g/cm3,具有良好的跟隨性.由于拖曳水池尺度較大,如何在拖曳水池進(jìn)行大范圍粒子均勻播撒是一項(xiàng)技術(shù)挑戰(zhàn),本試驗(yàn)應(yīng)用一套自主設(shè)計(jì)的粒子播撒裝置進(jìn)行示蹤粒子播撒,播撒裝置如圖7 所示.裝置上方裝有攪拌電機(jī)的水箱可以存儲(chǔ)大量的水并均勻攪拌粒子,打開水箱下方的閥門后,粒子沿下方不同管道播撒到水池中,下方每條管道都開有若干小孔,其潛深可以在水下50 cm 范圍內(nèi)自由調(diào)節(jié),完全符合潛艇的潛深要求,在不同深度下均能保證潛艇模型周圍存在適量粒子.每航次試驗(yàn)結(jié)束后進(jìn)行示蹤粒子補(bǔ)充,通過相機(jī)觀測窗口進(jìn)行粒子濃度判別,保證每次測量過程中平均粒子濃度約為150 g/m3即每個(gè)查問區(qū)域具有15～30 個(gè)粒子,這樣能夠增大互相關(guān)分析的精確度.

在試驗(yàn)過程中,脈沖激光照亮待測區(qū)域示蹤粒子,相機(jī)雙幀拍攝,采集頻率調(diào)至最大采樣頻率7.5 Hz,根據(jù)不同航速,調(diào)整激光(相機(jī))在相鄰兩幀之間的時(shí)間間隔,以保證在此間隔內(nèi)粒子最大位移不超過查詢區(qū)域的1/4,每次測量采集的圖片數(shù)控制在150～200 對之間以滿足數(shù)據(jù)收斂性的要求.試驗(yàn)過程如圖8 所示.

圖7 粒子播撒裝置Fig.7 The particle seeding device

圖8 工作中的SPIV 系統(tǒng)Fig.8 The SIV system in operation

采用Dantec Studio 6.2 軟件對試驗(yàn)圖片進(jìn)行批量處理與后續(xù)分析.自適應(yīng)互相關(guān)算法獲得流場矢量信息,查詢區(qū)域大小設(shè)置為64×64 像素,32×32 像素,16×16 像素,重疊率為50%.鄰域驗(yàn)證設(shè)置為5×5 矩陣,通過比較矩陣內(nèi)的向量并啟用替換,可以有效過濾掉偽向量并用鄰域向量的中值進(jìn)行插值.通過實(shí)時(shí)監(jiān)測迭代步數(shù)對結(jié)果收斂性的影響,迭代次數(shù)選擇為10,檢測最少有效粒子數(shù)設(shè)置為10,最終得到的相鄰速度矢量間距為2.45 mm.

吳鐵成[36]應(yīng)用該SPIV 系統(tǒng)在拖曳水池中進(jìn)行了大量船舶艉流場試驗(yàn),對該系統(tǒng)不確定度進(jìn)行分析,并對KCS 標(biāo)模多工況下伴流場進(jìn)行測量,與公開數(shù)據(jù)進(jìn)行對比,驗(yàn)證了該系統(tǒng)的準(zhǔn)確性.

3 模型試驗(yàn)準(zhǔn)確性驗(yàn)證

應(yīng)用泰勒水池的SUBOFF 實(shí)驗(yàn)數(shù)據(jù)作為參照,對本次模型試驗(yàn)及數(shù)值計(jì)算進(jìn)行驗(yàn)證.由于測試水池以及PIV 測試能力的限制,模型縮尺比略有不同,兩者主尺度參數(shù)見表2.測量航速由Fr相似換算得到,由于此時(shí)浸深較大,自由面影響較小,使用Fr相似來確定本試驗(yàn)航速不夠嚴(yán)謹(jǐn).受試驗(yàn)條件中模型長度及拖曳航速的限制,在今后的試驗(yàn)研究中應(yīng)盡量滿足雷諾數(shù)相似條件,使用更大的模型及航速.

表2 模型差異Table 2 Model differences

提取0.6R圓周范圍內(nèi)若干點(diǎn)上的潛艇槳盤面軸向伴流速度與標(biāo)準(zhǔn)實(shí)驗(yàn)值進(jìn)行對比.取點(diǎn)原則為周向間隔5?取一個(gè)點(diǎn),在0?,90?,180?的±10?鄰域內(nèi)角度間隔加密為2?.其中,周向角0?規(guī)定指向艇體正上方(十二點(diǎn)鐘位置),從艇艉前視,順時(shí)針旋轉(zhuǎn)為正方向.最終,軸向速度使用u/U作無量綱化表示,其中U為拖曳航速.軸向速度對比結(jié)果如圖9 所示.槳盤面上方速度值徑向分布如圖10 所示.

圖9 槳盤面0.6R 處無量綱軸向速度值對比圖Fig.9 Comparison diagram of dimensionless axial velocity value at 0.6R of propeller disk

圖10 槳盤面上方速度徑向分布圖Fig.10 Radial velocity distribution above the propeller disk

PIV 試驗(yàn)得到的軸向速度值普遍小于泰勒水池實(shí)驗(yàn)值,將其歸結(jié)于模型大小引起的尺度效應(yīng).本試驗(yàn)采用的潛艇模型尺度較小,因此相同F(xiàn)r下的Re較小,導(dǎo)致邊界層稍厚,測量得到的速度要小一些.以泰勒水池的實(shí)驗(yàn)值作為參照計(jì)算本試驗(yàn)測量的誤差值,最大誤差在5%以內(nèi).對比發(fā)現(xiàn)CFD 計(jì)算的軸向速度值偏小,其原因?yàn)榇蟛糠志W(wǎng)格位于自由液面區(qū)域,艇體周圍網(wǎng)格只能達(dá)到四五百萬量級,另外模型太小,模型本身產(chǎn)生的流場特征也較小,經(jīng)過一些耗散實(shí)際測量到的量就會(huì)小很多,而數(shù)值模擬條件比較理想化,模擬出與試驗(yàn)環(huán)境相同的真實(shí)物理環(huán)境還存在差異,即CFD 結(jié)果過高估計(jì)了馬蹄渦的強(qiáng)度,使得馬蹄渦的影響區(qū)域增大,造成了槳盤面處偏低的軸向速度.通過以上流場對比結(jié)果來看,本次試驗(yàn)和數(shù)值計(jì)算結(jié)果與泰勒水池標(biāo)準(zhǔn)實(shí)驗(yàn)結(jié)果吻合良好,因此本試驗(yàn)結(jié)果能夠應(yīng)用于后續(xù)的研究.

4 試驗(yàn)結(jié)果分析

4.1 不同潛深下螺旋槳盤面速度場分析

圖11 給出Fr=0.3～0.7 范圍內(nèi)不同潛深工況下槳盤面軸向速度場,左為1.5D潛深結(jié)果,中間為2.5D潛深結(jié)果,右為4D潛深結(jié)果.其中空白區(qū)域?qū)?yīng)槳軸遮擋處的無數(shù)據(jù)區(qū)域.

從圖11 中可以看出,本次PIV 試驗(yàn)大致描繪出了潛艇槳盤面處流場的幾項(xiàng)重要特征.首先,速度值呈現(xiàn)出沿徑向連續(xù)增大的圓圈式形狀;其次,4 個(gè)十字舵角度(0?,90?,180?,270?) 最外側(cè)的速度等值線向外凸出,速度等值線0.8～0.95 范圍表現(xiàn)出了十字舵后端由于速度削減效應(yīng)產(chǎn)生的低速區(qū)的存在;另外,內(nèi)側(cè)的速度等值線朝向槳軸中心凹陷,對應(yīng)速度等值線范圍為0.4～0.8,這是由于主附體交接部位馬蹄渦作用造成動(dòng)量交換,使得中心位置動(dòng)量增大,兩側(cè)動(dòng)量降低,從而形成近似“V”型的速度等值線結(jié)構(gòu).在徑向上表現(xiàn)為馬蹄渦將高動(dòng)量的流體卷攜至徑向內(nèi)部,低動(dòng)量的流體卷攜至徑向外部.大潛深工況下,隨著Fr增大,流場結(jié)構(gòu)特征基本一致.

在1.5D潛深工況下,隨著Fr增大,潛艇尾流場依舊存在上述幾項(xiàng)重要特征.進(jìn)一步比較Fr對軸向速度的影響,可以發(fā)現(xiàn)隨著Fr增大潛艇上方外側(cè)速度等值線有不斷向下的趨勢,最外層0.95 等值線輪廓位于Y/R=1 處.與此同時(shí),上方等值線整體趨于扁平化,由指揮臺圍殼脫落的馬蹄渦而形成的速度等值線凸起越來越不明顯.這是由于潛艇近水面航行時(shí),不同F(xiàn)r下興波高度不同,自由面與艇體伴流相互作用.

圖11 不同潛深工況下槳盤面軸向速度云圖(左側(cè)為1.5D 潛深工況,中間為2.5D 潛深工況,右側(cè)為4D 潛深工況)Fig.11 Axial velocity of propeller disk under different submergence conditions(1.5D submergence depths on the left,2.5D submergence depths in the middle and 4D submergence depths on the right

圖11 不同潛深工況下槳盤面軸向速度云圖(左側(cè)為1.5D 潛深工況,中間為2.5D 潛深工況,右側(cè)為4D 潛深工況)(續(xù))Fig.11 Axial velocity of propeller disk under different submergence conditions(1.5D submergence depths on the left,2.5D submergence depths in the middle and 4D submergence depths on the right(continued)

不同潛深工況下潛艇尾流場流動(dòng)特征基本類似.對比4D潛深、2.5D潛深與1.5D潛深槳盤面軸向速度分布,可以發(fā)現(xiàn)在低Fr時(shí),1.5D潛深的上部分速度分布呈現(xiàn)較為飽和的半圓形分布,4D潛深的上部分速度有壓縮現(xiàn)象,2.5D潛深介于1.5D潛深和4D潛深之間.隨著Fr增大,自由面主要影響流場的上半?yún)^(qū)域,三種潛深工況下的速度均有壓縮趨勢,而1.5D潛深的壓縮現(xiàn)象更為明顯,2.5D潛深次之.在試驗(yàn)過程中,觀察到高速工況對應(yīng)的興波越發(fā)明顯,速度等值線體現(xiàn)出的特征應(yīng)與興波密切相關(guān),自由面產(chǎn)生了一種擠壓作用使得上半?yún)^(qū)域流場結(jié)構(gòu)整體下移.由于2.5D潛深工況整體特性介于1.5D潛深和4D潛深之間,因此下文僅對1.5D潛深和4D潛深工況展開分析.

4.2 不同航速及潛深下潛艇尾部脈動(dòng)速度場分析

由一系列的瞬時(shí)速度云圖分析得到脈動(dòng)速度,并進(jìn)行無量綱化處理,得到uRMS/U脈動(dòng)速度云圖.4D潛深與1.5D潛深工況不同航速下脈動(dòng)速度云圖分別如圖12 和圖13 所示.

從圖12 和圖13 中可以看出,脈動(dòng)速度云圖與軸向速度云圖較為相似,為一圈圈向外延伸的等值線結(jié)構(gòu),不同點(diǎn)在于向外的脈動(dòng)速度值不斷降低,在外流場達(dá)到0.02～0.03 的量級,同時(shí)在槳盤面內(nèi)半徑區(qū)域以及尾舵梢部附近即Y/R=1 出現(xiàn)較大的脈動(dòng)速度,這對應(yīng)于馬蹄渦及梢渦的位置.

對比不同潛深下的脈動(dòng)速度值,可以看出1.5D潛深下的脈動(dòng)速度云圖受自由面影響較大,在局部上體現(xiàn)出脈動(dòng)速度極大值點(diǎn),同時(shí)隨著Fr增大,槳盤面上方速度等值線明顯下移,表明受自由面的影響,較大的脈動(dòng)速度結(jié)構(gòu)整體下移,體現(xiàn)出與軸向速度等值線相似的規(guī)律.

圖12 1.5D 潛深不同航速下脈動(dòng)速度云圖Fig.12 Fluctuating velocity at different speeds of 1.5D depth

圖13 4D 潛深不同航速下脈動(dòng)速度云圖Fig.13 Fluctuating velocity at different speeds of 4D depth

進(jìn)一步提取0.6R范圍內(nèi)的脈動(dòng)速度值,如圖14所示,由于槳軸反光區(qū)域遮擋,60?～120?范圍數(shù)據(jù)缺失.從圖中可以看出不同潛深工況下槳盤面脈動(dòng)速度隨Fr增大所體現(xiàn)出的趨勢是相同的.首先,在每個(gè)90?范圍內(nèi),脈動(dòng)速度呈現(xiàn)出先下降后上升的趨勢,即對應(yīng)每個(gè)尾舵后方都會(huì)出現(xiàn)較大的速度脈動(dòng),具有一定周期性.另外上半部分(0?～90?范圍)的脈動(dòng)速度值略大于下半部分(90?～180?范圍)的脈動(dòng)速度值,這就表明槳軸上方的流場會(huì)受到指揮臺尾流的影響,從而產(chǎn)生更大的脈動(dòng)速度.相比于4D潛深工況,1.5D潛深時(shí)局部上出現(xiàn)脈動(dòng)速度極大值點(diǎn),且不同F(xiàn)r間的脈動(dòng)速度差值也有所增大,因此可知自由面的引入增大了流動(dòng)的脈動(dòng)特性.

4.3 對于PIV 測試結(jié)果的CFD 補(bǔ)充與擴(kuò)展

圖14 槳盤面0.6R 處脈動(dòng)速度對比Fig.14 Fluctuating velocity at 0.6R of propeller disk

由于PIV 測試難以實(shí)現(xiàn)對于潛艇興波的精確探測,僅能夠?qū)Υ笾虏ㄐ芜M(jìn)行觀察以完成粗略的對比,而CFD 計(jì)算可以精確地模擬出波浪的細(xì)節(jié)信息.本小節(jié)基于CFD 模擬結(jié)果,對潛艇近水面航行所興起的波浪進(jìn)行量化分析,并給出中縱切面的波浪高度,從而在興波的角度對潛艇繞流場隨航速與潛深的變化規(guī)律進(jìn)行解釋.

通過數(shù)值模擬發(fā)現(xiàn),在4D潛深時(shí),興波高度隨航速增大不斷提高,但不同航速對應(yīng)的波浪峰值差距僅為毫米量級,完全可以忽略自由表面的影響,因此下文僅對1.5D潛深工況進(jìn)行分析.

從圖15 可以看出,不同航速的潛艇興波具有以下規(guī)律:首先潛艇首部后方存在一個(gè)波峰,且航速越高,波幅越大.由于波長隨Fr增大而增大,在航速增大的情況下,興波波峰的起始位置會(huì)向艇后方移動(dòng).同時(shí)在潛艇尾肩部存在一波谷,隨著航速變大,波谷位置同樣向后移.因此在槳盤面位置總體表現(xiàn)為隨航速增大而降低的液面高度.

圖15 1.5D 潛深不同速度下潛艇中縱剖面波浪高度對比Fig.15 Comparison of wave height of longitudinal section in centre plane of submarine at different speeds of 1.5D

給出不同潛深及不同航速下的潛艇中縱切面速度云圖,如圖16 所示,左為1.5D潛深工況,右為4D潛深工況.

圖16 不同航速下潛艇中縱剖面速度云圖(左為1.5D,右為4D)Fig.16 Velocity nephogram of longitudinal section in centre plane of submarine at different speeds(1.5D submergence depths on the left and 4D submergence depths on the right)

通過對比不同潛深工況下速度云圖發(fā)現(xiàn),自由液面存在對流場的最大影響體現(xiàn)在上半?yún)^(qū)域.同一航速下,近水面航行時(shí),艇體與自由面間流場速度明顯增大,特別在波谷區(qū)域(對應(yīng)槳盤面區(qū)域),流場速度明顯提升.這可以由文丘里效應(yīng)合理地進(jìn)行解釋,即受限流動(dòng)在經(jīng)過縮小的過流斷面時(shí),流體流速會(huì)增大.由前文可知,隨著Fr提升,槳盤面處的液面高度越來越低,這就導(dǎo)致了槳盤面位置出現(xiàn)更大的流體速度,而流速的增大勢必會(huì)導(dǎo)致伴流場結(jié)構(gòu)的變化.

5 結(jié)論

本文應(yīng)用SPIV 系統(tǒng)對潛艇尾部流場進(jìn)行測試,同時(shí)輔以數(shù)值模擬對試驗(yàn)結(jié)果進(jìn)行補(bǔ)充與擴(kuò)展.通過對不同潛深工況不同航速的分析得到以下結(jié)論:

(1)在1.5D潛深工況下,隨著Fr增大,由于自由面產(chǎn)生的擠壓作用使得槳盤面上方速度等值線有不斷向下的趨勢,上方等值線整體趨于扁平化,由指揮臺圍殼脫落的馬蹄渦形成的速度等值線凸起越來越不明顯.

(2)每個(gè)尾舵后方都會(huì)出現(xiàn)較大的脈動(dòng)速度,且具有一定周期性,在周期內(nèi),整體脈動(dòng)速度呈現(xiàn)先下降后上升的趨勢.受指揮臺尾流影響,槳盤面上方流場體現(xiàn)出更大的脈動(dòng)速度.自由面會(huì)增大流體的脈動(dòng)特性,在1.5D潛深工況的脈動(dòng)速度云圖中存在局部脈動(dòng)速度極大值點(diǎn),同時(shí)隨著Fr增大,槳盤面上方0.04 以上的速度等值線明顯下移,表明受自由面的影響,較大的脈動(dòng)速度結(jié)構(gòu)整體下移.

(3) 由于波長隨Fr增大而變大,隨著航速增大,興波波峰的起始位置向艇后方移動(dòng).同時(shí)在潛艇尾肩部存在一波谷,隨著航速增大大,波谷位置同樣向后移.因此在槳盤面處總體表現(xiàn)為自由面高度隨航速增大而降低.

(4)自由面存在時(shí),艇體與自由面間流場速度明顯增大,特別是波谷區(qū)域,也恰好對應(yīng)槳盤面區(qū)域,流場速度明顯提升.隨著Fr增大,槳盤面處的液面高度越來越低,這就導(dǎo)致了槳盤面位置出現(xiàn)更大的流體速度,即造成了槳盤面伴流場擠壓現(xiàn)象的發(fā)生.