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

3D reconstruction and quantification of macropores using X-ray computed tomography and image analysis

Axial X-ray Computed tomography (CT) scanning provides a convenient means of recording the three-dimensional form of soil structure. The technique has been used for nearly two decades, but initial development has concentrated on qualitative description of images. More recently, increasing effort has been put into quantifying the geometry and topology of macropores likely to contribute to preferential now in soils. Here we describe a novel technique for tracing connected macropores in the CT scans. After object extraction, three-dimensional mathematical morphological filters are applied to quantify the reconstructed structure. These filters consist of sequences of so-called erosions and/or dilations of a 32-face structuring element to describe object distances and volumes of influence. The tracing and quantification methodologies were tested on a set of undisturbed soil cores collected in a Swiss pre-alpine meadow, where a new earthworm species (Aporrectodea nocturna) was accidentally introduced. Given the reduced number of samples analysed in this study, the results presented only illustrate the potential of the method to reconstruct and quantify macropores. Our results suggest that the introduction of the new species induced very limited chance to the soil structured for example, no difference in total macropore length or mean diameter was observed. However. in the zone colonised by, the new species. individual macropores tended to have a longer average length. be more vertical and be further apart at some depth. Overall, the approach proved well suited to the analysis of the three-dimensional architecture of macropores. It provides a framework for the analysis of complex structures, which are less satisfactorily observed and described using 2D imaging. (C) 2002 Elsevier Science B.V. All rights reserved.

Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

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.

Comparative surface accessibility of a pore-lining threonine residue (T6') in the glycine and GABA(A) receptors

The substituted cysteine accessibility method was used to probe the surface exposure of a pore-lining threonine residue (T6') common to both the glycine receptor (GlyR) and gamma-aminobutyric acid, type A receptor (GABAAR) chloride channels. This residue lies close to the channel activation gate, the ionic selectivity filter, and the main pore blocker binding site. Despite their high amino acid sequence homologies and common role in conducting chloride ions, recent studies have suggested that the GlyRs and GABA(A)Rs have divergent open state pore structures at the 6' position. When both the human alpha1(T6'C) homomeric GlyR and the rat alpha1(T6'C)beta1(T6'C) heteromeric GABA(A)R were expressed in human embryonic kidney 293 cells, their 6' residue surface accessibilities differed significantly in the closed state. However, when a soluble cysteine-modifying compound was applied in the presence of saturating agonist concentrations, both receptors were locked into the open state. This action was not induced by oxidizing agents in either receptor. These results provide evidence for a conserved pore opening mechanism in anion-selective members of the ligand-gated ion channel family. The results also indicate that the GABA(A)R pore structure at the 6' level may vary between different expression systems.

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.

Comparative surface accessibility of a pore-lining threonine residue (T6') in the glycine and GABA(A) receptors

The substituted cysteine accessibility method was used to probe the surface exposure of a pore-lining threonine residue (T6’) common to both the glycine receptor (GlyR) and GABAA receptor (GABAAR) chloride channels. This residue lies close to the channel activation gate, the ionic selectivity filter and the main pore blocker binding site. Recent studies have suggested that the GlyRs and GABAARs have divergent open state pore structures at the 6’ position. When both the human a1T6’C homomeric GlyR and the rat a1T6’Cb1T6’C heteromeric GABAAR were expressed in HEK293 cells, their 6’ residue surface accessibilities differed significantly in the closed state. However, when a soluble cysteine-modifying compound was applied in the presence of saturating agonist concentrations, both receptors were locked into the open state. This action was not induced by oxidising agents in either receptor. These results provide evidence for a conserved pore opening mechanism in anion-selective members of the ligand-gated ion channel family. The results also indicate that the GABAAR pore structure at the 6’ level may vary between different expression systems.

Site-directed mutagenesis of dimethylsulfoxide reductase from Rhodobacter capsulatus: The critical roles of tyrosine-114 and tryptophan-116

The crystal structure of six functionally-distinct enzymes of the DMSO reductase family of molybdenum enzymes has revealed that the tertiary structure of the polypeptide that binds the bis(MGD)Mo cofactor is highly conserved. Differences in the catalytic properties of enzymes of this family are almost certainly dependent upon differences in the structure ofthe MO active site. In DMSO reductase from Rhodobacter species tryptophan- 116 (W 116) hydrogen-bonds to an 0x0 group coordinated to the MO ion. In addition a second amino acid side chain from tyrosine-114 (Y 114) is in close proximity to the 0x0 group. We have investigated the role of Y 114 and W 116 in DMSO reductase using site-directed mutagenesis,

Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity

Exact analytical solutions of the critical Rayleigh numbers have been obtained for a hydrothermal system consisting of a horizontal porous layer with temperature-dependent viscosity. The boundary conditions considered are constant temperature and zero vertical Darcy velocity at both the top and bottom of the layer. Not only can the derived analytical solutions be readily used to examine the effect of the temperature-dependent viscosity on the temperature-gradient driven convective flow, but also they can be used to validate the numerical methods such as the finite-element method and finite-difference method for dealing with the same kind of problem. The related analytical and numerical results demonstrated that the temperature-dependent viscosity destabilizes the temperature-gradient driven convective flow and therefore, may affect the ore body formation and mineralization in the upper crust of the Earth. Copyright (C) 2003 John Wiley Sons, Ltd.

Innate immune surveillance of spontaneous B cell lymphomas by natural killer cells and gamma delta T cells

Few studies have demonstrated that innate lymphocytes play a major role in preventing spontaneous tumor formation. We evaluated the development of spontaneous tumors in mice lacking beta-2 microglobulin (beta2m; and thus MHC class I, CD1d, and CD16) and/or perform, since these tumor cells would be expected to activate innate effector cells. Approximately half the cohort of perform gene-targeted mice succumbed to spontaneous disseminated B cell lymphomas and in mice that also lacked beta2m, the lymphomas developed earlier (by more than 100 d) and with greater incidence (84%). B cell lymphomas from perforin/beta2m gene-targeted mice effectively primed cell-mediated cytotoxicity and perform, but not IFN-gamma, IL-12, or IL-18, was absolutely essential for tumor rejection. Activated NK1.1(+) and gammadeltaTCR(+) T cells were abundant at the tumor site, and transplanted tumors were strongly rejected by either, or both, of these cell types. Blockade of a number of different known costimulatory pathways failed to prevent tumor rejection. These results reflect a critical role for NK cells and gammadeltaTCP(+) T cells in innate immune surveillance of B cell lymphomas, mediated by as yet undetermined pathway(s) of tumor recognition.

Infrarenal aortic diameter predicts all-cause mortality

Objective - To assess the relationship between infrarenal aortic diameter and subsequent all-cause mortality in men aged 65 years or older. Methods and Results - Aortic diameter was measured using ultrasound in 12 203 men aged 65 to 83 years as part of a trial of screening for abdominal aortic aneurysms. A range of cardiovascular risk factors was also documented. Mortality over the next 3 to 7 years was assessed using record linkage. Initial aortic diameter was categorized into 10 intervals, and the relationship between increasing diameter and subsequent mortality was explored using Cox proportional hazard models. Median diameter increased from 21.4 mm in the youngest men to 22.1 mm in the oldest men. The cumulative all-cause mortality increased in a graded fashion with increasing aortic diameter. Using the diameter interval 19 to 22 mm as the reference, the adjusted hazard ratio for all-cause mortality increased from 1.26 (95% CI: 1.09, 1.44; P = 0.001) for aortic diameters of 23 to 26 mm to 2.38 (95% CI: 1.22, 4.61; P = 0.011) for aortic diameters of 47 to 50 mm. Analysis of causes of death indicated that cardiovascular disease was an important contributor to this increase. Conclusion - Infrarenal aortic diameter is an independent marker of subsequent all-cause mortality.