1、大型邮轮泳池形状对其晃荡性能影响研究大型邮轮泳池形状对其晃荡性能影响研究吴思莹，丁悦，黄毅铭(上海外高桥造船有限公司，上海200137)摘 要:针对大型邮轮泳池晃荡的问题，基于 CFD 软件 Fluent 对不同形状泳池进行晃荡分析，研究影响泳池晃荡的因素。为获得准确结果，对邮轮水动力及泳池晃荡研究分别基于势流理论与粘流理论。研究结果表明：通过结合势流粘流理论可以高效准确的模拟泳池晃荡；3 级海况下，对于矩形、圆形/椭圆以及棱形的泳池，横向宽度越大泳池晃荡引起的横向力越大，液面最大高度较高；纵横向过度越平缓引起横向力越小；相同海况下，棱形与椭圆形截面泳池晃荡幅度最小，可提高恶劣海况时泳池的舒适
2、性。关键词：邮轮；泳池；舒适性；CFD中图分类号：U675文献标识码：A文章编号：16727649(2023)12003104doi：10.3404/j.issn.16727619.2023.12.006Research on the influence of swimming pool shape of large cruise ship onits sloshing performanceWUSi-ying,DINGYue,HUANGYi-ming(ShanghaiWaigaoqiaoShipbuildingCo.,Ltd.,Shanghai200137,China)Abstract:Ai
3、mingatthesloshingproblemoflargecruiseswimmingpool,thesloshingofswimmingpoolswithdiffer-entshapesisanalyzedbasedonCFDsoftwareFluent,andthefactorsaffectingthesloshingofswimmingpoolarestudied.Inordertoobtainaccurateresults,potentialflowtheoryandviscousflowtheoryareusedforcruisehydrodynamicforcesandswim
4、mingpoolsloshingrespectively.Theresultsshowthatthesloshingofswimmingpoolcanbesimulatedefficientlyandaccuratelybycombiningthepotentialflowandviscousflowtheory.Forrectangular,circularandellipticalswimmingpools,thegreaterthetransversewidth,thegreaterthetransverseforcecausedbyswimmingpoolsloshing,andthe
5、higherthemax-imumheightofliquidlevel.Underthesameseaconditions,theswimmingpoolswithprismaticandellipticalsectionshavethesmallestsloshingrange,whichcanimprovethecomfortoftheswimmingpoolinbadseaconditions.Key words:cruise；swimmingpool；comfort；CFD0引言近年来，随着人民生活质量和消费水平不断提高，邮轮旅游已成为我国新型的消费产业，预计到2030 年，我国邮轮
6、出入境游客人数将达 1750 万人次1，成为邮轮产业最大的消费市场之一。同时，旅游方式也将从旅游观光为主向休闲娱乐的转变，因此邮轮旅游的舒适性与多样性变得尤为重要2。其中，邮轮泳池多布置在船首以及尾部，并附带多功能如烧烤台，海滩派对，日光浴场等，是游客休闲放松的主要区域。然而，泳池的自由液面会影响船舶稳性降低舒适性，同时泳池的晃荡也会对乘客造成不适与危险，所以一般船舶泳池布置在船舶中心线附近3，其尺寸也会受到限制。许多邮轮还规定在 3 级海况时禁止开放游泳池。为提高乘客舒适性，减少高海况对乘客旅游观光的不利影响，本文将对不同形状泳池进行晃荡分析，研究影响泳池晃荡的因素，提出舒适性更高的泳池形状
7、。对于晃荡问题，国内外专家学者开展了大量的分析研究。早期的液舱晃荡问题大多依赖于模型试验的方法。Nasar 等4对波浪中单液舱浮箱进行试验研究，设置不同波高、周期等参数和载液量，研究晃荡现象。卫志军等5开展了单个矩形液舱晃荡的大尺度模第 45卷第12期舰船科学技术Vol.45,No.122023年6月SHIPSCIENCEANDTECHNOLOGYJun.,2023收稿日期:20220711基金项目:工信部高技术船舶科研项目（KD1-223）作者简介:吴思莹(1991)，女，硕士，工程师，研究方向为船舶总体性能及设计。型试验，研究不同载液量下舱内液体晃荡冲击荷载的特性。随着理论方法与数值方法的
8、不断完善，基于数值计算的研究方法不断被采用。Jiang 等6针对某型双液舱船舶，采用 VOF 法(volumeoffluid)和 IRF 法(im-pulseresponsefunction)求解船舶与液舱晃荡间的耦合响应，研究了不同载液工况下船舶运动响应与液舱晃荡间的耦合效应。丁刚7采用非结构混合网格，基于Fluent 软件中的 VOF 模型、Segregated，以及压力速度祸合的 PISO 算法，模拟三维条件下具有自由液面的、不互溶的、油一水一气三相不可压缩非定常流动。徐博等8基于 CFD 软件 STAR-CCM+建立一种二维黏性流体数值模型，研究了矩形液舱在特定水平余弦激励下舱内液体的
9、晃动问题。Stephen 等9建立考虑液舱耦合效应的时域运动方程，求解了规则波作用下矩形浮体的液舱晃荡问题，研究不同载液量对浮体水动力响应的影响。DongmingLiu 等10基于 CFD 方法，采用大涡模拟（LES）研究了三维液舱晃荡液面变化情况，采用 VOF 追踪自由液面，计算结果与试验吻合度高。目前研究液舱晃荡问题多以模型试验结合数值仿真为主，经验证 CFD 方法能够与试验结果很好吻合，尤其对于非线性现象能够准确预测。基于势流理论的方法可在时域范围内求解液舱与浮体耦合运动，在线性条件下可以很好预测响应结果且效率较高。当浮体出现抨击：液面破碎等强非线性现象时势流理论将不再适用。因此，根据这
10、一结论，本文在求解大型邮轮水动力响应时采用势流理论，将得到的结果作为液舱的强迫运动，基于 CFD 方法模拟液舱晃荡，得到更多流动细节。在此基础上对不同形状泳池进行晃荡分析，研究影响泳池晃荡的因素，提出舒适性更高的泳池形状。1邮轮水动力分析 1.1 模型建立以某大型邮轮为例，其主尺度参数如图 1 所示。采用势流软件 Ansys-AQWA 对其水动力性能进行分析，由于船舶排水量与泳池重量相差甚大，故计算忽略泳池对邮轮运动响应的影响。参照目前研究方法，邮轮三维模型及网格图如图 1 所示，总网格数为7902，坐标系原点建立在船体中纵剖面与水线面的交线尾部。考虑到此类船舶纵向摇动远小于横向摇动，本文仅考
11、虑横向最危险工况，即重点研究横浪作用下船舶在固有周期时的横摇及其造成的液舱晃荡。根据经验公式，横摇惯性半径选取 0.36B，按规范增加人工阻尼。计算邮轮横浪作用时的横摇 RAO，为节约计算时间去除主甲板以上上层建筑。1.2 计算结果考虑横浪作用下，船舶横摇响应幅值算子如图 2所示。图 2 为邮轮横摇 RAO 曲线。可看出单位波幅下邮轮最大响应幅值可达 3.3，此时波浪激励力频率为0.5rad/s。考虑到大型邮轮尺寸较大，单位波幅作用时满足线性微幅波条件，当波高提升至 1.25m（对应3 级最高海况）时，邮轮最大响应幅值可达 2.3，对应频率仍为 0.5rad/s。因此，考虑最危险工况，即邮轮横
12、向受频率为 0.5rad/s 的规则波，则邮轮上任一点绕重心旋转角度为：=0.04cos(0.5t)rad/s。(1)将 式（1）作 为 泳 池 绕 重 心 的 强 迫 运 动 输 入CFD 软件即可得到 3 级海况下不同形状泳池内液面升高及压力变化情况。表 1 主尺度参数Tab.1Mainscaleparameters参数数值垂线间长/m304.00型宽/m40.00设计吃水/m8.00排水量/m3130000泳池尺寸量/m13101.5图1船体模型Fig.1Modelofshiphull图2横浪作用下邮轮横摇 RAOFig.2RollingRAOundertransversewaves32
13、舰船科学技术第45卷2泳池晃荡分析为深入研究泳池晃荡性能，基于 CFD 软件 Flu-ent 计算不同形状泳池 3 级海况下液面变化情况，为泳池设计及提高舒适性措施提供参考。CFD 具有丰富的物理模型，功能强大，精度比较高，用来捕捉液面变化规律具有较强的优势。2.1 基础理论2.1.1流体的运动控制方程根据连续介质假设，流体质量守恒方程（又称为连续性方程）可表达为：t+(vx)x+(vy)y+(vz)z=0，(2)tvxvyvzxyz式中：为流体的密度；为时间；，分别为流体速度矢量沿 轴，轴，轴的分量。动量守恒方程（Navier-Stokes 方程）可表达为：|(vx)t+(vx)=px+xx
14、x+yxy+zxz+fx，(vy)t+(vy)=py+xyx+yyy+zyz+fy，(vz)t+(vz)=pz+xzx+yzy+zzz+fz。(3)pfg式中：表示流体微团所受的压力；为 3 个坐标轴方向的外力分量，若只考虑重力，即仅有垂直方向的分量为，其他方向的分量均为 0。2.1.2强迫运动方法1）指定加速度。Fluent 中可以指定流体运动加速度，并将加速度以体积力的形式施加到计算域中的流体上。然而 Fluent 并没有提供变加速度的直接添加，若要计算变加速度情况，则需要手工分段计算。2）指定计算域速度。将加速度或位移数据转化为速度添加到计算域上。可以通过 DEFINE_ZONE_MO-
15、TION 宏或 PROFILE 文件的方式进行指定。2021 版 An-sys-Fluent 可通过函数直接定义域的运动此方式要比加速度方式更加灵活，本文采用此方法。2.2 数值计算2.2.1模型建立以矩形泳池为例，参数设置与实际情况一致，四周采用固壁边界条件，顶部设置大气压力出口，计算模型如图 3 所示。网格总数为 190 万（经无关性验证），自由液面处加密。启用 movingmesh，设置泳池绕重心的强迫运动，为保证计算精度采用 PISO 压力、流速耦合算法。保证面积不变，改变泳池形状采用相同设置，如图 3 所示。2.2.2计算结果k采用 VOF 模型，时间步长经验证可兼顾效率与准确性，R
16、NG湍流模型基于 PISO 算法模拟泳池 3 级海况下液面变化及泳池受力情况，计算结果如图 4 所示。图4矩形泳池计算结果Fig.4Calculationresultsofrectangularpool可知，由于重心距泳池垂向距离较小，泳池内液面升高不明显，泳池两侧受到的压力交替变化可比较明显地看出晃荡带来的影响，横向力变化频率与波浪激励力保持一致符合实际。监测点位于泳池壁面液面在 1.4801.522m 之间，游客受到影响较小。3不同形状泳池晃荡对比为进一步提高泳池舒适性，对目前市内泳池形状进行总结，通过上述方法研究不同形状泳池晃荡性能，通过对比分析最终确定有利于乘客舒适性的截面形状。不同形状泳池横向力曲线如图 5 所示。图3计算域网格模型Fig.3Computationaldomain第45卷吴思莹，等：大型邮轮泳池形状对其晃荡性能影响研究33图 6 为不同形状泳池晃荡结果。可以看出：1）受力曲线总体上变化规律一致，均在晃动开始68s 后达到最大值。稳定后矩形、棱形与椭圆受到横向力较小，圆形最大。2）液面变化曲线波动规律性较差，这主要是因为泳池整体晃荡较小任何一点小的扰动都会引起液