Services on Demand
Journal
Article
Indicators
 Cited by SciELO
 Access statistics
Related links
 Cited by Google
 Similars in SciELO
 Similars in Google
Share
Brazilian Journal of Chemical Engineering
Print version ISSN 01046632
Braz. J. Chem. Eng. vol.29 no.3 São Paulo July/Sept. 2012
http://dx.doi.org/10.1590/S010466322012000300017
SEPARATION PROCESSES
Effective transport properties for the pyridinegranular activated carbon adsorption system
S. A. BazRodríguez^{I}; R. OcampoPérez^{II}; J. P. RuelasLeyva^{I}; C. G. AguilarMadera^{III, * }
^{I}División de Ciencias Básicas e Ingenierías, Universidad Autónoma MetropolitanaIztapalapa, Apartado Postal 55534, México D.F. 09340, México
^{II}Centro de Investigación y Estudios de Posgrado, Facultad de Ciencias Químicas, Universidad Autónoma de San Luis Potosí, Av. Dr. M. Nava No. 6, San Luis Potosí, SLP 78210, México
^{III}División de Ciencias Básicas e Ingenierías, Universidad Autónoma MetropolitanaIztapalapa, Apartado Postal 55534, México D.F. 09340, México. Email: carlos_aguilarmadera@hotmail.com
ABSTRACT
In this work, the kinetics of pyridine adsorption onto granular activated carbon was studied from the point of view of an upscaling process by using the method of volume averaging. The pore and surface effective diffusivities were estimated by supposing simple microscale geometries (ordered media of cylinders and spheres) and those of images processed from SEM (Scanning Electron Microscopy) micrographs. In addition, as a rough estimate, the point surface diffusivity is reported. The results revealed that the upscaled diffusional model satisfactorily interpreted the concentration decay curves and the effective diffusivity was found to be an increasing function of the concentration, mainly due to the contribution of surface diffusion. In general, the diffusivity coefficients involved in the adsorption system are related through the expression molecular diffusivity = 22 ï' point surface diffusivity = 5/2 x' pore effective diffusivity = 1/12 x ' surface effective diffusivity.
Keywords: Pore and surface effective diffusion; Method of volume averaging; Granular activated carbon; Pyridine adsorption.
INTRODUCTION
Pyridine is a colorless and flammable liquid with an unpleasant odor, which is widely used in the manufacturing of fine chemicals (i.e., drugs, vitamins, dyes), agrochemicals, among others (Sabah and Celik, 2002). Due to this, it is commonly found in industrial wastewater, leading to environmental pollution problems. This has motivated several efforts to develop methods to remove pollutants from aqueous solutions such as: biological degradation (Mathur et al., 2008; Qiao and Wang, 2010), ozonation (Andreozzi, 1991), adsorption (Lataye et al., 2008; OcampoPérez et al., 2010) and others.
In adsorption systems, granular activated carbon (GAC) is one of the most efficient adsorbents used for removing organic compounds from aqueous solutions (Liu et al., 2009). GAC has been considered in the U.S. by the EPA (Environmental Protection Agency) to be one of the best available technologies for controlling several aromatic pollutants (volatile organic compounds, phenols, pesticides and herbicides) (Zhang and Wang, 2009) due to its chemical and morphological properties (MorenoCastilla et al., 1995). Moreover, it has been reported that GAC can efficiently remove pyridine from aqueous solution (Mohan et al., 2004; Lataye et al., 2008).
The global adsorption rate on porous materials like GAC involves three consecutive steps: external mass transport, intraparticle diffusion and adsorption on active sites of the solid matrix. The intraparticle diffusion can be due to pore volume diffusion, surface diffusion, or a combination of both mechanisms (LeyvaRamos and Geankoplis, 1994; OcampoPérez et al., 2010). Pore volume diffusion refers to the movement of the adsorbate due to concentration gradients in the fluidphase (i.e., molecular mechanisms), but affected by the geometry of the porous matrix. For diluted solutions, the pore volume diffusion depends only on the molecular diffusivity and the spatial distribution of phases at the microscale. The latter concept defines the geometrical properties of the porous medium such as the porosity and the tortuosity. Surface diffusion refers to the movement of the adsorbate through the solid surface and is also influenced by the distribution of phases; the main driving force is the surface concentration gradient. It has been found that surface diffusion depends on the surface concentration and the sorptive affinity between the molecules of adsorbate and adsorbent (Suzuki and Fujii, 1982; Furuya et al., 1996).
In many systems presenting a hierarchy of scales, like GAC, one of the main challenges is the theoretical description on a convenient scale of observation. Thus, for design and modeling macroscale average quantities are enough, but in order to take into account the physics at lower size scales an upscaling method is required. One standard procedure is the method of volume averaging (Whitaker, 1999) which has been used in many types of transport problems. For instance, OchoaTapia et al. (1993) carried out the theoretical upscaling of the diffusive and adsorptive transport of solute considering a linear isotherm in a homogeneous twophase porous medium and reported the associated closure problems for prediction of the effective properties. Subsequently, Guelli U. de Souza et al. (2007) incorporated the effect of bulk and surface reactions in their study of the dyeing of packed cotton threads and also reported the associated closure problems for prediction of the effective properties. Nonetheless, the influence of reactive conditions on effective properties like the effective diffusivity is unclear and some efforts have been devoted in the literature to gaining insight into this issue (Borges da Silva et al., 2007; ValdésParada and AlvarezRamírez, 2010; ValdésParada et al., 2011).
The closure problems arising theoretically in the method of volume averaging must be solved in representative domains of the microstructure in order to predict the effective transport properties. Thus, simple geometries have commonly been used (Chang, 1983, Whitaker, 1999). Even some more complicated domains from real material images have been employed (Lux et al., 2006).
In this work, the problem of adsorption of pyridine onto GAC was studied from the point of view of an upscaling process and the corresponding effective transport properties were determined. In order to model mass transport, we retook the theoretical upscaling process carried out by OchoaTapia et al. (1993) for the phenomenon on the microscale via the volume averaging method. Thus, governing equations on the macroscale are obtained that are expressed in terms of effective properties. These equations were used to interpret experimental data and to predict the pore and surface effective diffusivities. To this end, the microstructure was assumed to have simple geometries (ordered media of cylinders and spheres) or were obtained from images processed from SEM (Scanning Electron Microscopy) micrographs. Furthermore, as a rough estimate, the point surface diffusivity is reported.
UPSCALING OF MASS TRANSPORT
The system under study is presented in Figure 1. This consists of a fluid phase in which the adsorbate is diluted and of particles of the porous adsorbent media. In order to clearly identify each homogenous portion of the system (i.e., a portion with positionindependent effective properties), two different regions were defined: the ï¨region represents the homogenous fluid outside the particles of adsorbent; such particles are named the region. Since the region is composed of two phases (the solid matrix and the fluid within the pores), it is necessary to focus first on the derivation of the governing equation for mass transport in this region.
Mass Transport Description on the Microscale in the Porous Particle
The upscaling process begins by defining the mass balance for the adsorbate on the porescale (microscale). A single fluid phase (phase) is considered that completely saturates the interstices of the solid matrix (kphase), comprising the porous region (see Figure 1). In this type of system, it is appropriate to neglect the mass transport by convection within the pores as long as the associated Reynolds number is less than 1, which implies that mass transport mainly occurs by diffusive processes. Thus, the mass balance of pyridine (identified by the subscript A) inside the pores is given by:
In the phase
which is valid for dilute solutions. In the above equation represents the adsorbate concentration on the microscale, is the molecular diffusivity and is the spatial differential operator. The surface mass transport at the fluidsolid interface (k) is expressed as (Slattery, 1990, section 5.2.1),
At the k interface
where C_{As} is the surface adsorbate concentration at the k interface, is the unit normal vector pointing from the fluid phase to the solid phase, D_{As} represents the point surface diffusivity of the adsorbate and, in this case, is the surface differential operator [see Eq. (1.11) in the work of OchoaTapia et al. (1993) for its definition]. In Eq. it has been explicitly pointed out that the surface mass balance for the adsorbed pyridine encompasses accumulation, mass flux from or to the pore bulk and transport by a diffusive mechanism. As a matter of fact, the term allows quantifying the mass of pyridine adsorbed or desorbed at a certain position of the fluidsolid interface.
The problem of mass transport on the microscale is completed by setting the initial conditions for the pore and surface concentrations and the boundary conditions applying at the limits of the system. It is worth pointing out that the different concentrations handled in Eqs. and differ in the units, so has units of mass/volume and C_{As} of mass/area. In addition, D_{Ay} and D_{As} are physical properties related to the microscale and it is important to clearly distinguish from the effective transport properties involved in the upscaled models (explained later), which are related to the volume averaging scale.
Macroscale Description within the Adsorbent Medium
Effective Medium Equation
In principle, it is possible to describe completely the pyridine transport by numerically solving the microscale governing equations. However, it is a highly demanding computational task and possibly unnecessary, given that, for practical purposes, it is enough to describe the mass transport by using upscaled models (effective medium equations) which are expressed in terms of average variables and effective transport properties. Generally, such upscaled models are obtained by: 1) postulating the form of the equations, 2) applying some upscaling method to the microscale governing equations. In this work, the theoretical results of OchoaTapia et al. (1993) are used in which the method of volume averaging was applied for the upscaling process. These researchers derived the effective medium equation for solute mass transport in the porous particle (region):
where is the porosity, a_{v} is the interfacial area per unit volume, D_{eff} is the effective diffusivity tensor and is the intrinsic average concentration in the porous particle, defined as,
Here represents the portion of the averaging volume occupied by the fluidphase (see Figure 1). Notice that the form of Eq.(3) is similar to the microscale balance equation, Eq. , but with the term for accumulation modified by the adsorption equilibrium and the porosity. In addition, in Eq. (3) K_{eq} is the adsorption equilibrium constant, assuming the following linear constitutive expression:
It must be mentioned that expression is a strong assumption since it has been demonstrated that the pyridineGAC adsorption system obeys a PrausnitzRadke isothermal model (OcampoPérez et al., 2010). As a matter of fact, one can think of the coefficient K_{eq} as being the derivative of the adsorption isotherm at a given value of as is suggested in Figure 2. With this in mind, a decreasing dependency of K_{eq} with concentration is expected, which is corroborated later in the work.
Eq. (3) is valid at every point within the porous bulk, shown in Figure 1 and its validity is reduced as we approach to the boundary. Strictly, in order to use Eq. close to the particle porousfluid boundary, more detailed work about the jump boundary condition that applies at such boundary is required (ValdésParada et al., 2009); however, the development and implementation of this type of boundary conditions are beyond the scope of this paper. In addition, it should be stressed that, in obtaining Eq. (3), one of the main assumptions is that the adsorption equilibrium follows a linear behavior.
Closure
Eq. (3) is in closed form since it does not contain unknown (or empirical) parameters and it is explicitly related to the geometrical and physical microscale properties. Thus, the effective diffusivity tensor is given by (OchoaTapia et al., 1993):
where I is the identity tensor, represents the geometric place of the fluidsolid interface, f is the socalled closure vector and is a function defined as,
where H represents the mean curvature radius of the k interface (Slattery, 1990, section A5.3). In cases without adsorption, K_{eq}, D_{As}=0 and the expression valid for passive diffusion is given by (Ryan et al., 1981):
To differentiate from the effective diffusivity (D_{eff} ), D_{p} is used as the pore effective diffusivity tensor which is only a function of the molecular diffusivity and of the microstructure (quantified through the integral) commonly associated with the tortuosity. Several works have been addressed to obtain theoretical and experimental values of D_{p} for various porous materials (Currie, 1960; Ryan et al., 1981; among others).
The closure variable f satisfies the following steady boundary value problem (closure problem).
Periodicity conditions
The closure problem is solved in representative zones of the homogeneous porous medium namely unit cell. The periodicity conditions are applied at the external boundaries of the unit cell. The formal derivation of such a closure problem implies the assumptions of several lengthscale and time restrictions that must be considered during the application step. The presentation and discussion of these restrictions is vast and is beyond the scope of this work. Thus, for a detailed discussion of the validity and theoretical development of the closure problem we refer the reader to the work of OchoaTapia et al. (1993).
Note that the definition of the effective diffusivity, Eq. (6), cannot be divided in two contributions for pore and surface diffusion because the variables f and depend on the values of K_{eq} and D_{As} . Moreover, following the developments of OchoaTapia et al. (1993) it is possible to carry out the decomposition [see Eq. (4.43) therein],
In this case, the closure problem for prediction of D_{p} is obtained by setting b and g equal to zero in Eqs. (9)(11). In Eq. (14) D_{s} is the surface effective diffusivity, which is a function of the physical properties (K_{eq}, D_{As}, D_{Ay}) and the geometrical parameters (a_{v}, E_{y}) of the microscale. In this way, knowing the individuals values for D_{eff} and D_{p} , the corresponding one for D_{s} can be calculated. Using Eq. (13), the scaled model (3) can be rewritten as,
This equation presents a form similar to those models usually used for experimentally determining the values of the surface effective diffusion coefficient (OcampoPérez et al., 2010), but the substantial difference lies in the available clear physical meaning of the coefficients involved.
Solution of the Closure Problem
The domain in which the closure problem is solved must contain enough information about geometrical features of the porous material in order to accurately capture the distribution of phases on the microscale. With this idea in mind, SEM images of the GAC, filtered with an image processing software, are one of the best options. Nevertheless, the solution of the closure problem on simplified geometries like those shown in Figure 3 (geometrical models 1 and 2) is generally enough to obtain acceptable predictions of the effective properties (Ryan et al., 1981; Nozad et al., 1985) and to elucidate the effect of variables (such as the porosity and the equilibrium constant). Thus, in this work, the three geometrical models shown in Figure 3 are used and the results are compared.
OchoaTapia et al. (1993) analytically solved the closure problem in 2D and 3D Chang cells (units cells 1 and 2 shown in Figure 3) (Chang, 1983). Using the 2D cell, the longitudinal component of the effective diffusivity is given by:
where
while for the 3D cell:
with
For the pyridineGAC adsorption system, geometrical and physical properties (like the interfacial area and the molecular diffusivity) are available in the literature (OcampoPérez et al., 2010). Nevertheless, to the best of our knowledge, there are no reported values for the point coefficient D_{As} . In this regard, Eqs. (16) and allow us to move in either of two directions:

All the microscale parameters (D_{Ay}, K_{eq}, D_{As}, distribution of phases) being available, we are able to predict all the macroscale ones (D_{eff}, E_{y}, a_{v}) or,

From the known values of all the macroscale and microscale parameters except one of the physical properties, we are in a position to predict the remaining one.
Thus, it is possible from comparisons with experimental data to obtain the effective diffusivity and, with the known values of porosity, molecular diffusivity, equilibrium constant and interfacial area, and assuming the microstructure, we can infer the value of D_{As} .
Macroscale Model for the Entire System
The intraparticle mass transport governed by Eq.(3) is coupled with the corresponding mass balance equation in the fluid región (nregion shown in Figure 1) given by:
where V is the volume of the external fluid; is the intrinsic average concentration of adsorbate in the nregion, defined in a similar form to Eq.(4) ; m is the mass of adsorbent in the system; S represents the contact area per mass of adsorbent between the GAC and the external fluid ( boundary sketched in Figure 1); k_{L} is the mass transfer coefficient for the external fluid, and _{¨} represents the vector locating the boundary. In addition, it should be stressed that Eq. (20) assumes a perfectly mixed external fluid.
The governing equations for each homogeneous region are matched with the following boundary condition,
where is the unit normal vector pointing from the  to the nregion. Eq.(21) implies continuity for the mass flux through the n boundary and, therefore, surface accumulative effects (excess quantities) have been discarded. A study focused on measuring excess quantities is beyond the scope of this work, and the problematic is treated in the following way. For systems where resistance to mass transfer is greater in the external fluid, the right hand of Eq. can be simplified according to the expression (LeyvaRamos and Geankoplis, 1994; OcampoPérez et al., 2010),
Note that Eqs. (15) and (21) can be applied to whatever coordinate system. Thus, assuming spherical porous particles, the model can be rewritten as:
where R_{p} is the particle radius. In these equations gradients in the angular directions have been discarded, assuming that significant changes take place only in the radial direction. The model is complemented by the symmetry condition
and the initial conditions
where C_{A0} is the initial concentration of adsorbate in the external fluid.
METHODOLOGY
Experimental Procedure
Adsorbent
The granular activated carbon was manufactured from a bituminous carbon by Calgon, Inc. (Pittsburgh, PA) and commercially available as F400. The GAC was sieved to an average particle diameter of 1.02 mm, washed several times with deionized water, dried at 110°C for 24 h and stored under isolated conditions. The surface area, pore volume and average pore diameter of GAC were determined by the N_{2}BET method using a physisorption equipment (Micromeritics Model ASAP 2010), with S_{V} = 925 m^{2}/g, V_{p} = 0.534 cm^{3}/g, and d_{p} = 2.2 nm, respectively. The solid density was p_{s} = 2.320 g/cm^{3} determined by the helium displacement method using a helium pycnometer (Micromeritics Model Accupic 1330). The particle density and void fraction were p_{p} = 1.036 g/cm^{3} and E_{p} = 0.554, respectively. The molecular diffusivity of pyridine in aqueous solution at 25 °C is 7.7 × 10^{6} cm^{2}/s (OcampoPérez et al., 2010).
Adsorption Equilibrium Data
An Erlenmeyer flask of 500 mL was used as a batch adsorber for obtaining the experimental adsorption equilibrium data. A Nylon mesh bag with a known load of GAC was introduced into the flask and then the pyridine solution was added. The adsorption process was isothermal and continuously stirred. The experimental adsorption equilibrium data for pyridine on GAC were obtained as follows. A Nylon mesh bag with 1 g of GAC and 480 mL of a pyridine solution, with a known initial concentration at pH=10, were added to the batch adsorber. The adsorption equilibrium of pyridine on CAG at different pH and temperature values was previously studied by OcampoPerez et al. (2010) and they showed that, at pH 10, GAC presents the highest adsorption capacity to remove pyridine from aqueous solution. Therefore, this pH was chosen to study the pyridine adsorption rate on GAC.
The initial pyridine concentration was varied from 20 to 1000 mg/L. The pyridine solution was maintained in contact with the GAC particles until equilibrium was reached. The solution pH was measured periodically and kept constant by adding 0.01 and 0.1 M solutions of HCl and NaOH, as appropriate. After reaching equilibrium, the pyridine concentration in the solution was determined spectrophotometrically and the mass of pyridine adsorbed at equilibrium was calculated by making a mass balance of pyridine. In previous runs without GAC, it was demonstrated that pyridine did not adsorb on the Nylon bag.
Determination of Pyridine Concentration in Water Solution
The pyridine concentration in aqueous solution was determined by UVVisible spectroscopy. The absorbance of a pyridine solution was measured using a Shimadzu model UV160 spectrophotometer at a wavelength of 249.5 nm. Thus, the pyridine concentration of a sample was estimated with a calibration curve (absorbance vs. concentration) prepared with five standard pyridine solutions in the concentration range of 1050 mg/L.
Method to Obtain the Rate of Adsorption Data
A rotating basket batch adsorber was used to obtain the concentration profiles for pyridine adsorption on GAC. The adsorber was composed of a 1 L threeneck reaction flask and an impeller with its blades replaced with stainless steel baskets. A pyridine solution was poured into the adsorber and the GAC particles were placed in the stainless steel mesh baskets, which were attached to a shaft connected to a variable speed motor. The adsorber was partially immersed in a constant temperature water bath.
A pH=10 solution (800980 mL) was prepared by appropriately mixing 0.01N HCl and NaOH solutions and added to the adsorber. The baskets containing GAC (15 g) were subject to several rotating speed conditions (100, 150, 200 rpm) according to the experiment. The pH of the solution was controlled by appropriately adding 0.01 and 0.1N NaOH solutions. The solution and GAC were left in contact until the temperature and pH remained constant. Subsequently, the stirrer was turned off and an aliquot of a solution of known concentration of pyridine at pH = 10 was expeditiously poured into the adsorber solution to obtain the desired initial concentration. After the addition, the total volume of the solution was 1 L. Right away, the impeller motor and the timer were turned on. The solution pH was monitored with a pHmeter and adjusted as mentioned above. The total volume of NaOH solution added was always less than 2 mL and this volume represented 0.1% of the total volume. The solution was periodically sampled (5 mL) and analyzed to determine the pyridine concentration. Sampling times were 0, 1, 3, 5, 10, 15, 20, 25, 30, 40, 50, 60, 90, 120, 150, and 180 min. The total volume of the adsorber solution (1 L) was kept constant by adding 5 mL of a makeup solution immediately after sampling. The concentration of the makeup solution was the average of the initial concentration and the final concentration at equilibrium. The purpose of this addition was to replace the mass of pyridine withdrawn in the sample. Errors due to sampling and adding makeup solution were calculated by performing a mass balance and were always less than 2 % in all experiments.
Image Processing and Numerical Solution
The morphology of GAC was examined by means of a scanning electron microscope (LeicaCambridge, model S420I) using a magnification factor of 2400x. The images were processed by using the Image Processing Toolbox included in the commercial software Matlab. The procedure was as follows. Firstly, the contrast between dark and light zones was improved by translation to grayscale [from 0 (black) to 255 (white)]. Since dark zones visually indicate profundity, it was assumed that these represent the portion occupied by pores, whereas light zones are related to the solid matrix (GAC). With this idea in mind, the grayscale value of each pixel was tracked (from the black color) and accumulated to give the volumetric fraction ε_{γ} of the darker pixels. Secondly, all the pixels with a value equal or less than a threshold value were converted to white (space occupied by pores), whereas pixels with higher values to black (solid matrix). Finally, in order to enhance the definition of the contours, filters were applied. This procedure allows us to obtain a final binary image in which the ratio between white and total number of pixels equals the porosity of the GAC. Indeed, the processed images from the SEM micrographs provided the domains in which the closure problem was solved. For this purpose, the finite element solver COMSOL Multiphysics^{TM} 3.5a was used to obtain the field solution of the closure problem associated with the effective tensor D_{p}.
The effective diffusivity (D_{eff}) was determined by comparison between experimental and numerical dynamic profiles of the pyridine concentration in the solution. The numerical profiles were obtained by solving the diffusional model given in Equations (28) using COMSOL Multiphysics^{TM} 3.5a. The values reported in this work are those calculated when the experimental and numerical profiles matched within an average relative error of 0.05. The average relative error is defined as follows:
with n being the total number of experimental data for each experiment. Knowing the values of D_{eff} and D_{p}, the surface effective diffusivity was then calculated using Eq. (14). Moreover, D_{s} and D_{As} were calculated from the solution of the closure problems [Eqs. (16)(19)] in cylindrical and spherical Chang unit cells.
RESULTS AND DISCUSSION
As a first step to predict the effective properties of the pyridineGAC adsorption system, the diffusional model given by Eqs. (23)(28) was solved and fitted with the set of experimental data. The dynamical results are depicted in Figure 4, in which means the nondimensional concentration in the solution. For all the cases shown, the diffusional model interprets the experimental data reasonably well, with the following average relative errors: 0.05, 0.032, 0.0084, 0.014, and 0.0084 for = 41.4, 66.3, 80.8, 99.9 and 129.4 mg/L, respectively. As expected, the higher the equilibrium concentration, the more time necessary to reach equilibrium, though in practical terms 180 min were enough in each experiment.
The fitted diffusional model provides the effective transport properties. The effective diffusivity vs. the equilibrium pyridine concentration in the solution is presented in Figure 5, where D_{eff} is found to be an increasing function of the concentration. In this case, as stated in Eq. (14), the effective diffusivity encompasses two different contributions: the pore and the surface effective diffusivities. Since the first one is independent of concentration (i.e., only depends on the microstructure and the molecular diffusivity), the dependency of the effective diffusivity shown in Figure 5 is due to the surface contribution. This idea is in agreement with previous behavior found elsewhere (Furuya et al., 1996; OcampoPérez et al., 2010), as will be discussed later.
It should be stressed that, in the present work, three different geometrical models of the microstructure were used (see Figure 3). From the analytical results of OchoaTapia et al. (1993) (Chang's cells) and the numerical solution of the closure problem in domains derived from SEM micrographs, the pore effective diffusivities were calculated. The values of D_{p} are shown in Table 1. The results for Chang cells are significantly different than those from SEM micrographs. Since the latter is a result obtained from more "realistic" representations of the microstructure of the porous medium, it is concluded that the Chang cells are not appropriate for predicting this effective transport property for the case under study. Because Chang cells assume an isotropic medium, an analysis of the transverse and longitudinal components of the effective pore diffusivity tensor in SEM micrographs is useful to clarify the mismatch of results in Table 1.
The value reported in Table 1 for SEM micrographs is an arithmetic average of the longitudinal and transverse diffusivity values for the three processed samples shown in Figure 6; the standard deviation of that average value is 0.02035. Additionally, in Figure 6 the xx and yycomponents of the effective pore diffusivity tensor are shown. It can be noted that the longitudinal components are larger than their transverse counterparts. This suggests that the tensor D_{p} is slightly anisotropic and suggests that the pores within the GAC are, in general, oriented in a given direction. Indeed, this anisotropy is related to the laminar microstructure experimentally observed in the GAC (Marsh and RodriguezReynoso, 2006).
By subtracting the pore contribution to the effective diffusivity, the value of D_{s} is directly obtained for images processed from SEM micrographs. In Table 2, the values of D_{s} for the three geometrical models are shown as a function of the equilibrium concentration in the solution. By comparing with the D_{p} values in Table 1, it can be noted that the contribution of the surface effects is larger than the pore contribution for the effective diffusivity. It leads to similar values of D_{s} for the analyzed geometries. In particular, a rough estimate indicates that D_{s}=12D_{Ay} .
The values of k_{L} for the experimental runs are reported in Table 2; this parameter ranges from 0.009 to 0.015 cm/s and has a weak and decreasing dependency with concentration. This is consistent with the trend of values reported elsewhere (OcampoPérez et al., 2010). Finally, the fitted values of K_{eq} are also shown in Table 2. This parameter exhibits a clearly decreasing functionality with concentration and, as mentioned above, this is consistent because K_{eq} is directly related to the derivative of the nonlinear adsorption isotherm.
To conclude this work, it is opportune to take advantage of the analytical approaches given by OchoaTapia et al. (1993) [Eqs. (16)(19)] to calculate the local physical property D_{As}. Note that this is the only unknown parameter and it is required to find the roots satisfying Eqs. (16) and in order to determine D_{As}. In Figure 7 the values of the local surface diffusivity are reported as a function of the equilibrium concentration for cylindrical and spherical Chang cells. As long as the equilibrium concentration increases, the local surface diffusivity increases. For all the cases, the results with the 3D cell are larger than those for the 2D cell. In fact, the difference between the two geometries increases with concentration.
It should be stressed that the values reported in Figure 7 represent a first attempt to give a rough estimated of the local surface diffusivity. To carry this out, the simple geometries involved in the Chang cells were used. By comparison with the microstructures in the SEM micrographs in Figure 6, these are obviously poor representations of the microstructure of the GAC. However, as the effective parameter D_{s} was found to be weakly dependent on the local geometry, the same behavior can be inferred for D_{As}. It is recognized that this might be a serious assumption and a more detailed study must be done. Such a study would involve solving the closure problem given in Eqs. (9)(13) in more realistic geometries like those shown in Figure 6, which in turn represents a challenging task due to the necessity of knowing the curvature H at each point in the fluidsolid matrix interface.
CONCLUSIONS
The upscaled diffusional model given by Eqs. (28) satisfactorily interpreted the concentration decay curves of pyridine absorption onto granular activated carbon. The total effective diffusivity coefficient was found to be an increasing function of the pyridine concentration at equilibrium, mainly due to the larger contribution of the effective surface diffusion in comparison with the pore volume diffusion. In addition, a rough estimate of the point surface diffusivity was obtained by using simple geometrical models of the microstructure. In general, for the adsorption system under study, the following expression can be accepted for the different diffusivity coefficients involved:
NOMENCLATURE
REFERENCES
Andreozzi, R., Insola, A., Caprio, V. and D'Amore, M. G., Ozonation of pyridine in aqueous solution: echanistic and kinetic aspects. Water Research, 25, 655659 (1991). [ Links ]
Borges de Silva, E. A., Souza, D. P., Ulson de Souza, A. A. and Guelli U. de Souza, S. M. A., Prediction of effective diffusivity tensors for bulk diffusion with chemical reactions in porous media, Brazilian Journal of Chemical Engineering, 24, (1), 4760 (2007). [ Links ]
Chang, H. C., Effective diffusion and conduction in twophase media: a unified approach. AIChE Journal, 29, 846853 (1983). [ Links ]
Currie, J. A., Gaseous diffusion in porous media: Nonsteady state method. British Journal of Applied Physics, 11, 314320 (1960). [ Links ]
Furuya, E. G., Chang, H. T., Miura, Y., Yokomura, H., Tajima, S., Yamashita, S. and Noll, K. E., Intraparticle mass transport mechanism in activated carbon adsorption of phenols. Journal of Environmental Engineering, 122, 909916 (1996). [ Links ]
Guelli Ulson de Souza, S. M. A., Pacheco de Souza, D., Borges de Silva, E. A. and Ulson de Souza, A. A., Modelling of the dyeing process of packed cotton threads using reactive dyes. Transport in Porous Media, 68, 341363 (2007). [ Links ]
Lataye, D. H., Mishra, I. M. and Mall, I. D., Pyridine sorption from aqueous solution by rice husk ash (RHA) and granular activated carbon (GAC): Parametric, kinetic, equilibrium and thermodynamic aspects. Journal of Hazardous Materials, 154, 858870 (2008). [ Links ]
LeyvaRamos, R. and Geankoplis, C. J, Diffusion in liquidfilled pores of activated carbon. I. Pore volume diffusion. Canadian Journal of Chemical Engineering, 72, 262271 (1994). [ Links ]
Liu, G., Ma, J., Li, X. and Qin, Q., Adsorption of bisphenol A from aqueous solution onto activated carbons with different modification treatments. Journal of Hazardous Materials, 164, 12751280 (2009). [ Links ]
Lux, J., Ahmadi, A., Gobbé, C. and Delisée, C., Macroscopic thermal properties of real fibrous materials: Volume averaging method and 3D image analysis. International Journal of Heat and Mass Transfer, 49, 19581973 (2006). [ Links ]
Marsh, H. and RodriguezReynoso, F., Activated carbon, 1st Ed. Elsevier, (2006). [ Links ]
Mathur, A. K., Majumder, C.B., Chatterjee, S. and Roy, P., Biodegradation of pyridine by the new bacterial isolates S. putrefaciens and B. sphaericus. Journal of Hazardous Materials, 157, 335343 (2008). [ Links ]
Mohan, D., Singh, K. P., Sinha, S. and Gosh, D., Removal of pyridine from aqueous solution using low cost activated carbons derived from agricultural waste materials. Carbon, 42, 24092421 (2004). [ Links ]
MorenoCastilla, C., RiveraUtrilla, J., LópezRamón, M. V. and CarrascoMarín, F., Adsorption of some substituted phenols on activated carbons from a bituminous coal. Carbon, 33, 845851 (1995). [ Links ]
Nozad, I., Carbonell, R. G. and Whitaker, S., Heat conduction in multiphase systems I. Theory and experiment for twophase systems. Chemical Engineering Science, 40, 843855 (1985). [ Links ]
OcampoPérez, R., LeyvaRamos, R., AlonsoDavila, P., RiveraUtrilla, J. and SanchezPolo, M., Modeling adsorption rate of pyridine onto granular activated carbon. Chemical Engineering Journal, 165, 133141 (2010). [ Links ]
OchoaTapia, J. A., del Rio, J. A. and Whitaker, S., Bulk and surface diffusion in porous media: an application of the surfaceaveraging theorem. Chemical Engineering Science, 48, 20612082 (1993). [ Links ]
Qiao, L. and Wang, J., Microbial degradation of pyridine by Paracoccus sp. isolated from contaminated soil. Journal of Hazardous Materials, 176, 220225 (2010). [ Links ]
Ryan, D., Carbonell, R. G. and Whitaker, S., A theory of diffusion and reaction in porous media. AIChE Symposium, 71, 4662 (1981). [ Links ]
Sabah, E. and Celik, M. S., Interaction of pyridine derivatives with sepiolite. Journal of Colloid and Interface Science, 251, 3338 (2002). [ Links ]
Slattery, J. C., Interfacial Transport Phenomena. Ed. Springer, New York, (1990). [ Links ]
Suzuki, M. and Fujii, T., Concentration dependence of surface diffusion coefficient of propionic acid in activated carbon particles. AIChE Journal, 28, 380385 (1982). [ Links ]
ValdésParada, F. J., AlvarezRamírez, J., Goyeau, B. and OchoaTapia, J. A., Jump condition for diffusive and convective mass transfer between a porous medium and a fluid involving adsorption and chemical reaction. Transport in Porous Media, 78, 459476 (2009). [ Links ]
ValdésParada, F. J. and AlvarezRamírez, J., On the effective diffusivity under chemical reaction in porous media. Chemical Engineering Science, 65, 41004104 (2010). [ Links ]
ValdésParada, F. J., AguilarMadera, C. G. and AlvarezRamírez, J., On diffusion, dispersion and reaction in porous media, Chemical Engineering Science, 66, 21772190 (2011). [ Links ]
Whitaker, S., The Method of Volume Averaging. Kluwer Academic Publishers, Netherlands, (1999). [ Links ]
Zhang, C. and Wang, Y., Removal of dissolved organic matter and phthalic acid esters from landfill Leachate through a complexation  flocculation process. Waste Management, 29, 110116 (2009). [ Links ]
Submitted: July 5, 2011
Revised: December 21, 2011
Accepted: January 28, 2012
* To whom correspondence should be addressed