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.

Polyethylene flow prediction with a differential multi-mode Pom-Pom model

We report the first steps of a collaborative project between the University of Queensland, Polyflow, Michelin, SK Chemicals, and RMIT University; on simulation, validation and application of a recently introduced constitutive model designed to describe branched polymers. Whereas much progress has been made on predicting the complex flow behaviour of many - in particular linear - polymers, it sometimes appears difficult to predict simultaneously shear thinning and extensional strain hardening behaviour using traditional constitutive models. Recently a new viscoelastic model based on molecular topology, was proposed by McLeish and Larson (1998). We explore the predictive power of a differential multi-mode version of the pom-pom model for the flow behaviour of two commercial polymer melts: a (long-chain branched) low-density polyethylene (LDPE) and a (linear) high-density polyethylene (HDPE). The model responses are compared to elongational recovery experiments published by Langouche and Debbaut (1999), and start-up of simple shear flow, stress relaxation after simple and reverse step strain experiments carried out in our laboratory.

Comparison of DEM and experiment for a scale model SAG mill

Predictions of flow patterns in a 600-mm scale model SAG mill made using four classes of discrete element method (DEM) models are compared to experimental photographs. The accuracy of the various models is assessed using quantitative data on shoulder, toe and vortex center positions taken from ensembles of both experimental and simulation results. These detailed comparisons reveal the strengths and weaknesses of the various models for simulating mills and allow the effect of different modelling assumptions to be quantitatively evaluated. In particular, very close agreement is demonstrated between the full 3D model (including the end wall effects) and the experiments. It is also demonstrated that the traditional two-dimensional circular particle DEM model under-predicts the shoulder, toe and vortex center positions and the power draw by around 10 degrees. The effect of particle shape and the dimensionality of the model are also assessed, with particle shape predominantly affecting the shoulder position while the dimensionality of the model affects mainly the toe position. Crown Copyright (C) 2003 Published by Elsevier Science B.V. All rights reserved.

A unified model of flames as gasdynamic discontinuities

Viewed on a hydrodynamic scale, flames in experiments are often thin so that they may be described as gasdynamic discontinuities separating the dense cold fresh mixture from the light hot burned products. The original model of a flame as a gasdynamic discontinuity was due to Darrieus and to Landau. In addition to the fluid dynamical equations, the model consists of a flame speed relation describing the evolution of the discontinuity surface, and jump conditions across the surface which relate the fluid variables on the two sides of the surface. The Darrieus-Landau model predicts, in contrast to observations, that a uniformly propagating planar flame is absolutely unstable and that the strength of the instability grows with increasing perturbation wavenumber so that there is no high-wavenumber cutoff of the instability. The model was modified by Markstein to exhibit a high-wavenumber cutoff if a phenomenological constant in the model has an appropriate sign. Both models are postulated, rather than derived from first principles, and both ignore the flame structure, which depends on chemical kinetics and transport processes within the flame. At present, there are two models which have been derived, rather than postulated, and which are valid in two non-overlapping regions of parameter space. Sivashinsky derived a generalization of the Darrieus-Landau model which is valid for Lewis numbers (ratio of thermal diffusivity to mass diffusivity of the deficient reaction component) bounded away from unity. Matalon & Matkowsky derived a model valid for Lewis numbers close to unity. Each model has its own advantages and disadvantages. Under appropriate conditions the Matalon-Matkowsky model exhibits a high-wavenumber cutoff of the Darrieus-Landau instability. However, since the Lewis numbers considered lie too close to unity, the Matalon-Matkowsky model does not capture the pulsating instability. The Sivashinsky model does capture the pulsating instability, but does not exhibit its high-wavenumber cutoff. In this paper, we derive a model consisting of a new flame speed relation and new jump conditions, which is valid for arbitrary Lewis numbers. It captures the pulsating instability and exhibits the high-wavenumber cutoff of all instabilities. The flame speed relation includes the effect of short wavelengths, not previously considered, which leads to stabilizing transverse surface diffusion terms.

Stability of planar flames as gasdynamic discontinuities

The stability of a steadily propagating planar premixed flame has been the subject of numerous studies since Darrieus and Landau showed that in their model flames are unstable to perturbations of any wavelength. Moreover, the instability was shown to persist even for very small wavelengths, i.e. there was no high-wavenumber cutoff of the instability. In addition to the Darrieus-Landau instability, which results from thermal expansion, analysis of the diffusional thermal model indicates that premixed flames may exhibit cellular and pulsating instabilities as a consequence of preferential diffusion. However, no previous theory captured all the instabilities including a high-wavenumber cutoff for each. In Class, Matkowsky & Klimenko (2003) a unified theory is proposed which, in appropriate limits and under appropriate assumptions, recovers all the relevant previous theories. It also includes additional new terms, not present in previous theories. In the present paper we consider the stability of a uniformly propagating planar flame as a solution of the unified model. The results are then compared to those based on the models of Darrieus-Landau, Sivashinsky and Matalon-Matkowsky. In particular, it is shown that the unified model is the only model to capture the Darrieus-Landau, cellular and pulsating instabilities including a high-wavenumber cutoff for each.

Numerical investigation of turbulent flow around a rotating stepped cylinder for corrosion study

Direct numerical simulation has been carried out for turbulent flow set up by a rotating cylinder with two backward-facing steps axisymmetrically mounted in the circumferential direction. This flow geometry creates a qualitatively similar flow pattern as observed near, a sudden, pipe expansion or a plane backward-facing step, characterized by flow separation and reattachment. A region of intense turbulence intensity and high wall-shear-stress fluctuations is formed in,the recirculating I region downstream of the step, where high mass-transfer capacity was also experimentally observed. Since, corrosion is frequently mass-transfer., controlled, our findings, put forward this apparatus as useful tool for future corrosion research.

Effect of cold wire inclination on temperature fluctuation measurements in a heated jet

Hot-wire anemometers at low operating currents are used as fast response resistance thermometers for the study of heated turbulent flows. Simultaneous measurement of temperature and velocity is generally performed with multi-wire arrays. In order to give good spatial resolution a new layout has been tested which uses an inclined temperature wire positioned parallel to the nearest inclined velocity wire. This leads to an asymmetric wire arrangement relative to the mean flow direction. As expected, a reduction in thermal interference from the velocity wires results when compared with an array containing a temperature wire placed normal to the flow. However, measurement of higher order moments of fluctuating quantities in an axisymmetric jet shows considerable distortion of radial distributions which is traced to alteration of the temperature field sensed by the temperature wire. When inclined velocity sensitive wires contain a temperature component, the latter may be affected by the same phenomenon.

Stratigraphy, facies architecture and tectonic implications of the Upper Devonian to Lower Carboniferous Campwyn Volcanics of the northern New England Fold Belt

Upper Devonian to Lower Carboniferous strata of the Campwyn Volcanics of east central Queensland preserve a substantial sequence of first-cycle volcaniclastic sedimentary and coeval volcanic rocks that record prolonged volcanic activity along the northern New England Fold Belt. The style and scale of volcanism varied with time, producing an Upper Devonian sequence of mafic volcano-sedimentary rocks overlain by a rhyolitic ignimbrite-dominated sequence that passes upward into a Lower Carboniferous limestone-bearing sedimentary sequence. We define two facies associations for the Campwyn Volcanics. A lower facies association is dominated by mafic volcanic-derived sedimentary breccias with subordinate primary mafic volcanic rocks comprising predominantly hyaloclastite and peperite. Sedimentary breccias record episodic and high energy, subaqueous depositional events with clastic material sourced from a mafic lava-dominated terrain. Some breccias contain a high proportion of attenuated dense, glassy mafic juvenile clasts, suggesting a syn-eruptive origin. The lower facies association coarsens upwards from a lithic sand-dominated sequence through a thick interval of pebble- to boulder-grade polymict volcaniclastic breccias, culminating in facies that demonstrate subaerial exposure. The silicic upper facies association marks a significant change in eruptive style, magma composition and the nature of eruptive sources, as well as the widespread development of subaerial depositional conditions. Crystal-rich, high-grade, low- to high-silica rhyolite ignimbrites dominate the base of this facies association. Biostratigraphic age controls indicate that the ignimbrite-bearing sequences are Famennian to lower-mid Tournaisian in age. The ignimbrites represent extra-caldera facies with individual units up to 40 m thick and mostly lacking coarse lithic breccias. Thick deposits of pyroclastic material interbedded with fine-grained siliceous sandstone and mudstone (locally radiolarian-bearing) were deposited from pyroclastic flows that crossed palaeoshorelines or represent syn-eruptive, resedimented pyroclastic material. Some block-bearing lithic-pumice-crystal breccias may also reflect more proximal subaqueous silicic explosive eruptions. Crystal-lithic sandstones interbedded with, and overlying the ignimbrites, contain abundant detrital volcanic quartz and feldspar derived from the pyroclastic deposits. Limestone is common in the upper part of the upper facies association, and several beds are oolitic (cf. Rockhampton Group of the Yarrol terrane). Overall, the upper facies association fines upward and is transgressive, recording a return to shallow-marine conditions. Palaeocurrent data from all stratigraphic levels in the Campwyn Volcanics indicate that the regional sediment-dispersal direction was to the northwest, and opposed to the generally accepted notion of easterly sediment dispersal from a volcanic arc source. The silicic upper facies association correlates in age and lithology to Early Carboniferous silicic volcanism in the Drummond (Cycle 1) and Burdekin Basins, Connors Arch, and in the Yarrol terranes of eastern Queensland. The widespread development of silicic volcanism in the Early Carboniferous indicates that silicic (rift-related) magmatism was not restricted to the Drummond Basin, but was part of a more substantial silicic igneous province.

Numerical simulations of wind and temperature structure within an Alpine lake basin, Lake Tekapo, New Zealand

High-resolution numerical model simulations have been used to study the local and mesoscale thermal circulations in an Alpine lake basin. The lake (87 km(2)) is situated in the Southern Alps, New Zealand and is located in a glacially excavated rock basin surrounded by mountain ranges that reach 3000 m in height. The mesoscale model used (RAMS) is a three-dimensional non-hydrostatic model with a level 2.5 turbulence closure scheme. The model demonstrates that thermal forcing at local (within the basin) and regional (coast-to-basin inflow) scales drive the observed boundary-layer airflow in the lake basin during clear anticyclonic summertime conditions. The results show that the lake can modify (perturb) both the local and regional wind systems. Following sunrise, local thermal circulations dominate, including a lake breeze component that becomes embedded within the background valley wind system. This results in a more divergent flow in the basin extending across the lake shoreline. However, a closed lake breeze circulation is neither observed nor modelled. Modelling results indicate that in the latter part of the day when the mesoscale (coast-to-basin) inflow occurs, the relatively cold pool of lake air in the basin can cause the intrusion to decouple from the surface. Measured data provide qualitative and quantitative support for the model results.

Integrated control of nitrate recirculation and external carbon addition in a predenitrification system

The integrated control of nitrate recirculation and external carbon addition in a predenitrification biological wastewater treatment system is studied. The proposed control structure consists of four feedback control loops, which manipulate the nitrate recirculation and the carbon dosage flows in a highly coordinated manner such that the consumption of external carbon is minimised while the nitrate discharge limits (based on both grab and composite samples) are met. The control system requires the measurement of the nitrate concentrations at the end of both the anoxic and the aerobic zones. Distinct from ordinary control systems, which typically minimise the variation in the controlled variables, the proposed control system essentially maximises the diurnal variation of the effluent nitrate concentration and through this maximises the use of influent COD for denitrification, thus minimising the requirement for external carbon source. Simulation studies using a commonly accepted simulation benchmark show that the controlled system consistently achieves the designated effluent quality with minimum costs.

Water table waves in an unconfined aquifer: Experiments and modeling

[1] Comprehensive measurements are presented of the piezometric head in an unconfined aquifer during steady, simple harmonic oscillations driven by a hydrostatic clear water reservoir through a vertical interface. The results are analyzed and used to test existing hydrostatic and nonhydrostatic, small-amplitude theories along with capillary fringe effects. As expected, the amplitude of the water table wave decays exponentially. However, the decay rates and phase lags indicate the influence of both vertical flow and capillary effects. The capillary effects are reconciled with observations of water table oscillations in a sand column with the same sand. The effects of vertical flows and the corresponding nonhydrostatic pressure are reasonably well described by small-amplitude theory for water table waves in finite depth aquifers. That includes the oscillation amplitudes being greater at the bottom than at the top and the phase lead of the bottom compared with the top. The main problems with respect to interpreting the measurements through existing theory relate to the complicated boundary condition at the interface between the driving head reservoir and the aquifer. That is, the small-amplitude, finite depth expansion solution, which matches a hydrostatic boundary condition between the bottom and the mean driving head level, is unrealistic with respect to the pressure variation above this level. Hence it cannot describe the finer details of the multiple mode behavior close to the driving head boundary. The mean water table height initially increases with distance from the forcing boundary but then decreases again, and its asymptotic value is considerably smaller than that previously predicted for finite depth aquifers without capillary effects. Just as the mean water table over-height is smaller than predicted by capillarity-free shallow aquifer models, so is the amplitude of the second harmonic. In fact, there is no indication of extra second harmonics ( in addition to that contained in the driving head) being generated at the interface or in the interior.

On the evaluation of some differential formulations for the pom-pom constitutive model

In recent years, progress has been made in modelling long chain branched polymers by the introduction of the so-called pompom model. Initially developed by McLeish and Larson (1998), the model has undergone several improvements or alterations, leading to the development of new formulations. Some of these formulations however suffer from certain mathematical defects. The purpose of the present paper is to review some of the formulations of the pom-pom constitutive model, and to investigate their possible mathematical defects. Next, an alternative formulation is proposed, which does not appear to exhibit mathematical defects, and we explore its modelling performance by comparing the predictions with experiments in non-trivial rheometric flows of an LDPE melt. The selected rheometric flows are the double step strain, as well as the large amplitude oscillatory shear experiments. For LAOS experiments, the comparison involves the use of Fourier-transform analysis.