Bubbly Flow Structure in Hydraulic Jump

In an open channel, a hydraulic jump is the rapid transition from super- to sub-critical flow associated with strong turbulence and air bubble entrainment in the mixing layer. New experiments were performed at relatively large Reynolds numbers using phase-detection probes. Some new signal analysis provided characteristic air-water time and length scales of the vortical structures advecting the air bubbles in the developing shear flow. An analysis of the longitudinal air-water flow structure suggested little bubble clustering in the mixing layer, although an interparticle arrival time analysis showed some preferential bubble clustering for small bubbles with chord times below 3 ms. Correlation analyses yielded longitudinal air-water time scales Txx*V1/d1 of about 0.8 in average. The transverse integral length scale Z/d1 of the eddies advecting entrained bubbles was typically between 0.25 and 0.4, irrespective of the inflow conditions within the range of the investigations. Overall the findings highlighted the complicated nature of the air-water flow

Analysis of stirred mill performance using DEM simulation: Part 2 - Coherent flow structures, liner stress and wear, mixing and transport

Stirred Mills are becoming increasingly used for fine and ultra-fine grinding. This technology is still poorly understood when used in the mineral processing context. This makes process optimisation of such devices problematic. 3D DEM simulations of the flow of grinding media in pilot scale tower mills and pin mills are carried out in order to investigate the relative performance of these stirred mills. In the first part of this paper, media flow patterns and energy absorption rates and distributions were analysed to provide a good understanding of the media flow and the collisional environment in these mills. In this second part we analyse steady state coherent flow structures, liner stress and wear by impact and abrasion. We also examine mixing and transport efficiency. Together these provide a comprehensive understanding of all the key processes operating in these mills and a clear understanding of the relative performance issues. (C) 2006 Elsevier Ltd. All rights reserved.

Turbulent flow around a rotating stepped cylinder

Direct numerical simulation (DNS) of turbulent flow around a rotating cylinder with two backward-facing steps axisymmetrically mounted in the circumferential direction was performed and compared with DNS of plane backward-facing step flow (PBSF) of Le [J. Fluid Mech. 330, 349 (1997)]. The original motivation of this work stemmed from the efforts to design a simple device which can generate flows of high turbulence intensity at low cost for corrosion researchers. It turned out that the current flow shows flow structures quite similar to those of PBSF downstream of the step, even though configurations of the two flows are totally different from one another. The stepped cylinder appears to be a cost-effective tool in the generation of flow structures similar to those of PBSF. (C) 2002 American Institute of Physics.

Air Entrainment Processes in a Full-Scale Rectangular Dropshaft at Large Flows

Rectangular dropshafts, commonly used in sewers and storm water systems, are characterised by significant flow aeration. New detailed air-water flow measurements were conducted in a near-full-scale dropshaft at large discharges. In the shaft pool and outflow channel, the results demonstrated the complexity of different competitive air entrainment mechanisms. Bubble size measurements showed a broad range of entrained bubble sizes. Analysis of streamwise distributions of bubbles suggested further some clustering process in the bubbly flow although, in the outflow channel, bubble chords were in average smaller than in the shaft pool. A robust hydrophone was tested to measure bubble acoustic spectra and to assess its field application potential. The acoustic results characterised accurately the order of magnitude of entrained bubble sizes, but the transformation from acoustic frequencies to bubble radii did not predict correctly the probability distribution functions of bubble sizes.

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.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Analysis of stirred mill performance using DEM simulation: Part 1 - Media motion, energy consumption and collisional environment

Stirred mills are becoming increasingly used for fine and ultra-fine grinding. This technology is still poorly understood when used in the mineral processing context. This makes process optimisation of such devices problematic. 3D DEM simulations of the flow of grinding media in pilot scale tower mills and pin mills are carried out in order to investigate the relative performance of these stirred mills. Media flow patterns and energy absorption rates and distributions are analysed here. In the second part of this paper, coherent flow structures, equipment wear and mixing and transport efficiency are analysed. (C) 2006 Published by Elsevier Ltd.

Turbulence in Hydraulic Jumps: Experimental Measurements

A hydraulic jump is the transition from a supercritical open channel flow to a subcritical regime. It is characterised by a highly turbulent flow with macro-scale vortices, some kinetic energy dissipation and a bubbly two-phase flow structure. New air-water flow measurements were performed in hydraulic jump flows for a range of inflow Froude numbers. The experiments were conducted in a large-size facility using two types of phase-detection intrusive probes: i.e., single-tip and double-tip conductivity probes. These were complemented by some measurements of free-surface fluctuations using ultrasonic displacement meters. The present study was focused on the turbulence characteristics of hydraulic jumps with partially-developed inflow conditions. The void fraction measurements showed the presence of an advective diffusion shear layer in which the void fractions profiles matched closely an analytical solution of the advective diffusion equation for air bubbles. The present results highlighted some influence of the inflow Froude number onto the air bubble entrainment process. At the largest Froude numbers, the advected air bubbles were more thoroughly dispersed vertically, and larger amount of air bubbles were detected in the turbulent shear layer. In the air-water mixing layer, the maximum void fraction and bubble count rate data showed some longitudinal decay function in the flow direction. Such trends were previously reported in the literature. The measurements of interfacial velocity and turbulence level distributions provided new information on the turbulent velocity field in the highly-aerated shear region. The present data suggested some longitudinal decay of the turbulence intensity. The velocity profiles tended to follow a wall jet flow pattern. The air–water turbulent time and length scales were deduced from some auto- and cross-correlation analyses based upon the method of CHANSON (2006,2007). The results provided the integral turbulent time and length scales of the eddy structures advecting the air bubbles in the developing shear layer. The experimental data showed that the auto-correlation time scale Txx was larger than the transverse cross-correlation time scale Txz. The integral turbulence length scale Lxz was a function of the inflow conditions, of the streamwise position (x-x1)/d1 and vertical elevation y/d1. Herein the dimensionless integral turbulent length scale Lxz/d1 was closely related to the inflow depth: i.e., Lxz/d1 = 0.2 to 0.8, with Lxz increasing towards the free-surface. The free-surface fluctuations measurements showed large turbulent fluctuations that reflected the dynamic, unsteady structure of the hydraulic jumps. A linear relationship was found between the normalized maximum free-surface fluctuation and the inflow Froude number.

Finite element modelling of dissipative structures for nonequilibrium chemical reactions in fluid-saturated porous media

We use the finite element method to model and predict the dissipative structures of chemical species for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. In particular, we explore the conditions under which dissipative structures of the species may exist in the Brusselator type of nonequilibrium chemical reaction. Since this is the first time the finite element method and related strategies have been used to study the chemical instability problems in a fluid-saturated porous medium, it is essential to validate the method and strategies before they are put into application. For this purpose, we have rigorously derived the analytical solutions for dissipative structures of chemical species in a benchmark problem, which geometrically is a square. Comparison of the numerical solutions with the analytical ones demonstrates that the proposed numerical method and strategy are robust enough to solve chemical instability problems in a fluid-saturated porous medium. Finally, the related numerical results from two application examples indicate that both the regime and the magnitude of pore-fluid flow have significant effects on the nature of the dissipative structures that developed for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. The motivation for this study is that self-organization under conditions of pore-fluid flow in a porous medium is a potential mechanism of the orebody formation and mineralization in the upper crust of the Earth. (C) 2000 Elsevier Science S.A. All rights reserved.

Hydraulic tracer studies in a pilot scale subsurface flow constructed wetland

Current design procedures for Subsurface Flow (SSF) Wetlands are based on the simplifying assumptions of plug flow and first order decay of pollutants. These design procedures do yield functional wetlands but result in over-design and inadequate descriptions of the pollutant removal mechanisms which occur within them. Even though these deficiencies are often noted, few authors have attempted to improve modelling of either flow or pollutant removal in such systems. Consequently the Oxley Creek Wetland, a pilot scale SSF wetland designed to enable rigorous monitoring, has recently been constructed in Brisbane, Australia. Tracer studies have been carried out in order to determine the hydraulics of this wetland prior to commissioning it with sealed sewage. The tracer studies will continue during the wetland's commissioning and operational phases. These studies will improve our understanding of the hydraulics of newly built SSF wetlands and the changes brought on by operational factors such as biological films and wetland plant root structures. Results to date indicate that the flow through the gravel beds is not uniform and cannot be adequately modelled by a single parameter, plug flow with dispersion, model. We have developed a multiparameter model, incorporating four plug flow reactors, which provides a better approximation of our experimental data. With further development this model will allow improvements to current SSF wetland design procedures and operational strategies, and will underpin investigations into the pollutant removal mechanisms at the Oxley Creek Wetland. (C) 1997 IAWQ. Published by Elsevier Science Ltd.

Knowledge-based system on optimum design of liquid retaining structures with genetic algorithms

This paper delineates the development of a prototype hybrid knowledge-based system for the optimum design of liquid retaining structures by coupling the blackboard architecture, an expert system shell VISUAL RULE STUDIO and genetic algorithm (GA). Through custom-built interactive graphical user interfaces under a user-friendly environment, the user is directed throughout the design process, which includes preliminary design, load specification, model generation, finite element analysis, code compliance checking, and member sizing optimization. For structural optimization, GA is applied to the minimum cost design of structural systems with discrete reinforced concrete sections. The design of a typical example of the liquid retaining structure is illustrated. The results demonstrate extraordinarily converging speed as near-optimal solutions are acquired after merely exploration of a small portion of the search space. This system can act as a consultant to assist novice designers in the design of liquid retaining structures.

A coupled knowledge-based expert system for design of liquid-retaining structures

This paper describes a coupled knowledge-based system (KBS) for the design of liquid-retaining structures, which can handle both the symbolic knowledge processing based on engineering heuristics in the preliminary synthesis stage and the extensive numerical crunching involved in the detailed analysis stage. The prototype system is developed by employing blackboard architecture and a commercial shell VISUAL RULE STUDIO. Its present scope covers design of three types of liquid-retaining structures, namely, a rectangular shape with one compartment, a rectangular shape with two compartments and a circular shape. Through custom-built interactive graphical user interfaces, the user is directed throughout the design process, which includes preliminary design, load specification, model generation, finite element analysis, code compliance checking and member sizing optimization. It is also integrated with various relational databases that provide the system with sectional properties, moment and shear coefficients and final member details. This system can act as a consultant to assist novice designers in the design of liquid-retaining structures with increase in efficiency and optimization of design output and automated record keeping. The design of a typical example of the liquid-retaining structure is also illustrated. (C) 2003 Elsevier B.V All rights reserved.

An expert system on design of liquid-retaining structures with blackboard architecture

The design of liquid-retaining structures involves many decisions to be made by the designer based on rules of thumb, heuristics, judgement, codes of practice and previous experience. Structural design problems are often ill structured and there is a need to develop programming environments that can incorporate engineering judgement along with algorithmic tools. Recent developments in artificial intelligence have made it possible to develop an expert system that can provide expert advice to the user in the selection of design criteria and design parameters. This paper introduces the development of an expert system in the design of liquid-retaining structures using blackboard architecture. An expert system shell, Visual Rule Studio, is employed to facilitate the development of this prototype system. It is a coupled system combining symbolic processing with traditional numerical processing. The expert system developed is based on British Standards Code of Practice BS8007. Explanations are made to assist inexperienced designers or civil engineering students to learn how to design liquid-retaining structures effectively and sustainably in their design practices. The use of this expert system in disseminating heuristic knowledge and experience to practitioners and engineering students is discussed.

Timing and chemistry of fluid-flow events in the Lawn Hill Platform, Northern Australia

Mudrocks and carbonates of the Isa superbasin in the Lawn Hill platform in northern Australia host major base metal sulfide mineralization, including the giant strata-bound Century Zn-Pb deposit. Mineral paragenesis, stable isotope, and K-Ar dating studies demonstrate that long-lived structures such as the Termite Range fault acted as hot fluid conduits several times during the Paleoproterozoic and Mesoproterozoic in response to major tectonic events. Illite and chlorite crystallinity studies suggest the southern part of the platform has experienced higher temperatures (up to 300 degrees C) than similar stratigraphic horizons in the north. The irregular downhole variation of illite crystallinity values provides further information oil the thermal regime in the basin and shows that clay formation was controlled not only by temperature increase with depth but also by high water/rock ratios along relatively permeable zones. K-Ar dating of illite, in combination with other data, may indicate three major thermal events in the central and northern Lawn Hill platform Lit 1500, 1440 to 1400, and 1250 to 1150 Ma. This study did not detect the earlier Century base metal mineralizing event at 1575 Ma. 1500 Ma ages are recorded only in the south and correspond to the age of the Late Isan orogeny and deposition of the Lower Roper superbasin. They may reflect exhumation of a provenance region. The 1440 to 1300 Ma ages are related to fault reactivation and a thermal pulse at similar to 1440 to 1400 Ma possibly accompanied by fluid flow, with subsequent enhanced cooling possibly due to thermal relaxation or further crustal exhumation. The youngest thermal and/or fluid-flow event at 1250 to 1150 Ma is recorded mainly to the cast of the Tern-lite Range fault and may be related to the assembly of the Rodinian supercontinent. Fluids in equilibrium with illite that formed over a range of temperatures, at different times in different parts of the platform. have relatively uniform oxygen isotope compositions and more variable hydrogen isotope compositions (delta O-18 = 3.5-9.7 parts per thousand V-SMOW; delta D = -94 to -36 parts per thousand V-SMOW). The extent of the 180 enrichment and the variably depleted hydrogen isotope compositions suggest the illite interacted with deep-basin hypersaline brines that were composed of evaporated seawater and/or highly evolved meteoric water. Siderite is the most abundant iron-rich gangue phase in the Century Zn-Pb deposit, which is surrounded by all extensive ferroan carbonate alteration halo. Modeling suggests that the ore siderite formed at temperatures of 120 degrees to 150 degrees C, whereas siderite and ankerite in the alteration halo formed at temperatures of 150 degrees to 180 degrees C. The calculated isotopic compositions of the fluids are consistent with O-18-rich basinal brines and mixed inorganic and organic carbon Sources (6180 = 3-10 parts per thousand V-SMOW, delta C-13 = -7 to -3 parts per thousand V-PDB). in the northeast Lawn Hill platform carbonate-rich rocks preserve marine to early diagenetic carbon and oxygen isotope compositions, whereas ferroan carbonate cements in siltstones and shales in the Desert Creek borehole are O-18 and C-13 depleted relative to the sedimentary carbonates. The good agreement between temperature estimates from illite crystallinity and organic reflectance (160 degrees-270 degrees C) and inverse correlation with carbonate delta O-18 values indicates that organic maturation and carbonate precipitation in the northeast Lawn Hill platform resulted from interaction with the 1250 to 1150 Ma fluids. The calculated isotopic compositions of the fluid are consistent with evolved basinal brine (delta O-18 = 5.1-9.4 parts per thousand V-SMOW; delta C-13 = -13.2 to -3.7 parts per thousand V-PDB) that contained a variable organic carbon component from the oxidation and/or hydrolysis of organic matter in the host sequence. The occurrence of extensive O-18- and C-13-depleted ankerite and siderite alteration in Desert Creek is related to the high temperature of the 1250 to 1150 Ma fluid-flow event in the northeast Lawn Hill platform, in contrast to the lower temperature fluids associated with the earlier Century Zn-Pb deposit in the central Lawn Hill platform.