Abstract

The uptake data of CO2 adsorption at 273 K by coconut shell chars prepared at various carbonization temperatures from 250 to 550°C were used for characterizing pore texture of chars as well as the analysis of CO2 adsorption kinetics. The equilibrium isotherms were used to determine the porous texture of chars, employing the DR equation and GCMC simulation. It was found that all the test chars contain micropores of a size range from 0.8 to 2.2 nm with the pore size distribution becoming wider for char prepared at a higher carbonization temperature. Porous properties of chars, including surface area, total pore volume, and the average pore size, appear to increase with an increasing carbonization temperature. The analysis of CO2 uptake during the transient measurement of isotherms revealed that the kinetics of CO2 adsorption is governed by the internal diffusional transport of the adsorptive molecules. The effective pore diffusivity characterizing this transport process increases with increasing CO2 loading and passes through a maximum at a certain loading. This maximum pore diffusivity shifts to a higher value as the carbonization temperature is increased. A semiempirical equation was developed to correlate the effective pore diffusivity of CO2 with the equilibrium adsorption loading and its predictive capability is satisfactory.

1. Introduction

Pyrolysis is a simple and attractive thermal process for converting biomass into various fuel sources, including low calorific-value gases, liquid biooil, and solid carbon char, as well as a variety of useful chemical species [14]. The solid char in particular can be alternatively used as a primary low-cost adsorbent [5, 6] or further activated to produce a highly porous activated carbon which finds many applications in the separation and purification processes [7, 8]. The yield, physicochemical properties, and pore characteristics of the char product are normally determined as part of the evaluation of the pyrolysis process. In general, the nitrogen gas adsorption is the most widely used method for the characterization of porous adsorbents [9, 10]. However, this technique is not so practical and may lead to erroneous results when it is applied for biomass chars; for example, the decrease of the adsorbed amount with increasing pressure could be observed or the time to reach adsorption equilibrium is extremely long (the so-called activated diffusion effect) [11]. This is partly explained by the low diffusion rate of the N2 probe molecule at 77 K inside the narrow pores of the poorly connected pore structure of char. To solve this problem, carbon dioxide adsorption at 273 K is a preferred choice [12] because the higher adsorption temperature of CO2 enables the adsorptive molecules to enter more quickly into the narrow micropores of the char. Apart from the required data on equilibrium isotherms for characterizing porous properties of chars, information on the transport of carbon dioxide inside the internal pores is also important for the design of a thermal conversion system [13]. In addition, the diffusion of carbon dioxide in chars with different pore structures can have a significant influence on the reactivity of char gasification in the production of syngas and activated carbon [14, 15].

This work presents a standard analysis of CO2 uptake data derived from gravimetric adsorption equipment for the simultaneous determination of both pore texture of coconut shell chars and the kinetics of CO2 adsorption. Chars prepared from the carbonization of coconut shell in an inert atmosphere of nitrogen at temperatures over the range from 250 to 550°C were investigated.

2. Materials and Methods

2.1. Raw Materials

Coconut shell was used as a raw material in this work. The as-received coconut shell was first dried in an electric oven at 110°C overnight and the dry sample was crushed and sieved to obtain a screen size fraction of 12 × 16 mesh, with the average particle size of 1.44 mm. The prepared sample was kept for proximate, ultimate, and thermal analyses, as well as for char preparation.

2.2. Coconut Shell Characterization

The coconut shell sample was analyzed for the proximate analysis on a dry basis, namely, the volatile matter, the fixed carbon content, and the ash content, by employing a thermogravimetric analyzer (TGA) (model SDT 2960, TA Instruments, USA) with the analysis procedure being reported elsewhere [14]. The elemental or ultimate analysis of the coconut shell was analyzed by a CHNS analyzer (model PE 2400 Series II, Perkin Elmer, USA).

Thermal decomposition behavior of the coconut shell was determined using the same thermogravimetric analyzer by heating the sample under a constant flow of nitrogen gas from room temperature to the final temperature of 750°C and the sample weight loss was monitored.

2.3. Char Preparation

About 30 g of coconut shell sample was heated in a horizontal tube furnace (model CTF 12/75/700, Carbolite, UK) under a constant flow of N2 at the rate of 100 cm3/min. In all runs, the heating temperature was programmed to increase from room temperature to the desired carbonization temperature at the heating rate of 5°C/min and held at this temperature for 2 h. After that, the furnace was turned off and the char was cooled inside the furnace to room temperature under the flow of N2. The carbonization temperature studied was varied in the range from 250 to 550°C.

2.4. Char Characterization

The derived char products were characterized for the proximate analysis and their images of surface morphology were ascertained with a scanning electron microscopy instrument (model JEOL JSM-6400, JEOL Ltd., Japan). The true densities of the prepared chars were measured with a helium pycnometer (model AccuPyc 1330, Micromeritics, USA).

2.5. Carbon Dioxide Adsorption

Investigation of CO2 adsorption by solid chars from coconut shell was performed in an Intelligent Gravimetric Analyzer (IGA) (Hiden Analytical, UK). Initially, the char sample weighing about 0.2 g was outgassed at 300°C for 10 h at a vacuum pressure of 75 mmHg. Next, the adsorption of CO2 was carried out at 273 K by dosing CO2 at a specified pressure and the increasing amount of CO2 adsorbed was recorded as a function of time until the equilibrium was established. The measurement was then repeated for a number of increasing pressures up to the maximum pressure of 500 kPa. The uptake data obtained during the transient period and at equilibrium were subsequently used for the analyses of CO2 adsorption kinetics and pore texture of the chars, respectively.

3. Results and Discussion

3.1. Proximate and Ultimate Analyses of Coconut Shell

Table 1 presents the proximate and ultimate analyses of the raw coconut shell used in this work. The results indicate that coconut shell contains a high content of volatile matter with low compositions of ash and sulfur, and carbon and oxygen are the major elements. The fixed carbon content of 19.07% is comparable to those of some biomass materials, for example, longan seed (19.61%) [16], oil palm fiber (18.23%) [17], palm kernel cake (16.74%) [18], palm shell (19.20%) [19], and eucalyptus wood (18.30%) [20].

3.2. TGA and DTG Analysis of Coconut Shell

Nonisothermal thermogravimetric analysis was used to determine the thermal degradation behavior of coconut shell and the results are presented in Figure 1, as weight percent remaining (TG curve) and its corresponding first derivative (DTG curve). The small first peak of the DTG curve which occurs at a temperature below 120°C should result from the liberation of inherent moisture, giving the weight loss of about 2%. It is then followed by two distinct peaks where the main devolatilization occurs, with the first peak appearing in the range from 200 to 315°C and the second peak from 315 to 385°C. The maximum decomposition rates occur at 265°C and 350°C for the first peak and the second peak, respectively. Generally, the thermal decomposition of biomass results from the decomposition of its lignocellulosic components, including hemicellulose, cellulose, and lignin. Lignin usually decomposes first at a low temperature around 160°C and continues over a wide temperature range up to 900°C [21]. Hemicellulose and cellulose decompose over the temperature range from 200 to 380°C and 250 to 400°C, respectively [22]. Thus, it could be inferred that the appearance of the first peak of DTG curve is attributed mainly to the decomposition of hemicellulose component, giving the weight loss of about 28.5%, followed by the decomposition of cellulose showing the second peak with additional weight loss of 26.6%.

Surface morphology, proximate analysis, and true density were determined for the char products derived from coconut shell.

3.3. Surface Morphology of Chars

Figure 2 shows the SEM images of the raw coconut shell and those of the derived chars prepared at different carbonization temperatures of 250, 350, 450, and 550°C. It appears that the coconut shell surface consists of a stack of parallel solid flakes with relatively smooth surfaces. As the carbonization temperature is increased, the char surface shows an increase in the surface roughness. This is probably caused by the increased rupture of the biomass structure due to the release of more volatile species as the heating temperature is increased.

3.4. Proximate Analysis of Chars

The effect of the carbonization temperature on chemical compositions of solid chars is displayed in Figure 3. The results show that the volatile content decreases continuously from 68.35% to 15.50% as the carbonization temperature is increased from 250 to 550°C which results from the release of noncarbon elements such as oxygen, nitrogen, hydrogen, and sulfur as volatile gaseous products [23]. As a consequence, the fixed carbon content is increased accordingly from 31.41% to 81.1%. Since the actual mass of ash remains constant over the range of carbonization temperatures studied, the decrease of char yield with an increasing carbonization temperature would cause the ash composition to increase from 0.23% to 1.16%, as shown in Figure 3.

3.5. Char Density

Chars prepared at carbonization temperatures of 250, 350, 450, and 550°C possess the true densities of 1.19, 1.28, 1.39, and 1.51 g/cm3, respectively. Since the yield of char decreases with the increasing carbonization temperature, the increase in the true density of char indicates a certain degree of particle shrinkage. A simple calculation based on char yield and true density data shows that the average spherical diameter of a char particle decreases from 1.16 mm to 0.92 mm (a drop of about 20%) when the carbonization temperature is increased from 250 to 550°C.

3.6. Adsorption Isotherms of CO2 by Coconut Shell Chars

Figure 4 shows isotherms of CO2 adsorption at 273 K by coconut shell chars prepared at various carbonization temperatures from 250 to 550°C. Over the applied pressure up to 500 kPa (5 bars), the isotherms exhibit the initial part of Type II isotherm according to the IUPAC classification [24]. With respect to the pressure increasing, two regions of the isotherm curve can be observed, that is, a sharp increase of the amount adsorbed over the initial low pressure range from zero to 10.5 kPa ( = 0.003, where = saturation vapor pressure of CO2 at 273 K = 3500 kPa), followed by a continuous increase in the adsorbed amounts with a gradually decreasing slope of the isotherm curve.

The characteristic shape of the Type II isotherm curve suggests that there is a distribution of pore sizes even for chars with a relatively small pore volume. The adsorption of CO2 at low pressure range possibly involves the adsorption in small narrow micropores of a few molecular dimensions of the adsorptive molecules, known as primary micropore filling [25]. Our simulation study of CO2 adsorption in porous carbons using a slit-pore model [26] has revealed that the adsorption by the primary pore filling effect occurred for pore width less than 0.85 nm, in comparison with the collision diameter of CO2 of 0.303 nm. Carbon dioxide adsorption at pressures higher than 10.5 kPa occurs in wider supermicropores (1.4 to 2.0 nm) by quasi-multilayer adsorption due to adsorbate-adsorbate interaction, known as the secondary micropore filling mechanism [25]. It is also interesting to note that the Type II isotherm found for CO2 adsorption in coconut shell chars was also observed by a number of investigators for CO2 adsorption in activated carbons having much larger internal surface area [2730]. This indicates that the basic interaction forces between CO2 and carbon surfaces in solid chars and activated carbons are fundamentally similar, involving the attraction-repulsion dispersive forces and Coulombic forces induced by the quadrupole moments of CO2 molecules [31].

An attempt was also made in this work to fit the experimental isotherm data with available adsorption isotherm models. Two popularly used empirical isotherm equations were tested, namely, the two-parameter Freundlich equation and the three-parameter Sips equation [32]. The Freundlich isotherm equation was developed based on the assumption that the surface is heterogeneous; that is, the adsorbent-adsorbate interaction energy is distributed and the surface is grouped into various patches with each patch having its own adsorption energy.

The Freundlich isotherm has the following form:where is amount of CO2 adsorbed, mmol/g, is pressure, kPa, is affinity constant, and is surface heterogeneity constant.

This isotherm equation has no Henry law (linear) behavior at low pressure and no saturation limit at high pressure.

Sips equation takes the formwhere and are amount adsorbed and monolayer adsorbed capacity in mmol/g, respectively, is pressure, kPa, is affinity constant, and is heterogeneity constant.

This isotherm equation combines the characteristics of the Langmuir and Freundlich equations, thus showing a finite adsorption limit when the pressure is sufficiently high.

Table 2 shows the best fitted parameters of the Freundlich and Sips isotherm equations for CO2 adsorption at 273 K by coconut shell chars prepared at various carbonization temperatures. As seen from Figure 4 and the values of regression coefficient () in Table 2, the isotherm data are very well described by both adsorption models. It is interesting to observe that the surface heterogeneity constants, and , tend to increase with increasing carbonization temperatures. This result is possibly attributed to the increasing number of surface active sites containing unpaired electrons caused by the increasing degree of bond rupture during the devolatilization of chars [33]. The adsorbed capacity constants, and , also increase with the increase of carbonization temperature, indicating the consequent increase of available surface area for CO2 adsorption. However, the affinity constant, , of the Sips equation shows a decrease as the carbonization temperature is increased.

3.7. Pore Characterization of Coconut Shell Chars

Gas adsorption of a single adsorptive molecule is still the most widely used method for the characterization of porous solids. The required information is the surface area and pore size distribution which can be achieved by the application of adsorption theories to the derived isotherm data. For microporous materials like solid chars or activated carbons, the micropore volume and surface area can be assessed, for example, by the well-known equations such as the Dubinin-Radushkevich (DR) equation and the Dubinin-Astakhov (DA) equation or by employing the techniques of -plot and -plot [10]. The pore size distribution can be determined by the MP method [34] and the Horvath-Kawazoe (HK) method [35] or by more modern techniques of the density functional theory (DFT) [36] or the grand canonical Monte Carlo simulation (GCMC) [37]. In this work only the DR equation and GCMC simulation were used to determine the porous properties of coconut shell chars derived at different carbonization temperatures based on CO2 adsorption isotherms at 273 K.

The DR equation was developed based on the potential theory of Polanyi [38] to account for the adsorption in microporous solids by the micropore filling mechanism. The cumulated volume of the adsorbed phase is a function of the adsorption potential and this relationship is called the characteristic curve. To obtain the isotherm equation, Dubinin [39] has proposed the specific empirical functional form of the characteristic curve for microporous adsorbents such as activated carbon. The final form of DR equation readsor in linear form,where is the amount adsorbed at pressure , mmol/g adsorbent, is the amount adsorbed at the saturation pressure , mmol/g, is gas constant = 8.314 × 10−3 kJ/mol-K, is adsorption temperature, K, is characteristic energy, kJ/mol, is scaling factor that brings the characteristic curves for different adsorbates into a single curve, taken as 0.35 for carbon adsorbents.

The plotting of (4) from isotherm data should yield a straight line and the values of and are obtained from the slope and intercept, respectively. The micropore volume () of the adsorbent (coconut shell chars) is obtained by converting into the volume of adsorbate in liquid state (liquid density of CO2 is 0.0207 mol/cm3) and the micropore area () in m2/g can be calculated from by the equationwhere is the micropore adsorbed capacity (mol/g), is the cross-sectional area of CO2 (2.1 × 10−9 m2/molecule at 273 K), and is the Avogadro number (6.02 × 1023 molecules/mol).

The calculation method used for the characterization of porous chars by the GCMC simulation is that proposed by Herrera et al. [40]. The set of simulated isotherms are generated using the technique of grand canonical Monte Carlo simulation for carbon dioxide adsorption at 273 K, based on the slit-pore model (parallel graphitic surfaces) with pore width varying from 0.55 to 3 nm. The experimental isotherms and those of simulated results are then matched by using an optimization function of the MATLAB code to arrive at a pore volume for each pore size as well as its surface area. The pore size distribution, the average pore size (), the cumulated surface area, and the total pore volume are then calculated.

Figure 5 shows the DR plots for CO2 adsorption at 273 K by coconut shell chars prepared at various carbonization temperatures. The linear fitting is satisfactory over a low pressure range where adsorption takes place by micropore filling in small micropores (pore size < 1.4 nm). However, at higher pressures the experimental isotherms tend to deviate from linearity showing an upward turn toward the -axis and this characteristic corresponds to adsorption in larger supermicropores or some mesopores [41]. The pressure at which the isotherm points start to bend upward appears to increase with increasing carbonization temperatures, having the values of about 100 kPa for the 250°C char and 400 kPa for the 550°C char. The increase of this threshold pressure indicates that the char of a higher carbonization temperature contains a larger volume of small micropores, as compared to that of the lower carbonization temperature char, thus enabling the adsorption to proceed to a higher pressure. Evidence to support this view is seen from Figure 7 for the pore size distribution of chars, as will be discussed later on.

Table 3 and Figure 6 show the porous properties of coconut shell chars estimated by the DR equation and the GCMC simulation. The results on the micropore volume and micropore area of coconut shell chars derived from DR equation and GCMC simulation agree reasonably well, differing on average by about 10% except for the lowest carbonization temperature char that gives 20% and 70% differences for the micropore volume and micropore area, respectively. Porous properties of chars increase linearly with the increase of carbonization temperature from 250 to 550°C. This suggests that over this temperature range the creation of internal pores basically involves the primary decomposition of biomass components of coconut shell, probably without secondary reactions such as thermal cracking, repolymerization, and recondensation of devolatilization products that could occur at higher temperatures [18].

The DR equation gives the largest range of micropore volume and micropore area varying from 0.064 to 0.217 cm3/g and 167 to 566 m2/g, respectively, when the carbonization temperature is increased from 250 to 550°C. It is also noted that over this temperature range the average pore diameter of chars as determined by GCMC simulation method increases from 1.25 to 1.47 nm, indicative of pore enlargement by the effect of increasing carbonization temperature.

Results on pore size distribution of coconut shell chars estimated by the GCMC simulation are presented in Figure 7. It is seen that the pore size distributions of all chars are in the range from 0.8 to 2.2 nm. There is a tendency for the distribution to become wider as the carbonization temperature is increased but the difference is not pronounced, being in the ranges 0.8 to 1.6 nm, 0.9 to 2.0 nm, 0.9 to 2.1 nm, and 1.0 to 2.2 nm for chars carbonized at 250, 350, 450, and 550°C, respectively. All pores of the 250°C char are micropores (pore size < 2 nm), while chars of higher carbonization temperatures contain mostly micropores with around 5% being in the mesopore size range (2 to 50 nm). It is noted that the volume of micropores smaller than 1.4 nm, which is responsible for the major adsorption by micropore filling, tends to increase with increasing carbonization temperature, in support of the previous discussion on the deviation from linearity of the DR plots.

3.8. Adsorption Kinetics of CO2 and Diffusivity Equations

The kinetic analysis for CO2 adsorption by coconut shell char is obtained by writing a mass balance equation for the transport of an adsorbate molecule inside a single spherical adsorbent particle [42], as shown in where is the adsorbate concentration, is the effective pore diffusivity, and is the radial distance. Upon solving (6) using appropriate initial and boundary conditions, the final solution gives the concentration of adsorbate as a function of time and radial position inside the particle, .

The solution equation for the fractional uptake of the adsorbate () is then given by the equation where is average concentration of adsorbate in the particle at time given as follows: where is particle radius = 0.72 mm, is initial concentration of adsorbate, taken as zero for a fresh adsorbent, and is adsorbate concentration at equilibrium, obtained from an isotherm equation.

To facilitate the calculation, only the first two terms of the series solution in (7) are used. This gives

Typical kinetic data obtained during the transient measurement by the IGA are shown plotted in Figure 8. These data were fitted with (9) using a nonlinear regression algorithm to obtain the best fit of the effective pore diffusivity () in the equation as a function of applied pressure.

The correlation between and the equilibrium adsorption loading of CO2 () can now be obtained by converting the pressure to via the Freundlich isotherm equation and the results are delineated in Figure 9. As seen, char prepared at a higher carbonization temperature has a higher value of diffusivity for all adsorption loading, due principally to its larger pore size and pore volume. Further, the pore diffusivity increases with increasing CO2 loading and passes through a maximum () at an optimum value of (). The values of and the corresponding are 5.40 × 10−11 m2/s and 0.70 mmol/g, 1.68 × 10−10 m2/s and 1.54 mmol/g, 2.44 × 10−10 m2/s and 1.87 mmol/g, and 3.19 × 10−10 m2/s and 2.10 mmol/g for chars prepared at the carbonization temperatures of 250, 350, 450, and 550°C, respectively. The drop of the pore diffusivity is obviously caused by the blocking of pore space when the adsorption loading is sufficiently high and this pore blocking effect persists to a higher optimum loading as the carbonization temperature is increased. It should be noted that the value of (3.19 × 10−10 m2/s) for CO2 adsorption in the 550°C char is about two orders of magnitude smaller than the CO2 diffusivity at 308 K (2.61 × 10−8 m2/s) in BPL activated carbon [42], indicative of lower pore diffusional resistance of activated carbon due to its more developed and well-connected pore structure.

An attempt is now made to estimate the pore size at which the diffusivity starts to decline () for all chars. By assuming a cylindrical pore of length and an average diameter , the pore volume () is given byThe adsorbed volume corresponding to () can be calculated from the product of a cross-sectional area of the adsorbed layer and the pore length, so thatEliminating from (10) and (11) and rearranging the terms givewhere is and is the density of liquid carbon dioxide.

By employing (12), the values of are estimated to be 0.86, 0.95, 0.99, and 1.07 nm, for the carbonization temperatures of 250, 350, 450, and 550°C, respectively. If we consider that the variation of is not so great, the average value of about 1 nm indicates that the pore resistance comes into effect when the adsorption continues until the pore size reduces to approximately three times the collision diameter (size) of the CO2 molecule.

Finally, it is worthwhile to have a mathematical function to describe the relationship between the effective pore diffusivity and equilibrium loading of CO2. This is first accomplished by writing down the following three related expressions:(1)The adsorption rate equation based on the linear driving force (LDF) model [43],where is the adsorbed-phase concentration in equilibrium with the adsorbate concentration in the bulk phase.(2)The adsorption kinetic equation of (7), considering only the first term of the series, that is,(3)The equation for calculating the average adsorbed phase concentration which is derived by integrating (8), assuming a parabolic concentration profile of the adsorbate in the adsorbent particle (). Thus,where and are equation constants.

By combining (13) to (15) and considering the condition at equilibrium, , and , the diffusivity equation as a function of CO2 loading is finally obtained,

However, this functional equation predicts a continued increase of with . To account for the fall of at a certain CO2 loading, (16) is modified by multiplying with an empirical equation,where and are constants.

As a result, the final equation for describing the effective pore diffusivity as a function of equilibrium CO2 loading for the adsorption of CO2 at 273 K by coconut shell chars is in m2/s and in mmol/g.

The four constants in (18), being dependent on the carbonization temperature, are then determined by fitting the equation against the diffusivity data employing a nonlinear regression technique and the results are listed in Table 4. The predictive capability of the proposed equation is shown in Figure 9. Overall, the agreement between the computed and experimental values is satisfactory but the equation fails to predict the experimental results for CO2 loading less than about 0.32 mmol/g of char from which the adsorption takes place in very small pores.

4. Conclusions

Porous properties of coconut shell chars prepared at carbonization temperatures in the range from 250 to 550°C were determined from the isotherm data of CO2 adsorption at 273 K and pressures up to 500 kPa, employing the DR equation and GCMC simulation. It was found that all chars contain largely micropores with a small percent being in the mesopore size range. The range of pore size distribution tends to increase with increasing carbonization temperature and varies from 0.8 to 2.2 nm. The micropore volume, micropore area, and average pore size appear to increase with the increase of carbonization temperature, being in the range of 0.064 to 0.217 cm3/g, 167 to 566 m2/g, and 1.25 to 1.47 nm, respectively, for carbonization temperatures varying from 250 to 550°C. The adsorption kinetics of CO2 in the pores of a coconut shell char is well described by the diffusional transport model and the characteristic effective pore diffusivity increases with increasing CO2 loading, showing a maximum at an optimum adsorption loading. The maximum pore diffusivity shifts to a higher loading as the carbonization temperature is increased, giving the values of 5.4 × 10−11 m2/s at 0.70 mmol/g loading for 250°C char to 3.19 × 10−10 m2/s at 2.10 mmol/g loading for 550°C char. A semiempirical correlation was proposed to predict the variation of the pore diffusivity with the equilibrium loading of CO2 adsorption.

Competing Interests

The authors declare that there is no conflict of interests concerning the publication of this paper.

Acknowledgments

The authors gratefully acknowledge the financial support for this work from the Research and Development Fund, Suranaree University of Technology.