海底管线附近底床侵蚀平衡剖面的湍流数值模拟.pdfVIP

45卷增刊 中吲造船 V01．45 Special 2004年lO月 SHIPBUILDINGOFCHINA Oct～2004 文章编号：11300—4882(2004)03—0209—8 海底管线附近底床侵蚀平衡剖面的 湍流数值模拟 李玉成吕林 (大连理工大学海岸与近海工程国家重点实验室) 摘 要 本文建立了模拟海底管线附近底床侵蚀平衡剖面的垂向二维数学模型。对其中的三个关键性问题—— 复杂边界下的湍流场计算、泥沙底床侵蚀动边界平衡剖面的确定——分别采用有限元数值方法求解改进的 膏1双方程湍流模式、预网格剖分的边界拟合技术和基于壁面律计算剪切应力结合泥沙临界起动切应力准 则进行判别的方法进行处理。为检验数学模型的可靠性，文中将数值模拟结果与实验资料进行了对比，结 果显示：本文建立的数学模型能够比较准确的模拟出侵蚀平衡剖面的形态并揭示了侵蚀发生的机制。同时 也表明经Shih等改进的七1湍流模式适合于海底管线附近湍流场的数值计算；预网格边界拟合技术是处 理复连通域下活动边界问题的有效方法；管线附近泥沙底床的侵蚀平衡剖面主要受底床附近切应力的控 制，临界起动切应力判别准则在该问题中的应用是合适。 关键词：海底管线；数值模拟：平衡剖面；湍流；侵蚀 1 前言 海底管线作为一种相对高效、经济的输运方式在海洋工程中得到了J一泛的应用，特别是对近海 油气的开发和输运有着重要的意义。针对海底管线附近底床在水流和管线的共同作用下可能会发生 局部侵蚀，进而影响管线安全这一问题，本文基于湍流理论和有限元数值方{．，}：建立了模拟底床局部 冲刷平衡剖面的数学模型，其中流场的正确模拟和冲刷平衡剖面的确定是两个最重要的问题。 活动边界问题是一个广泛存在的自然现象，也是各种生产活动中经常会遇到的一个难题，因此 长期以来引起了人们的普遍关注，关于这一领域早期研究成果的总结可见于Crattk的著作”’。目前使 用最普遍的三类基本方法是：主要针对流动问题的欧拉方法，如VOF、MAC等；拉格朗H类方法则 被广泛用结构分析问题中；任意拉格朗日一欧拉方法(ALE)是近年发展起来的一种新型方法，是 上述两种方法的混合形式，兼备二者的优点。但是这些方法都有各自的缺陷，难以直接用于追踪海 底管线附近底床冲刷活动边界，主要集中在：对复杂边界的适应性、难以定义有效地体积函数、数 值稳定性以及计算量巨大等。另外，在泥沙输运、底床变形问题中，基于流场数值计算结果结合沉 积物浓度扩散方程和底床变形连续方程进行预测分析也是一种有效的手段，但是对于推移质输沙问 题扩散方程失效；同时，在本问题中由于管线的存在使流场成为一个复连通区域，此时采用深度平 均、引入水深函数等常规方法都将遇到诸多难题。因此，要对管线附近局部冲刷平衡剖面进行模拟、 预测，必须发展新型有效的底床冲刷变形追踪技术。 210 中 同 造 船 学术论文 自20世纪80年代以来，一些学者已经对海底管线附近底床侵蚀问题开展了数值模拟工作。Van Beek和Wind【21采用标准七吨湍流模式结合沉积物浓度扩散和底床变形连续方程针对悬移质输沙进 行了数值研究，虽然数值模拟结果在总体趋势上与实测资料比较符合，但是管线后方底床形态与实 Li【4】都认为造成这种差异的根本原因在于原始的七吨 测资料存在本质上的差异。Sumer等口1和Fangjun Li口1分别采用势流理论和大涡模拟的方法 湍流模型不能捕捉到管线后方的涡泄现象。最近，Fangjun 对之进行了数值计算，但是前者效果不理想，后者计算量的要求较大。