桥墩-冰塞-局部冲刷相关问题研究进展

桥梁建设改变了河流的边界条件、水流条件和河床泥沙的运动状态,冬季河流中有冰塞出现时则产生相互影响作用.近年来,基于冰塞稳定性力学分析,开展了桥墩影响下冰塞稳定性研究,所得桥墩影响下冰塞稳定性判别公式计算结果和实测资料能够较好得吻合;研究了桥墩对河道卡封以及临界流凌密度的影响,计算得到的临界流凌密度与实验值较为接近;介绍了明流、冰盖下桥墩周围局部冲刷等方面的相关研究工作,得出了冰盖条件下桥墩局部最大冲刷深度的回归经验公式,结果表明该公式具有很好的适用范围.对冰塞条件下桥墩附近局部冲刷问题探索性地开展了试验研究,研究表明,冰塞条件下的局部冲刷与冰盖条件下有着很大不同,不仅冲刷深度要大的多,而且桥墩周围局部冲刷过程与冰塞厚度变化相关.在冰塞条件下,桥墩处的最大冲刷深度和最大冲刷长度均随流速增大而增大,冲刷坑尾部淤积沙丘高度也受到流速的影响.就相关研究难点和今后需要进一步开展的研究工作进行了分析和展望.      

一种新的往复流试验平台及其泥沙冲刷试验研究

本文提出一种可以模拟潮流运动的自动化控制水槽。该水槽可以时控制往复流两端的水位,实现流的往复运动,并可以按计算机要求实现水位的稳定或渐变。本文同时利用该水槽进行了圆柱开桥墩周围的测试数据、照片为基础、分析了不同水流条件下圆柱形桥墩周围的床面冲刷情况。试验发现在桥墩附近,水流的旋滚运动及旋涡系是造成床面冲刷的主要原因,另外对于往复流情况,桥墩附近床面受到的冲刷比固定水位冲刷时要严重,也更为复杂 。     

床面上直立圆柱的三维湍流数值模拟

从数值预报桥墩等结构物床面局部冲刷的角度发展绕直立圆柱的三维湍流的数值模拟技术. 基于Wilcox的k-ω两方程湍流模式，采用基于有限体积法的压力修正SIMPLE算法, 计算了绕床面直立圆柱的三维湍流流场，分析了光滑和粗糙床面两种情况下的流动情况. 通过系列的验证计算，表明该计算模型能够比较准确地反映不同外来流条件下绕直立圆柱的流场. 计算结果揭示了床面粗糙度对绕圆柱的湍流流动的影响.     

非牛顿体不稳定流的研究

本文通过试验和分析研究了非牛顿体明渠不稳定流的机理,阐明了不稳定流的产生和发展是由于流体具有屈服应力,利用这一机理,解释了高含沙洪水的波动和“浆河”现象、冰川流和泥流的不稳定流和泥石流的阵流现象。     

城市河网洪水过程的一维数值模拟

基于邻接矩阵和邻接表概念,提出了新的拓扑结构,即利用net、edge和node对象实现河网拓扑信息及河道断面和结点信息的储存;基于虚设河段法,推导出城市河网内边界即集中入流(雨水管出口)、闸堰的追赶系数,从而实现了基于GIS平台的城市河网非恒定流一维数值模拟程序的通用性和可移植性.最后,将温州市区“9.4”暴雨期间的河网洪水过程的计算值与实测值进行了对比,结果表明本模型可较好地模拟城市河网的洪水过程.     

热局部非平衡下一维半无限长饱和多孔柱体的温度场分析

对于端部受温度载荷的一维半无限长多孔介质柱体,给出了热局部非平衡下固相和流相温度场在Laplace变换域中的解析表达式.对于冲击温度载荷的情况,获得了温度场在短时间内的Laplace逆变换渐近解析解.数值分析了流、固两相热扩散系数之比以及热交换系数对固相和流相温度场的影响,比较了热局部非平衡下加权温度与热局部平衡下温度之间的差别.     

近床面水平圆柱局部冲刷二维数值模拟

为研究近床面水平圆柱的局部冲刷问题,基于N-S方程和有限体积法,在FLUENT中通过二次开发建立了局部冲刷二维数值模型。模型采用标准κ-ε紊流模型来计算水平圆柱周围的流场,同时借助FLUENT软件中的自定义函数功能提取床面剪应力参数来计算该时刻的推移质输沙率及床面节点位移变化值,然后运用动网格技术来模拟床面地形的变化,通过物理模型试验来进行模型验证。结果表明,计算结果与试验结果基本一致,从而证实了冲刷模型的准确性。     

Lattice Boltzmann方法应用实例--一维溃坝波模拟

通过分析Ｌａｔｔｉｃｅ Ｂｏｌｔｚｍａｎｎ（ＬＢ）方法的多尺度方程和一维浅水波方程,建立了计算一有渠非恒定流的ＬＢ模型,对一维溃坝波动这一典型应用实例进行了模拟分析,结果与Ｐｒｅｉｓｓｍａｎｎ隐格式及解析解十分一致,说明将ＬＢ方法引入实际的一维明渠非恒定流计算是可行的。     

非恒定孔口出流的几个新公式

非恒定孔口出流的几个新公式高进（湘潭矿业学院勘建系，湘潭４１１２０１）由水动力学知道，具有等横截面的柱形容器的非恒定孔口出流，放空容器所需的时间等于定水头孔口出流排出等量流体所需时间的两倍￣［１］．事实上，非恒定孔口出流的规律并不止于此。本文的研究表...     

管壳式换热器内气液两相流冲刷垂直管束的阻力特性研究

本文对垂直布置的管壳式换热器中，气液两相流冲刷管束时的流动阻力进行了实验研究。实验段采用管壳式换热器模型，其圆形外壳内的管束由４９根管子组成，并沿长度方向用三块折流板将管束分成四个冲刷流程。当气液两相流自下而上冲刷管束时，分别测量摩擦阻力和局部阻力，并用分相流动模型进行分析，得到良好的计算关联式。     

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.     

Numerical modelling of ice jam resistance to main channel flow

In northern countries, subfreezing temperatures during the winter season result in the formation of ice covers on most rivers. Towards the end of the winter season, during the spring break-up period, stationary ice covers become weak in strength and break up. The resulting broken ice pieces or ice floes are significantly larger in thickness and have a rougher undersurface relative to sheet ice and impose higher hydraulic resistance. The downstream movement of the ice floes may be arrested under conditions such as an intact ice cover, bridge piers or channel constrictions, among others, thereby initiating a break-up ice jam. These ice jams most often have been observed to cause very high water stages. Detrimental effects caused by these high water levels encompass those of operational and design-related problems such as the flooding of communities due to ice-jam-induced backwater, flood risk assessments, altering of the open water flow regime, bed scour and flooding of bridges. The ability to predict the influence of an ice jam on the main flow is of considerable importance in river engineering and can be viewed upon by its effects on the variation in the water surface levels. All other information is dependent on the foregoing. The ice jam influence on the main flow can be regarded with respect to local and global standpoints. The primary objective of this study is to formulate the influence of the ice jam on the main channel flow. The formulation is then coupled with a two-dimensional numerical model for the simulation of the water flow regime. The data from different laboratory experiments on ice jams are reproduced numerically. Various simulations are then carried out to compute the water surface levels and velocities in channels under ice jam conditions. The numerical results are then compared with the laboratory data. Results show that the mathematical formulation developed to predict the water surface levels and velocities along the ice jam length as well as upstream and downstream of its leading and trailing edges respectively gives satisfactory predictions.     

Zero mass error using unsteady wetting–drying conditions in shallow flows over dry irregular topography

A wetting–drying condition (WDC) for unsteady shallow water flow in two dimensions leading to zero numerical error in mass conservation is presented in this work. Some applications are shown which demonstrate the effectiveness of the WDC in flood propagation and dam break flows over real geometries. The WDC has been incorporated into a cell centred finite volume method based on Roe's approximate Riemann solver across the edges of both structured and unstructured meshes. Previous wetting–drying condition based on steady‐state conditions lead to numerical errors in unsteady cases over configurations with strong variations on bed slope. A modification of the wetting–drying condition including the normal velocity to the cell edge enables to achieve zero numerical errors. The complete numerical technique is described in this work including source terms discretization as a complete and efficient 2D river flow simulation tool. Comparisons of experimental and numerical results are shown for some of the applications. Copyright © 2004 John Wiley & Sons, Ltd.     

Theoretical analyses on bed topography responses in large depth-to-width ratio river bends with constant curvatures

Bed morphology is the result of a dynamic response to a complex meandering river system. It is an important factor for the further development of river. Based on the meandering river characterized by a large depth-to-width ratio, a theoretical model is established with the coupling of Navier-Stokes (N-S), sediment transport, and bed deformation equations. The flow characteristics and bed response of river are obtained with the perturbation method. The research results show that, under the effect of two-dimensional flow disturbance, the bars and pools present the regular response. For a given sinuousness, the amplitude of the bed response can be used as a criterion to judge the bedform stability. The effects of the Reynolds number, disturbance wavenumber, sinuousness, and bed morphology gradient on the bed response development are described.     

Bathymetry reconstruction based on the zero-inertia shallow water approximation

The knowledge of the channel bed topography is paramount in modeling the hydrodynamics of open channel flows. Indeed, flow models based on the Shallow Water Approximation require prior information on the channel bed topography to accurately capture the flow features in natural rivers, estuaries, and flood plains. We present here a numerical technique for reconstructing the channel bed topography from given free surface elevation data for steep open channel flows for which the zero-inertia shallow water approximation holds. In this context, the shallow water equations are modified by neglecting inertia terms while retaining the effects of the bed slope and friction terms. We show in this work that by algebraic manipulation, we can recast the governing equations into a single first-order partial differential equation which describes the inverse problem which consists in finding the bed topography from known free surface elevation data. Interestingly, the analysis shows that the inverse problem does not require the knowledge of the bed roughness. The forward problem is solved using MacCormack’s explicit numerical scheme by considering unsteady modified shallow water equations. However, the inverse problem is solved using the method of characteristics. The results of the inverse and the forward problem are successfully tested against each other on two different test cases.     

河型发育过程的二维数值模拟研究

改进已有二维平面水沙模型,进一步考虑河床粗化、弯道二次流及河床结构对泥沙输运的影响;并针对河湾形态对崩岸的影响提出新的非粘性土崩岸模拟方法,使其能模拟边滩的形成。利用改进后模型分别对Yen弯道水槽泥沙冲淤试验及Friendkin室内弯曲小河塑造试验进行模拟;模拟结果与试验观测资料吻合较好,表明考虑弯道对泥沙输运及河岸崩塌影响的模型能够更准确反映天然河流演化过程,为研究河型演化机理及各控制因素变化规律提供新手段。     

Flow Field Around Single and Tandem Piers

The present study provides a comparison between the flow pattern around two circular piers in tandem and a single pier set up on a moderately rough flat bed in a laboratory flume. Velocities are measured by an Acoustic Doppler Velocimeter (ADV). The contours of the time-averaged velocity components, Reynolds shear stresses, turbulence intensities and turbulence kinetic energy at different planes are presented. Streamlines and vectors are used to study the flow features. The analysis of power spectra around the piers is also presented. The results show that the presence of downstream pier changes the flow structure to a great extent, particularly in the near-wake region. Within the gap between the two piers, a stronger and substantial upflow is shaped. However, a weaker transverse-deflection is formed in comparison with that in the single pier. Near the bed, the velocity of flow approaching the downstream pier decreases to 0.2–0.3 times of the approach mean velocity (U ₀) due to the sheltering effect of the upstream pier. In the wake of downstream pier, the flow structure is completely different from the one in the wake of single pier. In comparison with the single pier, the values of turbulence kinetic energy and turbulence intensities show a considerable decrease around the tandem piers. In tandem piers, the high values of turbulence characteristics are found near the downstream pier. There is a recirculation zones just upstream of the sheltered pier close to the bed and another behind that pier near the free surface. The results show a decrease in the strength of vortical structure in the wake of tandem piers in comparison with single pier. It is shown that the formation of flow with different Reynolds number along the flow depth due to the effect of bed roughness, as well as pier spacing, can influence the type of flow regime of tandem case. In addition to enhancing the flow structure indulgence, the present detailed measurements can also be used for verification of numerical models.     

圆柱绕流局部冲刷机制的实验研究

采用粒子图像测速技术和粒子图像识别技术研究了冲刷坑的演化过程及局部冲刷的流体动力机制。首先研究了动床条件下圆柱周围的床面演化过程，给出垂直对称面输沙率随时间变化以及圆柱下游垂直对称面的沙颗粒浓度和泥沙分布规律；其次针对典型冲刷坑形状固定底床，测量了圆柱下游的瞬时流场，并给出平均流场。结果表明尾流区平均剪应力较低，但垂向涡量和湍流强度较高。分析指出平均剪应力较低的尾流区域，局部冲刷的主要动力机制是尾涡和湍流的综合作用。