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.

Influence of Mg content on the microstructure and solid solution chemistry of Al-7%Si-Mg casting alloys during solution treatment

The solution treatment stage of the T6 heat-treatment of Al-7%Si-Mg foundry alloys influences microstructural features such as Mg2Si dissolution, and eutectic silicon spheroidisation and coarsening. Microstructural and microanalytical studies have been conducted across a range of Sr-modified Al-7%Si alloys, with an Fe content of 0.12% and Mg contents ranging from 0.3-0.7wt%. Qualitative and quantitative metallography have shown that, in addition to the above changes, solution treatment also results in changes to the relative proportions of iron-containing intermetallic particles and that these changes are composition-dependent. While solution treatment causes a substantial transformation of pi phase to beta phase in low Mg alloys (0.3-0.4%), this change is not readily apparent at higher Mg levels (0.6-0.7%). The pi to beta transformation is accompanied by a release of Mg into the aluminum matrix over and above that which arises from the rapid dissolution of Mg2Si. Since the level of matrix Mg retained after quenching controls an alloy's subsequent precipitation hardening response, a proper understanding of this phase transformation is crucial if tensile properties are to be maximised.

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.

Change in the rate of shell deposition and shell microstructure in response to shell borers in the abalone Haliotis rubra

Gastropod shells consist of two crystal types of calcium carbonate, an outer, prismatic calcite layer and an inner nacreous layer made of aragonite. In cross-section, the nacre of the nacreous layer appears to have a regular brick-like microstructure composed of thin laminae of aragonite crystals, separated by very thin sheets of protein (Lutz and Rhoads, 1980; Nakahara, 1983). In abalone (Genus, Haliotis) and other gastropods, thin layers of non-lamellar pigmented material occur within the nacre and have been termed alternatively, fine lines, growth rings or growth lines (Shepherd et al., 1995). It has been suggested that these pigmented layers are small, prismatic, calcite layers (Shepherd and Avalos-Borja, 1997; Zaremba et al., 1996) but investigations using a Raman laser in Haliotis rubra show that they contain aragonite rather than calcite (Hawkes et al, 1996). Day and Fleming (1992) suggest that the occurrence of pigmented layers is correlated with regular exogenous cues such as reproduction or temperature changes and indeed in some species, pigmented layers in the shell can be used to age abalone (review: Shepherd and Triantafillos, 1997). However, McShane and Smith (1992) suggest that pigmented layers can occur irregularly and therefore may be unreliable indicators of age.

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.

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.

Eutectic modification and microstructure development in Al-Si alloys

Recent increasing applications for cast Al-Si alloys are particularly driven by the need for lightweighting components in the automotive sector. To improve mechanical properties, elements such as strontium, sodium and antimony can be added to modify the eutectic silicon from coarse and plate-like to fine and fibrous morphology. It is only recently being noticed that the morphological transformation resulting from eutectic modification is also accompanied by other, equally significant, but often unexpected changes. These changes can include a 10-fold increase in the eutectic grain size, redistribution of low-melting point phases and porosity as well as surface finish, consequently leading to variations in casting quality. This paper shows the state-of-the-art in understanding the mechanism of eutectic nucleation and growth in Al-Si alloys, inspecting samples, both quenched and uninterrupted, on the macro, micro and nano-scale. It shows that significant variations in eutectic nucleation and growth dynamics occur in AI-Si alloys as a function of the type and amount of modifier elements added. The key role of AIP particles in nucleating silicon is demonstrated. (c) 2005 Elsevier B.V. All rights reserved.

The microstructure of mechanically alloyed MoSi2 with Ni and Al additions

The effect of Ni and Al additions on grain boundary silica in mechanically alloyed and hot isostatically pressed (HTPed) MoSi2 was investigated. Mechanical alloying Mo and Si in the absence of Al produced finely dispersed silica within a fine grained structure. Mechanically alloyed and HIPed Mo and Si with Ni and Al partially transformed the silica to crystalline oxide phases, including Al2O3. An improvement in high temperature properties is not expected due to the retention of a grain boundary silica film. Rapid grain growth resulted during HIPing, possibly due to the formation of a Ni/Fe/Al liquid phase.