基于图像处理的明渠紊流近壁区条带结构试验

基于图像处理技术,对明渠紊流近壁区的水流低速条带结构进行了测量和分析。条带结构采用天然沙进行示踪,使用高清相机拍摄高分辨率图像,通过形态学方法提取其几何特征。分析条带的宽度及间距的变化特点,得到结果:① 低速条带宽度从水槽中线到两侧壁逐渐变小,沿中线基本呈对称分布;② 摩阻雷诺数在276~465范围时,条带宽度和间距量纲一值的变化范围分别为20.8~29.8和 99.3~126.7,条带面积占图像总面积百分比为23.4%~15.9%;③ 随摩阻雷诺数增大,低速条带的宽度和间距加大,而条带面积百分比逐渐减小。试验结果表明,低速条带结构随水流摩阻雷诺数的变化非常复杂,而图像处理技术提供了一种高效率、高精度的测量手段。     

明渠紊流中涡结构的运动规律

明渠紊流中随水流迁移的涡结构在传质过程中具有重要作用。为观测涡结构在整个生存期内的运动过程,将两套高频粒子图像测速系统的测量范围沿水流方向拼接,实现流场的大范围测量。将上述方法用于测量摩阻雷诺数为490~870的明渠紊流的纵垂面流场,并结合涡结构跟踪算法,得到了涡的生存时间、运动速度、位移及强度变化等参数。分析表明,涡结构的水平运动与水流同步,但整体有向水面抬升的趋势,在低雷诺数条件下,最大生存时间约10倍涡周转时间,最大垂向位移约0.4倍水深;随着雷诺数的增大,涡结构的量纲一运动速度未发生明显变化,但极限强度和强度变化速率均增大,且极限强度的增加幅度相对较小,使得涡结构的生存时间缩短,纵、垂向位移相应减小。研究结果揭示了生存时间是影响涡结构运动规律的关键变量,初步得到了涡结构运动规律随水流条件的变化趋势。     

粗糙透水床面明渠水流的垂线流速分布

为了进一步揭示粗糙透水床面明渠水流运动特性,针对垂线流速分布研究存在的问题,根据边界层理论,推导了含有摩阻流速、理论床面流速和原点位移等3个参数的修正对数公式;通过水槽试验,运用激光多普勒测速仪,分别对用直径1 cm玻璃珠构成的粗糙不透水和透水床面明渠水流的垂线流速分布进行了测量。结果表明:推导的修正对数公式与实测符合很好;相同水流条件下,透水床面的摩阻流速要大于不透水床面,各自的阻力系数保持基本不变;理论床面流速是主流平均流速的0.35~0.45倍,且床面相对流速随着雷诺数的增大而略有减小。     

渗流边界上推移质输沙率

天然河道边界上常由于地下水和地表水的水位不同而产生渗流,引起水流的非静水压力分布,进而影响推移质输移规律,也使得常用的推移质输沙率公式不再适用.利用水槽试验测定沙波的尺寸及运动速度与床面渗流水力梯度的关系,结果表明沙波的高度及运动速度随着向下渗流强度的加大而增加,随着向上渗流强度的加大而减少;进而结合理论分析,利用沙波运动特点估算推移质输沙率,并推导出推移质输沙率随床面渗流的引入而变化的计算公式,分析渗流引起非静水压力分布下推移质输沙率的变化.     

泥沙知识痤讲 (三)推移质运动

河流泥沙运动的主要形式有沿床面滚动、跳跃和悬浮。此外在个别条件下可能发生滑动(当床面光滑时)和床面成层移动(当水流强度很大时)。悬浮运动的泥沙称为悬移质;滚动、滑动、跳跃、层移的泥沙称为推移质。推移质运动是由底部水动力和颗粒碰撞所引起,只发生在床面附近的区域(底层)内。1.泥沙起动当水流速度由小增大到一定数值后,床面上原来静止的个别泥沙,或因其粒径较细,或由于其位置突出,再碰上较大的纵向瞬时流速,就开始滚动。此时这种泥沙处于起动状态;而相应的水流时均速度则称为该类泥沙的起动流速。泥沙起动流速关系到开始测验推移质的时机,不冲渠道的设计、河床冲刷的极限深度,以及截流和抢险抛石的尺寸等等,故有重要的实用意义。泥沙起动时所受的作用力有:水流正面推力 P=     

非均匀沙运动特性试验

天然河流床沙通常为非均匀沙,准确把握非均匀沙颗粒运动规律是模拟和预测天然河流河床演变的基础。开展了恒定均匀流条件下的非均匀沙推移质运动水槽试验,床沙粒径范围为0.10~20 mm。利用摄像机从顶部拍摄了粗化条件下的推移质颗粒运动,获取大量非均匀沙颗粒的运动轨迹,提取了颗粒运动速度、走停时间等基本运动参数,推移质运动颗粒粒径范围为0.74~8.19 mm。试验结果表明,非均匀沙床面聚集体或大颗粒使推移质颗粒运动方向发生改变,与均匀沙成果相比,非均匀沙推移质颗粒的纵向运动速度减小,横向运动速度增大;推移质颗粒纵向运动速度遵循指数分布,单次运动速度遵循Γ分布,横向运动速度及运动速度矢量角则遵循正态分布。     

基于水下摄影的床面泥沙运动特性试验研究

为提高推移质试验的观测精度与效率,基于水下摄影和粒子跟踪（UP/PTV）技术,实现明槽流床面泥沙运动状态的精细试验观测。通过剔除床面颗粒震颤干扰、设置颗粒临界运动阈值和多重滤波筛选程序,提高采集样本数据的有效性。利用多组低强度推移质试验,提取床面颗粒运动轨迹、速度、单步时长等数据,并进行粒子运动的Lagrange过程分析和概率密度分布（PDF）研究。研究表明:粒间碰撞和近底紊流扫荡的影响使粒子速度在单步步长内呈现先急剧增加再缓慢衰减的变化特征;速度PDF曲线显示细尾Gamma函数特性,同时受粒子震颤效应影响,加速度PDF曲线则具有拉普拉斯分布特征;粒子速度与摩阻流速之间关系密切,保持3.4~3.5的比值;粒子单步时长与步长的联合分布呈幂函数变化规律,拟合曲线指数一般为1.25~1.3。幂律指数大小受推移质输沙强度与床面粒子异质性的影响。     

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

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

双峰型非均匀沙推移运动特性及输移规律

采用图像识别与推移质动态监测技术,开展基于双峰型非均匀推移质的系列水槽试验.通过引入反映床面粗糙度、粘性底层特性与颗粒非均匀度η(粗细比)的综合水流强度函数Ψ_b、特征弗劳德数Fr_b,系统研究了不同水流强度与床沙组成条件下的推移质输移特性以及颗粒非均匀度对输沙率的影响.通过对关键因子的辨识与量纲分析,提出了双峰型非均匀推移质输移模式,建立了基于近壁特征因子的水流强度Ψ_b与非均匀推移质输移强度Φ'的函数关系.对双峰型底沙输移机理的分析表明,非均匀沙的组成特征使得η成为影响Φ'的重要参量;正是细粒对粗粒的解怙作用对粗沙运动产生重要影响,使推移质输移率与颗粒非均匀度间呈现驼峰关系,峰值对应的粗细比η_c约为3∶7.     

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

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

On the relationship between flow and suspended sediment transport over the crest of a sand dune, Río Paraná, Argentina

The links between large‐scale turbulence and the suspension of sediment over alluvial bedforms have generated considerable interest in the last few decades, with past studies illustrating the origin of such turbulence and its influence on flow resistance, sediment transport and bedform morphology. In this study of turbulence and sediment suspension over large sand dunes in the Río Paraná, Argentina, time series of three‐dimensional velocity, and at‐a‐point suspended sediment concentration and particle‐size, were measured with an acoustic Doppler current profiler and laser in situ scattering transmissometer, respectively. These time series were decomposed using wavelet analysis to investigate the scales of covariation of flow velocity and suspended sediment. The analysis reveals an inverse relationship between streamwise and vertical velocities over the dune crest, where streamwise flow deceleration is linked to the vertical flux of fluid towards the water surface in the form of large turbulent fluid ejections. Regions of high suspended sediment concentration are found to correlate well with such events. The frequencies of these turbulent events have been assessed from wavelet analysis and found to concentrate in two zones that closely match predictions from empirical equations. Such a finding suggests that a combination and interaction of vortex shedding and wake flapping/changing length of the lee‐side separation zone are the principal contributors to the turbulent flow field associated with such large alluvial sand dunes. Wavelet analysis provides insight upon the temporal and spatial evolution of these coherent flow structures, including information on the topology of dune‐related turbulent flow structures. At the flow stage investigated, the turbulent flow events, and their associated high suspended sediment concentrations, are seen to grow with height above the bed until a threshold height (ca 0·45 flow depth) is reached, above which they begin to decay and dissipate.     

Quantifying velocity and turbulence structure in depositing sustained turbidity currents across breaks in slope

Two-dimensional experiments investigating sediment transport and turbulence structure in sustained turbidity currents that cross breaks in slope are presented as analogue illustrations for natural flows. The results suggest that in natural flows, turbulence generation at slope breaks may account for increased sand transport into basins and that the formation of a hydraulic jump may not be necessary to explain features such as the occurrence of submarine plunge pools and the deposition of coarser-grained beds in the bottomsets of Gilbert-type fan deltas. Experimental flows were generated on 0°, 3°, 6° and 9° slopes of equal length which terminated abruptly on a horizontal bed. Two-component velocities were measured on the slope, at the slope break and downstream of the slope break. Flows were depositional and non-uniform, visibly slowing and thickening with distance downstream. One-dimensional continuous wavelet transforms of velocity time series were used to produce time-period variance maps. Peaks in variance were tested against a background red-noise spectrum at the 95% level; a significant period banding occurs in the cross-wavelet transform at the slope break, attributed to increased formation of coherent flow structures (Kelvin–Helmholtz billows). Variance becomes distributed at progressively longer periods and the shape of the bed-normal-velocity spectral energy distribution changes with distance downstream. This is attributed to a shift towards larger turbulent structures caused by wake stretching. Mean velocity, Reynolds shear stress and turbulent kinetic energy profiles illustrate the mean distribution of turbulence through the currents. A turbulent kinetic energy transfer balance shows that flow non-uniformity arises through the transfer of mean streamwise slowing to mean bed-normal motion through the action of Reynolds normal stresses. Net turbulence production through the action of normal stresses is achieved on steeper slopes as turbulence dissipation due to mean bed-normal motion is limited. At the slope break, an imbalance between the production and dissipation of turbulence occurs because of the contrasting nature of the wall and free-shear boundaries at the bottom and top of the flows, respectively. A rapid reduction in mean streamwise velocity predominately affects the base of the flows and steeper proximal slope flows have to slow more at the break in slope. The increased turbulent kinetic energy, limited bed-normal motion and strong mixing imposed by steep proximal slopes means rapid slowing enhances turbulence production at the break in slope by focusing energy into coherent flow structures at a characteristic period. Thus, mean streamwise slowing is transferred into turbulence production at the slope break that causes increased transport of sediment and a decrease in deposit mass downstream of the slope break. The internal effects of flow non-uniformity therefore can be separated from the external influence of the slope break.     

含植物河道紊流结构研究进展

自然河流中,水生植物的存在改变了河道中的紊流结构.与普通明渠紊流相比,含植物河道水体掺混更剧烈,紊流结构更复杂.由于紊流结构直接影响河流中泥沙、污染物和营养物的输移,含植物河道紊流结构研究越来越受到重视.从紊流平均结构、紊动特性、拟序结构和紊流类型几方面归纳并总结了国内外研究进展,综合分析了影响紊流结构的因素,提出利用场与点相结合的综合分析方法,对拟序结构的形成及发展过程进行定量刻画,探究植物条件下紊流机理是今后的重要研究方向.     

浅水非线性波作用下沙纹床面底层流动特性试验研究

通过水槽试验研究浅水非线性波作用下沙纹床面底层流动特性,利用CCD图像技术观测分析非对称沙纹的形成和演化规律。利用声学多普勒测速仪(ADV)测量非对称沙纹底床上的流场,得到了不同波高、周期、水深条件下的沙纹峰顶和谷底断面的瞬时速度。试验结果分析表明,浅水非线性波作用下床面上形成非对称沙纹,其近底流速具有较强紊动特性,随着距床面距离的增大紊动强度逐渐减弱。在水流方向改变时,沙纹背部具有明显漩涡运动。沙纹背后形成的漩涡能起到维持沙纹的作用。浅水非线性波作用下,沙纹的形成原因主要是床面泥沙颗粒在非对称流动和床面近壁粘性底层中漩涡结构动力作用下,作受迫摆动、推移所致。     

泥沙起动概率的高速摄影测量方法

为了通过实验测量泥沙起动概率,基于封闭槽道内开展的泥沙起动实验,系统论述了基于高速摄影的起动概率测量方法所涉及的起动概率定义、相机采样频率、图像处理方法、起动比例与运动比例的关系等关键问题.根据泥沙在水流紊动猝发作用下起动的物理机制,定义实测起动概率为单个紊动猝发周期内,从一定面积床面上冲刷外移的泥沙比例.根据推移质运动概念模型,提出高速图片拍摄间隔应小于推移质运动的中间时间尺度,以捕获泥沙的每一次起动.在高速图像处理时,通过连续两帧图片之间的灰度差实现运动比例的无偏差提取.进一步分析表明,运动泥沙可处于起动、止动或滑翔状态,运动比例等于起动比例、止动比例和滑翔比例之和,起动比例与运动比例之比等于两张图片之间的时间间隔除以推移质运动的中间时间尺度.基于上述关键技术,将高速摄影测量方法用于泥沙起动概率的测量,结果表明,临界起动水流条件时,床沙处于少量起动现象,对应的临界起动概率与窦国仁所推导的2.28%接近.     

推移质翻越低坝输移特性的试验研究

对于山区河流低坝而言，平时淤积在坝前的推移质粗沙可能会在洪水期集中翻越坝顶，形成高强度输沙。本文开展水槽试验，研究推移质粗沙自上游起动、推进、再翻越坝顶后向下游输移的过程，分析了输沙参数的变化特性及数理规律，描述了翻坝输沙模式及运动特征，揭示了输沙规律与河床形态之间的自然联系。取得如下认识：①输沙量随时间大致以幂函数规律增长。②低坝附近区域河床形态终将趋于稳定，上游和下游均形成相对稳定的曲面斜坡淤积体。③在不同的水流强度下推移质翻坝输移模式存在差异。对于一般水流强度工况，上游淤积体曲面斜坡表面泥沙颗粒以滚动或滑动模式起动，推移至接近坝顶位置时再跃移翻坝，后向下游输移；对于更高水流强度工况，后期的翻坝输沙模式可能发生显著转变，周期性边壁漩涡成为翻坝输沙的主要动力来源。     

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

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

The importance of suppression of turbulence by stratification on the estuarine turbidity maximum

A simple numerical model demonstrates that the reduction in turbulence due to stratification greatly enhances the trapping of suspended sediment that occurs at the estuarine turbidity maximum. In moderately and highly stratified estuaries the turbulent diffusivity decreases markedly between the region upstream of the salinity intrusion, where the turbulence is uninhibited by salt stratification, and the stratified regime within the salinity intrusion, where turbulence is reduced by the inhibitory influence of salt stratification. This reduction in turbulent diffusion results in a reduction in the quantity of sediment that can be carried by the flow, causing sediment to be trapped near the landward limit of the salinity intrusion. This trapping process occurs at the same location as that due to the estuarine convergence, but it appears to be many times more effective at trapping silt-size particles. A model is formulated that is similar to Festa and Hansen's (1978) model of the estuarine turbidity maximum, with the addition of a stratification-dependent eddy diffusivity. For silt-size sediment particles, the model indicates as much as a 20-fold increase in the trapping rate with inclusion of the stratification effect. it is likely that this mechanism is important in many partially mixed and highly stratified estuaries.     

The vertical turbulence structure of experimental turbidity currents encountering basal obstructions: implications for vertical suspended sediment distribution in non‐equilibrium currents

Large roughness features, caused by erosion of the sea floor, are commonly observed on the modern sea floor and beneath turbidite sandstone beds in outcrop. This paper aims to investigate the effect of such roughness elements on the turbulent velocity field and its consequences for the sediment carrying capacity of the flows. Experimental turbidity currents were run through a rectangular channel, with a single roughness element fixed to the bottom in some runs. The effect of this roughness element on the turbulent velocity field was determined by measuring vertical profiles of the vertical velocity component in the region downstream of the basal obstruction with the Ultrasonic Doppler Velocity Profiling technique. The experiments were set up to answer two research questions. (i) How does a single roughness element alter the distribution of vertical turbulence intensity? (ii) How does the altered profile evolve in the downstream direction? The results for runs over a plane substrate are similar to data presented previously and show a lower turbulence maximum near the channel floor, a turbulence minimum associated with the velocity maximum, and a turbulence maximum associated with the upper flow interface. In the runs in which the flows were perturbed by the single roughness element, the intensity of the lower turbulence maximum was increased between 41% to 81%. This excess turbulence dissipated upwards in the flow while it travelled further downstream, but was still observable at the most distal measurement location (at a distance ca 39 times the roughness height downstream of the element). All results point towards a similarity between the near bed turbulence structure of turbidity currents and free surface shear flows that has been proposed by previous authors, and this proposition is supported further by the apparent success of a shear velocity estimation method that is based on this similarity. Theory of turbulent dispersal of suspended sediment is used to discuss how the observed turbulent effects of a single large roughness element may impact on the suspended sediment distribution in real world turbidity currents. It is concluded that this impact may consist of a non‐equilibrium net‐upwards transport of suspended sediment, counteracting density stratification. Thus, erosive substrate topography created by frontal parts of natural turbidity flows may super‐elevate sediment concentrations in upper regions above equilibrium values in following flow stages, delay depletion of the flow via sedimentation and increase their run‐out distance.     

Flow, sediment transport and bedform dynamics over the transition from dunes to upper-stage plane beds: implications for the formation of planar laminae

Preliminary results are reported from an experimental study of the interaction between turbulence, sediment transport and bedform dynamics over the transition from dunes to upper stage plane beds. Over the transition, typical dunes changed to humpback dunes (mean velocity 0–8 ms^-1, depth 01 m, mean grain size 0.3 mm) to nominally plane beds with low relief bed waves up to a few mm high. All bedforms had a mean length of 0.7–0.8 m. Hot film anemometry and flow visualization clearly show that horizontal and vertical turbulent motions in dune troughs decrease progressively through the transition while horizontal turbulence intensities increase near the bed on dune backs through to a plane bed. Average bedload and suspended load concentrations increase progressively over the transition, and the near-bed transport rate immediately downstream of flow reattachment increases markedly relative to that near dune crests. This relative increase in sediment transport near reattachment appears to be due to suppression of upward directed turbulence by increased sediment concentration, such that velocity close to the bed can increase more quickly downstream of reattachment. Low-relief bedwaves on upper-stage plane beds are ubiquitous and give rise to laterally extensive, mm-thick planar laminae; however, within such laminae are laminae of more limited lateral extent and thickness, related to the turbulent bursting process over the downstream depositional surface of the bedwaves.