Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses

We present finite element simulations of temperature gradient driven rock alteration and mineralization in fluid saturated porous rock masses. In particular, we explore the significance of production/annihilation terms in the mass balance equations and the dependence of the spatial patterns of rock alteration upon the ratio of the roll over time of large scale convection cells to the relaxation time of the chemical reactions. Special concepts such as the gradient reaction criterion or rock alteration index (RAI) are discussed in light of the present, more general theory. In order to validate the finite element simulation, we derive an analytical solution for the rock alteration index of a benchmark problem on a two-dimensional rectangular domain. Since the geometry and boundary conditions of the benchmark problem can be easily and exactly modelled, the analytical solution is also useful for validating other numerical methods, such as the finite difference method and the boundary element method, when they are used to dear with this kind of problem. Finally, the potential of the theory is illustrated by means of finite element studies related to coupled flow problems in materially homogeneous and inhomogeneous porous rock masses. (C) 1998 Elsevier Science S.A. All rights reserved.

Modelling dispersal behaviour on a fractal landscape

We use a spatially explicit population model to explore the population consequences of different habitat selection mechanisms on landscapes with fractal variation in habitat quality. We consider dispersal strategies ranging from random walks to perfect habitat selectors for two species of arboreal marsupial, the greater glider (Petauroides volans) and the mountain brushtail possum (Trichosurus caninus). In this model increasing habitat selection means individuals obtain higher quality territories, but experience increased mortality during dispersal. The net effect is that population sizes are smaller when individuals actively select habitat. We find positive relationships between habitat quality and population size can occur when individuals do not use information about the entire landscape when habitat quality is spatially autocorrelated. We also find that individual behaviour can mitigate the negative effects of spatial variation on population average survival and fecundity. (C) 1998 Elsevier Science Ltd. All rights reserved.

Experimental study of phase equilibria in the system PbO-ZnO-SiO2

An experimental study on the ternary system PbO-ZnO-SiO2, in air by high-temperature equilibration and quenching techniques followed by electron probe X-ray microanalysis was carried out as part of the wider research program on the six-component system PbO-ZnO-SiO2-CaO-FeO-Fe2O3, which combines experimental and thermodynamic computer modeling techniques to characterize zinc and lead industrial slags. Liquidus and solidus data were reported for all primary phase fields in the system PbO-ZnO-SiO2 in the temperature range 640 degrees C to 1400 degrees C (913 to 1673 K).

Observation of a cubic phase in mesoporous silicate films grown at the air/water interface

X-Ray diffraction is reported from mesoporous silicate films grown at the air/water interface. The films were studied both as powdered films, and oriented on silicon or mica sheets. At early stages of growth we observe Bragg diffraction from a highly ordered cubic phase, with both long and short d-spacing peaks. We have assigned this as a discontinuous micellar Pm3n phase in which the silica is partly ordered. Later films retain only the known hexagonal p6m peaks and have lost any order both at short d-spacings and the longer d-spacing Bragg peaks characteristic of the cubic structure. The silica framework is considerably expanded from that in bulk amorphous silica, average Si Si distances are some 30% greater. Incorporation of glycerol or polyethylene glycol preserves the earlier cubic structure. To be consistent with earlier, in situ, X-ray and neutron reflectivity data we infer that both structures are produced after a phase transition from a less-ordered him structure late in the induction phase. The structural relations between the film Pm3n and p6m phase(s) and the known bulk SBA-1 and MCM-41 phases are briefly discussed.

Modelling the large deformations in stratified media - the Cosserat continuum approach

Methods employing continuum approximation in describing the deformation of layered materials possess a clear advantage over explicit models, However, the conventional implicit models based on the theory of anisotropic continua suffers from certain difficulties associated with interface slip and internal instabilities. These difficulties can be remedied by considering the bending stiffness of the layers. This implies the introduction of moment (couple) stresses and internal rotations, which leads to a Cosserat-type theory. In the present model, the behaviour of the layered material is assumed to be linearly elastic; the interfaces are assumed to be elastic perfectly plastic. Conditions of slip or no slip at the interfaces are detected by a Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformation analysis. The model is incorporated into the finite element program AFENA and validated against analytical solutions of elementary buckling problems in layered medium. A problem associated with buckling of the roof and the floor of a rectangular excavation in jointed rock mass under high horizontal in situ stresses is considered as the main application of the theory. Copyright (C) 1999 John Wiley & Sons, Ltd.

Phase diagram for a class of spin-1/2 Heisenberg models interpolating between the square-lattice, the triangular-lattice, and the linear-chain limits

We study the spin-1/2 Heisenberg models on an anisotropic two-dimensional lattice which interpolates between the square lattice at one end, a set of decoupled spin chains on the other end, and the triangular-lattice Heisenberg model in between. By series expansions around two different dimer ground states and around various commensurate and incommensurate magnetically ordered states, we establish the phase diagram for this model of a frustrated antiferromagnet. We find a particularly rich phase diagram due to the interplay of magnetic frustration, quantum fluctuations, and varying dimensionality. There is a large region of the usual two-sublattice Neel phase, a three-sublattice phase for the triangular-lattice model, a region of incommensurate magnetic order around the triangular-lattice model, and regions in parameter space where there is no magnetic order. We find that the incommensurate ordering wave vector is in general altered from its classical value by quantum fluctuations. The regime of weakly coupled chains is particularly interesting and appears to be nearly critical. [S0163-1829(99)10421-1].

Phase diagram of a Heisenberg spin-Peierls model with quantum phonons

Using a new version of the density-matrix renormalization group we determine the phase diagram of a model of an antiferromagnetic Heisenberg spin chain where the spins interact with quantum phonons. A quantum phase transition from a gapless spin-fluid state to a gapped dimerized phase occurs at a nonzero value of the spin-phonon coupling. The transition is in the same universality class as that of a frustrated spin chain, to which the model maps in the diabatic limit. We argue that realistic modeling of known spin-Peierls materials should include the effects of quantum phonons.

Effect of disorder on quantum phase transitions in anisotropic XY spin chains in a transverse field

We present some exact results for the effect of disorder on the critical properties of an anisotropic XY spin chain in a transverse held. The continuum limit of the corresponding fermion model is taken and in various cases results in a Dirac equation with a random mass. Exact analytic techniques can then be used to evaluate the density of states and the localization length. In the presence of disorder the ferromagnetic-paramagnetic or Ising transition of the model is in the same universality class as the random transverse field Ising model solved by Fisher using a real-space renormalization-group decimation technique (RSRGDT). If there is only randomness in the anisotropy of the magnetic exchange then the anisotropy transition (from a ferromagnet in the x direction to a ferromagnet in the y direction) is also in this universality class. However, if there is randomness in the isotropic part of the exchange or in the transverse held then in a nonzero transverse field the anisotropy transition is destroyed by the disorder. We show that in the Griffiths' phase near the Ising transition that the ground-state energy has an essential singularity. The results obtained for the dynamical critical exponent, typical correlation length, and for the temperature dependence of the specific heat near the Ising transition agree with the results of the RSRODT and numerical work. [S0163-1829(99)07125-8].

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.

Absorption spectrum of a strongly driven atom in a detuned squeezed vacuum

We present numerical and analytical results for the Mollow probe absorption spectrum of a coherently driven two-level system in a narrow bandwidth squeezed vacuum field. The spectra are calculated for the case where the Rabi frequency of the driving field is much larger than the natural linewidth and the squeezed vacuum carrier frequency is detuned from the driving laser frequency. The driving laser is on resonance. We show that in a detuned squeezed vacuum the standard Mellow features are each split into triplets. The central components of each triplet are weakly dependent on the squeezing phase but the sidebands strongly depend on the phase and can have dispersive or absorptive/emissive profiles. We also derive approximate analytical expressions for the spectral features and find that the multi-peak structure of the spectrum can be interpreted either via the eigenfrequencies of a generalized Floquet Hamiltonian or in terms of three-photon transitions between dressed stales involving a probe field photon and a correlated photon pair from the squeezed vacuum field.

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.

Incorporating phase retardation effects into radiofrequency resonator models: application to magnetic resonance microscopy

A method is presented for including path propagation effects into models of radiofrequency resonators for use in magnetic resonance imaging. The method is based on the use of Helmholtz retarded potentials and extends our previous work on current density models of resonators based on novel inverse finite Hilbert transform solutions to the requisite integral equations. Radiofrequency phase retardation effects are most pronounced at high field strengths (frequencies) as are static field perturbations due to the magnetic materials in the resonators themselves. Both of these effects are investigated and a novel resonator structure presented for use in magnetic resonance microscopy.

Role of the crystallite phase of TiO2 in heterogeneous photocatalysis for phenol oxidation in water

A series of TiO2 samples with different anatase-to-rutile ratios was prepared by calcination, and the roles of the two crystallite phases of titanium(IV) oxide (TiO2) on the photocatalytic activity in oxidation of phenol in aqueous solution were studied. High dispersion of nanometer-sized anatase in the silica matrix and the possible bonding of Si-O-Ti in SiO2/TiO2 interface were found to stabilize the crystallite transformation from anatase to rutile. The temperature for this transformation was 1200 degrees C for the silica-titania (ST) sample, much higher than 700 degrees C for Degussa P25, a benchmarking photocatalyst. It is shown that samples with higher anatase-to-rutile ratios have higher activities for phenol degradation. However, the activity did not totally disappear after a complete crystallite transformation for P25 samples, indicating some activity of the rutile phase. Furthermore, the activity for the ST samples after calcination decreased significantly, even though the amount of anatase did not change much. The activity of the same samples with different anatase-to-rutile ratios is more related to the amount of the surface-adsorbed water and hydroxyl groups and surface area. The formation of rutile by calcination would reduce the surface-adsorbed water and hydroxyl groups and surface area, leading to the decrease in activity.