Study on diffusion and flow of benzene, n-hexane and CCl4 in activated carbon by a differential permeation method

In this paper the diffusion and flow of carbon tetrachloride, benzene and n-hexane through a commercial activated carbon is studied by a differential permeation method. The range of pressure is covered from very low pressure to a pressure range where significant capillary condensation occurs. Helium as a non-adsorbing gas is used to determine the characteristics of the porous medium. For adsorbing gases and vapors, the motion of adsorbed molecules in small pores gives rise to a sharp increase in permeability at very low pressures. The interplay between a decreasing behavior in permeability due to the saturation of small pores with adsorbed molecules and an increasing behavior due to viscous flow in larger pores with pressure could lead to a minimum in the plot of total permeability versus pressure. This phenomenon is observed for n-hexane at 30degreesC. At relative pressure of 0.1-0.8 where the gaseous viscous flow dominates, the permeability is a linear function of pressure. Since activated carbon has a wide pore size distribution, the mobility mechanism of these adsorbed molecules is different from pore to pore. In very small pores where adsorbate molecules fill the pore the permeability decreases with an increase in pressure, while in intermediate pores the permeability of such transport increases with pressure due to the increasing build-up of layers of adsorbed molecules. For even larger pores, the transport is mostly due to diffusion and flow of free molecules, which gives rise to linear permeability with respect to pressure. (C) 2002 Elsevier Science Ltd. All rights reserved.

Surface diffusion of strong adsorbing vapours on porous carbon

Surface diffusion of strongly adsorbing hydrocarbon vapours on activated carbon was measured by using a constant molar flow method (D.D. Do, Dynamics of a semi-batch adsorber with constant molar supply rate: a method for studying adsorption rate of pure gas, Chem. Eng. Sci. 50 (1995) 549), where pure adsorbate is introduced into a semi-batch adsorber at a constant molar flow rate. The surface diffusivity was determined from the analysis of pressure response versus time, using a linear mathematical model developed earlier. To apply the linear theory over the non-linear range of the adsorption isotherm, we implement a differential increment method on the system which is initially equilibrated with some pre-determined loading. By conducting the experiments at different initial loadings, the surface diffusivity can be extracted as a function of loading. Propane, n-butane, n-hexane, benzene, and ethanol were used as diffusing adsorbate on a commercial activated carbon. It is found that the surface diffusivity of these strongly adsorbing vapours increases rapidly with loading, and the surface diffusion flux contributes significantly to the total flux and cannot be ignored. The surface diffusivity increases with temperature according to the Arrhenius law, and for the paraffins tested it decreases with the molecular weight of the adsorbate. (C) 2002 Elsevier Science Ltd. All rights reserved.

Difference equations in Banach spaces

Difference equations which discretely approximate boundary value problems for second-order ordinary differential equations are analysed. It is well known that the existence of solutions to the continuous problem does not necessarily imply existence of solutions to the discrete problem and, even if solutions to the discrete problem are guaranteed, they may be unrelated and inapplicable to the continuous problem. Analogues to theorems for the continuous problem regarding a priori bounds and existence of solutions are formulated for the discrete problem. Solutions to the discrete problem are shown to converge to solutions of the continuous problem in an aggregate sense. An example which arises in the study of the finite deflections of an elastic string under a transverse load is investigated. The earlier results are applied to show the existence of a solution; the sufficient estimates on the step size are presented. (C) 2003 Elsevier Science Ltd. All rights reserved.

A mechanistic model for carbon dioxide corrosion of mild steel in the presence of protective iron carbonate films - Part 3: Film growth model

A model of iron carbonate (FeCO3) film growth is proposed, which is an extension of the recent mechanistic model of carbon dioxide (CO2) corrosion by Nesic, et al. In the present model, the film growth occurs by precipitation of iron carbonate once saturation is exceeded. The kinetics of precipitation is dependent on temperature and local species concentrations that are calculated by solving the coupled species transport equations. Precipitation tends to build up a layer of FeCO3 on the surface of the steel and reduce the corrosion rate. On the other hand, the corrosion process induces voids under the precipitated film, thus increasing the porosity and leading to a higher corrosion rate. Depending on the environmental parameters such as temperature, pH, CO2 partial pressure, velocity, etc., the balance of the two processes can lead to a variety of outcomes. Very protective films and low corrosion rates are predicted at high pH, temperature, CO2 partial pressure, and Fe2+ ion concentration due to formation of dense protective films as expected. The model has been successfully calibrated against limited experimental data. Parametric testing of the model has been done to gain insight into the effect of various environmental parameters on iron carbonate film formation. The trends shown in the predictions agreed well with the general understanding of the CO2 corrosion process in the presence of iron carbonate films. The present model confirms that the concept of scaling tendency is a good tool for predicting the likelihood of protective iron carbonate film formation.

Selecting methods to solve multi-well master equations

There are several competing methods commonly used to solve energy grained master equations describing gas-phase reactive systems. When it comes to selecting an appropriate method for any particular problem, there is little guidance in the literature. In this paper we directly compare several variants of spectral and numerical integration methods from the point of view of computer time required to calculate the solution and the range of temperature and pressure conditions under which the methods are successful. The test case used in the comparison is an important reaction in combustion chemistry and incorporates reversible and irreversible bimolecular reaction steps as well as isomerizations between multiple unimolecular species. While the numerical integration of the ODE with a stiff ODE integrator is not the fastest method overall, it is the fastest method applicable to all conditions.

Effect of vehicle pretreatment on the flux, retention, and diffusion of topically applied penetrants in vitro

Purpose. The flux of a topically applied drug depends on the activity in the skin and the interaction between the vehicle and skin. Permeation of vehicle into the skin can alter the activity of drug and the properties of the skin barrier. The aim of this in vitro study was to separate and quantify these effects. Methods. The flux of four radiolabeled permeants (water, phenol, diflunisal, and diazepam) with log K-oct/water values from 1.4 to 4.3 was measured over 4 h through heat-separated human epidermis pretreated for 30 min with vehicles having Hildebrand solubility parameters from 7.9 to 23.4 (cal/cm(3))(1/2). Results. Enhancement was greatest after pretreatment with the more lipophilic vehicles. A synergistic enhancement was observed using binary mixtures. The flux of diazepam was not enhanced to the same extent as the other permeants, possibly because its partitioning into the epidermis is close to optimal (log K-oct 2.96). Conclusion. An analysis of the permeant remaining in the epidermis revealed that the enhancement can be the result of either increased partitioning of permeant into the epidermis or an increasing diffusivity of permeants through the epidermis.

An approach to verifying and debugging simulation models governed by ordinary differential equations: Part 1. Methodology for residual generation

For dynamic simulations to be credible, verification of the computer code must be an integral part of the modelling process. This two-part paper describes a novel approach to verification through program testing and debugging. In Part 1, a methodology is presented for detecting and isolating coding errors using back-to-back testing. Residuals are generated by comparing the output of two independent implementations, in response to identical inputs. The key feature of the methodology is that a specially modified observer is created using one of the implementations, so as to impose an error-dependent structure on these residuals. Each error can be associated with a fixed and known subspace, permitting errors to be isolated to specific equations in the code. It is shown that the geometric properties extend to multiple errors in either one of the two implementations. Copyright (C) 2003 John Wiley Sons, Ltd.

An approach to verifying and debugging simulation models governed by ordinary differential equations: Part 2. Residuals analysis and a case study

In Part 1 of this paper a methodology for back-to-back testing of simulation software was described. Residuals with error-dependent geometric properties were generated. A set of potential coding errors was enumerated, along with a corresponding set of feature matrices, which describe the geometric properties imposed on the residuals by each of the errors. In this part of the paper, an algorithm is developed to isolate the coding errors present by analysing the residuals. A set of errors is isolated when the subspace spanned by their combined feature matrices corresponds to that of the residuals. Individual feature matrices are compared to the residuals and classified as 'definite', 'possible' or 'impossible'. The status of 'possible' errors is resolved using a dynamic subset testing algorithm. To demonstrate and validate the testing methodology presented in Part 1 and the isolation algorithm presented in Part 2, a case study is presented using a model for biological wastewater treatment. Both single and simultaneous errors that are deliberately introduced into the simulation code are correctly detected and isolated. Copyright (C) 2003 John Wiley Sons, Ltd.

Surface diffusion of strongly adsorbing vapors in activated carbon by a differential permeation method

Conventional methods to determine surface diffusion of adsorbed molecules are proven to be inadequate for strongly adsorbing vapors on activated carbon. Knudsen diffusion permeability (B-k) for strongly adsorbing vapors cannot be directly estimated from that of inert gases such as helium. In this paper three models are considered to elucidate the mechanism of surface diffusion in activated carbon. The transport mechanism in all three models is a combination of Knudsen diffusion, viscous flow and surface diffusion. The collision reflection factor f (which is the fraction of molecules undergoing collision to the solid surface over reflection from the surface) of the Knudsen diffusivity is assumed to be a function of loading. It was found to be 1.79 in the limit of zero loading, and decreases as loading increases. The surface diffusion permeability increases sharply at very low pressures and then starts to decrease after it has reached a maximum (B(mum)s) at a threshold pressure. The initial rapid increase in the total permeability is mainly attributed to surface diffusion. Interestingly the B(mum)s for all adsorbates appear at the same volumetric adsorbed phase concentration, suggesting that the volume of adsorbed molecules may play an important role in the surface diffusion mechanism in activated carbon. (C) 2003 Elsevier Ltd. All rights reserved.

A new wavelet-based adaptive method for solving population balance equations

A new wavelet-based adaptive framework for solving population balance equations (PBEs) is proposed in this work. The technique is general, powerful and efficient without the need for prior assumptions about the characteristics of the processes. Because there are steeply varying number densities across a size range, a new strategy is developed to select the optimal order of resolution and the collocation points based on an interpolating wavelet transform (IWT). The proposed technique has been tested for size-independent agglomeration, agglomeration with a linear summation kernel and agglomeration with a nonlinear kernel. In all cases, the predicted and analytical particle size distributions (PSDs) are in excellent agreement. Further work on the solution of the general population balance equations with nucleation, growth and agglomeration and the solution of steady-state population balance equations will be presented in this framework. (C) 2002 Elsevier Science B.V. All rights reserved.

Hydrodynamic origin of diffusion in nanopores

We study the transport of a subcritical Lennard-Jones fluid in a cylindrical nanopore, using a combination of equilibrium and nonequilibrium as well as dual control volume grand canonical molecular dynamics methods. We show that all three techniques yield the same value of the transport coefficient for diffusely reflecting pore walls, even in the presence of viscous transport. We also demonstrate that the classical Knudsen mechanism is not manifested, and that a combination of viscous flow and momentum exchange at the pore wall governs the transport over a wide range of densities.

Analysis of multicomponent adsorption kinetics on activated carbon

An integrated mathematical model for the kinetics of multicomponent adsorption on microporous carbon was developed. Transport in this bidisperse solid is represented by balance equations in the macropore and micropore phases, in which gas-phase diffusion dominates the mass transfer in the macropores, with the phenomenological diffusivities represented by the generalized Maxwell-Stefan (GMS) formulation. Viscous flow also contributes to the macropore fluxes and is included in the MS expressions. Diffusion of the adsorbed phase controls the mass transfer in the micro ore phase, p which is also described in a similar way by the MS method. The adsorption isotherms are represented by a new heterogeneous modified vacancy solution theory formulation of adsorption, which has proved to be a robust method for adsorption on activated carbons. The model is applied to the coadsorption and codesorption of C2H6 and C3H8 on Ajax and Norit carbon, as well as the displacement on Ajax carbon. The effect of the viscous flow in the macropore phase is not significant for the cases studied. The model accurately predicts the overshoot behavior and rollup of C2H6 during coadsorption. The prediction for the heavier compound C3H8 is always satisfactory, though at higher C3H8 mole fraction, the overshoot extent of C2H6 is overpredicted, possibly due to neglect of heat effects.

Modeling the Diffusion of Energy Performance Contracting

The Portuguese Energy Policy considers the development of a commercially viable and competitive market for energy performance contracting (EPC) as a main mechanism to achieve the objectives of energy efficiency improvement. This paper proposes a study to investigate how to achieve widespread adoption of energy performance contracting by means of system dynamics modelling and simulation. To explore and gather insights on this question, a system dynamics model representing the system of the Portuguese EPC market at industry level will be created. The simulation of that model will provide a helpful basis for analysing and explaining the development of key variables, and accelerating learning on the managerial, organizational and political adaptation processes that foster the diffusion of EPC adoption. The first phase of this research project aims at identifying and analysing the key factors and critical cause-effect relations that drive the adoption of EPC. With this purpose, a qualitative content analysis on relevant documents was performed and a set of interviews was conducted. That data was analysed to capture the critical variables and its interrelation to formulate a preliminary representation of the system structure as stock and flow diagrams.