對(duì)轉(zhuǎn)槳推進(jìn)的高速水下航行體實(shí)尺度自航計(jì)算與分析

黃永生，楊晨俊，董小倩

1 上海交通大學(xué) 海洋工程國家重點(diǎn)實(shí)驗(yàn)室，上海 200240

2 高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心，上海 200240

0 引 言

快速性預(yù)報(bào)是水面船舶及水下航行體設(shè)計(jì)的重要環(huán)節(jié)之一，其傳統(tǒng)做法是在阻力、敞水及自航模型試驗(yàn)的基礎(chǔ)上對(duì)阻力、伴流及敞水性能的尺度效應(yīng)采用1978 ITTC規(guī)程等經(jīng)驗(yàn)方法進(jìn)行修正。對(duì)于常規(guī)螺旋槳推進(jìn)的船舶，采用該方法雖然可靠性好，但成本較高。高速水下航行體一般采用對(duì)轉(zhuǎn)槳（CRP）或泵噴推進(jìn)，推進(jìn)器部件間的相互作用給尺度效應(yīng)修正帶來了一定的困難（主要是伴流）。另一方面，隨著計(jì)算流體力學(xué)（CFD）的快速發(fā)展，粘流CFD計(jì)算技術(shù)在快速性預(yù)報(bào)方面的應(yīng)用日益增多，并有逐步取代模型試驗(yàn)的趨勢(shì)。目前，阻力和敞水性能的數(shù)值計(jì)算精度已基本達(dá)到工程實(shí)用要求，自航計(jì)算的精度也在逐步提高。與模型試驗(yàn)相比，粘流CFD計(jì)算技術(shù)的另一個(gè)優(yōu)勢(shì)在于能夠進(jìn)行實(shí)尺度模擬，從而無需進(jìn)行尺度效應(yīng)修正。實(shí)尺度模擬對(duì)計(jì)算機(jī)容量和速度的要求較高，其主要困難在于船體，而推進(jìn)器則相對(duì)容易實(shí)現(xiàn)。

自航計(jì)算需要同時(shí)模擬船體與推進(jìn)器，網(wǎng)格量和計(jì)算量均較大；從理論上看，需要進(jìn)行非定常模擬，但所需計(jì)算時(shí)間會(huì)更長(zhǎng)。因此，一般采用定?；驕?zhǔn)定常的方法模擬船體與推進(jìn)器的相互作用。為了降低網(wǎng)格生成難度，節(jié)省計(jì)算時(shí)間，進(jìn)行定常計(jì)算時(shí)可以采用體積力模型來代替螺旋槳的作用。體積力模型分描述型和迭代型2種，前者的體積力徑向分布形式是基于經(jīng)驗(yàn)公式，后者則是通過螺旋槳升力面、面元等勢(shì)流方法計(jì)算得到?；诿枋鲂腕w積力模型，Choi等［1］對(duì)VLCC進(jìn)行了自航模擬和快速性預(yù)報(bào)；傅慧萍［2］針對(duì)KCS進(jìn)行了帶自由液面的自航模擬，研究了扭矩對(duì)計(jì)算精度的影響；呂曉軍等［3］開展了潛艇自航計(jì)算研究。Kim等［4］的研究表明，由于迭代型體積力模型是根據(jù)槳葉的真實(shí)幾何及實(shí)效伴流分布來計(jì)算體積力的分布，所以它對(duì)艉部流場(chǎng)的模擬比描述型體積力模型的結(jié)果更加準(zhǔn)確。雖然采用體積力方法對(duì)自航因子的預(yù)報(bào)結(jié)果較好，但畢竟是用力場(chǎng)代替了真實(shí)槳葉，對(duì)于槳葉對(duì)周圍流體的排擠效應(yīng)和擾動(dòng)難以精確體現(xiàn)，從而限制了計(jì)算精度和流場(chǎng)細(xì)節(jié)。2005年以后，采用螺旋槳真實(shí)幾何的數(shù)值自航模擬方法逐步發(fā)展起來。Lübke［5］應(yīng)用CFX軟件對(duì)KCS進(jìn)行了自航模擬；Choi等［6］對(duì)一系列船型進(jìn)行數(shù)值自航模擬和快速性預(yù)報(bào)，充分驗(yàn)證了其計(jì)算方法的精度和可靠性；熊鷹等［7］則針對(duì)傅汝德數(shù)較低的排水型船舶提出了一種簡(jiǎn)化的自由面處理方法，用以加快計(jì)算速度；楊琴等［8］針對(duì)配有七葉槳的SUBOFF潛艇進(jìn)行了自航模擬，并分析給出了自航因子。上述自航模擬均是在模型尺度下進(jìn)行，仍需通過修正尺度效應(yīng)來預(yù)報(bào)實(shí)船快速性。隨著計(jì)算機(jī)硬件能力的提高，近年來，有關(guān)實(shí)尺度自航計(jì)算的研究逐步增多。Castro等［9］以KCS為對(duì)象開展了實(shí)尺度自航模擬，并研究了推進(jìn)因子的尺度效應(yīng)；李亮等［10］也對(duì)KCS開展了類似的研究。

目前，針對(duì)高速水下航行體推進(jìn)性能方面的數(shù)值計(jì)算研究相對(duì)較少。何文生［11］較早地采用螺旋槳的真實(shí)幾何對(duì)高速水下航行體進(jìn)行了RANS模擬，但受計(jì)算機(jī)能力的限制，簡(jiǎn)化了附體。張濤等［12-13］針對(duì)對(duì)轉(zhuǎn)槳的敞水性能，開展了RANS計(jì)算研究。但是，針對(duì)高速水下航行體的自航計(jì)算研究，在模型和實(shí)尺度方面均未見公開發(fā)表的文獻(xiàn)。

為此，本文擬開展對(duì)轉(zhuǎn)槳推進(jìn)的高速水下航行體實(shí)尺度自航CFD計(jì)算研究，基于阻力、敞水及自航模擬結(jié)果，分析自航因子。將計(jì)算處于深潛狀態(tài)的航行體，而不考慮自由液面和空化的影響。

1 計(jì)算方法

本文采用求解RANS方程的方法，來開展對(duì)轉(zhuǎn)槳推進(jìn)的高速水下航行體自航性能實(shí)尺度計(jì)算研究。假定流體不可壓縮，流動(dòng)為全湍流，采用標(biāo)準(zhǔn)k-ε湍流模型和非平衡壁面函數(shù)?？刂品匠贪ㄟB續(xù)性方程、動(dòng)量輸運(yùn)方程及湍流輸運(yùn)方程，此處不再贅述。網(wǎng)格生成采用ICEMCFD及GAM?BIT軟件，數(shù)值計(jì)算采用FLUENT軟件。

1.1 計(jì)算域及其離散

如圖1所示（圖中D0為航行體平行中體段的直徑；L為航行體長(zhǎng)度），阻力及自航計(jì)算域?yàn)榘鼑叫畜w的回轉(zhuǎn)體，但阻力計(jì)算域不包含前、后槳葉。將計(jì)算域劃分為4個(gè)子域：Fore-body和Rudder/fin為靜止域；Rotor-F和Rotor-A在阻力計(jì)算中為靜止域，在自航計(jì)算中為旋向相反的旋轉(zhuǎn)域。4個(gè)子域之間依次設(shè)有Int.mid，Int.aft和Int.CRP三對(duì)交界面。

子域Fore-body和Rudder/fin采用結(jié)構(gòu)化單元離散，對(duì)航行體頭部及附體等幾何變化劇烈的區(qū)域進(jìn)行了網(wǎng)格加密。用圓柱面CS將子域Rotor-F和Rotor-A分割為2個(gè)部分，其中CS的外部用結(jié)構(gòu)化單元離散；用通過尾端面的橫截面XS進(jìn)一步分割子域Rotor-A，其中CS內(nèi)部和XS下游部分用棱柱單元離散。在自航計(jì)算中，CS內(nèi)部以及包含前、后槳葉的區(qū)域采用四面體及五面體非結(jié)構(gòu)化單元離散；在阻力計(jì)算中，原前、后槳葉所在區(qū)域用結(jié)構(gòu)化單元離散。圖2所示為阻力、自航計(jì)算關(guān)鍵區(qū)域的物面/截面部分區(qū)域面網(wǎng)格。

如圖3所示，對(duì)轉(zhuǎn)槳的敞水計(jì)算采用圓柱形全流道計(jì)算域，圖中，DF為前槳直徑。前、后槳之間的交界面Int.CRP將計(jì)算域分為了2個(gè)旋轉(zhuǎn)方向相反的子域，每個(gè)旋轉(zhuǎn)域又進(jìn)一步劃分為3個(gè)子域。其中子域Far-F和Far-A采用六面體結(jié)構(gòu)化單元離散；包圍前槳及轂帽的子域Rotor-F和包圍后槳的子域Rotor-A采用四面體及五面體非結(jié)構(gòu)化單元離散；子域UpS和DownS采用棱柱體單元離散。為了保證槳葉間網(wǎng)格的周期性，對(duì)于所有子域，在相鄰槳葉之間生成周期面，將全流道沿周向等分成Z（Z為槳葉數(shù)）個(gè)子流道，通過復(fù)制子流道網(wǎng)格、合并周期面網(wǎng)格，得到全流道網(wǎng)格。圖4所示為子域UpS，Rotor-F，Rotor-A及DownS的子流道幾何。

需要說明的是，本文的模型尺度及實(shí)尺度的計(jì)算均采用壁面函數(shù)，因此，需要根據(jù)計(jì)算對(duì)象的雷諾數(shù)，對(duì)航行體及槳葉、槳轂表面第1層單元采用合適的高度，以使壁面y+值處于可適用壁面函數(shù)的范圍，y+的具體數(shù)值見各算例。在阻力驗(yàn)證計(jì)算中，模型尺度與實(shí)尺度的計(jì)算網(wǎng)格拓?fù)潢P(guān)系以及各子域的單元類型均完全相同；而在敞水驗(yàn)證計(jì)算中，實(shí)尺度的計(jì)算網(wǎng)格由模型尺度網(wǎng)格按縮尺比（1∶1.5）放大得到。

1.2 計(jì)算設(shè)置

將航行體、推進(jìn)器表面分別設(shè)為靜止和旋轉(zhuǎn)坐標(biāo)系中的光滑無滑移壁面；計(jì)算域的下游邊界面設(shè)為壓力出口，其余邊界面設(shè)為速度入口。動(dòng)量方程和湍流輸運(yùn)方程采用二階迎風(fēng)格式離散，壓力方程采用標(biāo)準(zhǔn)格式離散，壓力與速度的耦合采用SIMPLE算法。

阻力計(jì)算采用定常模型；敞水及自航計(jì)算分別采用準(zhǔn)定常模型和非定常（滑移網(wǎng)格）模型，在非定常計(jì)算的每個(gè)時(shí)間步前、后槳葉各轉(zhuǎn)動(dòng)1°。

實(shí)尺度自航計(jì)算在給定的航速下進(jìn)行，通過調(diào)整轉(zhuǎn)速（前、后槳相等），直至螺旋槳推力與航行體阻力之差小于后者的1%，即認(rèn)為達(dá)到了航行體的自航點(diǎn)，結(jié)束計(jì)算。在自航的非定常計(jì)算中，推力、阻力均為一個(gè)槳葉旋轉(zhuǎn)周期（360°）的時(shí)間平均值。敞水和自航的非定常計(jì)算收斂準(zhǔn)則是，推力在連續(xù)2個(gè)旋轉(zhuǎn)周期內(nèi)的時(shí)間平均值之差小于0.5%。

2 計(jì)算精度驗(yàn)證

基于類似的航行體及對(duì)轉(zhuǎn)槳的模型試驗(yàn)結(jié)果，分別驗(yàn)證阻力和敞水性能的數(shù)值計(jì)算精度；同時(shí)，進(jìn)行實(shí)尺度計(jì)算，定性分析相關(guān)結(jié)果的合理性。因缺乏試驗(yàn)數(shù)據(jù)，對(duì)高速水下航行體的自航計(jì)算精度目前還無法驗(yàn)證。

2.1 航行體阻力

以圖5所示的航行體為對(duì)象，驗(yàn)證阻力計(jì)算精度。在網(wǎng)格依賴性檢驗(yàn)的基礎(chǔ)上，先進(jìn)行模型尺度的計(jì)算，對(duì)比試驗(yàn)數(shù)據(jù)以驗(yàn)證計(jì)算精度；然后，進(jìn)行實(shí)尺度計(jì)算，并根據(jù)模型試驗(yàn)結(jié)果對(duì)實(shí)尺度阻力進(jìn)行預(yù)報(bào)，間接驗(yàn)證實(shí)尺度阻力的RANS計(jì)算精度。

2.1.1 網(wǎng)格依賴性檢驗(yàn)

為保證壁面y+值基本相同，采用相同的單元徑向尺度，沿航行體軸向和周向按2∶2∶1的尺度比由疏到密建立3套網(wǎng)格M1，M2和M3，然后在模型尺度下進(jìn)行了2個(gè)航速的阻力計(jì)算，雷諾數(shù)Rem分別為1.220×107和2.034×107，壁面y+的平均值均約為50。表1所示為3套網(wǎng)格參數(shù)的比較。航行體總阻力系數(shù)Ctm隨網(wǎng)格尺度的變化如圖6所示。隨著網(wǎng)格的加密，Ctm趨于收斂；M2網(wǎng)格與M3網(wǎng)格Ctm的相對(duì)差僅為0.4%。因此，后續(xù)阻力計(jì)算采用M2網(wǎng)格的尺度，以兼顧計(jì)算精度和效率。

表1 3套阻力計(jì)算網(wǎng)格參數(shù)的比較Table 1 Comparison of the parameters for three sets of grids of resistance computation

2.1.2 模型和實(shí)尺度阻力計(jì)算精度驗(yàn)證

首先，驗(yàn)證模型尺度的阻力計(jì)算精度。表2所示為RANS計(jì)算結(jié)果與模型試驗(yàn)結(jié)果的比較。由于附體與殼體相比面積很小，因此分析形狀因子時(shí)未區(qū)分殼體與附體，即

對(duì)模型試驗(yàn)，摩擦阻力系數(shù)Cfm按1957 ITTC公式計(jì)算：

對(duì)于RANS計(jì)算，Cfm直接采用計(jì)算結(jié)果。在計(jì)算雷諾數(shù)的范圍內(nèi)，模型總阻力計(jì)算值總體上低于試驗(yàn)值，但誤差不超過3%；摩擦阻力的計(jì)算誤差不超過3.5%，且隨著航速的增高誤差減小，這可能是因?yàn)橛?jì)算是在全湍流狀態(tài)下進(jìn)行，而模型試驗(yàn)雖然采取了激流措施，但在航速較低時(shí)仍受到層流及過渡區(qū)的影響；形狀因子的計(jì)算誤差不超過4%。

實(shí)尺度網(wǎng)格拓?fù)渑c模型尺度相同，計(jì)算單元總數(shù)為1 985萬。雷諾數(shù)Res=1.682×108，壁面y+平均值約為50。表3所示為實(shí)尺度航行體阻力的RANS計(jì)算結(jié)果與基于模型試驗(yàn)的預(yù)報(bào)結(jié)果間的比較，用于對(duì)實(shí)尺度阻力計(jì)算精度的間接驗(yàn)證。在預(yù)報(bào)中，忽略了形狀因子的尺度效應(yīng)，實(shí)尺度的形狀因子(1+k)s為表2中模型試驗(yàn)雷諾數(shù)范圍內(nèi)(1+k)m的算術(shù)平均值，Cfs按1957 ITTC公式計(jì)算，Cts等于(1+k)s與Cfs的乘積，為航行體總阻力系數(shù)。在RANS計(jì)算結(jié)果分析中，Cfs和Cts均為直接計(jì)算值，(1+k)s為Cts與Cfs的比值。由表3可見，實(shí)尺度航行體總阻力的RANS計(jì)算值略低于基于模型試驗(yàn)的預(yù)報(bào)值，誤差為3.3%。

表2 模型尺度航行體阻力RANS計(jì)算與試驗(yàn)的比較Table 2 Comparison of RANS-simulated model-scale resistance of the underwater vehicle with experimental data

表3 實(shí)尺度航行體阻力的RANS計(jì)算結(jié)果與基于模型試驗(yàn)的預(yù)報(bào)結(jié)果間的比較Table 3 Comparison between RANS-simulated full-scale resistance ofthe underwater vehicle and that predicted from model test results

2.2 對(duì)轉(zhuǎn)槳敞水性能

以某高速航行體對(duì)轉(zhuǎn)槳為對(duì)象，驗(yàn)證敞水性能的RANS計(jì)算精度。該對(duì)轉(zhuǎn)槳前槳11葉、后槳9葉，前、后槳的直徑之比為1.062。首先，進(jìn)行模型尺度的網(wǎng)格依賴性分析，確定合適的網(wǎng)格尺度；然后，進(jìn)行模型尺度和實(shí)尺度敞水性能的計(jì)算與試驗(yàn)比較。

2.2.1 網(wǎng)格依賴性檢驗(yàn)

計(jì)算在模型尺度進(jìn)行，進(jìn)速系數(shù)J=1.0和1.4時(shí)，如式（3）所示的0.7R（R為槳直徑）處雷諾數(shù)分別為8.45×105和9.05×105。

式中：VA為進(jìn)速；b0.7R，n和D分別為槳葉0.7R處的弦長(zhǎng)、轉(zhuǎn)速與直徑；ν為水的運(yùn)動(dòng)粘性系數(shù)。由疏到密建立了3套網(wǎng)格G1，G2和G3，其單元尺度之比為2∶2∶1，面單元和體單元的密度分布及增長(zhǎng)率保持不變。表4所示為3套網(wǎng)格的參數(shù)比較，其中面單元相對(duì)尺度為面單元尺度與前槳直徑的比值。

表5所示為網(wǎng)格尺度不同引起推力系數(shù)KT和扭矩系數(shù)KQ計(jì)算結(jié)果的相對(duì)變化。表中KT和KQ均以前槳的直徑無量綱化，下標(biāo)F和A分別表示前、后槳。分別用R1，R2和R3代表G1，G2，G3的KT或KQ計(jì)算結(jié)果，則表中C21=R2/R1-1，表示R2相對(duì)R1的變化，C32=R3/R2-1，表示R3相對(duì)R2的變化。從表中可以看出，前、后槳的推力和扭矩基本上是隨網(wǎng)格尺度的減小而減小；3套網(wǎng)格計(jì)算結(jié)果之差均在1.2%以內(nèi)。經(jīng)權(quán)衡精度與效率，后續(xù)的計(jì)算采用G2網(wǎng)格的尺度。

表4 3套敞水計(jì)算網(wǎng)格參數(shù)的比較Table 4 Comparison of the parameters for three sets of grids of open-water computation

表5 敞水性能計(jì)算結(jié)果隨網(wǎng)格尺度的相對(duì)變化Table 5 Influence of grid size on computed open-water performance

2.2.2 敞水性能計(jì)算

用于驗(yàn)證的對(duì)轉(zhuǎn)槳模型與實(shí)槳的尺度比為1∶1.5，實(shí)槳的計(jì)算網(wǎng)格由槳模網(wǎng)格放大得到。在計(jì)算工況范圍內(nèi)，槳模與實(shí)槳在0.7R處的雷諾數(shù)分別為8.45×105～9.05×105和3.17×106～3.39×106，槳葉表面的y+平均值分別為35～50和130～200，均在壁面函數(shù)的適用范圍內(nèi)。當(dāng)前、后槳的葉數(shù)相互不構(gòu)成整倍數(shù)關(guān)系時(shí)，兩者相互作用的非定常性很弱［14］，非定常計(jì)算結(jié)果的時(shí)間平均值與準(zhǔn)定常結(jié)果很接近，因此采用準(zhǔn)定常計(jì)算以節(jié)約計(jì)算時(shí)間。

圖7所示為敞水性能的模型和實(shí)尺度RANS計(jì)算結(jié)果與模型試驗(yàn)結(jié)果比較，圖中η0為敞水效率，下標(biāo)F和A分別表示前、后槳。在計(jì)算進(jìn)速系數(shù)范圍內(nèi)，總推力系數(shù)和總扭矩系數(shù)的計(jì)算誤差分別小于2%和4%，表明模型尺度RANS計(jì)算具有較高的精度。比較實(shí)尺度與模型尺度的RANS計(jì)算結(jié)果發(fā)現(xiàn)，KT的變化很小，而KQ則有較明顯的降低，導(dǎo)致實(shí)槳的敞水效率比模型槳的高2.5%左右，該結(jié)果與單槳敞水性能尺度效應(yīng)定性一致，表明對(duì)轉(zhuǎn)槳實(shí)尺度敞水性能的RANS計(jì)算方法比較可靠。

3 實(shí)尺度自航計(jì)算與分析

基于上述計(jì)算方法，對(duì)采用對(duì)轉(zhuǎn)槳推進(jìn)的某高速水下航行體進(jìn)行實(shí)尺度阻力、敞水及自航計(jì)算，并分析其推進(jìn)因子。對(duì)轉(zhuǎn)槳的前槳9葉、右旋，后槳7葉、左旋，二者轉(zhuǎn)速相同。

3.1 航行體阻力計(jì)算

阻力計(jì)算單元總數(shù)約1 127萬。計(jì)算航速范圍為設(shè)計(jì)航速的64%～109%，雷諾數(shù)為8.07×107～13.83×107，航行體壁面y+的平均值約為50。圖8所示為該航行體實(shí)尺度阻力的RANS計(jì)算模型，圖9所示為拖航阻力系數(shù)計(jì)算結(jié)果。圖10所示為航行體摩擦阻力系數(shù)的RANS結(jié)果與1957 ITTC公式計(jì)算結(jié)果的比較，由圖可見，前者略高，兩者間的差異隨航速的增大而減小，在設(shè)計(jì)航速下相差1.9%。圖11所示為形狀因子(1+k)s隨航速的變化，雖然形狀因子隨航速的增高而變大，但在計(jì)算航速范圍內(nèi)其差值小于1%。

3.2 對(duì)轉(zhuǎn)槳敞水性能計(jì)算

敞水計(jì)算單元總數(shù)約為4 634萬，槳葉0.7R處雷諾數(shù)范圍為2.38×106～2.57×106，槳葉表面y+的平均值約為210。分別采用準(zhǔn)定常和非定常方法模擬前、后槳之間的相互作用，其中準(zhǔn)定常方法用于模擬槳葉旋轉(zhuǎn)過程中某一瞬間的流動(dòng)，而非定常方法則是在時(shí)域中模擬槳葉旋轉(zhuǎn)過程中的流動(dòng)。圖12和圖13所示分別為實(shí)尺度對(duì)轉(zhuǎn)槳敞水性能的RANS計(jì)算模型和計(jì)算結(jié)果，其中非定常方法模擬的結(jié)果為槳葉旋轉(zhuǎn)一周的平均值，圖中扭矩系數(shù)KQ的第2下標(biāo)0表示敞水，以與后文中的自航數(shù)據(jù)相區(qū)別。從圖13中可以看出，采用準(zhǔn)定常方法與非定常方法模擬的敞水性能結(jié)果非常接近，差值在1%以內(nèi)。

3.3 自航計(jì)算結(jié)果

自航計(jì)算在設(shè)計(jì)航速下進(jìn)行，計(jì)算單元總數(shù)為4 827萬，航行體雷諾數(shù)為12.69×107，壁面y+的平均值約為50。分別采用準(zhǔn)定常和非定常方法進(jìn)行計(jì)算，自航轉(zhuǎn)速根據(jù)1.2節(jié)所述方法確定。圖14所示為自航的RANS計(jì)算模型。為了節(jié)約計(jì)算時(shí)間，非定常計(jì)算采用準(zhǔn)定常計(jì)算得到的自航轉(zhuǎn)速及流場(chǎng)作為初始值。圖15所示為非定常計(jì)算的軸向力收斂歷程。圖中，KQBF和KQBA分別為前、后槳的扭矩系數(shù)。其中推力系數(shù)和扭矩系數(shù)的脈動(dòng)幅值分別為各自平均值的0.2%與0.1%，而Cts的脈動(dòng)幅值則約為平均值的2.5%。Cts的脈動(dòng)主要源于由槳葉尾流引起的航行體尾端面壓力脈動(dòng)。

自航的非定常與準(zhǔn)定常計(jì)算結(jié)果間的比較見表6，表中Ns為航行體推力和阻力平衡時(shí)前、后槳的轉(zhuǎn)速（兩者相等）。由于航行體伴流不均勻性的影響，自航非定常計(jì)算與準(zhǔn)定常計(jì)算結(jié)果間的差別比敞水情況下的大，主要表現(xiàn)為非定常計(jì)算的后槳負(fù)荷比準(zhǔn)定常的高。從理論上講，非定常計(jì)算模型更接近于物理實(shí)際，因此在前、后槳扭矩平衡要求較高的設(shè)計(jì)場(chǎng)合，采用非定常計(jì)算相對(duì)可靠。

表6 自航計(jì)算的準(zhǔn)定常與非定常計(jì)算結(jié)果的相對(duì)差Table 6 Relative difference between quasi-steady and unsteady results of self-propulsion simulation

圖16所示為自航狀態(tài)基于Q準(zhǔn)則的尾部渦流形態(tài)（準(zhǔn)定常結(jié)果）。從圖中可以看出，在槳葉邊緣、附體的尾緣和航行體尾端面均出現(xiàn)了明顯的渦流，較合理地反映了槳葉對(duì)周圍流場(chǎng)的擾動(dòng)。由于前、后槳旋轉(zhuǎn)域間交界面上的網(wǎng)格無法做到相同，因此前槳尾渦在進(jìn)入后槳域后因網(wǎng)格耗散幾乎消失了。

圖17所示為在尾端面下游0.05倍前槳直徑處，橫截面內(nèi)切向速度沿半徑的分布。其中，切向速度為沿圓周線的平均值，并用當(dāng)?shù)匕霃綐~的旋轉(zhuǎn)線速度進(jìn)行無量綱化。所謂當(dāng)?shù)匕霃?，是指在航行體中縱剖面上，將通過前、后槳盤半徑外端點(diǎn)的線段延長(zhǎng)到該縱向位置，線段終點(diǎn)與航行體軸線的距離。圖17中的橫坐標(biāo)為徑向相對(duì)坐標(biāo)2r/DF，其中r為徑向坐標(biāo)值。0.65R至葉梢的切向速度幾乎為零，但內(nèi)半徑的切向速度仍有殘留。總體來看，在設(shè)計(jì)航速下，該航行體尾流基本沒有旋轉(zhuǎn)，表明前、后槳的扭矩平衡較好。

3.4 自航因子分析

綜合阻力、敞水及自航計(jì)算結(jié)果，采用等推力法進(jìn)行自航因子分析，結(jié)果如表7所示。表中：J0為根據(jù)KT在圖13對(duì)應(yīng)的推力系數(shù)曲線上插值得到的進(jìn)速系數(shù)；η0為J0所對(duì)應(yīng)的效率曲線上的值；KQ的第2下標(biāo)B表示自航狀態(tài)，用于與敞水狀態(tài)的0相區(qū)別；ηR為相對(duì)旋轉(zhuǎn)效率；ηD為推進(jìn)效率；下標(biāo)F和A分別表示前槳和后槳。

表7 自航分析結(jié)果Table 7 Results of self-propulsion analysis

分析比較表明：

1）自航因子計(jì)算值處于合理的范圍。參考文獻(xiàn)［15］，水下航行體的伴流分?jǐn)?shù)w和推力減額系數(shù)（1-t）的范圍分別為0.10～0.25和0.82～0.90，本文中航行體的計(jì)算結(jié)果也處于此范圍內(nèi)。

2）前槳的伴流分?jǐn)?shù)高、后槳的伴流分?jǐn)?shù)低，分析原因認(rèn)為主要為：一是前槳更靠近尾附體，附體的尾流速度虧缺量更大；二是前、后槳轂徑之比為1.31，而直徑之比僅為1.06，這就意味著前槳盤面有較多的面積處于槳轂邊界層中。此外，槳葉的抽吸作用會(huì)導(dǎo)致槳盤面處的實(shí)效伴流低于標(biāo)稱伴流，如果這種作用在槳盤前較弱、槳盤后較強(qiáng)，也會(huì)成為原因之一，不過該問題還有待進(jìn)一步的研究。

3）推進(jìn)因子及推進(jìn)效率的準(zhǔn)定常與非定常計(jì)算結(jié)果間的誤差在2%以內(nèi)。若采用準(zhǔn)定常計(jì)算方法得到的推力減額較低、推進(jìn)效率較高，則有可能導(dǎo)致設(shè)計(jì)槳偏重，預(yù)報(bào)的航速偏高。

4 結(jié) 論

本文通過求解RANS方程，對(duì)采用對(duì)轉(zhuǎn)槳推進(jìn)的高速水下航行體開展了實(shí)尺度自航計(jì)算與分析研究。通過與模型試驗(yàn)結(jié)果的比較，驗(yàn)證了在航行體阻力和對(duì)轉(zhuǎn)槳敞水性能方面，采用本文方法計(jì)算精度良好；間接的驗(yàn)證及定性分析表明，實(shí)尺度阻力和敞水的計(jì)算結(jié)果合理、計(jì)算方法可靠，但實(shí)尺度自航計(jì)算仍有待驗(yàn)證。對(duì)某航行體實(shí)尺度的自航計(jì)算與分析表明：

1）采用準(zhǔn)定常和非定常方法計(jì)算得到的推進(jìn)因子其數(shù)值比較接近，也在合理范圍內(nèi)，可望為高速水下航行體對(duì)轉(zhuǎn)槳設(shè)計(jì)提供較可靠的輸入，從而提高設(shè)計(jì)精度、減少模型試驗(yàn)，縮短設(shè)計(jì)周期。

2）在自航模擬中，準(zhǔn)定常計(jì)算方法具有精度合理、計(jì)算量小的優(yōu)點(diǎn)，比較適合工程應(yīng)用；但當(dāng)航速預(yù)報(bào)精度以及前、后槳的扭矩平衡精度要求較高時(shí)，采用非定常計(jì)算方法更合理。