Carbon Dioxide Adsorption in Nanopores of Coconut Shell Chars for Pore Characterization and the Analysis of Adsorption Kinetics.
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 [1-4]. 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) . This is partly explained by the low diffusion rate of the [N.sub.2] 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  because the higher adsorption temperature of C[O.sub.2] 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 . 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 C[O.sub.2] uptake data derived from gravimetric adsorption equipment for the simultaneous determination of both pore texture of coconut shell chars and the kinetics of C[O.sub.2] adsorption. Chars prepared from the carbonization of coconut shell in an inert atmosphere of nitrogen at temperatures over the range from 250 to 550[degrees]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[degrees]C overnight and the dry sample was crushed and sieved to obtain a screen size fraction of 12 x 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 . 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[degrees]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 [N.sub.2] at the rate of 100 [cm.sup.3]/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[degrees]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 [N.sub.2]. The carbonization temperature studied was varied in the range from 250 to 550[degrees]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 C[O.sub.2] 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[degrees]C for 10 h at a vacuum pressure of 75 mmHg. Next, the adsorption of C[O.sub.2] was carried out at 273 K by dosing C[O.sub.2] at a specified pressure and the increasing amount of C[O.sub.2] 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 C[O.sub.2] 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%) , oil palm fiber (18.23%) , palm kernel cake (16.74%) , palm shell (19.20%) , and eucalyptus wood (18.30%) .
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[degrees] 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[degrees]C and the second peak from 315 to 385[degrees]C. The maximum decomposition rates occur at 265[degrees] C and 350[degrees] 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[degrees]C and continues over a wide temperature range up to 900[degrees]C . Hemicellulose and cellulose decompose over the temperature range from 200 to 380[degrees]C and 250 to 400[degrees]C, respectively . 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[degrees]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[degrees]C which results from the release of noncarbon elements such as oxygen, nitrogen, hydrogen, and sulfur as volatile gaseous products . 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[degrees]C possess the true densities of 1.19,1.28,1.39, and 1.51 g/[cm.sup.3], 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[degrees] C.
3.6. Adsorption Isotherms of C[O.sub.2] by Coconut Shell Chars. Figure 4 shows isotherms of C[O.sub.2] adsorption at 273 K by coconut shell chars prepared at various carbonization temperatures from 250 to 550[degrees]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 . 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 (P/[P.sup.o] = 0.003, where [P.sup.o] = saturation vapor pressure of C[O.sub.2] 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 C[O.sub.2] at low pressure range possibly involves the adsorption in small narrow micropores of a few molecular dimensions of theadsorptivemolecules, known as primarymicropore filling . Our simulation study of C[O.sub.2] adsorption in porous carbons using a slit-pore model  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 C[O.sub.2] 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 . It is also interesting to note that the Type II isotherm found for C[O.sub.2] adsorption in coconut shell chars was also observed by a number of investigators for C[O.sub.2] adsorption in activated carbons having much larger internal surface area [27-30]. This indicates that the basic interaction forces between C[O.sub.2] 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 C[O.sub.2] molecules .
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 . 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:
[mathematical expression not reproducible], (1)
where [q.sub.F] is amount of C[O.sub.2] adsorbed, mmol/g, P is pressure, kPa, [K.sub.F] is affinity constant, and [n.sub.F] 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 form
[mathematical expression not reproducible], (2)
where [q.sub.S] and [q.sub.m] are amount adsorbed and monolayer adsorbed capacity in mmol/g, respectively, P is pressure, kPa, b is affinity constant, and ns 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 C[O.sub.2] adsorption at 273 K by coconut shell chars prepared at various carbonization temperatures. As seen from Figure 4 and the values of regression coefficient ([R.sup.2]) in Table 2, the isotherm data are very well described by both adsorption models. It is interesting to observe that the surface heterogeneity constants, [n.sub.F] and [n.sub.S], 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 . The adsorbed capacity constants, [K.sub.F] and [q.sub.m], also increase with the increase of carbonization temperature, indicating the consequent increase of available surface area for C[O.sub.2] adsorption. However, the affinity constant, b, 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 t-plot and [[alpha].sub.s]-plot . The pore size distribution can be determined by the MP method  and the Horvath-Kawazoe (HK) method  or by more modern techniques of the density functional theory (DFT)  or the grand canonical Monte Carlo simulation (GCMC) . 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 C[O.sub.2] adsorption isotherms at 273 K.
The DR equation was developed based on the potential theory of Polanyi  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  has proposed the specific empirical functional form of the characteristic curve for microporous adsorbents such as activated carbon. The final form of DR equation reads
V/[V.sub.o] = exp [- [[RT ln([P.sup.o]/P)].sup.2]/[beta][E.sub.o.sup.2]] (3)
or in linear form,
log V = log [V.sub.o] - 2.303 [(RT/[beta][E.sub.o]).sup.2] [log.sup.2] ([P.sup.o]/P), (4)
where V is the amount adsorbed at pressure P, mmol/g adsorbent, [V.sub.o] is the amount adsorbed at the saturation pressure [P.sup.o], mmol/g, R is gas constant = 8.314 x [10.sup.-3] kJ/mol-K, T is adsorption temperature, K, E0 is characteristic energy, kJ/mol, [beta] 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 [E.sub.o] and [V.sub.o] are obtained from the slope and intercept, respectively. The micropore volume ([V.sub.mic]) of the adsorbent (coconut shell chars) is obtained by converting [V.sub.o] into the volume of adsorbate in liquid state (liquid density of C[O.sub.2] is 0.0207 mol/[cm.sup.3]) and the micropore area ([A.sub.mic]) in [m.sup.2]/g can be calculated from [V.sub.o] by the equation
[A.sub.mic] = [V.sub.o][a.sub.m][N.sub.a], (5)
where [V.sub.o] is the micropore adsorbed capacity (mol/g), [a.sub.m] is the cross-sectional area of C[O.sub.2] (2.1 x [10.sup.-9] [m.sup.2]/molecule at 273 K), and [N.sub.a] is the Avogadro number (6.02 x [10.sup.23] molecules/mol).
The calculation method used for the characterization of porous chars by the GCMC simulation is that proposed by Herrera et al. . 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 3nm. 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 ([D.sub.av]), the cumulated surface area, and the total pore volume are then calculated.
Figure 5 shows the DR plots for C[O.sub.2] 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 y-axis and this characteristic corresponds to adsorption in larger supermicropores or some mesopores . 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[degrees]C char and 400 kPa for the 550[degrees]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[degrees]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 .
The DR equation gives the largest range of micropore volume and micropore area varying from 0.064 to 0.217 [cm.sup.3]/g and 167 to 566[m.sup.2]/g, respectively, when the carbonization temperature is increased from 250 to 550[degrees]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[degrees]C, respectively. All pores of the 250[degrees]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 C[O.sub.2] and Diffusivity Equations. The kinetic analysis for C[O.sub.2] 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 , as shown in
[partial derivative]q/[partial derivative]t = [D.sub.e] [[[partial derivative].sup.2]/[partial derivative][r.sup.2] + 2[[partial derivative].sup.q]/r[partial derivative]r], (6)
where q(t) is the adsorbate concentration, [D.sub.e] is the effective pore diffusivity, and r 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, q = q(r, t).
The solution equation for the fractional uptake of the adsorbate (F) is then given by the equation
[mathematical expression not reproducible], (7)
where [[bar.q].sub.(t)] is average concentration of adsorbate in the particle at time t given as follows:
[mathematical expression not reproducible], (8)
where [R.sub.P] is particle radius = 0.72 mm, [q.sub.(0)] is initial concentration of adsorbate, taken as zero for a fresh adsorbent, and [q.sub.e] 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
[mathematical expression not reproducible]. (9)
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 ([D.sub.e]) in the equation as a function of applied pressure.
The correlation between [D.sub.e] and the equilibrium adsorption loading of C[O.sub.2] (qe) can now be obtained by converting the pressure to qe 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 C[O.sub.2] loading and passes through a maximum ([D.sub.e,max]) at an optimum value of [q.sub.e] ([q.sub.op]). The values of De max and the corresponding [q.sub.op] are 5.40 x [10.sup.-11] [m.sup.2]/s and 0.70 mmol/g, 1.68 x [10.sup.-10] [m.sup.2]/s and 1.54 mmol/g, 2.44 x [10.sup.-10] [m.sup.2]/s and 1.87 mmol/g, and 3.19 x [10.sup.-10] [m.sup.2]/s and 2.10 mmol/g for chars prepared at the carbonization temperatures of 250, 350, 450, and 550[degrees]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 [D.sub.e,max] (3.19 x [10.sup.-10] [m.sup.2]/s) for C[O.sub.2] adsorption in the 550[degrees]C char is about two orders of magnitude smaller than the C[O.sub.2] diffusivity at 308 K (2.61 x [10.sup.-8] [m.sup.2]/s) in BPL activated carbon , 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 ([D.sub.op]) for all chars. By assuming a cylindrical pore of length L and an average diameter [D.sub.av], the pore volume ([V.sub.T]) is given by
[V.sub.T] = ([pi]/4) [D.sub.av.sup.2] L. (10)
The adsorbed volume corresponding to [q.sub.op] ([V.sub.op]) can be calculated from the product of a cross-sectional area of the adsorbed layer and the pore length, so that
[V.sub.op] = ([pi]/4) ([D.sub.av.sup.2] - [D.sub.op.sup.2]) L. (11)
Eliminating L from (10) and (11) and rearranging the terms give
[D.sub.op] = [D.sub.av] [(1 - [V.sub.op]/[V.sub.T]).sup.1/2], (12)
where [V.sub.op] is [q.sub.op]/[[rho].sub.L] and [[rho].sub.L] is the density of liquid carbon dioxide.
By employing (12), the values of [D.sub.op] are estimated to be 0.86, 0.95, 0.99, and 1.07 nm, for the carbonization temperatures of 250, 350, 450, and 550[degrees]C, respectively. If we consider that the variation of [D.sub.op] 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 C[O.sub.2] molecule.
Finally, it is worthwhile to have a mathematical function to describe the relationship between the effective pore diffusivity and equilibrium loading of C[O.sub.2]. 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 ,
[partial derivative][bar.q]/[partial derivative]t = 15[D.sub.e] ([q.sup.*] - q)/[R.sub.P.sup.2], (13)
where [q.sup.*] 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 q which is derived by integrating (8), assuming a parabolic concentration profile of the adsorbate in the adsorbent particle (q = [C.sub.1] + [C.sub.2][r.sup.2]). Thus,
[bar.q] = [C.sub.1] + 0.6[C.sub.2][R.sub.P.sup.2], (15)
where [C.sub.1] and [C.sub.2] are equation constants.
By combining (13) to (15) and considering the condition at equilibrium, t = [t.sub.e], and [q.sup.*] = [q.sub.e], the diffusivity equation as a function of C[O.sub.2] loading is finally obtained,
[D.sub.e] = (5.25 x [10.sup.-8]/[t.sup.e]) x ln [[q.sub.e] - [C.sub.1] - 3.11 x [10.sup.-3][C.sub.2]/0.4[q.sub.e]], [m.sup.2]/s. (16)
However, this functional equation predicts a continued increase of [D.sub.e] with [q.sub.e] To account for the fall of [D.sub.e] at a certain C[O.sub.2] loading, (16) is modified by multiplying with an empirical equation,
[mathematical expression not reproducible], (17)
where [C.sub.3] and [C.sub.4] are constants.
As a result, the final equation for describing the effective pore diffusivity as a function of equilibrium C[O.sub.2] loading for the adsorption of C[O.sub.2] at 273 K by coconut shell chars is
[mathematical expression not reproducible], (18)
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 C[O.sub.2] loading less than about 0.32 mmol/g of char from which the adsorption takes place in very small pores.
Porous properties of coconut shell chars prepared at carbonization temperatures in the range from 250 to 550[degrees]C were determined from the isotherm data of C[O.sub.2] 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 [cm.sup.3]/g, 167 to 566 [m.sup.2]/g, and 1.25 to 1.47 nm, respectively, for carbonization temperatures varying from 250 to 550[degrees]C. The adsorption kinetics of C[O.sub.2] 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 C[O.sub.2] 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 x [10.sup.-11] [m.sup.2]/s at 0.70 mmol/g loading for 250[degrees]C char to 3.19 x [10.sup.-10] [m.sup.2]/s at 2.10 mmol/g loading for 550[degrees]C char. A semiempirical correlation was proposed to predict the variation of the pore diffusivity with the equilibrium loading of C[O.sub.2] adsorption.
The authors declare that there is no conflict of interests concerning the publication of this paper.
The authors gratefully acknowledge the financial support for this work from the Research and Development Fund, Suranaree University of Technology.
 T. Kan, V. Strezov, and T. J. Evans, "Lignocellulosic biomass pyrolysis: a review of product properties and effects of pyrolysis parameters," Renewable and Sustainable Energy Reviews, vol. 57, pp. 1126-1140, 2016.
 A. Sharma, V Pareek, and D. Zhang, "Biomass pyrolysis--a review of modelling, process parameters and catalytic studies," Renewable and Sustainable Energy Reviews, vol. 50, pp. 1081-1096, 2015.
 T. Aysu and M. M. Kucuk, "Biomass pyrolysis in a fixed-bed reactor: effects of pyrolysis parameters on product yields and characterization of products," Energy, vol. 64, pp. 1002-1025, 2014.
 F.-X. Collard and J. Blin, "A review on pyrolysis of biomass constituents: mechanisms and composition of the products obtained from the conversion of cellulose, hemicelluloses and lignin," Renewable and Sustainable Energy Reviews, vol. 38, pp. 594-608, 2014.
 Y. Shen, "Chars as carbonaceous adsorbents/catalysts for tar elimination during biomass pyrolysis or gasification," Renewable and Sustainable Energy Reviews, vol. 43, pp. 281-295, 2015.
 Y. W. Lee, J. Park, C. K. Ryu et al., "Comparison of biochar properties from biomass residues produced by slow pyrolysis at 500[degrees]C," Bioresource Technology, vol. 148, pp. 196-201, 2013.
 A. Alabadi, S. Razzaque, Y. Yang, S. Chen, and B. Tan, "Highly porous activated carbon materials from carbonized biomass with high C[O.sub.2] capturing capacity," Chemical Engineering Journal, vol. 281, pp. 606-612, 2015.
 N. M. Nor, L. C. Lau, K. T. Lee, and A. R. Mohamed, "Synthesis of activated carbon from lignocellulosic biomass and its applications in air pollution control--a review," Journal of Environmental Chemical Engineering, vol. 1, no. 4, pp. 658-666, 2013.
 J. W. Patricks, Porosity in Carbons: Characterization and Applications, Edward Arnold Publisher, 1995.
 S. J. Gregg and K. S. W. Sing, Adsorption, Surface Area and Porosity, Academic Press, 2nd edition, 1982.
 H. Marsh and F. Rodriguez-Reinoso, Activated Carbon, Elsevier, 2006.
 D. Lozano-Castello, D. Cazorla-Amoros, and A. Linares-Solano, "Usefulness of C[O.sub.2] adsorption at 273 K for the characterization of porous carbons," Carbon, vol. 42, no. 7, pp. 1233-1242, 2004.
 D. Ozcimen and A. Ersoy-Mericboyu, "Carbonization kinetics of various biomass sources," Journal of Chemical Engineering Research Updates, vol. 1, no. 1, pp. 20-28, 2014.
 C. Tangsathitkulchai, S. Junpirom, and J. Katesa, "Comparison of kinetic models for C[O.sub.2] gasification of coconut-shell chars: carbonization temperature effects on char reactivity and porous properties of produced activated carbons," Engineering Journal, vol. 17, no. 1, pp. 13-28, 2013.
 D. Lopez-Gonzolez, M. Fernandez-Lopez, J. L. Valverde, and L. Sanchez-Silva, "Gasification of lignocellulosic biomass char obtained from pyrolysis: kinetic and evolved gas analyses," Energy, vol. 71, pp. 456-467, 2014.
 S. Junpirom, C. Tangsathitkulchai, and M. Tangsathitkulchai, "Thermogravimetric analysis of longan seed biomass with a two-parallel reactions model," Korean Journal of Chemical Engineering, vol. 27, no. 3, pp. 791-801, 2010.
 A. C. Lua and J. Guo, "Preparation and characterization of chars from oilpalm waste," Carbon, vol. 36, no. 11, pp. 1663-1670,1998.
 P. Weerachanchai, C. Tangsathitkulchai, and M. Tangsathitkulchai, "Characterization of products from slow pyrolysis of palm kernel cake and cassava pulp residue," Korean Journal of Chemical Engineering, vol. 28, no. 12, pp. 2262-2274, 2011.
 P. Luangkiattikhun, C. Tangsathitkulchai, and M. Tangsathitkulchai, "Non-isothermal thermogravimetric analysis of oilpalm solid wastes," Bioresource Technology, vol. 99, no. 5, pp. 986-997, 2008.
 Y. Ngernyen, C. Tangsathitkulchai, and M. Tangsathitkulchai, "Porous properties of activated carbon produced from Eucalyptus and Wattle wood by carbon dioxide activation," Korean Journal of Chemical Engineering, vol. 23, no. 6, pp. 1046-1054, 2006.
 H. Yang, R. Yan, H. Chen, D. H. Lee, and C. Zheng, "Characteristics of hemicellulose, cellulose and lignin pyrolysis," Fuel, vol. 86, no. 12-13, pp. 1781-1788, 2007.
 D. Vamvuka, E. Kakaras, E. Kastanaki, and P. Grammelis, "Pyrolysis characteristics and kinetics of biomass residuals mixtures with lignite," Fuel, vol. 82, no. 15-17, pp. 1949-1960, 2003.
 R. C. Bansal and M. Goyal, Activated Carbon Adsorption, Taylor & Francis, 2005.
 S. Lowell, J. E. Shields, M. A. Thomas, and M. Thomas, Characterization of Porous Solids and Powders: Surface Area, Pore Size and Density, Kluwer Academic, Dordrecht, the Netherlands, 2004.
 F. Rouquerol, J. Rouquerol, and K. Sing, Adsorption by Powders and Porous Solids, Academic Press, 1999.
 A. Wongkoblap, W. Intomya, W. Somrup, S. Charoensuk, S. Junpirom, and C. Tangsathitkulchai, "Pore size distribution of carbon with different probe molecules," Engineering Journal, vol. 14, no. 3, pp. 45-56, 2010.
 S. Deng, H. Wei, T. Chen, B. Wang, J. Huang, and G. Yu, "Superior C[O.sub.2] adsorption on pine nut shell-derived activated carbons and the effective micropores at different temperatures," Chemical Engineering Journal, vol. 253, pp. 46-54, 2014.
 J. Srenscek-Nazzal, U. Narkiewicz, A. W. Morawski, R. J. Wrobel, and B. Michalkiewicz, "Comparison of optimized isotherm models and error functions for carbon dioxide adsorption on activated carbon,," Journal of Chemical & Engineering Data, vol. 60, no. 11, pp. 3148-3158, 2015.
 R.-L. Tseng, F.-C. Wu, and R.-S. Juang, "Adsorption of C[O.sub.2] at atmospheric pressure on activated carbons prepared from melamine-modified phenol-formaldehyde resins," Separation and Purification Technology, vol. 140, pp. 53-60, 2015.
 B. B. Saha, S. Jribi, S. Koyama, and I. I. El-Sharkawy, "Carbon dioxide adsorption isotherms on activated carbons," Journal of Chemical & Engineering Data, vol. 56, no. 5, pp. 1974-1981, 2011.
 Y. Liu and J. Wilcox, "Molecular simulation studies of C[O.sub.2] adsorption by carbon model compounds for carbon capture and sequestration applications," Environmental Science & Technology, vol. 47, no. 1, pp. 95-101, 2013.
 D. D. Do, Adsorption Analysis: Equilibria and Kinetics, Imperial College Press, 1998.
 W. Li, K. Yang, J. Peng, L. Zhang, S. Guo, and H. Xia, "Effects of carbonization temperatures on characteristics of porosity in coconut shell chars and activated carbons derived from carbonized coconut shell chars," Industrial Crops and Products, vol. 28, no. 2, pp. 190-198, 2008.
 R. S. Mikhail, S. Brunauer, and E. E. Bodor, "Investigations of a complete pore structure analysis: I. Analysis of micropores," Journal of Colloid and Interface Science, vol. 26, no. 1, pp. 45-53, 1968.
 G. Horvath and K. Kawazoe, "Method for the calculation of effective pore size distribution in molecular sieve carbon," Journal of Chemical Engineering of Japan, vol. 16, no. 6, pp. 470-475, 1983.
 P. I. Ravikovitch, G. L. Haller, and A. V. Neimark, "Density functional theory model for calculating pore size distributions: pore structure of nanoporous catalysts," Advances in Colloid and Interface Science, vol. 76-77, pp. 203-226, 1998.
 S. Samios, A. K. Stubos, N. K. Kanellopoulos, R. F. Cracknell, G. K. Papadopoulos, and D. Nicholson, "Determination of micropore size distribution from grand canonical Monte Carlo simulations and experimental C[O.sub.2] isotherm data," Langmuir, vol. 13, no. 10, pp. 2795-2802, 1997.
 M. Polanyi, "Theories of the adsorption of gases," Transactions of the Faraday Society, vol. 28, pp. 316-333, 1932.
 M. M. Dubinin, "Adsorption in micropores," Journal of Colloid and Interface Science, vol. 23, no. 4, pp. 487-499, 1967.
 L. F. Herrera, C. Fan, V. Nguyen, D. D. Do, T. Horikawa, and D. Nicholson, "A self-consistent method to determine accessible volume, area and pore size distribution (APSD) of BPL, Norit and AX-21 activated carbon," Carbon, vol. 50, no. 2, pp. 500-509, 2012.
 N. Kanellopoulos, Nanoporous Materials: Advanced Techniques for Characterization, Modelingand Processing, CRC Press, 2011.
 J. A. Delgado, V. I. Agueda, M. A. Uguina, J. L. Sotelo, P. Brea, and C. A. Grande, "Adsorption and diffusion of [H.sub.2], CO, C[H.sub.4], and C[O.sub.2] in BPL activated carbon and 13X zeolite: evaluation of performance in pressure swing adsorption hydrogen purification by simulation," Industrial & Engineering Chemistry Research, vol. 53, no. 40, pp. 15414-15426, 2014.
 E. Glueckauf, "Theory of chromatography. Part 10.--formula for diffusion into spheres and their application to chromatography," Transactions of the Faraday Society, vol. 51, pp. 1540-1551, 1955.
Chaiyot Tangsathitkulchai, Supunnee Junpirom, and Juejun Katesa
School of Chemical Engineering, Institute of Engineering, Suranaree University of Technology, Muang District, Nakhon Ratchasima 30000, Thailand
Correspondence should be addressed to Chaiyot Tangsathitkulchai; email@example.com
Received 20 August 2016; Accepted 27 November 2016
Academic Editor: Ajayan Vinu
Caption: FIGURE 1: TGA and DTG curves of coconut shell used in this study.
Caption: FIGURE 2: SEM images of coconut shell and coconut shell chars prepared at various carbonization temperatures.
Caption: FIGURE 3: Effect of carbonization temperature on the proximate analysis of coconut shell chars.
Caption: FIGURE 4: Isotherms of C[O.sub.2] adsorption at 273 K by coconut shell chars prepared at various carbonization temperature. The solid lines represent isotherms predicted by the Freundlich and Sips isotherm equations.
Caption: FIGURE 5: DR plots for C[O.sub.2] adsorption at 273 K by coconut shell chars prepared at different carbonization temperatures ([P.sup.o] = 3500kPa for C[O.sub.2] at 273 K).
Caption: FIGURE 6: Effect of carbonization temperature on the porous properties of coconut shell chars, calculated by DR equation and GCMC simulation.
Caption: FIGURE 7: Pore size distribution of coconut shell chars prepared at various carbonization temperatures calculated by GCMC simulation technique.
Caption: FIGURE 8: Typical uptake of C[O.sub.2] adsorption at 273 K by coconut shell char as a function of adsorption time and applied pressure measured by the IGA instrument.
Caption: FIGURE 9: Effect of C[O.sub.2] loading on effective pore diffusivity ([D.sub.e]) of C[O.sub.2] at 273 K in coconut shell chars prepared at various carbonization temperatures.
Table 1: Proximate and ultimate analyses of coconut shell used in this study. Proximate analysis (dry basis), wt% Volatile matter 80.80 Fixed carbon 19.07 Ash content 0.13 Ultimate analysis, wt% C 49.75 H 5.60 O (by diff.) 44.28 N 0.35 S 0.02 TABLE 2: Fitted parameters of Freundlich and Sips isotherm equations for C[O.sub.2] adsorption at 273 K by coconut-shell chars prepared at various carbonization temperatures. Freundlich equation Carbonization temperature [K.sub.F] [n.sub.F] [R.sup.2] [q.sub.m] ([degrees]C) 250 0.0327 1.43 0.9991 30.4 350 0.1130 1.83 0.9988 32.3 450 0.2192 2.17 0.9987 47.4 550 0.3577 2.43 0.9968 97.7 Sips equation Carbonization temperature b [n.sub.s] [R.sup.2] ([degrees]C) 250 7.4 x [10.sup.-5] 1.37 0.9985 350 4.8 x [10.sup.-5] 1.73 0.9981 450 1.3 x [10.sup.-5] 2.08 0.9987 550 2.0 x [10.sup.-6] 2.36 0.9980 TABLE 3: Calculated porous properties of coconut-shell chars prepared at various carbonization temperatures based on DR equation and GCMC simulation method. DR equation Carbonization [V.sub.mic] [A.sub.mic] [E.sub.o] temperature ([cm.sup.3]/g) ([m.sup.2]/g) (kJ/mol) ([degrees]C) 250 0.064 167 28.7 350 0.122 320 29.0 450 0.169 443 29.7 550 0.217 566 29.9 GCMC simulation Carbonization [V.sub.mic] [A.sub.mic] [E.sub.o] temperature ([cm.sup.3]/g) ([m.sup.2]/g) (kJ/mol) ([degrees]C) 250 0.079 285 1.25 350 0.125 356 1.45 450 0.153 418 1.46 550 0.189 503 1.47 TABLE 4: Values of constants in (18) for predicting the effective pore diffusivity of C[O.sub.2] ([D.sub.e]) at 273 K in various coconut shell chars. Carbonization [t.sub.e] [C.sub.1] [C.sub.2] [C.sub.3] temperature ([degrees]C) 250 0.38 -0.06 80.5 4 x [10.sup.-6] 350 228 -0.37 182 2.92 450 168 0.25 -15.5 3.61 550 137 -0.28 132 4.57 Carbonization [C.sub.4] [R.sup.2] temperature ([degrees]C) 250 0.60 0.6805 350 4.21 0.9733 450 5.44 0.9725 550 4.94 0.8923
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Tangsathitkulchai, Chaiyot; Junpirom, Supunnee; Katesa, Juejun|
|Publication:||Journal of Nanomaterials|
|Date:||Jan 1, 2017|
|Previous Article:||Surface Plasmon Enhanced Light Trapping in Metal/Silicon Nanobowl Arrays for Thin Film Photovoltaics.|
|Next Article:||Cell Growth Inhibition Effect of DsiRNA Vectorised by Pectin-Coated Chitosan-Graphene Oxide Nanocomposites as Potential Therapy for Colon Cancer.|