Processing math: 100%

结合Savitsky方法和重叠网格技术的滑行艇阻力数值计算与分析

郭军, 扈喆, 朱子文, 陈作钢, 崔连正, 李贵斌

郭军, 扈喆, 朱子文, 等. 结合Savitsky方法和重叠网格技术的滑行艇阻力数值计算与分析[J]. 中国舰船研究, 2022, 17(3): 126–134. DOI: 10.19693/j.issn.1673-3185.02417
引用本文: 郭军, 扈喆, 朱子文, 等. 结合Savitsky方法和重叠网格技术的滑行艇阻力数值计算与分析[J]. 中国舰船研究, 2022, 17(3): 126–134. DOI: 10.19693/j.issn.1673-3185.02417
GUO J, HU Z, ZHU Z W, et al. Numerical calculation and analysis of resistance performance of planing craft combining Savitsky method and overset grid technology[J]. Chinese Journal of Ship Research, 2022, 17(3): 126–134. DOI: 10.19693/j.issn.1673-3185.02417
Citation: GUO J, HU Z, ZHU Z W, et al. Numerical calculation and analysis of resistance performance of planing craft combining Savitsky method and overset grid technology[J]. Chinese Journal of Ship Research, 2022, 17(3): 126–134. DOI: 10.19693/j.issn.1673-3185.02417
郭军, 扈喆, 朱子文, 等. 结合Savitsky方法和重叠网格技术的滑行艇阻力数值计算与分析[J]. 中国舰船研究, 2022, 17(3): 126–134. CSTR: 32390.14.j.issn.1673-3185.02417
引用本文: 郭军, 扈喆, 朱子文, 等. 结合Savitsky方法和重叠网格技术的滑行艇阻力数值计算与分析[J]. 中国舰船研究, 2022, 17(3): 126–134. CSTR: 32390.14.j.issn.1673-3185.02417
GUO J, HU Z, ZHU Z W, et al. Numerical calculation and analysis of resistance performance of planing craft combining Savitsky method and overset grid technology[J]. Chinese Journal of Ship Research, 2022, 17(3): 126–134. CSTR: 32390.14.j.issn.1673-3185.02417
Citation: GUO J, HU Z, ZHU Z W, et al. Numerical calculation and analysis of resistance performance of planing craft combining Savitsky method and overset grid technology[J]. Chinese Journal of Ship Research, 2022, 17(3): 126–134. CSTR: 32390.14.j.issn.1673-3185.02417

结合Savitsky方法和重叠网格技术的滑行艇阻力数值计算与分析

基金项目: 福建省中青年教师教育科研项目(JAT200243,JAT190333);福建省自然科学基金资助项目(2021J05163)
详细信息
    作者简介:

    郭军,男,1992年生,博士,副教授。研究方向:船舶计算流体力学。E-mail:guojun6049@163.com

    扈喆,男,1988年生,博士,教授。研究方向:船舶与海洋结构物研发与应用。E-mail:zhehu@jmu.edu.cn

    朱子文,男,1991年生,博士,讲师。研究方向:计算流体力学。E-mail:zwzhu@jmu.edu.cn

    陈作钢,男,1967年生,博士,研究员。研究方向:船舶计算流体力学,风洞循环水槽研发与应用。E-mail:zgchen@sjtu.edu.cn

    通讯作者:

    陈作钢

  • 中图分类号: U661.31+1

Numerical calculation and analysis of resistance performance of planing craft combining Savitsky method and overset grid technology

知识共享许可协议
结合Savitsky方法和重叠网格技术的滑行艇阻力数值计算与分析郭军,采用知识共享署名4.0国际许可协议进行许可。
  • 摘要:
      目的  为提高数值预报精度,对滑行艇的静水阻力高精度数值模拟方法进行研究。
      方法  应用计算流体动力学(CFD)方法,结合Savitsky方法和重叠网格技术,对滑行艇在静水中的三维黏性流场进行数值模拟,并对不同载荷系数和航速下滑行艇的流场特性进行分析。
      结果  结果显示,滑行艇的阻力、升沉及纵倾角等计算结果与试验结果吻合良好,艇底的喷溅现象及水气分布模拟正常,表明采用所提方法可以准确、有效地预报滑行艇的阻力性能;随着载荷系数的增加,龙骨线压力系数的峰值增加,压力中心位置逐渐前移;随着航速的增加,龙骨线压力系数的峰值减小,压力中心位置逐渐后移,驻点线与中纵剖面的夹角减小,艇后“空穴”的深度减小、长度增大。
      结论  所做研究可为滑行艇阻力预报提供一种准确、有效的数值计算方法,能为滑行艇水动力性能数值研究提供技术支撑。
    Abstract:
      Objective  In this paper, the hydrostatic resistance of a planing craft is studied using the high-precision numerical simulation method to improve the numerical prediction accuracy.
      Methods  The three-dimensional viscous flow field of a planing craft in calm water is numerically simulated using the computational fluid dynamics (CFD) method combined with the Savitsky method and overset grid technique, and the flow field characteristics of the craft under different load coefficients and speeds are analyzed.
      Results  The calculated results of the resistance, sinkage and trim angle of the planing craft are in good agreement with the experimental results, and the spray phenomenon and distribution of water and air on the bottom of the craft are simulated normally, which shows that this method can accurately and effectively predict the resistance performance of planing craft. With the increase in the load coefficient, the peak value of the pressure coefficient on the keel increases, and the position of the pressure center moves forward. With the increase in speed, the peak value of the pressure coefficient on the keel decreases, the position of the pressure center gradually moves towards the stern, the angle between the stagnation line and the longitudinal section in the center plane decreases, the depth of the cavity behind the transom decreases, and the length of the cavity increases.
      Conclusions  This study provides an accurate and effective numerical calculation method for the resistance prediction of planing craft, and can provide technical support for the numerical study of the hydrodynamic performance of such craft.
  • 滑行艇是依靠流体的动力升力将艇体托出水面而使其处于滑行状态的一种高速艇,具有航速高、机动性好等特点。凭借这些优越的性能,滑行艇已被应用于导弹艇、鱼雷艇、侦察艇、交通艇和游艇等船型,在军用和民用领域具有广阔的应用前景。滑行艇在整个航行过程中航态变化大、阻力变化剧烈,如何快速、准确地预报滑行艇的阻力性能,是研究人员普遍关心的问题。

    早期的滑行艇阻力性能预报方法主要有半理论半经验公式方法(如Savitsky法[1]等)、根据系列滑行艇试验资料整理得出的图谱方法,以及滑行艇阻力模型试验。其中,半经验半理论方法和图谱方法的预报精度受限于计算船型与母型船的相似程度,而滑行艇的模型试验虽然能够直观地反映出艇体的受力情况以及流场特性,但成本高、耗时长。随着计算流体动力学(CFD)计算精度的不断提高,采用CFD技术不仅可以预报滑行艇的受力及运动情况,还可以捕捉艇体周围的流场细节,是研究滑行艇水动力问题的一种经济、有效的方法。近年来,越来越多的研究人员采用CFD技术对滑行艇进行数值模拟,并取得了一定的研究成果。Brizzolara和Serra[2]对棱柱型滑行艇流场的CFD数值模拟精度进行了研究,结果表明CFD技术可以用于滑行艇水动力问题的分析与计算。曹洪建[3]基于FLUENT软件预报了滑行艇在固定姿态下的实船阻力,结果显示计算误差约为10%,认为采用FLUENT软件对滑行艇进行阻力预报可行,但对升力的预报则不理想。Ghadimi等[4]采用RANS方法对滑行艇的姿态进行了数值模拟,结果显示RANS方法可以用于滑行艇的初步设计。马伟佳等[5]采用混合网格对滑行艇的阻力进行了数值模拟,结果显示数值计算结果与试验结果有一定的误差。蒋一[6]对某超高速三体滑行艇的快速性进行了数值计算,结果显示在速度较低时阻力的计算值与试验值吻合较好,但随着航速的增加,误差会逐渐增大,在最大航速下误差约为30%。Lotfi等[7]基于 CFX软件对某断级滑行艇的水动力性能进行了研究,结果显示其阻力预报精度约为5%。Frisk等[8]采用FLUENT软件和STAR-CCM+软件对滑行艇的阻力与姿态进行了数值计算,结果显示滑行状态下数值预报的阻力偏低,在半滑行状态下倾角的误差达32%。De Marco等[9]采用重叠网格技术和变形网格技术对滑行艇的水动力进行了分析,发现重叠网格技术表现较好。孙华伟等[10]采用STAR-CCM+软件探讨了船体表面第1层网格节点高度、船体表面网格尺度以及网格节点分布系数这3个因素对棱柱型滑行艇阻力计算精度、收敛速度及计算稳定性的影响。邵文勃等[11] 采用STAR-CCM+软件探究了时间步长分析与近壁面网格划分对滑行艇阻力计算精度的影响,结果显示时间步长和近壁面网格划分均会对滑行艇的摩擦阻力与艇底的水气分布造成影响。魏子凡等[12]采用NUMECA软件对翼滑艇进行了数值模拟,发现加水翼和防飞溅条能改善滑行艇的阻力性能。丁江明等[13]研究了在RANS计算中网格类型、网格尺寸等因素对滑行艇艇底水气分布、兴波形状、艇体姿态以及阻力的影响,并与船模阻力试验结果进行了比较。易文彬等[14]采用FLUENT,CFX 以及 STAR-CCM+软件探究了时间步长及网格数对预报结果的影响,发现采用CFX和STAR-CCM+软件模拟的滑行艇艇底表面有部分水体积分数约为0.5~0.8,艇底的水体积分数预报偏低,艇底水气分布异常,并认为较小的时间步长以及较细致的网格分布能够改善艇底水气分布的效果。李屺楠等[15]采用重新建模法对滑行艇绕流场进行了数值模拟,结果显示重新建模方法相比重叠网格方法其计算精度和计算效率更高。王慧等[16]采用重叠网格技术对滑行艇的兴波进行了模拟,结果显示采用CFD方法能够准确模拟滑行艇在静水中的航行姿态。

    综上所述,虽然现在采用CFD技术可以预报滑行艇的阻力性能,但在半滑行状态和滑行状态下,滑行艇艇体与流体之间存在剧烈的相互作用,在数值模拟中常存在艇底水气分布异常以及喷溅区域较难真实模拟的现象,且当滑行艇处于高速航行状态时其阻力预报与模型试验之间仍存在较大的误差,因此需进一步研究半滑行状态和滑行状态下的滑行艇阻力高精度数值预报方法。为此,本文拟基于CFD方法,首先采用Savitsky方法对滑行艇的航行姿态进行预估,再结合重叠网格技术对滑行艇的阻力性能进行预报,以提供一种准确、有效的滑行艇阻力数值预报方法,并分析滑行艇航态预估对阻力预报的影响;然后,在3种载荷工况及6个航速下对滑行艇的阻力性能进行数值预报,并与试验结果进行对比;最后,分析滑行艇的流场特性。

    本文将以Fridsma[17]于1969年进行的滑行艇系列模型试验中的一种模型为研究对象。该艇长1.143 m,宽0.2286 m,型深0.143 m,排水量7.26 kg,底部斜升角为10°。滑行艇的艇体型线和几何模型如图1所示。沿艇长方向,将滑行艇等分为20站,从第2站到20站之间的型线相同,第0站到第2站之间的型线如图1(a)所示,图中数值表示型线图的站位。

    图  1  滑行艇模型
    Figure  1.  Geometry model of planing craft

    滑行艇的无量纲速度为体积弗劳德数Fr

    Fr=Ug1/3 (1)

    式中:U为滑行艇航速;g为重力加速度;为排水体积。

    按照Fr的不同,滑行艇的运动可划分3种典型的航态:排水型状态(Fr<1.0)、半滑行状态(1.0<Fr<3.0)和滑行状态(Fr>3.0)。

    滑行艇的载荷系数C定义为:

    C=B3 (2)

    式中,B为滑行艇艇宽。

    滑行艇的不同载荷工况如表1所示。

    表  1  滑行艇的载荷工况
    Table  1.  Load conditions of planing craft
    工况载荷系数吃水/m重心纵向位置/m重心垂向位置/m
    10.3040.0260.4000.067
    20.6080.0420.4570.067
    30.9120.0570.4570.067
    下载: 导出CSV 
    | 显示表格

    对滑行艇水动力性能的数值模拟采用商业软件STAR-CCM+,不可压缩流体的湍流基本方程由连续性方程和动量方程组成[18]

    连续性方程:

    uixi=0 (3)

    动量方程:

    uit+xj(uiuj)=1ρpxi+υxj(uixj+ujxi)+xj(¯uiuj) (4)

    式中:uiuj为速度时均量;xixj为坐标;t为时间;ρ为流体密度; p为压力时均量;υ为流体运动黏性系数;¯uiuj为雷诺应力。

    本文的湍流模型采用SST k-ω湍流模型,该模型常被用于船舶水动力性能的数值计算,其表达形式为[18]

    ρkt+ρkuixi=xj[(μ+μtσk)kxj]+˜GkYk (5)
    ρωt+ρωuixi=xj[(μ+μtσω)ωxj]+GωYω+Dω (6)

    式中:μμt分别为流体的动力黏性系数和湍流黏性系数;σkσω分别为湍动能k和耗散率ω的湍流普朗特数;˜GkGω分别为kω的衍生项;YkYω分别为kω的耗散项;Dωω的交叉扩散项。

    在滑行艇阻力数值计算中,采用流体体积(VOF)[19]方法和高分辨率界面捕捉格式(high resolution interface capturing,HRIC)对自由液面进行捕捉,求解器采用基于分离流的黏性求解器,压力与速度的耦合运用SIMPLE法,对流项离散格式为二阶迎风格式。采用动态流体−固体相互作用模型(dynamic fluid-body interaction,DFBI)模拟滑行艇的六自由度运动,在计算中具有升沉和纵倾的自由度,滑行艇的释放时间为0.5 s,缓冲时间为2.5 s。采用隐式不定常格式,国际拖曳水池ITTC推荐的船舶阻力数值计算时间步长Δt与垂线间长LPP及航速U的关系为Δt =(0.005~0.01)LPP/U[20],本文针对滑行艇阻力计算采用的时间步长为0.005LPP/U

    由于模型几何及其绕流场关于中纵剖面基本对称,故数值模拟对象为滑行艇的半个绕流场。数值计算的计算域和边界条件如图2所示。计算域由静止的背景区域和随艇体运动的重叠区域组成,采用重叠网格法处理滑行艇高速航行时艇体姿态变化较大的问题。计算域的长度、宽度和高度分别为13Lpp,5Lpp和5Lpp,速度入口位于艇前2Lpp处,计算域顶部距离滑行艇1.0Lpp。将计算域入口、侧面、顶部和底部设为速度入口边界,计算域的出口设为压力出口边界,艇体表面设为无滑移壁面边界。在计算域的出口前设置消波阻尼以消除传至出口处的船行波,为避免波的反射,阻尼长度为3Lpp[21, 22]

    图  2  计算域和边界条件
    Figure  2.  Computational domain and boundary conditions

    网格类型为切割体非结构化网格,在滑行艇周围和自由面等区域进行了网格细化。在滑行艇周围设置3个开尔文波系形状的网格加密区以进行局部加密,在设计水线处设置3个长方体的网格加密区以对自由面网格进行局部加密,在背景区域与重叠区域之间设置1个网格加密区以进行网格过渡,以使背景区域和重叠区域的网格尺寸处于相同量级。滑行艇的表面网格尺寸为1.75%LPP ,滑行艇边界层网格使用棱柱形网格进行划分,共15层,增长率为1.2,边界层的厚度δ为:

    δ=0.382LWRe1/5LW (7)

    式中: LW为滑行艇的平均浸湿长度;Re为雷诺数。

    计算域的自由面网格和滑行艇表面网格的分布如图3所示。背景区域和重叠区域的网格量分别为90万和141万,最终的计算域网格总数为231万。

    图  3  计算域网格
    Figure  3.  Grid of computational domain

    滑行艇在排水型状态下的艇体姿态与静浮时差别不大,而半滑行状态和滑行状态下的艇体姿态则与静浮时差别较大。由于重叠区域内的网格是随滑行艇姿态的变化而做整体运动,此时艇体周围及自由面等区域需要更大的网格加密区,但这样会极大地增加网格数量,降低计算效率。因此在开始计算时,调整滑行艇的航态以使最终计算的航态接近于初始航态,这样不仅能降低网格数量,提高计算效率,还能改善计算的收敛性。此外,由于网格的方向与来流方向较为接近,故还能改善艇底的水气分布异常现象[23],从而提高计算精度。本文将采用半理论半经验公式Savitsky方法[24-25]对滑行艇的航态进行估算。

    CLβ=FLβ0.5ρU2B2=CL00.0065βC0.6L0 (8)

    式中:CLβ为升力系数;FLβ为升力(滑行状态时约等于重力);β为斜升角;CL0为零斜升角时的升力系数,

    CL0=θ1.1(0.012λ0.5W+0.0055λ2.5W/Fr2B) (9)

    式中:θ为滑行艇纵倾角;λW为平均浸湿长宽比;FrB为艇宽弗劳德数,FrB = U/(gB)0.5

    xgλWB=0.7515.21Fr2B/λ2W+2.39 (10)

    式中,xg为重心距方尾的纵向距离。

    开始计算时,在不同载荷系数下设置的滑行艇初始航态如图4所示。

    图  4  滑行艇的航态估算
    Figure  4.  Attitude estimation of planing craft

    滑行艇的总阻力系数Ct定义为:

    Ct=Rt0.5ρAwetU2 (11)

    式中:Rt为滑行艇总阻力;Awet为滑行艇的动态湿表面积。

    在载荷系数C= 0.608及航速U = 4 m/s(Fr= 2.904)的工况下,对滑行艇进行数值计算,分析航态预估对滑行艇阻力预报的影响,计算结果与试验结果的比较如表2(表中,σ为滑行艇的深沉,Lwet为浸湿长度)所示。由表可见,与单纯采用重叠网格相比,结合Savitsky经验公式进行航态预估后,滑行艇航态、阻力和浸湿长度的预报精度均有所提高。

    表  2  航态预估对阻力预报的影响
    Table  2.  Influence of trim estimation on resistance prediction
    方法σ/Bθ/(°)CtLwet /B
    试验结果 (EFD)0.0644.150.006823.41
    重叠网格0.0583.940.007523.54
    Savitsky+重叠网格0.0604.070.007253.47
    下载: 导出CSV 
    | 显示表格

    在3种载荷工况及6个航速(U = 1,2,3,4,5,6 m/s)下,对滑行艇的阻力性能进行数值计算,并与阻力试验结果进行比较,以验证数值方法的准确性及有效性。

    滑行艇升沉计算结果与试验结果的对比如图5所示。由图5可见,升沉的计算结果稍微偏小,3种载荷工况下升沉的平均误差分别为0.12%,0.31%和−5.06%。滑行艇纵倾角计算结果与试验结果的对比如图6所示。由图6可见,纵倾角的计算结果稍微偏小,3种载荷工况下纵倾角的平均误差分别为−0.62%,−2.11%和−5.86%。滑行艇总阻力计算结果与试验结果的对比如图7所示。由图7可见,阻力的数值计算结果偏大,3种载荷工况下阻力的平均误差分别为9.03%,4.09%和6.93%。滑行艇平均浸湿长度计算结果与试验结果的对比如图8所示。由图8可见,3种载荷工况下浸湿长度的平均误差分别为−0.63%,1.23%和1.61%。整体上,滑行艇的升沉、纵倾角、总阻力和平均浸湿长度的计算结果均与试验值吻合良好,说明本文所提的数值计算方法有效,数值计算结果准确可靠,可用于半滑行状态和滑行状态下滑行艇的阻力预报。

    图  5  升沉计算值与试验值对比
    Figure  5.  Comparison of heaving between CFD and EFD results
    图  6  纵倾计算值与试验值对比
    Figure  6.  Comparison of trim between CFD and EFD results
    图  7  阻力计算值与试验值对比
    Figure  7.  Comparison of resistance between CFD and EFD results
    图  8  平均浸湿长度计算值与试验值对比
    Figure  8.  Comparison of average wet length between CFD and EFD results

    滑行艇压力系数CP的定义为:

    CP=PP00.5ρU2 (12)

    式中:P为滑行艇表面的绝对压力;P0为参考点压力(1个大气压)。

    当航速U = 4 m/s时,在不同的载荷工况下,滑行艇艇体表面的压力分布规律如图9所示,其中图9(a)为滑行艇龙骨线的压力分布图,图9(b)为滑行艇艇体表面压力分布云图。图中,x为沿船长方向的位置。从中可以看出,龙骨线的压力系数在与水接触的位置显著增加,在驻点处时达到峰值,然后减小并达到稳定,在尾部逐渐减小至负值,至艉封板时达到最小;当C = 0.304,0.608时,压力系数的峰值比较接近,导致滑行艇的纵倾角也比较接近;当C = 0.912时,压力系数的峰值增加,压力中心位置前移,导致滑行艇的纵倾角增加。

    图  9  不同载荷工况下的压力分布
    Figure  9.  Pressure distribution at different load conditions

    当载荷系数C = 0.912时,不同航速下滑行艇艇体表面的压力分布规律如图10所示。从中可以看出,当Fr=1.357时,滑行艇处于半滑行状态的初始阶段,在喷溅区域其压力系数是先快速增加至最大,随后迅速减小至常压,然后缓慢增加,在艉部又迅速减小至负值;当Fr = 2.714时,滑行艇处于半滑行状态的结束阶段,此时龙骨线压力分布的变化规律与滑行状态的变化规律基本相同;随着航速的增加,驻线压力系数的峰值减小,驻点线逐渐向艉部方向移动,驻点线与中纵剖面的夹角逐渐减小。

    图  10  不同航速下的压力分布
    Figure  10.  Pressure distribution at different speeds

    当载荷系数C = 0.912时,不同航速下滑行艇中纵剖面的自由面波高对比如图11(图中,z为自由面的波高位置)所示。滑行艇的自由面波形如图12所示,其中左图为整体波形,右图为喷溅区域波形。由图可见,随着航速的增加,船行波的角度逐渐减小,艏部兴波的高度逐渐降低,艇后鸡尾状水丘的高度逐渐增大,且逐渐向后移动,艇后“空穴”的深度也逐渐减小,而“空穴”的长度则逐渐增大;当Fr = 1.357时,在滑行艇艏部发生了明显的喷溅现象,水流冲击滑行艇艏部后向前翻卷,喷溅的浸湿面积较大;当Fr = 2.714时,在艇体底部发生了喷溅现象,喷溅区域的浸湿形状接近于三角形;当Fr = 4.071时,喷溅区域的浸湿形状已完全为三角形。

    图  11  不同航速下y = 0剖面的自由面波高对比
    Figure  11.  Wave pattern of y = 0 at different speeds
    图  12  不同航速下的自由波面分布
    Figure  12.  Free wave surface distribution at different speeds

    当载荷系数C = 0.912时,不同航速下滑行艇表面的水气分布如图13所示。由图可见,滑行艇艇底与水接触部分水相的体积分数接近于1.0,说明本文艇底水气分布预报正常,未出现艇底水气分布模拟异常的现象,这从另一个角度也说明了本文滑行艇阻力数值模拟的准确性。

      13  不同航速下的滑行艇表面的水气分布
      13.  Volume fraction of water of planing craft at different speeds

    本文在求解RANS方程的基础上,首先采用Savitsky方法对滑行艇的航行姿态进行了预估,然后结合重叠网格技术对滑行艇的阻力性能进行预报,并提供了一种滑行艇阻力高精度数值预报方法,接着,在此基础上对3种载荷工况下滑行艇的阻力性能进行数值预报,并与试验结果进行了对比,最后,分析了滑行艇的流场特性。通过研究,主要得出以下结论:

    1) 随着载荷系数的增加,龙骨线压力系数的峰值增加,压力中心位置前移,滑行艇的纵倾角增加。

    2) 随着航速的增加,龙骨线压力系数的峰值减小,驻点线逐渐后移,驻点线与中纵剖面的夹角减小;随着航速的增加,滑行艇船行波的角度逐渐减小,艏部兴波逐渐减弱,艉部鸡尾流的高度逐渐增大,艇后“空穴”的深度逐渐减小,但长度逐渐增大。

    3) 本文模拟的艇底水气分布正常,能够模拟出滑行艇的喷溅现象,且滑行艇阻力数值计算结果与试验结果吻合良好,这些均说明本文数值方法准确可靠,可为滑行艇水动力数值研究提供有力的技术支撑。

    本文得到了中国船舶及海洋工程设计研究院喷水推进技术重点实验室的支持,在此深表谢意!

  • 图  1   滑行艇模型

    Figure  1.   Geometry model of planing craft

    图  2   计算域和边界条件

    Figure  2.   Computational domain and boundary conditions

    图  3   计算域网格

    Figure  3.   Grid of computational domain

    图  4   滑行艇的航态估算

    Figure  4.   Attitude estimation of planing craft

    图  5   升沉计算值与试验值对比

    Figure  5.   Comparison of heaving between CFD and EFD results

    图  6   纵倾计算值与试验值对比

    Figure  6.   Comparison of trim between CFD and EFD results

    图  7   阻力计算值与试验值对比

    Figure  7.   Comparison of resistance between CFD and EFD results

    图  8   平均浸湿长度计算值与试验值对比

    Figure  8.   Comparison of average wet length between CFD and EFD results

    图  9   不同载荷工况下的压力分布

    Figure  9.   Pressure distribution at different load conditions

    图  10   不同航速下的压力分布

    Figure  10.   Pressure distribution at different speeds

    图  11   不同航速下y = 0剖面的自由面波高对比

    Figure  11.   Wave pattern of y = 0 at different speeds

    图  12   不同航速下的自由波面分布

    Figure  12.   Free wave surface distribution at different speeds

    13   不同航速下的滑行艇表面的水气分布

    13.   Volume fraction of water of planing craft at different speeds

    表  1   滑行艇的载荷工况

    Table  1   Load conditions of planing craft

    工况载荷系数吃水/m重心纵向位置/m重心垂向位置/m
    10.3040.0260.4000.067
    20.6080.0420.4570.067
    30.9120.0570.4570.067
    下载: 导出CSV

    表  2   航态预估对阻力预报的影响

    Table  2   Influence of trim estimation on resistance prediction

    方法σ/Bθ/(°)CtLwet /B
    试验结果 (EFD)0.0644.150.006823.41
    重叠网格0.0583.940.007523.54
    Savitsky+重叠网格0.0604.070.007253.47
    下载: 导出CSV
  • [1]

    SAVITSKY D. Hydrodynamic design of planing hulls[J]. Marine Technology, 1964, 1(4): 71–95.

    [2]

    BRIZZOLARA S, SERRA F. Accuracy of CFD codes in the prediction of planing surfaces hydrodynamic characteristics[C]//2nd International Conference on Marine Research and Transportation, 2007.

    [3] 曹洪建. 基于FLUENT的滑行艇阻力计算研究[D]. 哈尔滨: 哈尔滨工程大学, 2008.

    CAO H J. The computation and research on resistance of planing craft on the software FLUENT[D]. Harbin: Harbin Engineering University, 2008 (in Chinese).

    [4]

    GHADIMI P, MIRHOSSEINI S H, DASHTIMANESH A, et al. RANS simulation of dynamic trim and sinkage of a planing hull[J]. Applied Mathematics and Physics, 2013, 1(1): 6–10.

    [5] 马伟佳, 庞永杰, 孙华伟, 等. 混合网格在滑行艇阻力数值模拟中的应用[J]. 船舶工程, 2013, 35(4): 8–10, 58.

    MA W J, PANG Y J, SONG H W, et al. Application of mixed grid in numerical simulation of planning-hull resistance[J]. Ship Engineering, 2013, 35(4): 8–10, 58 (in Chinese).

    [6] 蒋一. 基于CFD的超高速三体滑行艇快速性分析[D]. 哈尔滨: 哈尔滨工程大学, 2013.

    JIANG Y. CFD-based analysis on speed performance of high-speed trimaran planing boat[D]. Harbin: Harbin Engineering University, 2013 (in Chinese).

    [7]

    LOTFI P, ASHRAFIZAADEH M, ESFAHAN R K. Numerical investigation of a stepped planing hull in calm water[J]. Ocean Engineering, 2015, 94: 103–110. doi: 10.1016/j.oceaneng.2014.11.022

    [8]

    FRISK D, TEGEHALL L. Prediction of high-speed planing hull resistance and running attitude[D]. Gothenburg, Sweden: Chalmers University of Technology, 2015.

    [9]

    DE MARCO A, MANCINI S, MIRANDA S, et al. Experimental and numerical hydrodynamic analysis of a stepped planing hull[J]. Applied Ocean Research, 2017, 64: 135–154. doi: 10.1016/j.apor.2017.02.004

    [10] 孙华伟, 马伟佳, 朱江波. 影响滑行艇阻力数值计算的网格因素研究[J]. 中国造船, 2015, 56(2): 170–178.

    SUN H W, MA W J, ZHU J B. Research on grid factor in numerical calculation of planing craft resistance[J]. Shipbuilding of China, 2015, 56(2): 170–178 (in Chinese).

    [11] 邵文勃, 马山, 段文洋, 等. 基于CFD技术的滑行艇静水阻力计算[J]. 船舶工程, 2019, 41(9): 41–45,137.

    SHAO W B, MA S, DUAN W Y, et al. Calm-water resistance calculation of planing craft based on CFD method[J]. Ship Engineering, 2019, 41(9): 41–45,137 (in Chinese).

    [12] 魏子凡, 井升平, 杨松林. 新型高速艇的CFD模拟和对比分析[J]. 中国舰船研究, 2016, 11(4): 22–28.

    WEI Z F, JING S P, YANG S L. CFD simulation and comparison analysis of a new type high-speed boat[J]. Chinese Journal of Ship Research, 2016, 11(4): 22–28 (in Chinese).

    [13] 丁江明, 江佳炳, 秦江涛, 等. 高速滑行艇阻力性能RANS计算中网格影响因素[J]. 哈尔滨工程大学学报, 2019, 40(6): 1065–1071.

    DING J M, JIANG J B, QIN J T, et al. Influencing mesh factors in the calculation of the resistance performance of high-speed planing crafts through RANS[J]. Journal of Harbin Engineering University, 2019, 40(6): 1065–1071 (in Chinese).

    [14] 易文彬, 王永生, 彭云龙, 等. 滑行艇阻力数值预报若干影响因素研究[J]. 华中科技大学学报(自然科学版), 2017, 45(9): 120–126.

    YI W B, WANG Y S, PENG Y L, et al. Research on several influence factors in numerical prediction of planning craft resistance[J]. Journal of Huazhong University of Science and Technology (Natural Science Edition), 2017, 45(9): 120–126 (in Chinese).

    [15] 李屺楠, 秦江涛, 周利兰. 基于重新建模法的滑行艇阻力数值计算[J]. 船舶工程, 2020, 42(1): 42–46, 121.

    LI Q N, QIN J T, ZHOU L L. Numerical calculation of planing boat resistance based on remesh method[J]. Ship Engineering, 2020, 42(1): 42–46, 121 (in Chinese).

    [16] 王慧, 朱仁传, 杨云涛, 等. 基于CFD的滑行艇兴波与姿态模拟分析[J]. 中国造船, 2020, 61(3): 1–14.

    WANG H, ZHU R C, YANG Y T, et al. Simulation and analysis of wave-making and attitudes of planing hull by CFD[J]. Shipbuilding of China, 2020, 61(3): 1–14 (in Chinese).

    [17]

    FRIDSMA G. A systematic study of the rough-water performance of planing boats. irregular waves-part II[R]. New Jersey: Davidson Laboratory Stevens Institute of Technology, 1969.

    [18] 王福军. 计算流体动力学分析——CFD软件原理与应用[M]. 北京: 清华大学出版社, 2004: 7-12.

    WANG F J. Computational fluid dynamics analysis—principles and applications of CFD software[M]. Beijing: Tsinghua University Press, 2004: 7-12 (in Chinese).

    [19]

    HIRT C W, NICHOLS B D. Volume of fluid (VOF) method for the dynamics of free boundaries[J]. Journal of Computational Physics, 1981, 39(1): 201–225. doi: 10.1016/0021-9991(81)90145-5

    [20]

    ITTC. The specialist committee on CFD in marine hydrodynamics-recommended procedures and guidelines: practical guidelines for ship CFD applications[C]//Proceeding of the 26th ITTC. Rio de Janeiro, Brazil: ITTC, 2011.

    [21]

    GUO J, CHEN Z G, DAI Y X. Numerical study on self-propulsion of a waterjet propelled trimaran[J]. Ocean Engineering, 2020, 195: 106655. doi: 10.1016/j.oceaneng.2019.106655

    [22]

    GUO J, ZHANG Y, CHEN Z G, et al. CFD-based multi-objective optimization of a waterjet-propelled trimaran[J]. Ocean Engineering, 2020, 195: 106755. doi: 10.1016/j.oceaneng.2019.106755

    [23] 李昆鹏, 魏成柱, 梁晓锋. 多面体网格在滑行艇数值仿真计算中的应用[J]. 舰船科学技术, 2020, 42(2): 33–37.

    LI K P, WEI C Z, LIANG X F, et al. Application of polyhedral mesh in numerical simulations of planing hulls[J]. Ship Science and Technology, 2020, 42(2): 33–37 (in Chinese).

    [24] 孙源, 卢晓平, 李井煜, 等. 滑行艇阻力计算方法对比研究[J]. 中国舰船研究, 2019, 14(1): 27–32.

    SUN Y, LU X P, LI J Y, et al. Comparative study on simulation of resistance for planing craft[J]. Chinese Journal of Ship Research, 2019, 14(1): 27–32 (in Chinese).

    [25]

    FALTINSEN O M. Hydrodynamics of high-speed marine vehicles[M]. New York: Cambridge University Press, 2005: 516-530.

  • 期刊类型引用(1)

    1. 程宏佳,闫冬花,沈恒,徐云涛,周爱春,张莉. 基于Savitsky法的平底滑行艇阻力估算影响因素研究. 现代制造技术与装备. 2023(12): 23-25 . 百度学术

    其他类型引用(2)

  • 其他相关附件

图(14)  /  表(2)
计量
  • 文章访问数:  939
  • HTML全文浏览量:  246
  • PDF下载量:  98
  • 被引次数: 3
出版历程
  • 收稿日期:  2021-06-15
  • 修回日期:  2021-09-13
  • 网络出版日期:  2021-09-14
  • 刊出日期:  2022-06-29

目录

/

返回文章
返回