Analysis of Siloxane Adsorption Characteristics Using Response Surface Methodology
Analysis of Siloxane Adsorption Characteristics Using Response Surface Methodology
Environmental Engineering Research. 2012. Jun, 17(2): 117-122
Copyright ©2012, Korean Society of Environmental Engineering
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : March 03, 2012
  • Accepted : May 05, 2012
  • Published : June 30, 2012
Export by style
Cited by
About the Authors
Jin-Kyu Park
Department of Environmental Engineering, Anyang University, Anyang 430-714, Korea
Gyeung-Mi Lee
Department of Environmental Engineering, Anyang University, Anyang 430-714, Korea
Chae-Young Lee
Department of Civil Engineering, The University of Suwon, Suwon 445-743, Korea
Kwang-Beom Hur
Korea Electric Power Corporation, Daejeon 305-706, Korea
Nam-Hoon Lee
Department of Environmental Engineering, Anyang University, Anyang 430-714, Korea
A central composite design and response surface methodology were applied to investigate the optimum conditions for maximum adsorption capacity in activated alumina as an adsorbent. The optimized conditions were determined for adsorption capacity using variables of flow rate and temperature. It was found that flow rate and temperature greatly influenced the adsorption capacity, as determined by analysis of variance analysis of these variables. Statistical checks indicated that second order polynomial equations were adequate for representing the experimental values. The optimum conditions for adsorption capacity were 0℃ and 2,718 mL/min, with the estimated maximum adsorption capacity of 17.82%. The experimental adsorption capacity was 17.75% under these optimum conditions, which was in agreement with the predicted value of 17.82%.
1. Introduction
Worldly consumption of fossil fuel has multiplied greenhouse gas emission, which has stimulated researchers’ interest in development of renewable energy instead of traditional energy sources. Biogas production from municipal solid wastes or organic wastes is one of the renewable energy sources and is getting more attention because of the attractiveness of simultaneous waste treatment and electrical generation [1 , 2] .
However, to use biogas as an energy source, several trace compounds such as hydrogen sulfide, ammonia and siloxanes should be removed [3] . Among these compounds, siloxanes have negative influence on the capability and life duration of generators [4] . Siloxanes are a subgroup of silicones, containing Si-O bonds with organic radicals (mostly methyl-groups) bound to Si [5] . The structure of the siloxane is linear or cyclic. Several technical methods (e.g., adsorption, gas chilling, and absorption) are commercially available to remove the siloxanes today [6] . In adsorption, the most widely used method, various adsorbents and different conditional variables (e.g., flow rate, temperature, and concentration) have been applied to find an optimal condition for increasing adsorption capacity [6 - 8] . For the traditional searching procedure, one conditional variable at a time is changed without any changes to other variables for finding an optimal condition. This method may produce a certain optimum stage, however, interdependent influence or relations for various conditional variables are not considered in general.
Therefore, an alternative searching procedure involving a statistical tool (e.g., factorial experimental design and response surface methodology [RSM]) should be adopted to solve this complexity [9] . RSM is a collection of mathematical and statistical techniques used to model and analyze problems in which a response of interest is influenced by several variables [10] .
In this study, characteristics of the siloxane adsorption with various flow rates and temperatures were analyzed using alumina as an adsorbent of the siloxane. To improve the conventional test method which searches for the effect of each variable, RSM was applied to investigate the individual and interactive effects of flow rate and temperature on adsorption capacity.
2. Materials and Methods
- 2.1. Adsorption Test
The experimental equipment for this study is presented in Fig. 1 . Activated alumina (Duksan Co., Ansan, Korea), an adsorbent of siloxane, was circular-shaped with a white color and its density was 0.6 g/cm 3 . The setup consisted of an adsorption column (12 mm, diameter; 120 mm, length) filled with 5 g of dried granular activated alumina (65.7% porosity) heated until 300℃ for approximately 2 hr. The gas used in this column was adsorbed
PPT Slide
Lager Image
Schematic diagram of siloxane adsorption apparatus.
PPT Slide
Lager Image
Central composite design (CCD): 4 cube points, 4 axial points, and 5 replicates at the central point (0,0) [12].
in three impingers containing 6 mL of methanol. The inside of an ice box was steadily kept at a lower temperature of 4 ± 2℃. The gas flow rate into the impingers was 100 mL/min constantly, with discharging of extra gas to the outside. In the last part of the impingers, a wet gas meter was installed to keep the total gas input constant. A gas chromatography/mass spectrometry (GC-MS) system (QP5050A; SIMADZU, Kyoto, Japan) equipped with a ZB-5 column (0.25 μm × 0.25 mm × 30 m; Zebron, Newport Beach, CA, USA) was operated to analyze the siloxane stored in three impingers.
In this study, octamethylcyclotetrasiloxane (D4; Acros Organics, Fair Lawn, NJ, USA) and decamethylcyclopentasiloxane (D5; Merck-Schuchardt, Hohenbrunn, Germany) were used as
Central composite design
PPT Slide
Lager Image
Central composite design
siloxane test materials and are major components of siloxanes in biogas. To generate constant input concentration of the siloxane without any dependence of input gas rate and gas temperature, nitrogen gas is injected into two parts separately. One part of nitrogen gas injection was fed to a glass bottle containing a standard siloxane solution (D4, 100 ppm; D5, 70 ppm), and it was mixed with the other nitrogen gas injection before being introduced into the adsorption column. The siloxane solution was vaporized into the nitrogen gas. The average concentration of the injected siloxane mixture to the adsorption column was 250 μg/m 3 . Teflon pipe (5 m) was installed between mixing chamber and adsorption columns in a water bath for constant gas temperature.
- 2.2. Central Composite Design and Data Analysis
A central composite design is a second order factorial design including a central point and axial points [11] . If the number of independent variables is k, the points of factorial experiments are 2 k and axial points are 2k. Therefore, if the number of repeats at a central point is n, the total number of experiments is 2 k +2k+n. As an example, Fig. 2 shows four (2 k = 4) points as cube points, four (2k = 2) axial points, and one center point (0, 0). Therefore if k is equal to 2, the minimum number of experimental points is 9 for central composite design.
In this research, the applied flow rate and adsorption temperature as independent variables with the siloxane adsorption capacity as a dependent variable were examined with 13 experiments ( Table 1 ). A central point was replicated 5 times in order to investigate experimental error.
Factorial design levels were codified from -1 to +1. The usage of axial point ( α ) is for indicating rotatability. In central composite design, α is presented as [2 k ] 1/4 and calculated as 1.414 in this experiment. The ranges for adsorption temperature and flow rate were 0-40℃ and 700-3,400 mL/min (linear velocity, 0.1-0.5 m/sec), respectively.
The experimental design was analyzed using the RSM. The variables, X i , were coded as xi according to Eq. (1), such that X 0 corresponded to the central value [11] :
  • xi= (Xi- X*i)/ΔXi(1)
  • xi: dimensionless coded value of an independent variable
  • Xi: actual value of an independent variable for thei-th test
  • X*i: actual value of an independent variable at the central point
  • ΔXi: step change
Regression analysis for the experimental results was performed using Minitab software ver. 14 (Minitab Inc., State College, PA, USA), and Eq. (2) was applied to search the optimum conditions.
  • Y=β0+ β1X1+ β2X2+ β11X12+ β22X22+ β12X1X2(2)
  • Y: predicted response
  • X1, X2: independent variables
  • β0: constant coefficient
  • β1, β2: linear coefficients
  • β11, β22: quadratic coefficients
  • β12: interaction coefficient
3. Results and Discussion
- 3.1. Adsorption Capacity of Siloxane
In this study, as an index of the adsorption performance, the adsorption capacity (the weight of the adsorbed siloxanes/the weight of the adsorbent) was calculated [7] . The results of the siloxane adsorption capacity depending on various flow rates and temperatures in the adsorption column filled with activated alumina are shown in Fig. 3 . The 1,080 mL/min of flow rate with two separate temperatures, 6℃ and 34℃, produced adsorption capacities of 15.42% and 11.93%, respectively. The temperature of 6℃ produced better adsorption capacity. When flow rate was fixed to 2,050 mL/min with various temperatures (i.e., 0℃, 20℃, and 40℃), adsorption capacities were 17.5% (0℃) and 13.41% (40℃), respectively. The results of the adsorption capacity when flow rate was 3,000 mL/min with temperatures of 6℃ and 34℃ were 16.90% and 14.31%, respectively. From these experimental results, we clearly see that reduced temperature provides better conditions for the adsorption capacity with fixed flow rate. Ricaurte Ortega and Subrenat [13] reported that the increase in the temperature reduced the adsorption capacity. The reason for decreased adsorption capacity when temperatures were increased might be that the adsorption process between siloxane and activated alumina relies more on physical than chemical adsorption.
PPT Slide
Lager Image
Scatter plot for adsorption capacity as a function of temperature and flow rate.
Alternatively, we tried to study the pattern of fixed temperature with various flow rates. When temperature was fixed to 20℃ with diverse flow rates (700 mL/min, 2,050 mL/min, and 3,400 mL/min) the adsorption capacities were increased to 12.69%, 14.75-14.99%, and 15.7%, respectively.
Adsorbates must be transported through the boundary layer that surrounds the adsorbents when gas is flowing past them. The distance of transport is determined by the flow rate past the adsorbent. Higher flow rate results in shorter distance. Therefore, when gas flow rate was increased, the siloxane was transported more quickly through the boundary layer. After passing through this boundary layer, siloxane was transported more quickly through the adsorbent’s pores to the available adsorption sites. This rapid diffusion into pores of adsorbent might generate more reactions on the adsorbent with increased adsorption capacity [14] .
- 3.2. Model Fitting
Based on the above results, second order polynomial equations for the adsorption capacity fitted in terms of coded and uncoded factors were obtained, as shown in Eqs. (3) and (4):
  • Parameter estimate from coded data
  • Y=14.832-1.5828X1+1.1483X2+0.2151X12-0.2049X22-0.1372X1X2(3)
  • Parameter estimate from uncoded data
  • Y=13.6726-0.1341X1+0.00233X2+0.00108X12-2.2497×10-7X22-1.0167×10-6X1X2(4)
Where, X 1 and X 2 represent temperature and flow rate, respectively.
Table 2 shows experimental results and predicted values with ranges of 11.93-17.50% and 12.24-17.64%, respectively.
- 3.3. Validation of the Model
In this study, we applied two different methods (i.e., graphical and numerical methods) to validate the model. Graphical methods utilize residual analysis between the values of dependent variables and predicted values from a regression equation. Fig. 4
PPT Slide
Lager Image
Residuals vs. fitted values for adsorption capacity.
PPT Slide
Lager Image
Residuals vs. observation orders of data for adsorption capacity.
Central composite design and the results obtained
PPT Slide
Lager Image
Central composite design and the results obtained
Analysis of variance results for the adsorption capacity
PPT Slide
Lager Image
Analysis of variance results for the adsorption capacity
PPT Slide
Lager Image
Normal probability plots for adsorption capacity.
presents randomly distributed residuals and predicted values in a scatter diagram within ± 0.8. In Fig. 5 , each residual is plotted against an index of observation orders of data, which was used to check for any drift in the process. As previously shown in Figs. 4 and 5 , we are not able to see any pattern of residual distribution. From this result, we may say thatdual is independent. Fig. 6 is a graph to check the normality distribution of the residuals. In this figure, experimental points are close to the center cross line, which means the test follows the normal distribution relatively well.
The model was then checked using a numerical method employing analysis of variance (ANOVA). The results of ANOVA are shown in Table 3 . The value of p > F less than 0.05 indicates that this model is significant. For this model, the value of p > F (0.000) shows statistically strong significance. The fit of the model was also expressed by the coefficient of determination R 2 , which was found to be 0.9946, which indicates a good agreement between experimental and predicted values. The lack of fit measures the
PPT Slide
Lager Image
Two dimensional response surface plots for adsorption capacity.
PPT Slide
Lager Image
Response optimization plot.
Estimation of the second-order response surface parameters (coded unit)
PPT Slide
Lager Image
Estimation of the second-order response surface parameters (coded unit)
Confirmation experiments at optimum conditions( ): standard deviation.
PPT Slide
Lager Image
Confirmation experiments at optimum conditions ( ): standard deviation.
variation of the data around the fitted model. When the p -value for the lack of fit is more than 0.05, there is a statistically significant lack-of-fit at the 95% confidence level, which means that the model is adequate [15] . In this study, the p -value (0.108) for the lack of fit showed that the model was valid.
- 3.4. Optimization and Verification of the Model
The p -values of the linear, quadratic and interaction terms, shown in Table 4 with 95% confidence, showed that all the terms were significant. The effects of the linear terms from temperature (X 1 ) and flow rate (X 2 ) presented were significant. In terms of temperature, negative effect was presented. In the non-linear effects of temperature (X 1 ) × temperature (X 1 ) and flow rate (X 2 ) × flow rate (X 2 ), quadratic terms showed that it was significant in the model. Flow rate (X 2 ) × flow rate (X 2 ) has a negative effect. Finally, the interaction effect with 0.025 as the p -value presented significance of the adsorption capacity attributed to the mixture effects from temperature and flow rate.
In Fig. 7 , a two contour line graph shows the effects on adsorption capacity based on flow rate and temperature. The optimum status (i.e., temperature below 5℃ and flow rate over 1,500 mL/min) for adsorption capacity was presented in the two contour line graph. In the case of response surface graph, depending on the shape of the response surface, we can find the optimum conditions. However, this study had some difficulties in finding the optimum conditions under the given conditions because the contour lines are raised continuously. It seems more research is necessary to find the optimum conditions under 0℃.
The optimum conditions under the given conditions were predicted using Minitab software. Fig. 8 shows the optimum values at the top of the graph at 0℃ and 2,718 mL/min, with a maximum adsorption capacity (17.82%).
In order to verify the model, experiments were performed to compare the experimental results with corresponding predicted values. Table 5 shows that the experimental adsorption capacity of 17.75% and the predicted value of 17.82% corresponded well. The experiment results were close to the model prediction. Additionally, this study confirmed that the RSM method is appropriate for optimizing the adsorption capacity.
4. Conclusions
This research applied response surface methodology to find the optimum conditions of temperature and flow rate for siloxane adsorption capacity in activated alumina as an adsorbent. The results are shown below: 1) With faster flow rate, increased temperature minimizes siloxane adsorption capacity on activated alumina because the adsorption mechanism has a preference for physical adsorption instead of chemical adsorption relying on activate energy; 2) Under the conditions of constant temperature, flow rate and adsorption capacity have a positive relation. When gas flow rate was increased, siloxane was transported more quickly through the boundary layer. After passing through this boundary layer, siloxane was transported more quickly through the adsorbent’s pores to the available adsorption sites. This rapid diffusion into pores of adsorbent might generate more reactions on the adsorbent with increased adsorption capacity; 3) From statistical analysis on siloxane adsorption capacity, second order polynomial equations were confirmed by the ANOVA result, with a p-value lower than 0.05; and 4) In numerical optimization, optimal conditions (i.e., 0℃ and 2,718 mL/min for temperature and flow rate, respectively) were obtained with 17.82% adsorption capacity. Additionally, the three verification tests showed small error rate (0.37%) using the obtained model.
This work was supported by the Power Generation & Electricity Delivery of the Korea Institute of Energy Technology Evaluation and Planning (KETEP) grant funded by the Korea Government Ministry of Knowledge Economy (Development of Biogas Pretreatment Fuel System for 5MW-GT) (no. 2009T100200127).
Münster M , Lund H. 2010 Comparing waste-to-energy technologies by applying energy system analysis. Waste Manag. 30 1251 - 1263
Holm-Nielsen JB , Al Seadi T , Oleskowicz-Popiel P. 2009 The future of anaerobic digestion and biogas utilization. Bioresour. Technol. 100 5478 - 5484
Arnold M , Kajolinna T. 2010 Development of on-line measurement techniques for siloxanes and other trace compounds in biogas. Waste Manag. 30 1011 - 1017
Ohannessian A , Desjardin V , Chatain V , Germain P. 2008 Volatile organic silicon compounds: the most undesirable contaminants in biogases. Water Sci. Technol. 58 1775 - 1781
Dewil R , Appels L , Baeyens J , Buczynska A , Van Vaeck L. 2007 The analysis of volatile siloxanes in waste activated sludge. Talanta 74 14 - 19
Ajhar M , Travesset M , Yüce S , Melin T. 2010 Siloxane removal from landfill and digester gas - a technology overview. Bioresour. Technol. 101 2913 - 2923
Matsui T , Imamura S. 2010 Removal of siloxane from digestion gas of sewage sludge. Bioresour. Technol. 101 (Suppl 1) S29 - 32
Schweigkofler M , Niessner R. 2001 Removal of siloxanes in biogases. J. Hazard. Mater. 83 183 - 196
Hong C , Haiyun W. 2010 Optimization of volatile fatty acid production with co-substrate of food wastes and dewatered excess sludge using response surface methodology. Bioresour. Technol. 101 5487 - 5493
Montgomery DC. 2005 Design and analysis of experiments. 6th ed. John Wiley & Sons Hoboken
Molinuevo-Salces B , García-González MC , González-Fernández C , Cuetos MJ , Morán A , Gómez X. 2010 Anaerobic co-digestion of livestock wastes with vegetable processing wastes: a statistical analysis. Bioresour. Technol. 101 9479 - 9485
Trinh TK , Kang LS. 2010 Application of response surface method as an experimental design to optimize coagulation tests. Environ. Eng. Res. 15 63 - 70
Ricaurte Ortega D , Subrenat A. 2009 Siloxane treatment by adsorption into porous materials. Environ. Technol. 30 1073 - 1083
Kim JS , Kim MC , Kang EJ , Kim MS. 2002 H2S adsorption characteristics of KIO3impregnated activated carbon. J. Korean Oil Chem. Soc. 20 72 - 79
An SW , Yoo JY , Choi JY , Park JW. 2009 Adsorption characterization of Cd by activated carbon containing hydroxyapatite using response surface methodology (RSM). J. Korean Soc. Water Qual. 25 943 - 950