低雷诺数下圆柱强迫振动特性光滑粒子流体动力学模拟

基于光滑粒子流体动力学(SPH)方法,开发了能够准确描述水流作用下圆柱强迫振动特性的数学模型。通过引入适合于无网格粒子法的开边界算法,来模拟出入流边界条件,建立了具有造流功能的SPH数值水槽。圆柱及计算域的上下边界均采用修正的动力边界条件进行模拟。借助于粒子位移矫正和压力修正算法,避免了圆柱周围流体粒子压力大幅震荡以及结构下游区域出现空腔等非物理性现象。使用典型的圆柱绕流数据,验证了所建SPH模型的计算性能,研究了固定圆柱在低雷诺数情况下的尾涡脱落模式和升阻力变化规律。明确了低雷诺数下强迫振动圆柱在频率锁定以及非锁定区间内的升力变化规律,量化了升力与外界激励频率之间的关系。     

亚临界区雷诺数下圆柱绕流场电磁力控制数值研究

在亚临界区雷诺数下,采用脱体涡模拟方法对弱电解质中电磁力作用下圆柱绕流场及其升阻力特性进行了数值模拟与分析。结果表明,电磁力可以提高圆柱体边界层内的流体动能,延缓圆柱体近壁面流动分离,减弱绕流场中流向和展向大尺度涡的强度,减小圆柱体阻力及其升力脉动幅值;当电磁力作用参数大于某个临界值后,流动分离角消失,在圆柱体尾部产生射流现象,电磁力产生净推力作用,出现负阻力现象,而且升力脉动幅值显著减小且接近于零。     

各种间隙比下海底管线绕流二维数值模拟与实验研究

通过数值模拟方法与物理模型实验,对不同间隙比情况下的海底管线绕流特性进行研究。数值模拟采用有限差分法对标准k-ε湍流方程和控制方程进行离散,利用GMRES算法求解离散方程,模拟绕流流场。物理模型实验中采用超声测速仪ADV系统测取不同断面的流速分布情况,分别分析不同间隙比情况下的漩涡运动特性。同时分析各垂直断面不同间隙比时尾流流速变化情况,数值模拟结果与实验结果吻合较好,本文成果为进一步研究海底管线的防护措施打下基础。     

亚临界雷诺数下V型沟槽圆柱减阻数值模拟

为了深入研究亚临界雷诺数下光滑和沟槽圆柱绕流阻力及流场特性。采用SSTk-w模型对雷诺数Re分别为2×104,6×104,1×105和1.4×105的光滑圆柱和V型沟槽圆柱绕流进行三维流场数值模拟,首先以光滑圆柱绕流为算例与相关实验研究结果进行了对比,验证了SSTk-w方法的准确性,然后在不同雷诺数条件下,开展了V型沟槽圆柱绕流数值模拟研究;与光滑圆柱相比,V型沟槽圆柱具有较好的减阻效果,最大减阻率达28.4%;在相同雷诺数下,沟槽圆柱表面边界层分离点与光滑圆柱表面分离点相比更靠后。V型沟槽织构可有效改变尾流区旋涡脱落频率和尾流区湍流的发展,为隔水管减阻和提高隔水管柱系统的稳定性提供参考。     

张力腿平台绕流数值模拟分析

在一定来流条件下,张力腿平台(tension leg platform,简称TLP)的立柱后缘出现周期性的交替旋涡脱落,致使立柱受到垂直于来流方向的升力和平行于来流方向的阻力作用,导致TLP产生大幅度往复运动,显著增加平台结构和系泊系统的负载。目前,关于单柱、多柱结构绕流问题的研究较多,但对于TLP绕流特性的研究较少,机理尚存不明确的地方。为研究TLP的绕流力变化情况和流场特征,开展了数值模拟分析。利用计算流体动力学数值模拟软件,基于雷诺平均(RANS)法和分离涡模拟(DES)法对TLP绕流场进行仿真分析,揭示了TLP的绕流特性。结果表明,在3种来流角度和多个折合速度V_r下,TLP绕流的流体力系数存在明显差异,升力系数时域曲线呈现脉动性。TLP的上、下游立柱间存在明显的相互作用,影响了旋涡的形成与发展。TLP的旋涡脱落大多集中在平台固有频率附近,且在V_r=7,来流角度为0°时,升力系数频谱峰值最大,旋涡脱落集中。     

基于MPS方法模拟低雷诺数方柱绕流问题

柱体绕流问题是流体力学领域一个非常经典的问题。当流体流经柱体时,由于黏性的存在,会发生许多复杂的流动现象,如流动分离、涡旋周期性生成与脱落等,经常被作为标准验证算例。同时,柱体绕流广泛存在于实际工程中,并在一定工况下可能对工程产生巨大危害,因此对柱体绕流进行深入研究具有重要意义。研究中,拟将一种无网格类方法——半隐式移动粒子方法(moving particle semi-implicit method,简称MPS)引入到柱体绕流问题的数值研究中,并对不同雷诺数下二维方柱绕流问题进行数值模拟。首先,使用基于MPS方法自主开发的MLParticle-SJTU求解器,结合入口边界条件和出口边界条件,模拟了雷诺数Re分别为40、200和1 000时均匀来流条件下的方柱绕流。随后,将模拟的绕流结果与文献中试验和数值计算结果进行了对比,结果吻合较好,并且在雷诺数为200和1 000时,可以清晰地捕捉到方柱尾流中的卡门涡街现象,验证了MPS方法在柱体绕流问题模拟上的有效性和适用性。     

不同雷诺数下倾斜圆柱绕流三维数值模拟研究

基于计算流体力学(CFD)开源代码OpenFOAM开展了不同雷诺数(Re=100、1500和3900)和倾斜角度(-60°≤α≤60°)工况下倾斜圆柱绕流流场的三维数值模拟,研究倾斜圆柱绕流的三维瞬时及时均尾流流场、流线拓扑、升阻力系数与旋涡脱落频率随雷诺数及倾斜角度变化的规律,探讨在顺流向及逆流向情况下独立性原则对倾斜圆柱绕流的适用性。研究结果表明:随着圆柱倾角的增大,倾斜圆柱尾流产生较为明显的轴向流,尾流旋涡脱落受到明显抑制,细碎旋涡逐渐消失,尾流宽度随之减小;随着雷诺数的增大,圆柱尾流涡管发生明显的变形,展向掺混使得大量细碎旋涡产生,呈现出明显的三维特性。在不同雷诺数下,阻力系数均值、升力系数均方根及无量纲涡脱频率在一定倾角范围内符合独立性原则。     

近壁圆柱绕流水动力特性数值模拟与实验研究

通过数值模拟和物理模型实验,对距壁面一定高度的圆柱绕流水动力特性进行了研究。数值模拟采用有限体积法对标准k-ε模式方程进行离散,采用SIMPLE算法进行求解,模拟绕流流场。在物理模型实验中,将PVC圆管制作的实验模型安放在水槽内,在圆管的跨中沿表面周向均匀布置水下压力传感器,用于测量绕流圆柱体表面动水压力分布。通过改变Re数和间隙比来分析它们对近壁圆柱绕流水动力特性的影响。基于数值流动显示技术,给出了近壁绕流流场的尾流流态分析。通过数值结果与实验结果的对比,对近壁绕流圆柱体的升力系数及其表面动水压力分布进行了研究,对比结果显示了较好的一致性。     

雷诺数为3 900时三维圆柱绕流的大涡模拟

自从人们对层流的圆柱绕流现象有了系列研究及清楚的认识后,人们逐渐把目光投向湍流的圆柱绕流,但相关研究主要关注于模拟湍流方法的数值格式和精度问题,而忽略了对高雷诺下圆柱绕流流场本身的认识及规律的总结。基于开源代码Open FOAM的大涡模拟方法以连续方程和Navier-Stokes方程为控制方程,选用Smargorinsky模式为亚格子应力模型,采用有限体积法和一次预测两次修正的PISO算法,对Re=3 900时三维圆柱绕流问题进行了数值模拟研究,并着重分析了其尾流特征和性质。数值计算结果表明:大涡模拟方法可以模拟出细致的流场结构,该雷诺数下的圆柱绕流具有很强的三维及湍流效应,在圆柱后方约一倍直径的范围内存在回流区域,在靠近圆柱壁面的尾流区域的速度剖面呈"U"型,远离壁面的速度剖面呈"V"型。瞬时速度剖面始终围绕着时均速度的周围脉动,且距离圆柱越远瞬时速度场的脉动范围越大。     

后挖沟深度对深海海底管道屈曲影响数值分析

鉴于海底管道的服役水深越来越深,主要采用犁式挖沟机对预铺设于海床之上的海底管道采取后挖沟的方式将海底管道埋设于海床之下,以保护其免受不必要的损伤。针对后挖沟深度H是海底挖沟机的重要设计参数,也是影响管道悬跨的重要因素的问题,对SMD(UK)犁式挖沟机展开参数优化,确保作业过程中悬跨段管道在外部静水压力作用下,海底管道不会发生屈曲破坏。采用ABAQUS软件,分别建立了作业前和作业中两种工况下的悬跨模型,分析机械手对接触部分管道的损伤,结果显示,作业中的机械手对悬跨管道的损伤更大;同时,建立了作业中不同管径下,后挖沟深度对管道损伤的安全裕量关系曲线。进一步,结合作业中不同挖沟深度下的管跨段屈曲数值模型,对处于外部静水压力作用下的悬跨管的屈曲失效展开分析,结果显示,随着后挖沟深度的加大,不同管径下的悬跨段管道局部出现塑性压溃的临界压力值不断降低;管道外径的增大,降低了同一后挖沟深度下发生屈曲失效的压力值。最后,在后挖沟深度与外部静水压力组成的区域内,建立屈曲失效临界关系曲线,并划分出工作区和压溃区,为深海管道后挖沟埋管的施工提供工程参考。     

Numerical simulation of a partially buried pipeline in a permeable seabed subject to combined oscillatory flow and steady current

Hydrodynamic forces exerting on a pipeline partially buried in a permeable seabed subjected to combined oscillatory flow and steady current are investigated numerically. Two-dimensional Reynolds-Averaged Navier-Stokes equations with a k−ω turbulent model closure are solved to simulate the flow around the pipeline. The Laplace equation is solved to calculate the pore pressure below the seabed with the simulated seabed hydrodynamic pressure as boundary conditions. The numerical model is validated against the experimental data of a fully exposed pipeline resting on a plane boundary under various flow conditions. Then the flow with different embedment depths, steady current ratios and KC numbers is simulated. The amplitude of seepage velocity is much smaller than the amplitude of free stream velocity as expected. The normalized Morison inertia, drag and lift coefficients based on the corresponding force coefficients of a fully exposed pipeline are investigated. The normalized Morison force coefficients reduce almost linearly with the increase of embedment depth and that the KC only has minor effect on the normalized Morison coefficients. It is also found that the permeable seabed condition causes a slight increase on the inline force and has a little effect on the lift force, compared with corresponding conditions in an impermeable bed.     

Experimental Study on Free Spanning Submarine Pipeline Under Dynamic Excitation

Seismic load has a significant effect on the response of a free spanning submarine pipeline when the pipeline is con-structed in a seismically active region. The model experiment is performed on an underwater shaking tahle to simulate the response of submarine pipelines under dynamic input. In consideration of the effects of the terrestrial and submarine pipeline , water depth, support condition, distance from seabed, empty and full pipeline, and span on dynamic response, 120 groups of experiments are conducted. Affecting factors are analyzed and conclusions are drawn for reference. For the con-trol of dynamic response, the span of a submarine pipeline is by far more important than the other factors. Meanwhile, the rosponse difference between a submarine pipeline under sine excitation and that under random excitation exists in ex-periments.     

海底滑坡作用下滩海管道结构安全分析

海底滑坡作为常见的海洋地质灾害,对海洋油气工程安全产生巨大威胁。海床土体失稳引起滑坡体滑动,会对海底管道产生拖曳作用。基于计算流体动力学方法(CFD)建立海底滑坡体对管道作用的评估模型,采用H-B模型描述块状滑坡体并与试验比较验证,分析不同海床倾斜度滑坡对管道的作用并拟合表达式;研究了海底管道在滑坡作用下的力学响应,并采用极限状态方法开展海底滑坡作用下管道结构极限安全分析,探讨了管道埋地状态时的极限安全界限,建立滑坡作用下管道结构安全分析方法。研究表明:滑坡对管道作用力与海床倾角呈现正相关,而覆土层厚度对作用力影响较小;随着不排水抗剪强度的减小,允许的滑坡宽度和速度均增加,表明土体不排水抗剪强度与引起的拖曳力呈正相关;滑坡土体宽度对极限安全速度影响较大。     

海底滑坡对置于海床表面管线作用力的CFD模拟

海底管线是海洋工程中用于传输原油和天然气等的重要通道,通常放置在海床表面或处于悬跨状态。本文采用计算流体动力学CFD法模拟了不同冲击角度下海底滑坡对置于海床表面的海底管线的作用,得到了管线所受的轴向荷载和法向荷载与滑坡冲击角度之间的关系。同时分析了沿冲击方向管线截面形状与管线所受阻力之间的关系。对已有研究进行拓展延伸,丰富了不同工况下轴向阻力系数和法向阻力系数的计算成果,得出了海底滑坡对置于海床表面管线冲击力的计算公式。     

Submarine debris flow impact on pipelines — Part I: Experimental investigation

Estimating the impact forces exerted by a submarine debris flow on a pipeline is a challenge, and there is room for considerably more work to advance the state of the art. To this end, an experimental program was performed to investigate the impact on two pipeline installation scenarios: 1) suspended pipeline and 2) laid-on-seafloor pipeline. The results and observations from the experimental investigation are discussed. The definition of Reynolds number was modified for non-Newtonian fluids and an ad hoc method was developed to estimate the drag force exerted by an impact perpendicular to the pipe axis. The method may be used in prototype situations to estimate the drag force from submarine debris flow impact on pipelines. The experimental program was complemented by Computational Fluid Dynamics (CFD) analyses, the details of which are discussed in the accompanying paper.     

Parametric study of the wave-induced residual liquefaction around an embedded pipeline

Wave-induced liquefaction in a porous seabed around submarine pipeline may cause catastrophic consequences such as large horizontal displacements of pipelines on the seabed, sinking or floatation of buried pipelines. Most previous studies in relation to the wave and seabed interactions with embedded pipeline dealt with the wave-induced instaneous seabed response and possible resulting momentary liquefaction (where the soil is liquefied instantaneously during the passage of a wave trough), using theory of poro-elasticity. Studies for the interactions between a buried pipeline and a soil undergoing build-up of pore pressure and residual liquefaction have been comparatively rare. In this paper, this complicated process was investigated by using a new developed integrated numerical model with RANS (Reynolds averaged Navier–Stokes) equations used for governing the incompressible flow in the wave field and Biot consolidation equations used for linking the solid–pore fluid interactions in a porous seabed with embedded pipeline. Regarding the wave-induced residual soil response, a two-dimensional poro-elastoplastic solution with the new definition of the source term was developed, where the pre-consolidation analysis of seabed foundation under gravitational forces including the body forces of a pipeline was incorporated. The proposed numerical model was verified with laboratory experiment to demonstrate its accuracy and effectiveness. The numerical results indicate that residual liquefaction is more likely to occur in the vicinity of the pipeline compared to that in the far-field. The inclusion of body forces of a pipeline in the pre-consolidation analysis of seabed foundation significantly affects the potential for residual liquefaction in the vicinity of the pipeline, especially for a shallow-embedded case. Parametric studies reveal that the gradients of maximum liquefaction depth with various wave and soil characteristics become steeper as pipeline burial depth decreases.     

波浪作用下孔隙海床-管线动力相互作用分析

波浪作用下海床中的孔隙水压力与有效应力是影响海底管线稳定性的主要因素。然而,在目前的海床响应分析中一般将管线假定为刚性,并不能合理地考虑海床与管线之间的相互作用效应,同时也没有考虑土体和管线加速度对海床动力响应的惯性影响,从而无法确定由此所引起的管线内应力。为此考虑管线的柔性,分别采用饱和孔隙介质的Biot动力固结理论和弹性动力学理论列出了海床与管线的控制方程,进而采用摩擦接触理论考虑海床与管线之间的相互作用效应,基于有限元方法建立了海床-管线相互作用的计算模型及其数值算法。通过变动参数对比计算讨论了管线几何尺寸、海床土性参数对波浪所引起的管线周围海床孔隙水压力和管线内应力的影响。     

基于OLGA的海底管道水合物的生成和影响因素分析

海底管道是海洋油气输送的重要纽带。为解决海底管道面临的水合物生成和堵塞问题,文章结合海底多相流管道水合物生成的数学模型,采用OLGA对海底管道不同含水率、气油比和流量下水合物的生成情况进行数值模拟。研究结果表明：在某海底管道的工艺参数条件下,水合物生成风险随着含水率和流量的增大而降低,随着气油比的增大而增加;海底水平管路和立管都有可能生成水合物,尤其立管常是水合物最大生成量的位置;模拟结果可为海底管道水合物的防治和保障海底管道的安全运行提供参考。     

Interaction between submarine landslides and suspended pipelines with a streamlined contour

Recently, the security and stability of submarine pipelines have attracted much attention in ocean engineering. In this paper, pipelines with a streamlined contour (wedge, airfoil, double-ellipse, and arc-angle hexagon) are designed in hopes of defending against the impact of submarine landslides, and the computational fluid dynamics (CFD) approach is used to investigate the interaction between submarine landslides and streamlined pipelines. The results show that the peak interactional force is more representative of the hazard level of pipelines imposed by submarine landslides. It is also found that the streamlined pipelines possess a significant advantage in reducing the drag force and lift force of landslide–pipeline interaction with a maximum lessening percentage of 66.32 and 40.17%, compared with a conventional circular pipeline. In addition, the influence of applying streamlined pipelines to engineering is briefly discussed, and the empirical equation for estimating the drag force and lift force of streamlined pipelines induced by landslides is recommended based on the numerical test results.     

Numerical analysis of the impact forces exerted by submarine landslides on pipelines

Submarine pipelines that transport crude oil and natural gas are often in a complex marine geological environment and may become unstable and fail upon impact by submarine landslides. Previous research has mostly focused on the impact forces exerted by submarine landslides on suspended pipelines, but the impact of submarine landslides on pipelines laid on the seafloor at various impact angles, θ, have been relatively infrequently discussed, and the effects of suspended height, H, on the impact forces exerted by submarine landslides on pipelines have not been thoroughly investigated. In this study, based on the Herschel–Bulkley model, the impact forces exerted by a submarine landslide on laid-on or suspended pipelines at various impact angles θ were simulated using the computational fluid dynamics (CFD) approach. Equations for calculating the axial and normal drag coefficients of a submarine pipeline were proposed. The CFD numerical simulation results were rearranged based on the soil mechanics approach. By comparing the parameters, an essentially corresponding relationship was found between the soil mechanics and CFD approaches when the equations were used to calculate the impact forces exerted by a submarine landslide on a pipeline. In addition, a semi-analytical expression for the failure envelope was provided. Furthermore, the effects of H on the forces on a pipeline were discussed, and an equation for calculating the acting forces on a pipeline along the flow direction of a submarine landslide that comprehensively accounts for the effects of θ and H was proposed. The lift force was discussed preliminarily and the results provide a basis for further investigation. The achievement of this study is applicable for selecting locations of submarine pipeline routes and for designing submarine pipelines.