Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

Oxygen transport capacity in the air-breathing fish, Megalops cyprinoides: compensations for strenuous exercise

Tarpon have high resting or routine hematocrits (Hct) (37.6+/-3.4%) and hemoglobin concentrations (120.6+/-7.3 g 1(-1)) that increased significantly following bouts of angling-induced exercise (51.9+/-3.7% and 142.8+/-13.5 g 1(-1), respectively). Strenuous exercise was accompanied by an approximately tenfold increase in blood lactate and a muscle metabolite profile indicative of a high energy demand teleost. Routine blood values were quickly restored only when this facultative air-breathing fish was given access to atmospheric air. In vitro studies of oxygen transport capacity, a function of carrying capacity and viscosity, revealed that the optimal Hct range corresponded to that observed in fish under routine behaviour. During strenuous exercise however, further increase in viscosity was largely offset by a pronounced reduction in the shear-dependence of blood which conformed closely to an ideal Newtonian fluid. The mechanism for this behaviour of the erythrocytes appears to involve the activation of surface adrenergic receptors because pre-treatment with propranolol abolished the response. High levels of activity in tarpon living in hypoxic habitats are therefore supported by an elevated Hct with adrenergically mediated viscosity reduction, and air-breathing behaviour that enables rapid metabolic recovery. (C) 2002 Elsevier Science Inc. All rights reserved.

α′ phase stability in Nd-Li-sialon systems

The formability and stability of the alpha-sialon (alpha') phase was investigated in multi-cation Nd-Li-sialon systems. Four samples were prepared, ranging from a pure Nd-sialon to a pure Li-sialon, with two intermediate samples being prepared with either lithium or neodymium replacing the other alpha'-stabilising additive by 20 eq.%, as to maintain an equivalent design composition in all samples. After sintering, all samples were subsequently heat treated up to 192 h at 1450 and 1300 degreesC. While significant quantities of the beta'-sialon (beta' phase were found in most samples, the high-lithium Li-Nd-sialon sample was found to be almost pure a' phase after sintering. Furthermore, the long-term stability of the a' phase on heat treatment was also found to be superior in both multi-cation samples than in either of the single-alpha'-stabilising-cation samples. This is thought to be related to improved retention of the lithium in the multi-cation systems, as much of the lithium was found to volatilise during sintering in the neodymium-free sample. (C) 2002 Elsevier Science Ltd. All rights reserved.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.

Special issue: Selected papers from the sixth IFIP 2.6 working conference on visual database systems: Guest editors' introduction

While multimedia data, image data in particular, is an integral part of most websites and web documents, our quest for information so far is still restricted to text based search. To explore the World Wide Web more effectively, especially its rich repository of truly multimedia information, we are facing a number of challenging problems. Firstly, we face the ambiguous and highly subjective nature of defining image semantics and similarity. Secondly, multimedia data could come from highly diversified sources, as a result of automatic image capturing and generation processes. Finally, multimedia information exists in decentralised sources over the Web, making it difficult to use conventional content-based image retrieval (CBIR) techniques for effective and efficient search. In this special issue, we present a collection of five papers on visual and multimedia information management and retrieval topics, addressing some aspects of these challenges. These papers have been selected from the conference proceedings (Kluwer Academic Publishers, ISBN: 1-4020- 7060-8) of the Sixth IFIP 2.6 Working Conference on Visual Database Systems (VDB6), held in Brisbane, Australia, on 29–31 May 2002.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.

Designing teams for first-of-a-kind, complex systems using the initial phases of cognitive work analysis: Case study

We present a technique for team design based on cognitive work analysis (CWA). We first develop a rationale for this technique by discussing the limitations of conventional approaches for team design in light of the special characteristics of first-of-a-kind, complex systems. We then introduce the CWA-based technique for team design and provide a case study of how we used this technique to design a team for a first-of-a-kind, complex military system during the early stages of its development. In addition to illustrating the CWA-based technique by example, the case study allows us to evaluate the technique. This case study demonstrates that the CWA-based technique for team design is both feasible and useful, although empirical validation of the technique is still necessary. Applications of this work include the design of teams for first-of-a-kind, complex systems in military, medical, and industrial domains.

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.

Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity

Exact analytical solutions of the critical Rayleigh numbers have been obtained for a hydrothermal system consisting of a horizontal porous layer with temperature-dependent viscosity. The boundary conditions considered are constant temperature and zero vertical Darcy velocity at both the top and bottom of the layer. Not only can the derived analytical solutions be readily used to examine the effect of the temperature-dependent viscosity on the temperature-gradient driven convective flow, but also they can be used to validate the numerical methods such as the finite-element method and finite-difference method for dealing with the same kind of problem. The related analytical and numerical results demonstrated that the temperature-dependent viscosity destabilizes the temperature-gradient driven convective flow and therefore, may affect the ore body formation and mineralization in the upper crust of the Earth. Copyright (C) 2003 John Wiley Sons, Ltd.

Characterisation of the three-dimensional structure of earthworm burrow systems using image analysis and mathematical morphology

The aim of this work was to exemplify the specific contribution of both two- and three-dimensional (31)) X-ray computed tomography to characterise earthworm burrow systems. To achieve this purpose we used 3D mathematical morphology operators to characterise burrow systems resulting from the activity of an anecic (Aporrectodea noctunia), and an endogeic species (Allolobophora chlorotica), when both species were introduced either separately or together into artificial soil cores. Images of these soil cores were obtained using a medical X-ray tomography scanner. Three-dimensional reconstructions of burrow systems were obtained using a specifically developed segmentation algorithm. To study the differences between burrow systems, a set of classical tools of mathematical morphology (granulometries) were used. So-called granulometries based on different structuring elements clearly separated the different burrow systems. They enabled us to show that burrows made by the anecic species were fatter, longer, more vertical, more continuous but less sinuous than burrows of the endogeic species. The granulometry transform of the soil matrix showed that burrows made by A. nocturna were more evenly distributed than those of A. chlorotica. Although a good discrimination was possible when only one species was introduced into the soil cores, it was not possible to separate burrows of the two species from each other in cases where species were introduced into the same soil core. This limitation, partly due to the insufficient spatial resolution of the medical scanner, precluded the use of the morphological operators to study putative interactions between the two species.

Parity measurement of one- and two-electron double well systems

We outline a scheme to accomplish measurements of a solid state double well system (DWS) with both one and two electrons in nonlocalized bases. We show that, for a single particle, measuring the local charge distribution at the midpoint of a DWS using a SET as a sensitive electrometer amounts to performing a projective measurement in the parity (symmetric/antisymmetric) eigenbasis. For two-electrons in a DWS, a similar configuration of SET results in close-to-projective measurement in the singlet/triplet basis. We analyze the sensitivity of the scheme to asymmetry in the SET position for some experimentally relevant parameter, and show that it is experimentally realizable.