Analysis of bedload transport in laminar flows

Bedload transport generally depends on the bed shear stress and Reynolds number. Many studies conducted for the condition of turbulent flows have revealed the dependence of the transport rate on the bed shear stress, while knowledge of the Reynolds number effect on the transport rate is very limited. As an extreme case to reflect the viscous effect on sediment transport, sediment transport in laminar flows is considered in this paper. A stochastic approach is adopted to explore how the transport rate can be associated with characteristics of laminar flows. First, the probability of erosion in the absence of turbulence is assumed to depend only on the randomness of bed particles. The probability is then applied to formulate the sediment transport rate, of which the derivation is made largely based on Einstein’s bedload theory. The theoretical result indicates that the dimensionless transport rate for laminar flows is dependent on the dimensionless shear stress and dimensionless particle diameter or the shear Reynolds number. Comparisons are finally made between the derived formula and an empirical correlation available in the literature.     

Power-law velocity profile in turbulent boundary layers: An integral reynolds-number dependent solution

Geophysical flows of practical interest encompass turbulent boundary layer flows. The velocity profile in turbulent flows is generally described by a log- or a power-law applicable to certain zones of the boundary layer, or by wall-wake law for the entire zone of the boundary layer. In this study, a novel theory is proposed from which the power-law velocity profile is obtained for the turbulent boundary layer flow. The new power-law profile is based on the conservation of mass and the skin friction within the boundary layer. From the proposed theory, analytical expressions for the power-law velocity profile are presented, and their Reynolds-number dependency is highlighted. The velocity profile, skin friction coefficient and boundary layer thickness obtained from the proposed theory are validated by the reliable experimental data for zero-pressure gradient turbulent boundary layers. The expressions for Reynolds shear stress and eddy viscosity distributions across the boundary layer are also obtained and validated by the experimental data.     

The wave-induced boundary layer under long internal waves

The boundary layer formed under the footprint of an internal solitary wave is studied by numerical simulation for waves of depression in a two-layer model of the density stratification. The inviscid outer flow, in the perspective of boundary-layer theory, is based on an exact solution for the long wave-phase speed, yielding a family of fully nonlinear solitary wave solutions of the extended Korteweg–de Vries equation. The wave-induced boundary layer corresponding to this outer flow is then studied by means of simulation employing the Reynolds-averaged Navier–Stokes (RANS) formulation coupled with a turbulence closure model validated for wall-bounded flows. Boundary-layer characteristics are computed for an extensive range of environmental conditions and wave amplitudes. Boundary-layer transition, identified by monitoring the eddy viscosity, is correlated in terms of a boundary-layer Reynolds number. The frictional drag is evaluated for laminar, transitional, and turbulent cases, and correlations are presented for the friction coefficient plus relevant measures of the boundary-layer thickness.     

THE THEORETICAL STUDY ON THE LAWS OF DRAG REDUCTION BY AERATION IN OPEN CHANNEL 1

ＴＨＥＴＨＥＯＲＥＴＩＣＡＬＳＴＵＤＹＯＮＴＨＥＬＡＷＳＯＦＤＲＡＧＲＥＤＵＣＴＩＯＮＢＹＡＥＲＡＴＩＯＮＩＮＯＰＥＮＣＨＡＮＮＥＬ１ＺＨＥＮＧＹｏｎｇｇａｎｇ２,ＦＡＮＧＤｕｏ３ａｎｄＬＩＧｕｉｆｅｎ４ＡＢＳＴＲＡＣＴＴｈｅｌａｗｓｏｆｄｒａｇｒｅ...     

INITIATION OF SEDIMENT MOTION IN LAMINAR OVERLAND FLOW

IINTRODUCTIONWhileriverflowsareusuallydeepandturbulent,overlandflowisextremelyshallowandcanbelaminar,transitionalandturbulent.Becauseoftheshallownessoftheflolw,overlandflowhydraulicsisgreatlyaffectedbysurfaceroughness,raindropimpact,andinthecaseoflaminarflow,flui(Iviscosity.Theinitiationofsedimentmovementinoverlandflowisthereforeexpectedtodifferfromthatinriverflows.InriverstUdies,bedshearStressgbhastraditionallybeenusedtocharacterizethecriticalflowconditionatwhichsedimentbeginstomove.At…     

A reduced 3D hydrodynamic model of a shallow,long, and weakly curved stream

A reduced three-dimensional mathematical model of a free-flow stream in nondeformable channels (a model of a long shallow flow), proposed earlier, has been studied analytically and numerically. The reduced model has been verified by direct numerical simulation of the flow by full hydrodynamic models in COMSOL finite-element software complex for laminar and turbulent flows of a viscous fluid. The obtained results show that the proposed reduced model of a shallow weakly curved stream flow adequately describes its hydrodynamics, so it can be used in systems of complex simulation of the ecology of water objects and the use of water resources.     

Structure of numerically simulated katabatic and anabatic flows along steep slopes

Direct numerical simulation (DNS) is applied to investigate properties of katabatic and anabatic flows along thermally perturbed (in terms of surface buoyancy flux) sloping surfaces in the absence of rotation. Numerical experiments are conducted for homogeneous surface forcings over infinite planar slopes. The simulated flows are the turbulent analogs of the Prandtl (1942) one-dimensional laminar slope flow. The simulated flows achieve quasi-steady periodic regimes at large times, with turbulent fluctuations being modified by persistent low-frequency oscillatory motions with frequency equal to the product of the ambient buoyancy frequency and the sine of the slope angle. These oscillatory wave-type motions result from interactions between turbulence and ambient stable stratification despite the temporal constancy of the surface buoyant forcing. The structure of the mean-flow fields and turbulence statistics in simulated slope flows is analyzed. An integral dynamic similarity constraint for steady slope/wall flows forced by surface buoyancy flux is derived and quantitatively verified against the DNS data.     

Turbulence structures in non-uniform flows

This study investigates turbulence structures in steady and non-uniform flows. Equations of Reynolds shear stress and turbulent velocity fluctuations are derived and their physical interpretations are explained. The theoretical results show that, different from previous studies, the variation of water surface can generate the wall-normal velocity, resulting in deviations of Reynolds shear stress and turbulence intensities from those in uniform flows. A self-similarity relationship is found between the Reynolds shear stress and turbulence intensities in non-uniform flows. The existence of self-similarity indicates that the effect of non-uniformity does not influence the mixing length. An empirical equation has been proposed to express the relationship based on experimental data available in the literature. Good agreement is achieved between the measured and predicted turbulence intensities by applying the self-similarity relationship.     

Experimental investigation on local shear stress and turbulence intensities over a rough non-uniform bed with and without sediment using 2D particle image velocimetry

The current work focuses on locally resolving velocities,turbulence,and shear stresses over a rough bed with locally non-uniform character.A nonporous subsurface layer and fixed interfacial sublayer of gravel and sand were water-worked to a nature-like bed form and additionally sealed in a hydraulic flume.Two-dimensional Particle Image Velocimetry(2 D-PIV) was applied in the vertical plane of the experimental flume axis.Runs with clear water and weak sediment transport were done under slightly supercritical flow to ensure sediment transport conditions without formation of considerable sediment deposits or dunes.The study design included analyzing the double-averaged flow parameters of the entire measurement domain and investigating the flow development at 14 consecutive vertical subsections.Local geometrical variabilities as well the presence of sediment were mainly reflected in the vertical velocity component.Whereas the vertical velocity decreased over the entire depth in presence of sediment transport,the streamwise velocity profile was reduced only within the interfacial sublayer.In the region with decelerating flow conditions,however,the streamwise velocity profile systematically increased along the entire depth extent.The increase in the main velocity(reduction of flow resistance)correlated with a decrease of the turbulent shear and main normal stresses.Therefore,effects of rough bed smoothening and drag force reduction were experimentally documented within the interfacial sublayer due to mobile sediment.Moreover,the current study leads to the conclusion that in nonuniform flows the maximum Reynolds stress values are a better predictor for the bed shear stress than the linearly extrapolated Reynolds stress profile.This is an important finding because,in natural flows,uniform conditions are rare.     

Comparing different methods of bed shear stress estimates in simple and complex flow fields

Bed shear stress is a fundamental variable in river studies to link ?ow conditions to sediment transport. It is, however, dif?cult to estimate this variable accurately, particularly in complex ?ow ?elds. This study compares shear stress estimated from the log pro?le, drag, Reynolds and turbulent kinetic energy (TKE) approaches in a laboratory ?ume in a simple boundary layer, over plexiglas and over sand, and in a complex ?ow ?eld around de?ectors. Results show that in a simple boundary layer, the log pro?le estimate is always the highest. Over plexiglas, the TKE estimate was the second largest with a value 30 per cent less than the log estimate. However, over sand, the TKE estimate did not show the expected increase in shear stress. In a simple boundary layer, the Reynolds shear stress seems the most appropriate method, particularly the extrapolated value at the bed obtained from a turbulent pro?le. In a complex ?ow ?eld around de?ectors, the TKE method provided the best estimate of shear stress as it is not affected by local streamline variations and it takes into account the increased streamwise turbulent ?uctuations close to the de?ectors. It is suggested that when single‐point measurements are used to estimate shear stress, the instrument should be positioned close to 0·1 of the ?ow depth, which corresponds to the peak value height in pro?les of Reynolds and TKE shear stress. Copyright © 2004 John Wiley & Sons, Ltd.     

An analysis of the characteristics of rough bed turbulent shear stresses in an open channel

Entrainment of sediment particles from channel beds into the channel flow is influenced by the characteristics of the flow turbulence which produces stochastic shear stress fluctuations at the bed. Recent studies of the structure of turbulent flow has recognized the importance of bursting processes as important mechanisms for the transfer of momentum into the laminar boundary layer. Of these processes, the sweep event has been recognized as the most important bursting event for entrainment of sediment particles as it imposes forces in the direction of the flow resulting in movement of particles by rolling, sliding and occasionally saltating. Similarly, the ejection event has been recognized as important for sediment transport since these events maintain the sediment particles in suspension. In this study, the characteristics of bursting processes and, in particular, the sweep event were investigated in a flume with a rough bed. The instantaneous velocity fluctuations of the flow were measured in two-dimensions using a small electromagnetic velocity meter and the turbulent shear stresses were determined from these velocity fluctuations. It was found that the shear stress applied to the sediment particles on the bed resulting from sweep events depends on the magnitude of the turbulent shear stress and its probability distribution. A statistical analysis of the experimental data was undertaken and it was found necessary to apply a Box-Cox transformation to transform the data into a normally distributed sample. This enabled determination of the mean shear stress, angle of action and standard error of estimate for sweep and ejection events. These instantaneous shear stresses were found to be greater than the mean flow shear stress and for the sweep event to be approximately 40 percent greater near the channel bed. Results from this analysis suggest that the critical shear stress determined from Shield's diagram is not sufficient to predict the initiation of motion due to its use of the temporal mean shear stress. It is suggested that initiation of particle motion, but not continuous motion, can occur earlier than suggested by Shield's diagram due to the higher shear stresses imposed on the particles by the stochastic shear stresses resulting from turbulence within the flow.     

Turbulence modeling of compound open-channel flows with and without vegetation on the floodplain using the Reynolds stress model

A Reynolds stress model for the numerical simulation of compound open-channel flows with vegetation on the floodplain is described. The Reynolds stress model consists of various sub-models such as Speziale et al.’s model, Mellor and Herring’s model, and Rotta’s model for the pressure–strain correlation term, the turbulent diffusion term, and the dissipation term, respectively. For validation of the model, plain compound open-channel flows are simulated. The computed results were compared with measured data by [Tominaga A, Nezu I. Turbulent structure in compound open-channel flows. J Hydraul Eng, ASCE 1991;117(1):21–41] and the results show that the Reynolds stress model successfully simulates the mean flow and turbulence structure of plain compound channel flows. The model was then applied to compound open-channel flows with vegetated floodplains. Good agreement between the simulated results and data from an algebraic stress model by [Naot D, Nezu I, Nakagawa H. Hydrodynamic behavior of partly vegetated open channels. J Hydraul Eng, ASCE 1996;122(11):625–33] was found. However, it was shown that the RSM is capable of predicting the velocity dip and lateral shift in the maximum streamwise velocity, which were not observed in the data from algebraic stress modeling. Finally, a depth-averaged analysis of the streamwise momentum equation was performed to investigate the lateral momentum transfer in compound channel flows with vegetated floodplains. Compared with components by the secondary currents and Reynolds stress, the drag force due to the presence of vegetation appears to be a factor in reducing the bottom shear stress in both main channel and floodplain.     

ON THE NECESSITY OF APPLYING A ROTATION TO INSTANTANEOUS VELOCITY MEASUREMENTS IN RIVER FLOWS

In studies on river channel flow turbulence, it is often the case that the measured mean vertical velocity is different from zero, indicating that the frame of reference of the current meter is not parallel to the flow streamline. This situation affects the estimate of Reynolds shear stress in the streamwise and vertical planes and consequently the analysis of the flow turbulent structure. One way to solve this problem is to correct data by applying a rotation and this is reviewed in the first part of the paper. However, in fluvial geomorphology, the studied flow is often complex and streamlines may exhibit significant changes from one point of measurement to the other. In this context, applying a rotation complicates the situation more than it simplifies it. The second part of this paper examines the question of velocity data correction in complex flows using a field example of the turbulent boundary layer over a very rough gravel bed and a laboratory example taken from flow at a river channel confluence. In both cases, velocity vectors are spatially variable. In the first case, errors in the Reynolds shear stress estimates are relatively low (ranging from −13 to 7 per cent/deg) while in the second case, they are much larger (−200 to 164 per cent/deg). The significance of these errors on the interpretation of turbulence statistics in river channel flows is discussed. We propose that corrections should be applied in all clear cases of sensor misalignment and when the frame of reference changes spatially and temporally. However, no corrections should be used where different flow velocity vector orientations, not sensor misalignment, are responsible for the mean vertical velocity differing from zero.     

LABORATORY EXPERIMENTS ON NAVIGATION-INDUCED BED SHEAR STRESSES AND SEDIMENT RESUSPENSION

IINTRODUCTIONDependingonflowandoperatingconditions,navigationtrafficmaycausesignificantresuspensionofdepositedsediment.Jnanumberofsituationsresuspensionofdepositedsedimentcanhavesevereenvironmentalrepercussions.Forinstance,ifthesedimentcontainscontaminants,thecontaminantsmaybereentrainedwiththesediment,taintingthewaterquality(Erdmannetal.,1994).Inothersituations,..evedincreasesintheamountofcleansuspendedsedimentcanbedetrimentalforaquaticplantsandanimals(Garcfaetal.,1998).Inordertoassesst…     

Homogeneous planar and two-dimensional mean-field antidynamo theorems with zero mean flow

In an electrically conducting fluid, two types of turbulence with a preferred direction are distinguished: planar turbulence, in which every velocity in the turbulent ensemble of flows has no component in the given direction; and two-dimensional turbulence, in which every velocity in the turbulent ensemble is invariant under translation in the preferred direction. Under the additional assumptions of two-scale and homogeneous turbulence with zero mean flow, the associated magnetohydrodynamic alpha- and beta-effects are derived in the second-order correlation approximation (SOCA) when the electrically conducting fluid occupies all space. Limitations of the SOCA are well known, but alpha- and beta-effects of a turbulent flow are useful in interpreting the dynamo effects of the turbulence. Two antidynamo theorems, which establish necessary conditions for dynamo action, are shown to follow from the special structures of these alpha- and beta-effects. The theorems, which are analogues of the laminar planar velocity and two-dimensional antidynamo theorems, apply to all turbulent ensembles with the prescribed alpha- and beta-effects, not just the planar and two-dimensional ensembles. The mean magnetic field is general in the planar theorem but only two-dimensional in the two-dimensional theorem. The two theorems relax the previous restriction to turbulence which is both two-dimensional and planar. The laminar theorems imply decay of the total magnetic field for any velocity of the associated turbulent ensemble. However, the mean-field theorems are not fully consistent with the laminar theorems because further conditions beyond those arising from the turbulence must be imposed on the beta-effect to establish decay of the mean magnetic field. In particular, negative turbulent magnetic diffusivities must be restricted. It is interesting that there is no inconsistency in the alpha-effects. The failure of the SOCA with the two-scale approximation to simply preserve the laminar antidynamo theorems at the beta-effect level is a further demonstration of the restricted validity of the theory and shows that negative diffusivity effects derived by approximation methods must be treated cautiously.     

Variable and conflicting shear stress estimates inside a boundary layer with sediment transport

This paper presents a comparison between two methods for estimating shear stress in an atmospheric internal boundary layer over a beach surface under optimum conditions, using wind velocities measured synchronously at 13 heights over a 1.7 m vertical array using ultrasonic anemometry. The Reynolds decomposition technique determines at‐a‐point shear stresses at each measurement height, while the Law‐of‐the‐Wall yields a single boundary layer estimate based on fitting a logarithmic velocity profile through the array data. Analysis reveals significant inconsistencies between estimates derived from the two methods, on both a whole‐event basis and as time‐series. Despite a near‐perfect fit of the Law‐of‐the‐Wall, the point estimates of Reynolds shear stress vary greatly between heights, calling into question the assumed presence of a constant stress layer. A comparison with simultaneously measured sediment transport finds no relationship between transport activity and the discrepancies in shear stress estimates. Results do show, however, that Reynolds shear stress measured nearer the bed exhibits slightly better correlation with sand transport rate. The findings serve as a major cautionary message to the interpretation and application of single‐height measurements of Reynolds shear stress and their equivalence to Law‐of‐the‐Wall derived estimates, and these concerns apply widely to boundary layer flows in general. © 2015 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Numerical model for the generation of the ensemble of lithospheric plates and their penetration through the 660-km boundary

In the kinematic theory of lithospheric plate tectonics, the position and parameters of the plates are predetermined in the initial and boundary conditions. However, in the self-consistent dynamical theory, the properties of the oceanic plates (just as the structure of the mantle convection) should automatically result from the solution of differential equations for energy, mass, and momentum transfer in viscous fluid. Here, the viscosity of the mantle material as a function of temperature, pressure, shear stress, and chemical composition should be taken from the data of laboratory experiments. The aim of this study is to reproduce the generation of the ensemble of the lithospheric plates and to trace their behavior inside the mantle by numerically solving the convection equations with minimum a priori data. The models demonstrate how the rigid lithosphere can break up into the separate plates that dive into the mantle, how the sizes and the number of the plates change during the evolution of the convection, and how the ridges and subduction zones may migrate in this case. The models also demonstrate how the plates may bend and break up when passing the depth boundary of 660 km and how the plates and plumes may affect the structure of the convection. In contrast to the models of convection without lithospheric plates or regional models, the structure of the mantle flows is for the first time calculated in the entire mantle with quite a few plates. This model shows that the mantle material is transported to the mid-oceanic ridges by asthenospheric flows induced by the subducting plates rather than by the main vertical ascending flows rising from the lower mantle.     

High-resolution numerical simulation of turbulence in natural waterways

We develop an efficient and versatile numerical model for carrying out high-resolution simulations of turbulent flows in natural meandering streams with arbitrarily complex bathymetry. The numerical model solves the 3D, unsteady, incompressible Navier-Stokes and continuity equations in generalized curvilinear coordinates. The method can handle the arbitrary geometrical complexity of natural streams using the sharp-interface curvilinear immersed boundary (CURVIB) method of Ge and Sotiropoulos (2007) [1]. The governing equations are discretized with three-point, central, second-order accurate finite-difference formulas and integrated in time using an efficient, second-order accurate fractional step method. To enable efficient simulations on grids with tens of millions of grid nodes in long and shallow domains typical of natural streams, the algebraic multigrid (AMG) method is used to solve the Poisson equation for the pressure coupled with a matrix-free Krylov solver for the momentum equations. Depending on the desired level of resolution and available computational resources, the numerical model can either simulate, via direct numerical simulation (DNS), large-eddy simulation (LES), or unsteady Reynolds-averaged Navier-Stokes (URANS) modeling. The potential of the model as a powerful tool for simulating energetic coherent structures in turbulent flows in natural river reaches is demonstrated by applying it to carry out LES and URANS in a 50-m long natural meandering stream at resolution sufficiently fine to capture vortex shedding from centimeter-scale roughness elements on the bed. The accuracy of the simulations is demonstrated by comparisons with experimental data and the relative performance of the LES and URANS models is also discussed.     

LES of turbulent surface shear stress and pressure-gradient-driven flow on shallow continental shelves

Turbulent shear flows on shallow continental shelves (here shallow means that the interaction with the solid, no-slip bottom is important) are of great importance because tide- and wind-driven flows on the shelf are drivers of the transfer of momentum, heat, and mass (gas) across the air–sea interface. These turbulent flows play an important role because vertical mixing and current are vectors for the transport of sediment and bioactive material on continental shelves. Understanding the dynamics of this class of flows presents complications because of the presence of a free surface and also because the flow can be driven by a pressure gradient (a tidal current), a stress at the free surface (a wind-driven current), or a combination of both. In addition, the flow can be modified by the presence of a wave field that can induce Langmuir circulation (Langmuir, Science 87:119–123, 1938). Large eddy simulation is used to quantify the effects of pressure gradient and wind shear on the distinctive structures of the turbulent flow. From these computations, an understanding of the physics governing the turbulence of pressure-driven and wind-driven flows, how they can interact in a normal or a tangential direction, and the effect of wave forcing on these flows is obtained.