FPSO耐波性运动和液舱晃荡全耦合数值模拟.PDFVIP

第六屆海峽兩岸水動力學研討會 FPSO耐波性运动与液舱晃荡全耦合数值模拟 * 庄园 万德成 (上海交通大学 船舶海洋与建筑工程学院 海洋工程国家重点实验室 高新船舶与深海开发装备协同创新中心, 上海 200240) *通信作者：dcwan@sjtu.edu.cn 摘要 带有液舱的船舶运动时不同自由度对液舱晃荡的影响不同，而被影响的液舱晃荡会 反过来改变船舶运动性能。本文采用 CFD 求解器 naoe-FOAM-SJTU 对带有液舱的 FPSO 进行全 耦合数值模拟，对放开横荡自由度和不放开横荡自由度下的船舶运动和液舱晃荡进行了模拟 研究；并对这两种状态下的船舶运动曲线和液舱抨击压力进行了比较。 關鍵字 液舱晃荡；全耦合；naoe-FOAM-SJTU 求解器；横荡运动 1 引言 对于带有自由液面的 LNG 液舱的船舶，在运动时会给予 LNG 液舱外部激励，LNG 液舱 受到激励产生晃荡，晃荡的流体将会同样影响船舶的运动。因此，研究船舶内外流场耦合现 象既可以有效预报船舶的水动力性能，又可以防止过于剧烈的舱内晃荡造成舱壁结构破坏。 对于船舶在波浪中运动的内外流场耦合情况，国内外许多研究者采用不同的方法进行了研 究。Nam 等人[1]采用实验方法对简化的 FPSO 模型进行了不同冲水率下的耦合效应研究，并应 用数值模拟的方法进行对比验证。随着计算机的发展，数值方法逐渐成为应用研究的热点。 而 CFD 方法因其可以模拟物体大幅度运动、流体的翻卷破碎等现象，成为数值模拟中不可或 缺的方法。近年来，模拟船舶内外流场耦合现象的数值方法，从应用势流理论转变为粘流方 法与势流方法的结合，即内流场采用 CFD 求解器进行求解，再将获得的力和力矩传递给船体， 船体在波浪中的运动则采用势流方法求解。Lee 和 Kim[2]，Li 和 Zhu[3]，Jiang 等人[4]等均采用 内流场（液舱）应用粘流方法求解，外流场（船舶运动）采用势流理论预报。Shen 和 Wan[5] 以及 Zhuang 和 Wan[6]应用基于 OpenFOAM 的自主开发的 CFD 求解器 naoe-FOAM-SJTU，进行 了船舶内外流场耦合数值模拟，实现了内外流场同时求解的全耦合计算。然而以上研究对于 船舶运动均只放开了三个自由度，即垂荡、横摇和纵摇。Eswaran 等人[7]仅对液舱进行了平面 运动激励以观察液舱的液面情况；Seo 等人[8]应用势流理论对带有液舱的船舶迎浪下放开纵荡 自由度进行了模拟研究；Mitra 等人[9]将船舶放开六个自由度，并考虑了风和流的影响，着重 考虑了横摇运动响应。Rognebakke 和 Faltinsen[10]进行了实验和数值模拟，对只放开横荡 自由 度的船舶与液舱耦合晃荡进行了研究。Zhao 等人[11]应用势流理论和实验模拟了二维 FLNG 和 - 1 - 第六屆海峽兩岸水動力學研討會 液舱晃荡在放开横荡自由度下的液舱的非线性自由面，Nasar 等人[11]采用实验方法将液舱固 定在驳船中，放开横摇垂荡和横荡自由度，观察其耦合现象。这些研究均表明，水平面运动 将会对液舱产生激励，且会使液舱产生剧烈晃荡，并改变船舶运动性能。 本文基于 Zhuang 和 Wan[6]的研究基础上，放开横荡自由度，并与之前所做三自由度运动 情况进行对比，观察横荡与横摇、垂荡和纵摇耦合情况下船舶内外流场的耦合情况。同时检 测了舱内的抨击压力，并与三自由度下的液舱抨击压力进行对比，分析舱内流体晃荡情况。 最后通过比较船舶运动与抨击压力的相位差，分析液舱晃荡对于船舶的减摇作用。 2 数值方法 为了充分考虑流体的粘性效应，本文采用求解不可压 Navier-Stocks 方程进行数值模拟，控 制方程如下： ∇U = 0