Effects of surface heterogeneity on the adsorption of nitrogen on graphitized thermal carbon black

In this paper, we study the surface heterogeneity and the surface mediation on the intermolecular potential energy for nitrogen adsorption on graphitized thermal carbon black (GTCB). The surface heterogeneity is modeled as the random distribution of effective carbonyl functional groups on the graphite surface. The molecular parameters and the discrete charges of this carbonyl group are taken from Jorgensen, et al. (J. Am. Chem. Soc., (1984) 106, 6638) while those for nitrogen (dispersive parameters and discrete charges) are taken from Murthy et al. (Mol. Phys., (1983) 50, 531) in our Grand Canonical Monte Carlo (GCMC) simulation. The solid surface mediation in the reduction of intermolecular potential energy between two fluid molecules was taken from a recent work by Do et al. (Langmuir, (2004) 20, 7623). Our simulation results accounting for the surface heterogeneity and surface mediation on intermolecular potential energy were compared with the experimental data of nitrogen at 77 and 90 K. The solid-fluid dispersive parameters are determined from the Lorentz-Berthelot (LB) rule. The fraction of the graphite surface covered with carbonyl functional groups was then derived from the consideration of the Henry constant, and for the data of Kruk et al. (Langmuir, (1999) 15, 1435) we have found that 1% of their GTCB surface is covered with effective carbonyl functional groups. The damping constant, due to surface mediation, was determined from the consideration of the portion of the adsorption isotherm where the first layer is being completed, and it was found to take a value of 0.0075. With these parameters, we have found that the GCMC simulation results describe the data over the complete range of pressure substantially better than any other MC models in the literature. The implication of this work is demonstrated with local adsorption isotherms of 10 and 20 A slit pores. One was obtained without allowance for surface mediation, while the other correctly accounts for these factors. The two local isotherms differ substantially, and the implication is that if we used incorrect local isotherms (i.e. without the surface mediation) the pore size distribution would be incorrectly derived.

Adsorption of argon on homogeneous graphitized thermal carbon black and heterogeneous carbon surface

In this paper we investigate the effects of surface mediation on the adsorption behavior of argon at different temperatures on homogeneous graphitized thermal carbon black and on heterogeneous nongraphitized carbon black surface. The grand canonical Monte Carlo (GCMC) simulation is used to study the adsorption, and its performance is tested against a number of experimental data on graphitized thermal carbon black (which is known to be highly homogeneous) that are available in the literature. The surface-mediation effect is shown to be essential in the correct description of the adsorption isotherm because without accounting for that effect the GCMC simulation results are always greater than the experimental data in the region where the monolayer is being completed. This is due to the overestimation of the fluid–fluid interaction between particles in the first layer close to the solid surface. It is the surface mediation that reduces this fluid–fluid interaction in the adsorbed layers, and therefore the GCMC simulation results accounting for this surface mediation that are presented in this paper result in a better description of the data. This surface mediation having been determined, the surface excess of argon on heterogeneous carbon surfaces having solid–fluid interaction energies different from the graphite can be readily obtained. Since the real heterogeneous carbon surface is not the same as the homogeneous graphite surface, it can be described by an area distribution in terms of the well depth of the solid–fluid energy. Assuming a patchwise topology of the surface with patches of uniform well depth of solid–fluid interaction, the adsorption on a real carbon surface can be determined as an integral of the local surface excess of each patch with respect to the differential area. When this is matched against the experimental data of a carbon surface, we can derive the area distribution versus energy and hence the geometrical surface area. This new approach will be illustrated with the adsorption of argon on a nongraphitized carbon at 87.3 and 77 K, and it is found that the GCMC surface area is different from the BET surface area by about 7%. Furthermore, the description of the isotherm in the region of BET validity of 0.06 to 0.2 is much better with our method than with the BET equation.

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.

Kinetic study of the thermal degradation of high density polyethylene

The thermal degradation of high density polyethylene has been modelled by the random breakage of polymer bonds, using a set of population balance equations. A model was proposed in which the population balances were lumped into representative sizes so that the experimentally determined molecular weight distribution of the original polymer could be used as the initial condition. This model was then compared to two different cases of the unlumped population balance which assumed unimolecular initial distributions of 100 and 500 monomer units, respectively. The model that utilised the experimentally determined molecular weight distribution was found to best describe the experimental data. The model fits suggested a second mechanism in addition to random breakage at slow reaction rates. (c) 2005 Elsevier Ltd. All rights reserved.

Thermal Cycling Stability of Silica Membranes for Gas Separation

Hydrogen is being seen as an alternative energy carrier to conventional hydrocarbons to reduce greenhouse gas emissions. High efficiency separation technologies to remove hydrogen from the greenhouse gas, carbon dioxide, are therefore in growing demand. Traditional thermodynamic separation systems utilise distillation, absorption and adsorption, but are limited in efficiency at compact scales. Molecular sieve silica (MSS) membranes can perform this separation as they have high permselectivity of hydrogen to carbon dioxide, but their stability under thermal cycling is not well reported. In this work we exposed a standard MSS membrane and a carbonised template MSS (CTMSS) membrane to thermal cycling from 100 to 450°C. The standard MSS and carbonised template CTMSS membranes both showed permselectivity of helium to nitrogen dropping from around 10 to 6 in the first set of cycles, remaining stable until the last test. The permselectivity drop was due to small micropore collapse, which occurred via structure movement during cycling. Simulating single stage membrane separation with a 50:50 molar feed of H2:CO2, H2 exiting the permeate stream would start at 79% and stabilise at 67%. Higher selectivity membranes showed less of a purity drop, indicating the margin at which to design a stable membrane separation unit for CO2 capture.

Thermal effects in the evolution of initially layered mantle material

A simplified model for anisotropic mantle convection based on a novel class of rheologies, originally developed for folding instabilities in multilayered rock (MUHLHAUS et al., 2002), is extended ¨ through the introduction of a thermal anisotropy dependent on the local layering. To examine the effect of the thermal anisotropy on the evolution of mantle material, a parallel implementation of this model was undertaken using the Escript modelling toolkit and the Finley finite-element computational kernel (DAVIES et al., 2004). For the cases studied, there appears too little if any effect. For comparative purposes, the effects of anisotropic shear viscosity and the introduced thermal anisotropy are also presented. These results contribute to the characterization of viscous anisotropic mantle convection subject to variation in thermal conductivities and shear viscosities.

Growth of Bose-Einstein condensates from thermal vapor

We report on a quantitative study of the growth process of 87Rb Bose-Einstein condensates. By continuous evaporative cooling we directly control the thermal cloud from which the condensate grows. We compare the experimental data with the results of a theoretical model based on quantum kinetic theory. We find quantitative agreement with theory for the situation of strong cooling, whereas in the weak cooling regime a distinctly different behavior is found in the experiment.

XSophe-Sophe-XeprView (R). A computer simulation software suite (v. 1.1.3) for the analysis of continuous wave EPR spectra

The XSophe-Sophe-XeprView((R)) computer simulation software suite enables scientists to easily determine spin Hamiltonian parameters from isotropic, randomly oriented and single crystal continuous wave electron paramagnetic resonance (CW EPR) spectra from radicals and isolated paramagnetic metal ion centers or clusters found in metalloproteins, chemical systems and materials science. XSophe provides an X-windows graphical user interface to the Sophe programme and allows: creation of multiple input files, local and remote execution of Sophe, the display of sophelog (output from Sophe) and input parameters/files. Sophe is a sophisticated computer simulation software programme employing a number of innovative technologies including; the Sydney OPera HousE (SOPHE) partition and interpolation schemes, a field segmentation algorithm, the mosaic misorientation linewidth model, parallelization and spectral optimisation. In conjunction with the SOPHE partition scheme and the field segmentation algorithm, the SOPHE interpolation scheme and the mosaic misorientation linewidth model greatly increase the speed of simulations for most spin systems. Employing brute force matrix diagonalization in the simulation of an EPR spectrum from a high spin Cr(III) complex with the spin Hamiltonian parameters g(e) = 2.00, D = 0.10 cm(-1), E/D = 0.25, A(x) = 120.0, A(y) = 120.0, A(z) = 240.0 x 10(-4) cm(-1) requires a SOPHE grid size of N = 400 (to produce a good signal to noise ratio) and takes 229.47 s. In contrast the use of either the SOPHE interpolation scheme or the mosaic misorientation linewidth model requires a SOPHE grid size of only N = 18 and takes 44.08 and 0.79 s, respectively. Results from Sophe are transferred via the Common Object Request Broker Architecture (CORBA) to XSophe and subsequently to XeprView((R)) where the simulated CW EPR spectra (1D and 2D) can be compared to the experimental spectra. Energy level diagrams, transition roadmaps and transition surfaces aid the interpretation of complicated randomly oriented CW EPR spectra and can be viewed with a web browser and an OpenInventor scene graph viewer.

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

User guide for shock and blast simulation with the OctVCE code (version 3.5+), Department of Mechanical Engineering Report 2007/13

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel as it is simple to code and sufficient for practical engineering design problems. This also makes the code much more ‘user-friendly’ than structured grid approaches as the gridding process is done automatically. The CFD methodology relies on a finite-volume formulation of the unsteady Euler equations and is solved using a standard explicit Godonov (MUSCL) scheme. Both octree-based adaptive mesh refinement and shared-memory parallel processing capability have also been incorporated. For further details on the theory behind the code, see the companion report 2007/12.

Flow and axial dispersion simulation for traveling axisymmetric Taylor vortices

Numerical experiments using a finite difference method were carried out to determine the motion of axisymmetric Taylor vortices for narrow-gap Taylor vortex flow. When a pressure gradient is imposed on the flow the vortices are observed to move with an axial speed of 1.16 +/- 0.005 times the mean axial flow velocity. The method of Brenner was used to calculate the long-time axial spread of material in the flow. For flows where there is no pressure gradient, the axial dispersion scales with the square root of the molecular diffusion, in agreement with the results of Rosen-bluth et al. for high Peclet number dispersion in spatially periodic flows with a roll structure. When a pressure gradient is imposed the dispersion increases by an amount approximately equal to 6.5 x 10(-4) (W) over bar(2)d(2)/D-m, where (W) over bar is the average axial velocity in the annulus, analogous to Taylor dispersion for laminar flow in an empty tube.

Improved methods for predicting individual leaf area and leaf senescence in maize (Zea mays)

The ability to predict leaf area and leaf area index is crucial in crop simulation models that predict crop growth and yield. Previous studies have shown existing methods of predicting leaf area to be inadequate when applied to a broad range of cultivars with different numbers of leaves. The objectives of the study were to (i) develop generalised methods of modelling individual and total plant leaf area, and leaf senescence, that do not require constants that are specific to environments and/or genotypes, (ii) re-examine the base, optimum, and maximum temperatures for calculation of thermal time for leaf senescence, and (iii) assess the method of calculation of individual leaf area from leaf length and leaf width in experimental work. Five cultivars of maize differing widely in maturity and adaptation were planted in October 1994 in south-eastern Queensland, and grown under non-limiting conditions of water and plant nutrient supplies. Additional data for maize plants with low total leaf number (12-17) grown at Katumani Research Centre, Kenya, were included to extend the range in the total leaf number per plant. The equation for the modified (slightly skewed) bell curve could be generalised for modelling individual leaf area, as all coefficients in it were related to total leaf number. Use of coefficients for individual genotypes can be avoided, and individual and total plant leaf area can be calculated from total leaf number. A single, logistic equation, relying on maximum plant leaf area and thermal time from emergence, was developed to predict leaf senescence. The base, optimum, and maximum temperatures for calculation of thermal time for leaf senescence were 8, 34, and 40 degrees C, and apply for the whole crop-cycle when used in modelling of leaf senescence. Thus, the modelling of leaf production and senescence is simplified, improved, and generalised. Consequently, the modelling of leaf area index (LAI) and variables that rely on LAI will be improved. For experimental purposes, we found that the calculation of leaf area from leaf length and leaf width remains appropriate, though the relationship differed slightly from previously published equations.