基于CFD的双艉客船阻力特性分析

赵丙乾，方昭昭，喻露，陈庆任

1中国船级社武汉规范研究所，湖北武汉430022

2中国舰船研究设计中心，湖北武汉430064

基于CFD的双艉客船阻力特性分析

赵丙乾1，方昭昭2，喻露1，陈庆任1

1中国船级社武汉规范研究所，湖北武汉430022

2中国舰船研究设计中心，湖北武汉430064

［目的］随着绿色船舶规范的实施及推进，船舶能效指数，特别是船舶航速功率指标越来越受到业界的关注。为了精确获得双艉船型的航速功率指标及船舶流场信息，［方法］基于计算流体动力学（CFD）理论，提出一种双艉船型航行阻力预报及流场捕捉的工程方法。采用粘性流计算软件FINE/MARINE，对某双艉客船的航行阻力进行计算，探讨船舶第1层网格节点高度、船舶航速以及附体对计算结果的影响，并将不同航速下的阻力预报及其结果与试验数据进行对比。［结果］结果显示，预报误差基本在3%以内。［结论］研究表明，所提方法计算效率高且易于实现，预报精度能满足工程需要，具有较强的工程实用性。

双艉客船；计算流体动力学；阻力预报；FINE/MARINE

0 引 言

随着绿色船舶规范的实施以及内河标准船型建设的推进，内河船舶的能效指标越来越受到业界的关注。交通运输部更是要求［1］：对于新建、改建的内河运输船舶，应当满足内河运输船舶标准船型指标体系中的强制性指标要求。强制性指标［2］包括船舶主尺度系列标准、CO2排放指标和燃料消耗指标，其中CO2排放指标和燃料消耗指标均与船舶能效指数相关联。船舶能效指数评估包括能效指数前期评估和建造后评估2类［3-4］。用于能效设计指数前期评估的功率曲线一般基于水池试验。在船型设计初期，可以通过水池试验评估船舶快速性能，以获得优良的船型。但进行船模试验需要制作船体模型，费用高、周期长，且不易对船型进行多方案调整。因此，在船舶设计初期，如何精确预报船舶阻力，获得流场信息，以掌握船舶快速性性能指标并方便后续阶段的船舶优化改进是船舶设计和能效指标验证的关键。

随着计算机技术的飞速发展和计算数学理论的不断完善，计算流体动力学（Computational Fluid Dynamics，CFD）得到了蓬勃发展，成为研究船舶航行性能的重要手段之一。相比于模型试验方法，采用CFD方法不仅费用低、周期短，而且还可以准确获得船舶周围流场信息。目前，国内一些单位采用基于RANS方程加自由面的处理方法对船舶快速性能进行的预报在工程上基本能满足其精度要求［5-8］，但以上这些研究大多是基于海船船型。受吃水及航道的限制，内河船舶多采用双艉或者三艉设计，船型宽大而型深较小，主要船型系数与海船相比差异较大，而目前有关双艉船型阻力预报及流场捕捉方面的研究还较少。

FINE/MARINE是船舶与海洋工程水动力专业使用的粘性流计算软件［9］，主要包括网格生成器HEXPRESS、后处理工具CFView和不可压粘性流场求解器ISIS-CFD。其特色在于，HEXPRESS为全六面体非结构网格生成工具，采用体到面的网格生成方式可以快速生成高质量的网格，并且其对船舶海洋工程问题的模拟可以通过界面设置来实现，无需进行二次开发，方便快捷。

本文将基于计算流体动力学理论，提出一种双艉船型航行阻力预报及流场捕捉的工程方法。采用FINE/MARINE软件，对某双艉客船的航行阻力进行计算，探讨船舶第1层网格节点高度、船舶航速以及附体对计算结果的影响。通过不同航速下的阻力预报，并将数值结果与试验数据进行比较，以验证所提方法的正确性。

1 数学模型

本文对水面船舶自由面扰流问题的数值模拟是将自由面流动作为两相流来处理，自由面为水和空气的交界面，采用流体体积法（Volume of Fluid，VOF）求解。水面船舶自由面扰流问题数学模型的控制方程包括：连续性方程、体积分数方程、动量方程和湍流模型方程。

1.1 控制方程

不可压缩粘性流体的连续性方程和动量方程分别为：

式中：t为时间；ρ为密度；v为控制体；A为控制体面积；n为法矢量；Ud为控制体表面法矢量方向上的速度；U为速度矢量；p为压力；Ui为xi坐标轴方向上的平均速度分量；τij为粘性应力张量；gi为重力矢量；Ii和Ij为方向向量。

1.2 湍流模型

本文采用SSTk-ω湍流模型，其湍动能k和湍流耗散率ω的输运方程如下：

式中：μ为分子粘度；Uj为xj坐标轴方向上的平均速度分量；μt为湍流涡粘度；tij为湍流雷诺应力张量；Sij为平均应变率张量；FI为辅助混合函数；Pω为ω的导出项；β为湍流应力常数；β*，σk和σω2

均为为湍流模型常数。

1.3 自由面捕捉算法

自由液面捕捉法是将空气和水作为单一流体同时计算。该单一流体的性能（质量体积和粘性因数）在空间的变化取决于构成函数ci。在自由液面计算中，ci在空气中的取值为0，在水中的取值为1，通过求解下式确定：

自由液面捕捉算法具有更好的灵活性和适应性，可以较好地模拟破碎波等复杂的自由液面。

1.4 离散格式

控制方程的离散项采用隐式有限体积法，直接求解三维粘性不可压缩多相流的雷诺平均方程，其具有二阶空间和时间精度。动量方程离散采用GDS格式，时间离散采用时间步进算法。自由液面捕捉采用BRICS可压缩型离散格式，能减小自由液面附近构成函数的数值扩散。

2 研究对象

研究对象为航行于三峡库区的某双艉客船，该客船为双机、双桨、双舵，船舶线型为双艉线型、V型船艏。该船由于船型特殊，缺乏相关的资料，为准确预报其船舶阻力，在武汉理工大学的拖曳水池进行了阻力测试。本文对船模进行了CFD模拟计算，主要船型参数如表1所示，试验船模及布置图如图1所示。

表1 双艉客船船型主要参数Table 1 Main parameters of the twin-skeg passenger ship

图1 船模及试验布置图Fig.1 Model ship and test arrangement

3 数值计算

3.1 计算区域及边界条件

考虑到船型及流动的对称性，进行数值模拟时只计算一半区域即可。计算区域如图2所示，其中，计算边界定义如下：

1）前端——模型首部前约1.0LPP处；

2）后端——模型尾部后约3.0LPP处；

3）侧边界——模型侧方约1.5LPP处；

4）上边界——水线以上约0.5LPP处；

5）下边界——水线以下约1.5LPP处；

6）对称面——模型中纵剖面的延展面。

英国的体育课程内容以基础教育为例，就是学习领会国家课程标准所规定的具体内容，即：舞蹈、游戏、体操、游泳和水上安全、田径、户外、探险，由于不同的学段所要达到的目标不同，因此教师在进行教学时应根据学生所处的学段选择相应的内容。课程标准规定了体育课程学习领域，即：学生必须学习舞蹈、游戏、体操、游泳和水上安全、田径和户外探险活动等6大范畴内容，其中，学段一学生（5～7 岁）对游泳和水上安全选修，不学田径和户外探险活动，其他学段要求必修上述六大领域课程内容。

图2 数值计算区域及网格划分Fig.2 Numerical calculation domain and mesh scheme

计算域的上、下边界取为压力边界，前、后截面及远离船体的一侧取为远流场边界，靠近船体的一侧取为对称边界，船体甲板面为滑移物面，其他表面为不可滑移物面。因在求解过程中存在着时间偏导项，故船体从静止到某一航速，给定加速时间，然后船舶按照设计航速航行直至收敛。

3.2 网格划分

计算网格质量影响数值计算的正确性和精确性。本文采用软件自带的前处理器HEXPRESS生成全六面体非结构化网格。通过将网格粗化、细化、吸附以及在自由面附近进行细化形成贴体网格。其中，第1层边界层网格间距及边界层层数根据无量纲特征数Y+确定。分别对双艉客船带舵和不带舵这2种状态下的阻力进行计算，其中带舵的计算网格总数为82万，不带舵的网格总数为76万。网格质量通过正交性来保证，即绝大部分网格的正交性为90°。其中，不带舵网格的最小正交性≥22.5°，对于带舵的网格，由于舵尾缘尺寸较小，故最小正交性≥8°，网格质量满足数值计算要求。船体表面网格划分如图3所示。

图3 船体表面网格划分Fig.3 Mesh scheme of the hull

4 计算结果分析

为方便对结果的处理，引入了无量纲系数：傅汝德数Fr、雷诺数Re、船舶总阻力系数CT和摩擦阻力系数CF_ITTC，分别定义如下：

式中：Vm为船模速度；g为重力加速度；ρw和μw分别为水的密度与粘度；Aw为船舶静水湿表面积；R为船模总阻力。

4.1 对带舵和不带舵的影响分析

双艉客船的水池试验考虑带舵的影响。在进行CFD数值计算时，为便于比较分析舵对船艉流动的影响，分别进行了裸船体和带舵船体2种状态下的多航速计算，其中Fr=0.099～0.217。总阻力系数的数值计算结果及其与试验数据的比较如表2所示。表中：CT_EXP为水池试验结果，考虑了带舵的情况；CT_CAL1为裸船体阻力计算结果；CT_CAL2为带舵船体阻力计算结果。定义误差：

表2 双艉客船不同航速下阻力计算结果及其与试验数据的对比Table 2 Comparison of numerical results and experimental data for total resistance of the twin-skeg passenger ship at different speeds

由表2可以看出，不带舵和带舵2种状态下的阻力数值计算结果与试验结果均吻合良好，误差在3%以内。在不考虑舵影响的计算中，也即裸船体时，虽然网格质量和计算效率都很高，但计算结果与试验结果相比普遍偏小。相比于裸船体的计算，考虑带舵计算时网格数量会有所增加，网格划分难度也较大，但其计算结果与试验结果更为吻合。图4和图5所示分别为裸船体和带舵船体在Fr=0.178下航行时的船艉流场对比。从图中可以看出：一方面，舵作为附体会增加船舶阻力；另一方面，也会改变船艉的流线分布。为更真实地反映水池试验状态，在后续的算例中，将均采用带舵船体进行计算分析。

图4 Fr=0.178时裸船体艉部流线图Fig.4 Flow visualization of the ship stern without rudder at Fr=0.178

图5 Fr=0.178时带舵船体艉部流线图Fig.5 Flow visualization of the ship stern with rudder at Fr=0.178

在进行船舶CFD绕流问题模拟时，通常将船体表面视为无滑移壁面。根据流动特性的不同［10］，壁面边界层可以分为3个区域：粘性底层、过渡层和惯性子层。边界层的位置通常用无量纲特征数Y+来表示。无量纲分析和试验工作表明，在边界层特定的Y+范围内，具有特定的流动特点。对大多数高雷诺数流动来说，一般在粘性底层内使用壁面函数，以节约计算资源。在惯性子层区域，由于湍流应力近乎常数，故主流方向的流速分量呈对数分布。因此，为捕捉壁面附近流动的物理特性，可在物理特征变化较大的地方分布大量网格点，尤其是第1层网格点应落在此区域，且必须布置足够的网格数以保证分辨所有的流动梯度。这就要求正确布置边界层第1层网格点的位置以及边界层层数。对于无量纲特征数Y+，可采用下式进行估算：

式中：ywall为第1层网格距离；Vref为来流的参考速度；Lref为参考长度；υ为粘性系数。

针对本文所研究的双艉船舶，以Vm=1.721 m/s为例，对应的Fr=0.237，Re=9.781×106，该航速下对应的水池试验结果为CT_EXP=4.247×10-3。研究此航速下船舶第1层网格节点高度以及网格边界层层数对计算结果的影响，其中船舶初始网格设置及网格细化准则均与前述带舵航行时船舶阻力计算的设置相同。分别设置船舶边界层网格层数为6～10层，网格变化率为1.2，逐步改变船舶第1层网格节点高度和边界层层数，以查看边界层网格方案对计算结果的影响。计算结果如表3所示。

表3 不同网格方案计算结果Table 3 Simulation results of different mesh scheme

由表3可以看出，随着边界层第1层网格点位置的减小，阻力计算结果有所增加，而Y+则是逐步减小。当第1层网格节点距离相同时，不同的边界层网格数对计算结果的影响也不一样，较多的边界层网格数能够更好地捕捉边界层内速度梯度的变化，计算结果将更加精确。对于船体近壁面流动计算，推荐的Y+取值范围为20～60，边界层网格数约为8～10层。

4.3 不同航速下阻力的数值计算分析

双艉客船的设计航速范围较大，在进行船舶阻力模拟计算时，为进一步验证文中数值模拟方法的可靠性与可行性，对较高航速下的船舶总阻力也进行了计算，例如Fr取为0.237，0.257和0.276时。在高航速的数值计算中，考虑了最优Y+取值及边界层数量。表4给出了带舵双艉客船在全航速范围内（Fr=0.099～0.276）不同航速下的总阻力数值计算结果，并与试验值进行了比较分析。

表4 双艉客船阻力计算结果Table 4 Resistance calculated results of the twin-skeg passenger ship

由表4可以看出：在所计算的速度范围内，尽管航速变化较大，但和水池试验结果相比，数值计算结果的误差均在3%以内，说明本文中数值模拟精度较高，能满足工程需要。

不同航速下，船舶各阻力占总阻力的百分比不一样。一般对水面排水型舰船而言，低速航行时，兴波较小，摩擦和粘压阻力占主要成分。而当航速增加时，船体兴波成分会逐渐增加，船舶兴波阻力随之增大。图6所示为不同航速下船艏兴波高度的变化。从中可以看出，随着船舶航速的增加，船艏兴波高度也相应增大，因此在双艉客船高速航行时，应考虑船艏消波的必要性。

图6 不同航速下船艏兴波高度的变化Fig.6 Variation of wave height of the bow at different speeds

图7所示为不同航速下的自由液面波形图。从中可以看出，不同航速下船舶自由液面的波形差别较大。这同时也表明文中所提计算方法是正确的，能准确反映物理现象，可为后续船型的改进提供参考。

图7 不同航速下船舶自由面波形Fig.7 Free surface wave of the ship at different speeds

5 结 语

本文基于计算流体动力学理论，提出了一种双艉船型航行总阻力及流场信息捕捉的工程方法。采用粘性流计算软件FINE/MARINE，对某双艉客船的阻力特性进行预报分析，探讨了船体壁面第1层网格点的位置及边界层层数、船舶航速以及附体对总阻力计算结果的影响。最后，对某双艉客船带舵船体在全航速范围内的总阻力进行了数值计算，并与试验数据进行了比较，结果显示，误差均在3%以内，且阻力是随着航速的增加而增大。研究表明，本文所提方法计算效率较高，易于实现，预报精度较高，能满足工程需要，具有较强的工程实用性。

［1］交通运输部.内河运输船舶标准化管理规定［S］.北京：交通运输部，2014.

［2］交通运输部.内河运输船舶标准船型指标体系［S］.北京：交通运输部，2012.

［3］交通运输部.内河船舶能效设计指数（EEDI）评估指南2012［S］.北京：人民交通出版社，2012.

［4］中国船级社.绿色船舶规范［S］.北京：人民交通出版社，2012.

［5］黄少锋，李百齐，陈京谱.CFD在船舶工程中的实际运用［C］//第二十一届全国水动力学研讨会暨第八届全国水动力学术会议暨两岸船舶与海洋工程水动力学研讨会论文集.济南：中国力学学会，中国造船工程学会，2008：323-332.

［6］吴乘胜，赵峰，张志荣，等.基于CFD模拟的水面船功率性能预报研究［J］.中国造船，2013，54（1）：1-11.WU C S，ZHAO F，ZHANG Z R，et al.Predictionof powering performance for a surface shipbased on CFD simulations［J］.Shipbuilding of China，2013，54（1）：1-11（in Chinese）.

［7］方昭昭，赵丙乾，陈庆任.基于CFD的双艉船型阻力数值预报［J］.中国舰船研究，2014，9（4）：36-44，62.FANG Z Z，ZHAO B Q，CHEN Q R.CFD theory based resistance prediction for ships with twin-skeg［J］.Chi⁃nese Journal of Ship Research，2014，9（4）：36-44，62（in Chinese）.

［8］邓贤辉，方昭昭，赵丙乾.基于计算流体动力学的最小阻力船型自动优化［J］.中国舰船研究，2015，10（3）：19-25.DENG X H，FANG Z Z，ZHAO B Q.Automatedoptimi⁃zation of the ship hull form with minimum resistance⁃based on CFD［J］.Chinese Journal of Ship Research，2015，10（3）：19-25（in Chinese）.

［9］魏子凡，井升平，杨松林.新型高速艇的CFD模拟和对比分析［J］.中国舰船研究，2016，11（4）：22-28.WEI Z F，JING S P，YANG S L.CFD simulation and comparison analysis of a new type high-speed boat［J］.Chinese Journal of Ship Research，2016，11（4）：22-28（in Chinese）.

［10］刘应中，缪国平.高等流体力学［M］.2版.上海：上海交通大学出版社，2002.

Ananysis of resistance characteristics of twin-skeg passenger ship based on CFD

ZHAO Bingqian1，FANG Zhaozhao2，YU Lu1，CHEN Qingren1
1 Wuhan Rule and Regulation Research Institute，China Classification Society，Wuhan 430022，China
2 China Ship Development and Design Center，Wuhan 430064，China

With the implementation and promotion of Rules for Green Ships,more and more attention has been paid by the maritime industry to the ship energy efficiency index.In order to evaluate the function relationship between speed and effective horsepower,and acquire the accurate flow field information of twin-skeg ships,a method for resistance prediction of twin-skeg ships based on the Computational Fluid Dynamics(CFD)theory is outlined.Numerical simulation and resistance prediction by FINE/MARINE for a twin-skeg passenger ship advancing at different speeds are carried out respectively.The effects of the mesh characteristics,calculating speed range,and existence of appendages such as rudders on resistance are discussed.Comparisons between present results and corresponding experimental are made,showing that the prediction error can be controlled below 3%.The results demonstrate that the proposed method is efficient enough to be realized and will meet the engineering requirements.

twin-skeg passenger ship；Computational Fluid Dynamics（CFD）；resistance prediction；FINE/MARINE

U661.31+1

：ADOI：10.3969/j.issn.1673-3185.2017.03.004

http://kns.cnki.net/kcms/detail/42.1755.TJ.20170512.1258.034.html期刊网址：www.ship-research.com

赵丙乾，方昭昭，喻露，等.基于CFD的双艉客船阻力特性分析［J］.中国舰船研究，2017，12（3）：23-28.

ZHAO B Q，FANG Z Z，YU L，et al.Ananysis of resistance characteristics of twin-skeg passenger ship based on CFD［J］.Chinese Journal of Ship Research，2017，12（3）：23-28.

2016-12-20< class="emphasis_bold">网络出版时间

时间：2017-5-12 12:58

中国舰船研究设计中心研发基金资助项目

赵丙乾（通信作者），男，1985年生，硕士，工程师。研究方向：船舶阻力快速性。E-mail：zhaobingqian@ccs.org.cn

方昭昭，女，1986年生，硕士，工程师。研究方向：舰船水动力学计算。E-mail：281988856@163.com

喻露，女，1990年生，硕士，工程师。研究方向：船舶三维模型及船型优化。E-mail：luy@ccs.org.cn

陈庆任，男，1984年生，博士，高级工程师。研究方向：船舶与海洋工程水动力学。E-mail：qrchen@ccs.org.cn