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.

Analysis of steady-state heat transfer through mid-crustal vertical cracks with upward throughflow in hydrothermal systems

We conduct a theoretical analysis of steady-state heat transfer problems through mid-crustal vertical cracks with upward throughflow in hydrothermal systems. In particular, we derive analytical solutions for both the far field and near field of the system. In order to investigate the contribution of the forced advection to the total temperature of the system, two concepts, namely the critical Peclet number and the critical permeability of the system, have been presented and discussed in this paper. The analytical solution for the far field of the system indicates that if the pore-fluid pressure gradient in the crust is lithostatic, the critical permeability of the system can be used to determine whether or not the contribution of the forced advection to the total temperature of the system is negligible. Otherwise, the critical Peclet number should be used. For a crust of moderate thickness, the critical permeability is of the order of magnitude of 10(-20) m(2), under which heat conduction is the overwhelming mechanism to transfer heat energy, even though the pore-fluid pressure gradient in the crust is lithostatic. Furthermore, the lower bound analytical solution for the near field of the system demonstrates that the permeable vertical cracks in the middle crust can efficiently transfer heat energy from the lower crust to the upper crust of the Earth. Copyright (C) 2002 John Wiley Sons, Ltd.

Scope for further analytical solutions for constant flux infiltration into a semi-infinite soil profile or redistribution in a finite soil profile

[1] We attempt to generate new solutions for the moisture content form of the one-dimensional Richards' [1931] equation using the Lisle [1992] equivalence mapping. This mapping is used as no more general set of transformations exists for mapping the one-dimensional Richards' equation into itself. Starting from a given solution, the mapping has the potential to generate an infinite number of new solutions for a series of nonlinear diffusivity and hydraulic conductivity functions. We first seek new analytical solutions satisfying Richards' equation subject to a constant flux surface boundary condition for a semi-infinite dry soil, starting with the Burgers model. The first iteration produces an existing solution, while subsequent iterations are shown to endlessly reproduce this same solution. Next, we briefly consider the problem of redistribution in a finite-length soil. In this case, Lisle's equivalence mapping is generalized to account for arbitrary initial conditions. As was the case for infiltration, however, it is found that new analytical solutions are not generated using the equivalence mapping, although existing solutions are recovered.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

Finite-difference time-domain-based studies of MRI pulsed field gradient-induced eddy currents inside the human body

In modern magnetic resonance imaging (MRI), patients are exposed to strong, rapidly switching magnetic gradient fields that, in extreme cases, may be able to elicit nerve stimulation. This paper presents theoretical investigations into the spatial distribution of induced current inside human tissues caused by pulsed z-gradient fields. A variety of gradient waveforms have been studied. The simulations are based on a new, high-definition, finite-difference time-domain method and a realistic inhomogeneous 10-mm resolution human body model with appropriate tissue parameters. it was found that the eddy current densities are affected not only by the pulse sequences but by many parameters such as the position of the body inside the gradient set, the local biological material properties and the geometry of the body. The discussion contains a comparison of these results with previous results found in the literature. This study and the new methods presented herein will help to further investigate the biological effects caused by the switched gradient fields in a MRI scan. (C) 2002 Wiley Periodicals, Inc.

Conductivity, NMR and crystallographic study of N,N,N,N-tetramethylammonium dicyanamide plastic crystal phases: an archetypal ambient temperature plastic electrolyte material

N,N,N,N-Tetramethylammonium dicyanamide (Me(4)NDCA) has been examined via differential scanning calorimetry (DSC), thermogravimetric analysis, conductivity, single crystal X-ray diffraction and H-1 nuclear magnetic resonance (NMR) analyses, and was found to be highly conductive in the solid state (sigma = 10(-3) S cm(-2) at 420 K) and to also exhibit unusual plastic crystal behaviour. To investigate the correlation between such behaviour and the occurrence of molecular rotations in the crystal, H-1 NMR second moment measurements are compared with calculated values predicted from the crystal structure. While DSC analysis indicates a number of solid-solid transitions at ambient temperatures, subsequent H-1 NMR analysis of the Me4N+ cation shows that a variety of rotational motions become active at low (

Maturation temperature and rainfall influence seed dormancy characteristics of annual ryegrass (Lolium rigidum)

The role of temperature and rainfall during seed development in modulating subsequent seed dormancy status was studied for Lolium rigidum Gaud. (annual ryegrass). Climatic parameters relating to geographic origin were compared with annual ryegrass seed dormancy characteristics for seeds collected from 12 sites across the southern Western Australian cropping region. Seed germination was tested soon after collection and periodically during subsequent after-ripening. Temperature in the year of seed development and long-term rainfall patterns showed correlations with aspects of seed dormancy, particularly the proportion of seeds remaining dormant following 5 months of after-ripening. Consequently, for one population the temperature (warm/cool) and water supply (adequate/reduced) during seed development were manipulated to investigate the role of maternal environment in the quantity and dormancy characteristics of seeds produced. Seeds from plants grown at warm temperatures were fewer in number, weighed less, and were less dormant than those from plants grown at cool temperature. Seeds that developed under both cool temperature and reduced moisture conditions lost dormancy faster than seeds from well-watered plants. Seed maturation environment, particularly temperature, can have a significant effect on annual ryegrass seed numbers and seed dormancy characteristics.