Different mineralization styles in a volcanic-hosted ore deposit: the fluid and isotopic signatures of the Mt Morgan Au-Cu deposit, Australia

Quantitative laser ablation (LA)-ICP-MS analyses of fluid inclusions, trace element chemistry of sulfides, stable isotope (S), and Pb isotopes have been used to discriminate the formation of two contrasting mineralization styles and to evaluate the origin of the Cu and Au at Mt Morgan. The Mt Morgan Au-Cu deposit is hosted by Devonian felsic volcanic rocks that have been intruded by multiple phases of the Mt Morgan Tonalite, a low-K, low-Al2O3 tonalite-trondhjemite-dacite (TTD) complex. An early, barren massive sulfide mineralization with stringer veins is conforming to VHMS sub-seafloor replacement processes, whereas the high-grade Au-Cu. ore is associated with a later quartz-chalcopyrite-pyrite stock work mineralization that is related to intrusive phases of the Tonalite complex. LA-ICP-MS fluid inclusion analyses reveal high As (avg. 8850 ppm) and Sb (avg. 140 ppm) for the Au-Cu mineralization and 5 to 10 times higher Cu concentration than in the fluids associated with the massive pyrite mineralization. Overall, the hydrothermal system of Mt Morgan is characterized by low average fluid salinities in both mineralization styles (45-80% seawater salinity) and temperatures of 210 to 270 degreesC estimated from fluid inclusions. Laser Raman Spectroscopic analysis indicates a consistent and uniform array Of CO2-bearing fluids. Comparison with active submarine hydrothermal vents shows an enrichment of the Mt Morgan fluids in base metals. Therefore, a seawater-dominated fluid is assumed for the barren massive sulfide mineralization, whereas magmatic volatile contributions are implied for the intrusive related mineralization. Condensation of magmatic vapor into a seawater-dominated environment explains the CO2 occurrence, the low salinities, and the enriched base and precious metal fluid composition that is associated with the Au-Cu. mineralization. The sulfur isotope signature of pyrite and chalcopyrite is composed of fractionated Devonian seawater and oxidized magmatic fluids or remobilized sulfur from existing sulfides. Pb isotopes indicate that Au and Cu. originated from the Mt Morgan intrusions and a particular volcanic strata that shows elevated Cu background. (C) 2002 Elsevier Science B.V. All rights reserved.

Drop penetration into porous powder beds

The kinetics of drop penetration were studied by filming single drops of several different fluids (water, PEG200, PEG600, and HPC solutions) as they penetrated into loosely packed beds of glass ballotini, lactose, zinc oxide, and titanium dioxide powders. Measured times ranged from 0.45 to 126 s and depended on the powder particle size,viscosity, surface tensions, and contact angle. The experimental drop penetration times were compared to existing theoretical predictions by M. Denesuk et al. (J. Colloid Interface Sci. 158, 114, 1993) and S. Middleman (Modeling Axisymmetric Flows: Dynamics of Films, Jets, and Drops, Academic Press, San Diego, 1995) but did not agree. Loosely packed powder beds tend to have a heterogeneous bed structure containing large macrovoids which do not participate in liquid flow but are included implicitly in the existing approach to estimating powder pore size. A new two-phase model was proposed where the total volume of the macrovoids was assumed to be the difference between the bed porosity and the tap porosity. A new parameter, the effective porosity (epsilon)eff, was defined as the tap porosity multiplied by the fraction of pores that terminate at a macrovoid and are effectively blocked pores. The improved drop penetration model was much more successful at estimating the drop penetration time on all powders and the predicted times were generally within an order of magnitude of the experimental results. (C) 2002 Elsevier Science (USA).

Instabilities and drop formation in cylindrical liquid jets in reduced gravity

The effects of convective and absolute instabilities on the formation of drops formed from cylindrical liquid jets of glycerol/water issuing into still air were investigated. Medium-duration reduced gravity tests were conducted aboard NASA's KC-135 and compared to similar tests performed under normal gravity conditions to aid in understanding the drop formation process. In reduced gravity, the Rayleigh-Chandrasekhar Equation was found to accurately predict the transition between a region of absolute and convective instability as defined by a critical Weber number. Observations of the physics of the jet, its breakup, and subsequent drop dynamics under both gravity conditions and the effects of the two instabilities on these processes are presented. All the normal gravity liquid jets investigated, in regions of convective or absolute instability, were subject to significant stretching effects, which affected the subsequent drop and associated geometry and dynamics. These effects were not displayed in reduced gravity and, therefore, the liquid jets would form drops which took longer to form (reduction in drop frequency), larger in size, and more spherical (surface tension effects). Most observed changes, in regions of either absolute or convective instabilities, were due to a reduction in the buoyancy force and an increased importance of the surface tension force acting on the liquid contained in the jet or formed drop. Reduced gravity environments allow better investigations to be performed into the physics of liquid jets, subsequently formed drops, and the effects of instabilities on these systems. In reduced gravity, drops form up to three times more slowly and as a consequence are up to three times larger in volume in the theoretical absolute instability region than in the theoretical convective instability region. This difference was not seen in the corresponding normal gravity tests due to the masking effects of gravity. A drop is shown to be able to form and detach in a region of absolute instability, and spanning the critical Weber number (from a region of convective to absolute instability) resulted in a marked change in dynamics and geometry of the liquid jet and detaching drops. (C) 2002 American Institute of Physics.

Close packed transitions in slit-shaped pores: Density functional theory study of methane adsorption capacity in carbon

An important feature of improving lattice gas models and classical isotherms is the incorporation of a pore size dependent capacity, which has hitherto been overlooked. In this paper, we develop a model for predicting the temperature dependent variation in capacity with pore size. The model is based on the analysis of a lattice gas model using a density functional theory approach at the close packed limit. Fluid-fluid and solid-fluid interactions are modeled by the Lennard-Jones 12-6 potential and Steele's 10-4-3, potential respectively. The capacity of methane in a slit-shaped carbon pore is calculated from the characteristic parameters of the unit cell, which are extracted by minimizing the grand potential of the unit cell. The capacities predicted by the proposed model are in good agreement with those obtained from grand canonical Monte Carlo simulation, for pores that can accommodate up to three adsorbed layers. Single particle and pair distributions exhibit characteristic features that correspond to the sequence of buckling and rhombic transitions that occur as the slit pore width is increased. The model provides a useful tool to model continuous variation in the microstructure of an adsorbed phase, namely buckling and rhombic transitions, with increasing pore width. (C) 2002 American Institute of Physics.

Physiological profiles of restricted endemic plants and their widespread congenors in the North Queensland wet tropics, Australia

Little is known about causes of endemic rarity in plants. This study pioneered an approach that determined environmental variables in the rainforest habitat and generated physiological profiles for light, water, and nutrient relations for three endemically restricted versus widespread congeneric species' pairs. We found no overall consistent differences in the physiological variables between the group of restricted species and the group of widespread species, and congeneric species pairs were therefore examined individually. Availability of soil nutrients did not differ between restricted-widespread species sites suggesting that species grow under comparable nutrient conditions. Under ambient and manipulated higher light conditions, widespread Gardenia ovularis had a greater photosynthetic activity than restricted Gardenia actinocarpa suggesting that the two species differ in their photosynthetic abilities. Differences between Xanthostemon species included lower photosynthetic activity, higher transpiration rate, and a higher foliar manganese concentration in restricted Xanthostemon formosus compared to widespread Xanthostemon chrysanthus. It is suggested that X. formosus is restricted by its high water use to its current rainforest creek edge habitat, while X. chrysanthus grows in a range of environments, although naturally found in riparian rainforest. Restricted Archidendron kanisii had higher electron transport rates, greater dissipative capacity for removal of excess light, and more efficient investment of nitrogen into photosynthetic components, than its widespread relative Archidendron whitei. These observations and previous research suggest that restricted Archidendron kanisii is in the process of expanding its range. Physiological profiles suggest a different cause of rarity for each species. This has implications for the conservation strategies required for each species. (C) 2002 Elsevier Science Ltd. 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.

The modeling of turbulent reactive flows based on multiple mapping conditioning

A new modeling approach-multiple mapping conditioning (MMC)-is introduced to treat mixing and reaction in turbulent flows. The model combines the advantages of the probability density function and the conditional moment closure methods and is based on a certain generalization of the mapping closure concept. An equivalent stochastic formulation of the MMC model is given. The validity of the closuring hypothesis of the model is demonstrated by a comparison with direct numerical simulation results for the three-stream mixing problem. (C) 2003 American Institute of Physics.

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.

Explicit/implicit partitioning and a new explicit form of the generalized alpha method

Most finite element packages use the Newmark algorithm for time integration of structural dynamics. Various algorithms have been proposed to better optimize the high frequency dissipation of this algorithm. Hulbert and Chung proposed both implicit and explicit forms of the generalized alpha method. The algorithms optimize high frequency dissipation effectively, and despite recent work on algorithms that possess momentum conserving/energy dissipative properties in a non-linear context, the generalized alpha method remains an efficient way to solve many problems, especially with adaptive timestep control. However, the implicit and explicit algorithms use incompatible parameter sets and cannot be used together in a spatial partition, whereas this can be done for the Newmark algorithm, as Hughes and Liu demonstrated, and for the HHT-alpha algorithm developed from it. The present paper shows that the explicit generalized alpha method can be rewritten so that it becomes compatible with the implicit form. All four algorithmic parameters can be matched between the explicit and implicit forms. An element interface between implicit and explicit partitions can then be used, analogous to that devised by Hughes and Liu to extend the Newmark method. The stability of the explicit/implicit algorithm is examined in a linear context and found to exceed that of the explicit partition. The element partition is significantly less dissipative of intermediate frequencies than one using the HHT-alpha method. The explicit algorithm can also be rewritten so that the discrete equation of motion evaluates forces from displacements and velocities found at the predicted mid-point of a cycle. Copyright (C) 2003 John Wiley Sons, Ltd.

Skin-friction measurements in high-enthalpy hypersonic boundary layers

Skin-friction measurements are reported for high-enthalpy and high-Mach-number laminar, transitional and turbulent boundary layers. The measurements were performed in a free-piston shock tunnel with air-flow Mach number, stagnation enthalpy and Reynolds numbers in the ranges of 4.4-6.7, 3-13 MJ kg(-1) and 0.16 x 10(6)-21 x 10(6), respectively. Wall temperatures were near 300 K and this resulted in ratios of wall enthalpy to flow-stagnation enthalpy in the range of 0.1-0.02. The experiments were performed using rectangular ducts. The measurements were accomplished using a new skin-friction gauge that was developed for impulse facility testing. The gauge was an acceleration compensated piezoelectric transducer and had a lowest natural frequency near 40 kHz. Turbulent skin-friction levels were measured to within a typical uncertainty of +/-7%. The systematic uncertainty in measured skin-friction coefficient was high for the tested laminar conditions; however, to within experimental uncertainty, the skin-friction and heat-transfer measurements were in agreement with the laminar theory of van Driest (1952). For predicting turbulent skin-friction coefficient, it was established that, for the range of Mach numbers and Reynolds numbers of the experiments, with cold walls and boundary layers approaching the turbulent equilibrium state, the Spalding & Chi (1964) method was the most suitable of the theories tested. It was also established that if the heat transfer rate to the wall is to be predicted, then the Spalding & Chi (1964) method should be used in conjunction with a Reynolds analogy factor near unity. If more accurate results are required, then an experimentally observed relationship between the Reynolds analogy factor and the skin-friction coefficient may be applied.

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.

Anomalous water absorption in porous materials

The absorption of fluid by unsaturated, rigid porous materials may be characterized by the sorptivity. This is a simple parameter to determine and is increasingly being used as a measure of a material's resistance to exposure to fluids (especially moisture and reactive solutes) in aggressive environments. The complete isothermal absorption process is described by a nonlinear diffusion equation, with the hydraulic diffusivity being a strongly nonlinear function of the degree of saturation of the material. This diffusivity can be estimated from the sorptivity test. In a typical test the cumulative absorption is proportional to the square root of time. However, a number of researchers have observed deviation from this behaviour when the infiltrating fluid is water and there is some potential for chemo-mechanical interaction with the material. In that case the current interpretation of the test and estimation of the hydraulic diffusivity is no longer appropriate. Kuntz and Lavallee (2001) discuss the anomalous behaviour and propose a non-Darcian model as a more appropriate physical description. We present an alternative Darcian explanation and theory that retrieves the earlier advantages of the simple sorptivity test in providing parametric information about the material's hydraulic properties and allowing simple predictive formulae for the wetting profile to be generated.

Method for determining the shear stress in cylindrical systems

We develop a method for determining the elements of the pressure tensor at a radius r in a cylindrically symmetric system, analogous to the so-called method of planes used in planar systems [B. D. Todd, Denis J. Evans, and Peter J. Daivis, Phys. Rev. E 52, 1627 (1995)]. We demonstrate its application in determining the radial shear stress dependence during molecular dynamics simulations of the forced flow of methane in cylindrical silica mesopores. Such expressions are useful for the examination of constitutive relations in the context of transport in confined systems.