1. Introduction
Adsorption provides an efficient means of removing organic and inorganic contaminants from aqueous solutions [1–6]. The equilibrium and kinetic characteristics of a given adsorption system are usually determined from experimental data obtained from batch contactors. In the case of batch kinetic studies, simple surface reaction models are frequently used by investigators for data analysis. Notable examples include the pseudo first-order or Lagergren kinetic equation and the pseudo second-order kinetic equation. The popularity of these two kinetic equations is due, in part, to the simplicity of their mathematical forms which allow the equations to be linearized for parameter estimation by means of linear regression. The two pseudo reaction rate equations have however encountered criticisms on several fronts. For example, the speciousness of a linearized form of the pseudo second-order equation for parameter estimation has been discussed at some length by a number of investigators [7–9]. Methodological biases and flaws exist in the comparison of the modeling abilities of these two kinetic equations [10–12]. From a theoretical viewpoint, they should be considered as being purely phenomenological with no sound physicochemical basis [13]. The physical significance of their rate coefficients has been called into question [14–16]. Also, they are merely curve fitting tools and should not be used to identify the actual uptake mechanism [17, 18]. Two recent reviews [19, 20] as well as a recent textbook on adsorption [21] have summarized the various shortcomings of the two pseudo reaction rate equations.
For porous adsorbent particles the rate of adsorption is often controlled by intraparticle diffusion [22, 23]. Thus, a more instructive approach in the modeling of adsorption kinetics is to use physically reasonable rate models based on the kinetics of mass transfer processes. Rate models accounting for an external film mass transfer resistance and one or more intraparticle diffusion mechanisms allow a realistic mathematical description of the kinetics of many adsorption systems. Because intraparticle diffusion coefficients are an intrinsic property of a given adsorption system, they are independent of the mode of operation. For this reason, diffusion coefficients determined from a batch contactor may be used to model the dynamic behavior of a fixed bed adsorber. In general, two distinct intraparticle transport models, pore diffusion and surface/solid diffusion, are commonly used to describe the kinetics of contaminant adsorption by porous particles. Both diffusion models have been used by investigators to fit batch kinetic data with equal accuracy. More sophisticated models with pore and surface diffusion acting in parallel have also been studied by some practitioners [24–26]. From a practical standpoint, it is more convenient to use solid-phase diffusion models based on a single intraparticle diffusion mechanism to model adsorption kinetics because they are computationally simpler.
The present paper deals with an important aspect of pore diffusion: determination of pore diffusion coefficients from experimental data. Batch contactors provide a simple and convenient way to acquire kinetic data. By fitting a pore diffusion model to such data, the pore diffusion coefficient of the model can be determined. However, for nonlinear systems the data fitting process is often complicated by the need to solve the model equations numerically [27, 28]. As an alternative, the present article explores the feasibility of using a simple kinetic model to extract pore diffusivities from batch kinetic data without resorting to numerical methods. The approximation is based on the use of the classic Langmuir kinetics model – which reduces to the Langmuir isotherm model at large values of time – to describe the rate of uptake. Although the Langmuir kinetics model is based on a surface reaction mechanism, its rate coefficient, unlike those of the pseudo first-order and pseudo second-order kinetic equations, can be converted to the pertinent pore diffusion coefficient [29–31].
For batch adsorption systems, it is possible to derive from the Langmuir kinetics model an analytical expression in the form of an algebraic equation. Fitting the algebraic equation to batch kinetic data yields the rate coefficient of the Langmuir kinetics model, from which the corresponding pore diffusion coefficient may be computed. The data fitting process can be easily accomplished using standard spreadsheet programs. Using data taken from the environmental adsorption literature, the potential of this simple modeling approach was evaluated here on its ability to provide accurate pore diffusion coefficients of fluoride ions. Elevated levels of fluoride in drinking water are a major public health concern in many parts of the world. Adsorption has been recognized as an effective technique for the removal of fluoride from water sources [32–34].
2. Theory
2.1. Film-Pore Diffusion Model
In this model the transport of a solute within a spherical particle is described by pore diffusion and is given by the following equation and initial and boundary conditions [22, 23]:
In these equations, c and C are respectively the pore fluid and bulk solution concentrations at time t, q is the adsorbed phase concentration at time t, ɛp is the particle porosity, ρp is the particle density, Dp is the pore diffusion coefficient, kf is the external mass transfer coefficient, r is the particle radial coordinate, and Rp is the particle radius. Eq. (1(c)) describes mass transfer through a boundary film surrounding an individual particle.
where m is the mass of adsorbent, V is the volume of solution, q̄ is the solute concentration in the adsorbent particle averaged over the particle volume, and Co is the initial concentration. The relationship between adsorbed phase and pore fluid concentrations is defined by the following Langmuir isotherm model [35]:
where qm and b are respectively the Langmuir saturation capacity and Langmuir constant. These equations are nonlinear and must, in general, be solved numerically.
2.2. Langmuir Kinetics Modeling Approach
According to the classic Langmuir kinetics model [36, 37], the interaction between a solute and an adsorption site is described by the following reaction rate equation:
where k1 is a second-order forward rate coefficient and k2 is a first-order reverse rate coefficient. At equilibrium the derivative dq̄/dt approaches zero and Eq. (4) reduces to the Langmuir isotherm model given by Eq. (3) with b = k1/k2. Eq. (4) may be rewritten as
where
For an experimental kinetic test, the parameters Co, m, and V are known while the parameters α and β may be computed with known Langmuir isotherm parameters qm and b, and so the rate coefficient k1 is the only fitting parameter which may be obtained by matching Eq. (6) to the experiment’s data. Note that Eq. (6) may be rearranged in the following manner to allow estimation of k1 by means of linear regression:
In adsorption studies k1 is often regarded as a “lumped” kinetic rate constant which includes all mass transfer resistances (film diffusion, surface diffusion, pore diffusion, and surface reaction resistance). Here, we restrict our interest to the film mass transfer and pore diffusion resistances which can be combined in the following manner [29–31]:
where,
2.3. Goodness-of-Fit Measure
The following coefficient of determination (COD) gives a good idea of the overall goodness-of-fit of the Langmuir kinetics model to measured data:
where
is the model prediction for observation j,
is the experimental data for observation j, and z̄exp is the mean of experimental data. A COD score of unity indicates a perfect fit.
3. Results and Discussion
Intraparticle diffusion coefficients (pore and surface) serve as important input parameters to solid-phase mass transfer models which are often used in fixed bed adsorber design and analysis. In addition, these diffusion coefficients may be used to compare the kinetic properties of different adsorbents. For a given contaminant removal process, it is important to select an adsorbent which exhibits not only optimal adsorption capacity but also fast uptake kinetics. As such, the subject matter of this work is of considerable practical interest.
Three literature batch data sets obtained for the uptake of fluoride by porous adsorbents have been used in this work to evaluate the performance of the proposed Langmuir kinetics modeling approach. In the selected literature reports the pore diffusion coefficient for each data set has been obtained by fitting either a numerical film-pore diffusion model or a numerical shrinking core model to the kinetic data.
An attractive feature of the Langmuir kinetics model is its ability to account for the effects of isotherm nonlinearity. As noted previously, at equilibrium the Langmuir kinetics model reduces to the Langmuir isotherm model, so the effects of isotherm nonlinearity are inherently included within the rate model. Since many experimental adsorption isotherms follow the Langmuir isotherm model, the Langmuir kinetics model can thus be used to describe the kinetics as well as the equilibrium characteristics of these systems. For rate models with a nonlinear isotherm, a numerical solution is in general required. The Langmuir kinetics model is an exception in this respect since it can be solved analytically. However, the Langmuir kinetics model cannot extract pore diffusion coefficients directly from experimental data; they are instead computed from Eq. (8).
3.1. Case 1: Adsorption of Fluoride by Bone Char
Nigri et al. [41] measured the uptake of fluoride by a commercial bone char adsorbent as a function of time using Erlenmeyer flasks containing fluoride solutions and bone char which were placed on a thermostatic orbital shaker. The resultant kinetic data set is shown in Fig. 1. The authors extracted Dp from the Fig. 1 data by means of a pore diffusion model which was solved numerically. The numerical solution is based on the explicit finite difference method of lines. The optimal value of Dp so obtained is 7.15 × 10−6 cm2/s. In their modeling approach, Nigri et al. [41] assumed negligible resistance to external film mass transfer. As a result, the pore diffusion model defined by Eq. (1) reduces to their pore diffusion model by omitting the external mass transfer mechanism.
In the Langmuir kinetics modeling approach, determination of Dp is based on two distinct steps: estimation of k1 by fitting Eq. (6) to batch kinetic data followed by computation of Dp from Eq. (8) using the resultant k1 value. To extract k1 from the Fig. 1 data via Eq. (6), one requires values of the Langmuir isotherm parameters qm and b and two system parameters Co and m/V. These parameter values, given in the paper by Nigri et al. [41], are qm = 6.7 mg/g, b = 800 cm3/g, Co = 1 × 10−2 g/cm3, and m/V = 2 × 10−3 g/cm3.
Since the Langmuir isotherm parameters serve as input parameters to Eq. (6), it is important to check whether the Fig. 1 kinetic data are consistent with prediction of the Langmuir isotherm model. This consistency will ensure a good fit of Eq. (6) to the data points of Fig. 1 collected at large values of time when apparent equilibrium was established. The equilibrium concentration at large values of time, Ce, for a batch adsorber that follows the Langmuir isotherm may be computed from the following expression:
where
The value of Ce for the batch kinetic experiment of Nigri et al. [41] computed from Eq. (10) is 1.9 × 10−3 mg/cm3, which agrees with the observed Ce value of approximately 1.8 × 10−3 mg/cm3 (indicated by the last data point of Fig. 1). Thus, it is evident that the Langmuir parameters qm and b extracted from independent equilibrium experiments are able to describe the equilibrium state of the batch kinetic experiment with good accuracy.
Eq. (6) was fitted to the Fig. 1 data by a nonlinear regression procedure using the qm, b, Co, and m/V values given above, yielding k1 = 4.38 × 10−3 cm3/mg s and a COD score of 0.986. As seen in Fig. 1, the entire data set is fitted well by Eq. (6). Having determined k1, it is straightforward to compute Dp from Eq. (8), which has several other input parameters. Values of qm, b, and Co are listed above while values of the adsorbent radius Rp and adsorbent density ρp, as reported by Nigri et al. [41], are respectively 0.07 cm (mean) and 0.96 g/cm3. The remaining parameter of Eq. (8), the external mass transfer coefficient kf, is not available since Nigri et al. [41] assumed negligible external mass transfer in their modeling approach.
Because the initial uptake of solute is controlled only by the external mass transfer resistance, kf may be estimated from experimental values of C/Co obtained during the initial uptake period. The fluid phase concentration change during the initial period is given by [42, 43]:
The preceding equation implies that a plot of C/Co against t for small times gives a straight line with slope (3m/RpρpV)kf. Therefore, kf can be calculated from such a slope. We calculated kf as 1.01 × 10−3 cm/s by fitting Eq. (11) to the first two data points of Fig. 1.
Knowing all the relevant parameters, the pore diffusion coefficient can be computed from Eq. (8): Dp = 8.03 × 10−6 cm2/s. This Dp value obtained by the Langmuir kinetics modeling approach and the Dp value determined by Nigri et al. [41] differ only by 12%, as shown in Table 1. The good agreement confirms the validity of the proposed modeling approach.
3.2. Case 2: Adsorption of Fluoride by Bone Char
The uptake of fluoride by a commercial bone char adsorbent has also been investigated by Leyva-Ramos et al. [44] who employed a rotating basket batch adsorber to obtain kinetic data for a series of experimental runs with varying system and operational conditions. We will focus on a particular kinetic experiment (run 13 in Fig. 6 of their paper). The adsorbent size of this run is similar to that of case 1. The authors fitted the film-pore diffusion model given by Eq. (1) to the chosen kinetic data set, shown in Fig. 2, obtaining Dp = 3.31 × 10−6 cm2/s. The partial and ordinary differential equations of the film-pore diffusion model (Eqs. (1) and (2)) were numerically solved by using a specialized software package. The external mass transfer coefficient kf – required by the fitting procedure as one of the model inputs – was estimated from data points in the small time region which yielded kf = 3.07 × 10−3 cm/s. Other relevant system and operational parameters are as follows: Co = 5.66 × 10−3 g/cm3, m/V = 1 × 10−3 g/cm3, Rp = 0.065 cm, and ρp = 0.96 g/cm3. Note that the ρp value is taken from case 1 since it is not given in the paper by Leyva-Ramos et al. [44].
As in case 1, one should verify the ability of the Langmuir isotherm to predict the apparent equilibrium state of the kinetic data set of Fig. 2 before fitting Eq. (6) to the kinetic data. The final data point of Fig. 2 gives an experimental Ce value of approximately 2.4 × 10−3 mg/cm3. According to Leyva-Ramos et al. [44], the Langmuir isotherm parameters qm and b extracted from batch equilibrium data are 5.89 mg/g and 450 cm3/g, respectively. The theoretical Ce value calculated from Eq. (10) using these two isotherm parameters is 2.5 × 10−3 mg/cm3, confirming the validity of the Langmuir isotherm parameters. Knowing Co, m/V, qm, and b, the unknown rate coefficient k1 can be obtained by fitting Eq. (6) to the Fig. 2 data. The best-fit value of k1 is 2.19 × 10−3 cm3/mg s. A comparison between experimental and calculated curves is shown in Fig. 2 from which it may be seen that the Langmuir kinetics model provides a satisfactory fit of the experimental data (COD = 0.961).
We calculated kf as 6.17 × 10−4 cm/s by fitting Eq. (11) to the first two data points of Fig. 2. Note that our kf value is noticeably smaller than the kf value obtained by Leyva-Ramos et al. [44]. The difference may be due to the use of different initial data points to extract kf from the small time region. We were not able to extract accurate initial data points from the original figure because it contains several cluttered data points in the small time region. Also, as noted above we have used an assumed value for the adsorbent density ρp in the calculation of kf from Eq. (11). Knowing k1, kf, qm, b, Co, Rp, and ρp, the pore diffusion coefficient Dp can be computed from Eq. (8) as 3.04 × 10−6 cm2/s. The Dp value obtained here is in excellent agreement with that determined by Leyva-Ramos et al. [44], showing a deviation of 8% (see Table 1).
It is interesting to note that the literature Dp value for case 1 (7.15 × 10−6 cm2/s) is about two times larger than the literature Dp value for case 2 (3.31 × 10−6 cm2/s). This difference is somewhat unexpected given that the textural properties of the two commercial bone char adsorbents are quite similar, as shown in Table 2. In principle, adsorbents with similar textural properties should exhibit similar pore diffusivities.
3.3. Case 3: Adsorption of Fluoride by Laterite
Batch experiments were conducted by Maiti et al. [45] to obtain kinetic data of fluoride adsorption in a chemically treated laterite adsorbent. Three kinetic data sets obtained at different temperatures are available in their paper. The kinetic data set measured at 15 °C is presented here for analysis (Fig. 3). The experimental conditions are Co = 1 × 10−2 g/cm3 and m/V = 5 × 10−4 g/cm3. The properties of the adsorbent are Rp = 0.01 cm and ρp = 1.08 g/cm3. The authors adopted a shrinking core model to describe pore diffusion and external mass transfer. The advantage of this model is that it can be formulated in terms of ordinary differential equations, avoiding the need to solve the partial differential equation given by Eq. (1). The shrinking core model was numerically solved and fitted to the kinetic data of Fig. 3, yielding Dp = 2.3 × 10−6 cm2/s and kf = 8.3 × 10−2 cm/s.
By fitting Eq. (3) to independent equilibrium data, Maiti et al. [45] obtained qm = 36.31 mg/g and b = 230 cm3/mg. This maximum adsorption capacity is about five to six times larger than the qm values of the two bone char adsorbents analyzed in the previous two cases. Using known values of Co, m/V, qm, and b, Eq. (6) was fitted to the Fig. 3 data to determine k1. The resultant k1 is 6.54 × 10−3 cm3/mg s. As seen in Fig. 3, the calculated curve of the Langmuir kinetics model with a COD score of 0.957 is capable of fitting the data reasonably well. The external mass transfer coefficient kf was again determined by fitting Eq. (11) to the first two data points of Fig. 3. Knowing k1, kf, qm, b, Co, Rp, and ρp, the pore diffusion coefficient Dp can be calculated from Eq. (8) as 1.81 × 10−6 cm2/s. The Dp value obtained by our approach is consistent with that obtained by Maiti et al. [45]; they differ by 21% (see Table 1). It should be noted that the shrinking core model used by Maiti et al. [45] to estimate their Dp value is based on a number of assumptions and simplifications. As a result, their Dp value could differ from the Dp value of the formal film-pore diffusion model given by Eq. (1), to which our Dp should be compared.
Fig. 3 reveals that the fluoride-laterite adsorption system needed about 10 h to reach apparent equilibrium. By contrast, the contact time needed by the two fluoride-bone char adsorption systems to reach apparent equilibrium was more than 50 h, as may be seen in Figs. 1 and 2. This marked difference is due to the fact that the particle size of the laterite adsorbent (0.2 mm) is very much smaller than those of the bone char adsorbents (Table 2). This particle size effect implies that the intraparticle diffusion mechanism played a dominant role in the overall adsorption rates of fluoride in these porous adsorbents.
3.4. Errors in Fitted Dp Values
On the whole the parameter estimation method presented in this work has proven to be satisfactory, providing a simple and convenient way to extract Dp from batch kinetic data without resorting to numerical methods. The fitted Dp values differ from literature values by 8–21% (Table 1). These deviations are in line with errors reported in other similar studies. Differences in the range 3–23% have been reported by Chu [40] who used the Langmuir kinetics modeling approach to extract surface diffusion coefficients from batch kinetic data. Yao and Chen [46] developed a simplified method based on algebraic equations for estimating surface diffusion coefficients from batch kinetic data and reported errors in the range 11–24%. For the uptake of fluoride by bone char in batch adsorbers, Leyva-Ramos et al. [44] demonstrated that concentration decay curves predicted with best-fit Dp values are very similar to those predicted with non-optimal Dp values that differ from the best-fit values by 10–17% (see Fig. 4 of their paper). The magnitude of the errors in Dp reported here is therefore deemed tolerable.
4. Conclusions
From our evaluation of the Langmuir kinetics modeling approach considered here, we conclude that it affords a simple and convenient way to extract pore diffusion coefficients from batch kinetic data and, with limited computational effort, gives acceptable quantitative results. For the three case studies dealing with the uptake of fluoride by bone char and laterite adsorbents, the fitted pore diffusion coefficients are in good agreement with literature results. The Langmuir kinetics modeling approach requires accurate Langmuir isotherm parameters as well as external mass transfer coefficients which serve as important input parameters. The major drawback of this simple method is that it is restricted to adsorption systems in which the adsorption equilibrium behavior follows the Langmuir isotherm. Additional studies using adsorption systems dealing with other water contaminants such as heavy metals and dyes will be required to expand the evaluation presented here and to confirm the utility of the simple method as a rapid and reliable modeling tool for estimating pore diffusion coefficients from batch kinetic data.