彭建鋒 楊敬東

【摘 要】隨著CFD（Computational Fluid Dynamics）技術在船舶設計研究中的應用，數(shù)字化的船模水池，在民船與軍船水動力特性研究中得到廣泛的應用。本文介紹了以勢流理論為基礎的船舶設計軟件Shipflow，并介紹了其適用的船舶類型，分析了它在阻力預報中的優(yōu)越性。

【關鍵詞】CFD；勢流理論；阻力預報

中圖分類號： U675.91 文獻標識碼： A 文章編號： 2095-2457（2018）16-0194-002

DOI：10.19694/j.cnki.issn2095-2457.2018.16.088

【Abstract】With the application of CFD（Computational Fluid Dynamics）technology in ship design research， digitalized model vessel pools have been widely used in the study of hydrodynamic characteristics of civilian and military ships.

【Key words】CFD；Potential flow theory；Resistance prediction

21世紀，科技高速發(fā)展為船舶業(yè)帶來了巨大的機遇與挑戰(zhàn)，作為船舶設計者，要審時度勢，積極開發(fā)應用新技術，才能立足于行業(yè)前列。CFD是一種新興的計算方法，它作為一種新的優(yōu)化船型手段在國內外船舶設計中已得到了廣泛的應用。但由于受到當前理論水平和計算機的運算速度限制，CFD方法對于精準預報船舶水動力性能還些不足，但其計算結果仍能一定的精度要求，且能得到比船模試驗更多流場信息。目前主流的CFD設計軟件有：Fluent、CFX，Star-CCM、Comsol等，由于算法和軟件本身側重點不同，因此各有優(yōu)缺點。本文將介紹一種基于勢流理論的CFD計算軟件Shipflow，并分析其在特定船型設計中所具有的優(yōu)越性。

1 軟件簡介

Shipflow是以CDF為依托，處于領先地位的船舶水動力設計分析專用軟件，最初是由瑞典SSPA公司和Chalmers科技大學在八十年代聯(lián)合研發(fā)的，針對船體水動力數(shù)值模擬的專用軟件[1]。在國內，Shipflow應用主要集中在高校等科研院所，側重于應用研究[2]；在國外，該軟件側重應用于多學科設計優(yōu)化，例如，船型參數(shù)化軟件Friendship已將其作為系統(tǒng)默認求解器[3]。除此之外國外學者還十分重視軟件的檢驗與驗證，Lu[4]等對Shipflow軟件用于油船仿真的網(wǎng)格劃分、計算設置等進行了相關性分析；Mierlo[5]對Shipflow軟件的網(wǎng)格劃分、興波阻力計算方法和收斂準則等開展了較為深入的研究，并經過系列船型計算，得到了可用于修正Shipflow計算結果的回歸公式。

2 計算優(yōu)越性分析

船舶拖曳水用船舶模型實驗方法，模擬船舶運動狀態(tài)，是船舶水動力實驗的重要設備。船模實驗是驗證船舶設計、了解船舶阻力性能的重要環(huán)節(jié)，其實驗數(shù)據(jù)至今仍在新船型開發(fā)中起著重要的作用[6]。但這種方法實驗費用高，周期長，且對船池本身依賴性較大。船舶CFD是通過模擬不可壓縮牛頓流體的湍流運動狀況來預報船舶航行時的水動力性能。結合具體的船型，CFD可以比較精確預報船舶快速性和流場特性，能比較不同流線方案的快速性優(yōu)劣。由于是在“數(shù)值水池”中進行試驗，尤其在設計初期，不用進行船模實驗，也可以獲得具有參考價值的設計參數(shù)，從而縮短設計周期，降低成本，這在船舶設計有著重要意義。Shipflow是造船行業(yè)廣泛認同的功能最全面，最完善的有效計算船舶阻力性能的CFD軟件之一，它相當于數(shù)字化的船模水池，在民船與軍船水動力特性研究中有廣泛的應用。Shipflow可以通過模擬給出波浪模式，壓力分布，速度矢量，空間流線和波浪增阻，粘性阻力，興波阻力，誘導阻力，升力，螺旋槳效率等船體特征性參數(shù)[7]。新開發(fā)的Shipflow版本還增加了增添了淺水效應研究，盤形推進器網(wǎng)格自動劃分和多船體研究等功能。

Shipflow將船體周圍的流場劃分為三個計算流域，分別為勢流區(qū)、薄邊界層區(qū)和湍流區(qū)，根據(jù)各個區(qū)域流體流動特性的不同，分別在每一區(qū)采用不同的數(shù)值計算方法。其中每一區(qū)使用的CFD方法如下：

2.1 勢流區(qū)

Shipflow在該區(qū)域采用Rankine源法，根據(jù)自由表面的線型、非線性邊界條件，進行高階面元法的勢流理論進行計算，主要是求解船舶的興波阻力。

其原理是認為此區(qū)域船舶在不可壓縮的理想流體中運動，忽略水的粘性，且不考慮尾部伴流影響，進行數(shù)值計算時，運用高階面元法，聯(lián)立Laplace方程以及自由液面邊界條件得出速度勢分布，然后根據(jù)伯努利方程得到壓力與速度的關系，最后算出船舶興波阻力系數(shù)Cw。

2.2 薄邊界層區(qū)

此區(qū)域主要應用邊界層理論。采用動量積分的方法，對船舶前三分之二的船體沿邊界層進行厚度積分，起點可以從駐點開始，也可以從選定的站點開始。將勢流壓力分布作為輸入?yún)?shù)，可計算出此區(qū)域的摩擦阻力。

2.3 湍流區(qū)

此區(qū)域可看做是對粘性流體進行求解。Shipflow采用基于有限體積法的RANS求解器——XCHAP求解模塊[8]來進行計算。

計算中采用時均化的N-S方程，即雷諾方程。

不可壓縮流體的雷諾方程：

在數(shù)值計算中采用有限體積法，通過k-？著湍流模型進行求解。其中，RANS方程的邊界條件可由勢流和邊界層計算得到，方程中的對流項采用Roe平均策略，取二階精度格式，其余各項采用中心差分進行離散，并對離散方程進行迭代求解[9]，由于此區(qū)域船體受到流體的粘性作用，對船體濕表面進行積分，所得運動方向上的分量即粘性阻力Rv，再通過RANS方程求解得到船體表面壓力分布。采用這種分區(qū)域分步驟的求解方式，得到船舶航行時受到的總阻力。整個計算方法既能保證計算的精度與準確度，同時也有效縮短了計算周期，提高了計算效率。

3 阻力預報方法

艦船船體阻力可以分為興波阻力與粘性阻力，阻力預報有兩種方法。

第一種，分別計算兩種阻力，其中粘性阻力利用Holtrop方法估算，形狀因子可用經驗公式，摩擦阻力系數(shù)可用柏-徐公式計算。興波阻力計算可以通過對船體表面勢流壓力進行積分得到；粘性阻力是摩擦阻力與粘壓阻力之和，摩擦阻力是對沿船體表面局部摩擦應力積分得到，前半體局部摩擦阻力通過邊界層方法得到，后半體局部摩擦阻力通過N-S方程中壁面剪應力得到，粘壓阻力由粘性引起的壓力積分得到。

第二種是排擠厚度法，排擠厚度是指粘性作用造成邊界層速度降低，相比理想流體有流量損失，相當于中心區(qū)理想流體的流通面積減少，計算時將平板表面上移一個厚度，此為排擠厚度[11]。在計算興波阻力時就等效的計入粘性的影響，獲得剩余阻力，然后估算出艦船的船體阻力。

采用和這兩種方法估算剩余阻力時，只有當剩余阻力系數(shù)Crs>0.32時預報準確性高，系數(shù)較小時，則存在一定的誤差。其中采用排擠厚度法比第一種估算阻力方法更為準確，其預報阻力誤差基本可控在7%內。計算興波阻力時引入排擠厚度等效算作粘性影響，進行修正，可以達到對船舶阻力的快速預報。

4 關于Shipflow應用的討論

近年來，有關CFD的研究進展較快，Shipflow的應用也越來越廣泛。但是，無論是網(wǎng)格處理、求解器設計、后處理，還是計算結果可靠性分析上，由于源代碼不公開而存在修正的困難[12]。大量的算例分析得出，Shipflow對于大Cb，低Fn船型有較好的計算結果，但是對于特殊的高速船，其計算結果不如Fluent；尤其當船舶航態(tài)處于半滑行狀態(tài)時，出現(xiàn)破波，其計算精度往往無法滿足設計要求，只能用船模試驗。綜上，對于船型幾何形狀及船體周圍流動的過于復雜時，數(shù)值計算結果只適用于對船舶性能做定性的分析。若能解決這幾方面問題，Shipflow將會得到更為廣泛的應用。

5 結束語

本文對目前應用較為廣泛的Shipflow軟件進行了介紹，分析了其分區(qū)域計算的方法，以及其在船舶阻力預報中的實際運用。通過對船舶繞流的分析，對船舶阻力的計算，可以對船舶設計進行線型優(yōu)化，對初步線型的阻力性能進行對比、評價，為修改線型優(yōu)化設計提供幫助[13]，提高船舶快速性，從而降低船舶制造成本，提高工作效率，縮短建造周期，節(jié)省船模加工、試驗成本。

【參考文獻】

【References】

[1]陳紅梅，蔡榮全，勢流計算在船舶型線優(yōu)化改型中的適用性研究[J].船舶工程，2012，S2：9-19.

[2]刑瑤，易宏，李英輝，改進型雙支柱SWATH的興波阻力研究[J]，船舶工程，2008，30（2）：16-19.

[3]CLAUS A，HARRIES S，WUHDERLICH，et al.Flexible tool integration for simulation-driven design using XML，generic and COM interfaces [C]//8th International Conference on Computer and IT Applications in the Maritime Industrise.Budapest，2009：461-473.

[4]LU Z，LARSSON L，ORYCH M.Verification and validation of CFD predictions for a manoeuvring tanker[C]//9th international Conference on Hydrodynamics.Shanghai，2010：438-445.

[5]MIERLO K J.Trend validation of SHIPFLOW based on the bare hull upright resistance of the delft series[D].Netherlands：Delft University of Technology，2006.

[6]周志勇，CFD在船型優(yōu)化中的應用[J]，船舶設計通訊，2005，02，46-49.

[7]吳曉蓮，基于CFD的船舶球艏-球艉低阻線型研究[D]. 哈爾濱：哈爾濱工程大學，2008，02：5-6.

[8]盧雨，胡安康，尹遜濱，田明琦，基于SHIPFLOW的某大型集裝箱船阻力預報與實驗驗證[J]，中國造船，2016，02，22-30.

[9]胡春平，基于Friendship的船體型線優(yōu)化技術研究[D]. 武漢：武漢理工大學，2012，5：48-50.

[10]陳偉，許輝，邱遼原，姜治芳，王曉喆，基于SHIPFLOW軟件的方尾艦船阻力快速預[J].中國艦船研究，2012，04，17-22.

[11]夏國澤，船舶流體力學[M]，武漢，華中科技大學出版社，2003.

[12]萬德成，船舶與海洋工程數(shù)值水池研究[M]，中國力學學會學術大會，2009.

[13]魏菲菲，李曉峰，49000DWT油船化學品船線型設計與優(yōu)化[J]，船舶設計通訊，2010，09，S2，16-20.