Turbulence Characterisation of High-Velocity Free-Surface Flows

In high-velocity free-surface flows, air is continuously being trapped and released through the free-surface. Such high-velocity highly-aerated flows cannot be studied numerically because of the large number of relevant equations and parameters. Herein an advanced signal processing of traditional single- and dual-tip conductivity probes provides some new information on the air-water turbulent time and length scales. The technique is applied to turbulent open channel flows in a large-size facility. The auto- and cross-correlation analyses yield some characterisation of the large eddies advecting the bubbles. The transverse integral turbulent length and time scales are related to the step height: i.e., Lxy/h ~ 0.02 to 0.2, and T.sqrt(g/h) ~ 0.004 to 0.04. The results are irrespective of the Reynolds numbers. The present findings emphasise that turbulent dissipation by large-scale vortices is a significant process in the intermediate zone between the spray and bubbly flow regions (0.3 < C < 0.7). Some self-similar relationships were observed systematically at both macroscopic and microscopic levels. The results are significant because they provide a picture general enough to be used to characterise the air-water flow field in prototype spillways.

Velocimetry and thermometry of supersonic flow around a cylindrical body

The supersonic flow around a cylindrical body has been studied using two optical techniques. For both sets of measurements, the cylinder was mounted from the side of the tunnel, allowing investigation of the bow shock region as well as in the wake. A new technique, laser-enhanced ionization flow tagging, was used for streamwise velocity determinations behind the body. From these measurements, it was found that the downstream velocity outside the wake was (1.90 +/- 0.06) km/s, whereas inside the wake the velocity was about 0-500 m/s in the upstream direction. Planar laser induced fluorescence of nitric oxide was employed for temperature determinations. It was established that the freestream temperature was (2120 +/- 100) K, decreasing to around (1550 +/- 400) K in the wake.

Hypoxia stimulates cerebral blood flow in the estuarine crocodile (Crocodylus porosus)

The effect of N-2 respiration on cerebral blood flow (CBF) velocity on the dorsal surface of cerebellum was examined in the estuarine crocodile, Crocodylus porosus, using epi-illumination microscopy. Twelve minutes of N-2 respiration resulted in a 126% increase in CBF velocity. N-2 respiration had no effect on blood pressure, indicating an underlying cerebral vasodilation. In addition, heart rate increased significantly. Systemic injections of aminophylline and the NO synthase (NOS) inhibitor nitro-L-arginine (L-NA) did not affect the hypoxia induced increase in CBF. We conclude that C. porosus responds to hypoxia with adenosine and nitric oxide (NO) independent cerebral vasodilation, and that this is likely to be a mechanism protecting the brain from energy deficiency during prolonged dives. (C) 1999 Elsevier Science Ireland Ltd. All rights reserved.

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

Model for gas-liquid flow through wet porous medium

A method involving bubbling of air through a fibrous filter immersed in water has recently been investigated (Agranovski et al. [1]). Experimental results showed that the removal efficiency for ultra-fine aerosols by such filters was greatly increased compared to dry filters. Nuclear Magnetic Resonance (NMR) imaging was used to examine the wet filter and to determine the nature of the gas flow inside the filter (Agranovski et al. [2]). It was found that tortuous preferential pathways (or flow tubes) develop within the filter through which the air flows and the distribution of air and water inside the porous medium has been investigated. The aim of this paper is to investigate the geometry of the pathways and to make estimates of the flow velocities and particle removal efficiency in such pathways. A mathematical model of the flow of air along the preferred pathways has been developed and verified experimentally. Even for the highest realistic gas velocity the flow field was essentially laminar (Re approximate to 250). We solved Laplace's equation for stream function to map trajectories of particles and gas molecules to investigate the possibility of their removal from the carrier.

Analytical solution to the zero-inertia problem for surge flow phenomena in nonprismatic channels

Surge flow phenomena. e.g.. as a consequence of a dam failure or a flash flood, represent free boundary problems. ne extending computational domain together with the discontinuities involved renders their numerical solution a cumbersome procedure. This contribution proposes an analytical solution to the problem, It is based on the slightly modified zero-inertia (ZI) differential equations for nonprismatic channels and uses exclusively physical parameters. Employing the concept of a momentum-representative cross section of the moving water body together with a specific relationship for describing the cross sectional geometry leads, after considerable mathematical calculus. to the analytical solution. The hydrodynamic analytical model is free of numerical troubles, easy to run, computationally efficient. and fully satisfies the law of volume conservation. In a first test series, the hydrodynamic analytical ZI model compares very favorably with a full hydrodynamic numerical model in respect to published results of surge flow simulations in different types of prismatic channels. In order to extend these considerations to natural rivers, the accuracy of the analytical model in describing an irregular cross section is investigated and tested successfully. A sensitivity and error analysis reveals the important impact of the hydraulic radius on the velocity of the surge, and this underlines the importance of an adequate description of the topography, The new approach is finally applied to simulate a surge propagating down the irregularly shaped Isar Valley in the Bavarian Alps after a hypothetical dam failure. The straightforward and fully stable computation of the flood hydrograph along the Isar Valley clearly reflects the impact of the strongly varying topographic characteristics on the How phenomenon. Apart from treating surge flow phenomena as a whole, the analytical solution also offers a rigorous alternative to both (a) the approximate Whitham solution, for generating initial values, and (b) the rough volume balance techniques used to model the wave tip in numerical surge flow computations.

Simulation of mixing in unsteady flow through a periodically obstructed channel

Fluid mixing in steady and unsteady Bow through a channel containing periodic square obstructions has been studied using a finite-difference simulation to determine fluid velocities, followed by the use of passive marker particle advection to look at fluid transport out of the cavities formed between each of the obstructions. The geometry and Bow conditions were chosen from the work by Perkins (1989, M.S. Thesis, Lehigh University; 1992, Ph.D. Thesis, Lehigh University); who investigated heat transfer enhancement due to unsteady flow through such an obstructed channel. Particle advection shows that Bow regimes which are predicted to give good mixing based on snapshots of instantaneous streamline contour plots were not necessarily able to efficiently mix fluid which started in the cavity regions throughout the channel. The use of Poincare sections shows regular regions existing under these conditions which inhibit efficient fluid transport. These regular regions are found to disappear when the unsteady Bow velocity is increased. (C) 1997 Elsevier Science Ltd.

Ionic strontium fluorescence as a method for flow tagging velocimetry

A flow tagging technique based upon ionic fluorescence in strontium is investigated for applications to velocity measurements in gas flows. The method is based upon a combination of two laser based spectroscopic techniques, i.e. resonantly-enhanced ionisation and laser-induced ionic fluorescence. Strontium is first ionised and then planar laser-induced fluorescence is utilised to give 2D 'bright images' of the ionised region of the flow at a given time delay. The results show that this method can be used for velocity measurements. The velocities were measured in two types of air-acetylene flames - a slot burner and a circular burner yielding velocities of 5.1 +/- 0.1 m/s and 9.3 +/- 0.2 m/s, respectively. The feasibility of the method for the determination of velocities in faster flows than those investigated here is discussed.

Experimental investigation of liquid flow shift due to gas cross flow in non-wetted packed beds

An experimental study has been carried out for the gas-liquid two-phase flow in a packed bed simulating conditions of the gas and liquid flows in the lower part of blast furnace. The localised liquid flow phenomenon in presence of gas cross flow, which usually occurs around the cohesive zone and raceway in blast furnace, was investigated in detail. Such liquid flow is characterised in terms of liquid shift distance or liquid shift angle that can effectively be measured by the experiments involved in the current study. It is found that liquid shift angle does not significantly increase or decrease with different packing depth. This finding supports the hypothesis of the force balance model where a vectorial relationship among acting forces, i.e. gas drag force, gravitational force and solid-liquid friction force, and liquid shift angle does exist. Liquid shift angle is inversely proportional to particle size and liquid density, and proportional to square of gas superficial velocity, but is almost independent on liquid flowrate and liquid viscosity. The gas-liquid drag coefficient, an important aspect for quantifying the interaction between gas and liquid flows, was conceptually modified based on the discrete feature of liquid flow through a packed bed and evaluated by the combined theoretical and experimental investigation. Experimental measurements suggest that the gas-liquid drag coefficient is approximately a constant (C-DG(')=5.4+/-1.0) and is independent on liquid properties, gas velocity and packing structure. The result shows a good agreement with previous experimental data and prediction of the existing liquid flow model.

Shear stress and sediment transport calculations for sheet flow under waves

A simple method is provided for calculating transport rates of not too fine (d(50) greater than or equal to 0.20 mm) sand under sheet flow conditions. The method consists of a Meyer-Peter-type transport formula operating on a time-varying Shields parameter, which accounts for both acceleration-asymmetry and boundary layer streaming. While velocity moment formulae, e.g.., = Constant x calibrated against U-tube measurements, fail spectacularly under some real waves (Ribberink, J.S., Dohmen-Janssen, C.M., Hanes, D.M., McLean, S.R., Vincent, C., 2000. Near-bed sand transport mechanisms under waves. Proc. 27th Int. Conf. Coastal Engineering, Sydney, ASCE, New York, pp. 3263-3276, Fig. 12), the new method predicts the real wave observations equally well. The reason that the velocity moment formulae fail under these waves is partly the presence of boundary layer streaming and partly the saw-tooth asymmetry, i.e., the front of the waves being steeper than the back. Waves with saw-tooth asymmetry may generate a net landward sediment transport even if = 0, because of the more abrupt acceleration under the steep front. More abrupt accelerations are associated with thinner boundary layers and greater pressure gradients for a given velocity magnitude. The two real wave effects are incorporated in a model of the form Q(s)(t) = Q(s)[theta(t)] rather than Q(S)(t) = Q(S)[u(infinity)(t)], i.e., by expressing the transport rate in terms of an instantaneous Shields parameter rather than in terms of the free stream velocity, and accounting for both streaming and accelerations in the 0(t) calculations. The instantaneous friction velocities u(*)(t) and subsequently theta(t) are calculated as follows. Firstly, a linear filter incorporating the grain roughness friction factor f(2.5) and a phase angle phi(tau) is applied to u(infinity)(t). This delivers u(*)(t) which is used to calculate an instantaneous grain roughness Shields parameter theta(2.5)(t). Secondly, a constant bed shear stress is added which corresponds to the streaming related bed shear stress -rho ($) over bar((u) over tilde(w) over tilde)(infinity) . The method can be applied to any u(infinity)(t) time series, but further experimental validation is recommended before application to conditions that differ strongly from the ones considered below. The method is not recommended for rippled beds or for sheet flow with typical prototype wave periods and d(50) < 0.20 turn. In such scenarios, time lags related to vertical sediment movement become important, and these are not considered by the present model. (C) 2002 Elsevier Science B.V. All rights reserved.

Vertical fluxes of sediment in oscillatory sheet flow

Time series of vertical sediment fluxes are derived from concentration time series in sheet flow under waves. While the concentrations C(z,t) vary very little with time for \z\ < 10d(50), the measured vertical sediment fluxes Q(zs)(z,t) vary strongly with time in this vertical band and their time variation follows, to some extent, the variation of the grain roughness Shields parameter 02,5(t). Thus, sediment distribution models based on the pickup function boundary condition are in some qualitative agreement with the measurements. However, the pickup function models are only able to model the upward bursts of sediment during the accelerating phases of the flow. They are, so far, unable to model the following strong downward sediment fluxes, which are observed during the periods of flow deceleration. Classical pickup functions, which essentially depend on the Shields parameter, are also incapable of modelling the secondary entrainment fluxes, which sometimes occur at free stream velocity reversal. The measured vertical fluxes indicate that the effective sediment settling velocity in the high [(0.3 < C(z,t) < 0.4] concentration area is typically only a few percent of the clear water settling velocity, while the measurements of Richardson and Jeronimo [Chem. Eng. Sci. 34 (1979) 1419], from a different physical setting, lead to estimates of the order 20%. The data does not support gradient diffusion as a model for sediment entrainment from the bed. That is, detailed modelling of the observed near-bed fluxes would require diffusivities that go negative during periods of flow deceleration. An observed general trend for concentration variability to increase with elevation close to the bed is also irreconcilable with diffusion models driven by a bottom boundary condition. (C) 2002 Published by Elsevier Science B.V.

Investigation of the flow field of a highly heated jet of air

Measurements of mean and fluctuating velocity and temperature and their self- and cross-products to the third-order are presented for a heated axisymmetric air jet. Froude numbers in the range of 3500 13,190, Reynolds numbers in the range of 3470-8500 and non-dimensional streamwise distances. X*, from 0.27 to 1.98 are covered by the data. It was found that turbulence intensity decreases for the heated jet in the region between the inertia dominated and the buoyancy dominated regions which is contrary to findings with helium jets mixing with ambient air to produce density fluctuations. The effects of heating on the turbulent kinetic energy budget and the temperature variance budget show small differences for the inertia dominated region and the intermediate region which help to explain the transition process to the far field plume region. Constants are evaluated for the isotropic eddy diffusivity and generalised gradient hypothesis models as well as the scalar variance model. No significant effect of heating on the dissipation time-scale ratio was found. A novel wire array with an inclined cold wire was used. Measurements obtained with this probe are found to lead to asymmetries in some of the higher-order products. Further investigation suggested that the asymmetries are attributable to an as yet unreported interference effect produced by the leading prong of the inclined temperature wire, The effect may also have implications for inclined velocity wires which contain a temperature component when used in heated flows. (C) 2002 Elsevier Science Inc. All rights reserved.

Air-water flows down stepped chutes: turbulence and flow structure observations

Interactions between turbulent waters and atmosphere may lead to strong air-water mixing. This experimental study is focused on the flow down a staircase channel characterised by very strong flow aeration and turbulence. Interfacial aeration is characterised by strong air-water mixing extending down to the invert. The size of entrained bubbles and droplets extends over several orders of magnitude, and a significant number of bubble/droplet clusters was observed. Velocity and turbulence intensity measurements suggest high levels of turbulence across the entire air-water flow. The increase in turbulence levels, compared to single-phase flow situations, is proportional to the number of entrained particles. (C) 2002 Elsevier Science Ltd. All rights reserved.

Unsaturated liquid percolation flow through nonwetted packed beds

The unsaturated flow of liquid through packed beds of large particles was studied using six different liquids, all with contact angles greater than 90degrees on the bed packing (wax spheres of 9, 15 and 19.4 mm diameter). The liquid flow was discrete in nature, as drops for low flow rates and rivulets for high flow rates. For unsaturated liquid flows, the actual percolation velocity, not superficial velocity, should be used to characterize the flow. The percolation velocity did not vary with packed-bed depth, but was a strong function of liquid flow rate, liquid and particle properties. Effects of liquid and particle properties (but not flow rate) are well captured by a simple correlation between the liquid-particle friction factor and Reynolds number based on actual percolation velocities. Liquid dispersion, characterized by the maximum dispersion angle, varies significantly with liquid and particle properties. The tentative correlation suggested here needs further validation for a wider range of conditions.

The interaction of water flow and nutrients on aquatic plant growth

A long-term experiment was conducted to compare the effects of flowing and still water on growth, and the relationship between water flow and nutrients, in Aponogeton elongatus, a submerged aquatic macrophyte. A. elongatus plants were grown for 23 weeks with three levels of nutrition (0, 0.5 and 1g Osmocote Plus(R) fertiliser pot(-1)) in aquaria containing stirred or unstirred water. Fertilized plants grew much better than non-fertilized. The highest fertilizer level produced 29% wider leaves and 58% higher total dry weight in stirred water. Stirred water increased leaf area by 40% and tuber size by 81%, but only with the highest level of nutrition. These results suggest that this plant depends on its roots for mineral uptake, rather than from the open water, and the major limitation to growth in still water is the supply of dissolved inorganic carbon. It was the combined effects of nutrient availability and stirring that produced the strongest response in plant growth, morphology and composition. This study provides some explanation for the observations of others that these plants grow best in creeks or river systems with permanently flowing water.