波浪作用下沙纹床面底层流动特性研究

涡动沙纹很大程度决定着床面阻力和波浪衰减,针对涡动沙纹作用下的绕流结构具有强二维特性,建立了立面二维曲线坐标下的数值模型.利用国外公开发表的试验数据进行验证,表明模型能模拟波浪作用下底层分离涡随时间、空间变化的动力特性,对底层流动的紊动特性、床面剪切应力的随时间、空间的变化以及涡动结构的尺度分别进行了讨论.     

波流边界层水沙运动数值模拟——I:水动力模拟

波流边界层水动力模拟对研究波流相互作用和泥沙运动具有重要的理论意义和实践价值。开发了波流边界层1DV垂向一维水动力数值模型,可用于模拟漩涡沙波床面和平底床面水动力特征。模型的构建基于边界层控制方程,平底床面采用k-ε模型,沙波床面采用双层模型,提出了漩涡层和紊动扩散层交界面紊动动能和紊动耗散率表达式。试验资料验证表明,模型较好地模拟了波浪-水流-床面共同作用下的边界层水动力特征,包括波周期内不同相位流速分布、紊动动能、剪切应力等以及波致时均流速分布和波流相互作用下的时均流速分布等。根据所建模型,讨论了不同床面和波流组合条件下的水动力特征。该模型可为研究波流边界层内水动力特征提供工具。     

海啸波作用下泥沙运动——Ⅱ.床沙组成变化分析

在波浪水槽实验的基础上,对海啸波作用下的床沙组成变化规律开展研究。实验采用1/10~1/20的组合坡度,选取N波作为入射波。实验对波高进行采集,对波浪的上爬、回落和水跃过程进行拍摄记录,对每个波作用后的地形进行测量,并对初始和最终的床面顶层泥沙进行采样筛分。实验结果表明,N波作用下泥沙在离岸区水跃发生区域堆积,淤积沙坝泥沙粒径呈粗化趋势。同时采用规则波和非规则波进行对比,波浪作用后形成岸滩也为沙坝剖面,淤积沙坝泥沙粒径呈细化趋势,结果均遵循Çelikoğlu提出的泥沙运动的基本规律,细颗粒泥沙会在强烈的紊动作用下从床面中被筛选出来,并被搬运到低紊动地区,此过程造成了剧烈紊动区泥沙的粗化。     

淹没植物明渠床面冲淤及其对水流运动的影响

植物的存在改变了河流水动力特性,造成独特的床面冲淤态势。利用实验室水槽模拟含淹没植物的河道,对床面形态和紊流统计特性参数进行测量,研究不同类型紊流作用下的床面冲淤特征以及床面起伏对流动的影响。结果表明:床面剪切紊流条件下,床面形态为马蹄坑-沙沟/沙脊与沙波复合分布,床面变形加剧了流速沿水深不均匀分布并促进水流动量交换;在自由剪切混合层紊流条件下,床面形态为植物根部马蹄形冲坑及其后方沙沟、沙脊交错分布,床面变形对流动的影响并不显著;“类二重紊流”条件下,床面形态同样表现为马蹄坑-沙沟/沙脊-沙波复合,床面变形促进植物层内部的水流动量交换、抑制紊动清扫,抑制植物层外部的动量交换、促进紊动喷射。     

模拟降雨条件下坡面流水动力学特性研究

为探明降雨条件下沙黄土坡面水流水动力学特性,以流体力学和泥沙运动力学理论为依据,通过5个坡度和5个雨强组合条件下室内模拟降雨试验,系统研究了坡面水流水力参数的变化规律。结果表明,降雨条件下坡面薄层水流雷诺数均小于580,处于层流失稳区;水流流型随雨强和坡度的变化而发生转捩,当坡度较缓、雨强较小时,床面形态处于低能态区和过渡区,床面出现沙纹现象,水流宏观上呈缓流,反之,坡度较陡、雨强较大时,床面由沙纹和沙垄向动平床过渡,宏观上多呈急流;并根据薄层水流阻力组成特点,推导出沙黄土坡面薄层水流阻力计算公式。验证结果表明,该式误差较小,可为坡面侵蚀预报模型的构建提供参考。     

一种改进的近岸波浪破碎数值模型

提出一种基于完全非线性Boussinesq方程的改进波浪破碎模型,用于模拟近岸浅水波浪破碎.模型借助水滚的概念,提出一种确定k方程紊动模式中紊动源项的计算公式,通过解k方程计算Boussinesq方程中的涡黏系数,实现对破碎波的模拟.岸边界采用窄缝法,使得模型可用于波浪爬坡的计算.用实验室实测波高和增减水资料对模型进行了验证,得到了一致的结果.紊动源、紊动动能以及紊动耗散率的计算结果表明:①在破波点处紊动源项值最大,随着波浪向岸边传播,逐渐减小;②破波点处,水平方向的对流和扩散在紊动能量输移中发挥重要作用;③岸边附近紊动源与紊动耗散接近平衡.     

破碎波作用下沙质海床床面形态变化试验

波浪在斜坡沙质海床上破碎会加剧泥沙输移导致海床形态变化,研究破碎波作用下沙质海床形态变化机制对于岸滩演变分析极为重要。在波浪水槽中采用中值粒径0.47 mm原型沙铺设1∶20坡度的底床模型进行试验研究,测量不同波浪条件下床面形态和沙坝顶端悬浮泥沙浓度变化。通过测量和计算破碎带输沙率、沙坝尺度和沙坝移动速度,分析破碎波作用下沙质斜坡海床上床面形态变化规律。试验结果表明,破碎带沙坝顶端的悬浮泥沙浓度与水深和底部床面密切相关,在形成沙坝和沙坝水平方向移动时,悬浮泥沙浓度较大;斜坡上沙坝前后来回运动的周期大小具有随机性,沙坝既有向岸又有离岸移动;在多组波长时间作用后沙坝尺度趋于稳定,底床净输沙量趋于0。     

破碎波作用下沙质海床床面形态变化试验

波浪在斜坡沙质海床上破碎会加剧泥沙输移导致海床形态变化,研究破碎波作用下沙质海床形态变化机制对于岸滩演变分析极为重要。在波浪水槽中采用中值粒径0.47 mm原型沙铺设1∶20坡度的底床模型进行试验研究,测量不同波浪条件下床面形态和沙坝顶端悬浮泥沙浓度变化。通过测量和计算破碎带输沙率、沙坝尺度和沙坝移动速度,分析破碎波作用下沙质斜坡海床上床面形态变化规律。试验结果表明,破碎带沙坝顶端的悬浮泥沙浓度与水深和底部床面密切相关,在形成沙坝和沙坝水平方向移动时,悬浮泥沙浓度较大;斜坡上沙坝前后来回运动的周期大小具有随机性,沙坝既有向岸又有离岸移动;在多组波长时间作用后沙坝尺度趋于稳定,底床净输沙量趋于0。     

过渡区动床明渠流的流速分布

利用激光流速仪测量了随水流强度加大松散床面,从定床到具有一定输沙强度的动床流动,试验水流条件均处于过渡区。试验资料分析表明,动床过渡区明渠流的卡门常数小于通常所采用的0.4;水流强度逐渐增大,随着推移质运动量增加,阻力系数有减小再增大的趋势;雷诺应力分布和紊动强度分布规律同光滑明渠流类似。     

30°分水角明渠水流紊动特性试验研究

灌区分水口多为直角, 附近普遍存在冲刷、淤积及结构破坏现象。天然河流分、汇流多为锐角, 为探索非直角分水口区域明渠水流紊动特性, 以30°分水角明渠为研究对象, 采用声学多普勒点式流速仪测量各断面上的三维瞬时流速, 分析典型断面上的时均流速、环流强度、紊动强度、紊动能和床面剪切应力分布规律。试验结果显示: 分水口处水流横向流速沿渠宽方向变化较大, 易形成环流, 导致局部冲淤现象明显; 靠近口门侧水流紊动强度变化剧烈且分布无规律, 最大值出现在口门上唇断面; 分水口处下层水体紊动能明显大于上层水体, 紊动能较大值多集中在口门附近; 口门断面下唇附近床面剪切应力较大, 易出现冲刷侵蚀, 破坏渠道稳定; 与直角分水口相比, 分水角为30°时, 泥沙颗粒不易进入侧渠道, 水流对渠底及边壁侵蚀速率较低。研究结果可为灌区渠系设计及运行提供参考。     

Flow and turbulence structure across the ripple–dune transition: an experiment under mobile bed conditions

Current knowledge of flow and turbulent processes acting across the sand bed continuum is still unable to unequivocally explain the mechanism(s) by which ripples become dunes. Understanding has been improved by comparative high-resolution studies undertaken over fixed bedforms at different stages in the continuum. However, these studies both ignore the role of mobile sediment and do not examine flow structure during the actual transition from ripples to dunes. The aims of the paper are: (i) to describe flow and turbulence characteristics acting above mobile bedforms at several stages across the transition; and (ii) to compare these data with those arising from experiments over fixed ripples and dunes. Laboratory experiments are presented that examine the turbulence structure across seven distinct stages of the transition from ripples to dunes. Single-point acoustic Doppler velocimeter sampling at three flow heights above a developing mobile boundary was undertaken. Time-averaged statistics and the instantaneous quadrant record reveal distinct changes in flow structure either side of the change from ripples to dunes. Initially, shear-related, high-frequency vortex shedding dominates turbulence production. This increases until two-dimensional (2D) dunes have formed. Thereafter, turbulence intensities and Reynolds stress decline and three-dimensional dunes exhibit values found over 2D ripples. This is the result of shear layer dampening which occurs when the topographically-accelerated downstream velocity increases at a faster rate than flow depth. Activity at reattachment increases due to high velocity fluid imparting high mass and momentum transfer at the bed and/or wake flapping. Suspended sediment may also play a role in turbulence dampening and bed erosion. Ejections dominate over sweeps in terms of event frequency but not magnitude. Strong relationships between inward interactions and sweeps, and ejections and outward interactions, suggest that mass and momentum exchanges are dependent upon activity in all four quadrants. The results contradict the notion present in most physical models that larger bedforms exhibit most shear layer activity. Consequently an improved model for the ripple–dune transition is proposed.     

A field study of turbulence and sediment dynamics over subaqueous dunes with flow separation

This study examines flow, turbulence and sand suspension over large dunes in Canoe Pass, a distributary channel of the Fraser River delta, Canada. Dune morphology is characterized by a symmetrical shape and steep leeside slopes over 30°. Velocity was measured with an electromagnetic current meter and suspended sand concentration with four optical backscatter (OBS) probes. The general patterns of time-averaged velocity and sand suspension are consistent with previous studies, including an increase in mean velocity and decrease in turbulence intensity and sand concentration with height above the bed, reversed flow with high turbulence intensity and high sand concentrations in the leeside flow separation zone and an increase in near-bed velocity and sand concentration along the stoss side of the dune. Frequency spectra of near-bed velocity and OBS records from leeside separation zones are composed of two distinct frequencies, providing field confirmation of scale relations based on flume experiments. The low-frequency spectral signal probably results from wake flapping and the high-frequency signal from vortex shedding. The wake-flapping frequency predominates outside the separation zone and is linked to turbulent structures that suspend sand. Predictions from a depth-scale Strouhal Law show good agreement with measured wake-flapping frequencies. Cross-correlations of OBS records reveal that turbulent sand suspension structures advect downstream at 23–25° above the horizontal. These advection angles are similar to coherent flow structures measured in flumes and to sand suspension structures visualized over large dunes in the field.     

Depositional processes,bedform development and hybrid bed formation in rapidly decelerated cohesive (mud–sand) sediment flows

Flows with high suspended sediment concentrations are common in many sedimentary environments, and their flow properties may show a transitional behaviour between fully turbulent and quasi‐laminar plug flows. The characteristics of these transitional flows are known to be a function of both clay concentration and type, as well as the applied fluid stress, but so far the interaction of these transitional flows with a loose sediment bed has received little attention. Information on this type of interaction is essential for the recognition and prediction of sedimentary structures formed by cohesive transitional flows in, for example, fluvial, estuarine and deep‐marine deposits. This paper investigates the behaviour of rapidly decelerated to steady flows that contain a mixture of sand, silt and clay, and explores the effect of different clay (kaolin) concentrations on the dynamics of flow over a mobile bed, and the bedforms and stratification produced. Experiments were conducted in a recirculating slurry flume capable of transporting high clay concentrations. Ultrasonic Doppler velocity profiling was used to measure the flow velocity within these concentrated suspension flows. The development of current ripples under decelerated flows of differing kaolin concentration was documented and evolution of their height, wavelength and migration rate quantified. This work confirms past work over smooth, fixed beds which showed that, as clay concentration rises, a distinct sequence of flow types is generated: turbulent flow, turbulence‐enhanced transitional flow, lower transitional plug flow, upper transitional plug flow and a quasi‐laminar plug flow. Each of these flow types produces an initial flat bed upon rapid flow deceleration, followed by reworking of these deposits through the development of current ripples during the subsequent steady flow in turbulent flow, turbulence‐enhanced transitional flow and lower transitional plug flow. The initial flat beds are structureless, but have diagnostic textural properties, caused by differential settling of sand, silt and cohesive mud, which forms characteristic bipartite beds that initially consist of sand overlain by silt or clay. As clay concentration in the formative flow increases, ripples first increase in mean height and wavelength under turbulence‐enhanced transitional flow and lower transitional plug‐flow regimes, which is attributed to the additional turbulence generated under these flows that subsequently causes greater lee side erosion. As clay concentration increases further from a lower transitional plug flow, ripples cease to exist under the upper transitional plug flow and quasi‐laminar plug flow conditions investigated herein. This disappearance of ripples appears due to both turbulence suppression at higher clay concentrations, as well as the increasing shear strength of the bed sediment that becomes more difficult to erode as clay concentration increases. The stratification within the ripples formed after rapid deceleration of the transitional flows reflects the availability of sediment from the bipartite bed. The exact nature of the ripple cross‐stratification in these flows is a direct function of the duration of the formative flow and the texture of the initial flat bed, and ripples do not form in cohesive flows with a Reynolds number smaller than ca 12 000. Examples are given of how the unique properties of the current ripples and plane beds, developing below decelerated transitional flows, could aid in the interpretation of depositional processes in modern and ancient sediments. This interpretation includes a new model for hybrid beds that explains their formation in terms of a combination of vertical grain‐size segregation and longitudinal flow transformation.     

Visualization of wave-induced suspension patterns over two-dimensional bedforms

High spatial and temporal resolution measurements of suspended sand concentration ( c ) over vortex ripples were collected with a three-transducer acoustic backscatter sensor (ABS) array, under irregular `natural' waves in a multidirectional wave basin. These measurements permit two-dimensional visualization of the movement of sediment-laden vortices over an individual vortex ripple under a series of waves. Patterns of sediment motion were tracked through consecutive zero-crossings in the horizontal velocity ( U ) record measured at 0·05 m above the ripple crest elevation. It was possible to trace the advection of individual sediment-laden vortices at the zero-crossings. During 73% of these events, shedding and advection of coherent suspension events occurred before the flow reversal associated with the zero-crossing. This may be caused by the bedforms retarding the near-bed flow inducing the eddy shedding before the zero-crossing. While at maxima in U , secondary suspension events with low c were observed to pass over the ripple crest moving with U measured at 0·05 m. This pattern is attributed to vortex shedding from adjacent bedforms and/or antecedent suspension events. The most energetic events appeared to persist for several wave cycles and reached heights of ≈0·20 m. These suspension events appeared to be more persistent when smaller waves follow larger waves, possibly as a result of weaker reversals in vorticity. Although the events appeared to be vertically coherent in the time series from the individual transducers, it is apparent through visualization that these events are associated with the pairing of antecedent and developing vortices.     

Sediment suspension under waves and currents: time scales and vertical structure

Field measurements of the vertical structure of near-bed suspended sediment concentrations were obtained from arrays of fast response optical backscatter suspended solids sensors to examine the time-dependent response of sediment resuspension to waves and currents and the constraints imposed by bedforms. Data were recorded from both a nonbarred, marine shoreface and a barred lacustrine shoreface, under both shoaling and breaking waves (significant heights of 0·25–1·50m; peak periods of 3 and 8 s) and in water depths of 0·5–5·0 m. Sediment concentrations are positively correlated with increasing elevation above the bed, but lagged in time. The time lag varies directly with separation distance between measurement locations and inversely with the horizontal component of the near-bed oscillatory velocity. Both the presence of wave groups and the settling velocities of the sediment particules in suspension influence the temporal changes in concentration at a given elevation. Sediment concentrations appear to respond more slowly to the incident wind-wave forcing with distance away from the bed as a result of two factors: (1) the sequential increase in concentration induced by a succession of large waves in a group; and (ii) the relative increase in finer sediments with smaller settling velocities. Bedforms interact with the near-bed horizontal currents to impose a distinct constraint upon the timing of suspension events relative to the phase of the fluid motion, and, therefore, the vertical structure of the suspended sediment concentration at a range of time scales. The near-bed concentrations appear to be strongly dependent upon the vertical convection of sediment associated with the ejection from the wave boundary layer of separation vortices generated in the lee of ripple crests. Concentration gradients in the presence of vortex ripples are large, as are the correlation between concentrations measured at different elevations within the fluid.     

Effect of velocity hiatuses in oscillatory flow on migration and geometry of ripples: wave‐flume experiments

A series of wave‐flume experiments was conducted to closely look at characteristics of geometry and migration of wave‐generated ripples, with particular reference to the effect of velocity ‘hiatuses’ during which the near‐bed flow velocity becomes much smaller than the threshold of sediment movement. Three types of wave patterns were generated: two types for simulating waves with intervening velocity hiatuses; and regular waves for comparison purposes. In the former two types, two different wavelengths of water waves were generated alternately in the course of a wave test: the wave with a longer wavelength was set large enough to mobilize the bottom sediment, whereas the wave with a shorter wavelength was set too small to mobilize the sediment. The former two types were designed to be different in sequence of convexity and concavity of wave patterns. The sequence with the convex–concave longer wave and successive convex–concave shorter wave was described as a ‘zero‐up‐crossing’ wave pattern, and the inverse sequence was described as a ‘zero‐down‐crossing’ wave pattern. The ripples developed under oscillatory flow with intervening hiatuses manifested the following characteristics in geometry and migration. (i) The morphological characteristics of ripples, namely wavelength, height and the ripple steepness, are unaffected by the intervening hiatuses of velocity. (ii) The directions of ripple migration under the zero‐up‐crossing and zero‐down‐crossing wave patterns corresponded well with the directions of the flow immediately before onset of the hiatuses. (iii) The observation of sand particle movement on the ripple surface indicated that, under the zero‐up‐crossing waves, the velocity hiatus prevents the entrained sediment cloud from being thrown onshore, and thus the sediment grains thrown onshore are fewer than those thrown offshore. As a result of the sediment movement over one wave‐cycle, the net sediment transport is directed offshore under the zero‐up‐crossing wave pattern. (iv) The velocity of ripple migration was highly correlated with acceleration skewness. Under most of the zero‐up‐crossing (zero‐down‐crossing) wave patterns, flow acceleration skewed negative (positive) and ripples migrated offshore (onshore).     

Response of sand ripples to change in oscillatory flow

Ripples take time to evolve to a new equilibrium state in response to a change in wave-generated oscillatory flow. The paper presents results from flow tunnel experiments designed to examine oscillatory flow transient ripple processes under controlled, full-scale laboratory conditions. The experiments include study of the growth of ripples from flat bed and the evolution of existing ripples to new equilibrium ripples in response to a step change in the flow. In general, ripples evolve through a combination of two main processes: (i) from a flat bed or from a bed consisting of ripples that are smaller than the equilibrium ripples through a combination of 'slide' and 'merge'; (ii) from a bed consisting of ripples that are larger than the equilibrium ripples through a combination of 'split' and 'merge'. The experimental results show that equilibrium ripple geometry is independent of initial bed morphology while the time to reach equilibrium is largely independent of the initial bed and the equilibrium ripple size. The time to reach equilibrium depends strongly on the mobility number, and a new empirical equation relating mobility number and the number of flow cycles to equilibrium is proposed. This equation is combined with a simple exponential function for ripple height growth or decay to produce a new empirical model for ripple height evolution, which gives a reasonably good overall agreement with the measurements. The model is based on experiments involving one sediment size only and further work is needed to develop the model for other sand sizes.