细颗粒泥沙净冲刷和输移的大涡模拟研究

在传统水沙输移数值模拟研究中一般采用雷诺时均模拟技术(Reynolds-averaged simulation,RANS).与RANS相比,大涡模拟技术(large eddy simulation,LES)能够更加精确反映细部流动结构,计算机的发展使得采用LES探讨水流和泥沙运动规律成为可能.本文尝试给出净冲刷条件下悬沙计算的边界条件,采用动态亚格子模式对循环槽道和长槽道中的水流运动和泥沙输移进行了三维大涡模拟研究.利用直接数值模拟(direct numerical simulation,DNS)结果对LES模型进行了率定,计算结果符合良好,在此基础上初步探讨了泥沙浓度、湍动强度和湍动通量等的分布特征.结果表明,净冲刷条件下输沙平衡时泥沙浓度符合Rouse公式分布,单向流动中泥沙浓度沿着流向逐渐增大.泥沙浓度湍动强度和湍动通量都在近底部达到最大值,沿着垂向迅速减小.湍动黏性系数和扩散系数基本上在水深中间处达到最大.湍动Schmidt数沿着水深方向不是常数,在近底部和自由水面附近较大,在水深中间处较小.     

河口混合与泥沙输运

根据振荡边界层理论和波流分解方法,导出了河口往复水流的流速垂向分布廓线,据此建立了河口垂向准二维水流、盐度、泥沙运动模型。对泥沙输运,完整地考虑了其对流、扩散、起动和沉降的动力学过程。模拟结果与实测资料进行了对比。应用该模型研究河口泥沙输运,分析了河口混合对泥沙输运的影响及最大浑浊带的时空变化规律。     

THE INTERACTIONS BETWEEN WAVE AND COHESIVE SEDIMENT

淤泥质海岸泥沙运动活跃,泥沙淤积问题严重,直接关系到岸滩演变和海岸防护. 该文主要综合国内外有关研究成果,对淤泥质海岸黏性泥沙的流变特性,波浪与底泥之间的相互作用,波浪衰减以及波浪作用下的底泥质量输移作简要总结,供进一步研究参考.     

风沙运动力学机理研究的历史、进展与趋势

回顾了风沙运动力学机理研究历史,阐述了中国在风沙运动力学机理研究方 面所取得的最新进展，提出了目前存在的问题和今后发展的趋势.     

明渠中悬移质的弥散-对流方程及悬浮机理

悬移质泥沙通常构成冲积河流总输沙量的主体, 研究悬移质的悬浮机理具有重要的意义. 以双流体模型为基础, 通过引入弥散速度的概念, 建立了悬移质泥沙的输沙方程以及泥沙扩散系数的本构关系. 应用该方程分析了二维明渠均匀流中悬移质泥沙浓度垂向分布规律, 并与Einstein 和Chien 的泥沙浓度实验资料及经典扩散理论进行了对比. 以此为基础, 分析了紊动扩散、颗粒自身的紊动、颗粒碰撞应力对泥沙悬浮的影响在垂向上的变化, 以及浓度、粒径等对这些因素的影响. 结果表明, 泥沙颗粒在明渠紊流中的扩散是浑水的紊动扩散、颗粒自身的紊动、颗粒碰撞应力3 部分不同机制共同作用的结果, 把泥沙颗粒的悬浮简单归因于水流的紊动是不全面的.     

Multiple time scales of fluvial processes — theory and applications

Fluvial processes comprise water flow, sediment transport and bed evolution, which normally feature distinct time scales. The time scales of sediment transport and bed deformation relative to the flow essentially measure how fast sediment transport adapts to capacity region in line with local flow scenario and the bed deforms in comparison with the flow, which literally dictates if a capacity based and/or decoupled model is justified. This paper synthesizes the recently developed multiscale theory for sediment-laden flows over erodible bed, with bed load and suspended load transport, respectively. It is unravelled that bed load transport can adapt to capacity sufficiently rapidly even under highly unsteady flows and thus a capacity model is mostly applicable, whereas a non-capacity model is critical for suspended sediment because of the lower rate of adaptation to capacity. Physically coupled modelling is critical for fluvial processes characterized by rapid bed variation. Applications are outlined on very active bed load sediment transported by flash floods and landslide dam break floods.     

一维河网泥沙输移模式研究

推导了一种显式的河网 泥沙输移递推计算方法, 该方法在水动力学计算结果所确定的汊点类型基础上, 直接计算每 个河段中逐个节点的不平衡泥沙输沙过程以及悬沙和床沙的级配调整, 克服了传统方法中必 须通过首末断面含沙量关系确定河网悬沙矩阵求解的困难. 将该泥沙输移模式的一维河网水 沙数学模型, 应用于长江三峡下游荆江-洞庭湖河网区域, 并取得了较好的计算结果.在模型 验证之后, 进一步讨论了采用新的泥沙输移计算模式和原有泥沙输移计算模式在计算 时间上和计算内容方面的差别.     

Three-dimensional mathematical model and its application in the neighborhood of the Three Gorges Reservoir dam in the Yangtze River

The calculation of flow and sediment transport is one of the most important tasks in river engineering. The task is particularly difficult because a number of complex physical phenomena should be accounted for more realistically in a model with a predictive power. Three-dimensional calculations of river flow and suspended sediment transport are performed in this paper with application in the Three Gorges Reservoir in the Yangtze River. A period of 76 years after the dam is built is simulated and the results are compared with laboratory measurements obtained by Tsinghua University whereby the model is verified and calibrated. Generally speaking, the calculated results agree well with the experiments, demonstrating that the present model can be used for flow and sediment transport prediction in rivers. The project supported by the National Natural Science Foundation of China (50009004)     

Comparison of Kinetic and Hybrid-Equilibrium Models Simulating Colloid-Facilitated Contaminant Transport in Porous Media

The presence of colloidal particles in groundwater can enhance contaminant transport by reducing retardation effects and carrying them to distances further than predicted by a conventional advective/dispersive equation with normal retardation values. When colloids exist in porous media and affect contaminant migration, the system can best be simulated as a three-phase medium. Mechanisms of mass transfer from one phase to another by colloids and contaminants can be kinetic or equilibrium-based, depending on the sorption–desorption reaction rate between the aqueous and solid phases. When the rate of sorption between the water phase and the solid phase(s) is not much greater than the rate of change in contaminant concentration in the water phase, kinetic sorption models may better describe the phenomenon. In some cases of modeling one or more mass transfer processes, a useful simplification may be to introduce the local equilibrium assumption. In this study, the local equilibrium assumption for sorption processes on colloidal surfaces (hybrid equilibrium model) was compared with kinetic-based models. Sensitivity analyses were conducted to deduce the effect of major parameters on contaminant transport. The results obtained from the hybrid equilibrium model in predicting the transport of colloid-facilitated groundwater contaminant are very similar to those of the kinetic model, when the point of interest is not at contaminant and colloid source vicinities and the time of interest is sufficiently long for imposed sources.     

页岩气和致密砂岩气藏微裂缝气体传输特性

页岩气和致密砂岩气藏发育微裂缝,其开度多在纳米级和微米级尺度且变化大,因此微裂缝气体传输机理异常复杂.本文基于滑脱流动和努森扩散模型,分别以分子之间碰撞频率和分子与壁面碰撞频率占总碰撞频率的比值作为滑脱流动和努森扩散的权重系数,耦合这两种传输机理,建立了微裂缝气体传输模型. 该模型考虑微裂缝形状和尺度对气体传输的影响. 模型可靠性用分子模拟数据验证.结果表明:(1)模型能够合理描述微裂缝中所有气体传输机理,包括连续流动,滑脱流动和过渡流动;(2)模型能够描述不同开发阶段,微裂缝中各气体传输机理对传输贡献的逐渐变化过程;(3)微裂缝形状和尺度影响气体传输,相同开度且宽度越大的微裂缝,气体传输能力越强,且在高压和微裂缝大开度的情况下表现更明显.      

Improvements on bed‐shear stress formulation for pier scour computation

This paper introduces an improved formula for the bed‐shear stress by applying the vorticity effect and its application in a 3D flow and sediment model to estimate scouring around bridge piers. Up to now, the sediment transport formulae used for computing pier scour were developed based on the general scouring in unobstructed flow. The capability for numerical models to predict local scour around bridge piers was severely restricted by the sediment transport formulae. The new formula introduced in this paper can take into account vortices that affect the local scour process by adding some terms into the classic bed‐shear stress equation. The 3D numerical model system used in this study consists of three modules: (a) an unsteady hydrodynamic module; (b) a sediment transport module; and (c) a Fation module. The hydrodynamic module is based on the 3D RANS equations. The sediment transport module is comprised of semi empirical models of suspended load and non‐equilibrium bed load. The bed‐deformation module is based on the mass balance for sediment. The model was used to simulate pier scour in tree different test cases: (1) a circular pier; (2) a square pier; and (3) a rectangular pier, by applying the ordinary sediment equation and the newly introduced sediment equation. Results of both numerical simulations were compared against laboratory measured data and also in case 1 with result of Olsen and Melaaen (J. Hydraul. Eng. 1993; 119 (9):1048–1054). Comparisons show that the new sediment formula could predict the scour more accurately than the ordinary one. Copyright © 2010 John Wiley & Sons, Ltd.     

方柱尾流流态对床面冲淤形态的影响

应用平面二维悬沙数学模型对方柱尾流区的泥沙输运及床面调整进行了数值模拟，并对重要物理参数做了分析说明。数值计算采用时间分裂一全隐式有限差分格式，流场的计算基于沿水深平均的RANS方程。通过系统的数值模拟，揭示了不同来流情况下，方柱尾流区的流态变化，及与之相应的床面变化规律，并对不同流态下的泥沙运动作了分析。计算结果表明在方柱尾流区，不同的流场流态对尾流区的床面调整有明显影响。流场较弱的情况下，尾流区中湍动强度相应较弱，此时悬浮泥沙由于流速减小而普遍落淤，床面应力的减小也致使床面冲刷量减小。随着流场强度的逐渐增大，尾流区中湍动强度相应增强，床面应力增大，同时湍流的强烈交换作用增强了对泥沙的输运作用，床面变化不再是简单的淤积状态，部分区域出现了冲刷。     

Empirical model for estimating vertical concentration profiles of re-suspended,sediment-associated contaminants

Vertical distribution processes of sediment con-taminants in water were studied by flume experiments. Experimental results show that settling velocity of sedi-ment particles and turbulence characteristics are the major hydrodynamic factors impacting distribution of pollutants, especially near the bottom where particle diameter is similar in size to vortex structure.Sediment distribution was uniform along the distance, while contaminant distribution slightly lagged behind the sediment.The smaller the initial sediment concentration was,the more time it took to achieve a uniform concentration distribution for suspended sediment. A con-taminants transportation equation was established depending on mass conservation equations.Two mathematical estima-tion models of pollutant distribution in the overlying water considering adsorption and desorption were devised based on vertical distribution of suspended sediment:equilibrium par-tition model and dynamic micro-diffusion model.The ratio of time scale between the sediment movement and sorption can be used as the index of the models.When this ratio was large,the equilibrium assumption was reasonable,but when it was small,it might require dynamic micro-diffusion model.     

泥沙输运模拟综述——现状及其发展趋势

随着流体力学中数值方法的飞速发展,计算模型已经成为研究流体运动,泥沙输运和不同环境(诸如河流,湖泊及沿海地区)中相应污染物归趋过程等方面非常有吸引力的工具,在过去的30多年里,发展了许多计算水动力学/泥沙输运模型.文章追溯当前具有代表性的(一维、二维、三维)模型的发展历程,描述他们各自的特点,优势及局限,力图作为对模型方面感兴趣读者的第一指南,同时也为大家讨论模型的局限性,未来的发展趋势和研究需求方面搭建一个平台.给出了模型的表达,时空特性,水动力学和沉积物的耦合方式,处理非恒定流,推移质和悬移质,泥沙交换过程,泥沙类型(粘性或非粘性)及非均匀泥沙输运的能力.总结了不同模型的应用实例,读者可以运用这些例子作为研究模型设置,模型率定及模型验证的参考.给出了选择泥沙输运模型应遵循的原则,模型输入及率定方面存在的问题及改进的途径.探讨了现有水动力学/泥沙输运模型在处理复杂湍流,泥沙携带,流动与输沙耦合,非均匀泥沙,离散和扩散系数,河岸来沙处理等方面的局限性及改进的方向.最后,对基于多相流思想的泥沙输运模型及其它一些交叉性问题作了评述与展望.     

Modelling sediment transport over partially non‐erodible bottoms

In‐depth‐averaged and cross‐section‐averaged morphodynamic models, based on explicit time integration, it may happen that the computed bed level becomes lower than the top level of a non‐erodible layer (e.g. concrete, bedrock or armoured layer). This is a standard pitfall, which has been addressed in different ways. In this paper, we present an original approach for avoiding computation of non‐physical bed levels, using an iterative procedure to correct the outward sediment fluxes. The procedure is shown to be computationally efficient while it achieves a high accuracy in terms of mass conservation. We compare our original approach with the existing Struiksma's method and with a reformulation of the problem in terms of mathematical optimization of a linear or nonlinear objective function under linear constraints. The new procedure has been incorporated into an existing finite volume morphodynamic model. It has been validated with several 1D benchmarks leading to configurations with sediment transport over non‐erodible bottom. The computation time has been verified not to increase by more than 15% compared with runs without non‐erodible bottom. Copyright © 2011 John Wiley & Sons, Ltd.     

GAS-WATER TWO-PHASE DIFFUSION MODEL FOR UNSATURATED LAYERED COVER SYSTEMS

为了研究垃圾填埋场内部生物降解产生的气体(以甲烷为主) 和氧气等在非饱和覆盖层中的运移规律,建立了气体在成层非饱和覆盖层中的一维扩散模型. 该模型考虑了覆盖层含水量瞬态变化和气体扩散的耦合作用. 采用有限元多物理场耦合分析软件COMSOL Multiphysics 4.3进行求解. 其计算结果与前人实验结果能符合得很好. 研究结果表明,对典型垃圾填埋场中厚度为1m的非饱和黏土覆盖层,含水量变化对覆盖层中气体污染物的扩散运移有显著影响. 与含水量稳态变化情况相比,20 d到100 d的运移时间内同一覆盖层深度处耦合含水量作用下,气体相对浓度可以相差接近8倍. 含水量是气体在非饱和覆盖层中不稳定运移的重要影响因素.     

INTERNAL SOLITARY WAVE LOADS EXPERIMENTS AND ITS THEORETICAL MODEL FOR A CYLINDRICAL STRUCTURE

在大型重力式密度分层水槽中, 对内孤立波与圆柱型结构的相互作用特性开展了系列实验. 基于两层流体中 内孤立波的KdV,eKdV和MCC理论, 建立了圆柱型结构内孤立波载荷的理论预报模型, 给出了该载荷理论预报模型中3类内孤立波理论的适用性条件.研究表明, 圆柱型结构内孤立波水平载荷包括水平Froude-Krylov力、附加质量力和拖曳力3个部分, 可以由Morison公式计算, 而内孤立波垂向载荷主要为垂向Froude-Krylov力, 可以由内孤立波诱导动压力计算.系列实验结果表明, 附加质量系数可以取为常数1.0, 拖曳力系数与内孤立波诱导速度场的雷诺数之间为指数函数关系, 而且基于理论预报模型的数值结果与系列实验结果吻合.     

Dynamic release process of pollutants during suspended sediment transport in aquatic system

Pollutants release is highly consistent with suspended sediment concentration （SSC） in water column, especially during re-suspension and transport events. The present research focuses on pollutant dynamic release from re-suspended sediment, especially the vertical distribution relationship between them. The sediment erosion experiments on a series of uniform flow are conducted in a circulate flume. Reactive tracer （phosphorus） is used as the contaminant in fine-grained sediments to identify the release characteristic length and time. Experimental results show that the flow condition near-bed depends on the sediment surface roughness. The region with high turbulent intensities corresponds to a high concentration sediment layer. In addition, the SSC decreases with the distance, water depth, and particle grain size. The sediment in a smaller grain size takes much more time to reach equilibrium concentration. Total phosphorus （TP） concentration changes along the water depth as SSC in the initial re-suspension stage, appearing in two obvious concentration regimes： the upper low-concentration layer and the high-concentration near-bottom layer. This layered phenomenon remains for about 3 hours until SSC distri- bution tends to be uniform. Longitudinal desorption plays an important role in long-way transport to reduce the amount of suspended sediment in water column.     

考虑损伤的内变量黏弹-黏塑性本构方程

基于Rice 不可逆内变量热力学框架,在约束构型空间中讨论材料的蠕变损伤问题. 通过给定具体的余能密度函数和内变量演化方程推导出考虑损伤的内变量黏弹-黏塑性本构方程. 通过模型相似材料单轴蠕变加卸载试验对一维情况下的本构方程进行参数辨识和模型验证,本构方程能很好地描述黏弹性变形和各蠕变阶段.不同的蠕变阶段具有不同的能量耗散特点. 受应力扰动后,不考虑损伤的材料系统能自发趋于热力学平衡态或稳定态. 在考虑损伤的整个蠕变过程中,材料系统先趋于平衡态再背离平衡态发展. 能量耗散率可作为材料系统热力学状态偏离平衡态的测度;能量耗散率的时间导数可用于表征系统的演化趋势;两者的域内积分值可作为结构长期稳定性的评价指标.      

基于颗粒离散元的沙纹演变大涡模拟研究

应用大涡模拟、"点球" 浸入边界法和基于"事件驱动" 模型和"弹簧-阻尼" 模型的颗粒离散元法, 数值模拟了明渠湍流中沙纹的演变. 通过对不同谢尔兹数下无沙纹床面的推移质输沙率进行计算, 并与经典输沙率公式进行对比, 验证了模型的精度和可靠性. 随后, 对明渠湍流中沙纹床面的演变过程进行了数值模拟, 计算了有沙纹床面的推移质输沙率、沙纹长度和高度、等效床面高度的最大值、最小值和平均值、沙纹形状阻力、体积流速随时间的变化曲线. 研究发现:初始平整的床面在较短的时间内(tU_b=h≈100) 发展出数条沙纹, 随后沙纹逐渐发展, 在tU_b=h为1 600~2 000 时, 沙纹发生合并. 在沙纹数量不变的条件下, 沙纹高度随时间近似呈线性增长, 而沙纹的长度的平均值却保持恒定. 随着沙纹高度的增大, 输沙率和体积流速逐渐降低, 沙纹形状阻力则逐渐增大;当沙纹发生合并时, 沙纹高度快速增加, 输沙率、体积流速和沙纹形状阻力也出现了大幅跳跃. 在同等的水流强度条件下, 有沙纹床面的输沙率小于平整床面的输沙率.