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.

Formal derivation of finite state machines for class testing

Previous work on generating state machines for the purpose of class testing has not been formally based. There has also been work on deriving state machines from formal specifications for testing non-object-oriented software. We build on this work by presenting a method for deriving a state machine for testing purposes from a formal specification of the class under test. We also show how the resulting state machine can be used as the basis for a test suite developed and executed using an existing framework for class testing. To derive the state machine, we identify the states and possible interactions of the operations of the class under test. The Test Template Framework is used to formally derive the states from the Object-Z specification of the class under test. The transitions of the finite state machine are calculated from the derived states and the class's operations. The formally derived finite state machine is transformed to a ClassBench testgraph, which is used as input to the ClassBench framework to test a C++ implementation of the class. The method is illustrated using a simple bounded queue example.

Effect of temperature and a crosslinking promoter on the gamma-radiolysis of a perfluoroelastomer

Methods of promoting the radiation-induced cross linking of poly(tetrafluoro-ethylene-co-perfluoromethyl vinyl ether) (TFE/PMVE) have been investigated. Greater control of the crosslinking and chain-scission reactions was achieved by varying the radiolysis temperature. This was attributed to temperature affecting the mobilities of reactive species such as polymeric free radicals. These reactive species are precursors to radiation-induced cross links and chain-ends. Analysis of the sol/gel behaviour, tensile properties and FTIR indicated that the optimum temperature for the radiation crosslinking of TFE/PMVE, at a dose of 150 kGy, was 263 K. This temperature was 10 K below the glass transition temperature. Incorporation of 1 wt% triallyl isocyanurate (TAIC) greatly amplified the radiation crosslinking of TFE/PMVE, The dose for gelation was decreased by 70%, and the additive imparted superior mechanical properties compared to the neat irradiated TFE/PMVE. Electron spin resonance (ESR) measurements showed higher radical yields at 77 K with the 1 wt% TAIC, indicating that the crosslinking promoter was acting as a radical trap. (C) 1999 Society of Chemical Industry.

Adaptation of sorghum: characterisation of genotypic flowering responses to temperature and photoperiod

Sorghum [Sorghum bicolor (L.) Moench] is an important cereal crop grown in a wide range of tropical and temperate environments. This study was conducted to characterise the photothermal flowering responses of sorghum genotypes and to examine relationships between photothermal characteristics and environment of origin in order to better understand the phenological basis of adaptation to environment in sorghum. Twenty-four germplasm accessions and one hybrid from 24 major sorghum-growing areas were grown in a wide range of environments varying in temperature and photoperiod in India, Kenya and Mall between 1992 and 1995. Times from sowing to flowering (f) were recorded, and the responsiveness of 1/f to temperature and photoperiod was quantified using photothermal models. Times from sowing to flowering were accurately predicted in a wide range of environments using a multiplicative rate photothermal model. Significant variation in the minimum time to flower (F-m) and photoperiod sensitivity (critical photoperiod, P-c, and photoperiod-sensitivity slope, P-s) was observed among the genotypes; in contrast there was little variation in base temperature (Tb) Adaptation of sorghum to the diverse environments in which it is grown was largely determined by photoperiod sensitivity and minimum time to flower; photoperiod sensitivity determines bread adaptation to latitude (daylength), while variation in the minimum time to flower determines specific adaptation within smaller ranges of latitude, e.g. within the humid and sub-humid tropics.

Variable temperature and pressure study of the aquation reactions of cobalt(III) and chromium(III) penta- and tetra-amines

Preparation of a series of specific penta- and tetra-amine derivatives of Co-III and Cr-III with a neutral leaving ligand has been carried out in order to accomplish a fine tuning of the associativeness/dissociativeness of their substitution reactions. Spontaneous aquation reactions of the neutral ligands have been studied at variable temperature and pressure. Although rate constants and thermal activation parameters show an important degree of scatter, the values determined for the activation volumes of the substitution process illustrate the mechanistic fine tuning that may be achieved for these reactions. In all cases, in the absence of important steric constraints in the molecule, electronic inductive effects seem to be the most important factor accounting for the dissociative shifts observed both for pentaamine (i.e.Delta V double dagger=+4.0 or +14.0 cm(3) mol(-1) and +5.2 or +16.5 cm(3) mol(-1) for the aquation of cis- or trans-[Co(MeNH2)(NH3)(4)(DMF)](3+) and cis- or trans-[CoL15(DMF)](3+) respectively, where L-15 represents a pentaamine macrocyclic ligand), and tetraamine systems (i.e.Delta V double dagger=+4.1 or +8.4 cm(3) mol(-1) and -10.8 or -7.4 cm(3) mol(-1) for the aquation of cis-[Co(NH3)(4)Cl(DMAC)](2+) (DMAC=dimethylacetamide) or cis-[Co(en)(2)Cl(DMAC)](2+) and cis-[Cr(NH3)(4)Cl(DMF)](2+) or cis -[Cr(en)(2)Cl(DMF)](2+)). From the results, clear evidence is obtained which indicates that, only when the situation is borderline I-a/I-d, or the steric demands are increased dramatically, dissociative shifts are observed; in all other cases electronic inductive effects seem to be dominant for such a tuning of the substitution process.

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.

Long range forecasts of the numbers of Helicoverpa punctigera and H-armigera (Lepidoptera : Noctuidae) in Australia using the Southern Oscillation Index and the Sea Surface Temperature

The use of long-term forecasts of pest pressure is central to better pest management. We relate the Southern Oscillation Index (SOI) and the Sea Surface Temperature (SST) to long-term light-trap catches of the two key moth pests of Australian agriculture, Helicoverpa punctigera (Wallengren) and H. armigera (Hubner), at Narrabri, New South Wales over 11 years, and for H. punctigera only at Turretfield, South Australia over 22 years. At Narrabri, the size of the first spring generation of both species was significantly correlated with the SOI in certain months, sometimes up to 15 months before the date of trapping. Differences in the SOI and SST between significant months were used to build composite variables in multiple regressions which gave fitted values of the trap catches to less than 25% of the observed values. The regressions suggested that useful forecasts of both species could be made 6-15 months ahead. The influence of the two weather variables on trap catches of H. punctigera at Turretfield were not as strong as at Narrabri, probably because the SOI was not as strongly related to rainfall in southern Australia as it is in eastern Australia. The best fits were again given by multiple regressions with SOI plus SST variables, to within 40% of the observed values. The reliability of both variables as predictors of moth numbers may be limited by the lack of stability in the SOI-rainfall correlation over the historical record. As no other data set is available to test the regressions, they can only be tested by future use. The use of long-term forecasts in pest management is discussed, and preliminary analyses of other long sets of insect numbers suggest that the Southern Oscillation Index may be a useful predictor of insect numbers in other parts of the world.

Soil carbon determination by high temperature combustion - A comparison with dichromate oxidation procedures and the influence of charcoal and carbonate carbon on the measured value

The measurement of organic carbon in soils has traditionally used dichromate oxidation procedures including the Wakley and Black and the Heanes methods. The measurement of carbon in soils by high temperature combustion is now widely used providing a rapid automated procedure without the use of toxic chemicals. This procedure however measures total carbon thus requiring some means of correction for soil samples containing carbonate and charcoal forms of carbon. This paper examines the effects of known additions of charcoal to a range of soil types on the results obtained by the Walkley and Black, Heanes and combustion methods. The results show, that while the charcoal carbon does not react under Walkley and Black conditions, some proportion does so with the Heanes method. A comparison of six Australian Soil and Plant Analysis Council reference soil samples by the three methods showed good agreement between the Heanes method, the combustion method and only slightly lower recoveries by the Walkley and Black procedure. Carbonate carbon will cause an overestimation of soil organic carbon by the combustion method thus requiring a separate determination of carbonate carbon to be applied as a correction. This work shows that a suitable acid pre-treatment of alkaline soils in the sample boats followed by a drying step eliminates the carbonate carbon prior to combustion and the need for an additional measurement. The measurement of carbon in soils by high temperature combustion in an oxygen atmosphere has been shown to be a rapid and reliable method capable of producing results in good agreement with one of the established dichromate oxidation procedures.

Prediction of rock alteration patterns: a potential tool in mineral exploration

Hydrothermal alteration of a quartz-K-feldspar rock is simulated numerically by coupling fluid flow and chemical reactions. Introduction of CO2 gas generates an acidic fluid and produces secondary quartz, muscovite and/or pyrophyllite at constant temperature and pressure of 300 degrees C and 200 MPa. The precipitation and/or dissolution of the secondary minerals is controlled by either mass-action relations or rate laws. In our simulations the mass of the primary elements are conserved and the mass-balance equations are solved sequentially using an implicit scheme in a finite-element code. The pore-fluid velocity is assumed to be constant. The change of rock volume due to the dissolution or precipitation of the minerals, which is directly related to their molar volume, is taken into account. Feedback into the rock porosity and the reaction rates is included in the model. The model produces zones of pyrophyllite quartz and muscovite due to the dissolution of K-feldspar. Our model simulates, in a simplified way, the acid-induced alteration assemblages observed in various guises in many significant mineral deposits. The particular aluminosilicate minerals produced in these experiments are associated with the gold deposits of the Witwatersrand Basin.

Resistance to temperature extremes between and within life cycle stages in Drosophila serrata, D-birchii and their hybrids: intraspecific and interspecific comparisons

Intraspecific Drosophila studies suggest that resistance to heal and cold stresses are largely independent and that correlations across life cycle stages are low whereas comparisons of Drosophila species indicate correlations between heat and cold resistance as well as between resistance levels in different life cycle stages. These inconsistent results may reflect differences in associations among traits at the interspecific and intraspecific levels or interspecific correlations arising because of correlated selection pressures. These alternatives were tested using Drosophila serrata, D. birchii and hybrids derived from these species. Variation among hybrid lines and families was used to test associations at the interspecific level while intraspecific variation was examined using isofemale lines of D. serrata. There was a significant association between adult heat knockdown time at 38 degreesC and adult cold resistance in one set of hybrid lines. An association between female knockdown resistance to heat and larval heat resistance was also evident in one set of hybrids. Resistance to heat anti cold at the larval stage were not correlated at either the intraspecific or interspecific levels. At the intraspecific level, lan al heat resistance and two measures of adult heat resistance were uncorrelated. Moreover, adult and larval cold resistance measures were not correlated at either the intraspecific or interspecific levels. These results suggest that there are no associations between resistance to heat and cold extremes and that extreme temperature resistance is largely independent across life cycle stages at both the intraspecific and interspecific levels. Species associations may therefore arise from correlated selection pressures rather than trait correlations. (C) 2000 The Linnean Society of London.

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.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

Semisolid casting of AlSi7Mg0.35 alloy produced by low-temperature pouring

AlSi7Mg0.35 alloy was cast into permanent moulds using different pouring temperatures (725 to 625degreesC). As the pouring temperature decreased, the as-cast microstructure changed from a coarse dendritic structure, through fine equiaxed grains to fine rosette-like grains. The as-cast materials were then partially remelted and isothermally held at 580degreesC prior to semisolid casting into a stepped die. The feedstock material cast from a high temperature filled only half the die, with severe segregation and other defects. The low-temperature-poured material completely filled the die with negligible porosity. The quality of semisolid castings is significantly affected by the microstructure of the semisolid feedstock material that arises from a combination of as-cast and subsequent thermal treatment conditions. The paper describes (a) the influence of pouring temperature on the microstructure of feedstock; (b) microstructure evolution through remelting and (c) the quality of semisolid castings produced with this material. For A17Si0.35Mg alloy, low temperature pouring in the range of 625-650degreesC followed by suitable isothermal holding treatment can result in good quality semisolid casting.