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.

m144, a murine cytomegalovirus (MCMV)-encoded major histocompatibility complex class I homologue, confers tumor resistance to natural killer cell-mediated rejection

Until now, it has been unclear whether murine cytomegalovirus (MCMV)-encoded protein m144 directly regulates natural killer (NK) cell effector function and whether the effects of m144 are only strictly evident in the context of MCMV infection. We have generated clones of the transporter associated with antigen processing (TAP)-2-deficient RMA-S T lymphoma cell line and its parent cell line, RMA, that stably express significant and equivalent levels of m144. In vivo NK cell-mediated rejection of RMA-S-m144 lymphomas was reduced compared with rejection of parental or mock-transfected RMA-S clones, indicating the ability of m144 to regulate NK cell-mediated responses in vivo. Significantly, the accumulation of NK cells in the peritoneum was reduced in mice challenged with RMA-S-m144, as was the lytic activity of NK cells recovered from the peritoneum. Expression of m144 on RMA-S cells also conferred resistance to cytotoxicity mediated in vitro by interleukin 2-activated adherent spleen NK cells. In summary, the data demonstrate that m144 confers some protection from NK cell effector function mediated in the absence of target cell class I expression, but that in vivo the major effect of m144 is to regulate NK cell accumulation and activation at the site of immune challenge.

Engineered detoxification confers resistance against a pathogenic bacterium

We generated transgenic sugarcane plants that express an albicidin detoxifying gene (albD), which was cloned from a bacterium that provides biocontrol against leaf scald disease. Plants with albicidin detoxification capacity equivalent to 1-10 ng of AlbD enzyme per mg of leaf protein did not develop chlorotic disease symptoms in inoculated leaves, whereas all untransformed control plants developed severe symptoms. Transgenic lines with high AlbD activity in young stems were also protected against systemic multiplication of the pathogen, which is the precursor to economic disease. We have shown that genetic modification to express a toxin-resistance gene can confer resistance to both disease symptoms and multiplication of a toxigenic pathogen in its host.

Monooxygenases play only a minor role in resistance to synthetic pyrethroids in the cattle tick, Boophilus microplus

We investigated the role of monooxygenases in resistance to synthetic pyrethroids (SPs) in the cattle tick, Boophilus microplus. We found that monooxygenases play only a minor role in resistance to SPs in both resistant and susceptible strains of B. microplus. We blocked the monooxygenases with piperonyl butoxide (PBO) and simultaneously applied the SPs, flumethrin and cypermethrin to larval B. microplus. PBO increased the effect of flumethrin (synergism ratios 2.7-8.9) more than it increased the effect of cypermethrin (synergism ratios 1.9-3.1). Of the four strains tested, Parkhurst, which is resistant to SPs, was the least affected by the addition of PBO (synergism ratios after cypermethrin was applied 1.9; after flumethrin 2.7) whereas N.R.F.S., the strain susceptible to SPs, was the most affected by synergism between PBO and SPs (synergism ratio after cypermethrin was applied 3.1; after flumethrin 8.9). We hypothesize that B. microplus lacks monooxygenases capable of conferring resistance to SPs because it and its recent ancestors were blood-feeders rather than herbivores.

Sensitivity of multiple frequency bioelectrical impedance analysis to changes in ion status

Bioelectrical impedance analysis has found extensive application as a simple noninvasive method for the assessment of body fluid volumes, The measured impedance is, however, not only related to the volume of fluid but also to its inherent resistivity. The primary determinant of the resistivities of body fluids is the concentration of ions. The aim of this study was to investigate the sensitivity of bioelectrical impedance analysis to bodily ion status. Whole body impedance over a range of frequencies (4-1012 kHz) of rats was measured during infusion of various concentrations of saline into rats concomitant with measurement of total body and intracellular water by tracer dilution techniques. Extracellular resistance (R-o), intracellular resistance (R-i) and impedance at the characteristic frequency (Z(c)) were calculated. R-o and Z(c) were used to predict extracellular and total body water respectively using previously published formulae. The results showed that whilst R-o and Z(c) decreased proportionately to the amount of NaCl infused, R-i increased only slightly. Impedances at the end of infusion predicted increases iu TBW and ECW of approximately 4-6% despite a volume increase of less than 0.5% in TBW due to the volume of fluid infused. These data are discussed in relation to the assumption of constant resistivity in the prediction of fluid volumes from impedance data.

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.

Isotropic singularities in shear-free perfect fluid cosmologies

We investigate barotropic perfect fluid cosmologies which admit an isotropic singularity. From the General Vorticity Result of Scott, it is known that these cosmologies must be irrotational. In this paper we prove, using two different methods, that if we make the additional assumption that the perfect fluid is shear-free, then the fluid flow must be geodesic. This then implies that the only shear-free, barotropic, perfect fluid cosmologies which admit an isotropic singularity are the FRW models.

Insect peptides with improved protease-resistance protect mice against bacterial infection

At a time of the emergence of drug-resistant bacterial strains, the development of antimicrobial compounds with novel mechanisms of action is of considerable interest. Perhaps the most promising among these is a family of antibacterial peptides originally isolated from insects. These were shown to act in a stereospecific manner on an as-yet unidentified target bacterial protein. One of these peptides, drosocin, is inactive in vivo due to the rapid decomposition in mammalian sera. However, another family member, pyrrhocoricin, is significantly more stable, has increased in vitro efficacy against Gram-negative bacterial strains, and if administered alone, as we show here, is devoid of in vitro or in vivo toxicity. At low doses, pyrrhocoricin protected mice against Escherichia call infection, but at a higher dose augmented the infection of compromised animals. Analogs of pyrrhocoricin were, therefore, synthesized to further improve protease resistance and reduce toxicity. A linear derivative containing unnatural amino acids at both termini showed high potency and lack of toxicity in vivo and an expanded cyclic analog displayed broad activity spectrum in vitro. The bioactive conformation of native pyrrhocoricin was determined by nuclear magnetic resonance spectroscopy, and similar to drosocin, reverse turns were identified as pharmacologically important elements at the termini, bridged by an extended peptide domain. Knowledge of the primary and secondary structural requirements for in vivo activity of these peptides allows the design of novel antibacterial drug leads.

Strain-dependent fluid flow defined through rock mass classification schemes

Strain-dependent hydraulic conductivities are uniquely defined by an environmental factor, representing applied normal and shear strains, combined with intrinsic material parameters representing mass and component deformation moduli, initial conductivities, and mass structure. The components representing mass moduli and structure are defined in terms of RQD (rock quality designation) and RMR (rock mass rating) to represent the response of a whole spectrum of rock masses, varying from highly fractured (crushed) rock to intact rock. These two empirical parameters determine the hydraulic response of a fractured medium to the induced-deformations The constitutive relations are verified against available published data and applied to study one-dimensional, strain-dependent fluid flow. Analytical results indicate that both normal and shear strains exert a significant influence on the processes of fluid flow and that the magnitude of this influence is regulated by the values of RQD and RMR.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Petrographic and isotope constraints on the origin of authigenic carbonate minerals and the associated fluid evolution in Late Permian coal measures, Bowen Basin (Queensland), Australia

Authigenic carbonate minerals are ubiquitous throughout the Late Permian coal measures of the Bowen Basin, Queensland, Australia. In the northern Bowen Basin, carbonates include the following assemblages: siderite I (delta O-18(SMOW) = +11.4 to + 17%, delta C-13(PDB) = - 5.3 to + 120), Fe-Mg calcite-ankerite-siderite II mineral association (delta O-18(SMOW) = +7.2 to + 10.20, delta C-13(PDB) = 10.9 to - 1.80 for ankerite) and a later calcite (delta O-18(SMOW) = +5.9 to + 14.60, delta C-13(PDB) = -11.4 to + 4.40). In the southern Bowen Basin, the carbonate phase consists only of calcite (delta O-18(SMOW) = +12.5 to + 14.80, delta C-13(PDB) = -19.4 to + 0.80), where it occurs extensively throughout all stratigraphic levels. Siderite I occurs in mudrocks and sandstones and predates all other carbonate minerals. This carbonate phase is interpreted to have formed as an early diagenetic mineral from meteoric waters under cold climate and reducing conditions. Fe-Mg calcite-ankerite-siderite Il occur in sandstones as replacement of volcanic rock fragments. Clay minerals (illite-smectite, chlorite and kaolinite) postdate Ca-Fe-Mg carbonates, and precipitation of the later calcite is associated with clay mineral formation. The Ca-Fe-Mg carbonates and later calcite of the northern Bowen Basin are regarded as having formed as a result of hydrothermal activity during the latest Triassic extensional tectonic event which affected this part of the basin, rather than deep burial diagenesis during the Middle to Late Triassic as previously reported. This hypothesis is based on the timing relationships of the authigenic mineral phases and the low delta O-18 values of ankerite and calcite, together with radiometric dating of illitic clays and recently published regional geological evidence. Following the precipitation of the Ca-Fe-Mg carbonates from strongly O-18-depleted meteoric-hydrothermal fluids, continuing fluid circulation and water-rock interaction resulted in dissolution of these carbonate phases as well as labile fragments of volcaniclastic rocks. Subsequently, the later calcite and day minerals precipitated from relatively evolved (O-18-enriched) fluids. The nearly uniform delta O-18 values of the southern Bowen Basin calcite have been attributed to very low water/rock ratio in the system, where the fluid isotropic composition was buffered by the delta O-18 values of rocks. (C) 2000 Elsevier Science B.V. All rights reserved.

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.

Hybrid modelling of coupled pore fluid-solid deformation problems

A hybrid formulation for coupled pore fluid-solid deformation problems is proposed. The scheme is a hybrid in the sense that we use a vertex centered finite volume formulation for the analysis of the pore fluid and a particle method for the solid in our model. The pore fluid formally occupies the same space as the solid particles. The size of the particles is not necessarily equal to the physical size of materials. A finite volume mesh for the pore fluid flow is generated by Delaunay triangulation. Each triangle possesses an initial porosity. Changes of the porosity are specified by the translations of the mass centers of particles. Net pore pressure gradients are applied to the particle centers and are considered in the particle momentum balance. The potential of our model is illustrated by means of a simulation of coupled fracture and fluid flow developed in porous rock under biaxial compression condition.

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.