,

(大連理工大學(xué)a.工業(yè)裝備結(jié)構(gòu)分析國家重點實驗室；b.船舶工程學(xué)院，遼寧 大連 116024)

近20年通過CFD技術(shù)預(yù)測螺旋槳水動力性能方面的應(yīng)用越來越多，計算精度越來越高。目前，采用RANS湍流模型預(yù)報螺旋槳的水動力性能，結(jié)果比較準(zhǔn)確，這為螺旋槳的水動力性能研究提供了新的途徑[1-7]。由于螺旋槳在水中做旋轉(zhuǎn)前進，采用多參考坐標(biāo)系模型不能模擬螺旋槳的真實運動，這種模型也只能用于螺旋槳定常水動力性能研究。滑移網(wǎng)格技術(shù)不僅能夠模擬螺旋槳實際旋轉(zhuǎn)，而且還能用于螺旋槳非定常水動力性能研究。本文通過Gambit軟件建立螺旋槳的三維模型，應(yīng)用Fluent軟件，選取3種不同的湍流模型預(yù)報螺旋槳的定常水動力性能，得出在不同進速系數(shù)時螺旋槳的敞水性能曲線。采用SM模型預(yù)報螺旋槳的非定常水動力性能，并對MRF和SM模型計算結(jié)果進行比效。

1 數(shù)學(xué)模型

1.1 控制方程

假定流體是不可壓縮的，則流場的連續(xù)方程和動量方程[8]分別為

(1)

(2)

式中：ui，uj——速度分量時均值(i，j=1,2,3)；

p——壓力時均值；

ρ——流體密度；

μ——流體粘性系數(shù)；

gi——重力加速度分量；

方程中的雷諾應(yīng)力項屬于新的未知量，因此，要使方程封閉，必須對該應(yīng)力項做某種假設(shè)，即建立應(yīng)力的表達式(或引進新的湍流模型方程)，通過表達式或湍流模型，把應(yīng)力項中的脈動值與時均值聯(lián)系起來。

1.2 湍流模型

1.2.1 標(biāo)準(zhǔn)k-ε模型

在湍動能k方程的基礎(chǔ)上，引入一個耗散率ε的方程[9]，形成了k-ε雙方程模型，稱為標(biāo)準(zhǔn)的k-ε模型。標(biāo)準(zhǔn)k-ε模型的輸運方程為

Gk+Gb-ρε-YM+Sk

(3)

(4)

式中:Gk——由于平均速度梯度引起的湍動能產(chǎn)生;

Gb——由于浮力影響引起的湍動能的產(chǎn)生;

Ym——可壓縮湍流脈動膨脹對總的耗散率的影響;

C1ε,C2ε,C3ε——經(jīng)驗常數(shù),Fluent中默認值為C1ε=1.44,C2ε=1.92,C3ε=0.09;

δk,δε——湍動能和湍動耗散率對應(yīng)的普朗特數(shù)，F(xiàn)luent中默認值為δk=1.0,δε=1.3;

μt——湍動粘度。

1.2.2 RNGk-ε模型

RNGk-ε模型是對瞬時N-S方程用重整化的數(shù)學(xué)方法推導(dǎo)出來的模型。模型中的常數(shù)與標(biāo)準(zhǔn)的模型不同,而且方程中也出現(xiàn)了新的函數(shù)或項。所得到的方程如下。

Gk+Gb-ρε-YM

(5)

(6)

式中:μeff=μ+μt,μt=ρCμ(k2/ε);

αk、αε——湍動能和湍動耗散率的有效的普朗特數(shù)的倒數(shù);

對于高雷諾數(shù)問題,經(jīng)驗常數(shù)Cμ=0.084 5。

1.2.3 Reliablek-ε模型

對于時均應(yīng)變力特別大的情形，標(biāo)準(zhǔn)的k-ε模型有可能導(dǎo)致負的正應(yīng)力。為使流動符合湍流物理定律，需要對正應(yīng)力進行某種數(shù)學(xué)約束。Reliablek-ε模型正是基于此原理提出的，其運輸方程為

Gk+Gb-ρε-YM

(7)

(8)

σk、σε——湍動能和湍動耗散率對應(yīng)的普朗特數(shù)，F(xiàn)luent中默認值為σk=1.0，σε=1.3;

μt——湍動粘度。

2 計算模型

2.1 幾何建模與計算域的創(chuàng)建

本文采用DTMB4119標(biāo)準(zhǔn)螺旋槳，其尺寸如下[10]。

幾何參數(shù)DTMB4119直徑/m0.304 8葉數(shù)3螺距比(0.7R)1.084轂徑比0.2縱傾角/(°)0側(cè)斜角/(°)0葉剖面NACA66mod

幾何建模中采用了CFD前處理器Gambit軟件，將螺旋槳表面數(shù)據(jù)點導(dǎo)入Gambit中，采用樣條曲線擬合各個型值點，從而建立光滑的三維螺旋槳表面外形，見圖1。

圖1 螺旋槳幾何模型

由于研究的是螺旋槳的水動力性能，因此需要把模型置于流場域中，計算域由兩部分組成，外部大圓柱和內(nèi)部小圓柱,計算域內(nèi)邊界取在槳榖和葉片表面上，其中槳榖被簡化為中間為圓柱面，兩端為橢球面；外邊界面取在大圓柱體表面。

2.2 網(wǎng)格劃分

采用分區(qū)網(wǎng)格劃分方案，螺旋槳表面采用三角形網(wǎng)格進行劃分，見圖2，外部大圓柱形區(qū)域采用結(jié)構(gòu)化六面體網(wǎng)格，在接近小圓柱區(qū)域采用局部加密，包含螺旋槳面的內(nèi)部小圓柱區(qū)域采用非結(jié)構(gòu)化四面體網(wǎng)格，內(nèi)外區(qū)之間通過定義interface來傳替，見圖3。

圖2 螺旋槳表面網(wǎng)格

圖3 流場體網(wǎng)格截面

2.3 邊界條件的設(shè)置

進口邊界條件處設(shè)置為速度進口條件，命名為in，在不考慮空泡條件下速度入口給定來流速度；出口邊界定義為壓力出口邊界，命名為out，設(shè)置為0;小圓柱區(qū)域定義為旋轉(zhuǎn)區(qū)域命名為fluid2,設(shè)置轉(zhuǎn)速為600 r/min；大圓柱區(qū)域定義為靜止區(qū)域命名為fluid1；fluid1和fluid2的交界面設(shè)置為interface面。螺旋槳面設(shè)置為無滑移的固壁條件，在近壁區(qū)采用標(biāo)準(zhǔn)壁面函數(shù)。

2.4 螺旋槳模擬的旋轉(zhuǎn)模型

在計算區(qū)域中存在運動區(qū)域和靜止區(qū)域，F(xiàn)luent可以采用多參考系模型(MRF)，和滑移網(wǎng)格模型(SM)。MRF模型用于定常流動的計算，是最簡單、最經(jīng)濟的模型。在使用MRF模型進行計算時，整個計算域被分成多個小的子域。SM模型用于處理非定常問題，這是與MRF模型最大的區(qū)別，在計算時，SM模型要求計算域至少包含兩個以上存在相對運動的子域。SM模型相比較于MRF模型需要使用的計算資源比較大，對內(nèi)存，CPU速度都有較高的要求[11]。

3 數(shù)值模擬

3.1 螺旋槳定常水動力性能的預(yù)報

本文計算了進速系數(shù)為0.500、0.700、0.833、0.900時，分別采用的3種湍流模型即，標(biāo)準(zhǔn)k-ε模型、RNGk-ε模型和Reliablek-ε模型，數(shù)值計算了螺旋槳的推力和轉(zhuǎn)矩，數(shù)值結(jié)果見表1。

表1 推力T和轉(zhuǎn)矩Q計算結(jié)果

表1數(shù)據(jù)表明,采用Reliablek-ε湍流模型模擬出來的數(shù)值結(jié)果最好，RNGk-ε湍流模型次之，標(biāo)準(zhǔn)的k-ε湍流模型最差。將采用Reliablek-ε湍流模型計算的結(jié)果見圖4。

圖4 Reliable k-ε湍流模型計算結(jié)果

由圖4可知，模擬計算所得到的推力系數(shù)和轉(zhuǎn)矩系數(shù)與試驗偏差分別為3.4%和7.2%，推力系數(shù)模擬得較好，轉(zhuǎn)矩系數(shù)模擬的誤差偏大，但整體來說基本吻合。在J=0.833時吻合得最理想。

3.2 螺旋槳的非定常水動力性能預(yù)報

以采用MRF模型，Reliablek-ε湍流模型模擬的結(jié)果作為初始場，然后采用滑移網(wǎng)格,單機并行計算非定常螺旋槳的水動力性能。分別模擬了4種工況下螺旋槳的推力系數(shù)和轉(zhuǎn)矩系數(shù)，并與試驗值進行比較。計算結(jié)果見表2，繪制螺旋槳的特性曲線見圖5。

表2 推力系數(shù)kt和轉(zhuǎn)矩系數(shù)10 kq計算結(jié)果

圖5 敞水性能試驗值與計算值對比

由圖5可見不同進速系數(shù)下螺旋槳的水動力性能與試驗值的差距。由計算結(jié)果與試驗值可以求出kt、kq的平均誤差分別為1.33%和4.89%。數(shù)值結(jié)果與試驗值吻合較好，在J=0.833處吻合得最理想。

3.3 螺旋槳在不同的旋轉(zhuǎn)模型下的水動力性能分析

將由MRF模型和SM模型計算出來的螺旋槳的推力系數(shù)、扭轉(zhuǎn)系數(shù)和試驗值一起繪制成曲線，見圖6。

圖6 兩種模型敞水性能曲線

圖6表明，采用滑移網(wǎng)格技術(shù)計算出來的螺旋槳的水動力曲線結(jié)果更接近試驗值。對比這兩種方法可以看出，MRF算法求解出來的流場是一個充分發(fā)展的流場，這個流場再以一定的速度轉(zhuǎn)動就可以得到實際的流場，換句話說MRF實現(xiàn)的螺旋槳的旋轉(zhuǎn)并不是螺旋槳真正的旋轉(zhuǎn)，此時所求解的僅僅是周期性轉(zhuǎn)動某一瞬時的情況，這樣就把一個非定長的問題轉(zhuǎn)化成一個定長問題來求解；而SM方法計算的出來的流場就是實際的流場，也就是螺旋槳真正的轉(zhuǎn)動。采用SM模擬螺旋槳水動力性能耗用時間是MRF方法的4～5倍，并且對電腦的配置要求較高。

3.4 螺旋槳槳葉表面壓力分布

通過上述分析，可以看出采用滑移網(wǎng)格技術(shù)模擬的螺旋槳的水動力性能與試驗值更加吻合，尤其在進速系數(shù)J=0.833時吻合得最為理想。接下來將以進速系數(shù)J=0.833時的工況為例，對螺旋槳槳葉表面的壓力進行分析。

圖7和圖8顯示了槳葉吸力面和壓力面的壓力云圖。

圖7 吸力面壓力云圖

圖8 壓力面壓力云圖

從圖中可以看出，在吸力面壓力分布的主要特點是中間相當(dāng)一部分壓力較小，導(dǎo)邊與隨邊處相對較大。在壓力面上，從葉根到葉梢，壓力先增大然后逐漸減小，在葉梢處達到最小值，壓力由隨邊向?qū)н呏饾u增大，在導(dǎo)邊達到最大。

4 結(jié)論

在網(wǎng)格的劃分上采用了結(jié)構(gòu)網(wǎng)格與非結(jié)構(gòu)網(wǎng)格相結(jié)合的混合方法，這種方法更加符合像螺旋槳這種外形復(fù)雜的模型。不僅減少了網(wǎng)格數(shù)量，還可以獲得比較理想的結(jié)果。

將RANS方程與三種不同的湍流模型相結(jié)合，并采用MRF模型對螺旋槳定常水動力性能進行數(shù)值求解，通過計算結(jié)果與試驗值比較，驗證了本文數(shù)值方法的正確性。同時也說明了Reliablek-ε湍流模型更適于螺旋槳定常水動力性能的模擬，并具有較高的計算精度。并將此作為初始場，采用SM模型、單機并行計算螺旋槳非定常水動力性能。通過與定常計算結(jié)果比較，表明此方法得到的數(shù)值結(jié)果更加精確，論證了此方法的正確性。同時說明了采用SM模型得到的數(shù)值結(jié)果更加符合螺旋槳的實際情況。

運用CFD方法對螺旋槳J=0.833時槳葉表面壓力分布進行了分析。定性和定量上預(yù)報螺旋槳流場性能為螺旋槳設(shè)計、噪聲、振動的研究提供了理論依據(jù)和支持。

本文通過采用CFD軟件對螺旋槳流場特性進行分析。對計算結(jié)果分析比較淺顯，對模型的建立，網(wǎng)格劃分以及計算參數(shù)的設(shè)置還有值得改進的地方，希望在以后的工作中不斷地完善和提高。

[1] 劉 丹，陳鳳馨.CFD在計算船舶螺旋槳水動力性能中的應(yīng)用研究[J].現(xiàn)代制造工程，2010(4):18-20.

[2] 馮雪梅,陳鳳明，蔡榮泉.使用fluent軟件的螺旋槳敞水性能計算和考察[J].船舶,2006(1):14-19.

[3] 王 超,黃 勝,解學(xué)參.基于CFD方法的螺旋槳水動力性能預(yù)報[J].海軍工程大學(xué)學(xué)報,2008,20(4):107-112.

[4] 王 偉.基于滑移網(wǎng)格與RNG湍流模型的螺旋槳水動力性能研究[J].中國海洋平臺,2008,23(6):8-12.

[5] 高富東，潘存云，蔡汶珊，等.基于CFD的螺旋槳敞水性能數(shù)值分析與驗證[J].機械工程學(xué)報,2010,46(8):133-139.

[6] 顧鋮璋，鄭百林.船用螺旋槳敞水性能與槳葉應(yīng)力的數(shù)值分析[J].力學(xué)季刊，2011,32(3):440-443.

[7] 常欣，鄒經(jīng)湘，郭春雨，等.螺距角和縱傾角對全方向推進器水動力性能的影響[J].船海工程，2010，39(2)：26-29.

[8] 王福軍.計算流體動力學(xué)分析-CFD軟件原理與應(yīng)用[M].北京：清華大學(xué)出版社，2010.

[9] 朱紅鈞，林元華，謝龍漢.Fluent 12流體分析及工程仿真[M].北京:清華大學(xué)出版社，2011.

[10] 黃 勝,王 超,王詩洋.不同湍流模型在螺旋槳水動力性能計算中的應(yīng)用與比較[J].哈爾濱工程大學(xué)學(xué)報，2011,30(5):481-485.

[11] 胡 健，馬 騁，黃 勝.導(dǎo)管對螺旋槳水動力性能的影響[J].武漢理工大學(xué)學(xué)報：交通科學(xué)與工程版，2009，33(5)：992-995.