基于cfd方法的长峰不规则波数值模拟

基于cfd方法的长峰不规则波数值模拟

船舶的垂直移动运动主要包括船舶的纵向移动和悬挂。严重的纵摇和垂荡将引起甲板上浪、砰击和失速等一系列后果,对船舶在海上的航行安全以及舰载武器的使用影响极大。因此,本文选取船舶的纵向运动作为研究对象,基于CFD方法对船模在长峰不规则波中顶浪运动进行数值模拟,具有迫切的现实意义。1浮体运动数值模拟设施要实现对船模在长峰不规则波中顶浪纵向运动的数值模拟,首要前提就是构建精度满足耐波性计算要求的长峰不规则波数值波浪水池。所谓的数值波浪水池,是指对非线性波浪水动力以及浮体运动数值模拟设施的统称。它能够通过实验观测为各种波浪理论的研究奠定坚实的物理基础,还能为海洋、船舶等工程设计提供可靠而高效的试验数据。本文基于粘性流理论构建长峰波数值波浪水池,采用有限体积法(FVM)对RANS方程和连续控制方程进行离散求解,利用Fluent软件的二次开发功能UDF,编写边界条件完成数值造波和阻尼消波功能。选用ITTC海浪谱作为目标海浪谱,构建长峰波数值波浪水池,并对其计算精度进行误差计算。1.1波水池价值1.1.1双转台波的入射速度数值造波是指用数值方法模拟波浪的生成过程,为数值模型实验提供各种形式的波浪环境条件。本文采用的数值造波方法是边界条件造波法,即在数值波浪水池入口边界设置3个方向的速度。X方向沿水池向下游为正,Z方向向上为正,Y轴与X轴和Z轴符合右手法则。3个方向的入射速度分别为:⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪U=∑i=1NωiAiekizcos(kix−ωit+εi),V=0,W=∑i=1NωiAiekizsin(kix−ωit+εi)。(1)[BFQ]{U=∑i=1ΝωiAiekizcos(kix-ωit+εi),V=0,W=∑i=1ΝωiAiekizsin(kix-ωit+εi)。(1)[BFQ]式中:η为波动水面相对于静止水面的瞬时高度;Ai,ki,ωi和εi分别为第i个组成波的波幅、波数、圆频率和初始相位,εi是在(0,2π)范围内的随机相位;X轴为波浪传播方向;U,V,W分别为波浪X轴、Y轴和Z轴的速度分量。1.1.2阻尼消波器计算当波浪到达水池末端开边界处,会引发水波的二次反射,反射波与入射波叠加,会导致波场的失真。因此,消除反射波的影响也是数值波浪水池的一项重要技术。数值波浪水池常用的消波技术主要有辐射边界条件消波、主动消波器消波、阻尼消波3种。本文选取的是阻尼消波。阻尼消波法是指在流场中添加人工粘性,因其对来波的频率和波长不敏感,可以有效地消除各种频率和波长的来波,因而被广泛采用。在计算域出口边界前设置1～2倍波长的阻尼消波段,利用Fluent中的UDF宏DEFINE_SOURCE(mom_source,cell,thread,dS,eqn)编程实现消波。在阻尼消波段内,动量方程写为:∂u∂t+u∂u∂x+v∂u∂y=−1ρ∂p∂x+v(∂2u∂x2+∂2u∂y2)−μ(x)u,(2)∂u∂t+u∂u∂x+v∂u∂y=-1ρ∂p∂x+v(∂2u∂x2+∂2u∂y2)-μ(x)u,(2)∂v∂t+u∂v∂x+v∂v∂y=g−1ρ∂p∂y+v(∂2v∂x2+∂2v∂y2)−μ(x)v∂v∂t+u∂v∂x+v∂v∂y=g-1ρ∂p∂y+v(∂2v∂x2+∂2v∂y2)-μ(x)v。(3)其中μ(x)为在阻尼段起点为0的单调递增函数,可以取为线性递增、指数递增等形式。取μ=105(x-Xmin_D)/(Xmax_D-Xmin_D)。(4)式中:Xmin_D和Xmax_D分别为消波区的最小、最大X坐标。1.2长峰不规则波值波池建设1.2.1垂直网格化社区的网格结构三维数值水槽的网格划分如图1所示,本文构建的长峰不规则波数值水池长18m,其中12～18m为消波区,宽3m,深2.5m,自由面以上1.5m,整个水槽的网格数为413820。垂直自由面方向的网格尺寸取为有义波高1/5。造波区沿X轴正方向网格尺寸与垂直于自由面Z方向的最小网格尺寸相同,从自由面到水槽顶部网格按1∶1.1的比例等比分布;从自由面到水池底部网格按1∶1.05的比列等比分布。网格基本上是离自由面越远尺寸越大。对于消波区的网格划分,垂直方向划分与造波区的一致,水平方向网格以造波区网格尺度为基准向右边界逐渐扩大。1.2.2计算示例描述选取ITTC海浪谱作为目标靶谱,对3种海况下的长峰不规则波进行数值模拟。波浪的目标参数见表1。1.2.3长峰不法波数值波浪水池波面时历曲线长峰不规则波数值模拟瞬时波面场和局部速度矢量如图2和图3所示。图4和图5是有义波高分别取H1/3=0.08m,H1/3=0.125m时,长峰不规则波数值波浪水池X=1m,X=6m,X=11m处的波面时历曲线对比图。从上述数值造波水池不同位置的波面时历曲线对比可看出,长峰不规则波在沿X轴正方向向下游传播过程中,波能有一定程度的衰减,这是由于波浪在传播过程中,高频子波衰减导致的。1.2.4目标谱的建立对上面监测得到长峰不规则波波面时历进行谱分析与目标谱对比如图6所示。分别从谱面积m0、谱峰频率ωp和有义波高H1/33个方面对上述数值模拟海浪谱进行误差分析,误差计算结果如表2和表3所示。2船舶模型中的行程曲线模型对长峰不规则波中的波浪垂直运动进行了数值模拟2.1海上模型几何模型选取具有球鼻首和方位的DTMB5512船模作为研究对象,该船模是ITTC(国际船模试验水池会议)推荐的瘦削型标准船模DTMB5415的全相似几何模型,以美国海军DDG-51型驱逐舰为模板,船型数据详实,而且与我海军舰船船型相似(见图7),适合军舰参考。表4为该船模及其所对应的实船尺度的主要数据。2.2数值模拟方案2.2.1垂荡和纵摇2个自由度的运动本文对DTMB5512船模在长峰不规则波中顶浪纵向运动数值模拟试验中,只考虑垂荡和纵摇2个自由度的运动。船模CFD耐波性数值模拟试验对计算资源的要求较高,考虑到本文研究所使用的计算机配置的实际情况,以及试验水池网格划分所带来的计算效率等问题,试验计划如表5所示。2.2.2耐波性数值波浪水池网格划分参照《水面船模耐波性实验规程》,将计算域设置成长方体形状,如图8所示。船模与计算域各边界的位置关系如下:入口距船首1倍船长,出口距船尾2倍船长,顶部边界距水线0.5倍船长,底部边界距水线1倍船长,左、右边界距船中纵剖面0.5倍船长。耐波性数值波浪水池分成5个区域进行网格划分,即近船体区域、近流场区域、自由面区域、上下远流场区域及消波区域、各区之间互不重叠,且连接处选择connected方式,如图9所示。近船体区网格划分如图10所示。2.3udf模拟水池试验过程根据船舶六自由度运动的控制方程(5),当波浪作用于船体时,其运动的速度、角速度以及位置、姿态等可以通过控制方程求解、积分得到。对于流浮耦合运动,波浪作用在船体上的力和力矩使船体产生运动,同时船体的运动又对其周围流场产生影响。因此本文在数值模拟中,分段计算流体与船体运动的耦合,步骤如下:1)将船模按初始浮态固定,原点与重心重合;2)对流场进行初始化,设定初始航速后造波;3)以时间步长Δt=0.001s步进;4)通过当前流场变量迭代求解流场的速度矢量;5)通过压力—速度耦合算法获得压力场;6)求解体积分数方程重构自由面;7)读取船体受力F和力矩M,根据控制方程mz⋅⋅=F,Iϕϕϕ⋅⋅=Mmz⋅⋅=F,Ιϕϕϕ⋅⋅=Μ求解船体运动要素,经积分得到船体x方向的阻力、垂向受到的垂向力合纵摇角速度及纵摇幅值;8)更新船体位置和浮态;9)返回第3步,求解改变浮态后各量,并根据计算再次改变浮态,按此迭代求解,直到方程组的残差小于设定值或迭代次数达到设定值;10)返回第2步并重复以下步骤,直到设定的时间步数计算完毕。mz⋅⋅=−z⋅⋅∫Lm′dx−z˙∫L(N′−vdm′dx)dx−2ρgz∫Lywdx+ϕ⋅⋅∫Lm′xdx+ϕ˙∫L(N′x−2m′v−xvdm′dx)dx+ϕ∫L(2ρgywx−N′v+v2dm′dx)dx+ζA[∫L2ρgywe−kTcossinkxdxmω∫L(N′−vdm′dx)ekTsincoskxdx−ω2∫Lm′ekTcossinkxdx]cossinωet。(5)[BFQ]mz⋅⋅=-z⋅⋅∫Lm′dx-z˙∫L(Ν′-vdm′dx)dx-2ρgz∫Lywdx+ϕ⋅⋅∫Lm′xdx+ϕ˙∫L(Ν′x-2m′v-xvdm′dx)dx+ϕ∫L(2ρgywx-Ν′v+v2dm′dx)dx+ζA[∫L2ρgywe-kΤsincoskxdxmω∫L(Ν′-vdm′dx)ekΤcossinkxdx-ω2∫Lm′ekΤsincoskxdx]sincosωet。(5)[BFQ]通过以上迭代、循环,可实现流体与船体运动的耦合。步骤2按照1.2节构建的数值波浪水池进行造波和消波。试验过程中可以将事先保存好的稳定流场导入耐波性数值模拟水池,进行数值模拟,这样可以提高计算效率,缩短试验时间。步骤4～6,按前文设置的数值方法进行计算。步骤7则通过UDF编程实现。步骤9残差标准可取软件默认设定值,迭代次数上限为40步。2.4纵向运动模拟系统主要解决了以下方面的问题对于DTMB5512船模由于带球鼻首、具有方尾,线性复杂,所以实现其耐波性的数值模拟,要比一般商用船模困难得多。可采用以下方法对实验进行改进:1)改进船模贴体网格质量、数量及分区网格的匹配;2)采用递增方法实现船模的纵向运动,即先对船模进行单自由度纵摇,稳定之后再增加垂荡的数值模拟;3)逐渐增加船模质量(将船模质量降低到原船模的一半进行数值模拟,当残差稳定后再逐渐增大质量直至原值);4)减小欠松弛因子,控制单元网格内速度的变化量。通过上述方法,解决了DTMB5512船模复杂的几何船形与波浪作用的数值问题。DTMB5512型船模在长峰不规则波中纵向运动数值模拟的压力场和速度场,如图11和图12所示。DTMB5512型船模在长峰不规则波中纵摇及垂荡的时历曲线如图13和图14所示。3顶浪运动时历曲线在计算精度验证中的应用由于本文对船模在长峰不规则波中顶浪纵向运动的数值模拟研究目前在国内外还处于起步阶段,未找到具体的水池实验数据。考虑到基于势流理论舰船六自由度计算,在理论上是成熟的,并在工程应用上取得了很多成果,得到了水动力学界认可。一般来说,SCFD方法由于考虑到流体粘性,计算精度应略高于势流理论计算结果,但没有本质上的差异。至于非线性摇荡则需要另作考虑,本课题局限于线性摇荡,所以用势流理论计算结果进行验证是可行的。由于势流理论切片法计算过程相对繁琐,可参考文献,本文直接给出DTMB5512型船模在算例波浪环境下纵摇和垂荡的响应方差。将数值模拟的DTMB5512型船模摇荡运动时历曲线运用线性谱分析方法进行分析,获得纵摇及垂荡的摇荡谱,如图15和图16所示。对上述摇荡谱进行积分即可得到本文数值模拟DTMB5512船模顶浪运动纵摇和垂荡的响应方差。对DTMB5512型船模在长峰不规则波中顶浪纵向运动的数值结果与势流理论计算结果进行对比得出相对误差如表6所示。4船模、船模间的布置要求本文运用SCFD方法对舰船在长峰不规则波中顶浪运动进行数值模拟,得到如下可供参考的经验:1)船舶摇荡、阻力、操纵、推进等课题的数值研究,技术细节上的一个主要不同,体现为网格的布设上。2)船舶耐波性研究的网格,必须将造波和船体网格、动网格三者进行很好的协调、匹配,做到三者的有机结合,否则会出现计算发散及非物理现象等不合理现象的发生。3)船舶多自由度摇荡试验目前的困难主要集中在计算资源上,由于其计算量巨大,PC机以及低端的工作站、服务器已经不能满足其正常情况下的计算需要。如果计算资源等硬件设施有限,则需要在离散方法、格式,控制方程、调节参数、UDF开发等“软件”上下功夫。4)在三维空间内完成船模非规则波中的摇荡试验,其流场及自由面要比船模在规则波中的复杂得多,因此在船模贴体网格的布设上,要求网格质量非常高,同时在interface交界面处,左右网格要尺寸一致,且在船模外表面尽可能多地使用结构性网格,特别是在阻力计算上,还要尽可能多地布设边界层,以提高计算的精度。5)耐波性流场的高度复杂性还表现在输入、输出及船模的响应及其变化率上。具体表现为1个网格上,至少要输入3个方向的线速度u,v,ω,压力P,参数k,ε,流体体积分数