Chinese Journal of Ship Research
CFD基于 的双艉客船阻力特性分析
1，方昭昭2，喻露1，陈庆任1赵丙乾1 430022中国船级社 武汉规范研究所，湖北 武汉2 430064中国舰船研究设计中心，湖北 武汉
摘 要：［目的］随着绿色船舶规范的实施及推进，船舶能效指数，特别是船舶航速功率指标越来越受到业界的关注。为了精确获得双艉船型的航速功率指标及船舶流场信息，［方法］基于计算流体动力学（CFD）理论，提出FINE/MARINE，对某双艉客船的航一种双艉船型航行阻力预报及流场捕捉的工程方法。采用粘性流计算软件1行阻力进行计算，探讨船舶第 层网格节点高度、船舶航速以及附体对计算结果的影响，并将不同航速下的阻力3%以内。［结论］研究表明，所提方法计预报及其结果与试验数据进行对比。［结果］结果显示，预报误差基本在算效率高且易于实现，预报精度能满足工程需要，具有较强的工程实用性。关键词：双艉客船；计算流体动力学；阻力预报；FINE/MARINE中图分类号：U661.31+1 文献标志码：A DOI：10.3969/j.issn.1673-3185.2017.03.004 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
Abstract：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. Key words：twin-skeg passenger ship；Computational Fluid Dynamics（CFD）；resistance prediction； 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 控制方程
不可压缩粘性流体的连续性方程和动量方程分别为： ¶ )× （1） v A ρ(U - U ndA =0 ρdv + ¶t d ¶ )× v A (U - U ndA= ρUi dv + ρUi ¶t d (τij )× A v （2） Ij - pIi ndA + ρg dv i式中：t为时间；ρ 为密度；v 为控制体；A为控制体面积；n 为法矢量；U 为控制体表面法矢量方d
向上的速度； U 为速度矢量； p 为压力； Ui 为 xi坐标轴方向上的平均速度分量；τij 为粘性应力张量；g 为重力矢量；Ii 和 I 为方向向量。i j
1.2 湍流模型
SST本文采用 k - ω湍流模型，其湍动能 k 和湍流耗散率ω的输运方程如下： ¶ρk ¶k + ¶ = tij S - β* ρωk（3） ρU k - ( μ +σ μ t) ¶t ¶xj j k ¶xj ij
¶ρω ¶ω + ¶ ρU ω - ( μ +σ μ t) =P - ρβω2 + ¶t ¶xj j k ¶xj ω ρσω2
¶k ¶ω （4）
2(1 - F ) I ω ¶xj ¶xj式中：μ 为分子粘度；Uj 为 xj 坐标轴方向上的平均速度分量；μ 为湍流涡粘度；t 为湍流雷诺应t ij
力张量；S 为平均应变率张量；F 为辅助混合函ij I数；Pω 为 ω 的导出项； β 为湍流应力常数； β* ， σk 和 σω2 均为为湍流模型常数。
1.3 自由面捕捉算法
自由液面捕捉法是将空气和水作为单一流体同时计算。该单一流体的性能（质量体积和粘性因数）在空间的变化取决于构成函数 ci 。在自由液面计算中，ci 在空气中的取值为0，在水中的取1，通过求解下式确定：值为
¶ )× vci Aci (U - U =0 （5） dv + ndA ¶t d自由液面捕捉算法具有更好的灵活性和适应性，可以较好地模拟破碎波等复杂的自由液面。
1.4 离散格式
控制方程的离散项采用隐式有限体积法，直接求解三维粘性不可压缩多相流的雷诺平均方程，其具有二阶空间和时间精度。动量方程离散GDS采用 格式，时间离散采用时间步进算法。自BRICS由液面捕捉采用 可压缩型离散格式，能减小自由液面附近构成函数的数值扩散。
2 研究对象
研究对象为航行于三峡库区的某双艉客船，该客船为双机、双桨、双舵，船舶线型为双艉线型、V型船艏。该船由于船型特殊，缺乏相关的资料，为准确预报其船舶阻力，在武汉理工大学的拖曳CFD水池进行了阻力测试。本文对船模进行了 模1拟计算，主要船型参数如表 所示，试验船模及布1置图如图 所示。
3 数值计算3.1 计算区域及边界条件考虑到船型及流动的对称性，进行数值模拟
2时只计算一半区域即可。计算区域如图 所示，其中，计算边界定义如下： 1）前端——模型首部前约1.0LPP处； 2）后端——模型尾部后约3.0LPP处； 3）侧边界——模型侧方约1.5LPP处； 4）上边界——水线以上约0.5LPP处； 5）下边界——水线以下约1.5LPP处； 6）对称面——模型中纵剖面的延展面。 计算域的上、下边界取为压力边界，前、后截面及远离船体的一侧取为远流场边界，靠近船体的一侧取为对称边界，船体甲板面为滑移物面，其他表面为不可滑移物面。因在求解过程中存在着时间偏导项，故船体从静止到某一航速，给定加速时间，然后船舶按照设计航速航行直至收敛。
3.2 网格划分
计算网格质量影响数值计算的正确性和精确HEXPRESS性。本文采用软件自带的前处理器 生成全六面体非结构化网格。通过将网格粗化、细化、吸附以及在自由面附近进行细化形成贴体网1格。其中，第 层边界层网格间距及边界层层数根据无量纲特征数Y 确定。分别对双艉客船带2舵和不带舵这 种状态下的阻力进行计算，其中82带舵的计算网格总数为 万，不带舵的网格总数76为 万。网格质量通过正交性来保证，即绝大部90°。其中，不带舵网格的最小分网格的正交性为≥22.5°，对于带舵的网格，由于舵尾缘尺正交性寸较小，故最小正交性≥8°，网格质量满足数值计3算要求。船体表面网格划分如图 所示。
4 计算结果分析
为方便对结果的处理，引入了无量纲系数：傅汝德数Fr、雷诺数Re、船舶总阻力系数 C 和摩擦T
阻力系数C ，分别定义如下： F_ITTC
式中：V 为船模速度；g 为重力加速度；ρ 和 μw m w 2 2由表 可以看出，不带舵和带舵 种状态下的阻力数值计算结果与试验结果均吻合良好，误差3%以内。在不考虑舵影响的计算中，也即裸船在体时，虽然网格质量和计算效率都很高，但计算结果与试验结果相比普遍偏小。相比于裸船体的计算，考虑带舵计算时网格数量会有所增加，网格划分难度也较大，但其计算结果与试验结果更为吻4 5合。图 和图 所示分别为裸船体和带舵船体在Fr=0.178下航行时的船艉流场对比。从图中可以看出：一方面，舵作为附体会增加船舶阻力；另一方面，也会改变船艉的流线分布。为更真实地反映水池试验状态，在后续的算例中，将均采用带舵船体进行计算分析。
4.2 对 的影响分析
Y+ CFD在进行船舶 绕流问题模拟时，通常将船［10］体表面视为无滑移壁面。根据流动特性的不同 ， 3壁面边界层可以分为 个区域：粘性底层、过渡层和惯性子层。边界层的位置通常用无量纲特征数+ Y 来表示。无量纲分析和试验工作表明，在边界+层特定的Y 范围内，具有特定的流动特点。对大 分别为水的密度与粘度； A 为船舶静水湿表面w
积；R为船模总阻力。
4.1 对带舵和不带舵的影响分析
双艉客船的水池试验考虑带舵的影响。在进CFD行 数值计算时，为便于比较分析舵对船艉流2动的影响，分别进行了裸船体和带舵船体 种状Fr=0.099~0.217态下的多航速计算，其中 。总阻力系数的数值计算结果及其与试验数据的比较如2表 所示。表中：C 为水池试验结果，考虑了T_EXP带舵的情况；C 为裸船体阻力计算结果； T_CAL1
C 为带舵船体阻力计算结果。定义误差： T_CAL2 ERR1 = (C -C T_EXP)/C T_CAL1 T_EXP ERR = (C - C T_EXP)/C 2 T_CAL2 T_EXP
多数高雷诺数流动来说，一般在粘性底层内使用壁面函数，以节约计算资源。在惯性子层区域，由于湍流应力近乎常数，故主流方向的流速分量呈对数分布。因此，为捕捉壁面附近流动的物理特性，可在物理特征变化较大的地方分布大量网格1点，尤其是第 层网格点应落在此区域，且必须布置足够的网格数以保证分辨所有的流动梯度。这1就要求正确布置边界层第 层网格点的位置以及+边界层层数。对于无量纲特征数Y ，可采用下式进行估算：
1式中：y 为第 层网格距离；V 为来流的参考wall ref速度；L 为参考长度；υ为粘性系数。ref =1.721 m/s针对本文所研究的双艉船舶，以V m Fr=0.237，Re=9.781×106为例，对应的 ，该航速下对= 4.247×10-3。研究此应的水池试验结果为C T_EXP 1航速下船舶第 层网格节点高度以及网格边界层层数对计算结果的影响，其中船舶初始网格设置及网格细化准则均与前述带舵航行时船舶阻力计算的设置相同。分别设置船舶边界层网格层数为6~10 1.2，逐步改变船舶第1层层，网格变化率为网格节点高度和边界层层数，以查看边界层网格3方案对计算结果的影响。计算结果如表 所示。 3 1由表 可以看出，随着边界层第 层网格点位+置的减小，阻力计算结果有所增加，而Y 则是逐1步减小。当第 层网格节点距离相同时，不同的边界层网格数对计算结果的影响也不一样，较多的边界层网格数能够更好地捕捉边界层内速度梯度的变化，计算结果将更加精确。对于船体近壁20~60，边+面流动计算，推荐的Y 取值范围为 界8~10层网格数约为 层。
4.3 不同航速下阻力的数值计算分析
双艉客船的设计航速范围较大，在进行船舶阻力模拟计算时，为进一步验证文中数值模拟方法的可靠性与可行性，对较高航速下的船舶总阻力 0.237，0.257 0.276也进行了计算，例如Fr取为 和+时。在高航速的数值计算中，考虑了最优Y 取值4及边界层数量。表 给出了带舵双艉客船在全航速范围内（Fr=0.099~0.276）不同航速下的总阻力数值计算结果，并与试验值进行了比较分析。 4由表 可以看出：在所计算的速度范围内，尽管航速变化较大，但和水池试验结果相比，数值计3%以内，说明本文中数值模拟算结果的误差均在精度较高，能满足工程需要。不同航速下，船舶各阻力占总阻力的百分比不一样。一般对水面排水型舰船而言，低速航行时，兴波较小，摩擦和粘压阻力占主要成分。而当航速增加时，船体兴波成分会逐渐增加，船舶兴波6阻力随之增大。图 所示为不同航速下船艏兴波高度的变化。从中可以看出，随着船舶航速的增加，船艏兴波高度也相应增大，因此在双艉客船高速航行时，应考虑船艏消波的必要性。
（f）Fr=0.276图6 不同航速下船艏兴波高度的变化Fig.6 Variation of wave height of the bow at different speeds 7图 所示为不同航速下的自由液面波形图。从中可以看出，不同航速下船舶自由液面的波形差别较大。这同时也表明文中所提计算方法是正确的，能准确反映物理现象，可为后续船型的改进提供参考。
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 boa［t J］. Chinese Journal of Ship Research，2016，11（4）：22-28 （in Chinese）. ［10］ 平. 学［M］. 2 版.刘应中，缪国 高等流体力 上海：上海交通大学出版社，2002.