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.

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.

Size of the first spring generation of Helicoverpa punctigera (Wallengren) (Lepidoptera : Noctuidae) and winter rain in central Australia

Serious infestations of Helicoverpa punctigera are experienced yearly in the eastern cropping regions of Australia. Regression analysis was used to determine whether the size of the first generation in spring (G(1)), which is comprised mostly of immigrants from inland Australia, was related to monthly rainfall in inland winter breeding areas. Data from two long series of light-trap catches at Narrabri in New South Wales (NSW) and Turretfield in South Australia (SA) were used in the analyses. The size of G1 at Narrabri in each year was significantly regressed on the amount of rainfall in western Queensland and NSW in May and June. The size of G1 at Turretfield each year was significantly regressed on the amount of rain in May, June and July in western Queensland and NSW and also in the desert of central Western Australia. Low r(2) values of the regressions suggest that rainfall data for more sites, as well as biological and other physical factors, such as temperature, evaporation, and prevailing wind systems, may need to be included to improve forecasts of the potential magnitude of the infestations in coastal cropping regions.

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 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.

A novel method for tailoring the pore-opening size of MCM-41 materials

A novel pore tailoring method is proposed by which the pore-opening sizes of MCM-41 materials can be finely tuned without significant loss in pore volume and surface area.

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.

Life-history consequences of divergent selection on egg size in Drosophila melanogaster

Life histories are generally assumed to evolve via antagonistic pleiotropy (negative genetic correlations) among traits, and trade-offs between life-history traits are typically studied using either phenotypic manipulations or selection experiments. We investigated the trade-off between egg size and fecundity in Drosophila melanogaster by examining both the phenotypic and genetic relationships between these traits after artificial selection for large and small eggs, relative to female body size. Egg size responded strongly to selection in both directions, increasing in the large-egg selected lines and decreasing in the small-egg selected lines. Phenotypic correlations between egg size and fecundity in the large-egg selected lines were negative, but no relationship between these traits occurred in either the control or small-egg selected lines. There was no negative genetic correlation between egg size and fecundity. Total reproductive allocation decreased in the small-egg selected lines but did not increase in the large-egg lines. Our results have three implications. First, our selection procedure may have forced females selected for large eggs into a physiological trade-off not reflected in a negative genetic correlation between these traits. Second, the lack of a negative genetic correlation between egg size and number suggests that the phenotypic trade-off frequently observed between egg size and number in other organisms may not evolve over the short term via a direct genetic trade-off whereby increases in egg size are automatically accompanied by decreased fecundity. Finally, total reproductive allocation may not evolve independently of egg size as commonly assumed.

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.

Intraspecific co-variation between egg and body size affects fertilisation kinetics of free-spawning marine invertebrates

Fertilisation of eggs of free-spawning marine invertebrates depends on factors affecting sperm concentration in the field and also on gamete characteristics such as egg size. In the free-spawning intertidal ascidian Pyura stolonifera mean egg size increased with maternal size in 2 separate populations. The largest ascidian produced eggs that were, on average, 50% greater in volume than the eggs produced by the smallest individual studied. There was no evidence to suggest that egg density varied with adult size and egg dry organic weight increased with maternal size. The fertilisation kinetics of this species were strongly affected by the variation in egg size, with the eggs of large individuals requiring much less concentrated sperm to achieve maximal levels of fertilisation success than the eggs of small individuals. We suggest that variation in egg size between individuals of different sizes and ages may be an important factor in determining fertilisation success for ascidians of this species.

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.

Muscle fibre types and size distribution in sub-antarctic notothenioid fishes

The presumptive tonic muscles fibres of Cottoperca gobio, Champsocephalus esox, Harpagifer bispinis, Eleginops maclovinus, Patagonothen tessellata, P. cornucola and Paranotothenia magellanica stained weakly or were unstained for glycogen, lipid, succinic dehydrogenase (SDHase) and myosin ATPase (mATPase) activity. Slow, intermediate and fast twitch muscle fibres, distinguished on the basis of the pH stability of their mATPases, showed intense, moderate and low staining activity for SDHase, respectively. Slow fibres were the major component of the pectoral fin adductor profundis muscle. The proportion of different muscle fibre types varied from the proximal to distal end of the muscle, but showed relatively little variation between species. The myotomes contained a lateral superficial strip of red muscle composed of presumptive tonic, slow twitch and intermediate fibres, thickening to a major wedge at the horizontal septum. All species also had characteristic secondary dorsal and ventral wedges of red muscle. The relative abundance and localization of muscle fibre types in the red muscle varied between species and with body size in the protandric hermaphrodite E. maclovinus. The frequency distribution of diameters for fast twitch muscle fibres, the major component of deep white muscle, was determined in fish of a range of body sizes. The absence of fibres <20 mu m diameter was used as a criterion for the cessation of muscle fibre recruitment. Fibre recruitment had stopped in P, tessellata of 13.8 cm L-T and E, maclovinus of 32.8 cm L-T, equivalent to 49 and 36.5% of their recorded maximum sizes respectively. As a result in 20-cm P. tessellata, the maximum fibre diameter was 300 mu m and 36% of fibres were in excess of 200 mu m The unusually large maximum fibre diameter, the general arrangement of the red muscle layer and the extreme pH lability of the mATPase of fast twitch fibres are all common characters of the sub-Antarctic and Antarctic Notothenioids, including Cottoperca gobio, the suggested sister group to the Notothenidae. (C) 2000 The Fisheries Society of the British Isles.

The line structure of balanced ternary designs with block size three, index three and rho(2) = 1, 2

In this paper necessary and sufficient conditions for a vector to be the fine structure of a balanced ternary design with block size 3, index 3 and rho(2) = 1 and 2 are determined with one unresolved case.

Characterization of pore size distributions of mesoporous materials from adsorption isotherms

In this article, a new hybrid model for estimating the pore size distribution of micro- and mesoporous materials is developed, and tested with the adsorption data of nitrogen, oxygen, and argon on ordered mesoporous materials reported in the literature. For the micropore region, the model uses the Dubinin-Rudushkevich (DR) isotherm with the Chen-Yang modification. A recent isotherm model of the authors for nonporous materials, which uses a continuum-mechanical model for the multilayer region and the Unilan model for the submonolayer region, has been extended for adsorption in mesopores. The experimental data is inverted using regularization to obtain the pore size distribution. The present model was found to be successful in predicting the pore size distribution of pure as well as binary physical mixtures of MCM-41 synthesized with different templates, with results in agreement with those from the XRD method and nonlocal density functional theory. It was found that various other recent methods, as well as the classical Broekhoff and de Beer method, underpredict the pore diameter of MCM-41. The present model has been successfully applied to MCM-48, SBA's, CMK, KIT, HMS, FSM, MTS, mesoporous fly ash, and a large number of other regular mesoporous materials.

Osmolarity effects on observed insect cell size after baculovirus infection are avoided using growth medium for sample dilution

Rates of cell size increase are an important measure of success during the baculovirus infection process. Batch and fed batch cultures sustain large fluctuations in osmolarity that can affect the measured cell volume if this parameter is not considered during the sizing protocol. Where osmolarity differences between the sizing diluent and the culture broth exist, biased measurements of size are obtained as a result of the cell osmometer response. Spodoptera frugiperda (Sf9) cells are highly sensitive to volume change when subjected to a change in osmolarity. Use of the modified protocol with culture supernatants for sample dilution prior to sizing removed the observed error during measurement.