Effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in deformable fluid-saturated porous media heated from below

In this paper, a solution method is presented to deal with fully coupled problems between medium deformation, pore-fluid flow and heat transfer in fluid-saturated porous media having supercritical Rayleigh numbers. To validate the present solution method, analytical solutions to a benchmark problem are derived for some special cases. After the solution method is validated, a numerical study is carried out to investigate the effects of medium thermoelasticity on high Rayleigh number steady-state heat transfer and mineralization in fluid-saturated media when they are heated from below. The related numerical results have demonstrated that: (1) medium thermoelasticity has a little influence on the overall pattern of convective pore-fluid flow, but it has a considerable effect on the localization of medium deformation, pore-fluid flow, heat transfer and mineralization in a porous medium, especially when the porous medium is comprised of soft rock masses; (2) convective pore-fluid flow plays a very important role in the localization of medium deformation, heat transfer and mineralization in a porous medium. (C) 1999 Elsevier Science S.A. All rights reserved.

A numerical study of pore-fluid, thermal and mass flow iu fluid-saturated porous rock basins

We present a numerical methodology for the study of convective pore-fluid, thermal and mass flow in fluid-saturated porous rock basins. lit particular, we investigate the occurrence and distribution pattern of temperature gradient driven convective pore-fluid flow and hydrocarbon transport in the Australian North West Shelf basin. The related numerical results have demonstrated that: (1) The finite element method combined with the progressive asymptotic approach procedure is a useful tool for dealing with temperature gradient driven pore-fluid flow and mass transport in fluid-saturated hydrothermal basins; (2) Convective pore-fluid flow generally becomes focused in more permeable layers, especially when the layers are thick enough to accommodate the appropriate convective cells; (3) Large dislocation of strata has a significant influence off the distribution patterns of convective pore;fluid flow, thermal flow and hydrocarbon transport in the North West Shelf basin; (4) As a direct consequence of the formation of convective pore-fluid cells, the hydrocarbon concentration is highly localized in the range bounded by two major faults in the basin.

Theoretical and numerical analyses of convective instability in porous media with upward throughflow

Exact analytical solutions have been obtained for a hydrothermal system consisting of a horizontal porous layer with upward throughflow. The boundary conditions considered are constant temperature, constant pressure at the top, and constant vertical temperature gradient, constant Darcy velocity at the bottom of the layer. After deriving the exact analytical solutions, we examine the stability of the solutions using linear stability theory and the Galerkin method. It has been found that the exact solutions for such a hydrothermal system become unstable when the Rayleigh number of the system is equal to or greater than the corresponding critical Rayleigh number. For small and moderate Peclet numbers (Pe less than or equal to 6), an increase in upward throughflow destabilizes the convective flow in the horizontal layer. To confirm these findings, the finite element method with the progressive asymptotic approach procedure is used to compute the convective cells in such a hydrothermal system. Copyright (C) 1999 John Wiley & Sons, Ltd.

Finite element modelling of reactive mass transport problems in fluid-saturated porous media

We use the finite element method to solve reactive mass transport problems in fluid-saturated porous media. In particular, we discuss the mathematical expression of the chemical reaction terms involved in the mass transport equations for an isothermal, non-equilibrium chemical reaction. It has turned out that the Arrhenius law in chemistry is a good mathematical expression for such non-equilibrium chemical reactions especially from the computational point of view. Using the finite element method and the Arrhenius law, we investigate the distributions of PH (i.e. the concentration of H+) and the relevant reactive species in a groundwater system. Although the main focus of this study is on the contaminant transport problems in groundwater systems, the related numerical techniques and principles are equally applicable to the orebody formation problems in the geosciences. Copyright (C) 1999 John Wiley & Sons, Ltd.

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.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Finite element modelling of three-dimensional convection problems in fluid-saturated porous media heated from below

We use the finite element method to model three-dimensional convective pore-fluid flow in fluid-saturated porous media when they are heated from below. In particular, we employ the particle-tracking technique to mimic the trajectories of particles in three-dimensional fluid flow problems. The related numerical results demonstrated that: (1) The progressive asymptotic approach procedure, which was previously developed for the finite element modelling of two-dimensional convective pore-fluid flow problems, is equally applicable to the finite element modelling of three-dimensional convective pore-fluid flow in fluid-saturated porous media heated from below. (2) The perturbation of gravity at different planes has a significant effect on the pattern of three-dimensional convective pore-fluid flow and therefore, may influence the pattern of orebody formation and mineralization in three-dimensional hydrothermal systems. Copyright (C) 2001 John Wiley & Sons, Ltd.

Synthesis of anatase TiO2 supported on porous solids by chemical vapor deposition

Coating anatase TiO2 onto three different particle supports, activated carbon (AC), gamma -alumina (Al2O3) and silica gel (SiO2), by chemical vapor deposition (CVD) was studied. The effect of the CVD synthesis conditions on the loading rate of anatase TiO2 was investigated. It was found that introducing water vapor during CVD or adsorbing water before CVD was crucial to obtain anatase TiO2 on the surface of the particle supports. The evaporation temperature of precursor, deposition temperature in the reactor, flow rate of carrier gas, and the length of coating time were also important parameters to obtain more uniform and repeatable TiO2 coating. High inflow precursor concentration, high CVD reactor temperature and long coating time tended to cause block problem. Coating TiO2 onto small particles by CVD involved both chemical vapor deposition and particle deposition. It was believed that the latter was the reason for the block problem. In addition, the mechanism of CVD process in this study included two parts, pyrolysis and hydrolysis, and one of them was dominant in the CVD process under different synthesis route. Among the three types of materials, silica gel, with higher surface hydroxyl groups and macropore surface area, was found to be the most efficient support in terms of both anatase TiO2 coating and photocatalytic reaction. (C) 2001 Elsevier Science B.V. All rights reserved.

Model for gas-liquid flow through wet porous medium

A method involving bubbling of air through a fibrous filter immersed in water has recently been investigated (Agranovski et al. [1]). Experimental results showed that the removal efficiency for ultra-fine aerosols by such filters was greatly increased compared to dry filters. Nuclear Magnetic Resonance (NMR) imaging was used to examine the wet filter and to determine the nature of the gas flow inside the filter (Agranovski et al. [2]). It was found that tortuous preferential pathways (or flow tubes) develop within the filter through which the air flows and the distribution of air and water inside the porous medium has been investigated. The aim of this paper is to investigate the geometry of the pathways and to make estimates of the flow velocities and particle removal efficiency in such pathways. A mathematical model of the flow of air along the preferred pathways has been developed and verified experimentally. Even for the highest realistic gas velocity the flow field was essentially laminar (Re approximate to 250). We solved Laplace's equation for stream function to map trajectories of particles and gas molecules to investigate the possibility of their removal from the carrier.

Porous clays and pillared clays-based catalysts. Part 2: A review of the catalytic and molecular sieve applications

Metal oxide pillared clay (PILC) possesses several interesting properties, such as large surface area, high pore volume and tunable pore size (from micropore to mesopore), high thermal stability, strong surface acidity and catalytic active substrates/metal oxide pillars. These unique characteristics make PILC an attractive material in catalytic reactions. It can be made either as catalyst support or directly used as catalyst. This paper is a continuous work from Kloprogge's review (J.T. Kloprogge, J. Porous Mater. 5, 5 1998) on the synthesis and properties of smectites and related PILCs and will focus on the diverse applications of clay pillared with different types of metal oxides in the heterogeneous catalysis area and adsorption area. The relation between the performance of the PILC and its physico-chemical features will be addressed.

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.

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.

Systemic Model for Diagnosis of the Micro, Small and Medium Enterprises from Two Cities from the Countryside of the State of So Paulo in Brazil

This article presents a proposal of a systemic model composed for the micro and small companies (MSE) of the region of Ribeiro Preto and the agents which influenced their environment. The proposed model was based on Stafford Beer`s (Diagnosing the system for organizations. Chichester, Wiley, 1985) systemic methodologies VSM (Viable System Model) and on Werner Ulrich`s (1983) CSH (Critical Systems Heuristics). The VSM is a model for the diagnosis of the structure of an organization and of its flows of information through the application of the cybernetics concepts (Narvarte, In El Modelo del Sistema Viable-MSV: experiencias de su aplicacin en Chile. Proyecto Cerebro Colectivo del IAS, Santiago, 2001). On the other hand, CSH focus on the context of the social group applied to the systemic vision as a counterpoint to the organizational management view considered by the VSM. MSE of Ribeiro Preto and Sertozinho had been analyzed as organizations inserted in systems that relate and integrate with other systems concerning the public administration, entities of representation and promotion agencies. The research questions: which are the bonds of interaction among the subsystems in this process and who are the agents involved? The systemic approach not only diagnosed a social group, formed by MSE of Ribeiro Preto and Sertozinho, public authorities and support entities, but could also delineate answers that aimed the clarification of obscure questions generating financial assistance to the formularization of efficient actions for the development of this system.

Machining of micro aspherical mould inserts

This paper reports a parametric investigation and development of grinding technologies for micro aspherical mould inserts using parallel grinding method. The parametric investigation revealed that at nanometric scale the undeformed chip thickness has little influence on the surface finish of ground inserts. The grinding trace spacing has a slightly larger influence on the surface finish. A new technique was developed to true and dress the resin bonded micro wheels with mesh size of #3000, which produced a satisfactory wheel form accuracy and relatively high grain packing density. A form error compensation technique was also developed, with which mould inserts of submicron form accuracy were consistently produced. Using the developed technologies, micro aspherical inserts of diameters ranging from 200 mu m to 1000 mu m with surface finish of around 10 nm and form error of similar to 0.2-0.4 mu m were successfully fabricated. (c) 2005 Elsevier Inc. All rights reserved.

Numerical modelling of chemical effects of magma solidification problems in porous rocks

The solidification of intruded magma in porous rocks can result in the following two consequences: (1) the heat release due to the solidification of the interface between the rock and intruded magma and (2) the mass release of the volatile fluids in the region where the intruded magma is solidified into the rock. Traditionally, the intruded magma solidification problem is treated as a moving interface (i.e. the solidification interface between the rock and intruded magma) problem to consider these consequences in conventional numerical methods. This paper presents an alternative new approach to simulate thermal and chemical consequences/effects of magma intrusion in geological systems, which are composed of porous rocks. In the proposed new approach and algorithm, 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 the proposed mass source and physically equivalent heat source. The major advantage in using the proposed equivalent algorithm is that a fixed mesh of finite elements with a variable integration time-step can be employed to simulate the consequences and effects of the intruded magma solidification using the conventional finite element method. The correctness and usefulness of the proposed equivalent algorithm have been demonstrated by a benchmark magma solidification problem. Copyright (c) 2005 John Wiley & Sons, Ltd.