Finite element analysis of heat transfer and mineralization in layered hydrothermal systems with upward throughflow

We use the finite element method to solve the coupled problem between convective pore-fluid flow, heat transfer and mineralization in layered hydrothermal systems with upward throughflow. In particular, we present the improved rock alteration index (IRAI) concept for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in the systems. To validate the numerical method used in the computation, analytical solutions to a benchmark problem have been derived. After the numerical method is validated, it is used to investigate the pattern of pore-fluid Aom, the distribution of temperature and the mineralization pattern of gold minerals in a layered hydrothermal system with upward throughflow. The related numerical results have demonstrated that the present concept of IRAI is useful and applicable for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in hydrothermal systems. (C) 2000 Elsevier Science S.A. All rights reserved.

The mechanical properties of Saccharomyces cerevisiae

Cell-wall mechanical properties play an integral part in the growth and form of Saccharomyces cerevisiae, In contrast to the tremendous knowledge on the genetics of S. cerevisiae, almost nothing is known about its mechanical properties. We have developed a micromanipulation technique to measure the force required to burst single cells and have recently established a mathematical model to extract the mechanical properties of the cell wall from such data, Here we determine the average surface modulus of the S, cerevisiae cell wall to be 11.1 +/- 0.6 N/m and 12.9 +/- 0.7 N/m in exponential and stationary phases, respectively, giving corresponding Young's moduli of 112 +/- 6 MPa and 107 +/- 6 MPa, This result demonstrates that yeast cell populations strengthen as they enter stationary phase by increasing wall thickness and hence the surface modulus, without altering the average elastic properties of the cell-wall material. We also determined the average breaking strain of the cell wall to be 82% +/- 3% in exponential phase and 80% +/- 3% in stationary phase, This finding provides a failure criterion that can be used to predict when applied stresses (e,g,, because of fluid flow) will lead to wall rupture, This work analyzes yeast compression experiments in different growth phases by using engineering methodology.

Prediction of rock alteration patterns: a potential tool in mineral exploration

Hydrothermal alteration of a quartz-K-feldspar rock is simulated numerically by coupling fluid flow and chemical reactions. Introduction of CO2 gas generates an acidic fluid and produces secondary quartz, muscovite and/or pyrophyllite at constant temperature and pressure of 300 degrees C and 200 MPa. The precipitation and/or dissolution of the secondary minerals is controlled by either mass-action relations or rate laws. In our simulations the mass of the primary elements are conserved and the mass-balance equations are solved sequentially using an implicit scheme in a finite-element code. The pore-fluid velocity is assumed to be constant. The change of rock volume due to the dissolution or precipitation of the minerals, which is directly related to their molar volume, is taken into account. Feedback into the rock porosity and the reaction rates is included in the model. The model produces zones of pyrophyllite quartz and muscovite due to the dissolution of K-feldspar. Our model simulates, in a simplified way, the acid-induced alteration assemblages observed in various guises in many significant mineral deposits. The particular aluminosilicate minerals produced in these experiments are associated with the gold deposits of the Witwatersrand Basin.

Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

A simple solution of the laminar dam break wave

The sudden release of a mass of fluid in a channel generates a highly unsteady flow motion, called dam break wave. While industrial fluids exhibit sometimes non-Newtonian behaviours, the viscous fluid flow assumption remains a useful approximation for simplified analyses. In this study, new solutions of laminar dam break wave are proposed for a semi-infinite reservoir based upon the method of characteristics. The solutions yield simple explicit expressions of the wave front location, wave front celerity and instantaneous free-surface profiles that compare favourably with experimental observations. Both horizontal and sloping channel configurations are treated. The simplicity of the equations may allow future extension to more complicated fluid flows.

Propagation of nonstationary curved and stretched premixed flames with multistep reaction mechanisms

The propagation speed of a thin premixed flame disturbed by an unsteady fluid flow of a larger scale is considered. The flame may also have a general shape but the reaction zone is assumed to be thin compared to the flame thickness. Unlike in preceding publications, the presented asymptotic analysis is performed for a general multistep reaction mechanism and, at the same time, the flame front is curved by the fluid flow. The resulting equations define the propagation speed of disturbed flames in terms of the properties of undisturbed planar flames and the flame stretch. Special attention is paid to the near-equidiffusion limit. In this case, the flame propagation speed is shown to depend on the effective Zeldovich number Z(f) , and the flame stretch. Unlike the conventional Zeldovich number, the effective Zeldovich number is not necessarily linked directly to the activation energies of the reactions. Several examples of determining the effective Zeldovich number for reduced combustion mechanisms are given while, for realistic reactions, the effective Zeldovich number is determined from experiments. Another feature of the present approach is represented by the relatively simple asymptotic technique based on the adaptive generalized curvilinear system of coordinates attached to the flame (i.e., intrinsic disturbed flame equations [IDFE]).

Molecular transport in nanopores

Simulation of the transport of methane in cylindrical silica mesopores have been performed using equilibrium and nonequilibrium molecular dynamics (NEMD) as well as dual control volume grand canonical molecular dynamics methods. It is demonstrated that all three techniques yield the same transport coefficient even in the presence of viscous flow. A modified locally averaged density model for viscous flow, combined with consideration of wall slip through a frictional condition, gives a convincing interpretation of the variation of the transport coefficient over a wide range of densities, and for various pore sizes and temperatures. Wall friction coefficients extracted from NEMD simulations are found to be consistent with momentum transfer arguments, and the approach is shown to be more meaningful than the classical slip length concept. (C) 2003 American Institute of Physics.

Rare earth element fractionation in authigenic illite–smectite from Late Permian clastic rocks, Bowen Basin, Australia: implications for physico-chemical environments of fluids during illitization

REE analyses were performed on authigenic illitic clay. minerals from Late Permian mudrocks, sandstones and bentonites from the Bowen Basin (Australia). The mixed-layer illite-smectite exhibit REE patterns with an obvious fractionation of the HREE from the LREE and MREE, which is an apparent function of degree of illitization reaction. The highly illitic (R greater than or equal to 3) illite-smectite from the northern Bowen Basin show a depletion of LREE relative to the less illitic (R=0 and 1) clays. In contrast, an enrichment of HREE for the illite-rich clays relative to less. illitic clays is evident for the southern Bowen Basin samples. The North American Shale Composite-normalized (La/Lu)(sn) ratios show negative correlations with the illite content in illite-smectite and positive correlations with the delta(18)O values of the clays for both the northern and southern Bowen Basin samples. These correlations indicate that the increasing depletion of LREE in hydrothermal fluids is a function of increasing water/rock ratios in the northern Bowen Basin. Good negative correlations between (La/Lu)(sn) ratios and illite content in illite-smectite from the southern Bowen Basin suggest the involvement of fluids with higher alkalinity and higher pH in low water/ rock ratio conditions. Increasing HREE enrichment with delta(18)O decrease indicates the effect of increasing temperature at low water/rock ratios in the southern Bowen Basin. Results of the present study confirm the conclusions of some earlier studies suggesting that REE in illitic clay minerals are mobile and fractionated during illitization and that this fact should be considered in studies of sedimentary processes and in identifying provenance. Moreover, our results show that REE systematic of illitic clay minerals can be applied as an useful technique to gain information about physico-chemical conditions during thermal and fluid flow events in certain sedimentary basins. (C) 2003 Elsevier Science B.V. All rights reserved.

Three-dimensional simulation of the rolling of a saturated fibro-porous media

This paper describes recent advances made in computational modelling of the sugar cane liquid extraction process. The saturated fibro-porous material is rolled between circumferentially grooved rolls, which enhance frictional grip and provide a low-resistance path for liquid flow during the extraction process. Previously reported two-dimensional (2D) computational models, account for the large deformation of the porous material by solving the fully coupled governing fibre stress and fluid-flow equations using finite element techniques. While the 2D simulations provide much insight into the overarching cause-effect relationships, predictions of mechanical quantities such as roll separating force and particularly torque as a function of roll speed and degree of compression are not satisfactory for industrial use. It is considered that the unsatisfactory response in roll torque prediction may be due to the stress levels that exist between the groove tips and roots which have been largely neglected in the geometrically simplified 2D model. This paper gives results for both two- and three-dimensional finite element models and highlights their strengths and weaknesses in predicting key milling parameters. (c) 2005 Elsevier B.V. All rights reserved.

Modeling of Flowing Film Concentrators - Part 1. Water Split Behaviour

A large number of mineral processing equipment employs the basic principles of gravity concentration in a flowing fluid of a few millimetres thick in small open channels where the particles are distributed along the flow height based on their physical properties and the fluid flow characteristics. Fluid flow behaviour and slurry transportation characteristics in open channels have been the research topic for many years in many engineering disciplines. However, the open channels used in the mineral processing industries are different in terms of the size of the channel and the flow velocity used. Understanding of water split behaviour is, therefore, essential in modeling flowing film concentrators. In this paper, an attempt has been made to model the water split behaviour in an inclined open rectangular channel, resembling the actual size and the flow velocity used by the mineral processing industries, based on the Prandtl's mixing length approach. (c) 2006 Elsevier B.V. All rights reserved.

Numerical simulation of magnetite segregation in a dense medium cyclone

Numerical simulations of turbulent driven flow in a dense medium cyclone with magnetite medium have been conducted using Fluent. The predicted air core shape and diameter were found to be close to the experimental results measured by gamma ray tomography. It is possible that the Large eddy simulation (LES) turbulence model with Mixture multi-phase model can be used to predict the air/slurry interface accurately although the LES may need a finer grid. Multi-phase simulations (air/water/medium) are showing appropriate medium segregation effects but are over-predicting the level of segregation compared to that measured by gamma-ray tomography in particular with over prediction of medium concentrations near the wall. Further, investigated the accurate prediction of axial segregation of magnetite using the LES turbulence model together with the multi-phase mixture model and viscosity corrections according to the feed particle loading factor. Addition of lift forces and viscosity correction improved the predictions especially near the wall. Predicted density profiles are very close to gamma ray tomography data showing a clear density drop near the wall. The effect of size distribution of the magnetite has been fully studied. It is interesting to note that the ultra-fine magnetite sizes (i.e. 2 and 7 mu m) are distributed uniformly throughout the cyclone. As the size of magnetite increases, more segregation of magnetite occurs close to the wall. The cut-density (d(50)) of the magnetite segregation is 32 gm, which is expected with superfine magnetite feed size distribution. At higher feed densities the agreement between the [Dungilson, 1999; Wood, J.C., 1990. A performance model for coal-washing dense medium cyclones, Ph.D. Thesis, JKMRC, University of Queensland] correlations and the CFD are reasonably good, but the overflow density is lower than the model predictions. It is believed that the excessive underflow volumetric flow rates are responsible for under prediction of the overflow density. (c) 2006 Elsevier Ltd. All rights reserved.

CFD simulations of hydrocyclones with an air core - Comparison between large eddy simulations and a second moment closure

CFD simulations of the 75 mm, hydrocyclone of Hsieh (1988) have been conducted using Fluent TM. The simulations used 3-dimensional body fitted grids. The simulations were two phase simulations where the air core was resolved using the mixture (Manninen et al., 1996) and VOF (Hirt and Nichols, 1981) models. Velocity predictions from large eddy simulations (LES), using the Smagorinsky-Lilly sub grid scale model (Smagorinsky, 1963; Lilly, 1966) and RANS simulations using the differential Reynolds stress turbulence model (Launder et al., 1975) were compared with Hsieh's experimental velocity data. The LES simulations gave very good agreement with Hsieh's data but required very fine grids to predict the velocities correctly in the bottom of the apex. The DRSM/RANS simulations under predicted tangential velocities, and there was little difference between the velocity predictions using the linear (Launder, 1989) and quadratic (Speziale et al., 1991) pressure strain models. Velocity predictions using the DRSM turbulence model and the linear pressure strain model could be improved by adjusting the pressure strain model constants.

Experimental and analytical study of the effect of contact angle on liquid convective heat transfer in microchannels

In this paper we examine the effect of contact angle (or surface wettability) on the convective heat transfer coefficient in microchannels. Slip flow, where the fluid velocity at the wall is non-zero, is most likely to occur in microchannels due to its dependence on shear rate or wall shear stress. We show analytically that for a constant pressure drop, the presence of slip increases the Nusselt number. In a microchannel heat exchanger we modified the surface wettability from a contact angle of 20 degrees-120 degrees using thin film coating technology. Apparent slip flow is implied from pressure and flow rate measurements with a departure from classical laminar friction coefficients above a critical shear rate of approximately 10,000 s(-1). The magnitude of this departure is dependant on the contact angle with higher contact angles surfaces exhibiting larger pressure drop decreases. Similarly, the non-dimensional heat flux is found to decrease relative to laminar non-slip theory, and this decrease is also a function of the contact angle. Depending on the contact angle and the wall shear rate, variations in the heat transfer rate exceeding 10% can be expected. Thus the contact angle is an important consideration in the design of micro, and even more so, nano heat exchangers. (c) 2006 Elsevier Ltd. All rights reserved.