孤立波与固定海洋结构物相互作用数值模拟-上海交通大学.pdfVIP

孤立波与固定海洋结构物相互作用数值模拟1) 曹洪建，刘远传，万德成2) （上海交通大学 船舶海洋与建筑工程学院 海洋工程国家重点实验室，上海 200240) 摘要：基于开源程序库OpenFOAM 数值求解N-S 方程，采用VOF 方法捕捉自由面，开发了粘性数值水池。并采 用模拟推板运动的方式生成孤立波，与孤立波理波形比较验证了本文方法的有效性。基于粘性数值水池，实现了 孤立波与大尺寸直立圆柱结构的相互作用的数值模拟，重点研究了孤立波对直立圆柱的冲击现象和爬高等现象。 通过数值模拟得到了一系列不同波高的孤立波在圆柱前的最大爬高值，并计算分析了结构物受到的波浪载荷，为 海洋结构物的设计提供了很好的支持。 关键词：孤立波；圆柱；爬高；波浪载荷；OpenFOAM 引 言 近年来人类开采海洋油气资源的范围逐渐由近海扩大到深海，使得研究恶劣海况下海洋工程结构物的 安全性、稳定性变得越发重要。当海洋结构物遭遇深海恶劣海浪的冲击作用时，巨大的冲击力将严重影响 海洋平台的安全性，然而如何确定极限海况下的波浪对结构物的载荷是海洋结构物设计和应用的关键。目 前对于非线性很强的极限波浪与大型结构物的相互作用，采用传统的势流理论难以得到准确的结果。因此， 通常都会进行水池模型实验获取需要的数据。随着CFD 技术的不断发展，应用数值模拟方法构建粘性数值 波浪水池研究波浪与海洋结构物相互作用问题将是有效的手段。国内外众多学者对基于实验造波技术以及 [1]-[6] 数值模拟方法对数值造波理论进行了大量研究 。 近几年，由地震引起的海啸给人类造成了巨大的灾害，对于近海岸的各种海洋平台结构物也造成了很 大的破坏。而孤立波是一种典型的强非线性长波，可以用来描述海啸和风暴引起的巨浪的某些特征，在以 往的研究中常以孤立波来表示海啸来研究极限波浪与海洋结构物的相互作用。国内有许多学者对孤立波与 [7] 结构物的相互作用问题进行了研究，万德成、缪国平 应用VOF 方法研究了表面孤立波翻越直立方柱流动 [8] 的全过程；刘长根、陶建华 用基于时变雷诺方程和方程为基础的数学模型模拟了表面孤立波与近海不同 [9] 淹没深度水平圆柱的相互作用过程；谢伟松、陶建华 用Level Set 算法模拟表面孤立波与前台阶的相互作 用等等。 为了研究强非线性孤立波与大尺寸海洋结构物的相互作用，为海洋工程结构物的设计提供支持，本文 基于开源程序库OpenFOAM 数值求解Navier-Stokes 方程，采用VOF 方法捕捉自由面，开发三维粘性数值 波浪水池。在此基础上，以具有大直径圆柱结构的Spar 平台为研究对象，数值模拟了孤立波的生成以及传 1) 本文工作得到国家自然科学基金项目(Grant No., 海洋工程国家重点实验室研究基金 (Grant No. GKZD010053-11)和上海东方 学者人才计划基金(Grant No. 2008007)资助。 2) E-mail ：dcwan@sjtu.edu.cn 播过程中遇到Spar 平台立柱时波面的变化过程，重点分析了孤立波对Spar 平台立柱的冲击和爬高现象。 1 数学模型 1.1 控制方程 本文基于 OpenFOAM 开源程序库中求解两相流问题的求解器，将空气和水作为不可压缩粘性流体处 理，采用有限体积法(FVM)求解其控制方程： u 0 (1) (u) (uu) p (u) g Fs (2) t 其中： ，u ，p 分别为流体密度、