波浪作用下管道冲刷暴发计算的边界元方法

管道冲刷暴发临界条件是海底管道设计和运营的重要参数。Sumer等给出了3个Keulegan-Carpenter(KC)数的波浪作用下冲刷暴发的试验曲线。但对于任意KC数,Sumer等的方法并不准确。因此,建立波浪作用下冲刷暴发条件计算方法,具有重要意义。因管道与床面相交形成了复杂的几何形状,管道与床面的交叉点是数值奇点,因而管道冲刷暴发临界条件计算较为复杂。本研究采用边界元方法计算波浪运动及渗流压力。边界元可准确地拟合海底与管壁边界,方便地处理管道与海床交叉点的数值奇点问题,并可直接计算出奇点处的垂向压力梯度,准确地计算冲刷暴发的临界条件。本研究计算结果与实测值符合良好。     

粘性泥沙在盐水中冲刷和沉降特性的试验研究

本文根据絮凝电化学理论,探讨了淤积物浓度对粘性泥沙起动切应力的影响,获得起动切应力τc与淤积物浓度S的相关式τ_o=KS~(2·33) 该公式得到众多试验资料的证实。冲刷试验表明,具有均匀容重的底床在不同水流切力下,冲刷率随时间基本成直线变化,自养密实底床的冲刷率随时间呈指数曲线变化。试验结果还表明均匀容重底床的冲刷率与水流流速近似成四次方关系,论文运用了窦国仁博士从能量原理获得的床面泥沙动密实系数概念,获得冲刷率的表达式为 E=(K_0)/(C_0~2)-(γsγ)/(γs-γ)(V-V_c)V~3/gdω     

后台阶下游床面冲刷过程数值模拟研究

以床面瞬时剪应力作为泥沙起动及输运的水动力机制,建立了结构物周围复杂流场下床面局部冲刷的数学模型。并应用大涡模拟方法对后台阶下游三维湍流流动进行数值模拟,得到台阶下游床面瞬时剪应力的分布规律。为了确定床面瞬时剪应力与泥沙上扬通量的关系,先应用数学模型对不同模型参数下,冲刷开始后5分钟时台阶下游床面形状进行试算。通过试算与实验结果的比较,确定床面瞬时剪应力与泥沙上扬通量关系中需要的模型参数。进一步对冲刷开始后30分钟内台阶下游床面演化规律进行模拟,模拟结果与实验结果相符合。     

桩周局部冲刷三维数值模拟

用数值模拟手段研究桩周冲刷越来越引起人们的重视.基于OpenFOAM开源CFD软件,采用VOF方法处理自由水面,通过动网格技术考虑床面地形随时间的变化,构建起桩周局部冲刷的动态三维数学模型.其中,通过改进沙滑模型以修正当生成的剖面坡度超过泥沙休止角时造成的模拟失真,同时克服数值计算的不稳定性.计算结果和水槽实验资料对比分析,二者符合较好.     

河口海岸泥沙输移的分组数学模型研究

根据河口海岸水沙输移的特点,建立了一个新的二维分组数学模型,用来预测该区域的水沙输移过程。该模型耦合了水动力模块、泥沙输移模块和床面演变模块。其中水动力模块基于浅水方程组,综合考虑了柯氏力、床面切应力以及表面风应力的影响,引入干湿判断法处理动边界。泥沙输移模块首先将泥沙按照粒径分组,针对不同泥砂性质,对各组泥沙分别进行建模求解。床面演变模块基于质量守恒方程,实时更新床面高程以及床沙级配变化,并传递给水动力模块,更新底部边界。该模型被应用在了英国塞汶(Severn)河口,其预测的泥沙浓度和实测数据以及不分组的模型的预测结果进行了比较,结果显示,文中建立的分组模型预测的结果要明显好于不分组模型。     

植被对波浪作用下床面切应力影响的数值模拟分析

本文基于OpenFOAM建立三维波浪数值水槽,模拟计算植被水域波浪作用下的床面切应力,分析了入射波高、植被密度、植被淹没高度、水流对植被水域波浪作用下床面切应力的影响。结果表明：纯波时,由于植被的阻水作用,植被水域床面切应力沿程衰减,其衰减程度与入射波高、植被密度及植被淹没高度呈现正相关;与纯波时相比,在波浪和同向流共同作用下正向床面切应力幅值增大,负向床面切应力幅值减小;弱水流对植被水域床面切应力的大小及分布无明显影响;强水流时,床面切应力在植被水域先增大后逐渐减小并在植被水域后显著降低。     

单向流定床床面上水平管道周围流场与切应力分布特性

采用简化PIV方法,对悬空比G为0～1.6的单向流定床床面上水平圆柱绕流进行了系统的试验研究.通过对试验资料的分析,给出了圆柱绕流时的流场图和无量纲的床面、横管顶部和横管底部切应力的变化图,讨论了切应力随悬空比、管径、流速和水深的变化趋势.作者的研究成果有利于进一步研究海底管道的冲刷机理和冲刷坑的形成,进而有针对性的进行管道保护,减少管道破坏失效.     

床面冲淤计算机处理方法

冲淤计算是河流，河口，海岸等水域中分析，预测床面变化的基础，利用非结构网络自动剖分技术形成覆盖计算区域的三角网，能较好地适应边界和特殊地形，基于高程内插建立三角网数字高程模型，采用垂直三棱柱体作为基本单位计算不同时段冲淤量，从而提高计算精度，通过网络剖分，利用在计算区域内任意选择子域进行淤量的计算和分析。同时将生成的非结构网格直接有于二维床面冲淤数学模型，便于对每个节点的计算结果进行检证以提高模型的可靠性。     

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

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

波浪作用下直立堤前床面边界层内传质速度

本文研究立波作用下的层流、紊流边界层内的流速和传质速度,并采用谱法将传质速度推广到不规则波情况.根据实验与计算分析,研究了床面边界层内传质速度在堤前冲刷中的应用.研究结果表明,在传质速度的峰值附近形成底沙的冲刷坑;由于不规则波作用下堤前传质速度的峰值沿程递减,因此相应冲刷坑的深度沿程减少.     

动、定砂床底部剪切力对比试验研究

在河口海岸工程中,常常会面临岸滩冲蚀、岸线演变、航道淤积、建筑物底部淘刷等涉及泥沙起动和输运的问题,而水流对底床的剪切力是研究泥沙起动与输运的重要参数。本文利用自行设计的底床剪切力测量装置,在不同流速的水流中,分别在固定砂床(定床)和可移动砂床(动床)上进行了底部剪切力的直接测量;同时,根据试验中声学多普勒流速仪(ADV)测得的流速信息,采用湍流动能法对底床剪切力进行了估算。结果显示:当比例系数取值0.19时,估算出的底床剪切力与测量值吻合较好。对测量结果进行分析后发现,流速较小、砂粒未起动时,动、定砂床上的底部剪切力大致相同;在有砂粒起动的情况下,动床上的底部剪切力比定床上的大,相对差值最大约20%。因此当涉及底床剪切力的问题时,需要先确认床面形式,然后再进行分析研究。     

A Lagrangian model for boundary layer growth and bed shear stress in the swash zone

A new model for the boundary layer development and associated skin friction coefficients and shear stress within the swash zone is presented. The model is developed within a Lagrangian reference frame, following fluid trajectories, and can be applied to both laminar flow and smooth turbulent flow. The model is based on the momentum integral approach for steady, flat-plate boundary layers, with appropriate modifications to account for the unsteady flow regime and flow history. The model results are consistent with previous measurements of bed shear stress and skin friction coefficients within the swash zone. These indicate strong temporal and spatial variation throughout the swash cycle, and a clear distinction between the uprush and backwash phase. This variation has been previously attributed the unsteady flow regime and flow history effects, both of which are accounted for in the new model. Fluid particle trajectories and velocity are computed using the non-linear shallow water wave equations and the boundary layer growth across the entire swash zone is estimated. Predictions of the bed shear stress and skin friction coefficients agree reasonably well with direct bed shear stress measurements reported by Barnes et al. (Barnes, M.P., O’Donaghue, T., Alsina, J.M., Baldock, T.E., 2009. Direct bed shear stress measurements in bore-driven swash. Coastal Engineering 56 (8), 853–867) and, for a given flow velocity, give stresses which are consistent with the bias toward uprush sediment transport which has consistently been observed in measurements. The data and modelling suggest that the backwash boundary layer is initially laminar, which results in the late development of significant bed shear during the backwash, with a transition to a turbulent boundary layer later in the backwash. A new conceptual model for the boundary layer structure at the leading edge of the swash is proposed, which accounts for both the no-slip condition at the bed and the moving wet–dry interface. However, further development of the Lagrangian Boundary Layer Model is required in order to include bore-generated turbulence and to account for variable roughness and mobile beds.     

Numerical Simulation of Wave Scour Around A Large-Scale Circular Cylinder

A numerical model is developed for estimation of local scour around a large circular cylinder under vvave action. The model includes wave diffraction around structures, bed shear stress calculation inside the vvave boundary layer and topo-graphical change model. The vvave model is based on the improved Boussinesq equations for varying depth. The vvave boundary layer is calculaled by solving the integrated momentum equation over the boundary layer. The bed shear stress due to streaming, an important factor affecting the sediment transport around a large-scale cylinder, is calculated. The Lagrangian drift velocity is included in calculation of the suspended sediment transport rates. The model is implemented by a finite element method and the results from the present model, which agree well with experimental data, are com-pared vvith those from other methods.     

Sand transport under combined current and wave conditions: A semi-unsteady,practical model

For the general purposes of morphodynamic computations in coastal zones, simple formula-based models are usually employed to evaluate sediment transport. Sediment transport rates are computed as a function of the bottom shear stress or the near bed flow velocity and it is generally assumed that the sediment particles react immediately to changes in flow conditions. It has been recognized, through recent laboratory experiments in both rippled and plane bed sheet flow conditions that sediment reacts to the flow in a complex manner, involving non-steady processes resulting from memory and settling/entrainment delay effects. These processes may be important in the cross-shore direction, where sediment transport is mainly caused by the oscillatory motions induced by surface short gravity waves.The aim of the present work is to develop a semi-unsteady, practical model, to predict the total (bed load and suspended load) sediment transport rates in wave or combined wave-current flow conditions that are characteristic of the coastal zone. The unsteady effects are reproduced indirectly by taking into account the delayed settling of sediment particles. The net sediment transport rates are computed from the total bottom shear stress and the model takes into account the velocity and acceleration asymmetries of the waves as they propagate towards the shore.A comparison has been carried out between the computed net sediment transport rates with a large data set of experimental results for different flow conditions (wave-current flows, purely oscillatory flow, skewed waves and steady currents) in different regimes (plane bed and rippled bed) with fine, medium and coarse uniform sand. The numerical results obtained are reasonably accurate within a factor of 2. Based on this analysis, the limits and validity of the present formulation are discussed.     

Direct bed shear stress measurements in bore-driven swash

Direct measurements of bed shear in the swash zone are presented. The data were obtained using a shear plate in medium and large-scale laboratory bore-driven swash and cover a wide range of bed roughness. Data were obtained across the full width of the swash zone and are contrasted with data from the inner surf zone. Estimates of the flow velocities through the full swash cycle were obtained through numerical modelling and calibrated against measured velocity data. The measured stresses and calculated flow velocities were subsequently used to back-calculate instantaneous local skin friction coefficients using the quadratic drag law. The data show rapid temporal variation of the bed shear stress through the leading edge of the uprush, which is typically two–four times greater than the backwash shear stresses at corresponding flow velocity. The measurements indicate strong temporal variation in the skin friction coefficient, particularly in the backwash. The general behaviour of the skin friction coefficient with Reynolds number is consistent with classical theory for certain stages of the swash cycle. A spatial variation in skin friction coefficient is also identified, which is greatest across the surf-swash boundary and likely related to variations in local turbulent intensities. Skin friction coefficients during the uprush are approximately twice those in the backwash at corresponding Reynolds number and cross-shore location. It is suggested that this is a result of the no-slip condition at the tip leading to a continually developing leading edge and boundary layer, into which high velocity fluid and momentum are constantly injected from the flow behind and above the tip region. Finally, the measured stress data are used to determine the asymmetry and cross-shore variation in potential sediment transport predicted by three forms of sediment transport formulae.     

Mechanism of local scour around submarine pipelines based on numerical simulation of turbulence model

1 IntroductionSubmarine pipelines are important ocean engi-neeringequipm ents, especiallyfortheoiland gasin-dustries. M anyresearchershavebeendedicatedtotheproblem of local scour around the submarinepipelines, butmostoftheirwork arelaboratory testsand foc…     

Numerical study on vertical structures of undertow inside and outside the surf zone

Nearshore shoaling and breaking waves can drive a complex circulation system of wave-induced currents. In the cross-shore direction, the local vertical imbalance between the gradient of radiation stress and that of pressure due to the setup drives an offshore flow near the bottom, called ‘undertow’, which plays a significant role in the beach profile evolution and the structure stability in coastal regions. A 1DV undertow model was developed based on the relationship between the turbulent shear stress and t...     

Research on Measurement of Bed Shear Stress Under Wave?Current Interaction

The movement of sediment in estuary and on coast is directly restricted by the bed shear stress. Therefore, the research on the basic problem of sediment movement by the bed shear stress is an important way to research the theory of sediment movement. However, there is not a measuring and computing method to measure the bed shear stress under a complicated dynamic effect like wave and current. This paper describes the measurement and test research on the bed shear stress in a long launder of direct current by the new instrument named thermal shearometer based on micro-nanotechnology. As shown by the research results, the thermal shearometer has a high response frequency and strong stability. The measured results can reflect the basic change of the bed shear stress under wave and wave?current effect, and confirm that the method of measuring bed shear stress under wave?current effect with thermal shearometer is feasible. Meanwhile, a preliminary method to compute the shear stress compounded by wave?current is put forward according to the tested and measured results, and then a reference for further study on the basic theory of sediment movement under a complicated dynamic effect is provided.     

Research on Measurement of Bed Shear Stress Under Wave-Current Interaction

The movement of sediment in estuary and on coast is directly restricted by the bed shear stress. Therefore, the research on the basic problem of sediment movement by the bed shear stress is an important way to research the theory of sediment movement. However, there is not a measuring and computing method to measure the bed shear stress under a complicated dynamic effect like wave and current. This paper describes the measurement and test research on the bed shear stress in a long launder of direct current by the new instrument named thermal shearometer based on micro-nanotechnology. As shown by the research results, the thermal shearometer has a high response frequency and strong stability. The measured results can reflect the basic change of the bed shear stress under wave and wave?current effect, and confirm that the method of measuring bed shear stress under wave?current effect with thermal shearometer is feasible. Meanwhile, a preliminary method to compute the shear stress compounded by wave?current is put forward according to the tested and measured results, and then a reference for further study on the basic theory of sediment movement under a complicated dynamic effect is provided.