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

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.

Characterization of pore wall heterogeneity in nanoporous carbons using adsorption: the slit pore model revisited

In this paper, we propose a new nonlocal density functional theory characterization procedure, the finite wall thickness model, for nanoporous carbons, whereby heterogeneity of pore size and pore walls in the carbon is probed simultaneously. We determine the pore size distributions and pore wall thickness distributions of several commercial activated carbons and coal chars, with good correspondence with X-ray diffraction. It is shown that the conventional infinite wall thickness approach overestimates the pore size slightly. Pore-pore correlation has been shown to have a negligible effect on prediction of pore size and pore wall thickness distributions for small molecules such as argon used in characterization. By utilizing the structural parameters (pore size and pore wall thickness distribution) in the generalized adsorption isotherm (GAI) we are able to predict adsorption uptake of supercritical gases in BPL and Norit RI Extra carbons, in excellent agreement with experimental adsorption uptake data up to 60 MPa. The method offers a useful technique for probing features of the solid skeleton, hitherto studied by crystallographic methods.

Improved comparison plot method for pore structure characterization of MCM-41

MCM-41 samples of various pore dimensions are synthesized. Plotting of nitrogen adsorption data at 77 K versus the statistical film thickness (comparison plot) reveals three distinct stages, with a characteristic of two points of inflection. The steep intermediate stage caused by capillary condensation occurred in the highly uniform mesopores. From the slopes of the sections before and after the condensation, the surface area of the mesopores is calculated. The linear portion of the last section is extrapolated to the adsorption axis of the comparison plot, and this intercept is used to obtain the volume of the mesopores. From the surface area and pore volume, average mesopore diameter is calculated, and the value thus obtained is in good agreement with the pore dimension obtained from powder X-ray diffraction measurements. The principle of the calculation as well as problems associated are discussed in detail.

Quantification of orientation of pore patterns in X-ray images of deformed clay

The use of X-ray imaging to assess variability in ceramic fabrication is common in archaeological studies aimed at examining ancient pottery technologies. In this paper, a method based on the measurement of individual pores orientation is presented. This method is successfully applied to ceramic specimens of known origin whose structure signified different deformation histories.

Pore structure characterisation of pillared clays using a modified MP method

The data of nitrogen adsorption on pillared clays (PILC) are converted to comparison plots (t-plots) to derive their pore size distribution (PSD). As in the MP method, the surface area of a group of pores having similar pore sizes is calculated from the slopes of tangent lines at two succeeding points on a comparison plot. By the modified MP method in this work, the tangent line is extrapolated to the adsorption axis on the t-plot, and the difference between intercepts is used to obtain the volume of the group of pores. From the information of surface area and pore volume, the average width of the pore group can be calculated and hence the PSDs of PILCs are obtained by carrying out such calculation procedures from high to low t. With this method, PSDs of several pillared clays are calculated over a wide pore size range, from micropores to mesopores. It is found that the modified MP method could result in the underestimation of the width of ultramicropores due to the enhancement in adsorption energy in these pores. Nevertheless, the method can be very useful in calculating the surface area and pore volume, as well as a mean width of these pores. For super-micropores and mesopores, pore size can also be underestimated, due to deviation of the pore shape from a slit. The principles of the improved MP method, as well as problems associated with it are thoroughly discussed in this paper. In general, this modified method provides practically meaningful results which are consistent with the pore dimension obtained from powder X-ray diffraction measurements, but involves no complicated theoretical treatment or assumptions.

A two-dimensional fuzzy random model of soil pore structure

A new conceptual model for soil pore-solid structure is formalized. Soil pore-solid structure is proposed to comprise spatially abutting elements each with a value which is its membership to the fuzzy set ''pore,'' termed porosity. These values have a range between zero (all solid) and unity (all pore). Images are used to represent structures in which the elements are pixels and the value of each is a porosity. Two-dimensional random fields are generated by allocating each pixel a porosity by independently sampling a statistical distribution. These random fields are reorganized into other pore-solid structural types by selecting parent points which have a specified local region of influence. Pixels of larger or smaller porosity are aggregated about the parent points and within the region of interest by controlled swapping of pixels in the image. This creates local regions of homogeneity within the random field. This is similar to the process known as simulated annealing. The resulting structures are characterized using one-and two-dimensional variograms and functions describing their connectivity. A variety of examples of structures created by the model is presented and compared. Extension to three dimensions presents no theoretical difficulties and is currently under development.

Estimation of the pore size of the large-conductance mechanosensitive ion channel of Escherichia coli

The open channel diameter of Escherichia coli recombinant large-conductance mechanosensitive ion channels (MscL) was estimated using the model of Hille (Hille, B. 1968. Pharmacological modifications of the sodium channels of frog nerve. J. Gen. Physiol. 51:199-219)that relates the pore size to conductance. Based on the MscL conductance of 3.8 nS, and assumed pore lengths, a channel diameter of 34 to 46 Angstrom was calculated. To estimate the pore size experimentally, the effect of large organic ions on the conductance of MscL was examined. Poly-L-lysines (PLLs) with a diameter of 37 Angstrom or larger significantly reduced channel conductance, whereas spermine (similar to 15 Angstrom), PLL19 (similar to 25 Angstrom) and 1,1'-bis-(3-(1'-methyl-(4,4'-bipyridinium)-1-yl)-propyl)-4,4'-bipyridinium (similar to 30 Angstrom) had no effect. The smaller organic ions putrescine, cadaverine, spermine, and succinate all permeated the channel. We conclude that the open pore diameter of the MscL is similar to 40 Angstrom, indicating that the MscL has one of the largest channel pores yet described. This channel diameter is consistent with the proposed homohexameric model of the MscL.

A modified pore-filling isotherm for liquid-phase adsorption in activated carbon

This article modifies the usual form of the Dubinin-Radushkevich pore-filling model for application to liquid-phase adsorption data, where large molecules are often involved. In such cases it is necessary to include the repulsive part of the energy in the micropores, which is accomplished here by relating the pore potential to the fluid-solid interaction potential. The model also considers the nonideality of the bulk liquid phase through the UNIFAC activity coefficient model, as well as structural heterogeneity of the carbon. For the latter the generalized adsorption integral is used while incorporating the pore-size distribution obtained by density functional theory analysis of argon adsorption data. The model is applied here to the interpretation of aqueous phase adsorption isotherms of three different esters on three commercial activated carbons. Excellent agreement between the model and experimental data is observed, and the fitted Lennard-Jones size parameter for the adsorbate-adsorbate interactions compares well with that estimated from known critical properties, supporting the modified approach. On the other hand, the model without consideration of bulk nonideality, or when using classical models of the characteristic energy, gives much poorer bts of the data and unrealistic parameter values.

Characterization of pore structure and coordination of titanium in TiO2 and SiO2-TiO2 sol-pillared clays

Titania sol-pillared clay (TiO2 PILC) and silica-titania sol-pillared clay (SiO2-TiO2 PILC) were synthesized by the sol-gel method. Supercritical drying (SCD) and treatment with quaternary ammonium surfactants were used to tailor the pore structure of the resulting clay. It was found that SCD approach increased the external surface area of the PILCs dramatically and that treatment with surfactants could be used to tailor pore size because the mesopore formation in the galleries between the clay layers follows the templating mechanism as observed in the synthesis of MCM-41 materials. Highly mesoporous solids were thus obtained. In calcined TiO2 PILC, ultrafine crystallites in anatase phase, which are active for photocatalytic oxidation of organics, were observed. In SiO2-TiO2 PILCs and their derivatives, titanium was highly dispersed in the matrix of silica and no crystal phase was observed. The highly dispersed titanium sites are good catalytic centers for selective oxidation of organic compounds. (C) 2001 Academic Press.