Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

The Heisenberg antiferromagnet on an anisotropic triangular lattice: linear spin-wave theory

We consider the effect of quantum spin fluctuations on the ground-state properties of the Heisenberg antiferromagnet on an anisotropic triangular lattice using linear spin-wave (LSW) theory. This model should describe the magnetic properties of the insulating phase of the kappa-(BEDT-TTF)(2)X family of superconducting molecular crystals. The ground-state energy, the staggered magnetization, magnon excitation spectra, and spin-wave velocities are computed as functions of the ratio of the antiferromagnetic exchange between the second and first neighbours, J(2)/J(1). We find that near J(2)/J(1) = 0.5, i.e., in the region where the classical spin configuration changes from a Neel-ordered phase to a spiral phase, the staggered magnetization vanishes, suggesting the possibility of a quantum disordered state. in this region, the quantum correction to the magnetization is large but finite. This is in contrast to the case for the frustrated Heisenberg model on a square lattice, for which the quantum correction diverges logarithmically at the transition from the Neel to the collinear phase. For large J(2)/J(1), the model becomes a set of chains with frustrated interchain coupling. For J(2) > 4J(1), the quantum correction to the magnetization, within LSW theory, becomes comparable to the classical magnetization, suggesting the possibility of a quantum disordered state. We show that, in this regime, the quantum fluctuations are much larger than for a set of weakly coupled chains with non-frustrated interchain coupling.

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 pore-size-dependent equation of state for multilayer adsorption in cylindrical mesopores

The chemical potential of adsorbed film inside cylindrical mesopores is dependent on the attractive interactions between the adsorbed molecules and adsorbent, the curvature of gas/adsorbed phase interface, and surface tension. A state equation of the adsorbed film is proposed to take into account the above factors. Nitrogen adsorption on model adsorbents, MCM-41, which exhibit uniform cylindrical channels, are used to verify the theoretical analysis. The proposed theory is capable of describing the important features of adsorption processes in cylindrical mesopores. According to this theory, at a given relative pressure, the smaller the pore radius is, the thicker the adsorbed film will be. The thickening of adsorbed films in the pores as the vapor pressure increases inevitably causes an increase in the interface curvature, which consequently leads to capillary condensation. Besides, this study confirmed that the interface tension depends substantially on the interface curvature in small mesopores. A quantitative relationship between the condensation pressure and the pore radius can be derived from the state equation and used to predict the pore radius from a condensation pressure, or vice versa.

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

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.

Approach to designing rotating drum bioreactors for solid-state fermentation on the basis of dimensionless design factors

The development of large-scale solid-stale fermentation (SSF) processes is hampered by the lack of simple tools for the design of SSF bioreactors. The use of semifundamental mathematical models to design and operate SSF bioreactors can be complex. In this work, dimensionless design factors are used to predict the effects of scale and of operational variables on the performance of rotating drum bioreactors. The dimensionless design factor (DDF) is a ratio of the rate of heat generation to the rate of heat removal at the time of peak heat production. It can be used to predict maximum temperatures reached within the substrate bed for given operational variables. Alternatively, given the maximum temperature that can be tolerated during the fermentation, it can be used to explore the combinations of operating variables that prevent that temperature from being exceeded. Comparison of the predictions of the DDF approach with literature data for operation of rotating drums suggests that the DDF is a useful tool. The DDF approach was used to explore the consequences of three scale-up strategies on the required air flow rates and maximum temperatures achieved in the substrate bed as the bioreactor size was increased on the basis of geometric similarity. The first of these strategies was to maintain the superficial flow rate of the process air through the drum constant. The second was to maintain the ratio of volumes of air per volume of bioreactor constant. The third strategy was to adjust the air flow rate with increase in scale in such a manner as to maintain constant the maximum temperature attained in the substrate bed during the fermentation. (C) 2000 John Wiley & Sons, Inc.

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.

Further results on the relationship between mu-invariant measures and quasi-stationary distributions for absorbing continuous-time Markov chains

This note considers continuous-time Markov chains whose state space consists of an irreducible class, C, and an absorbing state which is accessible from C. The purpose is to provide results on mu-invariant and mu-subinvariant measures where absorption occurs with probability less than one. In particular, the well-known premise that the mu-invariant measure, m, for the transition rates be finite is replaced by the more natural premise that m be finite with respect to the absorption probabilities. The relationship between mu-invariant measures and quasi-stationary distributions is discussed. (C) 2000 Elsevier Science Ltd. All rights reserved.

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.

Mesoionic pyridopyrimidinylium and pyridooxazinylium olates and non-mesoionic pyridopyrimidinones. Structures in the solid state, solution, and matrices

X-Ray crystal structures, C-13 NMR spectra and theoretical calculations (B3LYP/6-31G*) are reported for the mesoionic (zwitterionic) pyridopyrimidinylium- and pyridooxazinyliumolates 2a, 3a and 5a,b as well as the enol ether 11b and the enamine 11c. The 1-NH compounds like 1a, 2a and 3a exist in the mesoionic form in the crystal and in solution, but the OH tautomers such as 1b and 2b dominate in the gas phase as revealed by the Ar matrix IR spectra in conjunction with DFT calculations. All data indicate that the mesoionic compounds can be regarded as intramolecular pyridine-ketene zwitterions (cf. 16 --> 17) with a high degree of positive charge on the pyridinium nitrogen, a long pyridinium N-CO bond (ca. 1.44-1.49 Angstrom), and normal C=O double bonds (ca. 1.22 Angstrom). All mesoionic compounds exhibit a pronounced tilting of the olate C=O groups (the C=O groups formally derived from a ketene) towards the pyridinium nitrogen, giving NCO angles of 110-118 degrees. Calculations reveal a hydrogen bond with 6-CH, analogous to what is found in ketene-pyridine zwitterions and the C3O2-pyridine complex. The 2-OH tautomers of type 1b, 2b, and 11 also show a high degree of zwitterionic character as indicated by the canonical structures 11 12.

Signatures of the pair-coherent state

We explore in detail the possibility of generating a pair-coherent state in the nondegenerate parametric oscillator when decoherence is included. Such states are predicted in the transient regime in parametric oscillation where the pump mode is adiabatically eliminated. Two specific signatures are examined to indicate whether the state of interest has been generated, the Schrodinger cat state-like signatures, and the fidelity. Solutions in a transient regime reveal interference fringes which are indicative of the formation of a Schrodinger cat state. The fidelity indicates the purity of our prepared state compared with the ideal pair-coherent state.

An investigation of pretend play, mental state terms and false belief understanding: In search of a metarepresentational link

The aim of this study was to investigate the association between false belief comprehension, the exhibition of pretend play and the use of mental state terms in pre-school children. Ferry children, aged between 36 and 54 months were videotaped engaging in free play with each parent. The exhibit-ion of six distinct acts of pretend play and the expression of 16 mental sr:ate terms were coded during play. Each child was also administered a pantomime task and three standard false belief casks. Reliable associations were also found between false belief performance and the pretence categories of object substitution and role assignment, and the exhibition of imaginary object pantomimes. Moreover, the use of mental state terms was positively correlated with false belief and the pretence categories of object substitution, imaginary play and role assignment, and negatively correlated with the exhibition of body part object pantomimes. These findings indicate that the development of a mental state lexicon and some, bur not all, components of pretend play are dependent on the capacity for metarepresentational cognition.

Entangled coherent-state qubits in an ion trap

We show how entangled qubits can be encoded as entangled coherent states of two-dimensional center-of-mass vibrational motion for two ions in an ion trap. The entangled qubit state is equivalent to the canonical Bell state, and we introduce a proposal for entanglement transfer from the two vibrational modes to the electronic states of the two ions in order for the Bell state to be detected by resonance fluorescence shelving methods.