Effect of surface-perturbed intermolecular interaction on adsorption of simple gases on a graphitized carbon surface

In this paper, we investigate the effect of the solid surface on the fluid-fluid intermolecular potential energy. This modified fluid-fluid interaction energy due to the inducement of a solid surface is used in the grand canonical Monte Carlo (GCMC) simulation of various noble gases, nitrogen, and methane on graphitized thermal carbon black. This effect is such that the effective interaction potential energy between two particles close to surface is less than the potential energy if the solid substrate is not present. With this modification the GCMC simulation results agree extremely well with the experimental data over a wide range of pressures while the simulation results with the unmodified potential energy give rise to a shoulder near the neighborhood of monolayer coverage and the significant overprediction of the second and higher layer coverages. In particular the unmodified GCMC results exhibit very sharp change in those higher layers while the experimental data have a much gradual change in the uptake. We will illustrate this theory with adsorption data of argon, xenon, neon, nitrogen, and methane on graphitized thermal carbon black.

Capillary coexistence and criticality in mesopores: Modification of the Kelvin theory

The classical model of capillary equilibrium in cylindrical pores is modified here by the introduction of molecular concepts and the solid fluid interaction potential. The new approach accurately predicts capillary coexistence and criticality, with results quantitatively matching those from density functional theory for nitrogen adsorption, while also predicting condensation pressures in agreement with reported experimental findings for MCM-41. The larger critical pore size for nitrogen adsorption in these materials, however, suggests a modification of the potential function parameters, evaluated here from data for hydroxylated silica.

Non-additivity of attractive potentials in modeling of N-2 and Ar adsorption isotherms on graphitized carbon black and porous carbon by means of density functional theory

We present a new approach accounting for the nonadditivity of attractive parts of solid-fluid and fluidfluid potentials to improve the quality of the description of nitrogen and argon adsorption isotherms on graphitized carbon black in the framework of non-local density functional theory. We show that the strong solid-fluid interaction in the first monolayer decreases the fluid-fluid interaction, which prevents the twodimensional phase transition to occur. This results in smoother isotherm, which agrees much better with experimental data. In the region of multi-layer coverage the conventional non-local density functional theory and grand canonical Monte Carlo simulations are known to over-predict the amount adsorbed against experimental isotherms. Accounting for the non-additivity factor decreases the solid-fluid interaction with the increase of intermolecular interactions in the dense adsorbed fluid, preventing the over-prediction of loading in the region of multi-layer adsorption. Such an improvement of the non-local density functional theory allows us to describe experimental nitrogen and argon isotherms on carbon black quite accurately with mean error of 2.5 to 5.8% instead of 17 to 26% in the conventional technique. With this approach, the local isotherms of model pores can be derived, and consequently a more reliab * le pore size distribution can be obtained. We illustrate this by applying our theory against nitrogen and argon isotherms on a number of activated carbons. The fitting between our model and the data is much better than the conventional NLDFT, suggesting the more reliable PSD obtained with our approach.

Effects of the juxtaposition of carbonaceous slit pores on the overall transport behavior of adsorbed fluids

It is a common approximation in the modeling of adsorption in microporous carbons to treat the pores as slit pores, whose walls are considered to consist of an infinite number of graphitic layers. In practice, such an approximation is appropriate as long as the number of graphitic layers in the wall is greater than three. However, it is understood that pore walls in microporous carbons commonly consist of three or fewer layers. As well as affecting the solid-fluid interaction within a pore, such narrow walls permit the interaction of fluid molecules through the wall, with consequences for the adsorption characteristics. We consider the effect that a distributed pore-wall thickness model can have on transport properties. At low density we find that the only significant deviation in the transport properties from the infinite pore-wall thickness model occurs in pores with single-layer walls. For a model of activated carbons with a distribution of pore widths and pore-wall thicknesses, the transport properties are generally insensitive to the effects of finite walls, in terms of both the solid-fluid interaction within a pore and fluid-fluid interaction through the pore walls.

Impact of potential models on adsorption of linear molecules on carbon black

In this paper, we investigate the effects of potential models on the description of equilibria of linear molecules (ethylene and ethane) adsorption on graphitized thermal carbon black. GCMC simulation is used as a tool to give adsorption isotherms, isosteric heat of adsorption and the microscopic configurations of these molecules. At the heart of the GCMC are the potential models, describing fluid-fluid interaction and solid-fluid interaction. Here we studied the two potential models recently proposed in the literature, the UA-TraPPE and AUA4. Their impact in the description of adsorption behavior of pure components will be discussed. Mixtures of these components with nitrogen and argon are also studied. Nitrogen is modeled a two-site plus discrete charges while argon as a spherical particle. GCMC simulation is also used for generating simulation mixture isotherms. It is found that co-operation between species occurs when the surface is fractionally covered while competition is important when surface is fully loaded.

Modeling of adsorption on nongraphitized carbon surface: GCMC simulation studies and comparison with experimental data

We model nongraphitized carbon black surfaces and investigate adsorption of argon on these surfaces by using the grand canonical Monte Carlo simulation. In this model, the nongraphitized surface is modeled as a stack of graphene layers with some carbon atoms of the top graphene layer being randomly removed. The percentage of the surface carbon atoms being removed and the effective size of the defect ( created by the removal) are the key parameters to characterize the nongraphitized surface. The patterns of adsorption isotherm and isosteric heat are particularly studied, as a function of these surface parameters as well as pressure and temperature. It is shown that the adsorption isotherm shows a steplike behavior on a perfect graphite surface and becomes smoother on nongraphitized surfaces. Regarding the isosteric heat versus loading, we observe for the case of graphitized thermal carbon black the increase of heat in the submonolayer coverage and then a sharp decline in the heat when the second layer is starting to form, beyond which it increases slightly. On the other hand, the isosteric heat versus loading for a highly nongraphitized surface shows a general decline with respect to loading, which is due to the energetic heterogeneity of the surface. It is only when the fluid-fluid interaction is greater than the surface energetic factor that we see a minimum-maximum in the isosteric heat versus loading. These simulation results of isosteric heat agree well with the experimental results of graphitization of Spheron 6 (Polley, M. H.; Schaeffer, W. D.; Smith, W. R. J. Phys. Chem. 1953, 57, 469; Beebe, R. A.; Young, D. M. J. Phys. Chem. 1954, 58, 93). Adsorption isotherms and isosteric heat in pores whose walls have defects are also studied from the simulation, and the pattern of isotherm and isosteric heat could be used to identify the fingerprint of the surface.

Adsorption of carbon tetrachloride on graphitized thermal carbon black and in slit graphitic pores: Five-site versus one-site potential models

The performance of intermolecular potential models on the adsorption of carbon tetrachloride on graphitized thermal carbon black at various temperatures is investigated. This is made possible with the extensive experimental data of Machin and Ross(1), Avgul et al.,(2) and Pierce(3) that cover a wide range of temperatures. The description of all experimental data is only possible with the allowance for the surface mediation. If this were ignored, the grand canonical Monte Carlo (GCMC) simulation results would predict a two-dimensional (2D) transition even at high temperatures, while experimental data shows gradual change in adsorption density with pressure. In general, we find that the intermolecular interaction has to be reduced by 4% whenever particles are within the first layer close to the surface. We also find that this degree of surface mediation is independent of temperature. To understand the packing of carbon tetrachloride in slit pores, we compared the performance of the potential models that model carbon tetrachloride as either five interaction sites or one site. It was found that the five-site model performs better and describes the imperfect packing in small pores better. This is so because most of the strength of fluid-fluid interaction between two carbon tetrachloride molecules comes from the interactions among chlorine atoms. Methane, although having tetrahedral shape as carbon tetrachloride, can be effectively modeled as a pseudospherical particle because most of the interactions come from carbon-carbon interaction and hydrogen negligibly contributes to this.

Adsorption of benzene on graphitized thermal carbon black: Reduction of the quadrupole moment in the adsorbed phase

The performance of intermolecular potential models on the adsorption of benzene on graphitized thermal carbon black at various temperatures is investigated. Two models contain only dispersive sites, whereas the other two models account explicitly for the dispersive and electrostatic sites. Using numerous data in the literature on benzene adsorption on graphitized thermal carbon black at various temperatures, we have found that the effect of surface mediation on interaction between adsorbed benzene molecules must be accounted for to describe correctly the adsorption isotherm as well as the isosteric heat. Among the two models with partial charges tested, the WSKS model of Wick et at. I that has only six dispersive sites and three discrete partial charges is better than the very expensive all-atom model of Jorgensen and Severance.(2) Adsorbed benzene molecules on graphitized thermal carbon black have a complex orientation with respect to distance from the surface and also with respect to loading. At low loadings, they adopt the parallel configuration relative to the graphene surface, whereas at higher loadings (still less than monolayer coverage) some molecules adopt a slant orientation to maximize the fluid-fluid interaction. For loadings in the multilayer region, the orientation of molecules in the first layer is influenced by the presence of molecules in the second layer. The data that are used in this article come from the work of Isirikyan and Kiselev,(3) Pierotti and Smallwood,(4) Pierce and Ewing,(5) Belyakova, Kiselev, and Kovaleva,(6) and Carrott et al.(7)

Mechanisms influencing levitation and the scaling laws in nanopores: Oscillator model theory

We provide here a detailed theoretical explanation of the floating molecule or levitation effect, for molecules diffusing through nanopores, using the oscillator model theory (Phys. Rev. Lett. 2003, 91, 126102) recently developed in this laboratory. It is shown that on reduction of pore size the effect occurs due to decrease in frequency of wall collision of diffusing particles at a critical pore size. This effect is, however, absent at high temperatures where the ratio of kinetic energy to the solid-fluid interaction strength is sufficiently large. It is shown that the transport diffusivities scale with this ratio. Scaling of transport diffusivities with respect to mass is also observed, even in the presence of interactions.

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.

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.

Study of flow behaviour in a three-product cyclone using computational fluid dynamics

Simplicity in design and minimal floor space requirements render the hydrocyclone the preferred classifier in mineral processing plants. Empirical models have been developed for design and process optimisation but due to the complexity of the flow behaviour in the hydrocyclone these do not provide information on the internal separation mechanisms. To study the interaction of design variables, the flow behaviour needs to be considered, especially when modelling the new three-product cyclone. Computational fluid dynamics (CFD) was used to model the three-product cyclone, in particular the influence of the dual vortex finder arrangement on flow behaviour. From experimental work performed on the UG2 platinum ore, significant differences in the classification performance of the three-product cyclone were noticed with variations in the inner vortex finder length. Because of this simulations were performed for a range of inner vortex finder lengths. Simulations were also conducted on a conventional hydrocyclone of the same size to enable a direct comparison of the flow behaviour between the two cyclone designs. Significantly, high velocities were observed for the three-product cyclone with an inner vortex finder extended deep into the conical section of the cyclone. CFD studies revealed that in the three-product cyclone, a cylindrical shaped air-core is observed similar to conventional hydrocyclones. A constant diameter air-core was observed throughout the inner vortex finder length, while no air-core was present in the annulus. (c) 2006 Elsevier Ltd. All rights reserved.

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.