海底管道冲刷势流模型的改进

对Li&Cheng的势流模型在数值方法、床面平衡条件及冲刷床面调整技术做出了3方面的改进。采用边界元法代替差分法求解Laplace方程,前者可以准确地拟合地形与管道边界,因而可以准确反映固壁边界对流态的影响,此外,还具有数据准备简单,降低计算维数,计算速度快等优点。Li&Cheng模型以床面水流切应力等于泥沙起动切应力,τb=τc,作为床面平衡条件,这只适用于清水冲刷。以沿程输沙相同作为平衡剖面条件,理论上,将模型推广到了动床冲刷。此外,为了提高模型的收敛性,提出了最速下降法与牛顿迭代法相结合的床面调整技术。采用实验资料对模型进行了检验,计算的冲刷深度与实验结果符合良好。     

海流作用下海底管道局部冲刷数值分析

在任意拉格朗日-欧拉参考坐标系下,采用基于雷诺平均的Navier-Stokes方程组(RANS)、流线迎风有限元方法、泥沙输运模型以及底床变形方程,对海流引起的海底管道局部冲刷进行了数值模拟。着重讨论了均匀来流流速和海底管道直径对局部冲刷发展过程及平衡冲刷深度的影响作用。数值结果表明,在冲刷的初始阶段,冲刷深度随时间迅速增加,之后缓慢逼近极限平衡深度;在管径一定的情况下,管道附近的局部平衡冲刷深度与流速大致呈线性关系;当流速超过某一临界区域后,最大平衡冲刷深度出现的位置并不在管道正下方,而是随流速的增加向管道下游方向移动;另外,管道直径也会对平衡冲刷深度产生比较明显的影响,在相同流速下,平衡冲刷深度大致随管径呈线性增大。在本文的计算范围内,海底管道的相对局部平衡冲刷深度基本随雷诺数线性增加,但流速对冲刷深度的影响作用要比管径的影响作用更为明显。     

基于高阶边界元的三维数值波浪港池

初步建立了一个基于高阶边界元的三维数值波浪港池,港池具有造波和消波功能。采用高阶边界元16节点四边形单元和基于二阶显式泰勒展开的混合欧拉-拉格朗日时间步进求解带自由表面的完全非线性势流方程。模型中对于影响数值精度的问题作了细致的处理。数值计算结果表明本港池可以用来模拟非线性波浪的传播,具有很高的数值精度和稳定性。     

往复流作用下海底管道局部冲刷机制数值模拟

基于开源的计算流体力学模式REEF3D,建立了海底管道局部冲刷水槽数值模型,在验证单向流实验结果的基础上,进一步对往复流作用下的海底管道局部冲刷机制进行了研究,并作对比分析。研究表明,总体上往复流对管道所在海床的局部冲刷规模比单向流弱。当流向改变后,原先下游的堆积区转变成上游,优先受到冲刷,并填充到管道下方的冲刷坑,同时使水流在管道下方的作用减弱。这种回填过程使短周期下的往复流作用需要更长的冲淤平衡时间。在回填与冲刷的共同作用下,上下游的冲刷坑坡度会因流向变化而变化。     

舟山册子岛—宁波镇海海底原油管道冲刷的影响因素

局部冲刷引起裸露甚至悬空是造成管道失稳的主要原因之一。舟山册子岛—宁波镇海海底原油管道自2006年运行以来,其赋存状态由全线掩埋逐渐演变为局部裸露甚至悬空的严峻局面,并且有继续恶化的趋势。基于管道2006—2012年历年检测资料,从多个层面分析造成海底管道冲刷的原因,认为造成管道冲刷的基本条件是研究区的自然冲淤规律,极端天气条件、镇海侧的大型围垦工程、舟山连岛大桥建设、长江来沙减少等其他因素的叠加也不可小视。     

基于CFD-DEM方法对海底管道床面冲刷的研究

采用计算流体力学—离散元耦合方法（CFD-DEM）模拟海底管道床面的冲刷过程。经过模型验证,该方法的计算结果与前人的研究具有较好的一致性,证明其可以应用于海底管道周围的冲刷模拟计算。冲刷初期的结果增强了目前对启动阶段粒子运动机理的理解,即管前后压力梯度造成的渗流作用导致粒子运动。对完整冲刷过程的模拟中,发现冲刷分为冲刷启动阶段、间隙冲刷阶段和尾迹冲刷阶段。间隙冲刷阶段管道下方粒子具有较大速度,冲刷坑快速向下方发展。进入尾迹冲刷阶段后,管道后方出现周期性脱落的涡旋,沙丘上的粒子速度更大。同时利用DEM更具直观性的独特优势,首次得到了14个典型位置处颗粒的运动轨迹和运动速度,对于理解冲刷过程中粒子的运动情况具有较大帮助。     

基于高阶边界元的三维数值波浪港池--波浪破碎的模拟

在势流理论的框架内,采用高阶边界元方法和混合欧拉-拉格朗日法,实现了对三维波浪破碎过程的数值模拟.数值模型使用可调节时间步长的基于二阶显式泰勒展开的混合欧拉-拉格郎日时间步进来求解自由表面的演化过程.在所使用的边界元方法中,采用16节点三次滑移四边形单元来表示,这种单元在单元内具有高阶的精度同时在单元之间具有良好的连续性.给出了孤立波的传播和周期性非线性波浪沿缓坡传播的计算结果,表明数值模型具有良好的稳定性.     

基于IDDES方法的冲刷地形海底管道绕流数值模拟

近年来,海底管道对于近海油气工业中碳氢化合物的运输十分重要,预测海底管道周围水动力特征对于研究管道损坏机理具有重要意义。基于Spalart-Allmaras改进延迟分离涡模拟(improved detached eddy simulation,简称IDDES)湍流模型,利用开源OpenFOAM求解器计算单管冲刷过程中4个特征阶段管道周围水动力特征。在数值模拟中先计算了无圆管时均匀来流充分发展后的流场,依据该流场中进口流速分布设定有管道时的边界层及进口边界条件。通过将最终模拟结果与前人试验数据及数值模拟结果进行对比,表明了该模型计算结果相对于k-ε模型及LES模型更加贴合试验结果,验证了定义边界层方法的可靠性及IDDES模型在预测管道周围流场的适用性。模型计算结果显示,随着冲刷坑深度的逐渐增加,管道下游回流区会经历一个大幅缩减又逐渐扩大的过程,且逐渐关于圆管中心所在水平线对称;通过对不同冲刷阶段尾涡形态的分析表明,冲刷前期并没有涡脱产生,中期已经有产生涡脱的趋势,后期出现了明显的且较为稳定的涡脱。     

波浪问题中唯一可解的高阶边界元方法

本文就波浪与结构物相互作用问题，提出了一个适用于高阶边界元应用的新的积分方程，并利用修改积分区域的方法得了适用于本积分方程的不规则频率消除方法。最后，通过数值计算对附加区域的选择、单元的离散做了研究     

册子岛—镇海海底原油管道冲刷稳定性研究

根据对册子岛-镇海段海底原油管道2006-2011年连续多年的监测,结合区域地理位置、水文泥沙、地形地貌、海洋开发活动、极端气象天气等因素,提出从海洋区域综合环境的角度研究海底管道所在海床的冲淤稳定性.管道所在区域海床总体表现为冲刷,2005年完成埋设的海管在2008年出现了多处裸露,到2010和2011年甚至有多个区段出现悬空.分析认为,舟山跨海大桥尤其是金塘大桥的建设、宁波澉浦新泓口围垦工程的实施、极端气象天气或短期灾变因子等对研究区海床的冲淤稳定产生的影响较大,也是导致近年来海管出现危险状况的主要原因.     

Investigation on the Local Scour Beneath Piggyback Pipelines Under Clear-Water Conditions

A piggyback pipeline consists of two pipes such that the secondary line rides on the main pipe with a fixed distance between two pipes in length. The novel strategy is utilized in offshore areas instead of a single flow line. In this regard, there are only a handful of experimental and numerical studies investigating the effect of scour below a piggyback pipeline under steady current. Hence, this study focuses on examining the influential factors on scouring due to steady current including the pipe diameter and the gap between pipes through numerical simulations and experimental tests. Accordingly, at the first phase of the research, a single pipe was established and tested in laboratory to compare the results with those of an empirical equation. After finishing experimental verifications,piggyback pipelines were also assembled to study the scouring under steady current conditions. It was concluded that by increasing the gap distance between the pipes, the maximum scour depth decreases; however, an increase in the small pipe's diameter results in a larger maximum scour depth. Secondly, numerical simulations were carried out using the FLOW-3D software which was found to be a suitable tool for the numerical investigation of this study.Finally, the numerical results have been compared with the corresponding experimental data and a relatively good agreement was achieved between them.     

Local Scour Around Piles Under Wave Action

The model tests are performed with regular waves, and the effect of wave height, wave period, water depth, sediment size and pile diameter is evaluated. The shape and size of local scour around piles are studied. There are three typical scour patterns due to wave action. It is found that a relationship exists between the erosion depth and the wave number. An empirical formula of the maximum local scour is thus derived.     

Numerical Simulation of Local Scour Around A Large Circular Cylinder Under Wave Action

A horizontal two- dimensional numerical model is developed for estimation of sediment transport and sea bed change around a large circular cylinder under wave action. The wave model is based on an elliptic mild slope equation. The wave-induced current by the gradient of radiation stress is considered and a depth integrated shallow water equation is applied to the calculation of the current. The mass transport velocity and the bed shear stress due to streaming are considered, which are important factors affecting the sediment transport around a structure due to waves, especially in reflective areas. Wave-current interaction is taken into account in the model for computing the bed shear stress. The model is implemented by a finite element method. The results of this model are compared with those from other methods and agree well with experimental data.     

Three-Dimensional Boundary Element Method Applied to Nonlinear Wave Transformation

For higher accuracy in simulating the transformation of three dimensional waves,in consid-eration of the advantages of constant panels and linear elements,a combined boundary elements is appliedin this research.The method can be used to remove the transverse vibration due to the accumulation ofcomputational errors.A combined boundary condition of sponge layer and Sommerfeld radiation condi-tion is used to remove the reflected waves from the computing domain.By following the water particle onthe water surface,the third order Stokes wave transform is simulated by the numerical wave flume tech-nique.The computed results are in good agreement with theoretical ones.     

The Boundary Element Method for Pile-Soil Interaction

In this paper, several mathmatical models for the pile- soil interaction are outlined. The Boundary Element Method is one of the very effective methods for the reasonable models of elasticity and elastoplasticity. The major of this paper is concerned with the Boundary Element Method for the pile-soil interaction, including general methods and calculating formulation of static and dynamic analysis of the pile and pile groups. Some results of analysis are also given.     

四阶波动方程的有限体积元法

建立一类四阶色散耗散方程全离散有限体积元格式,并给出格式的收敛性和误差估计,同时还给出数值结果。     

Simplified Method for Calculating Standing Wave Pressure on Vertical Breakwater

Through numerical modeling, a kind of simplified calculating method for standing wave pressure on vertical face breakwater have been put forward. Not only the formulas proposed in this paper are simple in form and very easy in use, but also they possess continuity on the full range of standing wave. And more, the precision requiremennts of calculation can be satisfied to a certain extent in engineering practice.     

海洋平台桩基周围冲刷过程及冲刷机理分析

对埕岛油田典型平台周边多年的水深地形资料进行分析对比,研究平台桩基周围冲刷过程,探讨桩柱周围形成冲刷坑的冲刷深度、几何形态和分布特征,分析其冲刷机理。结果表明,平台建成初期桩基冲刷剧烈,周围海底形成以平台为中心的盆状地形,经历近1年的时间逐步达到冲淤平衡。冲刷坑形态特征主要由水动力条件控制。     

Characteristics of Parameters for Local Scour Depth around Submarine Pipelines in Waves

The purpose of this study is to investigate the characteristics of parameters for local scour depth around submarine pipes due to waves. Because it is impossible to consider all the factors that influence the scour, it is reasonable to investigate the correlations between parameters and scour. The experimental works about scour were performed for the conditions of 40, 60, 90 m pipe diameter, various wave periods, and wave height in the case of horizontal bed and 1/10 slope bed, respectively. Analyzing the results of experiments, the Reynolds number, Shields parameter, Keulegan-Carpenter number, Ursell number, and Modified Ursell number were estimated. The correlations between the relative scour depth, which is the maximum equilibrium local scour depth divided by the pipe diameter, and 5 parameters were analyzed. It was shown that there was hardly any correlation with the Reynolds number. In the case of a horizontal bed, the Keulegan-Carpent number had the highest correlation, but in the case of the slope bed the correlation was greatly reduced. The modified Ursell number showed a high correlation regardless of the type of bed.     

Scaled Boundary Finite Element Analysis of Wave Passing A Submerged Breakwater

The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique combining the advantage of the finite element method （FEM） and the boundary element method （BEM） with its unique properties. In this paper, the SBFEM is used for computing wave passing submerged breakwaters, and the reflection coeffcient and transmission coefficient are given for the case of wave passing by a rectangular submerged breakwater, a rigid submerged barrier breakwater and a trapezium submerged breakwater in a constant water depth. The results are compared with the analytical solution and experimental results. Good agreement is obtained. Through comparison with the results using the dual boundary element method （DBEM）, it is found that the SBFEM can obtain higher accuracy with fewer elements. Many submerged breakwaters with different dimensions are computed by the SBFEM, and the changing character of the reflection coeffcient and the transmission coefficient are given in the current study.