(必維船級(jí)社(中國(guó))有限公司 先進(jìn)技術(shù)研究中心，上海 200011)

船型優(yōu)化過程中，多數(shù)研究是將船體阻力作為優(yōu)化目標(biāo)，通常采用勢(shì)流理論方法快速預(yù)報(bào)船體興波阻力，結(jié)合ITTC推薦的摩擦阻力計(jì)算方法，進(jìn)而獲得船體總阻力[1-5]。其優(yōu)點(diǎn)是計(jì)算速度快，易與優(yōu)化算法結(jié)合，可自動(dòng)搜尋最優(yōu)船型。然而，船體尾部型線不僅與阻力性能相關(guān)，還與船體的自航性能相關(guān)，僅以阻力為優(yōu)化目標(biāo)，將會(huì)忽視船型變化對(duì)自航性能的影響[6]。對(duì)自航性能進(jìn)行預(yù)報(bào)，通常是進(jìn)行物理水池船模自航試驗(yàn)，也可以利用CFD技術(shù)構(gòu)建數(shù)值水池，開展自航模擬。通常有兩種方法模擬螺旋槳與船體之間的相互作用：①直接模擬螺旋槳旋轉(zhuǎn)②采用制動(dòng)盤方法；通過在螺旋槳盤面處加入體積力，模擬螺旋槳旋轉(zhuǎn)形成的流場(chǎng)。模擬螺旋槳旋轉(zhuǎn)評(píng)估船舶自航性能，計(jì)算成本很高，很難應(yīng)用于船型優(yōu)化過程。而采用制動(dòng)盤方進(jìn)行自航模擬，即可保證自航性能預(yù)報(bào)結(jié)果精度，又可以提高計(jì)算效率。為此，利用三維船體參數(shù)化建模軟件OptNAV進(jìn)行船型變化，在基于RANSE-VOF求解器 ISIS-CFD構(gòu)建的數(shù)值水池中，開展船體拖曳阻力和自航虛擬實(shí)驗(yàn)，對(duì)型線變化前后的船體阻力和自航性能進(jìn)行預(yù)報(bào)。對(duì)比分析船體艏艉型線變化對(duì)阻力和自航性能的影響，為船體型線設(shè)計(jì)優(yōu)化提供參考。

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

(1)

采用有限體積法(FVM)對(duì)流體輸運(yùn)方程進(jìn)行空間離散。通過求解動(dòng)量守恒方程獲得速度場(chǎng)，再通過求解由連續(xù)性方法或質(zhì)量守恒條件轉(zhuǎn)換成的壓力方程獲得壓力場(chǎng)。

求解氣液兩相流動(dòng)問題時(shí)，采用VOF方處理自由液面[9]。引入流體體積分?jǐn)?shù)α。α=1表示計(jì)算網(wǎng)格單元完全被水占據(jù)；α=0表示網(wǎng)格單元完全被空氣占據(jù)；0<α<1則表示該計(jì)算網(wǎng)格單元含有水和空氣的交界面。

ISIS-CFD求解器是針對(duì)船舶海洋工程水動(dòng)力問題的專用求解器，經(jīng)過了大量標(biāo)準(zhǔn)問題求解驗(yàn)證，廣泛應(yīng)用于船舶工程的實(shí)際項(xiàng)目，可以高精度預(yù)報(bào)船體水動(dòng)力性能，模擬研究復(fù)雜流體流動(dòng)問題。DTMB5415船模(5.72 m長(zhǎng))CFD數(shù)值計(jì)算結(jié)果與水池試驗(yàn)結(jié)果的對(duì)比見圖1。

圖1 數(shù)值水池與物理水池阻力結(jié)果對(duì)比

由圖1可見，兩者吻合的很好，預(yù)報(bào)精度滿足實(shí)際工程應(yīng)用的需要。

基于ISIS-CFD求解器構(gòu)建三維數(shù)值水池，開展阻力和自航數(shù)值試驗(yàn)。在自航數(shù)值試驗(yàn)中，采用致動(dòng)盤方法[10]模擬螺旋槳與船體之間的相互作用。對(duì)比船體阻力計(jì)算結(jié)果和自航計(jì)算結(jié)果，并結(jié)合螺旋槳敞水性能數(shù)據(jù)進(jìn)行后處理，從而推導(dǎo)出推進(jìn)系數(shù)(推力減額系數(shù)(1-t)、伴流系數(shù)(1-w)、螺旋槳效率ηp、船身效率ηh)，進(jìn)一步計(jì)算得到船體收到功率PD。

2 計(jì)算結(jié)果與分析

2.1 初始船型性能預(yù)報(bào)

某內(nèi)河船船長(zhǎng)140 m，吃水6.0 m，計(jì)算航速為10.0 kn，在水深H=9.0 m的淺水航道航行。

采用長(zhǎng)方體計(jì)算域建立數(shù)值水池模型?？紤]到船體直航時(shí)流場(chǎng)的對(duì)稱性，采用半船模型。計(jì)算域長(zhǎng)為6倍船長(zhǎng)，寬2倍船長(zhǎng)，底部邊界距離水面9 m，頂部邊界距離水面1倍船長(zhǎng)。船體位于計(jì)算算域中部，船艏距離右側(cè)入口邊界2倍船長(zhǎng)。

計(jì)算域右側(cè)入口、側(cè)面和頂部邊界均為速度入口邊界條件，左側(cè)為出口邊界條件，底部為可滑動(dòng)移壁面條件，船體表面為不可滑移壁面條件，船體中縱剖面所在邊界設(shè)置為對(duì)稱面。

采用HEXPRESS軟件進(jìn)行網(wǎng)格劃分，計(jì)算網(wǎng)格數(shù)量約320萬。為精確模擬淺水水底邊界對(duì)船體的影響，如船體吸底坐沉等[11-12]，船體和水底邊界之間區(qū)域需要進(jìn)行網(wǎng)格加密。此外，螺旋槳槳盤面區(qū)域也進(jìn)行了網(wǎng)格加密，見圖2。

圖2 計(jì)算網(wǎng)格

考慮船體縱搖和垂蕩2個(gè)自由度運(yùn)動(dòng)。船體由靜止逐漸加速至目標(biāo)航速，當(dāng)船體阻力曲線穩(wěn)定后，獲取數(shù)值試驗(yàn)結(jié)果，見表1。

表1 初始船型阻力和自航性能結(jié)果

采用制動(dòng)盤法進(jìn)行自航模擬計(jì)算得到的速度場(chǎng)見圖3。船體興波高度見圖4。

圖3 速度云圖(縱向剖面Y=5.0 m)

圖4 船體興波高度云圖

由圖3可以看到，在船艉螺旋槳位置后方，有明顯的螺旋槳推水形成的高流速區(qū)。此外，船體底部和水底之間的水流相對(duì)速度在船體航速的1.1倍以上，相對(duì)較高的流速導(dǎo)致船體有明顯的坐沉現(xiàn)象。結(jié)果顯示，船體沉降值可達(dá)到0.4 m，要比在深水環(huán)境中的沉降值大很多。

2.2 船體參數(shù)化變形

采用三維參數(shù)化建模軟件OptNav對(duì)初始船型進(jìn)行型線變化，生成新船型。

對(duì)給定船型進(jìn)行艏、艉變形，見圖5。橫剖面線視圖和水線視圖見圖6。型線變化主要是改變艉和平行中體過度區(qū)的型線，使局部區(qū)域變寬或變窄。

圖5 艏型線變形

圖6 艉型線變化

2.3 型線變化的影響分析

對(duì)初始船型艏艉分別進(jìn)行4種不同幅值的變形，并在數(shù)值水池中進(jìn)行阻力和自航性能預(yù)報(bào)。通過與初始船型進(jìn)行比較，分析參數(shù)化變形對(duì)船體水動(dòng)力性能的影響。采用OptNav軟件可以將初始船型的計(jì)算網(wǎng)格自動(dòng)應(yīng)用到變形后的船型上，確保生成相同的網(wǎng)格拓?fù)浣Y(jié)構(gòu)和網(wǎng)格節(jié)點(diǎn)布置。

對(duì)于艏變形，考慮到變形區(qū)域距離艉螺旋槳距離足夠遠(yuǎn)，認(rèn)為艏變形對(duì)自航性能沒有影響，只分析船型變化對(duì)阻力的影響。對(duì)于艉變形，則需要同時(shí)考慮船型變化對(duì)阻力和自航性能的影響。

艏變形幅值范圍為[-25，30]，在此范圍內(nèi)，船體阻力變化與幅值的關(guān)系見圖7。

圖7 與艏橫剖面變形對(duì)應(yīng)的阻力相對(duì)變化

變形幅值0對(duì)應(yīng)初始船型。變形幅值為負(fù)值表示艏底部變胖，水線上部變瘦，排水體積下移。 與初始船型相比，阻力可降低；如當(dāng)變形幅值為-25時(shí)，船體阻力可降低1.6%。相反，變形幅值為正值，表示船艏底部變瘦，上部變胖，對(duì)應(yīng)阻力增大，如當(dāng)變形幅值為30時(shí)，船體阻力將增加4.0%。

對(duì)比艏部變形前后的興波高度變化見圖8。

圖8 船體興波自由面高度云圖

由圖8可見，變形幅值-25對(duì)于的船體艏部興波較初始船型行波有明顯變化，有利于降低興波阻力。

對(duì)比艏部變形前后的船體表面壓力系數(shù)變化見圖9。變形幅值-25對(duì)應(yīng)的艏高壓區(qū)范圍明顯減小。

圖9 船艏表面動(dòng)壓力系數(shù)云圖

艉變形幅值為[-30，30]，不同變形幅值對(duì)應(yīng)船型的阻力和自航性能參數(shù)見表2。與初始船型相比，阻力和收到功率變化見圖10。結(jié)果顯示，當(dāng)變形幅值為15時(shí)，阻力較初始船型降低了1%，但收到功率PD僅降低了0.3%。這主要是由于艉變形改變了艉伴流系數(shù)和推力減額系數(shù)，從而減小了船身效率和螺旋槳效率，導(dǎo)致收到功率并沒有阻力收益那么高。

圖10 船艉變形對(duì)應(yīng)的阻力和收到功率的相對(duì)變化

與初始船型相比較，在給定的變形區(qū)域內(nèi)，船體略變寬些時(shí)，有利于降低阻力，如變形幅值15對(duì)應(yīng)的船型；而將船體略變窄些時(shí)，船身效率和螺旋槳效率會(huì)有所改善，如變形幅值-15所對(duì)應(yīng)的船型。因此，在對(duì)船艉進(jìn)行優(yōu)化時(shí)，僅以最小阻力為優(yōu)化目標(biāo)是不夠的，更要注意對(duì)自航性能的準(zhǔn)確預(yù)報(bào)。

表2 船艉變形幅值對(duì)應(yīng)的阻力和自航性能結(jié)果

由于艉部變形幅度不大，而且變化幅度較大的部分在距離水線面較深的位置，所以對(duì)艉部興波影響很小，主要改變變形區(qū)域的壓力分布，見圖11。

圖11 船艉表面壓力場(chǎng)云圖

3 結(jié)論

船體艏艉局部范圍均有進(jìn)一步優(yōu)化的空間；船艉變形對(duì)船體阻力和自航性能的影響并不完全一致；對(duì)船艉進(jìn)行優(yōu)化時(shí)，以自航計(jì)算獲得收到功率為優(yōu)化目標(biāo)更為合適。