波浪及波流边界层泥沙起动规律

基于波浪边界层理论及单向流泥沙起动Shields曲线,推证出波浪泥沙起动Shields曲线;基于波流边界层理论,提出表述波流边界层动力特征的波流比因子X及非线性作用因子Y,并建立了Y与X的相关关系;在此基础之上,结合单向流及波浪泥沙起动Shields曲线,推证出波流共同作用下泥沙起动Shields曲线。结果表明:波浪泥沙起动Shields曲线在层流区与单向流光滑紊流区曲线保持一致,粗糙紊流区与单向流粗糙紊流区曲线保持一致,过渡区线型为折线,由层流区及粗糙紊流区曲线延长交汇获得;X及Y能够合理地表征波流边界层动力对比特征及非线性作用特征;波流泥沙起动Shields曲线介于波浪及单向流泥沙起动Shields曲线之间,随着波流比因子X的不同,依据非线性作用因子Y,自动在波浪及单向流泥沙起动Shields曲线之间非线性过渡。建立的波流泥沙起动Shields曲线与试验结果吻合较好,且能够概括单向流、波浪及波流等不同动力及细沙、粗沙等不同粒径的泥沙起动条件。     

波浪及波流边界层泥沙起动规律

基于波浪边界层理论及单向流泥沙起动Shields曲线,推证出波浪泥沙起动Shields曲线;基于波流边界层理论,提出表述波流边界层动力特征的波流比因子X及非线性作用因子Y,并建立了Y与X的相关关系;在此基础之上,结合单向流及波浪泥沙起动Shields曲线,推证出波流共同作用下泥沙起动Shields曲线。结果表明:波浪泥沙起动Shields曲线在层流区与单向流光滑紊流区曲线保持一致,粗糙紊流区与单向流粗糙紊流区曲线保持一致,过渡区线型为折线,由层流区及粗糙紊流区曲线延长交汇获得;X及Y能够合理地表征波流边界层动力对比特征及非线性作用特征;波流泥沙起动Shields曲线介于波浪及单向流泥沙起动Shields曲线之间,随着波流比因子X的不同,依据非线性作用因子Y,自动在波浪及单向流泥沙起动Shields曲线之间非线性过渡。建立的波流泥沙起动Shields曲线与试验结果吻合较好,且能够概括单向流、波浪及波流等不同动力及细沙、粗沙等不同粒径的泥沙起动条件。     

层流区Shields曲线分布规律试验

为探明层流条件下Shields曲线的分布规律,运用电荷耦合组件(Charge Coupled Device,CCD)成像技术结合激光扫描及计算机图像处理技术,在水和甘油混合液中对无黏性均匀玻璃颗粒的起动拖曳力进行测量,并根据层流绕流理论建立了无黏性均匀颗粒的滚动起动模型,提出层流条件下无黏性均匀颗粒的量纲一起动拖曳力公式.试验观测发现在层流范围内希尔兹数呈规则的带状分布,流体作用引起颗粒床表面粗化并导致床面颗粒突起减小,使颗粒起动拖曳力增大1倍以上.结果表明颗粒床表面结构性状直接影响着颗粒起动拖曳力的大小,在层流区Shields曲线具有带状分布特性.     

淤泥质浅滩泥沙临界起动切应力剖面确定

为了确定淤泥质浅滩泥沙的临界起动切应力垂线剖面,采用音叉密度计在淤泥质连云港徐圩浅滩进行了定点密度垂线分布测量,并针对该海域的泥沙利用长水槽和环形槽开展了泥沙起动室内试验。淤泥密度现场结果表明,浅滩泥沙密度与深度满足对数型关系;室内试验得出密度1 050~1 400 kg/m³的泥沙临界起动切应力值为0.1~1.0 Pa,泥沙临界切应力与密度呈指数关系;进而确定了临界起动切应力与深度的关系即临界起动切应力剖面,该剖面关系式可供数学模型模拟淤泥质浅滩的泥沙起动过程参考。     

磨料两相圆形射流紊流方程的建立及其求解

根据一定合理假设从单相圆形层流边界方程的建立和求解入手，分别建立了两相圆形层流射流边界层方程，单相圆形紊流边界层方程和两相圆形紊流边界层方程并给出了各种情况下的流速分布等参数的解答，为研究磨料两相射流对材料的打击和切割作用奠定了理论基础。     

过渡粗糙平板紊流边界层统计特性试验研究

本文应用偏振差动式激光测速仪对过渡粗糙平板紊流边界层的统计特性进行了研究,得到紊动强度、紊动流速的三阶矩、四阶矩、概率密度分布形式、功率谱、自相关系数等特征.分析结果表明,可以根据高阶矩的符号对紊流边界层沿垂向进行分区,并给出了分区标准,同时发现,加糙可使边界层紊流趋于均匀和各向同性,即概率密度分布更接近于高斯分布.     

波流边界层水沙运动数值模拟——II.泥沙运动模拟

为解析波流边界层内泥沙运动,建立了基于水动力-泥沙-床面互馈过程的波流边界层1DV泥沙数学模型,可用于模拟不同床面形态下粉沙-沙的含沙量过程。床面形态模块提供床面形态类型和相应参数;给出了平底和沙波床面粗糙高度和泥沙扩散系数的确定方法;采用了适宜粉沙及沙的制约沉速、底部参考浓度和起动剪切应力等公式;引入含沙量层化效应和制约沉降反映水动力与泥沙之间的相互影响。水槽试验资料验证表明,建立的模型较好地模拟了不同床面不同波流组合条件下的含沙量剖面。在此基础上,讨论了不同床面含沙量剖面模拟方法的差异,指出床面形态是决定含沙量变化的重要因素之一,仅通过改变床面粗糙高度不足以反映漩涡沙波床面的含沙量剖面特征。该模型可为研究波流边界层内泥沙运动和物质输运提供工具。     

刚性植被区域床面泥沙起动特性

为研究植被环境中床面泥沙起动特性,以天然河道中的植被及沙床为原型,对刚性植被环境中床面形态变化及泥沙起动规律进行了试验研究,并采用单位床面空间平均流速表述植被环境中泥沙的起动流速。研究结果表明,淹没植被条件下动床床面泥沙运动可分为3个阶段,即静止、部分起动及全部起动阶段,定义第3个阶段为刚性植被环境中泥沙起动的标准;推导的淹没刚性植被环境中的泥沙起动流速半经验函数关系式可以表述植被环境中泥沙起动特性,其中植被因子F_v对泥沙起动的影响最明显;植被密度一定时,植株周围冲坑深与空间平均流速成线性关系。研究认为采用植被环境中单位床面空间平均流速表述泥沙起动特性的方法是可行的。     

波浪渗流力对泥沙起动的影响

海岸泥沙起动的力学机理十分复杂。首先理论分析了波浪作用下床面泥沙颗粒所受的主要作用力,然后将上举力和有效重力分别与渗流力进行比较分析。计算结果表明渗流力随着波高的变化其方向也发生向上或向下的变化,渗流力与上举力和有效重力之间存在着相位差,不同粒径的泥沙在一个周期内均有一时段渗流力远大于上举力;对于渗透系数小的细颗粒泥沙,渗流力远大于泥沙颗粒的上举力和有效重力。进一步通过水槽试验分析了不同波浪参数作用下孔隙水压力梯度对泥沙起动的影响,试验结果表明波浪渗流力对海床细颗粒泥沙的起动有着重要作用,特别是在分析长周期波浪作用下细颗粒泥沙起动规律时,需考虑波浪引起的渗流力影响。     

紊流对粘性细颗粒泥沙絮凝影响

以分形聚集生长模型为基础,通过三维模拟泥沙颗粒在布朗运动、重力沉降和紊流作用下的碰撞、粘结和破碎过程,探讨了紊流对粘性细颗粒泥沙絮凝速度、絮团平均粒径及分形维数的时空影响。试验结果表明:紊流作用下,泥沙絮凝及絮团平均粒径的变化可分为加速段、等速段和减速段,絮凝等速段的作用时间及粒径等速增长段的速度均随紊动强度呈先增后减的规律,且中部区域的絮凝速度和絮团平均粒径较小;泥沙絮团分形维数随水流紊动强度的增加逐渐增大,最终趋于稳定,随时间呈先减后增的规律。模拟空间下部区域分形维数最大,上部次之,中部最小。     

粘性土及粘性土夹沙的起动规律研究

从水流脉动的机理出发,以粘性土的宾汉应力作参数,建立了粘性土的临界起动剪切力公式.根据一系列实验结果,分析了不同情况下粘性土夹沙的起动条件,归纳出统一的表达式.此外,还对上部水流为浑水的情况下,粘性土及粘性土夹沙的起动条件作了定性分析.     

Sensitivity of incipient particle motion to fluid flow penetration depth within a packed bed

A discrete element method is applied to a three‐dimensional analysis related to sediment entrainment on a micro‐scale. Sediment entrainment is the process by which a fluid medium accelerates particles from rest and advects them upward until they are either transported as bedload or suspended by the flow. Modelling of the entrainment process is a critically important aspect for studies of erosion, pollutant resuspension and transport, and formation of bedforms in environmental flows. Previous discrete element method studies of sediment entrainment have assumed the flow within the particle bed to be negligible and have only allowed for the motion of the topmost particles. At the same time, micro‐scale experimental studies indicate that there is a small slip of the fluid flow at the top of the bed, indicating the presence of non‐vanishing fluid velocity within the topmost bed layers. The current study demonstrates that the onset of particle incipient motion, which immediately precedes particle entrainment, is highly sensitive to this small fluid flow within the topmost bed layers. Using an exponential decay profile for the inner‐bed fluid flow, the discrete element method calculations are repeated with different fluid penetration depths within the bed for several small particle Reynolds numbers. For cases with slip velocity corresponding to that observed in previous experiments with natural sediment, the predicted particle velocity is found to be a few percent of the fluid velocity at the top of the viscous wall layer, which is a reasonable range of velocities for observation of incipient particle motion. This method for prescribing the fluid flow within the particle bed allows for the current discrete element method to be extended in future studies to the analysis of sediment entrainment under the influence of events such as turbulent bursting. Additionally, predictions for the slip velocities and fluid flow profile within the bed suggest the need for further experimental studies to provide the data necessary for additional improvement of the discrete element method models.     

Emplacement of massive deposits by sheet flow

Turbidity current and coastal storm deposits are commonly characterized by a basal sandy massive (structureless) unit overlying an erosional surface and underlying a parallel or cross-laminated unit. Similar sequences have been recently identified in fluvial settings as well. Notwithstanding field, laboratory and numerical studies, the mechanisms for emplacement of these massive basal units are still under debate. It is well accepted that the sequence considered here can be deposited by waning-energy flows, and that the parallel-laminated units are deposited under transport conditions corresponding to upper plane bed at the dune–antidune transition. Thus, transport conditions that are more intense than those at the dune–antidune transition should deposit massive units. This study presents experimental, open-channel flow results showing that sandy massive units can be the result of gradual deposition from a thick bedload layer of colliding grains called sheet flow layer. When this layer forms with relatively coarse sand, the non-dimensional bed shear stress associated with skin friction, the Shields number, is larger than a threshold value approximately equal to 0·4. For values of the Shields number smaller than 0·4 the sheet flow layer disappeared, sediment was transported by a standard bedload layer one or two grain diameters thick, and the bed configuration was characterized by downstream migrating antidunes and washed out dunes. Parallel laminae were found in deposits emplaced with standard bedload transport demonstrating that the same dilute flow can gradually deposit the basal and the parallel-laminated unit in presence of traction at the depositional boundary. Further, the experiments suggested that two different types of upper plane bed conditions can be defined, one associated with standard bedload transport at the dune–antidune transition, and the other associated with bedload transport in sheet flow mode at the transition between upstream and downstream migrating antidunes.     

推移质输沙率的非线性研究

目的是从非线性科学角度探讨泥沙运动的规律。通过分析非均匀沙起动的影响因素得到了尖点突变模型的状态变量和控制变量,建立了能描述推移质运动的尖点突变模型。在尖点突变标准方程的基础上,应用尖点突变理论的坐标变换和拓扑变换,导出了输沙强度和水流参数的函数关系式。将输沙强度作为状态变量,水流参数和床沙密实系数作为控制变量。用水槽实验资料和其它推移质输沙率公式进行了对比验证。验证结果表明,计算值与其它推移质输沙率公式和水槽实验结果基本相符,误差一般在-90%～80%之间。说明建立的推移质输沙率公式是合理的,能够反映泥沙的起动和输移规律。     

The motion of sediment-water mixtures during intense bedload transport: computer simulations

A new model, which couples fluid and particle dynamics, has been developed to study the motion of the sediment-water mixture during intense bedload transport, including the velocity profiles of both sediment and water, the roughness length of an upper plane bed and the thickness of moving sediment layers. Standard mixing length theory is used to model the motion of water above the boundary between the overlying water and the sediment-water mixture. The turbulent flow within the moving sediment layers is described by a shear stress model, in which the effective viscosity of the flowing water is proportional to the velocity difference between the fluid and the sediment. The particle dynamics method, in which the equations of motion of each of many particles are solved directly, is applied to model the movement of sediment particles. The particle-fluid interaction is expressed by a velocity-squared fluid drag force exerted on each sediment particle. Both computer simulation results and theoretical analysis have shown that the velocities of both sediment and fluid during intense sediment transport decrease exponentially with depth in the top layers of a fast-moving sediment—water mixture. The thickness of the moving sediment layers, obtained from the computer simulation results, is proportional to the shear stress, which agrees with previous experimental observations.     

Applications of grain-pivoting and sliding analyses to selective entrapment of gravel and to flow-competence evaluations

Models of sediment threshold by grain pivoting or sliding over underlying particles are examined in order to explore their application to evaluations of selective entrainment of gravel by flowing water. Of special interest is whether such process-based models provide satisfactory evaluations of flow competence and the movement of large clasts by floods. A detailed derivation is undertaken, focusing first on the fluid flow and forces at the particle level. The resulting threshold equation for the particle-level velocity is then modified to yield the mean entrainment stress for the flow as a whole. This approach is appropriate for considerations of selective entrainment of grains of varying sizes within a deposit, the sorting being due to their relative projection distances above the bed and the dependence of their pivoting angles on grain size and shape. The resulting threshold equations contain a number of coefficients (e.g. drag and lift) whose values are poorly known, but can be constrained by requiring agreement with the Shields curve for the threshold of grains in uniform deposits. If pivoting coefficients based on laboratory measurements with tetrahedral arrangements of particles are used in the models, smaller degrees of selective sorting are predicted than found in the field measurements of gravel entrainment. However, if reasonable modifications of those coefficients are made for expected field conditions, then the models yield good agreement with the data. Sliding models, where sorting is due entirely to projection distances of the grains above the bed, yield somewhat poorer agreement with the field data; however, the sliding models may have support from laboratory experiments on gravel entrainment in that the data and theoretical curves have similar concave trends. The existing measurements lack documentation of the mechanisms of grain movement, so it is not possible to conclusively determine the relative importance of grain pivoting versus sliding. In spite of such uncertainties, the results are encouraging and it is concluded that pivoting and sliding models for grain entrainment do have potential for field computations of selective entrainment and flow competence.     

Threshold of motion of bivalve and gastropod shells under oscillatory flow in flume experiments

The threshold of motion of non-fragmented mollusc shells was studied for the first time under oscillatory flow. In this regard, flume experiments were used to investigate the threshold of motion of three bivalve and three gastropod species, two typical mollusc classes of coastal coquina deposits. The sieve diameters ranged from 2·0 to 15·9 mm. These experiments were performed on a flat-bottom setup under regular non-breaking waves (swell) produced by a flap-type wave generator. The critical Shields values for each species of mollusc were plotted against the sieve and nominal diameter. Moreover, the dimensionless Corey shape factor of the shells was evaluated in order to investigate the effect of mollusc shell shapes on the threshold of motion. According to their critical Shields parameter, the mollusc threshold data under oscillatory flow present smaller values than the siliciclastic sediments when considering their sieve diameter. In addition, the mollusc datasets are below the empirical curves built from siliciclastic grain data under current and waves. When considering the nominal diameter, the critical Shields parameter increases and the mollusc data are closer to siliciclastic sediments. Bivalves, which have a flat-concave shape (form factor: 0·27 to 0·37), have a higher critical Shields parameter for smaller particles and more uniform datasets than the gastropod scattered data, which have a rounded shape (form factor: 0·58 to 0·62) and have varied morphologies (ellipsoidal, conical and cubic). The comparison between previous current-driven threshold data of bioclastic sediment motion and the data of mollusc whole shells under oscillatory flow shows a fair correlation on the Shields diagram, in which all datasets are below the mean empirical curves for siliciclastic sediments. These findings indicate that the shape effect on the transport initiation is predominant for smaller shells. The use of the nominal diameter is satisfactory to improve the bioclastic and siliciclastic data correlation.     

Threshold studies of gravel size particles under the co-linear combined action of waves and currents

Gravel size sediment beds are tested under the combined influence of simulated wave action and co-linear currents in a laboratory flume. Critical current speed, at threshold, increases with increasing size. Superimposed wave energy causes a small reduction in the unidirectional current energy. For low values of wave-induced near bed current velocities, the resistance to erosion increases when the wave period decreases from 10 to 6 s. Finally, combined critical shear stresses are found to be lower than those predicted using the Shields curve, as modified for oscillatory flow. Grain protrusion is suggested as a mechanism to explain this divergence.