曲线坐标系二维带自由表面强紊动水流数值模拟

将追踪自由表面的流体体积(VOF)法应用于曲线坐标系下水流控制方程的求解中,计入流线弯曲对水流紊动特性的影响,建立了垂向二维强紊动水流的曲率修正的紊流模型,并对溢流坝反弧段的紊流特性进行了数值模拟。数值计算时,采用有限体积法离散水流的控制方程;物理变量,如:压力P、紊动参量κ、ε、γt等,采用交错方式排列(交错网格布置),用SIMPLEC算法求解离散方程。计算结果表明,得到的溢流坝反弧段的自由表面位置、速度场、压力场、剪应力分布和紊动能分布与实验结果吻合良好。     

温度和盐度分层流的数值模拟

以N-S方程和k-ε湍流模型为基础,针对水温和盐度分层流的流动特性,考虑了浮力对紊动的影响,建立了温度和盐度分层流的k-ε模型,对同时存在温度和盐度梯度的一类密度分层流进行了数值模拟计算,其中紊动平均量控制方程中的紊动项是由k-ε紊流模型计算,给出了速度、温度和盐度的分布。计算结果正确地反映了分层及紊动特征,与实测结果吻合较好。     

顺直河流横向紊动扩散系数

利用垂直紊动扩散系数及水利几何形态关系,借助于抛物线型断面形态方程,提出了顺直河道中局部水深沿横断面的分布,在此基础上确定了横向紊动扩散系数的断面分布及其平均值表达式,阐明了无量纲横向紊流扩散系数间的关系,计算的断面平均横向紊动扩散系数与138组试验资料吻合良好。比较结果表明,建立的顺直河道横向紊流扩散系数计算公式能给出与实测值最接近的预测值。与现有的其它横向紊流扩散系数计算公式相比,其公式在理论上更加合理,机理上更加清楚,并且具有最小的预测误差。     

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

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

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

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

采样频率及样本容量对明渠紊流统计值的影响

为了研究采样频率和样本容量对明渠紊流统计特征值的影响,利用粒子图像测速技术,对明渠恒定均匀水流的二维流速分布进行测量,获得了包含340535个样本的流场时间序列.对该大容量时间序列进行抽样,得到不同采样频率和容量条件下的一系列子序列.对比分析各子序列的时均流速、紊动强度及雷诺应力的变化特点,发现存在临界的采样频率和样本容量,紊动统计特征参数在临界值以下波动较大,而在临界值以上基本趋于稳定;要满足同一测量精度,时均流速对采样频率和样本容量的要求最低,紊动强度次之,雷诺应力最高.     

水动力学条件对复氧过程的影响研究

通过水槽实验对紊动水体复氧规律进行了研究.将氧通量理论应用于紊动水体复氧过程的研究,得到了紊动水体表面传质系数.数值模拟研究了实验流场分布和紊动特性,得出了紊动水体表面传质系数与流速和紊动动能的关系式.不同工况的实验资料验证表明了紊动复氧理论的合理性和可靠性.     

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

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

颗粒沉降的格子Boltzmann模拟与PIV实验验证

在格子Boltzmann方法中引入大涡模拟,对球形颗粒在静水中沉降引起的紊动流场进行了数值模拟。数值模拟沉速与理论值以及粒子图像测速系统(PIV)实验结果吻合,验证了模型的合理性。同时分析比较了颗粒沉降过程中尾部紊动流场分布以及尾流流速值,发现数值模拟结果与实测结果趋势、数值基本一致,进一步说明了利用格子Boltzmann方法与大涡模拟技术相结合可以合理模拟泥沙颗粒在紊流区的沉降。     

淹没刚性植被明渠紊流沿程流动特性差异

淹没刚性植被通过改变水流结构,造成时均流速、紊动强度、雷诺应力等水力参数垂线分布不均匀和沿程差异化。采用声学多普勒测速仪(ADV)测量3种淹没度(3.0、4.0、5.0)下的流速,结合统计学方法,系统分析植被段内及其上、下游过渡段流速和紊动特性差异。结果表明：植被显著增强了水流紊动,且紊动特性取值与淹没度正相关;植被段内流速差异在低淹没度下的植被层内和高淹没度下的自由流动层内更加显著,而紊动特性沿程增强,且垂线分布具有相似性,最大值点位于冠层顶部附近;当淹没度满足KH涡的形成和发展条件时,随淹没度的增大,植被段内紊动特性垂线分布出现转折点(临近此点梯度急剧减小并趋于0)的断面数量增多,经验证,在充分发展的紊流区此点可作为KH涡的上边界点。     

单向流边界层泥沙起动规律

Shields曲线常用于表示泥沙起动的临界条件,基于边界层理论,对Shields曲线各个流区的线型进行了推证;考虑粘结力的作用,对Shields参数及Shields曲线进行了修正,并给出修正Shields曲线表达式;在此基础之上,从边界层角度重新阐述了Shields曲线。结果表明:Shields曲线在光滑紊流及层流区呈直线分布,在过渡区与阻力系数线型保持一致,在粗糙紊流区呈水平直线分布;修正后Shields曲线与原始Shields曲线在形式上保持一致,修正Shields曲线表达式与实测数据吻合较好,适用于粗、细颗粒泥沙起动条件的计算;Shields曲线事实上代表了Shields参数与沙粒周围绕流流态的关系,同一颗粒处于不同流区起动时,其起动切应力不同。     

Aerodynamic entrainment threshold: effects of boundary layer flow conditions

Experimental data are presented demonstrating the influence of boundary layer flow conditions on aerodynamic entrainment of grains in the absence of intersaltation collisions. New methods are proposed for (1) the unambiguous determination of aerodynamic threshold for any grain population and (2) approximation of the probability density function (PDF) distributions of threshold shear velocity for aerodynamic entrainment. In wind tunnel experiments, the orderly spatial development of flow conditions within a developing boundary layer over the roughened surface of a flat plate constrains the aerodynamic threshold condition in terms of both mean and fluctuating values. Initial grain dislodgements and subsequent erosion from narrow strips of loose, finely fractionated ballotini were recorded photographically as wind speed was increased. Boundary layer parameters, including average threshold shear velocity (U_*t), were calculated using the momentum integral method. Direct observations show that sporadic oscillation of grains preceded dislodgement. At slightly higher velocities most grains rolled over their neighbours before entering saltation. Initial entrainment in spatially semi-organized flurries of 50 or more grains was followed by quiescent periods at airflow velocities close to threshold. These observations provide strong circumstantial evidence linking both the nature and spatial pattern of initial grain motions to sweep events during the fluid bursting process. For each grain fraction, values of U_*t were found to span an unexpectedly wide range and to decrease downwind from the leading edge of the plate as turbulence intensity increased. A probabilistic entrainment model is applied to the aerodynamic threshold condition so as to incorporate the effects of changing turbulent flow regimes over the plate. Analysis of strip erosion curves gives both an objective definition of the threshold condition and usable approximations of the PDF for U_*t required by the model and for future stochastic treatment of the threshold condition.     

波流边界层水沙运动数值模拟——Ⅲ.底部高含沙层模拟

为了解析底部高含沙层特征,研究了物理影响机制和时均含沙量剖面。采用波流边界层1DV泥沙数学模型,结合试验资料,通过设置不同的计算工况进行了影响因素敏感性分析;在此基础上,提出波浪相关的时均泥沙扩散系数分布,考虑主要影响机制,推导了基于边界层物理过程的波浪作用下时均含沙量剖面表达式。结果表明,底部高含沙层与波浪边界层密切相关,是受水动力和床面形态综合影响的结果。仅建立高含沙层与水动力或床面形态的单一关系是有局限性的。含沙量层化效应和制约沉速对底部高含沙层具有重要影响。提出的平底床面含沙量剖面表达式为幂函数-Rouse-指数分布,漩涡沙波床面为指数-幂函数-Rouse分布。预期可应用于二维和三维泥沙数值模拟。     

In situ acoustic measurements of marine gravel threshold and transport

Measurements of the nearbed turbulent current flow and the bedload transport of marine gravel have been made over three tidal cycles. The turbulence in the bottom boundary layer was measured using two electromagnetic current meters, and the gravel transport was measured using a passive acoustic system which monitored the interparticle collision noise of locally mobile material. Visual estimates of bedload were also obtained with an underwater TV camera. The acoustic technique, unlike a conventional bedload sampler, has allowed estimates of transport to be obtained with a temporal resolution comparable with the turbulence data collected. This has enabled a detailed comparison to be made between the turbulent flow and the sediment response to the instantaneous flow conditions. The results of the study show that of the turbulent bursting events which contribute towards the Reynolds stress, only the sweeps and outward interactions play a significant role in the transport of coarse sedimentary material. The measurements show that it is the instantaneous increases in the horizontal turbulent velocity fluctuations that generate excess shear stresses which drive the transport process.     

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

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

Understanding the turbulent structure of the atmospheric boundary layer: A diagnostic approach

In this paper, we have attempted a diagnostic study of the turbulence characteristics of the ABL by means of two one-dimensional models. The first model uses a first order non-local closure, based on the Transilient Turbulence Theory, for parameterizing turbulent fluxes. while the second model uses second order local closure for parameterizing these. The models have been applied to conduct case studies using the Kytoon data taken at Kharagpur, during 17th–21st June, 1990, as part of the MONTBLEX programme. Our findings bring out various interesting features regarding the non-local and local turbulent statistics such as kinematic fluxes, turbulence kinetic energy, vertical velocity variance, the contribution of the eddies of various sizes to the fluxes at different level and the mixing lengths. The one-dimensional anisotropy of the turbulent eddies has been revealed by the findings from the transilient model. The vertical variation of the turbulence kinetic energy, as computed directly by the second order model, is found to be strongly correlated with the vertical velocity variance. In particular, for stably stratified boundary layers, identification of two distinct zones of the turbulence kinetic energy and corresponding vertical velocity maxima is possible, which has been interpreted as positive evidence of patchy turbulence in the boundary layer.     

极浅水边界层的沉积环境效应

极浅水环境是水深远小于正常边界层厚度的环境,潮滩滩面和潮水沟就经常处于这种环境之中。来自江苏潮滩的观察资料表明,落潮后期滩面和潮水沟在薄层低速水流作用下可形成浅水波痕、平床和次级潮水沟形态,而边界层内的流速结构仍然得以维持,物质输运强度也与Von Kármán-Prandtl模型所定义的u100数值相一致。涨急时段中潮位附近形成的滩面涌潮是极浅水边界层的另一种动力学行为,它代表薄层高速水流作用下发生的底部边界层系统的崩溃,此时Von Kármán-Prandtl模型所刻画的流速结构不复存在。滩面涌潮高度是系统崩溃的临界水深,它可以定量地表示为Hb=4z0。极浅水边界层过程对潮滩沉积和地貌形成具有独特的作用。     

Mean flow and turbulence structure over fixed, two-dimensional dunes: implications for sediment transport and bedform stability

Detailed measurements of flow velocity and its turbulent fluctuation were obtained over fixed, two-dimensional dunes in a laboratory channel. Laser Doppler anemometry was used to measure the downstream and vertical components of velocity at more than 1800 points over one dune wavelength. The density of the sampling grid allowed construction of a unique set of contour maps for all mean flow and turbulence parameters, which are assessed using higher moment measures and quadrant analysis. These flow field maps illustrate that: (1) the time-averaged downstream and vertical velocities agree well with previous studies of quasi-equilibrium flow over fixed and mobile bedforms and show a remarkable symmetry from crest to crest; (2) the maximum root-mean-square (RMS) of the downstream velocity values occur at and just downstream of flow reattachment and within the flow separation cell; (3) the maximum vertical RMS values occur within and above the zone of flow separation along the shear layer and this zone advects and diffuses downstream, extending almost to the next crest; (4) positive downstream skewness values occur within the separation cell, whereas positive vertical skewness values are restricted to the shear layer; (5) the highest Reynolds stresses are located within the zone of flow separation and along the shear layer; (6) high-magnitude, high-frequency quadrant-2 events (‘ejections’) are concentrated along the shear layer (Kelvin-Helmholtz instabilities) and dominate the contribution to the local Reynolds stress; and (7) high-magnitude, high-frequency quadrant-4 events occur bounding the separation zone, near reattachment and close to the dune crest, and are significant contributors to the local Reynolds stress at each location. These data demonstrate that the turbulence structure associated with dunes is controlled intrinsically by the formation, magnitude and downstream extent of the flow separation zone and resultant shear layer. Furthermore, the origin of dune-related macroturbulence lies in the dynamics of the shear layer rather than classical turbulent boundary layer bursting. The fluid dynamic distinction between dunes and ripples is reasoned to be linked to the velocity differential across the shear layer and hence the magnitude of the Kelvin-Helmholtz instabilities, which are both greater for dunes than ripples. These instabilities control the local flow and turbulence structure and dictate the modes of sediment entrainment and their transport rates.     

Upper-regime parallel lamination as the result of turbulent sediment transport and low-amplitude bed forms

A series of laboratory experiments has been carried out in which parallel-laminated deposits were produced from an upper-regime plane bed. The laminae had thicknesses of a few mm and could be traced continuously over distances up to the length and width of the depositional area (0–3 m by 1–5m). Fluctuations in bed elevation were measured both during deposition and at equilibrium; much of the bed fluctuation occurs at time scales that are too long to be due directly to turbulence, as most theories for lamina formation would require. We suggest instead that extremely low-amplitude bed forms are present even on apparently flat beds and that the migration of these bed forms produces laterally continuous lamination. All the lamination produced in the laboratory experiments was normally graded. Using high-speed photography it was observed that the normal grading results from rapid deposition of a layer of loosely packed coarse sand several grain-diameters thick followed by the slow sieving-out of a well packed surface layer of finer sand. The initial deposition is the result of small-scale turbulent fluctuations in boundary shear stress. The sieving-out that follows results in a smooth surface whose low friction coefficient temporarily inhibits further deposition; we term this process ‘glazing’. The alignment of small-scale turbulent scour-and-fill structures along the paths traced by migrating bedform troughs produces laterally continuous parallel lamination.     

地球表层运动和变形的GPS描述

利用IERS所公布的分在全球各大构造板块上的 6 5 7个GPS、SLR和VLBI连续观测站在ITRF框架下的速度场资料 ,采用刚体板块运动 +板块整体均匀应变 +板块内局部不均匀应变的变形分析模型 ,研究了全球各主要板块的运动和变形。结果表明板块的整体变形在统计上均不显著。在一级近似上板块间表现出来的整体相对运动显著 ,根据这些运动参数定量研究了板块边界的相对运动的大小和性质。认为地球的双重不对称变形可能主要表现为南北、东西两半球所含的板块边界的运动方式不同所致。板块内的局部不均匀变形明显 ,为板块内部可能应划分成次一级的活动地块提供了佐证。由于观测点分布的密度和均匀性不足 ,本文未能就板内不均匀变形作进一步的深入讨论。