Longitudinal dispersion in laboratory and natural streams

This study concerns the longitudinal dispersion of fluid particles which are initially distributed uninformly over one cross section of a uniform, steady, turbulent open channel flow. The primary focus is on developing a method to predict the rate of dispersion in a natural stream.

Taylor's method of determining a dispersion coefficient, previously applied to flow in pipes and two-dimensional open channels, is extended to a class of three-dimensional flows which have large width-to-depth ratios, and in which the velocity varies continuously with lateral cross-sectional position. Most natural streams are included. The dispersion coefficient for a natural stream may be predicted from measurements of the channel cross-sectional geometry, the cross-sectional distribution of velocity, and the overall channel shear velocity. Tracer experiments are not required.

Large values of the dimensionless dispersion coefficient D/rU* are explained by lateral variations in downstream velocity. In effect, the characteristic length of the cross section is shown to be proportional to the width, rather than the hydraulic radius. The dimensionless dispersion coefficient depends approximately on the square of the width to depth ratio.

A numerical program is given which is capable of generating the entire dispersion pattern downstream from an instantaneous point or plane source of pollutant. The program is verified by the theory for two-dimensional flow, and gives results in good agreement with laboratory and field experiments.

Both laboratory and field experiments are described. Twenty-one laboratory experiments were conducted: thirteen in two-dimensional flows, over both smooth and roughened bottoms; and eight in three-dimensional flows, formed by adding extreme side roughness to produce lateral velocity variations. Four field experiments were conducted in the Green-Duwamish River, Washington.

Both laboratory and flume experiments prove that in three-dimensional flow the dominant mechanism for dispersion is lateral velocity variation. For instance, in one laboratory experiment the dimensionless dispersion coefficient D/rU* (where r is the hydraulic radius and U* the shear velocity) was increased by a factory of ten by roughening the channel banks. In three-dimensional laboratory flow, D/rU* varied from 190 to 640, a typical range for natural streams. For each experiment, the measured dispersion coefficient agreed with that predicted by the extension of Taylor's analysis within a maximum error of 15%. For the Green-Duwamish River, the average experimentally measured dispersion coefficient was within 5% of the prediction.

Eliminating discharge imbalances in predicting shallow water flows with shocks

The shallow water equations are widely used in modelling environmental flows. Being a hyperbolic system of differential equations, they admit shocks that represent hydraulic jumps and bores. Although the water surface can be solved satisfactorily with the modern shock-capturing schemes, the predicted flow rate often suffers from imbalances where shocks occur, eg the mass conservation is violated by failing to maintain a constant discharge rate at every cross-section in a steady open channel flow. A total-variation-diminishing Lax-Wendroff scheme is developed, and used to demonstrate how to achieve an exact flux balance. The performance of the proposed methods is inspected through some test cases, which include 1- and 2-dimensional, flat and irregular bed scenarios. The proposed methods are shown to preserve the mass exactly, and can be easily extended to other shock-capturing models.

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.

Effect of substrate surface roughness on electric current induced flow of liquid metals

Electric current can induce long-range flow of liquid metals over a conducting substrate. This work reports on the effect of the substrate surface roughness on the liquid metal-front velocity during such a flow. Experiments were conducted by passing electric current through liquid gallium placed over similar to 170 nm thick, 500 mu m wide gold and platinum films of varying roughness. The ensuing flow, thus, resembles micro-fluidics behavior in an open-channel. The liquid-front velocity decreased linearly with the substrate surface roughness; this is attributed to the reduction in the effective electric field along the liquid metal-substrate interface with the substrate surface roughness. (C) 2013 American Institute of Physics. http://dx.doi.org/10.1063/1.4790182]

Sidewall flow stability in the MHD channel flows

The velocity distribution between two sidewalls is; M-shaped for the MHD channel flows with rectangular cross section and thin conducting walls in a strong transverse magnetic field. Assume that the dimensionless numbers R(m) much less than 1, M, N much greater than 1, and sigma

The flow behaviour of shallow fluidised beds

The flow behaviour of shallow gas-fluidised beds was studied. experimentally using a rotational viscometer, and an inclined open channel. Initially, tests were carried out with the viscometer in order to establish qualitative trends in the flow properties of a variety of materials over a wide range of fluidising conditions. Also, a technique was developed which enabled quantitative viscosity data to be extracted from the experimental results. The flow properties were found to be sensitive to the size, size-range and density of the fluidised material, the type of distributor used, and the moisture content of the fluidising gas. Tests in beds up to 120 mm deep showed that the fluidity of the bed improves with reduction in depth; and indicated a range of flow behaviour from shear-thinning to Newtonian, depending chiefly on fluidising velocity .. Later, an apparatus was built which provided for a steady, continuous flow of fluidised material down an inclined open channel of 3m length x 0.15m square, up to a mass flowrate of 10 kg/s (35 ton/hr). This facility has enabled data to be obtained that is of practical value in industrial applications; which is otherwise difficult in view of the present limited understanding of the true mechanism of fluidised flow. A correlation has been devised, based on analogy with laminar liquid flow, which describes the channel flow behaviour with reasonable accuracy over the whole range of shear-rates used. 1he channeI results indicated that at low fluidiising velocities the flow was adversely affected by settlement of a stagnant layer of particles on to the distributor, which gave rise to increased flow resistance. Conversely, at higher fluidising velocities the resistance at the distributor appeared to be less than at the walls. In view of this, and also because of the disparity in shear-rates between the two types of apparatus, it is not possible as yet to predict exactly the flow behaviour in an open channel from small-scale viscometer tests.

Modeling ion channel dynamics through reflected stochastic differential equations

Ion channels are membrane proteins that open and close at random and play a vital role in the electrical dynamics of excitable cells. The stochastic nature of the conformational changes these proteins undergo can be significant, however current stochastic modeling methodologies limit the ability to study such systems. Discrete-state Markov chain models are seen as the "gold standard," but are computationally intensive, restricting investigation of stochastic effects to the single-cell level. Continuous stochastic methods that use stochastic differential equations (SDEs) to model the system are more efficient but can lead to simulations that have no biological meaning. In this paper we show that modeling the behavior of ion channel dynamics by a reflected SDE ensures biologically realistic simulations, and we argue that this model follows from the continuous approximation of the discrete-state Markov chain model. Open channel and action potential statistics from simulations of ion channel dynamics using the reflected SDE are compared with those of a discrete-state Markov chain method. Results show that the reflected SDE simulations are in good agreement with the discrete-state approach. The reflected SDE model therefore provides a computationally efficient method to simulate ion channel dynamics while preserving the distributional properties of the discrete-state Markov chain model and also ensuring biologically realistic solutions. This framework could easily be extended to other biochemical reaction networks. © 2012 American Physical Society.

Effect of wave parameters on flood wave subsidence

A detailed analy~is on the propagation of a sinusoidal flood wave in a wide prismatic open channel b.as hen made by numc? ii.~ll~integrating We govemins nondimenional equations of unsteady flow in an open chamei. EmpE:dsis has been laid on the effect of wave parmefen on th propagation of 6.8 sinusoidal wave. Results show that the amount of subsidence is more in the case of small wave anplltude and wave duration cases. Further, wave duration has been noticed to have a relatively Vier influence on subsidence than wave amplitude. The speed at which the peak of the wave moves is observed to be a function of only the wave amplitude.

Metamodeling approach to predict friction factor of alluvial channel

The flow resistance of an alluvial channel flow is not only affected by the Reynolds number and the roughness conditions but also the Froude number. Froude number is the most basic parameter in the case of the alluvial channel, thus effect of Froude number on resistance to flow should be considered in the formulation of the friction factor, which is not in the case of present available resistance equations. At present, no generally acceptable quantitative description of the effects of the Froude number on hydraulic resistance has been developed. Metamodeling technique, which is particularly useful in modeling a complex processes or where knowledge of the physics is limited, is presented as a tool complimentary to modeling friction factor in alluvial channels. Present work uses, a radial basis metamodel, which is a type of neural network modeling, to find the effect of Froude number on the flow resistance. Based on the experimental data taken from different sources, it has been found that the predicting capability of the present model is on acceptable level. Present work also tries in formulating an empirical equation for resistance in alluvial channel comprising all the three majorm, parameters, namely, roughness parameter, Froude number and Reynolds number. (C) 2009 Elsevier B.V. All rights reserved.

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.

Burst Detection in Turbulent Channel Flows Based on Large Eddy Simulation Databases

Reliable turbulent channel flow databases at several Reynolds numbers have been established by large eddy simulation (LES), with two of them validated by comparing with typical direct numerical simulation (DNS) results. Furthermore, the statistics, such as velocity profile, turbulent intensities and shear stress, were obtained as well as the temporal and spatial structure of turbulent bursts. Based on the LES databases available, the conditional sampling methods are used to detect the structures of burst events. A method to deterimine the grouping parameter from the probability distribution function (pdf) curve of the time separation between ejection events is proposed to avoid the errors in detected results. And thus, the dependence of average burst period on thresholds is considerably weakened. Meanwhile, the average burst-to-bed area ratios are detected. It is found that the Reynolds number exhibits little effect on the burst period and burst-to-bed area ratio.

Examination of the LBM in simulation of microchannel flow in transitional regime

The gas flows in micro-electro-mechanical systems possess relatively large Knudsen number and usually belong to the slip flow and transitional flow regimes. Recently the lattice Boltzmann method (LBM) was proposed by Nie et al. in Journal of Statistical Physics, vol. 107, pp. 279-289, in 2002 to simulate the microchannel and microcavity flows in the transitional flow regime. The present article intends to test the feasibility of doing so. The results of using the lattice Boltzmann method and the direct simulation Monte Carlo method show good agreement between them for small Kn (Kn = 0.0194), poor agreement for Kn = 0.194, and large deviation for Kn = 0.388 in simulating microchannel flows. This suggests that the present version of the lattice Boltzmann method is not feasible to simulate the transitional channel flow.

A non-equilibrium sediment transport model for rill erosion

Sediment transport in rill flows exhibits the characteristics of non-equilibrium transport, and the sediment transport rate of rill flow gradually recovers along the flow direction by erosion. By employing the concept of partial equilibrium sediment transport from open channel hydraulics, a dynamic model of rill erosion on hillslopes was developed. In the model, a parameter, called the restoration coefficient of sediment transport capacity, was used to express the recovery process of sediment transport rate, which was analysed by dimensional analysis and determined from laboratory experimental data. The values of soil loss simulated by the model were in agreement with observed values. The model results showed that the length and gradient of the hillslope and rainfall intensity had different influences on rill erosion. Copyright (c) 2006 John Wiley & Sons, Ltd.

Multiscale Method for Turbulent-Flow Computations

A main method of predicting turbulent flows is to solve LES equations, which was called traditional LES method. The traditional LES method solves the motions of large eddies of size larger than filtering scale An while modeling unresolved scales less than Delta_n. Hughes et al argued that many shortcomings of the traditional LES approaches were associated with their inabilities to successfully differentiate between large and small scales. One may guess that a priori scale-separation would be better, because it can predict scale-interaction well compared with posteriori scale-separation. To this end, a multi-scale method was suggested to perform scale-separation computation. The primary contents of the multiscale method are l) A space average is used to differentiate scale. 2) The basic equations include the large scale equations and fluctuation equations. 3) The large-scale equations and fluctuation equations are coupled through turbulent stress terms. We use the multiscale equations of n=2, i.e., the large and small scale (LSS) equations, to simulate 3-D evolutions of a channel flow and a planar mixing layer flow Some interesting results are given.