常曲率弯道二维流动稳定性与非线性演化规律分析

国内外对弯曲河流的研究从线性特性转移到不稳定性与非线性特性上来, 弯曲河流作为一个不稳定的系统, 主要是受到水流和边界的不稳定性的影响.针对弯曲河流中存在的水流动力的不稳定性与非线性特征,利用弱非线性理论与摄动法进行展开, 在微弯条件下建立了时间模式下的常曲率 弯道二维水流扰动幅值与扰动幅角的非线性演化方程.研究了在不同弯曲度下的扰动发展特征,探讨了扰动流场在弯曲度影响下的时空分布的规律,分析了扰动流场中出现的扰动漩涡的运动过程, 阐述了弯曲度对弯曲边界内水流不稳定性的影响, 具体表现为在微弯的情况下, 弯道的弯曲度增大会提高河湾内水流的稳定性, 扰动振幅与扰动幅角会随弯曲度的增大其衰减速度更快, 并且扰动流速的对称性在弯曲度较大时减弱, 逐渐向弯道的凸岸偏斜, 在中性状态附近 的弯道水流动力具有对流不稳定性和非线性衰减的特征.本文的研究成果为构建水流动力非线性与床面形态几何非线性相互作用的模型提供了思路.      

蜿蜒边界下平面流的线性流动稳定性

为便于数值分析,蜿蜒河流水动力和演变模型中一般隐性假设二次时均流-二次涡的关系与明渠流时均流-明渠湍流的关系相同,但由于高雷诺数下的DNS算力限制和实验尺度限制,这种隐含假设是否成立目前尚无相关湍流研究来支撑.文章试图通过分析明渠湍流和二次湍流发展初期的研究,侧面揭示其湍流结构的异同.通过对曲线正交坐标系下的平面二维NS方程使用双参数摄动的方法,建立了一种求解蜿蜒边界弱非线性层流的摄动解法,并推导得出一个适用于蜿蜒边界的EOS方程以及其特征值问题的解法.蜿蜒边界下弱非线性层流解为一系列蜿蜒谐波分量的叠加,其中线性部分使得两壁产生流速差,非线性部分随着雷诺数增大呈指数增长.水流的扰动增长率特征谱的第一模态与直道流相似,由3条曲线、4个波段合成,但其长波段和短波段的扰动流场与直道流不同,所有短波段的扰动流速近似于KH涡.蜿蜒边界对内部水流扰动有一定的选择性.偏角幅值越大扰动增长越快;蜿蜒波数的影响则为先增后减,有一个使扰动增长最快的蜿蜒波数.扰动流场由一个典型的TS波和一对波包形式的二次涡叠加而成,波包只有纵向流速分量,包络线由蜿蜒波数控制,波包内是与直道扰动波参数相同的TS波.     

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

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

Theoretical analyses on hydrodynamic instability in narrow deep river with variable curvature

Bending river is the most common river type in nature, and it is also a typical example for river evolution. The transform of the flow pattern can affect the development of the riverbed form. In return, the variation in the riverbed form can affect the hydrodynamic characteristics of the flow, thereby leading to the continuous evolution of the bending river. Based on this, a theoretical model for the bending river is established. The hydrodynamic instability characteristics of the laminar flow in the channel with a variable curvature, a typical model as the meandering river, are studied, and the variations of some parameters such as the curvature, the wave number, and the wave frequency are also discussed.     

摆动河槽水动力稳定性特征分析

河流形态与水动力结构息息相关,形态约束水动力结构,水动力结构则通过泥沙运动进一步塑造形态,在自然界河流中形成一对辩证互馈关系.天然河流形态形式多样,大致可以分为顺直、微弯、分叉和散乱游荡几种类型,其中微弯及多个弯曲构成的河型为河流动力演化中最重要的一环.多个弯曲构成的河型可用正弦派生曲线来描述,它也是天然河流主槽与水动力结构复杂相互作用的结果.作为探讨这一过程的力学作用机理,构建摆动槽道并研究槽道摆动与其内部流动结构的互馈关系,既是流体力学研究的热点内容,也是目前河流动力过程研究的基础内容.在此重点讨论这一互馈关系前一部分,即水流对摆动边界的响应.文中建立了随体坐标系下摆动河槽与内部水流动力响应理论模型,通过给定摆动弯曲槽道的不同特征参数,研究讨论了正弦派生型摆动边界下的槽道水流动力稳定性特征,明确了弯曲槽道摆动对其内部主流及扰动水流结构的影响,确定弯曲槽道摆动波数、摆动频率对扰动流发展影响的相应参数定量关系,得到了槽道弯曲度和摆动特征对其内部水流不同尺度扰动影响的阈值选择性范围.     

Time and spatially resolved measurements of interfacial waves in vertical annular flow

Film thickness measurements have been performed in a vertical air/water annular flow in a pipe of 0.05 m diameter. A sensor has been built which allows to measure the film thickness evolution in time at 320 positions, such that the interface of the vertical annular flow can be reconstructed. The large-scale structures moving on the interface are described statistically, with a special attention to the disturbance waves. Probability density functions and mean statistics are given for the height, length, velocity, frequency and spatial distribution of the disturbance waves. In particular, it is shown that the disturbance waves are three-dimensional structures with large height fluctuations in the circumferential and axial direction, giving a meandering path between the maximum height around the circumference. It is also shown that the disturbance waves can flow with a slight inclination with respect to the axial direction. Finally, the disturbance waves are shown to be located randomly in space, within a Gamma distribution whose order only depends on the liquid superficial velocity. Due to the nature of the Gamma distribution, it could indicate that the spatial distribution of the disturbance waves results from a cascade of coalescence processes between the original disturbance waves on the film.     

????????????????????????о?

通过动床试验,结合黄河下游细沙河床研究了大尺度缩窄河道时非恒定流桥墩的局部冲刷. 研究结果表明,非恒定流洪水过程中的桥墩局部冲刷与恒定流的情况相比,两者之间存在差 异而且有必然的联系. 通过影响因素和试验数据的分析,建立了两者之间的关系和相应的非 恒定流桥墩局部冲刷深度计算方法. 经与原型非恒定流洪水实测资料对比,结果合理且比较 符合实际情况.     

Determination of real-time predictors of the wind turbine wake meandering

The present work proposes an experimental methodology to characterize the unsteady properties of a wind turbine wake, called meandering, and particularly its ability to follow the large-scale motions induced by large turbulent eddies contained in the approach flow. The measurements were made in an atmospheric boundary layer wind tunnel. The wind turbine model is based on the actuator disc concept. One part of the work has been dedicated to the development of a methodology for horizontal wake tracking by mean of a transverse hot wire rake, whose dynamic response is adequate for spectral analysis. Spectral coherence analysis shows that the horizontal position of the wake correlates well with the upstream transverse velocity, especially for wavelength larger than three times the diameter of the disc but less so for smaller scales. Therefore, it is concluded that the wake is actually a rather passive tracer of the large surrounding turbulent structures. The influence of the rotor size and downstream distance on the wake meandering is studied. The fluctuations of the lateral force and the yawing torque affecting the wind turbine model are also measured and correlated with the wake meandering. Two approach flow configurations are then tested: an undisturbed incoming flow (modelled atmospheric boundary layer) and a disturbed incoming flow, with a wind turbine model located upstream. Results showed that the meandering process is amplified by the presence of the upstream wake. It is shown that the coherence between the lateral force fluctuations and the horizontal wake position is significant up to length scales larger than twice the wind turbine model diameter. This leads to the conclusion that the lateral force is a better candidate than the upstream transverse velocity to predict in real time the meandering process, for either undisturbed (wake free) or disturbed incoming atmospheric flows.     

Development and testing of a simple 2D finite volume model of sub‐critical shallow water flow

Many environmental applications of shallow water flow modelling can be characterized as only slowly varying and everywhere sub‐critical. A simplified finite volume model is therefore developed that is capable of describing pertinent shallow water flow processes more efficiently than the usual Godunov/ Riemann characteristics approaches. The model is tested against a number of analytical and numerical solutions to the governing equations. The model reproduces accurately flow round a circular bend, flow over topography, flow up an initially dry beach and floodwave propagation down a meandering river reach, with mass conservative solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical analysis of turbulent structure in compound meandering open channel by algebraic Reynolds stress model

The turbulent flow in a compound meandering channel with a rectangular cross section is one of the most complicated turbulent flows, because the flow behaviour is influenced by several kinds of forces, including centrifugal forces, pressure‐driven forces and shear stresses generated by momentum transfer between the main channel and the flood plain. Numerical analysis has been performed for the fully developed turbulent flow in a compound meandering open‐channel flow using an algebraic Reynolds stress model. The boundary‐fitted coordinate system is introduced as a method for coordinate transformation in order to set the boundary conditions along the complicated shape of the meandering open channel. The turbulence model consists of transport equations for turbulent energy and dissipation, in conjunction with an algebraic stress model based on the Reynolds stress transport equations. With reference to the pressure–strain term, we have made use of a modified pressure–strain term. The boundary condition of the fluctuating vertical velocity is set to zero not only for the free surface, but also for computational grid points next to the free surface, because experimental results have shown that the fluctuating vertical velocity approaches zero near the free surface. In order to examine the validity of the present numerical method and the turbulent model, the calculated results are compared with experimental data measured by laser Doppler anemometer. In addition, the compound meandering open channel is clarified somewhat based on the calculated results. As a result of the analysis, the present algebraic Reynolds stress model is shown to be able to reasonably predict the turbulent flow in a compound meandering open channel. Copyright © 2005 John Wiley & Sons, Ltd.     

Meandering jets in shallow rectangular reservoirs: POD analysis and identification of coherent structures

The effect of the shallowness on meandering jets in a shallow rectangular reservoir is investigated. Four meandering flows were investigated in an experimental shallow rectangular reservoir. Their boundary conditions were chosen to cover a large range of friction numbers (defined with the sudden expansion width). Due to the unsteady characteristics of the flows, a proper orthogonal decomposition (POD) of the fluctuating part of the surface velocity fields measured using Large-Scale Particle Image Velocity was used for discriminating the flow structures responsible for the meandering of the jet. Less than 1 % of the calculated POD modes significantly contribute to the meandering of the jet, and two types of instability are in competition in such a flow configuration. The sinuous mode is the dominant mode in the flow, and it induces the meandering of the flow, while the varicose mode is a source of local mixing and weakly participates to the flow. The fluctuating velocity fields were then reconstructed using the POD modes corresponding to 80 % of the total mean fluctuating kinetic energy, and the coherent structures were identified using the residual vorticity, their centres being localised using a topology algorithm. The trajectories of the structures centres emphasise that at high friction number the coherent structures are small and laterally paired in the near, middle and far fields of the jet, while with decreasing friction number, the structures merge into large horizontal vortices in the far field of the jet, their trajectories showing more variability in space and time. The analysis of the stability regime finally reveals that the sinuous mode is convectively unstable and may become absolutely unstable at the end of the reservoir when the friction number is small.     

Turbulent flow structure at concordant and discordant open-channel confluences

Models of flow at river-channel confluences that consist of two concordant confluent channels with avalanche faces dipping into a scour zone are limited because this morphology may be the exception rather than the rule in nature. In this paper the mean and turbulent flow structure in the streamwise and vertical directions at both concordant and discordant laboratory confluences were examined in order to determine the effect of bed discordance on the flow field, and to assess its influence on sediment transport. Instantaneous velocities were measured with a laser Doppler anemometer using a dense spatial sampling grid. The spatial distribution of normal stress varies with bed geometry as bed discordance generates a distortion of the mixing layer between the confluent streams. Turbulent shear stress is larger in the discordant bed case and its peak is associated with the position of the mixing layer whereas for concordant beds the zone of mixing is characterised by a decrease in the Reynolds shear stress. Quadrant analysis also revealed differential dominating quadrants between the two bed geometries which will influence sediment transport routing and, consequently, the resulting bed morphology. These results highlight the need for significant modifications to current models of confluence flow dynamics in order to account for the bed configuration.We would like to thank Phil Fields of the Earth Sciences workshop at Leeds for his skillful construction of the confluence model and continued help during the course of this research. The LDA facility has been funded through grants from the UFC, NERC and University of Leeds. PB thanks NSERC for financial support and the Fonds FCAR for funding an eight month study-visit held at the Department of Earth Sciences, University of Leeds. AGR acknowledges the financial support of NSERC whilst JLB is grateful for the award of a Nuffield Science Foundation Fellowship which facilitated the preparation of this paper. We thank two anonymous reviewers for their constructive criticisms that helped clarifying the paper.     

Identification and analysis of the meandering of a fin-tip vortex using Proper Orthogonal Decomposition (POD)

The meandering of a vortex exists in a broad range of engineering applications and can lead to flow instability and other undesirable characteristics. Compared to a static vortex, measurement of a meandering vortex can result in a ‘smeared’ mean-flow field and increased levels of turbulence at the centre of the vortex. A case study was performed here on the meandering nature of a fin-tip vortex generated by a manoeuvring submarine. From stereoscopic particle image velocimetry (SPIV) measurements, it is possible to remove the meandering by shifting each instantaneous velocity field so as to produce a common centre for the vortex. In this paper, a snapshot Proper Orthogonal Decomposition (POD) technique is used to capture the dominant large-scale coherent structures (from inspection of eigenvalue or energy distributions) and to improve vortex centre identification. The POD reconstructed velocity field using only the most energetic modes enabled the coherent structures of the flow to be clearly visualised, providing improved identification of the vortex centre and subsequent evaluation of the meandering effect on the turbulent statistics. The present findings suggest that the vortex meandering only has a small impact on the ensemble-averaged resultant velocity, while contributing up to a maximum of 28% for the fluctuating component. The meandering correction also leads to an overall decrease of turbulence intensity in the peak fluctuating region of the vortex core.     

单一水平轴风电机组尾迹的模拟方法与流动机理研究综述

为实现碳达峰、碳中和“3060”目标, 风能将在我国能源体系发挥重要作用. 风力机尾迹是影响风电性能和度电成本的关键因素, 需在风力机布置和控制设计中充分考虑. 本文首先介绍风力机尾迹的数值模拟方法, 包括解析模型、低阶模型、大涡模拟和来流湍流生成方法. 解析模型和低阶模型可快速计算风力机尾迹, 但依赖于模型参数, 且不能或不能准确预测尾迹湍流特性. 结合风力机参数化模型的大涡模拟可准确预测尾迹蜿蜒等湍流特征, 是流动机理研究的有力工具, 可为发展快速预测模型提供数据和理论支撑. 接着, 本文介绍了叶尖涡、中心涡和尾迹蜿蜒并讨论其产生机理. 对于湍流来流, 叶尖涡主要存在于近尾迹. 蜿蜒是远尾迹的主要特征, 影响下游风力机的来流特征. 尾迹蜿蜒的产生有两种机制: 来流大尺度涡和剪切层失稳. 数值和观测结果显示两种机制共同存在. 机舱和中心涡对尾迹蜿蜒有重要影响. 采用叶片和机舱的致动面模型可准确预测尾迹蜿蜒. 研究显示不同风力机尾迹间的湍流特征存在相似性, 为发展尾迹湍流的快速预测模型提供了理论依据. 当前研究多关注平坦地形上的风力机尾迹, 复杂地形和海洋环境下的大气湍流和风力机尾迹的机理复杂, 现有工程模型无法准确预测, 有待深入研究.      

Simple analytical model for depth-averaged velocity in meandering compound channels

A simple but applicable analytical model is presented to predict the lateral distribution of the depth-averaged velocity in meandering compound channels. The governing equation with curvilinear coordinates is derived from the momentum equation and the flow continuity equation under the condition of quasi-uniform flow. A series of experiments are conducted in a large-scale meandering compound channel. Based on the experimental data, a magnitude analysis is carried out for the governing equation, and two lower-order shear stress terms are ignored. Four groups of experimental data from different sources are used to verify the predictive capability of this model, and good predictions are obtained. Finally, the determination of the velocity parameter and the limitation of this model are discussed.     

The loss of stability of laminar flow in open channel and the mechanism of sand ripple formation

ＩｎｔｒｏｄｕｃｔｉｏｎＤｕｒｉｎｇｔｈｅｔｉｍｅｏｆｂｅｄｓｕｒｆａｃｅｓｅｄｉｍｅｎｔｔｒａｎｓｐｏｒｔａｔｉｏｎ,ｓｅｄｉｍｅｎｔｐａｒｔｉｃｌｅｓｏｆｂｅｄｌｏａｄｍａｋｅｔｈｅｃｏｌｌｅｃｔｉｖｅｍｏｔｉｏｎｏｆｖａｒｉｏｕｓｄｉｆｆｅｒｅｎｔｆｏｒｍｓｏｎｔｈｅｒｉｖｅｒｂｅｄｓｕｒｆａｃｅ,ｔｈｕｓｔｈｅｆｏｒｍｏｆｂｅｄｓｕｒｆａｃｅｃｈａｎｇｅｓｃｏｎｓｔａｎｔｌｙ .Ｔｈｉｓｋｉｎｄｏｆｃｏｌｌｅｃｔｉｖｅｍｏｔｉｏｎｏｆｓｅｄｉｍｅｎｔｐａｒｔｉｃｌｅｓｏｎｔｈｅｂｅｄｓｕｒｆ…     

Three‐dimensional numerical simulation of compound meandering open channel flow by the Reynolds stress model

Turbulent flow in a compound meandering open channel with seminatural cross sections is one of the most complicated turbulent flows as the flow pattern is influenced by the combined action of various forces, such as centrifugal force, pressure, and shear stresses. In this paper, a three‐dimensional (3D) Reynolds stress model (RSM) is adopted to simulate the compound meandering channel flows. Governing equations of the flow are solved numerically with finite‐volume method. The velocity fields, wall shear stresses, and Reynolds stresses are calculated for a range of input conditions. Good agreement between the simulated results and measurements indicates that RSM can successfully predict the complicated flow phenomenon. Copyright © 2008 John Wiley & Sons, Ltd.     

An Open Channel Flow Experimental and Theoretical Study of Resistance and Turbulent Characterization over Flexible Vegetated Linings

Hydraulic engineers and scientists working on river restoration recognize the needfor a deeper understanding of natural streams as a complex and dynamicsystem, which involves not only abiotic elements (flow, sediments) but alsobiotic or biological components. From this point of view, the role played byriverine vegetation in river dynamics and flow conditions becomesessential. Hydro-mechanic interaction between flow and flexible plantscovering a river bed is studied in this paper and some previous works arediscussed. Experimental tests and measurements of turbulence on the flow in anopen channelwere performed using plastic plants seeded in a gravel bed. Characterization of flowresistance (friction factors) due to vegetation flexible roughness fordifferent plant densities was attained; furthermore, measuring detailedturbulent velocityprofiles within and above submerged and flexed stems allowed us to distinguishdifferent turbulent regimes. Some interesting relationships wereobtained between the velocity field and the deflected height of the plants, suchas a linear fit between the non-dimensional flexural parameter and the relativedeflection of the plants. Turbulent stresses weremeasured showing two different regions: above and inside the vegetationdomain. The spectral interaction between the plant oscillations, their wakes and theturbulence at different heights, forces strongly anisotropic Reynolds tensors andin order to clarify turbulent processes and their complexstructure, theoretical concepts (Taylor, Kolmogorov's K41) and several dataanalysis (autocorrelation functions, integral scales) were applied.     

Disturbance growth rate in turbulent wall flows

Disturbance development in turbulent wall flows is numerically investigated. The flows in a circular tube and in a plane channel are considered. The Navier-Stokes equations subjected to the condition of periodicity along the main flow are integrated in time until a statistically stationary “turbulent” flow regime is attained. Then the solution is disturbed and the further evolution of the disturbance is determined by comparing the two solutions, i.e., with and without the disturbance, which are calculated in parallel. It is shown that in the linear stage on average the solutions diverge exponentially. The main result of the study is that the small disturbance growth rate normalized by the wall time scale turns out to be constant, that is, dependent on neither the Reynolds number on the range considered nor the type of the flow: λ⁺ ≈ 0.021. The estimate of the disturbance growth rate is consistent with the previously obtained results concerning downstream disturbance growth and the estimate for the highest Lyapunov exponent calculated for turbulent flow in a plane channel.     

磨西河流域泥石流流域形态的非线性特征

泥石流是磨西河流域地质环境演化的一个重要地表过程,研究该流域内泥石流的地貌发育及演化特征可为磨西河流域第四纪以来地质环境的演变提供重要信息。通过对流域51条泥石流沟地貌形态的统计分析发现,磨西河流域内泥石流的分布在空间上表现出明显的不均性,流域内泥石流沟的数量与对应的主沟长度、主沟纵比降及流域面积之间存在着明显的分形现象,具有较高的自相似性。结果表明,泥石流沟谷形态的非线性演化特征及现象是对磨西河流域冰川环境和地质构造环境响应的直接结果,主要体现在现代冰川对流域的强烈侵蚀和贡嘎山隆升而导致流域的强烈下切。研究结果可为流域内泥石流灾害的形成机理、防灾减灾及流域环境演化提供理论依据。