肥大型船伴流场网格划分方法的比较与分析.docxVIP

肥大型船伴流场网格划分方法的比较与分析 在使用fd方法求解流场时，首先需要对计算区域进行分散，即网格的划分。网格是研究对象的几何表达形式，也是数值模拟和分析的载体。此外，网格分离的质量直接影响到计算的可行性、容忍度、计算的精度、计算效率，甚至计算的成败。如何生成高质量的计算网格已成为当前ifd应用的一个重要课题。在目前发表的论文中，网格分离有几种方法，主要包括单个网格、多个网格和非海战网格，它们都有各自的优缺点和相应的应用领域。然而，到目前为止，还没有关于网格分割的文章具体讨论了不同网格划分方法对计算结果的影响，以及适应特定船型的成熟网格划分方法。 该文以某肥大型散货船为例,较为详细地探讨了适合肥大型船的网格划分方法,同时研究了不同的网格划分方式对计算结果所造成的不同影响.为此,采用了一种适应性更强的多重网格划分方式. 1 船受参数设置 以某散货船为研究对象.实船船长283.4 m,船宽45 m,吃水16.5 m,方型系数0.82,模型的缩尺比为1/38.18.船体艉部外形如图1所示. 由于船体关于中线面对称,并且数值计算的主要目的是获得船后螺旋桨桨盘面处的伴流分布,因此利用粘性流场的对称性,所建计算模型采用叠模法,只考虑半个船体的流场.该计算域从船艏向前延伸1个船长,从船尾向后延伸4个船长,沿宽度方向延伸5个半宽,深度方向延伸7个吃水.计算模型在前处理软件GAMBIT中生成. 2 网格的分段方法和比较 2.1 网格划分方式 划分网格首先要考虑两个问题,也就是网格尺度和网格数量问题.对于网格数量问题,应该以能否得到满足精度要求的计算结果为准.Tzabiras指出:在网格尺度达到某一下限值之后,解对网格精度的变化不再敏感,网格精度提高50%,其计算结果变化很小. 为控制网格数量,同时使计算域中流体流动受船体几何外形影响较大的区域有足够分辨率的网格,需要将计算域分块.在实际计算中,经多次尝试发现,流体流动受船体几何外形影响较大的区域分布在离船体表面大概半个船宽的距离内,如图2所示.据此可将计算域分为船体近流场控制域和远流场控制域,其中近流场控制域以设计水线面、中线面和中站面为基准,向外延伸半个船宽的距离,生成一个长方体,剩余部分为远流场控制域,其分块方式如图3所示. 在近流场控制域,最小和最大网格尺度可分别取螺旋桨直径的0.02倍和0.25倍.远流场控制域网格则以近流场控制域外表面网格为基础,并以1.2倍的增长比例外推生成结构化网格,整个流场网格划分如图4所示. 2.2 网格划分及计算结果 船体表面网格的尺度大小对网格总数起决定作用.较小尺度的网格能准确反映船体表面曲率的变化,精确地捕捉到由于曲率变化而产生的流动细节,当然网格数也成倍增加.较大尺度的网格虽然可以极大地降低网格数,但是在船体表面曲率较大的地方却不能反映出相应的流动细节.而数值模拟的目的就是为了准确地捕捉到一些具体的流动细节,因此有必要探讨舭部不同的网格划分方式对计算结果的影响. 肥大型船的特点是船体艏尾曲率大,船舯附近曲率小,近似为圆角矩形,为了比较船体表面网格尤其是舭部网格划分方法的不同对尾流场的影响,这里给出了2种典型的划分方式.这2种方式都以螺旋桨轴包套为基准向船艏方向进行网格加密,轴包套处网格尺度取0.2 m,船尾和船艏曲率较大的地方网格尺度取0.3 m,唯一不同之处是对船体舭部网格的处理.第1种方式忽略舭部曲率大的特点,将其和船舯一样处理,网格尺度均取0.8 m,划分的网格及相应的计算结果如图5和图6所示;第2种方式将舭部网格尺度取为0.3 m,且随着曲率变小,网格尺度逐渐增加到0.8 m,划分的网格及相应的计算结果如图7和图8所示.从图中可以看出,船体表面网格的疏密对计算结果有着极大的影响. 2.3 网格加密方式的确定 显然,对于肥大型船,舭涡的预报相当重要,而反映到桨盘面处,则是该处的伴流分数,因此,如果要捕捉到该处的流动细节,桨盘面处必须要有足够的网格分辨率,使用ICEM CFD生成网格时,可以通过局部加密网格来实现.加密网格必然会使整个流场网格数量增加,如何既能保证局部网格的分辨率,使数值模拟成功进行,又能使网格增加的数量控制在计算机可接受的范围内,成为研究人员必须解决的一个问题.该文对此做了初步研究,比较了2种典型的加密方式. 第1种是以桨盘面为基准,沿径向1.5倍螺旋桨直径,沿轴向1倍螺旋桨直径形成一个椭球型的加密区域,其网格形式及相应的计算结果如图9和图10所示.第2种方式是以桨盘面为基准,沿径向1.5倍螺旋桨直径,沿轴向向下游延伸1倍螺旋桨直径、向上游延伸4倍螺旋桨直径,形成一个圆柱型的加密区域,其最终的网格形式及相应的计算结果如图11、12所示.比较两者的计算结果,可以看到,第2种网格加密方式成功地捕捉到了肥大型船桨盘面处所特有伴流特征