PVT behaviour of fluids with potential models optimised for phase equilibria

NPT and NVT Monte Carlo simulations are applied to models for methane and water to predict the PVT behaviour of these fluids over a wide range of temperatures and pressures. The potential models examined in this paper have previously been presented in the literature with their specific parameters optimised to fit phase coexistence data. The exponential-6 potential for methane gives generally good prediction of PVT behaviour over the full range of temperature and pressures studied with the only significant deviation from experimental data seen at high temperatures and pressures. The NSPCE water model shows very poor prediction of PVT behaviour, particularly at dense conditions. To improve this. the charge separation in the NSPCE model is varied with density. Improvements for vapour and liquid phase PVT predictions are achieved with this variation. No improvement was found in the prediction of the oxygen-oxygen radial distribution by varying charge separation under dense phase conditions. (C) 2004 Elsevier B.V. All rights reserved.

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.

Adsorption of ethylene on graphitized thermal carbon black and in slit pores: A computer simulation study

In this paper, we studied vapor-liquid equilibria (VLE) and adsorption of ethylene on graphitized thermal carbon black and in slit pores whose walls are composed of graphene layers. Simple models of a one-center Lennard-Jones (LJ) potential and a two-center united atom (UA)-LJ potential are investigated to study the impact of the choice of potential models in the description of VLE and adsorption behavior. Here, we used a Monte Carlo simulation method with grand canonical Monte Carlo (GCMC) and Gibbs ensemble Monte Carlo ensembles. The one-center potential model cannot describe adequately the VLE over the practical range of temperature from the triple point to the critical point. On the other hand, the two-center potential model (Wick et al. J. Phys. Chem. B 2000, 104, 8008-8016) performs well in the description of VLE (saturated vapor and liquid densities and vapor pressure) over the wide range of temperature. This UA-LJ model is then used in the study of adsorption of ethylene on graphitized thermal carbon black and in slit pores. Agreement between the GCMC simulation results and the experimental data on graphitized thermal carbon black for moderate temperatures is excellent, demonstrating that the potential of the GCMC method and the proper choice of potential model are essential to investigate adsorption. For slit pores of various sizes, we have found that the behavior of ethylene exhibits a number of features that are not manifested in the study of spherical LJ particles. In particular, the singlet density distribution versus distance across the pore and the angle between the molecular axis and the z direction provide rich information about the way molecules arrange themselves when the pore width is varied. Such an arrangement has been found to be very sensitive to the pore width.

Application of density functional theory to capillary phenomena in cylindrical mesopores with radial and longitudinal density distributions

In this paper, we applied a version of the nonlocal density functional theory (NLDFT) accounting radial and longitudinal density distributions to study the adsorption and desorption of argon in finite as well as infinite cylindrical nanopores at 87.3 K. Features that have not been observed before with one-dimensional NLDFT are observed in the analysis of an inhomogeneous fluid along the axis of a finite cylindrical pore using the two-dimensional version of the NLDFT. The phase transition in pore is not strictly vapor-liquid transition as assumed and observed in the conventional version, but rather it exhibits a much elaborated feature with phase transition being complicated by the formation of solid phase. Depending on the pore size, there are more than one phase transition in the adsorption-desorption isotherm. The solid formation in finite pore has been found to be initiated by the presence of the meniscus. Details of the analysis of the extended version of NLDFT will be discussed in the paper. (C) 2004 American Institute of Physics.

Application of density functional theory to analysis of energetic heterogeneity and pore size distribution of activated carbons

A new approach based on the nonlocal density functional theory to determine pore size distribution (PSD) of activated carbons and energetic heterogeneity of the pore wall is proposed. The energetic heterogeneity is modeled with an energy distribution function (EDF), describing the distribution of solid-fluid potential well depth (this distribution is a Dirac delta function for an energetic homogeneous surface). The approach allows simultaneous determining of the PSD (assuming slit shape) and EDF from nitrogen or argon isotherms at their respective boiling points by using a set of local isotherms calculated for a range of pore widths and solid-fluid potential well depths. It is found that the structure of the pore wall surface significantly differs from that of graphitized carbon black. This could be attributed to defects in the crystalline structure of the surface, active oxide centers, finite size of the pore walls (in either wall thickness or pore length), and so forth. Those factors depend on the precursor and the process of carbonization and activation and hence provide a fingerprint for each adsorbent. The approach allows very accurate correlation of the experimental adsorption isotherm and leads to PSDs that are simpler and more realistic than those obtained with the original nonlocal density functional theory.

Research on the BET surface area and packing of molecules on the activated carbon

Adsorption of different aromatic compounds (two of them are electrolytes) onto an untreated activated carbon (F100) is investigated. The experimental isotherms are fitted into Langmuir homogenous and heterogeneous Model. Theoretical maximum adsorption capacities that are based on the BET surface area of the adsorbent cannot be close to the real value. The affinity and the heterogeneity of the adsorption system observed to be related to the pK(a) of the solutes. The maximum adsorption capacity (Q(max)) of activated carbon for each solute dependent on the molecular area as well as the type of functional group attached on the aromatic compound and also pH of solution. The arrangement of the molecules on the carbon surface is not face down. Furthermore, it is illustrated that the packing arrangement is most likely edge to face (sorbate-sorbent) with various tilt angles. For characterization of the carbon, the N-2 and CO2 adsorption were used. X-ray Photoelectron Spectroscopy (XPS) measurement was used to surface elemental analysis of activated carbon.

Adsorption of p-nitrophenol in untreated and treated activated carbon

The adsorption of p-nitrophenol in one untreated activated carbon (F100) and three treated activated carbons (H-2, H2SO4 and Urea treated F100) was carried out at undissociated and dissociated conditions. To characterize the carbon, N-2 and CO2 adsorption were used. X-ray Photoelectron Spectroscopy (XPS) was used to analyze the surface of the activated carbon. The experimental isotherms are fitted via the Langmuir homogenous model and Langmuir binary model. Variation of the model parameters with the solution pH is studied. Both Q(max) and the adsorption affinity coefficient (K-1) were dependent on the PZC of the carbons and solution pH. The Effect of pH must be considered due to its combined effects on the carbon surface and on the solute molecules. Adsorption of p-nitrophenol at higher pH was found to be dependent on the concentration of the anionic form of the solute.

Diffusion of linear paraffins in nanoporous silica

The diffusion of hexane, heptane, octane, and decane in nanoporous MCM-41 silica at various temperatures is investigated by the zero-length-column method. The diffusion coefficients are derived by a complete-time-range analysis of desorption curves at different purge flow rates and temperatures. The results show that the calculated low-coverage diffusivity values decrease monotonically, and the derived Henry's law constants increase, as the carbon number of paraffins increases. The study reveals that transport is strongly influenced by intracrystalline diffusion and dominated by the sorbate-sorbent interaction. The diffusion activation energy and adsorption isosteric heat at zero loading increase monotonically with the carbon number of linear paraffins, but their ratio is essentially constant for each adsorbate compound.

Three-dimensional simulation of the rolling of a saturated fibro-porous media

This paper describes recent advances made in computational modelling of the sugar cane liquid extraction process. The saturated fibro-porous material is rolled between circumferentially grooved rolls, which enhance frictional grip and provide a low-resistance path for liquid flow during the extraction process. Previously reported two-dimensional (2D) computational models, account for the large deformation of the porous material by solving the fully coupled governing fibre stress and fluid-flow equations using finite element techniques. While the 2D simulations provide much insight into the overarching cause-effect relationships, predictions of mechanical quantities such as roll separating force and particularly torque as a function of roll speed and degree of compression are not satisfactory for industrial use. It is considered that the unsatisfactory response in roll torque prediction may be due to the stress levels that exist between the groove tips and roots which have been largely neglected in the geometrically simplified 2D model. This paper gives results for both two- and three-dimensional finite element models and highlights their strengths and weaknesses in predicting key milling parameters. (c) 2005 Elsevier B.V. All rights reserved.

Effect of metal chlorides on the sintering and densification of hydroxyapatite adsorbent

This work is part of a series of studies dealing with the evaluation of the effects of major elements of solid waste, especially metallic oxides, nitrates, sulfates, and chlorides, on the sintering and the densification of calcium hydroxyapatite (Ca-HAP) adsorbent. The effects of chloride salts of potassium (KCl) and zinc (ZnCl2) on sintering and densification of Ca-HAP were studied using surface area reduction and shrinkage measurements. The addition of KCl (2% w/w) activated the sintering process by bringing a swift reduction in surface area and lowering the densification temperature. However, a low final densification was achieved. Increasing the amount of this additive to 10% w/w further lowered the final densification and lowered the densification temperature of hydroxyapatite by 150 degrees C. On the other hand, the addition of 2 wt % of ZnCl2 deactivated the sintering process by slowing down the densification process and raising the densification temperature. However, the reduction of surface area was comparable to that of Ca-HAP. The densification rate contained two or more rate maxima indicating the additives (salts) bring multiple speeds in the densification process.

The physical and surface chemical characteristics of activated carbons and the adsorption of methylene blue from wastewater

Adsorption of a basic dye, methylene blue, from aqueous solutions onto as-received activated carbons and acid-treated carbons was investigated. The physical and surface chemical properties of the activated carbons were characterized using BET-N-2 adsorption, X-ray photoelectron spectroscopy (XPS), and mass titration. It was found that acid treatment had little effect on carbon textural characteristics but significantly changed the surface chemical properties, resulting in an adverse effect on dye adsorption. The physical properties of activated carbon, such as surface area and pore volume, have little effect on dye adsorption, while the pore size distribution and the surface chemical characteristics play important roles in dye adsorption. The pH value of the solution also influences the adsorption capacity significantly. For methylene blue, a higher pH of solution favors the adsorption capacity. The kinetic adsorption of methylene blue on all carbons follows a pseudo-second-order equation. (c) 2004 Elsevier Inc. All rights reserved.

Local liquid holdups and hysteresis in a 2-D packed bed using X-ray radiography

An X-ray visualization technique has been used for the quantitative determination of local liquid holdups distribution and liquid holdup hysteresis in a nonwetting two-dimensional (2-D) packed bed. A medical diagnostic X-ray unit has been used to image the local holdups in a 2-D cold model having a random packing of expanded polystyrene beads. An aqueous barium chloride solution was used as a fluid to achieve good contrast on X-ray images. To quantify the local liquid holdup, a simple calibration technique has been developed that can be used for most of the radiological methods such as gamma ray and neutron radiography. The global value of total liquid holdup, obtained by X-ray method, has been compared with two conventional methods: drainage and tracer response. The X-ray technique, after validation, has been used to visualize and quantify, the liquid hysteresis phenomena in a packed bed. The liquid flows in preferred paths or channels that carry droplets/rivulets of increasing size and number as the liquid flow rate is increased. When the flow is reduced, these paths are retained and the higher liquid holdup that persists in these regions leads to the holdup hysteresis effect. Holdup in some regions of the packed bed may be an order of magnitude higher than average at a particular flow rate. (c) 2005 American Institute of Chemical Engineers

The drying of sewage sludge by immersion frying

The objective of this work was to dry sewage sludge using a fry-drying process. The frying experiments were carried out in commercial fryers modified by adding thermocouples to the setup. During frying, typical drying curves were obtained and it was verified that, in relation to the parameters: oil temperature, oil type and shape of the sample, the shape factor the most effect on the drying rate, at least within the range chosen for the variables studied. Oil uptake and calorific value were also analyzed. The calorific value of the samples increased with frying time, reaching values around 24MJ/kg after 600s of frying (comparable to biocombustibles such as wood and sugarcane bagasse). The process of immersion frying showed great potential for drying materials, especially sewage sludge, obtaining a product with a high energy content, thereby increasing its value as a combustible.

Process systems modelling and applications in granulation: A review

Granulation is one of the fundamental operations in particulate processing and has a very ancient history and widespread use. Much fundamental particle science has occurred in the last two decades to help understand the underlying phenomena. Yet, until recently the development of granulation systems was mostly based on popular practice. The use of process systems approaches to the integrated understanding of these operations is providing improved insight into the complex nature of the processes. Improved mathematical representations, new solution techniques and the application of the models to industrial processes are yielding better designs, improved optimisation and tighter control of these systems. The parallel development of advanced instrumentation and the use of inferential approaches provide real-time access to system parameters necessary for improvements in operation. The use of advanced models to help develop real-time plant diagnostic systems provides further evidence of the utility of process system approaches to granulation processes. This paper highlights some of those aspects of granulation. (c) 2005 Elsevier Ltd. All rights reserved.

Modelling nucleation in wet granulation

Nucleation is the first stage in any granulation process where binder liquid first comes into contact with the powder. This paper investigates the nucleation process where binder liquid is added to a fine powder with a spray nozzle. The dimensionless spray flux approach of Hapgood et al. (Powder Technol. 141 (2004) 20) is extended to account for nonuniform spray patterns and allow for overlap of nuclei granules rather than spray drops. A dimensionless nuclei distribution function which describes the effects of the design and operating parameters of the nucleation process (binder spray characteristics, the nucleation area ratio between droplets and nuclei and the powder bed velocity) on the fractional surface area coverage of nuclei on a moving powder bed is developed. From this starting point, a Monte Carlo nucleation model that simulates full nuclei size distributions as a function of the design and operating parameters that were implemented in the dimensionless nuclei distribution function is developed. The nucleation model was then used to investigate the effects of the design and operating parameters on the formed nuclei size distributions and to correlate these effects to changes of the dimensionless nuclei distribution function. Model simulations also showed that it is possible to predict nuclei size distributions beyond the drop controlled nucleation regime in Hapgood's nucleation regime map. Qualitative comparison of model simulations and experimental nucleation data showed similar shapes of the nuclei size distributions. In its current form, the nucleation model can replace the nucleation term in one-dimensional population balance models describing wet granulation processes. Implementation of more sophisticated nucleation kinetics can make the model applicable to multi-dimensional population balance models.