某雙體客船阻力特性數(shù)值仿真分析

石曉川,李春陽,孫光曉,仉永超,張祥,王玉帥

1.濰柴動力股份有限公司,山東 濰坊 261061;2.濰柴重機股份有限公司,山東 濰坊,261108

0 引言

雙體客船具有較好的穩(wěn)定性和操作靈活性,在對抗惡劣海況時表現(xiàn)優(yōu)異,得到廣泛應(yīng)用。船舶阻力是衡量船舶性能的重要指標之一,阻力不僅影響船舶的速度,同時也與船舶的使用性和經(jīng)濟性相關(guān)[1]。船舶在定速巡航工況下,主機的輸出功率主要用于抵抗航行中產(chǎn)生的阻力,因此,船舶阻力的預(yù)測有利于船主機的選型匹配,如何更加快速準確地預(yù)測船體阻力成為船舶設(shè)計過程中的關(guān)鍵。

船舶阻力預(yù)測的研究方法主要有理論研究、試驗研究和數(shù)值模擬。近年來,隨著計算機算力的提升,計算流體動力學(xué)(computational fluid dynamics,CFD)成為預(yù)測阻力的熱門手段[2-3]。趙丙乾等[4]使用CFD軟件計算了雙艉客船在靜水中的航行阻力,發(fā)現(xiàn)使用軟件對船體阻力特性的預(yù)測具有快速性和高精度的優(yōu)點,在計算范圍內(nèi)船舶總阻力的仿真與試驗結(jié)果的相對誤差在3%以內(nèi);Wang等[5]結(jié)合試驗研究和數(shù)值模擬的方法分析了高速平面雙體船在不同工況下的阻力,發(fā)現(xiàn)在高速工況下SSTk-ω湍流模型更接近試驗結(jié)果;文獻[6-7]對不同工況下的船舶阻力使用SSTk-ω湍流模型進行仿真,所得到的結(jié)果具備參考價值。仿真模型的不斷優(yōu)化使得計算結(jié)果的精度得到了進一步提高。

RANS法將流體的控制方程時均化處理,解決了流體在湍流工況下在時間與空間上的隨機性。文獻[8-10]通過RANS法結(jié)合數(shù)值模擬方法對三維船體阻力進行預(yù)測,發(fā)現(xiàn)仿真獲得的自由面興波趨勢與試驗得到的結(jié)果趨勢相同,具有較好的實際應(yīng)用價值。

本文中使用STAR-CCM+軟件,基于RANS法預(yù)測雙體船阻力,采用二因次換算法將船模阻力換算至實船阻力,并根據(jù)阻力預(yù)測結(jié)果對船主機進行選型匹配。

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

1.1 控制方程

將海水視為不可壓縮黏性流體,其連續(xù)性方程為:

(1)

式中: ▽為哈密頓算子,ux、uy、uz分別為x、y、z方向上的流體的平均速度,m/s。

流體的動量守恒方程為:

(2)

式中:ρ為流體密度,kg/m3;t為時間,s;p為微元流體壓力,Pa;ν為流體運動黏度,m2/s;fi為微元流體單位質(zhì)量力,m/s2;下標i、j分別代表笛卡爾坐標系橫坐標與縱坐標的方向。

對于任何形式的流體流動,均符合Navier-Stokes方程(N-S方程),但隨著雷諾數(shù)增大,流動由層流轉(zhuǎn)變?yōu)橥牧?過度流動使得N-S方程的初值與邊值問題產(chǎn)生分歧解。為解決該問題,目前常采用時間平均法對方程進行處理,即:

(3)

將式(3)分別帶入到式(1)(2)中可以得到RANS方程[11]

(4)

(5)

1.2 自由表面數(shù)值模型

船體在水面上航行,需要考慮氣-液兩相面對船舶前進阻力的影響[12]。自由面捕捉算法是將液相和氣相同時計算,具有較好適用性和靈活性,對于破碎波等自由液面具有更好的模擬效果。

單一流體的物性參數(shù)由構(gòu)成函數(shù)ci決定:

(6)

式中:V為微元體積,m3;A為微元表面積,m2;ud為微元表面法向矢量上的速度,m/s;n為法向矢量。計算過程中,氣相中ci=0,液相中ci=1。

2 物理模型

某雙體客船的長度為33.4 m,垂線間長30.9 m。將雙體船的實船尺寸與模型尺寸按10：1縮小進行靜水場航行阻力計算,實船與模型的主要參數(shù)如表1所示。

表1 雙體船主要參數(shù)

根據(jù)重力相似準則,實船與模型的2個流動系統(tǒng)重力必然相似,其弗勞德數(shù)相等。弗勞德數(shù)[13]

(7)

式中:v為特征速度,m/s;g為自由落體加速度,m/s2;L為船舶特征長度,m。

根據(jù)式(7)計算得到不同航速工況下的實船與模型的速度轉(zhuǎn)換結(jié)果如表2所示。

表2 實船與模型的速度轉(zhuǎn)換結(jié)果

使用計算流體力學(xué)軟件STAR-CCM+和SSTk-ω湍流模型進行仿真計算,模擬航速為26～32 kn。雙體船三維模型如圖1所示,計算時未考慮附體對船體的影響。

a)主視圖 b)左視圖

采用RANS法結(jié)合SSTk-ω湍流模型對虛擬繞流場進行求解,計算域采用長方體形狀,計算域的設(shè)置示意如圖2所示。由于船體及流動具有對稱性,所以以半船為計算對象,邊界計算定義如下:船前為船艏前1.5LPP處,速度入口;船后為船艉后5LPP處,壓力出口;船側(cè)為中線面?zhèn)确?LPP處,速度入口;船舯為船體中線面的延展面,對稱平面;頂部為設(shè)計水線以上LPP處,速度入口;底部為設(shè)計水線以下2LPP處,速度入口。

圖2 計算域與邊界條件的設(shè)置示意圖

計算域采用切割六面體網(wǎng)格,船體附近六面體網(wǎng)格的密度進行加密處理,在靜水面附近采用開爾文波形形狀的密度盒進行加密,計算域內(nèi)網(wǎng)格總數(shù)為220萬左右。計算域網(wǎng)格如圖3所示。

a)計算域內(nèi)網(wǎng)格 b)局部網(wǎng)格加密

3 模型驗證和結(jié)果分析

3.1 模型驗證

在進行模型數(shù)值模擬之前,需要預(yù)先檢驗湍流模型對于雙體船阻力計算的準確性。本文中參考文獻[14]進行湍流模型驗證,文獻[14]中使用船舶模型為雙體船,湍流模型使用修正Realizablek-ε模型。共選取7種航速工況結(jié)果進行對比,對比結(jié)果如表3所示。由表3可知:通過SSTk-ω模型模擬得到的結(jié)果相對誤差較小,相對誤差絕對值最大為5.9%,模擬誤差在可接受范圍內(nèi),SSTk-ω模型以及CFD軟件計算船舶阻力可靠。

表3 不同湍流模型下總阻力計算結(jié)果與試驗數(shù)據(jù)對比

3.2 結(jié)果與討論

使用STAR-CCM+對雙體客船進行不同航速下的靜水阻力仿真計算,航速為26～32 kn,該船在此航速階段航行時,船體姿態(tài)變化明顯。自由面的波形體現(xiàn)出船體的興波特性,數(shù)值計算收斂后,監(jiān)測到的自由面波形如圖4、5所示。

a)自由波形圖 b)船側(cè)興波云圖

a)自由波形圖 b)船側(cè)興波云圖

由圖4、5可知:相比于vc=26 kn,當(dāng)vc=32 kn時,隨著弗勞德數(shù)提高,船兩側(cè)波峰提高且范圍增加,興波范圍擴大,船艏波峰后移,船艉后興波波峰更加明顯,波峰高度有所提高。

不同工況下的船底壓力分布如圖6所示。由圖6可知:海水壓力主要作用于2個分離船體上且高壓力點集中在船體的中后方;隨著航速提高,船體受到的壓力逐漸增大,航速為26 kn時,最大船舶壓力為1 456.76 Pa,航速為32 kn時,最大船舶壓力為2 123.11 Pa,相比于航速為26 kn時增大約為45.7%。隨著海水壓力增大,作用于船體的阻力隨之增大。

圖6 雙體船在不同工況下的壓力分布對比

計算出模型阻力后,使用二因次換算法將船模阻力換算至實船,規(guī)定艉傾的縱傾角為正,艏傾的縱傾角為負;升沉用上升或下降高度與吃水的比值表示,下沉為負,上升為正。考慮到平板摩擦阻力因數(shù)區(qū)別于真正的摩擦阻力因數(shù),二因次法忽略了附體及風(fēng)阻影響,工程上通常采用的實船總阻力補償因數(shù)為10%[15]。實船有效功率

Pc=wcR′tc,

(8)

計算得到的實船總阻力及有效功率如表4所示,不同航速下實船總阻力及總阻力因數(shù)變化曲線如圖7所示。

表4 阻力計算結(jié)果

a)總阻力 b)總阻力系數(shù)

由表4及圖7可知:基于RANS法得到的實船總阻力隨實船航速的增大而增大;航行速度為32 kn時,補償后的實船總阻力達到122.77 kN,相較于航速為26 kn時的阻力提升了約22.4%。隨著航速提高,興波阻力對船體的影響顯著,在總阻力中影響占比增大;實船總阻力因數(shù)明顯隨航速增大而減小。

船舶航行工況復(fù)雜,船-機-槳達到最佳匹配狀態(tài)可以使整個系統(tǒng)的推進效果最優(yōu)[16]。船舶的有效功率

Pc=nη0δPs,

(9)

式中:Ps為主機額定功率,kW;n為發(fā)動機數(shù);η0為船機槳總效率;δ為發(fā)動機負荷率,是主機實際功率與額定功率之比。

船舶航速計算采用插值法,計算式為:

(10)

式中:uc為實際航速,kn,下標1、2表示不同航速對應(yīng)的工況點。

此雙體客船采用雙機、雙槳的推進系統(tǒng)配置,主機的額定功率為1 545 kW,船-機-槳總效率約為64%。

根據(jù)數(shù)值模擬計算得到船舶總阻力結(jié)果,匹配發(fā)動機的負載特性,可以計算得到不同發(fā)動機功率下的船舶航速。按照發(fā)動機70%～100%負荷率(1 082～1 545 kW),共計算7個發(fā)動機功率在上述條件下對應(yīng)的船舶航速,如表5所示。由表5可知:在發(fā)動機額定功率時,基本滿足客戶最大設(shè)計航速為32 kn的要求;該船主機與實船匹配效果較好,發(fā)動機負載可以覆蓋航行常用速度,可以達到節(jié)能增效的目的。

表5 不同發(fā)動機功率下所對應(yīng)的船舶航速

4 結(jié)論

通過RANS法對某雙體客船的船舶阻力進行預(yù)測,并應(yīng)用二因次轉(zhuǎn)換法獲得了航速為26～32 kn時的實船總阻力。

1)阻力隨著航速的提升而增大,當(dāng)航行速度為32 kn時,補償后的實船總阻力達到122.77 kN;總阻力因數(shù)隨著航速的提升而減小,當(dāng)航行速度為32 kn時,總阻力因數(shù)降低至2.91×10-3;在高航速工況下,雙體船具有較好的快速性能。

2)使用STAR-CCM+計算流體力學(xué)軟件,針對雙體船的阻力特性預(yù)測具有效率高、精度高和可靠性強等特點,模擬計算的阻力結(jié)果可以更好地進行船舶早期設(shè)計的船機匹配,確保整船效率得到提升。