sub-sshr-rsm对水沙运动的响应.docxVIP

sub-sshr-rsm对水沙运动的响应 1 相关研究问题的解决 自1980年代中期以来，国际水动力学研究院（iii）对潜水员的水流进行了广泛的数值模拟研究。Huang(1992)对于SUBOFF艇型进行了系统的风洞尾流场试验,采用BL模型、修正的BL模型与k-ε模型详细计算了桨盘面内的轴向无量纲速度。Atkins(1996)对装有围壳与四个尾翼的SUBOFF潜艇模型的流场进行了数值计算。Bull(1996)通过变化网格形式、网格数量、湍流模型、流动求解器,对于SUBOFF AFF-1(裸艇体)与AFF-8(全附体)的尾流场的数值模拟进行了系统的验证。Yang(2003)基于非结构化网格数值模拟了SUBOFF主体带指挥台围壳后的表面压力和摩擦应力,以及有无攻角下的流动结构。Rhee(2003)采用RANS方程结合k-ω湍流模型求解了SUBOFF主体带围壳和主体带尾翼的三向尾流场,并与试验结果进行了对比,发展了一种新的近壁面处理方法。Alin(2003)采用大涡模拟方法研究了SUBOFF全附体模型的绕流,计算给出了流线和壁面压力分布。Bensow(2004)采用大涡模拟的方法数值计算了SUBOFF的粘性流场,并与试验结果进行了对比,又利用鼓动盘模型模拟了潜艇模型的自航运动。在国内,中国船舶科学研究中心在粘流数值模拟潜艇流场与水动力方面做了大量工作,从回转体到翼板结合体、从裸艇体到全附体、从单独潜艇到艇桨耦合、从宏观流场到局部精细流场,进行了全面系统的计算与分析,这些工作为今后的研究提供了大量计算实例和技术支撑,积累的有益经验可资今后研究借鉴[8,9,10,11,12,13,14]。 对船体周围流动的试验研究一般是在拖曳水池或风洞中采用缩比模型进行试验测量分析,而模型尺度绕流与实船尺度绕流存在显著的差异,这样也就带来了从模型尺度向实船尺度外推的“尺度效应”问题,这是长期困扰流体力学界的一大难题。如果利用试验方式来研究尺度效应,势必造成巨大的人力和物力投入,成本很高,技术难度也大。目前由于CFD技术的工程实用性不断得到改善,对复杂流场的计算精度不断得到提高,所以人们已逐渐采用数值模拟的方法来研究船舶领域的尺度效应问题,主要集中于雷诺数对于船舶流场影响的研究。Bull与Watson(1998)利用SUBOFF潜艇模型,采用RANS方法,结合RNG k-ε、k-ω、SST k-ω湍流模型,研究了潜艇绕流随雷诺数的变化问题。网格数为50~100万。计算了带4只尾翼的SUBOFF潜艇(AFF4)在Re=1.2×107、1.2×108、1.2×109时的绕流,比较了不同雷诺数下桨盘面内速度波动图、速度等值线云图和湍流动能。Watson与Bull(1998)利用具有U型艉的HSVA油轮,采用RANS方法,结合RNG k-ε、k-ω、SST k-ω湍流模型,探讨了船舶绕流随雷诺数的变化问题。计算了Re=5×106、5×107、5×108、1×109时的绕流。网格数为80万。计算给出了艉部油流示踪显示和桨盘面内速度等值线云图,定性地比较了不同雷诺数下这些图形的变化。 在上述研究的基础上,近年来,欧盟在研究尺度效应方面做了巨大的投入,进行了更为深入、更为广泛的研究。欧盟建立了高雷诺数不可压缩流动研究计划(High Reynolds Number Impressible FlowEuclid 10.12)。这个计划的主要目的就是研发模拟高雷诺数船体绕流的计算平台,并且建立试验数据库以验证这个计算平台的可行性。这一计划选定了两艘军舰和一艘潜艇开展相应研究。采用激光多普勒测速仪测量实船螺旋桨上游的平均轴向和切向速度(Re=1.66~5.3×108),采用皮托管测量实船边界层。在拖曳水池(Re=5~8.35×106)和大型循环水槽(Re=0.674~5.67×107)中测量沿船体五个剖面内以及螺旋桨区域内的速度场,形成了丰富的数据库。在此基础上,评价了四个公司提供的七种不同的计算软件,比较了这些软件的能力、精度和可信度。上述计划说明,伴随着大规模并行计算技术和粘流数值模拟软件的快速发展,国际上已经开始专注于实尺度模型绕流的计算问题,以期为尺度效应问题的研究做些贡献。 2 计算值的方法 2.1 控制方程式 不可压缩流体连续性方程与RANS方程。 2.2 压力面及内载荷 L为潜艇总长。 速度入口:潜艇艏部向前1L,设定来流速度的大小与方向。 压力出口:潜艇艉部向后2L,设定相对于参考压力点的流体静压值。 壁面:潜艇外表面设定无滑移条件。 对称面:垂直于对称面的速度分量为零;平行于对称面的速度分量的法向导数为零。 外场:距潜艇表面约1L,速度为未受扰动的主流区速度。 2.3 湍流模型类型 本文共选用五种湍流模型。其中两方程湍流模型四种:Standard k-ε湍