Transport of simple fluids in nanopores: Theory and simulation

A theory is discussed of single-component transport in nanopores, recently developed by Bhatia and coworkers. The theory considers the oscillatory motion of molecules between diffuse wall collisions, arising from the fluid-wall interaction, along with superimposed viscous flow due to fluid-fluid interaction. The theory is tested against molecular dynamics simulations for hydrogen, methane, and carbon tetrafluoride flow in cylindrical nanopores in silica. Although exact at low densities, the theory performs well even at high densities, with the density dependency of the transport coefficient arising from viscous effects. Such viscous effects are reduced at high densities because of the large increase in viscosity, which explains the maximum in the transport coefficient with increase in density. Further, it is seen that in narrow pore sizes of less than two molecular diameters, where a complete monolayer cannot form on the surface, the mutual interference of molecules on opposite sides of the cross section can reduce the transport coefficient, and lead to a maximum in the transport coefficient with increasing density. The theory is also tested for the case of partially diffuse reflection and shows the viscous contribution to be negligible when the reflection is nearly specular. (c) 2005 American Institute of Chemical Engineers AIChE J, 52: 29-38, 2006.

Effects of graphene layer size on the adsorption of fluids on graphitized thermal carbon black. A computer simulation study

The adsorption of Lennard-Jones fluids (argon and nitrogen) onto a graphitized thermal carbon black surface was studied with a Grand Canonical Monte Carlo Simulation (GCMC). The surface was assumed to be finite in length and composed of three graphene layers. When the GCMC simulation was used to describe adsorption on a graphite surface, an over-prediction of the isotherm was consistently observed in the pressure regions where the first and second layers are formed. To remove this over-prediction, surface mediation was accounted for to reduce the fluid-fluid interaction. Do and co-workers have introduced the so-called surface-mediation damping factor to correct the over-prediction for the case of a graphite surface of infinite extent, and this approach has yielded a good description of the adsorption isotherm. In this paper, the effects of the finite size of the graphene layer on the adsorption isotherm and how these would affect the extent of the surface mediation were studied. It was found that this finite-surface model provides a better description of the experimental data for graphitized thermal carbon black of high surface area (i.e. small crystallite size) while the infinite- surface model describes data for carbon black of very low surface area (i.e. large crystallite size).

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.

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.

Pore size characterization of mesoporous materials by a thermodynamic approach: A curvature-dependent solid-fluid potential

A thermodynamic analysis of nitrogen adsorption in cylindrical pores of MCM-41 and SBA-15 samples at 77 K is presented within the framework of the Broekhoff and de Boer (BdB) theory. We accounted for the effect of the solid surface curvature on the potential exerted by the pore walls. The developed model is in quantitative agreement with the non-local density functional theory (NLDFT) for pores larger than 2 tun. This modified BdB theory accounting for the Curvature Dependent Potential (CDP-BdB) was applied to determine the pore size distribution (PSD) of a number of MCM-41 and SBA-15 samples on the basis of matching the equilibrium theoretical isotherm against the adsorption branch of the experimental isotherm. In all cases investigated the PSDs determined with the new approach are very similar to those determined with the non-local density functional theory also using the same basis of matching of theoretical isotherm against the experimental adsorption branch. The developed continuum theory is very simple in its utilization, suggesting that CDP-BdB could be used as an alternative tool to obtain PSD for mesoporous solids from the analysis of adsorption branch of adsorption isotherms of any sub-critical fluids.

Influence of adsorbate interaction on transport in confined spaces

This article provides a review of the recent theory of transport in nanopores developed in the author's laboratory. In particular the influence of fluid-solid interactions on the transport coefficient is examined, showing that such interactions reduce the value of the coefficient by almost an order of magnitude in comparison to the Knudsen theory for non-interacting systems. The activation energy and potential energy barriers for diffusion in smooth pores with a one-dimensional potential energy profile are also discussed, indicating the inadequacy of the commonly used assumption of proportionality between the activation energy and heat of adsorption or the minimum pore potential energy. A further feature affected by fluid-solid interactions is the nature of the reflection of fluid molecules colliding with a pore wall surface, varying from being nearly specular - such as in carbon nanotubes - to nearly diffuse for amorphous solids. Diffuse reflection leads to momentum loss and reduced transport coefficients. However, fluid-solid interactions do not affect the transport coefficient in the single-file diffusion regime when the surface reflection is diffuse, and the transport coefficient in this case is largely independent of the adsorbed density.

Dissolution of a solid sphere in an unbounded, stagnant liquid

An exact analytical solution is obtained for the transient dissolution of solid spheres in a diffusion-controlled environment. This result provides a useful reference point for drug testing in humans. The dimensionless solution is expressed in terms of a single parameter, which accounts for solubility, bulk flow, and stagnant fluid composition. A simple, explicit and exact expression was found to predict time-to-complete dissolution (TCD). An approximate solution was also found which tracks the exact case for low solubility conditions.

Analysis of pore-fluid pressure gradient and effective vertical-stress gradient distribution in layered hydrodynamic systems

A theoretical analysis is carried out to investigate the pore-fluid pressure gradient and effective vertical-stress gradient distribution in fluid saturated porous rock masses in layered hydrodynamic systems. Three important concepts, namely the critical porosity of a porous medium, the intrinsic Fore-fluid pressure and the intrinsic effective vertical stress of the solid matrix, are presented and discussed. Using some basic scientific principles, we derive analytical solutions and explore the conditions under which either the intrinsic pore-fluid pressure gradient or the intrinsic effective vertical-stress gradient can be maintained at the value of the lithostatic pressure gradient. Even though the intrinsic pore-fluid pressure gradient can be maintained at the value of the lithostatic pressure gradient in a single layer, it is impossible to maintain it at this value in all layers in a layered hydrodynamic system, unless all layers have the same permeability and porosity simultaneously. However, the intrinsic effective vertical-stress gradient of the solid matrix can be maintained at a value close to the lithostatic pressure gradient in all layers in any layered hydrodynamic system within the scope of this study.

Modeling the fluid-flow-induced stress and collapse in a dendritic network

An analytical approach to the stress development in the coherent dendritic network during solidification is proposed. Under the assumption that stresses are developed in the network as a result of the friction resisting shrinkage-induced interdendritic fluid flow, the model predicts the stresses in the solid. The calculations reflect the expected effects of postponed dendrite coherency, slower solidification conditions, and variations of eutectic volume fraction and shrinkage. Comparing the calculated stresses to the measured shear strength of equiaxed mushy zones shows that it is possible for the stresses to exceed the strength, thereby resulting in reorientation or collapse of the dendritic network.

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.

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.

Functional cross-talk between cytokine receptors revealed by activating mutations in the extracellular domain of the beta-subunit of the GM-CSF receptor

Several reports have suggested an interaction between the erythropoietin receptor (EpoR) and the shared signaling subunit (hbeta(c)) of the human granulocyte macrophage-colony stimulating factor (GM-CSF), interleukin (IL)-3, and IL-5 receptors, although the functional consequences of this interaction are unclear. We previously showed that in vivo expression of constitutively active extracellular (EC) mutants of hbeta(c) induces erythrocytosis and Epo independence of erythroid colony-forming units (CFU-E). This occurs despite an apparent requirement of these mutants for the GM-CSF receptor alpha-subunit (GMRalpha), which is not expressed in CFU-E. Here, we show that coexpression of hbeta(c) EC mutants and EpoR in BaF-B03 cells, which lack GMRalpha, results in factor-independent proliferation and JAK2 activation. Mutant receptors that cannot activate JAK2 fail to produce a functional interaction. As there is no detectable phosphorylation of hbeta(c). on intracellular tyrosine residues, EpoR displays constitutive tyrosine phosphorylation. These observations suggest that JAK2 activation mediates cross-talk between EC mutants of hbeta(c) and EpoR. The implications of these data are discussed as are our findings that activated hbeta(c) mutants can functionally interact with certain other cytokine receptors.

Epstein-Barr virus encephalitis in a renal allograft recipient diagnosed by polymerase chain reaction on cerebrospinal fluid and successfully treated with ganciclovir

Epstein–Barr virus (EBV) encephalitis has been reported rarely in the context of solid-organ and bone-marrow transplantation [1]. We report a case of a renal transplant recipient who developed EBV encephalitis following OKT3 therapy for acute allograft rejection. The diagnosis was expedited by the detection of EBV DNA in the cerebrospinal fluid (CSF) by nested polymerase chain reaction (PCR). Moreover, clinical recovery and clearance of CSF EBV DNA appeared to follow the institution of parenteral ganciclovir treatment.

Diverse solid-state and solution structures within a series of hexaamine dicopper(II) complexes

Mono- and dicopper(II) complexes of a series of potentially bridging hexaamine ligands have been prepared and characterized in the solid state by X-ray crystallography. The crystal structures of the following Cu-II complexes are reported: [Cu(HL3)](ClO4)(3), C11H31Cl3CuN6O12, monoclinic, P2(1)/n, a = 8.294(2) Angstrom, b = 18.364(3) Angstrom, c = 15.674(3) Angstrom, beta = 94.73(2)degrees, Z = 4; {[Cu-2(L-4)(CO3)](2)}(ClO4)(4). 4H(2)O, C40H100Cl4Cu4N12O26, triclinic, P (1) over bar, a = 9.4888(8) Angstrom, b=13.353(1) Angstrom,. c = 15.329(1) Angstrom, alpha = 111.250(7)degrees, beta = 90.068(8)degrees, gamma = 105.081(8)degrees, Z=1; [Cu-2(L-5)(OH2)(2)](ClO4)(4), C(13)H(36)Cl(4)Cu(2)Z(6)O(18), monoclinic, P2(1)/c, a = 7.225(2) Angstrom. b = 8.5555(5) Angstrom, c = 23.134(8) Angstrom, beta = 92.37(1)degrees, Z = 2; [Cu-2(L-6)(OH2)(2)](ClO4)(4). 3H(2)O, C14H44Cl4Cu2N6O21, monoclinic, P2(1)/a, a = 15.204(5) Angstrom, b = 7.6810(7) Angstrom, c = 29.370(1) Angstrom, beta = 100.42(2)degrees, Z = 4. Solution spectroscopic properties of the bimetallic complexes indicate that significant conformational changes occur upon dissolution, and this has been probed with EPR spectroscopy and molecular mechanics calculations.

Fluid shifts resulting from exercise in rats as detected by bioelectrical impedance

Purpose: This study was designed to investigate the immediate effect of exercise intensity and duration on body fluid volumes in rats throughout a 3-wk exercise program. Methods: Changes in the extracellular water (ECW) and total body water (TBW) volumes of rats were measured preexercise and postexercise using multiple frequency bioelectrical impedance analysis. Groups of rats were exercised at two intensities (6 m.min(-1) and 12 m.min(-1)) for two exercise times (60 min and 90 min) 5 d.wk(-1) during a 3-wk period. Changes in plasma electrolytes, glucose, and lactate resulting from the exercise were also measured on 3 d of each week. Results: Each group of animals showed significant losses in ECW and TBW as a direct result of daily exercise. The magnitude of fluid loss was directly related to the intensity of the exercise, bur not to exercise duration; although the magnitude of daily fluid loss at the higher intensity exercise (12 m.min(-1)) decreased as the study progressed, possibly indicating a training effect. Conclusion: At low-intensity exercise, there is a small bur significant loss in both TBW and ECW fluids, and the magnitude of these losses does not change throughout a 3-wk exercise program. At moderate levels of exercise intensity, there is a greater loss of both TBW and ECW fluids. However, the magnitudes of these losses decrease significantly during the 3-wk exercise program, thus demonstrating a training effect.