Stream power criterion and design of sand bed channels at incipient motion

The critical stream power criterion may be used to describe the incipient motion of cohesionless particles of plane sediment beds. The governing equation relating ``critical stream power'' to ``shear Reynolds number'' is developed by using the present experimental data as well as the data from several other sources. Simultaneously, a resistance equation, relating the ``particle Reynolds number'' to the``shear Reynolds number'' is developed for plane sediment beds in wide channels with little or no transport. By making use of these relations, a procedure is developed to design plane sediment beds such that any two of the four design variables, including particle size, energy/friction slope, flow depth, and discharge per unit width in the channel should be known to predict the remaining two variables. Finally, a straightforward design procedure using design tables/design curves and analytical methods is presented to solve six possible design problems.

Analysis of gradually and spatially varied flow in sand-bed channels

Seepage through a sand bed affects the channel hydrodynamics, which in turn alters channel stability. Thus, the effect of seepage on its hydrodynamic parameters needs to be ascertained. The present work analyses spatially varied flow of a sand-bed channel subjected to seepage in the downward direction through a sand bed. Numerically calculated flow profiles affected by seepage have been verified using experimental observations. The present work also analyses the friction slope, velocity and bed shear stress variations along the channel for both seepage and no-seepage conditions. It was found that the downward seepage-induced channel flow has larger friction slope and bed shear stress than that of no-seepage.

A numerical study of division of flow in open channels

A two-dimensional numerical model which employs the depth-averaged forms of continuity and momentum equations along with k-e turbulence closure scheme is used to simulate the flow at the open channel divisions. The model is generalised to flows of arbitrary geometries and MacCormack finite volume method is used for solving governing equations. Application of cartesian version of the model to analyse the flow at right-angled junction is presented. The numerical predictions are compared with experimental data of earlier investigators and measurements made as part of the present study. Performance of the model in predicting discharge distribution, surface profiles, separation zone parameters and energy losses is evaluated and discussed in detail. To illustrate the application of the numerical model to analyse the flow in acute angled offtakes and streamlined branch entries, a few computational results are presented.

Spatial B-jump at sudden channel enlargements with abrupt drop

An analytical and experimental study of the hydraulic jump in stilling basins with abrupt drop and sudden enlargement, called the spatial B-jump here, is carried out for finding the sequent depth ratio and resulting energy dissipation. The spatial B-jump studied has its toe downstream of the expansion section, and the stream lines at the toe are characterized by downward curvature. An expression is obtained for the sequent depth ratio based on the momentum equation with suitable assumptions for the extra pressure force term because of the abrupt drop in the bed and sudden enlargement in the basin width. Predictions compare favorably with experiments. It is shown that the spatial B-jump needs less tailwater depth, thereby enhancing the stability of the jump when compared either with spatial jump, which forms in sudden expanding channels, or with B-jump, which forms in a channel with an abrupt drop in bed. It is also shown that there is a significant increase in relative energy loss for the spatial B-jump compared to either the spatial jump or B-jump alone.

Bed Shear from Velocity Profiles: A New Approach

A new binary law of velocity distribution has been developed to describe the velocity profile for the entire flow region. The law is a combination of logarithmic law, valid in the wall (inner) region, and parabolic law, valid in the core (outer) region of the flow. The validity of the law has been established based on earlier data on flat plates, rough and smooth pipes and experimental data obtained from rigid-walled open channels with plane sand beds. A procedure of estimating bed shear stress from the proposed law of velocity distribution using the measured velocity profile has been evolved. Bed shear estimates made according to this procedure are in agreement with the values obtained from uniform flow analysis in the case of open channel flow over a sediment bed. The proposed method of estimating the bed shear stress from the observed velocity profiles is found to be particularly useful in cases where it is difficult to determine precisely the true bed level, such as in the case of flow over sediment beds.

Transition of turbulent pipe flow

The velocity profile in turbulent pipe flow is usually divided into two regions, a wall or inner region and a core or outer region. For the inner region, the viscosity and wall shear stress are the important parameters governing the velocity distribution whereas for the outer region, the wall reduces the velocity below the maximum velocity independent of viscosity. In the present work, a velocity model is proposed for turbulent flow in the wall region of a pipe covering the entire transition from smooth to rough flows. Coupling this model for the wall region with the power law velocity model for the core region, an equation for the friction factor is obtained. The model constants are evaluated by using Nikuradse's experiments in the fully smooth and rough turbulent flows. The model shows good agreement with the friction factor and the velocity profiles obtained by Nikuradse for the transition region of turbulent flow.

Water balance modelling in a tropical watershed under deciduous forest (Mule Hole, India): Regolith matric storage buffers the groundwater recharge process

Accurate estimations of water balance are needed in semi-arid and sub-humid tropical regions, where water resources are scarce compared to water demand. Evapotranspiration plays a major role in this context, and the difficulty to quantify it precisely leads to major uncertainties in the groundwater recharge assessment, especially in forested catchments. In this paper, we propose to assess the importance of deep unsaturated regolith and water uptake by deep tree roots on the groundwater recharge process by using a lumped conceptual model (COMFORT). The model is calibrated using a 5 year hydrological monitoring of an experimental watershed under dry deciduous forest in South India (Mule Hole watershed). The model was able to simulate the stream discharge as well as the contrasted behaviour of groundwater table along the hillslope. Water balance simulated for a 32 year climatic time series displayed a large year-to-year variability, with alternance of dry and wet phases with a time period of approximately 14 years. On an average, input by the rainfall was 1090 mm year(-1) and the evapotranspiration was about 900 mm year(-1) out of which 100 mm year(-1) was uptake from the deep saprolite horizons. The stream flow was 100 mm year(-1) while the groundwater underflow was 80 mm year(-1). The simulation results suggest that (i) deciduous trees can uptake a significant amount of water from the deep regolith, (ii) this uptake, combined with the spatial variability of regolith depth, can account for the variable lag time between drainage events and groundwater rise observed for the different piezometers and (iii) water table response to recharge is buffered due to the long vertical travel time through the deep vadose zone, which constitutes a major water reservoir. This study stresses the importance of long term observations for the understanding of hydrological processes in tropical forested ecosystems. (C) 2009 Elsevier B.V. All rights reserved.

Geometrically simple linear weirs using circular quadrants

A detailed theoretical analysis of flow through a quadrant plate weir is made in the light of the generalized theory of proportional weirs, using a numerical optimization procedure. It is shown that the flow through the quadrant plate weir has a linear discharge-head relationship valid for certain ranges of head. It is shown that the weir is associated with a reference plane or datum from which all heads are reckoned.Further, it is shown that the measuring range of the quadrant plate weir can be considerably enhanced by extending the tangents to the quadrants at the terminals of the quadrant plate weir. The importance of this weir (when the datum of the weir lies below its crest) as an outlet weir for grit chambers is highlighted. Experiments show excellent agreement with the theory by giving a constant average coefficient of discharge.

On the dual head-discharge characteristics of a modified chimney weir

This paper is concerned with the dual head-discharge characteristics of a modified chimney weir. It is shown by an optimization procedure that the modified chimney weir having an inward trapezoidal weir over an inverted V-notch can produce discharges proportional to both the logarithm of the head as well as the linear power of the head reckoned independently over two different reference planes, within a prefixed maximum allowable percentage error from the theoretical discharge. A new technique is adopted to optimize the weir parameters, in order to obtain the maximum ranges of measurement under logarithmic, linear as well as combined characteristics. In the case of linear weir it is shown that it is possible to enhance the linearity range of the chimney weir by more than 540% and for a weir with constant indication accuracy by more than 350%. In addition, about 86% of the overall depth of the designed linear weir and over 90% of the logarithmic weir is converted as the corresponding measurable ranges. Experiments with four typical weirs give consistent constant average coefficient of discharge for each weir confirming the theory. The practical application of the weir in minor irrigation, hydraulic and other engineering fields is highlighted.

Geometry of sand-bed channels with seepage

Literature suggests that apart from a significant percentage of water loss from seepage, downward seepage from the channel causes an increase in the mobility of channel bed materials and thus changes the channel stability. Consequently, regime conditions (which provide the relationship among hydraulic parameters) should also be revised by incorporating downward seepage as an additional parameter. In the present work, regime conditions for the prediction of channel geometry in alluvial channels affected by downward seepage have been formulated on the basis of experimental observations. (C) 2011 Elsevier B.V. All rights reserved.

Effect of a magnetic field on the flow and blood oxygenation in channels of variable cross-section

The effect of a magnetic field on the flow and oxygenation of an incompressible Newtonian conducting fluid in channels with irregular boundaries has been investigated. The geometric parameter δ, which is a ratio of the mean half width of the channel d to the characteristic length λ along the channel over which the significant changes in the flow quantities occur, has been used for perturbing the governing equations. Closed form solutions of the various order equations are presented for the stream function. The equations for oxygen partial pressure remain nonlinear even after perturbation, therefore a numerical solution is presented. The expressions for shear stress at a wall and pressure distributions are derived. Here the separation in the flow occurs at a higher Reynolds number than the corresponding non-magnetic case. It is found that the magnetic field has an effect on local oxygen concentration but has a little effect on the saturation length.

A fuzzy dynamic flood routing model for natural channels

A fuzzy dynamic flood routing model (FDFRM) for natural channels is presented, wherein the flood wave can be approximated to a monoclinal wave. This study is based on modification of an earlier published work by the same authors, where the nature of the wave was of gravity type. Momentum equation of the dynamic wave model is replaced by a fuzzy rule based model, while retaining the continuity equation in its complete form. Hence, the FDFRM gets rid of the assumptions associated with the momentum equation. Also, it overcomes the necessity of calculating friction slope (S-f) in flood routing and hence the associated uncertainties are eliminated. The fuzzy rule based model is developed on an equation for wave velocity, which is obtained in terms of discontinuities in the gradient of flow parameters. The channel reach is divided into a number of approximately uniform sub-reaches. Training set required for development of the fuzzy rule based model for each sub-reach is obtained from discharge-area relationship at its mean section. For highly heterogeneous sub-reaches, optimized fuzzy rule based models are obtained by means of a neuro-fuzzy algorithm. For demonstration, the FDFRM is applied to flood routing problems in a fictitious channel with single uniform reach, in a fictitious channel with two uniform sub-reaches and also in a natural channel with a number of approximately uniform sub-reaches. It is observed that in cases of the fictitious channels, the FDFRM outputs match well with those of an implicit numerical model (INM), which solves the dynamic wave equations using an implicit numerical scheme. For the natural channel, the FDFRM Outputs are comparable to those of the HEC-RAS model.

Waves in branched hydraulic pipes

The paper deals with the classical problem of axi-symmetric transmission of low amplitude waves through a circular pipe containing a viscous liquid. Exact governing equations are identified and solved, the radial as well as the axial component of the velocity being considered. Attention is drawn to certain fallacies underlying the conventional approach. The parameters required in the formulation of the transfer matrix for a pipe have been evaluated. In order to evaluate the response at the terminal point of a branched system for a sinusoidal input at one of the ends, a general algorithm has been developed.