规则波中船舶迎浪航行波浪增阻数值计算-欢迎光临万德成教授研究团队.pdfVIP

第十三届全国水动力学学术会议暨第二十六届全国水动力学研讨会文集 规则波中船舶迎浪航行波浪增阻数值计算 * 查若思，沈志荣，万德成 摘要：本文采用基于 RANS 方程的黏性数值模拟方法，对某 Wigley 改型船型在规则波 中迎浪航行的阻力性能进行了研究与分析。使用了基于开源 CFD 软件 OpenFOAM 开发的 多功能求解器 naoe-FOAM-SJTU ，对航速为Fn=0.2 工况下船体的运动响应和波浪增阻进行 了计算，并与试验数据进行了对比验证。考虑了不同波长和不同波高的规则波对船体阻力 增值的影响，并通过运动幅值响应算子、自由面波形等对波浪增阻的影响因素进行了分析。 结果表明航速一定时，在 λ/L=1.1 附近阻力和运动响应达到峰值，另外短波范围内船体运 动辐射增阻的影响不大。 关键词：波浪增阻；运动响应；规则波；naoe-FOAM-SJTU 求解器 1 引言 近年来，IMO 提出的船舶能效指数 EEDI 已经开始正式施行。出于提高经济效益和加 强环境保护的双重要求，针对船舶能效的研究已成为备受瞩目的热门课题。船体在波浪中 的阻力增加会导致额外的燃料消耗和失速现象，影响其快速性和耐波性，因此精确预报波 浪增阻具有十分重要的意义。 国内外已有很多学者利用不同方法研究了波浪增阻，包括船模试验（EFD ）、基于势流 理论的理论计算和基于全黏性流场计算的计算流体力学（CFD ）数值模拟等方法。Maruo[1] 基于线性势流理论，将流场总的速度势分为入射势、辐射势和绕射势三部分，通过沿船体 中心线分布奇点来表达船型。Faltinsen[2]提出了短波渐近理论，研究了波浪作用于船体并产 生反射（绕射效应）作用于船体上的力。近年来日本专家提出了一种考虑短波影响的计算 增阻的方法[3]，并在 MEPC 第 59 次会议上提出用该方法计算失速系数。该方法首先利用 Maruo 理论计算船舶由于运动导致的增阻，然后再加上短波时船舶绕射修正项，最后得到 船舶的波浪增阻。 随着计算机处理能力的提高，三维时域分析方法得到了广泛应用。Dawson[4]提出了基 - 1361 - 第十三届全国水动力学学术会议暨第二十六届全国水动力学研讨会文集 于 Rankine 源的线性兴波数值方法，利用分布在船体表面和自由面上的 Rankine 源来满足 线性自由面条件和物面条件。Kim 等[5]分别使用了 WISH 和 AEGIR 程序计算了 Wigley 船 型、系列 60 和 S175 集装箱船在波浪中阻力增值，计算结果表明波长在接近船长附近辐射 增阻成分要远大于绕射增阻成分。但在短波和长波范围内，绕射增阻成分占较大比例。 近年来，CFD 方法得到了很大的发展，能很好地模拟大幅运动响应以及波浪破碎等强 非线性现象。Orihara 和 Miyata[6]使用 WISDAM-X 程序，对 SR108 集装箱船在波浪中航行 的运动响应和阻力增值进行了 CFD 模拟计算，结果显示该船型突出的船艏能够有效减小波 浪中的阻力增值。Simonsen 等[7]对 KCS 集装箱船在不同航速下的波浪增阻和船体运动响应 进行了试验研究和 CFD 计算预报。沈志荣和万德成等[8]利用基于开源代码 OpenFOAM 工 具箱和数据结构开发而成的 naoe-FOAM-SJTU 求解器计算分析了 Wigley III 、DTMB5415、 S175 等船型在迎浪中的运动响应和波浪增阻。 2 数值方法 本文利用 naoe-FOAM-SJTU 求解器[9]进行数值模拟求解。其控制方程为两相不可压 RANS 方程，包括连续性方程和动量方程： U 0 （1） U  p g x    f f （2） ( (U U )U)     