绕海底管线悬空段肩部的三维湍流结构的数值的研究.pdfVIP

中 国 造 船 V0143 a 43卷增刊 Speci 2002年10月 SHIPBUILDINGOFCHINA oct，2002 文章编号：1000—4882(2002)S一00330—06 绕海底管线悬空段肩部的三维湍流 结构的数值研究 陈兵1李玉成1程亮2 (1大连理工大学2．西澳大利亚大学) 捅 要 用三维有限元方法模拟了恒定来流绕悬空的海底管线与海床相交处，即管线悬空段肩部的三维流动情况。基于管 线直径和来流流速的Re数为5xl 理，由涡量幅值的三维等值面可见，绕悬空段肩部在管线上游形成了一条涡管，它的一端绕过半埋于海床的管线段的 的特征值定义的涡核可进一步显示三维流动结构的细节。本文的研究为利用数值模型进一步研究海底管线周围局部冲 刷沿管线轴向发展打下了初步的基础。 关键词：海底管线；悬空段肩部；三维流动结构；数值模拟 由于海底管线在现代海上石油天然气开发中的广泛应用，对其受力和局部冲刷等关系到管线稳定 性的研究也得到越来越多的重视。当海底管线铺设在可侵蚀海床上时，由于海床的地形变化或管线F 方的局部冲刷的发展，很可能形成局部管线段的悬空状态。预估由管线悬空段的长度及其由于局部冲 刷引起的冲坑增长速率是管线稳定性设计中需要考虑的重要条件。显然，绕悬空的海底管线Lj海床相 交处，即管线肩部的三维流动情况，对理解管线肩部周围的冲刷机理有至关重要的作用。Sumer和 管线悬空段肩部上游形成的螺旋形涡管对三维冲坑的发展也有很大影响。然而，绕管线悬空段肩部的 三维流态究竟如何还没有清晰确定的了解(Sumer和FredsOe，1999)。 2000b及2001)、BrOrs1999。这些研究的结果令人鼓舞，甚至对涡脱落引起的随时间发展的冲刷过程 的模拟也相当的成功。尽管这些数值研究都局限于二维情况，但其数学模型却可以扩展到三维，且本 将原有工作扩展到湍流情况。 建立的模拟绕管线悬空段肩部的三维流动的数值模型，并将其进一步发展应用于湍流情况，研究管线 肩部附近的三维湍流结构和底床剪切应力分布。希望本文的研究能为进一步发展模拟海底管线周围的 三维局部冲刷模型打下基础。 (二)数学模型 1．基本方程 亚格子(SGS)应力模型来模拟湍流应力。本文用分裂步法求解用有限元守间离散的控制方程。即把 每个时间步长分成三步，在头两步中先基于初始的速度和压力显式地估算速度场，再通过求解压力 43卷增刊 陈兵等：绕海底管线悬空段肩部的三维湍流结构的数值研究 等，1995)求解器用来迭代求解压力Poisson方程。 2．初边值条件 湍流条件下，假定来流速度沿垂向成指数线分布： 1 u：f上]_ U。16J 式中U。为来流最大速度，Y是到固体壁面的距离，踞入口处边界层厚度。 在出口边界上，压力设为常数且保持不变，速度的各分量的法向导数为零。在海床和管线表 件，即仅设速度的法向导数为零。 初始时刻流体静止，速度场和压力场都为零，然后在入口和出口边界上施加一个突然的起动流速 (恒定来流速度)。 (三)流动和地形参数及网格系统 上述数学模型用来模拟恒定来流绕横跨斜坡的水平管线的流动。图1是管线和海底地形的示意 图，管线半埋段长度B为2D，D为管线直径。上游入口边界和计算域顶部边界到管线中心的距离是 数。本文研究Re=5×10’的湍流流动。入口处边界层厚度子设为O．5D。 在三维有限元模型中使用8节点的长方体单元。图2为本研究采用的网格系统，为了清楚地表示 出网格的形状，图中只给出了网格在计算域边界上的投影。靠近海底和管线处的网格经过了适当的加 个单元。