2. 江苏现代造船技术有限公司, 江苏 镇江 212003
2. Jiangsu Modern Shipbuilding Technology, Ltd., Zhenjiang 212003, China
在传统的船舶设计中,操纵性相较于船舶快速性和耐波性是次要考虑因素,然而船舶操纵性是关系到航行安全的一个重要因素,随着 IMO(International Maritime Organization)对操纵性标准的提高,船舶操纵性受到越来越多的关注[1]。自航模或约束模试验是预报操纵性最可靠的方法,然而船模试验需要特定的实验设备,且船模试验存在实验周期较长、实验费用高昂及尺度效应无法消除的缺点。
基于计算流体动力学(Computational Fluid Dynamic, CFD)技术的数值模拟方法是近年来出现的最先进的船舶操纵预报方法,该方法通过离散求解非定常操纵运动方程,得到相关船舶运动参数和运动轨迹,从而实现用纯数值计算的方法对船舶操纵性进行理论预报。与船模试验相比,数值仿真计算具有计算费用少、计算周期短,而且能够获得局部流动状态信息的优点。随着计算机技术的发展,CFD 技术在操纵性领域从最初的定常运动仿真发展至现在用于非定常运动仿真[2-4]。本文以 DTMB5415 为仿真对象,以横荡幅值为 0.064m 的试验数据为基准[5],设置不同边界条件进行数值仿真计算,对计算结果进行分析得到最佳边界条件设置。然后基于小振幅假设对低频率纯横荡及纯首摇运动进行数值仿真计算,求解船舶操纵相关的粘性流动,得到满意流体水动力以及水动力矩。
1 运动数学模型1)纯横荡运动
船模匀速前进的同时叠加一个横向正弦规律变化的位移[6-7]。
$\left\{ \begin{array}{*{35}{l}} y=-y\sin \left( \omega t \right), \\ v=-{{v}_{\max }}\cos \left( \omega t \right),{{v}_{\max }}={{y}_{\max }}\omega \\ \dot{v}={{{\dot{v}}}_{\max }}\sin \left( \omega t \right),\dot{v}={{y}_{\max }}{{\omega }^{2}}\sin \left( \omega t \right) \\ \psi =0 \\ \end{array} \right.$ | (1) |
式中:ymax为横荡运动幅值;ω为横荡运动频率;ν为横荡运动速度;
由于本文研究的是小振幅横荡运动,线性假设下船模进行纯横荡运动时所受的水动力和水动力力矩可表示为:
$\left\{ \begin{array}{*{35}{l}} Y={{Y}_{{\dot{\nu }}}}\dot{\nu }+{{Y}_{\nu }}\nu , \\ N={{N}_{{\dot{\nu }}}}\dot{\nu }+{{N}_{\nu }}\nu \\ \end{array} \right.$ | (2) |
将式(1)和式(2)合并有:
$\left\{ \begin{array}{*{35}{l}} Y=-\alpha {{\omega }^{2}}{{Y}_{{\dot{\nu }}}}\sin \omega t+\alpha \omega {{Y}_{\nu }}\cos \omega t\text{,} \\ N=-\alpha {{\omega }^{2}}{{N}_{{\dot{\nu }}}}\sin \omega t+\alpha \omega {{N}_{\nu }}\cos \omega t \\ \end{array} \right.。$ | (3) |
由式(3)可知,只要测得船舶横荡运动过程中所受的侧向力及转首力矩就可以获得
2)纯首摇运动
船模保持匀速前进的同时首向角进行正弦规律变化[6-7]。船模运动轨迹时刻与纵轴保持相切。运动规律为:
$\left\{ \begin{array}{*{35}{l}} \psi =-{{\psi }_{\max }}\cos \left( \omega t \right)\text{,} \\ r={{r}_{\max }}\sin \left( \omega t \right),{{r}_{\max }}={{\psi }_{\max }}\omega \text{,} \\ \dot{r}={{{\dot{r}}}_{\max }}\cos \left( \omega t \right),{{{\dot{r}}}_{\max }}={{\psi }_{\max }}\omega \\ \end{array} \right. 。$ | (4) |
式中:ψmax为首向角幅值;ω 为首摇运动频率;r 为首摇运动角速度;
由于本文研究的是小振幅首摇运动,线性假设下船模进行纯首摇运动时所受的水动力和水动力力矩可表示为:
$\left\{ \begin{array}{*{35}{l}} Y={{Y}_{{\dot{r}}}}\dot{r}+{{Y}_{r}}r, \\ [5pt]N={{N}_{{\dot{r}}}}\dot{r}+{{N}_{r}}r, \\ \end{array} \right.$ | (5) |
合并式(4)和式(5)有:
$\left\{ \begin{array}{*{35}{l}} Y=-a{{\omega }^{2}}{{Y}_{{\dot{r}}}}\sin \omega t+a\omega {{Y}_{r}}\cos \omega t\text{,} \\ [5pt]N=-a{{\omega }^{2}}{{N}_{{\dot{r}}}}\sin \omega t+a\omega {{N}_{r}}\cos \omega t \\ \end{array} \right.。$ | (6) |
由式(6)可知,只要测得船舶首摇运动过程中所受的侧向力及转首力矩就可以获得
1)网格划分
本文以 DTMB5415 为对象进行数值仿真计算[8],缩尺比为 1∶46.588。几何模型如图 1 所示。
在数值仿真计算过程中,计算域的划分对计算结果的精确度与计算时长都密切相关。计算域过大,会使得网格数目增多,造成计算资源浪费;计算域划分过小,会产生浅水效应和岸壁效应,使数值模拟与真实的运动状态不同[9]。文中计算域划分如图 2 所示。
网格质量是影响数值仿真计算的可行性、收敛性和计算精度的主要因素。在进行网格划分时考虑到船体表面曲率变化较大,需进行局部加密,同时对自由液面进行加密,以便捕捉液面的变化。网格划分完成以后需对网格质量进行检查,所划分的网格不能包含负体积及凹体积,负体积是无效的计算单元,凹体积将影响计算稳定性,同时推荐网格正交性大于 5。
2)CFD 边界条件设置
对纯横荡运动进行 3 种不同幅值的低频运动仿真;纯首摇运动进行 3 种不同首向角的低频运动仿真。计算工况见表 1。
图 3 中 Slip、Wall 代表数值仿真计算过程中甲板设置为滑移壁面条件和壁面函数条件,Exp 代表试验值。滑移壁面条件应用于壁面切向速度为 0,壁面函数法主要适用于低雷诺数层流流动。图 3为 2 种边界条件下的数值仿真计算结果与实验数据拟合图形。
甲板边界设置为滑移壁面时,数值仿真结果与试验数据之间相对误差均值为 0.25%,甲板采用壁面函数时,数值仿真结果与试验数据之间相对误差均值为 0.45%。从计算结果来看,考虑空气粘性影响后,数值计算结果较之前计算值略有变大,增大幅度不超过 0.5%。但考虑空气粘性影响后,网格需在滑移壁面网格基础上进行加密处理,从节省计算资源及提高效率的角度出发,本文后续仿真计算,忽略甲板空气粘性的影响,甲板的边界条件采用滑移边界条件。
图 4 中 bottom_pres 及 bottom_far 分别表示流域下壁面边界条件为指定压力边界(Prescribed Pressure)和远场边界(Far Field)条件,Exp 代表试验值。远场边界条件是通过速度分量的大小初始化边界上的速度,所给定的速度方向取决于具体的流动状态。指定压力边界条件属于 Dirichlet 条件,在计算初始化过程中指定边界压力。通过图 4 可以发现:底部壁面设置为指定压力时,数值仿真计算结果与试验值吻合度较好,二者之间的平均相对误差为 8.75%;底部壁面设置为速度远场时的数值仿真计算结果已经无法用于实际工程。
综上,对不同边界条件设置得到的计算结果进行分析可知,空气密度较小,甲板设置为滑移壁面或壁面函数时,2 种边界条件下的数值仿真计算结果与试验数据之间的相对误差均较小。为充分利用计算资源,本文进行后续数值仿真计算过程中甲板采用滑移壁面边界;流域底部壁面设置为远场边界条件时的数值仿真计算结果已无法用于实际工程,而流域底部壁面设置为指定压力时,数值仿真计算结果与试验数据之间相对误差较小。因此,在后续数值仿真计算过程中为了得到理想结果,甲板边界选择滑移壁面边界,流域底部壁面需选取指定压力边界条件。
3)纯横荡运动
纯横荡运动是通过船模恒定速度(U = 1.531 m/s)前进,同时沿横向叠加正弦规律变化的位移实现。为了便于将数值仿真计算结果与试验数据进行对比,数值仿真计算过程中设置水温与试验水温保持一致,设置为 13.5 ℃,流体特性参照 ITTC 推荐值进行设置,横荡运动幅值分别为 0.064 m,0.128 m,0.317 m。
图 5 中 Sim 0.064 等分别代表不同横荡幅值下一个周期内数值仿真计算结果,Exp0.064 等分别代表不同横荡幅值下一个周期内试验值。船模试验与数值仿真计算最大、最小值及相对误差见表 2。采用本文介绍的 CFD 设置及网格划分,计算出的最大误差为-20.9%,最小误差为-0.5%,本文计算结果相比于参考文献,纯横荡的侧向力及转首力矩的相对误差有一定的改善。考虑到船体形状和运动的复杂性,这样的误差在可以接受的范围内。而误差的来源可能是:数值仿真计算中忽略了运动过程中产生的兴波;实验过程中没有限制纵倾和升沉运动,然而数值仿真过程中没有考虑这 2 种运动。
4)纯首摇运动
纯首摇运动是通过船模以 1.531 m/s 的速度匀速前进,同时船模绕铅直轴进行正弦规律变化的运动实现的。运动过程中船模合速度与船模的中纵剖面始终相切。为了便于将数值仿真计算结果与试验数据进行对比,设置水温为 13.5 ℃,流体特性参照 ITTC 推荐值设置,首摇运动首向角幅值分别为1.7°,5.1°,10.2°
图 6 中 Sim 1.7 等分别代表不同横荡幅值下一个周期内数值仿真计算结果,Exp1.7 等分别代表不同首向角下一个周期内试验数据。船模试验与数值仿真计算最大值、最小值及相对误差见表 3。数值仿真计算结果与船模试验数据最大相对误差为 27.5%,最小相对误差-0.88%。考虑到船体形状和运动的复杂性,这样的误差在可以接受的范围内。而误差的来源可能是:数值仿真计算中忽略了运动过程中产生的兴波;实验过程中没有限制纵倾和升沉运动,然而数值仿真过程中没有考虑这 2 种运动。
本文利用商用软件 Fine/Marine 以 DTMB5415 为模型对纯横荡运动的一种工况试验数据为基准,采用不同边界条件进行数值仿真计算,分析边界条件对数值仿真计算结果的影响,从中找出适宜的边界条件进行纯横荡运动及纯首摇运动仿真。将纯横荡运动及纯首摇运动数值仿真计算结果与实验数值结果进行对比,数值仿真计算与船模试验结果曲线吻合较好。考虑船型的复杂及局部信息处理的不全面性以及数值仿真计算过程中忽略自由兴波的影响,计算结果在可接受范围内。本文数值仿真计算计算得到相对误差总体相对参考文献相对误差有一定改善,从而说明本文在网格划分及边界条件设定符合本文数值仿真计算要求。
[1] | YOON H, SIMONSEN C D, BENEDETTI L, et al. Benchmark CFD validation data for surface combatant 5415 in PMM maneuvers-part I: force/moment/motion measurements[J]. Ocean Engineering , 2015, 109 :705–734. DOI:10.1016/j.oceaneng.2015.04.087 |
[2] |
田喜民, 邹早建, 王化明. KVLCC2船模斜航运动粘性流场及水动力数值计算[J]. 船舶力学 , 2010, 14 (8) :834–840.
TIAN Xi-min, ZOU Zao-jian, WANG Hua-ming. Computation of the viscous flow and hydrodynamic forces on a KVLCC2 model in oblique motion[J]. Journal of Ship Mechanics , 2010, 14 (8) :834–840. |
[3] | SIMONSEN C D, STERN F, AGDRUP K. CFD with PMM test validation for manoeuvring VLCC2 tanker in deep and shallow water[C]//Proceedings of International Conference on Marine Simulation and Ship Maneuvering (MARSIM). Terschelling, the Netherlands: MARSIM, 2006. |
[4] | LIU H, MA N, GU X. Maneuvering prediction of a VLCC model based on CFD simulation for PMM tests by using a circulating water channel[C]//ASME 201534th international conference on Ocean, Offshore and Arctic Engineering. St. John's, Newfoundland, Canada: The American Society of Mechanical Engineers, 2015. http://cn.bing.com/academic/profile?id=2333204890&encoded=0&v=paper_preview&mkt=zh-cn |
[5] |
刘山. 基于CFD技术数值模拟平面运动机构试验[D]. 武汉: 武汉理工大学, 2012.
LIU Shan. Numerical simulation of planar motion mechanism test based on CFD technology[D]. Wuhan: Wuhan University of Technology, 2012. http://www.oalib.com/references/14732837 |
[6] |
吴秀恒, 刘祖源, 施生达, 等.
船舶操纵性[M]. 北京: 国防工业出版社, 2005 .
WU Xiu-heng, LIU Zu-yuan, SHI Sheng-da, et al. Ship maneuverability[M]. Beijing: National Defence Industry Press, 2005 . |
[7] |
范尚雍.
船舶操纵性[M]. 北京: 国防工业出版社, 1988 .
FAN Shang-yong. Ship Maneuverability[M]. Beijing: National Defence Industry Press, 1988 . |
[8] | http://www.simman2008.dk |
[9] |
郭然, 贾力平, 樊小莉, 等.
NUMECA系列教程[M]. 北京: 机械工业出版社, 2013 .
GUO Ran, JIA Li-ping, FAN Xiao-li, et al. NUMECA series tutorial[M]. Beijing: Machinery Industry Press, 2013 . |