舒敏驊，陳科，尤云祥，胡天群，劉恒

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

動(dòng)網(wǎng)格方法在螺旋槳非定常軸承力數(shù)值計(jì)算中的應(yīng)用

舒敏驊1，2，陳科1，2，尤云祥1，2，胡天群1，2，劉恒1，2

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

為了探究螺旋槳軸承力的非定常變化特性，研究了動(dòng)網(wǎng)格（DM）方法在螺旋槳非定常軸承力（UBF）數(shù)值計(jì)算中的適用性問題。首先，通過3套不同質(zhì)量網(wǎng)格的數(shù)值計(jì)算比較，進(jìn)行網(wǎng)格無關(guān)性檢驗(yàn)，選出合適的計(jì)算網(wǎng)格。進(jìn)一步采用所選網(wǎng)格DM數(shù)值模擬，并將結(jié)果與滑移網(wǎng)格（SM）方法、勢(shì)流（PM）方法及試驗(yàn)方法得到的結(jié)果進(jìn)行對(duì)比。分析結(jié)果顯示了DM方法得到的螺旋槳UBF結(jié)果與SM方法以及試驗(yàn)的結(jié)果吻合較好，且優(yōu)于PM方法。同時(shí)，DM方法與SM方法得到的流場分布也幾乎一致。研究表明：DM方法應(yīng)用于螺旋槳UBF數(shù)值計(jì)算是可行的。

艦船；三葉常規(guī)螺旋槳；非定常軸承力；動(dòng)網(wǎng)格；滑移網(wǎng)格；數(shù)值計(jì)算

網(wǎng)絡(luò)出版地址：http://www.cnki.net/kcms/detail/42.1755.tj.20160531.1104.010.html期刊網(wǎng)址：www.ship-research.com

引用格式：舒敏驊，陳科，尤云祥，等.動(dòng)網(wǎng)格方法在螺旋槳非定常軸承力數(shù)值計(jì)算中的應(yīng)用［J］.中國艦船研究，2016，11（3）：25-31.

SHU Minhua，CHEN Ke，YOU Yunxiang，et al.Application of the dynamic mesh method in the numerical simulation of unsteady bearing forces of a propelle［rJ］.Chinese Journal of Ship Research，2016，11（3）：25-31.

0　引　言

近年來，隨著艦船向大型化和快速化方向發(fā)展，對(duì)艦船的振動(dòng)和噪聲性能提出了嚴(yán)峻的挑戰(zhàn)。大量的實(shí)船測試結(jié)果表明，造成艦船總聲壓級(jí)指標(biāo)過高的主要原因是1/3倍頻程聲壓譜和線譜中低頻段的若干線譜或窄帶輻射噪聲的聲壓過高，而這些主要特征與艦船艉部的激振力特性有著密切的關(guān)系。因此，從低頻螺旋槳激振力的特性入手，尋求控制螺旋槳—軸系—艦船振動(dòng)輻射噪聲的有效方法已成為船舶海洋工程領(lǐng)域的一項(xiàng)重要研究課題。

螺旋槳的激振力可分為非定常軸承力（Unsteady Bearing Forces，UBF）和螺旋槳的誘導(dǎo)脈動(dòng)壓力，其中工作于艦船艉部非均勻伴流場中的螺旋槳產(chǎn)生的UBF是艦船振動(dòng)的直接激勵(lì)源，對(duì)其的精確預(yù)測是分析艦船振動(dòng)與噪聲性能的基礎(chǔ)與前提［1-3］。

自1978年Kerwin等［4］根據(jù)勢(shì)流理論提出用渦格升力面法計(jì)算螺旋槳的UBF以來，螺旋槳UBF數(shù)值計(jì)算方法得到了快速的發(fā)展，基于勢(shì)流理論的升力面法［5-6］及面元法［7-8］逐漸走向成熟，但畢竟勢(shì)流理論方法沒有考慮流體的粘性作用，與實(shí)際的流場特征還存在一定的偏差，而粘流技術(shù)正好彌補(bǔ)了勢(shì)流理論的這個(gè)缺點(diǎn)。目前，螺旋槳UBF的粘流數(shù)值計(jì)算主要集中在采用滑移網(wǎng)格（SM）技術(shù)研究不同的流場條件、與槳葉相關(guān)的基本參數(shù)等［9-13］對(duì)螺旋槳UBF的影響。然而，螺旋槳軸承力的非定常變化特性非常復(fù)雜，其脈動(dòng)幅值占平均量的比例相對(duì)很小，且螺旋槳隨艦船和軸系振動(dòng)而誘發(fā)的流體脈動(dòng)效應(yīng)會(huì)導(dǎo)致螺旋槳的軸承力非定常變化特性發(fā)生急劇變化，尤其是在艦船和軸系處于主共振狀態(tài)時(shí)，這種固體與流體的耦合效應(yīng)會(huì)變得更為顯著。

對(duì)于流體與固體之間的相互干涉作用，雖然流固耦合方法能夠?qū)ζ溥M(jìn)行計(jì)算，但需要的時(shí)間太長、計(jì)算難度太大，而SM方法通過整個(gè)計(jì)算域的運(yùn)動(dòng)來模擬流場的變化也很難捕捉到螺旋槳的這種流固干涉效應(yīng)。動(dòng)網(wǎng)格（DM）方法則是通過流場中網(wǎng)格的變化來模擬固體物理模型的運(yùn)動(dòng)與變形，從而可簡易地實(shí)現(xiàn)流固干涉耦合效應(yīng)，計(jì)算難度相對(duì)較低，且DM方法在離心泵、渦輪壓縮機(jī)等旋轉(zhuǎn)機(jī)械［14-15］的非定常流固耦合及艦船螺旋槳［16-17］的敞水和空泡性能的數(shù)值計(jì)算中已得到了應(yīng)用。

為了探究DM方法在螺旋槳UBF數(shù)值計(jì)算中的適用性，本文根據(jù)文獻(xiàn)［18］中螺旋槳的UBF及非均勻流場的試驗(yàn)測量數(shù)據(jù)，將分別采用DM和SM方法并結(jié)合RANS湍流模型模擬螺旋槳的非定常流場變化，經(jīng)過與SM方法、勢(shì)流（PM）方法及試驗(yàn)等結(jié)果的對(duì)比分析來驗(yàn)證DM在螺旋槳UBF數(shù)值計(jì)算中的適用性，為后續(xù)開展外部激勵(lì)作用對(duì)螺旋槳UBF影響的數(shù)值研究打下基礎(chǔ)。

1　數(shù)值計(jì)算方法

1.1控制方程

本文使用Fluent流體計(jì)算軟件，分別采用DM 和SM方法對(duì)螺旋槳的三維不可壓縮非定常流場進(jìn)行數(shù)值模擬，其流場數(shù)值計(jì)算的控制方程表達(dá)式為［19］

式中：u=（ux，uy，uz），為流場的瞬時(shí)速度矢量；ug=（ug，x，ug，y，ug，z），為網(wǎng)格節(jié)點(diǎn)的運(yùn)動(dòng)速度矢量；，為應(yīng)力張量；μ為流體的動(dòng)力學(xué)粘性系數(shù)；ρ為流體的密度；p為壓力；S為控制體積V的面積。

1.2DM方法

在用DM方法模擬螺旋槳流場的過程中，設(shè)數(shù)值計(jì)算的時(shí)間步長為Δt，當(dāng)螺旋槳開始工作時(shí)，螺旋槳將以給定的轉(zhuǎn)速作繞軸旋轉(zhuǎn)的剛體運(yùn)動(dòng)，螺旋槳周圍的網(wǎng)格節(jié)點(diǎn)將隨著螺旋槳的運(yùn)動(dòng)而進(jìn)行自我調(diào)節(jié)和更新，因此，在t=nΔt，…，（n+1）Δt的計(jì)算時(shí)間內(nèi)，網(wǎng)格節(jié)點(diǎn)所在控制體單元的體積將發(fā)生變化。

首先，在滿足網(wǎng)格守恒規(guī)律的基礎(chǔ)上，通過求解控制體各個(gè)面上的點(diǎn)積ug，j·Sj來求取控制體的體積時(shí)間導(dǎo)數(shù)dV/dt，即

總之，在高中語文教學(xué)過程中關(guān)注生命體驗(yàn)，不僅體現(xiàn)了語文學(xué)科的人文性特點(diǎn)，而且有助于學(xué)生形成正確的生命價(jià)值觀。語文教師應(yīng)將生命體驗(yàn)不斷地滲透教學(xué)過程中，培養(yǎng)學(xué)生的生命情懷，這樣既可以讓他們學(xué)到知識(shí)，實(shí)現(xiàn)既定的教學(xué)目標(biāo)，也能讓他們學(xué)會(huì)認(rèn)識(shí)生命的價(jià)值，獲得生命的成長，不斷向著理想的人生進(jìn)發(fā)。

式中：nf為控制體積上的面的數(shù)目；Sj為第j面積向量；δVj為控制體上第j個(gè)方面在時(shí)間步長Δt內(nèi)掃描出來的體積。

然后，通過第nΔt時(shí)間步長的控制體積Vn以及控制體的體積時(shí)間導(dǎo)數(shù)dV/dt來求得第（n+1）Δt時(shí)間步長的控制體積Vn+1，即

之后，采用一階向后差分的離散方法對(duì)式（2）中的時(shí)間項(xiàng)進(jìn)行離散，其表達(dá)式為

在Fluent軟件中，采用UDF自定義程序?qū)崿F(xiàn)螺旋槳的剛體運(yùn)動(dòng)，并通過DM算法對(duì)螺旋槳周圍的網(wǎng)格節(jié)點(diǎn)進(jìn)行調(diào)節(jié)和更新。DM算法可分為彈性光順法、局部重構(gòu)法和鋪層法3種，不同的DM算法對(duì)流場的網(wǎng)格類型要求不同，其中彈性光順法和局部重構(gòu)法主要用于四面體網(wǎng)格，鋪層法用于六面體網(wǎng)格或三菱柱體網(wǎng)格。由于螺旋槳葉片存在較大的曲率變化，在其周圍產(chǎn)生規(guī)則的六面體或三菱柱體網(wǎng)格難度較大，故通常采用四面體網(wǎng)格對(duì)流場計(jì)算域空間進(jìn)行離散。針對(duì)以上特點(diǎn)，本文在采用DM方法模擬螺旋槳流場的數(shù)值過程中，采用彈性光順與局部重構(gòu)相結(jié)合的方法對(duì)螺旋槳流場中的體網(wǎng)格節(jié)點(diǎn)進(jìn)行調(diào)節(jié)和控制。

1.3SM方法

根據(jù)SM方法的基本原理，通常將流場計(jì)算域劃分為運(yùn)動(dòng)和靜止2部分，這2部分計(jì)算域之間通過交界面實(shí)現(xiàn)網(wǎng)格的相互滑移并進(jìn)行數(shù)據(jù)的交換。與DM方法有所不同，在采用SM方法模擬螺旋槳流場的過程中，流場計(jì)算域內(nèi)部的網(wǎng)格控制體單元的體積不發(fā)生變化，即

在螺旋槳的流場數(shù)值計(jì)算過程中采用有限體積法離散連續(xù)性方程和動(dòng)量方程，其中對(duì)流項(xiàng)采用二階迎風(fēng)格式，擴(kuò)散項(xiàng)采用二階中心差分格式，壓強(qiáng)與速度的耦合迭代過程選用SIMPLE算法，對(duì)于湍流的模擬采用SST k-ω湍流模型［20］對(duì)控制方程進(jìn)行封閉。

2　計(jì)算模型

2.1物理模型與網(wǎng)格

圖1　螺旋槳流場計(jì)算域示意圖及相關(guān)參數(shù)Fig.1　Schematic and parameters of computational domain of the propeller

對(duì)螺旋槳流場計(jì)算域進(jìn)行空間離散，槳葉及槳轂表面鋪設(shè)三角形網(wǎng)格，并在其周圍布置多層棱柱形邊界層體網(wǎng)格。在槳葉的導(dǎo)邊、隨邊及葉梢處進(jìn)行網(wǎng)格的局部加密，內(nèi)圓柱體計(jì)算域采用四面體網(wǎng)格進(jìn)行填充，外圓柱體計(jì)算域采用全六面體網(wǎng)格。流場數(shù)值計(jì)算中螺旋槳的幾何模型見圖2（a），槳葉表面網(wǎng)格分布見圖2（b）。

圖2　DTMB 4119模型及槳葉表面網(wǎng)格分布Fig.2　DTMB 4119 model and grids on the blade surface

2.2邊界條件

為了與試驗(yàn)結(jié)果進(jìn)行對(duì)比分析，本文選取文獻(xiàn)［18］中的試驗(yàn)工況進(jìn)行數(shù)值計(jì)算，其中，螺旋槳的進(jìn)速系數(shù)J=0.833，轉(zhuǎn)速n=600 r/min，進(jìn)口來流采用UBF自定義程序給定九階角頻的非均勻流場U，周向遠(yuǎn)場邊界給定速度進(jìn)口，出口設(shè)定為壓力出口。

根據(jù)文獻(xiàn)［18］中測量的流場數(shù)據(jù)，可得到如圖3所示螺旋槳流場進(jìn)口處各半徑上的九階非均勻流場的周向分布圖，其中θ為周向角度。從圖中可以看出，在0.3R處，進(jìn)口流場周向的速度分布較為均勻，九階峰值不是很明顯，且最大峰值僅為最大谷值的1.05倍。隨著半徑的增大，流場的周向九階非均勻特性變得越來越明顯，在0.5R和0.7R處，進(jìn)口流速周向分布的九階來流速度的最大峰值與最大谷值的比值分別為1.22和1.37。當(dāng)進(jìn)口流場半徑為0.9R時(shí)，進(jìn)口流速的周向分布具有十分明顯的9個(gè)峰值，且九階角頻來流速度的最大峰值已為最大谷值的1.55倍。

圖3　進(jìn)口各半徑處九階非均勻流場的周向分布Fig.3　Circumferential distribution of nine-order non-uniform fields at different radius positions

3　結(jié)果分析

螺旋槳UBF由6個(gè)分量組成，其中軸向的推力Tx和扭矩Qx為主要成分。通過螺旋槳的流場模擬計(jì)算得到螺旋槳工作時(shí)的UBF，并定義螺旋槳UBF的推力系數(shù)KT和扭矩系數(shù)KQ分別為：

式中：ρ為流體密度，kg/m3；n為螺旋槳轉(zhuǎn)速，s-1；D為螺旋槳直徑，m。

3.1網(wǎng)格無關(guān)性分析

數(shù)值計(jì)算研究中，網(wǎng)格質(zhì)量是一個(gè)對(duì)計(jì)算結(jié)果有著顯著影響的因素，故對(duì)網(wǎng)格無關(guān)性進(jìn)行驗(yàn)證是數(shù)值計(jì)算研究中必要的工作。在此，首先使用3套不同質(zhì)量的計(jì)算網(wǎng)格，使用DM方法進(jìn)行螺旋槳流場的數(shù)值模擬，對(duì)網(wǎng)格無關(guān)性問題進(jìn)行驗(yàn)證，網(wǎng)格總數(shù)如表1所示。

表1　3套網(wǎng)格的網(wǎng)格總數(shù)Tab.1　Total number of three sets of grid

圖4給出了螺旋槳旋轉(zhuǎn)1周的時(shí)間間隔內(nèi)3套網(wǎng)格計(jì)算得到的螺旋槳UBF時(shí)域特性。從圖中可以看出，3套網(wǎng)格計(jì)算得到的螺旋槳UBF時(shí)域曲線的變化規(guī)律基本一致，三者之間吻合很好，均體現(xiàn)了明顯的九階峰值變化特性，但由3套網(wǎng)格計(jì)算得到的螺旋槳UBF時(shí)域曲線之間也存在著相對(duì)較小的差異，即隨著網(wǎng)格數(shù)的增多，螺旋槳UBF時(shí)域曲線的變化幅度稍有增大。

圖4　3套網(wǎng)格下的螺旋槳UBF時(shí)域特性Fig.4　Time series features of a propeller UBFs in three sets of grid

表2給出了3套網(wǎng)格計(jì)算得到的螺旋槳UBF的計(jì)算精度。從表中可以看出，在3套網(wǎng)格計(jì)算得到的螺旋槳UBF的各階諧頻中，三階諧頻（N=3）的計(jì)算精度相對(duì)較低，九階諧頻（N=9）的計(jì)算精度最高。對(duì)比分析3套網(wǎng)格的計(jì)算精度發(fā)現(xiàn)，網(wǎng)格2和網(wǎng)格3要明顯優(yōu)于網(wǎng)格1，且網(wǎng)格2與網(wǎng)格1的最大精度差為5.87%，網(wǎng)格3與網(wǎng)格1的最大精度差為6.12%，而網(wǎng)格2與網(wǎng)格3的計(jì)算精度則較為接近，兩者之間的最大精度差僅為1.36%。因此，從3套網(wǎng)格的數(shù)值計(jì)算結(jié)果來看，網(wǎng)格2已能夠滿足螺旋槳UBF數(shù)值計(jì)算的網(wǎng)格無關(guān)性要求，從提高計(jì)算效率的角度來考慮，網(wǎng)格2最適用于本文的研究工作。

3.2DM方法適用性分析

選擇網(wǎng)格2，根據(jù)文獻(xiàn)［18］中的試驗(yàn)工況分別采用DM和SM方法進(jìn)行螺旋槳的UBF數(shù)值計(jì)算，并與文獻(xiàn)中的PM方法計(jì)算結(jié)果及試驗(yàn)結(jié)果進(jìn)行對(duì)比分析。

圖5　采用DM和SM方法計(jì)算得到的螺旋槳UBF時(shí)域特性Fig.5　Time series features of the UBFs of a propeller obtained from DM and SM method

表2　3套網(wǎng)格計(jì)算的螺旋槳UBF各階諧頻下的相對(duì)誤差Tab.2　Relative error of each order harmonic of UBFs in three sets of grid

圖5給出了采用DM和SM方法計(jì)算得到的螺旋槳UBF的時(shí)域特性，即螺旋槳旋轉(zhuǎn)1周的時(shí)間間隔內(nèi)的UBF非定常變化特性。從圖中可以看出，由兩種粘流方法（DM和SM）計(jì)算得到的螺旋槳UBF時(shí)域曲線的整體變化趨勢(shì)吻合很好，即均出現(xiàn)了較為明顯的9個(gè)峰值，且9個(gè)峰值均可分為3組，這體現(xiàn)了非均勻來流的九階角頻特性和三葉常規(guī)螺旋槳的一階葉頻特性，計(jì)算結(jié)果與文獻(xiàn)［10］的計(jì)算結(jié)果也比較吻合。但兩種方法的UBF時(shí)域曲線也存在著細(xì)微的差異，即SM方法計(jì)算得到的時(shí)域曲線相對(duì)于DM方法較為光滑。從兩種粘流方法的時(shí)域曲線可以說明，就螺旋槳UBF數(shù)值計(jì)算而言，DM方法能夠達(dá)到與SM方法相近的計(jì)算精度。

圖6　采用DM和SM方法計(jì)算得到的3個(gè)不同盤面上的壓力分布Fig.6　Pressure distribution calculated by DM and SM method at three different discs

圖6給出了當(dāng)θ=0°時(shí)，由DM和SM方法計(jì)算得到的螺旋槳流場在3個(gè)不同盤面上的壓力分布。從3個(gè)不同盤面上的壓力分布的對(duì)比分析可看出，由DM和SM方法計(jì)算得到的壓力分布整體趨勢(shì)比較相似，即在x=-0.05 m盤面上，對(duì)應(yīng)于槳葉前方的位置為較大面積的負(fù)壓區(qū)，而對(duì)應(yīng)于槳葉流道的位置已逐漸轉(zhuǎn)變?yōu)檎龎簠^(qū)。在x=0 m盤面上，螺旋槳槳葉吸力面附近存在較大面積的負(fù)壓區(qū)，且隨著遠(yuǎn)離吸力面的距離的增大，壓力逐漸增大，吸力面與壓力面之間的流道內(nèi)存在著較大面積的正壓區(qū)，在壓力面附近，槳轂和葉梢附近的壓力相對(duì)較低。在x=0.05 m盤面上，螺旋槳的槳轂周圍存在相對(duì)較低的負(fù)壓區(qū)，而對(duì)應(yīng)于槳葉流道的位置存在較為明顯的正壓區(qū)。然而，由兩種粘流方法計(jì)算得到的流場之間也存在著較小的差別，即在流場計(jì)算域內(nèi)的交界面附近，相對(duì)于DM方法，SM方法計(jì)算得到的壓力分布較為零散，這可能是因?yàn)镾M方法是通過交界面上網(wǎng)格的相互滑移來實(shí)現(xiàn)動(dòng)靜計(jì)算域的數(shù)據(jù)交換所引起的。

圖7給出了當(dāng)θ=0°時(shí)，螺旋槳后x=0.05 m盤面上采用SM和DM方法計(jì)算得到的軸向速度矢量圖。從圖中可以看出，由SM和DM方法計(jì)算得到的軸向速度的大小及流動(dòng)方向均較為相似。但兩種方法計(jì)算得到的矢量圖中速度的節(jié)點(diǎn)分布稍有不同，這主要是由于DM方法需在數(shù)值計(jì)算過程中的每一時(shí)間步長進(jìn)行體網(wǎng)格節(jié)點(diǎn)的調(diào)節(jié)與更新，而SM方法則保持體內(nèi)網(wǎng)格節(jié)點(diǎn)不發(fā)生變化。

圖7　螺旋槳后x=0.05 m盤面上軸向速度矢量圖Fig7　Vector diagram of axial velocity on the disc x=0.05 m behind the propeller

表3給出了DM，SM和PM方法在螺旋槳UBF數(shù)值計(jì)算中的計(jì)算精度，其中PM方法的計(jì)算結(jié)果出自文獻(xiàn)［18］。從表中可以看出，對(duì)于螺旋槳的UBF數(shù)值計(jì)算，DM和SM方法的計(jì)算精度十分接近，除了在諧頻N=6時(shí)SM方法的計(jì)算精度稍優(yōu)于DM方法之外，其它諧頻下的DM方法的計(jì)算精度都稍高于SM方法。對(duì)比粘流方法與PM方法的計(jì)算結(jié)果，在諧頻N=0和N=3時(shí)，PM方法的計(jì)算精度要稍高于DM方法和SM方法，但在諧頻N=6 和N=9時(shí)，PM方法的最大誤差分別為69.77%和34.06%，其計(jì)算精度遠(yuǎn)遠(yuǎn)低于兩種粘流方法。通過上述對(duì)比分析，可以認(rèn)為DM方法在螺旋槳UBF數(shù)值計(jì)算中能夠滿足精度要求。

表3　DM，SM和PM方法計(jì)算得到的螺旋槳UBF各階諧頻下的相對(duì)誤差Tab.3　Relative error of every order harmonic of UBFs for DM，SM and PM

通過與SM，PM方法及試驗(yàn)等結(jié)果的對(duì)比分析，可以看出DM方法用于螺旋槳UBF的數(shù)值計(jì)算是可行的，計(jì)算精度能夠滿足要求，且DM方法通過網(wǎng)格的變化來實(shí)現(xiàn)物理模型的運(yùn)動(dòng)與變形，其計(jì)算過程更為接近物理模型的真實(shí)變化特性。

4　結(jié)　語

本文分別采用DM和SM方法數(shù)值模擬了非均勻來流條件下螺旋槳軸承力的非定常變化特性，并與PM方法及試驗(yàn)等結(jié)果進(jìn)行了對(duì)比分析。從計(jì)算結(jié)果來看，在螺旋槳UBF的數(shù)值計(jì)算中，DM和SM方法的計(jì)算結(jié)果與試驗(yàn)結(jié)果吻合很好，且兩者的計(jì)算精度均遠(yuǎn)遠(yuǎn)優(yōu)于PM方法。在螺旋槳流場的3個(gè)不同盤面上，由兩種粘流方法計(jì)算得到的壓力分布除了在交界面附近存在極小的差異之外，其它壓力分布趨勢(shì)基本一致，而螺旋槳后盤面上速度矢量的節(jié)點(diǎn)分布除因DM方法中體網(wǎng)格的調(diào)節(jié)與更新稍有不同之外，由兩種方法計(jì)算得到的速度值和流向均較為相似。上述分析結(jié)果說明，DM方法應(yīng)用于螺旋槳UBF數(shù)值計(jì)算是可行的。通過3套不同數(shù)量的網(wǎng)格對(duì)螺旋槳的流場進(jìn)行數(shù)值模擬，其結(jié)果均能夠滿足計(jì)算精度的要求，且網(wǎng)格2和網(wǎng)格3的最大精度差僅為1.36%，說明網(wǎng)格已滿足無關(guān)性要求。

［1］何友聲，王國強(qiáng).螺旋槳激振力［M］.上海：上海交通大學(xué)出版社，1987.

［2］ 鄭小龍，王超，張立新，等.非均勻流場中螺旋槳水動(dòng)力及噪聲特性預(yù)報(bào)研究［J］.武漢理工大學(xué)學(xué)報(bào)（交通科學(xué)與工程版），2014，38（6）：1300-1307.

ZHENG Xiaolong，WANG Chao，ZHANG Lixin，et al.Prediction research on propeller hydrodynamic and noise characteristics in non-uniform flow［J］.Journal of Wuhan University of Technology（Transportation Sci?ence&Engineering Edition），2014，38（6）：1300-1307.

［3］ BOSWELL R J，MILLER M L.Unsteady propellerloading-measurement，correlation with theory，and parametric study［R］.Washington：Naval Ship Re?search and Development Center，1968.

［4］KERWIN J E，LEE C S.Prediction of steady and un?steady marine propeller performance by numerical lift?ing surface theory［J］.Transactions of the Society of Naval Architects and Marine Engineers，1978，86：218-253.

［5］陳家棟.非定常螺旋槳水動(dòng)力升力面預(yù)報(bào)［J］.水動(dòng)力學(xué)研究與進(jìn)展（A輯），1992，7（4）：420-433.

CHEN Jiadong.Unsteady hydrodynamics prediction of propellers with a lifting surface method［J］.Journal of Hydrodynamics（Ser.A），1992，7（4）：420-433.

［6］ 馬騁，錢正芳，陳科，等.采用渦格法和面元法預(yù)報(bào)吊艙推進(jìn)器非定常水動(dòng)力性能［J］.船舶力學(xué)，2014，18（9）：1035-1043.

MA Cheng，QIAN Zhengfang，CHEN Ke，et al.Using vortex lattice and surface panel method to predict the unsteady hydrodynamic performance of podded propul?sors［J］.Journal of Ship Mechanics，2014，18（9）：1035-1043.

［7］HOSHINO T.Hydrodynamic analysis of propellers in unsteady flow using a surface panel method［J］.Jour?nal of the Society of Naval Architects of Japan，1993，174：71-87.

［8］張赫，崔和.面元法預(yù)估導(dǎo)管螺旋槳水動(dòng)力性能［J］.水雷戰(zhàn)與艦船防護(hù)，2013，21（3）：33-37.

ZHANG He，CUI He.Predicting hydrodynamic perfor?mance of ducted propeller using surface panel method ［J］.Mine Warfare and Ship Self-Defence，2013，21 （3）：33-37.

［9］KOUSHAN K，KRASILNIKOV V.Experimental and numerical investigation of open thrusters in oblique flow conditions［C］//Proceeding of the 27th ONR Sym?posium on Naval Hydrodynamics.Seoul，Korea：［s.n.］，2008.

［10］胡小菲，黃振宇，洪方文.螺旋槳非定常力的黏性數(shù)值分析［J］.水動(dòng)力學(xué)研究與進(jìn)展（A輯），2009，24 （6）：734-739.

HU Xiaofei，HUANG Zhenyu，HONG Fangwen.Un?steady hydrodynamics forces of propeller predicted with viscous CFD［J］.Journal of Hydrodynamics （Ser.A），2009，24（6）：734-739.

［11］GUO C Y，MA N，YANG C J.Numerical simulation of a podded propulsor in viscous flow［J］.Journal of Hy?drodynamics（Ser.B），2009，21（1）：71-76.

［12］黃振宇.螺旋槳的非定常力計(jì)算［C］//第九屆全國水動(dòng)力學(xué)學(xué)術(shù)會(huì)議暨第二十二屆全國水動(dòng)力學(xué)研討會(huì)文集.北京：中國力學(xué)學(xué)會(huì)，2009.

HUANG Zhenyu.Numerical simulation of unsteady force of marine propeller［C］//Proceeding of the 9th National Congress on Hydrodynamics&22th National Conference on Hydrodynamics.Beijing：Chinese Soci?ety of Theoretical and Applied Mechanics，2009.

［13］張瑞，王先洲，張志國，等.側(cè)斜變化對(duì)螺旋槳水動(dòng)力及變形振動(dòng)特性的影響［J］.中國艦船研究，2015，10（6）：87-94.

ZHANG Rui，WANG Xianzhou，ZHANG Zhiguo，et al.Effects of blade skew on the hydrodynamic and de?formation performance of propellers［J］.Chinese Jour?nal of Ship Research，2015，10（6）：87-94.

［14］陳黎卿，張棟，陳無畏.基于流固耦合的分動(dòng)器齒輪兩相流動(dòng)數(shù)值模擬與試驗(yàn)［J］.農(nóng)業(yè)工程學(xué)報(bào)，2014，30（4）：54-61.

CHEN Liqing，ZHANG Dong，CHEN Wuwei.Nu?merical simulation and test on two-phase flow inside shell of transfer case based on fluid-structure interac?tion［J］.Transactions of the Chinese Society of Agri?cultural Engineering，2014，30（4）：54-61.

［15］江帆，陳維平，王一軍，等.基于動(dòng)網(wǎng)格的離心泵內(nèi)部流場數(shù)值模擬［J］.流體機(jī)械，2007，35（7）：20-24.

JIANG Fan，CHEN Weiping，WANG Yijun，et al. Numerical simulation of flow field inside of centrifu?gal pump based on dynamics mesh［J］.Fluid Machin?ery，2007，35（7）：20-24.

［16］ 王順杰，王易川.基于動(dòng)網(wǎng)格的螺旋槳空化數(shù)值模擬［J］.艦船科學(xué)技術(shù)，2013，35（3）：15-18.

WANG Shunjie，WANG Yichuan.Numerical simula?tion of propeller cavitation based on dynamic mesh ［J］.Ship Science and Technology，2013，35（3）：15-18.

［17］繆宇躍，孫江龍.CFD敞水螺旋槳性能計(jì)算分析［J］.中國艦船研究，2011，6（5）：63-68.

MIAO Yuyue，SUN Jianglong.CFD analysis of hydro?dynamic performance of propeller in open water［J］. Chinese Journal of Ship Research，2011，6（5）：63-68.

［18］JESSUP S D.Measurement of multiple blade rate un?steady propeller forces［R］.Bethesda，MD：David Taylor Naval Ship R&D Center，1990.

［19］ANSYS FLUENT theory guide［S］.［S.l.］：ANSYS Inc，2011.

［20］MENTER F R.Two-equation eddy-viscosity turbu?lence models for engineering applications［J］.AIAA Journal，1994，32（8）：1598-1605.

［21］JESSUP S D.An experimental investigation of viscous aspects of propeller blade flow［D］.Washington：The Catholic University of America，1989.

［22］BROCKETT T.Minimum pressure envelopes for mod?ified NACA-66 sections with NACA a=0.8 Camber and BUSHIPS type 1 and type 2 sections［R］.Wash?ington：David Taylor Model Basin，1966.

Application of the dynamic mesh method in the numerical simulation of unsteady bearing forces of a propeller

SHU Minhua1，2，CHEN Ke1，2，YOU Yunxiang1，2，HU Tianqun1，2，LIU Heng1，2
1 School of Naval Architecture，Ocean and Civil Engineering，Shanghai Jiao Tong University，

Shanghai 200240，China
2 State Key Laboratory of Ocean Engineering，Shanghai Jiao Tong University，Shanghai 200240，China

The applicability of the Dynamic Mesh（DM）method in numerical simulation of propeller Un?steady Bearing Forces（UBF）is studied in this paper.The grid independence is initially verified by compar?ing three sets of grid of different qualities and the most suitable set of grid is selected for the work.Further numerical simulation is conducted by adopting DM method，and the corresponding results are compared with those obtained from the Sliding Mesh（SM）method，Potential Flow（PM）method and experiments.The analysis indicates that for the various characteristics of a propeller UBF，the results of the DM method are consistent with those of the SM method and experiments，and are better than those of the PM method.Mean?while，the flow field distribution obtained from the DM method agrees well with that from the SM method. Therefore，the application of the DM method in the numerical simulation of a propeller UBFs is feasible and effective.

ship；three-blade conventional propeller；Unsteady Bearing Force（UBF）；Dynamic Mesh （DM）；Sliding Mesh（SM）；numerical calculation

U664.33

A

10.3969/j.issn.1673-3185.2016.03.005

2015-09-23網(wǎng)絡(luò)出版時(shí)間：2016-5-31 11:04

國家級(jí)重大基礎(chǔ)研究項(xiàng)目

舒敏驊，男，1984年生，博士生。研究方向：螺旋槳水動(dòng)力學(xué)及振動(dòng)噪聲。

E-mail：shuminhua@126.com

尤云祥（通信作者），男，1963年生，博士，教授，博士生導(dǎo)師。研究方向：內(nèi)波水動(dòng)力學(xué)及其工程后效。E-mail：youyx@sjtu.edu.cn