Prediction of velocity distribution in straight open-channel flow with partial vegetation by singular perturbation method

A numerical analysis model based on two-dimensional shallow water differential equations is presented for straight open-channel flow with partial vegetation across the channel. Both the drag force acting on vegetation and the momentum exchange between the vegetation and non-vegetation zones are considered. The depth-averaged streamwise velocity is solved by the singular perturbation method, while the Reynolds stress is calculated based on the results of the streamwise velocity. Comparisons with the experimental data indicate that the accuracy of prediction is significantly improved by introducing a term for the secondary current in the model. A sensitivity analysis shows that a sound choice of the secondary current intensity coefficient is important for an accurate prediction of the depth-averaged streamwise velocity near the vegetation and non-vegetation interfaces, and the drag force coefficient is crucial for predictions in the vegetation zone.     

Flow characteristics of rectangular open channels with compound vegetation roughness

An idealized parallel flow caused by a lateral bed roughness difference due to the partial vegetation across a channel is investigated. Similar to the flow in a compound channel, there are mixing layers adjacent to the interface between the vegetation and the non-vegetation lanes, and a lateral momentum exchange occurs between the slow-moving water in the former lane and the fast-moving water in the latter lane. Under a uniform flow condition, the three-dimensional (3D) instantaneous velocities of two cases with different discharges and water depths are measured with a 16MHz acoustic Doppler velocimeter (micro ADV). The longitudinal variation of the streamwise velocity and the vertical variation of the Reynolds stress are analyzed. A quadrant analysis is carried out to investigate the outward and inward interaction, ejection, and sweep phenomenon caused by the vegetation variation across the channel. The results show that the flow characteristics in the vegetation lane are similar to those in an open channel fully covered with submerged vegetation, and the flow characteristics in the smooth non-vegetation lane are similar to those in a free open channel. For the cases studied here, the width of the mixing region is about 10% of the channel width, and the mixing region is mainly on the non-vegetation half.     

Destabilising effects of plant flexibility in air and aquatic vegetation canopy flows

With an analytical model coupling a mixing layer flow with an oscillating vegetation canopy through a drag force, we show that for both the cases of wind over a crop field and a water stream over aquatic plants, the dynamics of the plants is responsible for increasing the growth rate of the coupled instability in the lock-in range through two mechanisms. Because the flow and the vegetation canopy move in phase and thus minimise their interactions, the drag dissipation is decreased. The correlation between the two components of the perturbation flow velocity in the mixing layer is increased and makes the perturbation flow more efficient at extracting energy from the mean flow.     

Experimental Investigation on Turbulent Structure of Drag Reducing Channel Flow with Blowing Polymer Solution from the Wall

Turbulent drag reducing flow with blowing polymer solution from the channel wall was investigated experimentally using particle image velocimetry (PIV). Experiments were carried out with varying conditions of blowing polymer solution (e.g. weight concentration of polymer solution). Reynolds number based on the channel height and mean velocity was set to 20000 and 40000. When the polymer solution was blown from the channel wall, streamwise velocity fluctuation little increased, but wall-normal velocity fluctuation, Reynolds shear stress and correlation coefficient decreased significantly only near the blower wall. This behavior corresponds to the decrease of the ejection and sweep in the near-wall region observed by the investigation of instantaneous velocity map. On the contrary, this characteristic behavior was not observed at a position away from the blower wall (y/(H/2) > 0.4) and the scatter plot was almost the same as that of the water flow in this region. These results suggest that there are two regions in the drag reducing flow with blowing polymer solution from the wall; one is a non-Newtonian region which exists near the blower wall, and the other is a Newtonian region at a distance from the wall. The non-Newtonian region plays a key role in the drag reduction by the blowing polymer solution.     

Characteristic turbulent structure of a modified drag-reduced surfactant solution flow via dosing water from channel wall

The modification of the near-wall structure is very important for the control of wall turbulence. To ascertain the effect of near-wall modulation on the viscoelastic drag-reduced flow, the modified characteristics of a surfactant solution channel flow were investigated experimentally. The modulation was conducted on the boundary of the channel flow by injecting water from the whole surface of one side of the channel wall. The diffusion process of the injected water was observed by using the planar laser-induced fluorescence technique. The velocity statistics and characteristic structure including the spatial distributions of instantaneous streamwise velocity, swirling strength, and Reynolds shear stress were analyzed based on the velocity vectors acquired in the streamwise wall-normal plane by using the particle imaging velocimetry technique. The results indicated that the disturbance of the injected water was constricted within a finite range very near the dosing wall, and the Reynolds shear stress was increased in this region. However, the eventual drag reduction rate was found to be increased due to a relatively large decrement of viscoelastic shear stress in this near-wall region. Moreover, the flow structure under this modulation presented obvious regional characteristics. In the unstable disturbed region, the mixing of high-speed and low-speed fluids and the motions of ejection and sweep occurred actively. Many clockwise vortex cores were also found to be generated. This characteristic structure was similar to that in the ordinary turbulence of Newtonian fluid. Nevertheless, outside this disturbed region, the structure still maintained the characteristics of the drag-reduced flow with non-Newtonian viscoelastic additives. These results proved that the injected Newtonian fluid associated with the modified stress distribution creates a diverse characteristic structure and subsequent enhanced drag reduction. This investigation can provide the experimental basis for further study of turbulence control.     

The simulation of turbulent particle‐laden channel flow by the Lattice Boltzmann method

We perform direct numerical simulation of three‐dimensional turbulent flows in a rectangular channel, with a lattice Boltzmann method, efficiently implemented on heavily parallel general purpose graphical processor units. After validating the method for a single fluid, for standard boundary layer problems, we study changes in mean and turbulent properties of particle‐laden flows, as a function of particle size and concentration. The problem of physical interest for this application is the effect of water droplets on the turbulent properties of a high‐speed air flow, near a solid surface. To do so, we use a Lagrangian tracking approach for a large number of rigid spherical point particles, whose motion is forced by drag forces caused by the fluid flow; particle effects on the latter are in turn represented by distributed volume forces in the lattice Boltzmann method. Results suggest that, while mean flow properties are only slightly affected, unless a very large concentration of particles is used, the turbulent vortices present near the boundary are significantly damped and broken down by the turbulent motion of the heavy particles, and both turbulent Reynolds stresses and the production of turbulent kinetic energy are decreased because of the particle effects. We also find that the streamwise component of turbulent velocity fluctuations is increased, while the spanwise and wall‐normal components are decreased, as compared with the single fluid channel case. Additionally, the streamwise velocity of the carrier (air) phase is slightly reduced in the logarithmic boundary layer near the solid walls. Copyright © 2015 John Wiley & Sons, Ltd.     

Two dimensional analytical solution for a partially vegetated compound channel flow

The theory of an eddy viscosity model is applied to the study of the flow in a compound channel which is partially vegetated. The governing equation is constituted by analyzing the longitudinal forces acting on the unit volume where the effect of the vegetation on the flow is considered as a drag force item, The compound channel is divided into 3 sub-regions in the transverse direction, and the coefficients in every region＇s differential equations were solved simultaneously. Thus, the analytical solution of the transverse distribution of the depth-averaged velocity for uniform flow in a partially vegetated compound channel was obtained. The results can be used to predict the transverse distribution of bed shear stress, which has an important effect on the transportation of sediment. By comparing the analytical results with the measured data, the analytical solution in this paper is shown to be sufficiently accurate to predict most hydraulic features for engineering design purposes.     

Temporal and spatial transients in turbulent boundary layer flow over an oscillating wall

Direct numerical simulations have been performed to study the effect of an oscillating segment of the wall on a turbulent boundary layer flow. Two different oscillation amplitudes with equal oscillation period have been used, which allows a direct comparison between a relatively weak and strong forcing of the flow. The weaker forcing results in 18% drag reduction while the stronger forcing, with twice the amplitude, yields 29% drag reduction. The downstream development of the drag reduction is compared with earlier simulations and experiments. In addition, a simulation with identical oscillation parameters as in previous numerical and experimental investigations allows for an estimation of the effect of the Reynolds number on the drag reduction.Reductions in the Reynolds stresses and the important role that the edge of the Stokes layer has is explained.An estimation of the idealized power consumption shows that a positive energy budget is only possible for the weaker wall velocity case.Spatial and temporal transients are investigated and a transformation between spatial and temporal coordinates via a convection velocity is shown to facilitate a comparison between the two transients in a consistent manner. The streamwise shear exhibits a similar monotonic behavior in the spatial and temporal transients, while the non-monotinic temporal transient of the longitudinal Reynolds stress has no counterpart in the spatial development. Furthermore, the evolution in time of the spanwise Reynolds stress is very similar to previously reported channel flow data.The instantaneous spanwise velocity profile (only averaged in the homogeneous spanwise direction) will for the first time be presented from a boundary layer over an oscillating wall, and comparisons with the analytical solution to the laminar Navier–Stokes equations show very good agreement.     

Flow control over an undulating membrane

The flow along a flexible membrane forced to undulate in the form of a streamwise travelling wave pattern is studied experimentally in detail. Flow field and force measurements confirm that the form drag of the wavy wall is significantly reduced when starting the undulatory motion. A mechanical model of an undulating membrane was built, based on previous investigations described in literature, and placed in an open water channel. The motion pattern of the membrane was prescribed in such a way to achieve a downstream travelling wave with increasing amplitude. The exploratory focus laid on the identification of hydrodynamic mechanisms of drag reduction due to undulatory motion. The wave-speed c of the travelling wave was set proportional to the incoming flow velocity U, according to an optimum ratio identified by previous numerical and experimental investigations. Poisson’s equation for the pressure was used to calculate the 2D pressure field from the experimental data of the unsteady flow field. In addition, the integral drag force of the membrane, as a function of c/U, was measured with a force balance to compare with previous published numerical findings. Furthermore, the velocities close to the surface of the membrane were measured, and the boundary layer profiles were determined. The resulting normalised velocity profiles affirm an oscillation between laminar and turbulent flow over one period of the motion. The results are in good agreement with previous experimental and numerical findings. Additionally, the characteristics of the flow along a travelling wave with increasing amplitude are discussed in more detail.     

Velocity distribution of flow with submerged flexible vegetations based on mixing-length approach

By choosing a PVC slice to simulate flexible vegetation, we carried out experiments in an open channel with submerged flexible vegetation. A 3D acoustic Doppler velocimeter （micro ADV） was used to measure local flow velocities and Reynolds stress. The results show that hydraulic characteristics in non-vegetation and vegetation layers are totally different. In a region above the vegetation, Reynolds stress distribution is linear, and the measured velocity profile is a classical logarithmic one. Based on the concept of new-riverbed, the river compression parameter representing the impact of vegetation on river is given, and a new assumption of mixing length expression is made. The formula for time-averaged velocity derived from the expression requires less parameters and simple calculation, and is useful in applications.     

Experiments in Turbulent Pipe Flow with Polymer Additives at Maximum Drag Reduction

In this paper we report on (two-component) LDV experiments in a fully developed turbulent pipe flow with a drag-reducing polymer (partially hydrolyzed polyacrylamide) dissolved in water. The Reynolds number based on the mean velocity, the pipe diameter and the local viscosity at the wall is approximately 10000. We have used polymer solutions with three different concentrations which have been chosen such that maximum drag reduction occurs. The amount of drag reduction found is 60–70%. Our experimental results are compared with results obtained with water and with a very dilute solution which exhibits only a small amount of drag reduction. We have focused on the observation of turbulence statistics (mean velocities and turbulence intensities) and on the various contributions to the total shear stress. The latter consists of a turbulent, a solvent (viscous) and a polymeric part. The polymers are found to contribute significantly to the total stress. With respect to the mean velocity profile we find a thickening of the buffer layer and an increase in the slope of the logarithmic profile. With respect to the turbulence statistics we find for the streamwise velocity fluctuations an increase of the root mean square at low polymer concentration but a return to values comparable to those for water at higher concentrations. The root mean square of the normal velocity fluctuations shows a strong decrease. Also the Reynolds (turbulent) shear stress and the correlation coefficient between the stream wise and the normal components are drastically reduced over the entire pipe diameter. In all cases the Reynolds stress stays definitely non-zero at maximum drag reduction. The consequence of the drop of the Reynolds stress is a large polymer stress, which can be 60% of the total stress. The kinetic-energy balance of the mean flow shows a large transfer of energy directly to the polymers instead of the route by turbulence. The kinetic energy of the turbulence suggests a possibly negative polymeric dissipation of turbulent energy. This revised version was published online in July 2006 with corrections to the Cover Date.     

暖季强降雨对多年冻土南界斜坡路基稳定性影响分析

为阐明表面活性剂水溶液的减阻作用,使用LDV对零压梯度的二维湍流平板边界层中的CTAB 表面活性剂水溶液的湍流特性进行了实验研究. 结果表明：与牛顿流体相比,CTAB水溶液边 界层的粘性底层增厚;主流时均速度分布有被层流化的趋势,对数分布域上移;主流方向速 度湍动强度峰值减小,且远离壁面,在靠近边界层中部,出现第2峰值;垂直于主流方向的 速度湍动强度受到了大幅度抑制,雷诺应力沿着边界层厚度方向几乎为零. 结果说明CTAB 水溶液具有减弱湍流湍动各个成分相关度的作用,从而能够使雷诺应力降低、湍流能量生成 项减小最终降低流体的输送动力.     

固体颗粒对沟槽湍流边界层影响的实验研究

本文采用粒子图像测速技术(particles image velocimetry, PIV)研究固体颗粒对放置在平板湍流边界层中的平壁和沟槽壁面减阻效果的影响. 实验对清水和加入粒径为155 μm聚苯乙烯颗粒的流法向二维速度场信息进行采集, 对不同工况下的平均速度剖面、雷诺应力和湍流度等统计量进行对比, 分析流体在边界层中的行为. 运用空间局部平均结构函数提取了不同工况湍流边界层喷射?扫掠行为的空间拓扑结构并进行比较. 结果发现, 在不同的壁面条件下, 粒子加入后的对数律区中无量纲速度均略大于清水组, 雷诺切应力有所降低, 湍流度有所减弱. 对于不同流场速度下的沟槽而言, 颗粒的加入均降低了壁面附近的阻力, 而颗粒单独作用于光滑壁面的减阻效果并不明显. 加入粒子后的相干结构数目有所增加, 法向脉动速度下降. 沟槽壁面附近的相干结构数目有所增加, 法向脉动速度在自由来流速度较大时有所上升, 在速度较小时有所下降. 这表明不同减阻状况下的沟槽均能将大涡破碎成更多的涡, 并且粒子的加入强化了这种破碎作用.      

上海陆域古河道分布及对工程建设影响研究

采用直接数值模拟方法, 对槽道湍流中展向振动流向传播的波动壁面的流动 控制和减阻问题进行了研究, 讨论了流向参数k_{x}对Stokes层、湍流拟序结构、湍流猝 发事件以及壁面阻力的影响, 并对此类波动壁面的湍流控制和减阻机理进行了讨论. 结果表 明, 当此类波动壁面被用来调制近壁流动时, 仅低频波对湍流流场具有显著影响, 可导致湍 流猝发事件的频率和强度的显著变化; 波数k_{x}的增大对于湍流猝发事件的频率和强度增 减的影响并不同步, 存在一个最优的波数k_{x}, 在其调制下, 固有流场对诱导流场的影响 最弱, 而诱导流场对固有流场的影响显著, 减阻效果最好.     

Layered structure of turbulent plane wall jet

This paper investigates the layered structure of a turbulent plane wall jet at a distance from the nozzle exit. Based on the force balances in the mean momentum equation, the turbulent plane wall jet is divided into three regions: a boundary layer-like region (BLR) adjacent to the wall, a half free jet-like region (HJR) away from the wall, and a plug flow-like region (PFR) in between. In the PFR, the mean streamwise velocity is essentially the maximum velocity, and the simplified mean continuity and mean momentum equations result in a linear variation of the mean wall-normal velocity and Reynolds shear stress. In the HJR, as in a turbulent free jet, a proper scale for the mean wall-normal flow is the mean wall-normal velocity far from the wall and a proper scale for the Reynolds shear stress is the product of the maximum mean streamwise velocity and the velocity scale for the mean wall-normal flow. The BLR region can be divided into four sub-layers, similar to those in a canonical pressure-driven turbulent channel flow or shear-driven turbulent boundary layer flow. Building on the log-law for the mean streamwise velocity in the BLR, a new skin friction law is proposed for a turbulent wall jet. The new prediction agrees well with the correlation of Bradshaw and Gee (1960) over moderate Reynolds numbers, but gives larger skin frictions at higher Reynolds numbers.     

Considerations on shock wave/boundary layer interaction in undular hydraulic jumps in horizontal channels with a very high aspect ratio

It can be seen in the literature that the fundamental factors governing oblique shock wave development, typically in very large channels with straight sidewalls, have not yet been completely understood and remain at the level of indicating its presence and formation. In this study, in addition to an analysis of various properties of hydraulic jump behaviour in very large channels, some aspects of boundary layer development and its detachment from the channel lateral sidewall are also investigated. At the detachment point of the lateral shock waves, it was noted that the displacement thickness experiences a significant increase; this is accompanied by a significantly reduced gradient normal to the channel sidewalls of the flow velocity as well as the occurrence of a strong, sudden adverse pressure gradient. An analysis of the flow velocity distribution and the background turbulence intensity of both the streamwise and spanwise velocity components was also carried out. Furthermore, it is argued that the supersonic flow separation analogy with a supercritical free surface flow can be applied to this case study and that the behaviour of the supercritical flow during separation can be interpreted by the free interaction theory typically used in aerodynamics.     

带滩槽地形的连续弯道中纵向流速横向分布解析

本文基于沿水深积分的动量方程,假定二次流项和弯道附加应力项沿横断面呈线性分布,提出了预测弯道垂线平均纵向流速的解析计算方法,进一步提出了河槽区和河滩区垂线平均纵向流速沿断面分布的求解模式,并将其应用于带滩槽地形的反向连续弯道水槽中.根据实测数据率定计算参数,该模式可计算不同出口水深条件下断面垂线平均纵向流速分布,计算结...     

Development of a generalized multi‐layer model for 3‐D simulation of free surface flows

A mathematical model was developed for three‐dimensional (3‐D) simulation of free surface flows. In this model, the flow depth is divided into a number of layers and shallow water equations are integrated in each layer to derive the hydrodynamic equations. To give a general form to this model, each layer is assumed to be non‐horizontal with varying thickness in the flow domain. A non‐orthogonal curvilinear coordinate system is employed in the model, to allow for flexibility in dealing with the irregular geometry of natural watercourses. Due to the similarity in governing equations, two‐dimensional (2‐D) depth averaged programs can be developed into a multi‐layer model. The development for a depth averaged program and its numerical scheme is described in this paper. Experimental data and semi‐analytical solutions are used to evaluate the performance of the model. Three different cases of open channel flow are tested: 1‐flow in a straight open channel, 2‐the flow development region in a channel, and 3‐flow in a meandering channel. It is shown that the model has the capability to predict velocity distribution and secondary flows in complex 3‐D flow conditions. Copyright © 2004 John Wiley & Sons, Ltd.     

Spatio-temporal analysis of the turbulent flow in a ribbed channel

The turbulent flow in a channel with transverse ribs over one wall is studied experimentally. The height of the obstacles is about one tenth of the channel height, and the spacing is 10 times their height. The Reynolds number based on the channel hydraulic diameter and bulk flow velocity is 15,000. Velocity fields are obtained with high spatial and temporal resolution along the streamwise/wall-normal plane by means of time-resolved particle image velocimetry. Beside mean velocity and Reynolds stresses, the flow is investigated through two-point correlations, distributions of spanwise-swirling events, space–time velocity diagrams and power spectral density. Although the separated flow reattaches before the following obstacle is approached, a strong rib-to-rib interaction occurs. Spanwise vortices, 0.2 rib heights in size, are generated in the free shear layer, travel across the whole pitch, and may impact on the next rib. The large scale motions triggered by the separation grow in size until they reach the following obstacle. Flapping of the separated shear layer is observed at frequencies consistent with previous studies, causing the instantaneous reattachment point to fluctuate. The flapping initiates at the downstream edge of the obstacle tip, rather than at the upstream edge where the instantaneous separation occurs.     

Optimizing a fin ray for stiffness

We optimize the elastic constants of the supporting structure in flexible fish fins—the fin ray—to obtain the least deflection under loading. We first solve the problem numerically for rays with shear layers of uniform thickness. The optimal distributions of bending and shear moduli are nonzero on disjoint regions. The numerical solution suggests the form of the analytical solution, which we obtain using calculus of variations on two intervals with continuity conditions at the interface. The deflection of the optimal ray is less than half that of the uniform ray. When the shear layer thickness is allowed to vary, we find a tapered thickness distribution which gives zero deflection for a given positive loading.