JIANG Yichen;ZHAO Xiaojie;ZONG Zhi
【摘 要】为了研究带舭龙骨的船体横摇自由衰减运动,本文建立了基于OpenFOAM的粘性流场数学模型.模型采用六自由度运动方程与动网格技术计算船舶运动,应用阻尼消波技术消除远处的数值水池壁面产生的反射波,有效提高了数值预报的精度;通过对比模拟得到的衰减曲线和阻尼系数与试验所得的结果,分析了舭龙骨附近的涡量场,并探讨了舭龙骨减摇的机理.结果表明本文建立的数值模型可以准确地模拟船体横摇运动,成功预报了不同高度舭龙骨的减摇效果. 【期刊名称】《哈尔滨工程大学学报》
【年(卷),期】2019(040)001
【总页数】7页(P81-87)
【关键词】免摇启动器横摇衰减;舭龙骨;OpenFOAM;阻尼消波;动网格;六自由度;粘性流场;数值模拟
【作 者】JIANG Yichen;ZHAO Xiaojie;ZONG Zhi
光纤环网【作者单位】
锰氧化物【正文语种】中 文
【中图分类】U661.1
稳定性是船舶设计和航行过程中是重要因素，而横摇运动影响船舶稳定性，所以减少横摇运动是船舶安全航行的关键。在船体上安装舭龙骨是常用的减摇方法，它可以增加船体的横摇阻尼。Ikeda等[1]基于试验将横摇阻尼分成摩擦阻尼、兴波阻尼、旋涡阻尼、升力阻尼和舭龙骨阻尼。为了设计合理的舭龙骨，研究带有舭龙骨的船体横摇衰减模拟具有重要的意义。由于粘性对横摇运动具有重要的影响，横摇运动相比于其他自由度运动具有强烈的非线性，横摇运动的模拟比其他自由度运动更加困难。为了计算船舶的横摇运动，Ikeda等[1]基于试验结果得到半经验公式，但是精度不高。传统的势流理论方法不考虑粘性，仅仅可以计算出兴波阻尼，但是无法计算关于粘性的阻尼。如果船体加上舭龙骨，它会产生更强的粘性效应，传统的势流理论方法计算会更加困难。随着计算机技术的发展，计算流体
力学方法已经广泛地应用到船舶横摇的研究中。Wilson等[2]利用RANS方法研究DTMB5512船型在不同航速下的横摇衰减运动，并且展示了横摇运动产生的流场信息。Yang等[3]研究了船舶航行中的船舶横摇的问题，重点对强非线性进行了分析。罗天等[4]利用自主开发的naoe-FOAM-SJTU求解器研究了S60与DTMB5415船模的自由与强迫横摇运动，发现旋涡阻尼所占的比例最大。 Lavrov等[5]研究了二维剖面强迫横摇运动，重点分析了横摇中泄涡的现象。对于带有舭龙骨的横摇模拟，Xie等[6]基于FLUENT计算了带舭龙骨的柱体的横摇阻尼，并且得到了横摇过程中船体周围的压力分布图像，但没有对周围涡量场进行分析。蒋银等[7]结合重叠网格方法对三维带有舭龙骨船体的横摇运动进行了数值模拟研究，重点研究了阻尼系数的频率相关性。Irkal等[8]采用Flow3D软件研究了带有舭龙骨船体剖面的横摇阻尼运动，结合自己试验数据进行验证分析，但在数值模拟的后期出现了局部数值震荡。
本文基于OpenFOAM[9]中InterDyMFoam求解器，建立粘性流体数值模型，考虑计算域边界辐射波反射的影响，对带有舭龙骨的二维船体横截面的横摇运动进行数值模拟。利用流体体积分数方法(volume of fluid,VOF)方法捕捉自由面，利用六自由度运动模型和动网格技术计算船舶运动，由于本文只研究横摇运动，所以只考虑其中一个自由度。采用湍流模型
网联网来模拟横摇中的湍流，对不同高度舭龙骨下的横摇展开模拟，分析舭龙骨附近涡量场。
1 两相流数值模型
1.1 控制方程
采用两相流模型求解水和空气运动，利用VOF方法捕捉自由液面，假设水和空气不可压缩。连续性方程，动量方程，相方程表示为[10]：
连续性方程：
(1)
动量方程：
(2)
相方程：
(3)
式中：ρ表示流体密度；U表示速度矢量；Ug表示网格速度；μeff是考虑湍流的有效动力粘度系数，表达式为μeff=ρ(ν+νt)，其中ν为流体运动粘度系数，νt为湍流粘度系数；pd表示相对动压力，表达式为pd=p-ρgX，X为相对于参考高度的位置矢量；σ为流体表面张力系数，本文取0.728 kg/s2；κ表示自由表面曲率，它可以计算表面张力，计算公式为其中“·”表示求散度；γ为流体体积分数，由VOF相方程计算得出，如果γ=0，则表示网格内是空气，如果0<γ<1，则表示网格内是流体交界面，如果γ=1，则表示网格内是水；fs为阻尼项，用于消除横摇运动在远处产生的辐射波。
相方程中，Ur为用于界面压缩的速度场，只在自由面处起作用，压缩速度方向垂直于自由面，具体的表达式为其中表示自由面的法向方向，|U|为压缩速度的大小，表达式为表达式中φ表示速度通量；Ca为压缩因子，本文中取值1；Sf表示网格表面的面积矢量。
1.2 自由表面阻尼消波方法
本文采用二维剖面进行横摇数值模拟，计算模型示意图如图1。
图1 二维计算模型Fig.1 Two-dimensional calculation model
由于CFD模拟是在有限域中模拟，所以为了消除两边远处辐射波的反射所造成的数值误差，采用自由表面阻尼消波的方法，在动量方程里面添加阻尼项达到消波的目的：
fs=ρθU+ρθ′U
(4)
式中：θ与θ′分别用于右边和左边的消波系数。
(5)
式中：θ1为无因次阻尼系数，本文取θ1为用于确定消波长度，消波长度约为0.5倍波长，本文取1.0 m。
1.3 六自由度运动方程与动网格
为了计算横摇运动，通过求解六自由度运动方程来得到横摇的运动。
通过求解控制方程，可以得到船体周围流场的速度和压力，船体剖面受到的力F和力矩M可以通过边界表面积分得到：
(6)
式中：m为船体质量，r为力臂，τ为应力张量。
得到船体表面作用力之后，代入船体运动方程可得船体运动：
(7)
it运维系统详细设计式中：分别为3个坐标轴方向的加速度，分别为3个坐标轴方向的角加速度。X、Y、Z和K、M、N分别为力和力矩的3个分量。Ix、Iy、Iz为船体绕质心的惯性矩。
求得船体加速度后，对于时间进行积分可以得到时间段内船体的位移，进而得到船体下一个时刻的位置。时间积分步进方法是通过Newmark算法实现。
Newmark算法[11]的基本原理如下，下个时间步的加速度和速度为
(8)
式中δ和β是与精度和稳定性有关的参数，本文分别取0.5与0.25。
为了实现船体边界的运动，计算过程中采用动网格技术，本文采用的是类似于“随体网格”的方法，如图2所示。变形网格层用来确定网格变形的区域，在船体附近的那层网格同船体一起运动，运动方向和大小和船体一致，而远离船体的最外层网格固定不动，之间的网格按照线性插值决定网格变形的大小。
图2 动网格示意图Fig.2 Sketch of dynamic grids
值得注意的是要选取适当的距离避免求解过程网格过于形变造成计算发散，近距离值选用距离壁面相对近的位置，大概4～5层近壁面网格左右，而远距离值可以适当选用较大的距离。本文模拟计算选用近距离值为0.005 m，而远距离值选用0.4 m。
1.4 湍流模型
船体横摇过程中流场处于湍流状态，本文采用二方程的湍流模型进行处理，采用低雷诺数湍流模型k-ω模型进行计算。
k-ω双方程模型引入两个附加的传输方程，一个为湍流动能k，另一个为湍流耗散率ω。该模型由Wilcox等[12]提出，其中k与ω分别为:
(9)
(10)
湍流动能和湍流耗散率的约束方程为：
(11)
Pk是湍动能生成项，定义为：
(12)
式中μt为湍流粘度，计算公式为μt=ρk/ω。该模型的常数α=0.52,β=0.072,β*=0.09,σk=0.5,σω=0.5。
麻元友在运用湍流模型的时候要处理近壁面内的流动，近壁面处边界层内的流动可以分为粘性层，缓冲层和对数律层。网格高度的选取根据离壁面的无量纲距离y+确定。粘性层满足y+≤5。y+的定义为
(13)
式中：y为第一层网格节点高度;uτ为摩擦速度；ν为流体运动粘度系数。
本文采用壁面函数法，由于模型尺度较小，将第一层网格放置在粘性层，靠近壁面最小的网格在1 mm，这样可以应用于低雷诺数湍流模型。
1.5 边界条件和初始条件
由于本文采用二维边界，所以前后两个边界定义为empty边界，表示模型是二维。除此之外，其他的边界有左右两边的出流边界outflow1与outflow2，上面的大气边界atmosphere，下面的底部边界bottom，还有二维船体剖面边界hull。
关于数值边界条件，在OpenFOAM中是通过设置各个物理场的边界来实现。
对于速度场，在两边出流边界和大气边界定义为自由出流边界条件，在底部速度为0，在船体表面为物面速度。对于压力场，在出流和底部边界定义为零梯度边界，在大气边界处压力为一个大气压，在船体表面为针对压力的第二边界条件。因为采用VOF方法捕捉自由面，所以要设置相场边界条件，在边界处设置为零梯度边界。另外关于湍流场的边界，除了在壁面处如船体和底部处采用壁面函数边界，其他边界设置为零梯度边界。
对于初始条件，本文计算自由横摇运动，初始时刻将船体绕质心旋转一定角度，船体初始速度为0，通过setField工具设置相场来定义自由面，水的体积分数为1，空气的体积分数为0。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议