Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Response surface methodology (RSM) and its application for optimization of ammonium ions removal from aqueous solutions by pumice as a natural and low cost adsorbent

Response surface methodology (RSM) and its application for optimization of ammonium ions removal... This research was conducted to study the adsorption of ammonium ions onto pumice as a natural and low-cost adsorbent. The physico-chemical properties of the pumice granular were characterized by X-ray diffraction (XRD), Fourier transforms infrared spectroscopy (FTIR) and scanning electron microscopy (SEM). Modeling and optimization of a NH4+ sorption process was accomplished by varying four independent parameters (pumice dosage, initial ammonium ion concentration, mixing rate and contact time) using a central composite design (CCD) under response surface methodology (RSM). The optimum conditions for maximum removal of NH4+ (70.3%) were found to be 100 g, 20 mg/l, 300 rpm and 180 min, for pumice dosage, initial NH4+ ion concentration, mixing rate and contact time. It was found that the NH4+ adsorption on the pumice granular was dependent on adsorbent dosage and initial ammonium ion concentration. NH4+ was increased due to decrease the initial concentration of NH4 and increase the contact time, mixing rate and amount of adsorbent. Introduction The presence of nitrate and ammonium ions in effluents is a major concern due to their extreme toxicity, oxygen depletion in receiving water bodies and promotion of eutrophication (Arslan and Veli 2012). For these circumstances, it is necessary to develop efficient solutions to remove these pollutants from wastewater ensuring its reuse and/or decreasing the pollution to the acceptable standards for discharging into the receiving bodies of water. There are several different approaches based on the mechanical and physicochemical methods for removing ammonium ions from wastewaters, such as biological nitrification-denitrification, chemical treatment, precipitation, air stripping, selective ion exchange and adsorption (Vassileva and Voikova 2009). Among these techniques, adsorption seems to be one of the economical and the most effective methods because of simple operation and easy handling (Akbal 2005, Khosravi et al. 2014, Pirsaheb et al. 2014). Adsorption is the most versatile process and widely used for the pollutants removal from wastewaters. Several types of ammonia adsorbents have been tested so far, but most of them remain on a laboratorialscale application (Bandosz and Petit 2009). In recent years, the treatment of the ammonia found in liquid or gaseous effluents by using adsorption onto activated carbon or other adsorbent material has been studied by several researchers. Yuan and Kusuda (Yuan and Kusuda 2005) used a hydrogel of poly (N-isopropylacrylamide-co-chlorophyllin) to adsorb NH4- N and found that in a swollen state, the developed hydrogel had an adsorption capacity of 12.89 mg/g for NH4+N. Karapinar (Karapinar 2009) reported that zeolite can be used as an adsorbent for the removal of ammonium and phosphorus from the wastewaters, and ammonium uptake of zeolite was about 0.4 mg/g for initial NH4+ concentration of 40 mg/L and zeolite dosage of 8 g/L. Demir and coworkers (Demir et al. 2002) studied the capacity of Bigadiç clinoptilolite for removing ions from synthetic aqueous solutions under varying conditions. The actual ion-exchange capacity of the conditioned fine and coarse clinoptilolite was found to be 0.57 and 0.38 mg/g zeolite, respectively. Another study was carried out on adsorption behavior of NH4+-N from aqueous solution by activated carbons prepared from rice husk. The results obtained from this study Analysis and data calculation A Fourier transformed infrared (FTIR) was carried out using a spectrometer (WQF-510) in the transmittance mode with a resolution of 4 cm-1 in the range of 400­4000 cm-1 with KBr pellets technique. FTIR is an instrument that indicates the main functional groups of adsorbent and these groups play the main role in adsorbate removal from solutions. So that, any wave in the range of 400­4000 cm-1 indicates the particular functional group (Bardakçi and Bahçeli 2010). The chemical characteristics of the pumice samples were determined by a XRD (Shimadzu XRD-6000). The NH4+ ion concentration in the aqueous phase was determined using the standard Nesslerre agent method employing a Perkin-Elmer 402 UV spectrophotometer. A surface morphology of pumice granular was studied by scanning electron microscopy using a Philips XL30 microscope at an accelerating voltage of 10 kV. After oven-drying of the pumice granular for 12 hours, the sample was coated with a platinum layer using an SCDOOS sputter coater (BAL-TEC, Sweden) in an argon atmosphere. Subsequently, the sample was scanned and photomicrographs were obtained. Adsorption experiment Adsorption experiments were carried out in 500 mL Erlenmeyer flasks by using five different concentrations of NH4+ , mixing rate at temperature of 23 ±2°C and pH=7. A certain amount of pumice granular was mixed with 250 mL ammonium solutions of concentrations ranging from 20 to 100 mg/L. The amount of NH4+ adsorbed on the pumice granular was calculated from the difference between the initial concentration of NH4+ and the residual NH4+ concentration in the solution. The calculation can be described by the following equation (Wang et al. 2011): indicated that the developed activated carbon had an adsorption capacity of 1.3 mg/g for NH4+-N (Zhu et al. 2012). Commercial activated carbon and zeolites are the most common adsorbents that have been implemented in the process of ammonium ion removal. These products are relatively effective for removal of pollutants but are expensive. In recent years, studies have been devoted to low-cost materials such as local mineral sorbents for the elimination of pollutants from water. These sorbents can be used in natural or modified forms. One such low-cost material is pumice. Pumice is an amorphous, porous volcanic rock which is composed mainly of SiO2 and whose porous structure is formed by dissolved gases precipitated during the cooling as the lava hurtles through the air (Gunduz et al. 1998, Lura et al. 2004). After Italy, Greece and China, Iran is the fourth largest pumice producer in the world (Bolen 2003). Due to its micro porous structure, pumice has a high specific surface area and low apparent density (0.35­0.65 g/cm3) in aggregate form. Recently, many researchers used pumice for removal of cadmium (Panuccio et al. 2009), heavy metals (Moraci and Calabrò 2010), sulfur dioxide (Ozturk and Yildirim 2008) and disinfection by-products (Bekaroglu et al. 2010). However, based on our knowledge, it is seldom reported that the pumice stone is used in removal of NH4+ ion. Response surface methodology (RSM) is an effective statistical technique and also a critical tool that is used for modeling and analyzing the effects of multiple parameters on the process (Ba and Boyaci 2007, Mansouri et al. 2012). The main advantages of this method are the reduced number of experimental trials, calculating the complex interaction between the independent variables, analysis and optimization as well as the improvement of existing design. In other words, this statistical approach is more practical compared to conventional experimental work as it arises from experimental methodology which includes interactive effects among the variables and, finally, it shows the overall effects of the variables on the process (Box and Draper 1987). Also, there are many classes of response surface designs that are occasionally useful in practice, such as Central Composite design, Box-Behnken design, Hybrid design and Three-level Factorial design (Khuri and Cornell 1996, Mason et al. 2003). So, due to the several advantages of pumice stone and its accessibility in Iran, the aim of the present work was to investigate its effectiveness for the removal of NH4+ ion from aqueous solution. The Central composite design (CCD) and RSM were used to design the experiments, build models and determine the optimum conditions. qe (C 0 -C e )V (1) where qe (mg/g) is the adsorbent capacity, C0 (mg/L) is the initial concentration of NH4+, Ce (mg/L) is the final or equilibrium concentration of NH4+, V is the experimental solution volume (L), and M is the weight of pumice granular (g). Experimental design RSM along with the Design of Experiments remove systematic errors with an estimate of the experimental error and also reduce the number of experiments to obtain the optimum operating conditions. There are many classes of response surface designs that are occasionally useful in practice such as Central Composite design, Box-Behnken design, Hybrid design and Three-level Factorial design (Khuri and Cornell 1996, Mason et al. 2003). Among these methods the CCD is the most frequently used one under RSM. The RSM employed in the current study was a CCD involving four different factors; pumice dosage , NH4+ concentration , mixing rate (C) and contact time (D). The range investigated for these variables is presented in Table 1. The NH4+ removal efficiency and adsorption capacity of pumice granular were evaluated based on the CCD experimental plan (Table 2). The design consisted of 2k factorial points augmented by 2k axial points and a center point where k is the number of variables. The Materials and methods Materials Natural pumice sample used in this study was obtained from the Kourdistan (Ghorve) of Iran and has a density of 0.64 g cm-3. Prior to use in the experiments, pumice samples were crushed and sieved through a 65 × 65 mesh sieve (based on Taylor standard seive) to obtained the particles size less than 0.227 mm. NH4Cl supplied from Merck was used as adsorbate. 1000 ppm NH4Cl solution was prepared as a stock solution and then it was diluted to the required initial concentrations. Synthetic samples of aqueous solution were prepared by adding appropriate amounts of NH4Cl stock solution to deionized water to obtain ion concentrations of 20, 40, 60, 80 and 100 mg/l, respectively. coded values of four operating variables (A, B, C and D) were set at five levels: -1.5 (minimum), -1, 0 (central), +1, and +1.5 (maximum). These five levels were assessed based on the full face-centered CCD experimental plan. Accordingly, 30 experiments (=2k +2k +6), where k is the number of factors, were conducted with 25 experiments organized in a factorial design (including 16 factorial points, 8 axial points and 1 center point) and the remaining 5 involving the replication of the central point to get a good estimate of the experimental error. Repetition experiments were carried out followed by order of runs designed by Design of Experiments (DOE) as shown in Table 2. Table 1. Experimental range and levels of the independent variables Variables Pumice dosage, g/L Initial ammonium concentration, mg/L Mixing rate Contact time (min) Range and levels - (-1.5) 20 20 100 60 -1 40 40 150 90 0 60 60 200 120 1 80 80 250 150 + (1.5) 100 100 300 180 Table 2. Experimental conditions and results of central composite design Variables Factor 1 Run A: Pumice dosage g/l 60 60 20 60 100 100 60 60 20 100 20 100 100 60 20 20 80 100 60 100 40 100 60 20 60 20 60 60 60 20 Factor 2 B: Initial NH4+ Conc. mg/l 60 60 100 60 100 20 60 60 100 100 20 20 100 60 20 20 60 100 60 20 60 20 60 20 60 100 40 80 60 100 Factor 3 C: Mixing speed rpm 200 200 300 200 100 300 200 200 300 300 300 100 300 200 100 100 200 100 250 100 200 300 200 300 200 100 200 200 150 100 Factor 4 D: Contact time min 120 120 180 120 60 60 120 150 60 60 60 180 180 120 180 60 120 180 120 60 120 180 90 180 120 180 120 120 120 60 Actual % 45.3 40.2 24.12 37.6 32.77 66.25 43.1 47.69 20.7 37.88 37.79 61.3 44.5 43.6 39.5 32.64 51.2 39.5 46.3 57.18 29.88 71.4 35.13 43.9 39.4 19.5 41.67 31.16 39.7 18.3 % 40.64 40.64 26.33 40.64 33.38 64.69 40.64 42.05 20.70 39.33 38.48 64.37 44.97 40.64 38.16 32.52 46.24 39.01 42.13 58.73 35.03 70.32 39.23 44.12 40.64 20.37 46.03 35.25 39.15 14.74 Responses qe (180 min), mg/g Actual mg/g 0.140 0.125 0.330 0.115 0.090 0.050 0.120 0.133 0.273 0.110 0.117 0.040 0.114 0.117 0.100 0.105 0.114 0.113 0.140 0.037 0.120 0.040 0.102 0.110 0.106 0.270 0.080 0.128 0.108 0.246 mg/g 0.13 0.13 0.28 0.13 0.091 0.046 0.13 0.13 0.28 0.11 0.11 0.026 0.11 0.13 0.089 0.089 0.097 0.091 0.13 0.026 0.16 0.046 0.13 0.11 0.13 0.26 0.097 0.13 0.12 0.26 (moisture) which was adsorbed by the sample from the outside environment (Grim 1968). The water adsorbed from the outside environment also produces a deformation vibration peak around 1620 cm-1 (Ersoy et al. 2010). SEM images detailing the micro-structure (i.e., inner pore structure) of the pumice are displayed in Fig. 2. According to SEM images, the pumice sample has irregular or oval shaped and fibrous cavities (or pores), in addition, it may be said that these pores are either in closed or in open forms. The pore diameter of these cavities varied between 1 m and 30 m and these cavities generally did not intersect each other. It is reported that the connectivity of the pore structure of pumice may range from completely closed to completely open (Lura et al. 2004). Moreover, the glassy phase structure of the pumice can be seen in Fig. 2. According to the element analysis of the XRD, the biggest peaks were recorded as related to Si and Al. Model fitting and statistical analysis The design matrix in the term of coded factors and experimental results of the NH4+ removal and equilibrium adsorption capacity are presented in Table 2. The obtained results were then, subjected to RSM to evaluate the relationship between the pumice dosage , NH4+ concentration , mixing rate (C) and contact time (D). The values obtained from model fitting technique using the design expert software for both responses were seen to be sufficiently correlated to the observed values (Table 2). The response values slightly deviated from the experimental data. Multiple regression coefficients of a second-order polynomial models describing the NH4+ removal and equilibrium adsorption capacity are summarized in Table 2. For determination of coefficients' significance and the suitability of proposed models, the ANOVA technique was performed. The significance or insignificant of each coefficient was determined by P-value. The P-values lower than 0.05 indicated that coefficient of a model term is significant. The resulting ANOVA are listed in Table 3 for each of the model terms. Based on the P-values, the first order effects of all variables (A, B, D and C) and two level interaction of A and B (AB) are significant for both process responses. Other model terms are not significant (with probability values larger than Mathematical modeling After conducting the experiments, the relationship between the dependent and independent variables was calculated using the following equation (Box and Draper 1987): Y = 0 + i X i + j X j + ii X i + jj X j + ij X i X j + ... (2) Where Y, i, j, ß, X represent process response, linear coefficient, quadratic coefficient, regression coefficient and coded independent variables respectively. Model terms are accepted or rejected based on the probability of error (P) value with 95% confidence level. The results obtained from CCD were entirely examined by means of analysis of variance (ANOVA) by Design Expert software. Three-dimensional (3D) plots and their respective contour plots were obtained due to the effect of the two factors while other factors can be changed by default. Therefore, the results of CCD can be presented in 3D presentations with contours. This will help to study the simultaneous interaction of the variables on the responses. Results and discussion Pumice characterization The chemical analysis of pumice granular samples was determined using XRD. The results showed the composition as SiO2 (67%), Al2O3 (15.6%), K2O (5.4%), Fe2O3 (1.1%), Na2O (3.7%), CaO (3.9%), MgO (0.2%) and etc. As seen from the values, the two most important components of the sample are SiO2 and Al2O3. FTIR spectrum is the feature of a particular compound that gives the information about its functional groups, molecular geometry and inter/intra molecular interactions. Fig. 1 shows the FTIR spectra of pumice granular. The peaks at 467and 798 cm-1 in the FTIR spectrum of the pumice (Fig. 1) may have resulted from the Si­O bending strength vibrations of the SiO2 which constituted the structure. The strong peak around 1070 cm-1 may have resulted from Si­O stretching vibrations (Ersoy et al. 2010). The other important peak was detected around 3406 cm-1. This peak pointed out the OH stretching vibrations of the water Fig. 1. FTIR spectrum of the pumice Fig. 2. SEM images of the pumice sample scaled up by ×15 times, scaled up by×300 times 0.05). In order to simplify the model, these insignificant model terms (A2, B2, C2, D2, AC, AD, BC and BD) were eliminated. The analysis of variance (ANOVA) was employed to examine the significance of the obtained models under an extensive range of operating conditions. The results are depicted in Table 4 (numerically) and Figs. 3­5 (graphically). The ANOVA results revealed that the quadratic model was significant at the 5% confidence level respectively for NH4+ removal and equilibrium adsorption capacity since P values were less than 0.05 for both models. Fisher F-test value explains the distribution of the actual data around the fitted model. According to the calculated F-test values (86.66 for NH4+ removal and 92.38 for equilibrium adsorption capacity) along with a very low probability value for both responses (p model < 0.0001), it can be concluded that both models are highly significant. The fitting of models to the empirical data was tested by calculating determination coefficient (R2). The high values of R2 (0.95 for NH4+ removal and 0.94 for equilibrium adsorption capacity) designate the capability of developed models to satisfactorily describe the system behavior within the investigated range of operating parameters. Adequate Precision (AP) compares the range of the values at the design points to the average prediction error. Ratios greater than 4 indicate adequate model discrimination and can be used to navigate the design space defined by the CCD (Mason et al. 2003). In this study, the ratios of AP value were found to be 39.32 and 33.6 for NH4+ removal and equilibrium adsorption capacity respectively, which are greater than 4 and considered satisfied results. Other important information about fitting, reliability, adequacy, the correlation between the observed and data, and homogeneous and heterogeneous error variances of the models performance can be attained in the diagnostic plots (Figs. 3­5). These figures give an obvious sight about any deficiency of the models fitting to the experimental data Figs. 3a and 3b display the normal probability of the residuals for both responses, to confirm whether the standard deviations between the actual and the response values follow a normal distribution (Lee and Gilmore 2005). Points and points' clusters in Figs. 3a and 3b indicate that experimental values are distributed relatively near to the straight line and show satisfactory correlation between these values. Therefore, there are no serious violations in the hypothesis that errors are normally distributed and independent of each other, and also that the error variances are homogeneous and residuals are independent (Kusic et al. 2011). The prediction of Amunium ion removal efficincy and qe (Eqs. 2 and 3) was compared with the experimental values given in Tab. 2 and shown in Fig. 4. It can be seen from Fig4 that in the both models equations predictions are statistically a satisfactory match with the experimental values. Therefore, this model can be used to navigate the design space. In fact, the figures 4a and 4b depict the relationship between the values of R2 and adjusted R2. The R2 values of 0.95 and 0.94 are in reasonable agreement with the Adj R2 values of 0.93 and 0.91 for the NH4+ removal and equilibrium adsorption capacity, which shows a very good agreement between the and actual data. The plots of residual against responses (Figs. 5a and 5b) check the assumption of constant variance. It can be seen that all points of experimental runs are randomly distributed and all values lie within the range of -3 and 3. These findings indicate that the models proposed by RSM are satisfactory and that the constant variance assumptions were confirmed. Interactive effect of the variables on the NH4+ removal efficiency RSM is applied to build up mathematical relationships between the process variables and the responses of interest, using the simulation results. RSM is one of the best known empirical modeling techniques, which helps in carrying out the analysis of experiments with the least experimental effort. From the ANOVA indicated in Table 4, the main order effects of all variables (A, B, C and D) and the two level interaction between pumice dosage and initial concentration of NH4+ (AB) are the significant model terms on the both process responses. ANOVA results also indicate that the effects of A and B do not depend Table 3. Estimated regression coefficients and corresponding to ANOVA results from the data of central composite design experiments before elimination of insignificant model terms Model variables C B2 -15.82 -0.038 8.24 0.038 41.64 2.8 *10-4 0.011 57.27 0.011 5.07 41.59 0.6 0.3 Not Not Not 0.04 < 0.0001 0.36 Not 0.89 10.10 18.43 1.63 0.2 1.18 0.26 1.3 *10-3 2.88 41.64 1.3 *10-2 3.69 1.01 0.06 0.33 Not Not 3.2 *10-4 1.4 *10-3 1.58 1.4 *10-3 0.14 1.28 0.71 0.3 Not Not 18.43 2.88 57.27 10.10 1.58 0.038 0.004 0.004 8.40 8.40 0.74 0.40 Not sig 8.24 8.24 0.84 0.84 0.84 0.84 0.84 0.004 1.14 1.3 *10-4 1.14 1.3 *10-3 0.1 1.22 0.75 0.28 Not Not 0.86 Not sig 0.042 -0.027 -0.005 -0.005 10.52 4.16 -1.89 0.79 0.31 -0.72 -0.27 0.15 0.84 0.36 0.36 0.03 C2 D2 AB AC AD BC BD CD 2.98 0.83 146.41 146.41 12.96 0.003 Not Not sig 0.03 0.72 0.0003 0.93 5.59 0.13 11.60 0.007 1.2 *10-2 1.3 *10-4 130.99 0.007 1.5 *10-3 3.4 *10-5 146.41 0.08 0.004 0.038 0.83 8.24 0.001 0.014 2.82 0.68 D A2 Coefficient Estimate -10.78 -0.058 Standard Error Sum of squares (SS) Mean square (MS) F-value Prob> F < 0.0001 < 0.0001 < 0.0001 < 0.0001 Significant or insignificant model termes Table 4. Analysis of variance (ANOVA) for fit of Ammonium ion removal (%) and qe(mg/g) from central composite design after elimination of insignificant model terms Model Type of model Significant model terms F value Prob > F SD R2 Adj. R2 Pred. R2 Mean Probability for Adeq. PRESS lack of fit precision 39.3 309.15 0.4520 Ammonium Quadratic A, B, C, D, AB ion removal, model % Quadratic model A, B, C, AB < 0.0001 3.16 0.95 < 0.0001 0.018 0.94 R2: determination coefficient, Adj. R2: adjusted R2, Adeq. precision: adequate precision, SD: standard deviation, PRESS: residual error sum of squares Normal % Probability Normal % Probability -1.67 -0.80 -1.97 -0.68 Fig. 3. Normal probability plots of residual for NH4+ removal (%), qe (mg/g) Actual Actual Fig. 4. vs. actual values plots for NH4+ removal (%), qe (mg/g) on the level of factors C and D, whereas factors A and B do affect each other. Hence A and B are said to be interacting. The following regression equation is the reduced quadratic models in terms of implied factors for the both responses: Based on the results obtained from ANVA (Table 4), the response surface models constructed in this study for predicting NH4+­N removal efficiency and equilibrium adsorption capacity were considered reasonable. The final regression models, in terms of their coded factors as determined by the Design Experts software for both responses, are expressed by the following equations. Amonium ion removal efficincy (%)= =+40.64+11.21A-10.78B+2.98C +2.82D-1.89 AB (2) -1.50 -1.45 -3.00 -3.00 Run Number Fig. 5. Plots of residual versus responses for NH4+removal (%) and qe (mg/g) B: Ammonium ion conc. mg/l40 D: Contact time, min Mixing speed, rpm Fig. 6. Response surface plots for NH4+ removal (%) with respect to: . pumice dosage and ammonium ion concentration, mixing speed (rpm) and contact time(min) =+0.13­0.058 A+0.059 B+0.1 C­0.027 AB (3) The relationship between the responses and variables was visualized by the three-dimensional surface plots to see the influence of the parameters. Fig. 6a shows threedimensional contour plots of the model for variation in NH4+ removal, as a function of pumice dosage and initial NH4+ concentration at constant values of mixing rate (200 rpm) and contact time (60 min). It is clear from this figure that one of the most significant variables on the response was pumice dosage. The adsorption of NH4+ increased with increase in pumice dosage, presumably from the greater availability of active binding sites and to the presence of a greater surface area for adsorption. But the amount of NH4+ ions adsorbed per unit mass of adsorbent decreased with increase in adsorbent dose (Fig 6a). The phenomenon that the reduction of the adsorption unit of NH4 -N on the pumice granular was mainly due to the fact that the total adsorption sites increased with the increase of the dosage of pumice granular, but the total amount of NH4 -N in the solution still remained the same. As shown in Fig. 6a, an increase in the initial ammonium concentration from 20 to 100 mg/L caused an increase in the response. Increasing the initial NH4+ ions concentration would increase the mass transfer driving force and therefore the rate at which NH4+ ions pass from the solution to the particle surface. Figure 6b shows the three-dimensional plot of the model for NH4+ removal with respect to mixing rate and contact time the design space, with pumice dosage at its high level (100 g/L) and initial ammonium concentration at its low level (20 mg/l). As shown in this figure, the contact time (range 60­120 min) did not show any strong effect on the process, as only a small difference is observed in the results obtained with different contact time. However, the improvement of NH4+ removal efficiency (about 6%) was observed with increase of contact time from 60 to 120 min. Therefore, contact time of 60 min was chosen as optimal parameter in the work. These results also show that NH4+ adsorption by pumice granular was very fast at the initial period of contact time and then to become slow with the increase of contact time. This was caused by strong attractive forces between the NH4+ ions and the sorbent; fast diffusion onto the external surface was followed by a fast pore diffusion into the intra particle matrix to attain rapid equilibrium. As shown in Fig. 6b, an increase in the mixing rate increases the NH4+ removal. It can be stated that increase the mixing rate cause increases the contact between adsorbent and adsorbate, this event improved the access of NH4+ ions to the active binding sites of pumice granular as a result the NH4 increases. Perturbation plot (Fig. 7) also shows the comparative effects of the variables on NH4+ removal efficiency. In Fig. 7, steep lines in pumice dosage and ammonium concentration show that the response of NH4+ removal efficiency was very sensitive to these factors. The maximum modeled NH4+ removal was 70% at pumice dosage, initial NH4+ concentration and mixing rate of 100 g/l, 20 mg/l and 300 rpm, respectively. Meanwhile, the minimum response (20.6%) was obtained at pumice dosage, initial NH4+ concentration and mixing rate of 20 g/l, 100 mg/l and 100 rpm, respectively. Fig. 8 shows response surface and contour plots of the model for variation in qe, as a function of pumice dosage and initial ammonium concentration with three different mixing rates (100, 200 and 300). As can be seen in the figures (8a­c), the same trends were found as the mixing rate changed from 100 to 300. In Fig. 8, the response increased upon increasing the NH4+ concentration and decreasing the pumice dosage. It is demonstrated by perturbation plot (Fig. 9). The highest qe was attained when ammonium ion concentration and mixing rate were high, while the pumice dosage was low as can be seen from Fig. 8. It must be noted that qe values were calculated based on contact time of 180 min. It was observed that the equilibrium adsorption capacity decreased with increasing the amount of pumice dosage from 20 to 100 g/L. The optimum for the maximum equilibrium adsorption capacity was found to be 0.28 mg/g. According to obtained results, increases the adsorbent dosage can be attributed to increase the pumice surface area and the availability of more adsorption sites, nevertheless, the values of equilibrium adsorption capacity decreased with an increase in the adsorbent dosage. The primary factor explaining this characteristic is that adsorption sites remain unsaturated DC CD -0.973 -0.466 (c) Fig. 8. Response surface plots for qe (mg/g) with respect topumice dosage and ammonium ion concentration at constant value of mixing rate of 100, 200 and (c) 300 rpm Dev iation f rom Ref erence Point Fig. 7. Perturbation plot for NH4+ removal Bekaroglu, S.K., Yigit, N., Karanfil, T. & Kitis, M. (2010) The adsorptive removal of disinfection by-product precursors in a high-SUVA water using iron oxide-coated pumice and volcanic slag particles, Journal of Hazardous Materials, 183(1), pp. 389­394. Bolen, W.P. (2003) Pumice and pumicite. US Geological Survey Minerals Yearbook. Box, G.E. & Draper, N.R. (1987) Empirical model-building and response surfaces, Wiley New York 1987. Demir, A., Gunay, A. & Debik, E. (2002) Ammonium removal from aqueous solution by ion-exchange using packed bed natural zeolite, Water SA, 28(3), pp. 329­336. Ersoy, B., Sariisik, A., Dikmen, S. & Sariisik, G. (2010) Characterization of acidic pumice and determination of its electrokinetic properties in water, Powder Technology, 197(1), pp. 129­135. Grim, R.E. (1968) Clay mineralogy. McGraw-Hill Book Company, New York, 596. Gunduz, L., Sariisik, A., Tozacan, B., Davraz, M., Ugur, I. & Cankiran, O. (1998) Pumice technology. Skin 1. Karapinar, N. (2009) Application of natural zeolite for phosphorus and ammonium removal from aqueous solutions, Journal of hazardous materials, 170(2), pp. 1186­1191. Khosravi, R., Fazlzadehdavil, M., Barikbin, B. & Taghizadeh, A.A. (2014) Removal of hexavalent chromium from aqueous solution by granular and powdered Peganum Harmala, Applied Surface Science, 292, pp. 670­677. Khuri, A.I. & Cornell, J.A. (1996) Response surfaces: designs and analyses, CRC press 1996. Kusic, H., Koprivanac, N. & Bozic, A.L. (2011) Treatment of chlorophenols in water matrix by UV/ferrioxalate system: Part I. Key process parameter evaluation by response surface methodology, Desalination, 279(1), pp. 258­268. Lee, K.-M. & Gilmore, D.F. (2005) Formulation and process modeling of biopolymer (polyhydroxyalkanoates: PHAs) production from industrial wastes by novel crossed experimental design, Process Biochemistry, 40(1), pp. 229­246. Lura, P., Bentz, D.P., Lange, D.A., Kovler, K. & Bentur, A. (2004) Pumice aggregates for internal water curing, pp. 22­24, RILEM Publications SARL Evanston. Mansouri, Y., Zinatizadeh, A.A., Mohammadi, P., Irandoust, M., Akhbari, A. & Davoodi, R. (2012) Hydraulic characteristics analysis of an anaerobic rotatory biological contactor (AnRBC) using tracer experiments and response surface methodology (RSM), Korean Journal of Chemical Engineering, 29(7), pp. 891­902. Mason, R.L., Gunst, R.F. & Hess, J.L. (2003) Statistical design and analysis of experiments: with applications to engineering and science, John Wiley & Sons 2003. Moraci, N. & Calabrò, P.S. (2010) Heavy metals removal and hydraulic performance in zero-valent iron/pumice permeable reactive barriers, Journal of Environmental Management, 91(11), pp. 2336­2341. Ozturk, B. & Yildirim, Y. (2008) Investigation of sorption capacity of pumice for SO 2 capture, Process Safety and Environmental Protection, 86(1), pp. 31­36. Panuccio, M.R., Sorgonà, A., Rizzo, M. & Cacco, G. (2009) Cadmium adsorption on vermiculite, zeolite and pumice: batch experimental studies, Journal of Environmental Management, 90(1), pp. 364­374. Pirsaheb, M., Dargahi, A., Hazrati, S. & Fazlzadehdavil, M. (2014) Removal of diazinon and 2, 4-dichlorophenoxyacetic acid (2, 4-D) from aqueous solutions by granular-activated carbon, Desalination and Water Treatment, 52(22­24), pp. 4350­4355. Vassileva, P. & Voikova, D. (2009) Investigation on natural and pretreated Bulgarian clinoptilolite for ammonium ions removal from aqueous solutions, Journal of hazardous materials, 170(2), pp. 948­953. during the adsorption reaction whereas the number of sites available for adsorption increased by increasing the adsorbent dosage. This means that the higher ratio of the ammonium ion concentration to the surface binding site on the sorbent could be obtained at lower pumice dosage and higher initial NH4+ ion concentration. DC CD -0.973 -0.466 Dev iation f rom Ref erence Point Fig. 9. Perturbation plot for qe ( mg/g) Conclusion This study has demonstrated the possibility of response surface methodology application for studying the effect of process factors on removal of NH4+ ions from aqueous solution by pumice granular. Based on the response surface plots, among the studied variables, pumice dosage and initial ammonium concentration significantly affected the NH4+ removal. The maximum NH4+ was modeled to be 70.3% under pumice dosage, initial NH4+ ion concentration, mixing rate and contact time of 100 g/l, 20 mg/l, 300 rpm and 180 min, respectively, while the observed value was 71.4%. Based on the equilibrium study, the optimum for the maximum equilibrium adsorption capacity was found to be 0.28 mg/g. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Archives of Environmental Protection de Gruyter

Response surface methodology (RSM) and its application for optimization of ammonium ions removal from aqueous solutions by pumice as a natural and low cost adsorbent

Loading next page...
 
/lp/de-gruyter/response-surface-methodology-rsm-and-its-application-for-optimization-xov7UurDJi
Publisher
de Gruyter
Copyright
Copyright © 2016 by the
ISSN
2083-4810
eISSN
2083-4810
DOI
10.1515/aep-2016-0018
Publisher site
See Article on Publisher Site

Abstract

This research was conducted to study the adsorption of ammonium ions onto pumice as a natural and low-cost adsorbent. The physico-chemical properties of the pumice granular were characterized by X-ray diffraction (XRD), Fourier transforms infrared spectroscopy (FTIR) and scanning electron microscopy (SEM). Modeling and optimization of a NH4+ sorption process was accomplished by varying four independent parameters (pumice dosage, initial ammonium ion concentration, mixing rate and contact time) using a central composite design (CCD) under response surface methodology (RSM). The optimum conditions for maximum removal of NH4+ (70.3%) were found to be 100 g, 20 mg/l, 300 rpm and 180 min, for pumice dosage, initial NH4+ ion concentration, mixing rate and contact time. It was found that the NH4+ adsorption on the pumice granular was dependent on adsorbent dosage and initial ammonium ion concentration. NH4+ was increased due to decrease the initial concentration of NH4 and increase the contact time, mixing rate and amount of adsorbent. Introduction The presence of nitrate and ammonium ions in effluents is a major concern due to their extreme toxicity, oxygen depletion in receiving water bodies and promotion of eutrophication (Arslan and Veli 2012). For these circumstances, it is necessary to develop efficient solutions to remove these pollutants from wastewater ensuring its reuse and/or decreasing the pollution to the acceptable standards for discharging into the receiving bodies of water. There are several different approaches based on the mechanical and physicochemical methods for removing ammonium ions from wastewaters, such as biological nitrification-denitrification, chemical treatment, precipitation, air stripping, selective ion exchange and adsorption (Vassileva and Voikova 2009). Among these techniques, adsorption seems to be one of the economical and the most effective methods because of simple operation and easy handling (Akbal 2005, Khosravi et al. 2014, Pirsaheb et al. 2014). Adsorption is the most versatile process and widely used for the pollutants removal from wastewaters. Several types of ammonia adsorbents have been tested so far, but most of them remain on a laboratorialscale application (Bandosz and Petit 2009). In recent years, the treatment of the ammonia found in liquid or gaseous effluents by using adsorption onto activated carbon or other adsorbent material has been studied by several researchers. Yuan and Kusuda (Yuan and Kusuda 2005) used a hydrogel of poly (N-isopropylacrylamide-co-chlorophyllin) to adsorb NH4- N and found that in a swollen state, the developed hydrogel had an adsorption capacity of 12.89 mg/g for NH4+N. Karapinar (Karapinar 2009) reported that zeolite can be used as an adsorbent for the removal of ammonium and phosphorus from the wastewaters, and ammonium uptake of zeolite was about 0.4 mg/g for initial NH4+ concentration of 40 mg/L and zeolite dosage of 8 g/L. Demir and coworkers (Demir et al. 2002) studied the capacity of Bigadiç clinoptilolite for removing ions from synthetic aqueous solutions under varying conditions. The actual ion-exchange capacity of the conditioned fine and coarse clinoptilolite was found to be 0.57 and 0.38 mg/g zeolite, respectively. Another study was carried out on adsorption behavior of NH4+-N from aqueous solution by activated carbons prepared from rice husk. The results obtained from this study Analysis and data calculation A Fourier transformed infrared (FTIR) was carried out using a spectrometer (WQF-510) in the transmittance mode with a resolution of 4 cm-1 in the range of 400­4000 cm-1 with KBr pellets technique. FTIR is an instrument that indicates the main functional groups of adsorbent and these groups play the main role in adsorbate removal from solutions. So that, any wave in the range of 400­4000 cm-1 indicates the particular functional group (Bardakçi and Bahçeli 2010). The chemical characteristics of the pumice samples were determined by a XRD (Shimadzu XRD-6000). The NH4+ ion concentration in the aqueous phase was determined using the standard Nesslerre agent method employing a Perkin-Elmer 402 UV spectrophotometer. A surface morphology of pumice granular was studied by scanning electron microscopy using a Philips XL30 microscope at an accelerating voltage of 10 kV. After oven-drying of the pumice granular for 12 hours, the sample was coated with a platinum layer using an SCDOOS sputter coater (BAL-TEC, Sweden) in an argon atmosphere. Subsequently, the sample was scanned and photomicrographs were obtained. Adsorption experiment Adsorption experiments were carried out in 500 mL Erlenmeyer flasks by using five different concentrations of NH4+ , mixing rate at temperature of 23 ±2°C and pH=7. A certain amount of pumice granular was mixed with 250 mL ammonium solutions of concentrations ranging from 20 to 100 mg/L. The amount of NH4+ adsorbed on the pumice granular was calculated from the difference between the initial concentration of NH4+ and the residual NH4+ concentration in the solution. The calculation can be described by the following equation (Wang et al. 2011): indicated that the developed activated carbon had an adsorption capacity of 1.3 mg/g for NH4+-N (Zhu et al. 2012). Commercial activated carbon and zeolites are the most common adsorbents that have been implemented in the process of ammonium ion removal. These products are relatively effective for removal of pollutants but are expensive. In recent years, studies have been devoted to low-cost materials such as local mineral sorbents for the elimination of pollutants from water. These sorbents can be used in natural or modified forms. One such low-cost material is pumice. Pumice is an amorphous, porous volcanic rock which is composed mainly of SiO2 and whose porous structure is formed by dissolved gases precipitated during the cooling as the lava hurtles through the air (Gunduz et al. 1998, Lura et al. 2004). After Italy, Greece and China, Iran is the fourth largest pumice producer in the world (Bolen 2003). Due to its micro porous structure, pumice has a high specific surface area and low apparent density (0.35­0.65 g/cm3) in aggregate form. Recently, many researchers used pumice for removal of cadmium (Panuccio et al. 2009), heavy metals (Moraci and Calabrò 2010), sulfur dioxide (Ozturk and Yildirim 2008) and disinfection by-products (Bekaroglu et al. 2010). However, based on our knowledge, it is seldom reported that the pumice stone is used in removal of NH4+ ion. Response surface methodology (RSM) is an effective statistical technique and also a critical tool that is used for modeling and analyzing the effects of multiple parameters on the process (Ba and Boyaci 2007, Mansouri et al. 2012). The main advantages of this method are the reduced number of experimental trials, calculating the complex interaction between the independent variables, analysis and optimization as well as the improvement of existing design. In other words, this statistical approach is more practical compared to conventional experimental work as it arises from experimental methodology which includes interactive effects among the variables and, finally, it shows the overall effects of the variables on the process (Box and Draper 1987). Also, there are many classes of response surface designs that are occasionally useful in practice, such as Central Composite design, Box-Behnken design, Hybrid design and Three-level Factorial design (Khuri and Cornell 1996, Mason et al. 2003). So, due to the several advantages of pumice stone and its accessibility in Iran, the aim of the present work was to investigate its effectiveness for the removal of NH4+ ion from aqueous solution. The Central composite design (CCD) and RSM were used to design the experiments, build models and determine the optimum conditions. qe (C 0 -C e )V (1) where qe (mg/g) is the adsorbent capacity, C0 (mg/L) is the initial concentration of NH4+, Ce (mg/L) is the final or equilibrium concentration of NH4+, V is the experimental solution volume (L), and M is the weight of pumice granular (g). Experimental design RSM along with the Design of Experiments remove systematic errors with an estimate of the experimental error and also reduce the number of experiments to obtain the optimum operating conditions. There are many classes of response surface designs that are occasionally useful in practice such as Central Composite design, Box-Behnken design, Hybrid design and Three-level Factorial design (Khuri and Cornell 1996, Mason et al. 2003). Among these methods the CCD is the most frequently used one under RSM. The RSM employed in the current study was a CCD involving four different factors; pumice dosage , NH4+ concentration , mixing rate (C) and contact time (D). The range investigated for these variables is presented in Table 1. The NH4+ removal efficiency and adsorption capacity of pumice granular were evaluated based on the CCD experimental plan (Table 2). The design consisted of 2k factorial points augmented by 2k axial points and a center point where k is the number of variables. The Materials and methods Materials Natural pumice sample used in this study was obtained from the Kourdistan (Ghorve) of Iran and has a density of 0.64 g cm-3. Prior to use in the experiments, pumice samples were crushed and sieved through a 65 × 65 mesh sieve (based on Taylor standard seive) to obtained the particles size less than 0.227 mm. NH4Cl supplied from Merck was used as adsorbate. 1000 ppm NH4Cl solution was prepared as a stock solution and then it was diluted to the required initial concentrations. Synthetic samples of aqueous solution were prepared by adding appropriate amounts of NH4Cl stock solution to deionized water to obtain ion concentrations of 20, 40, 60, 80 and 100 mg/l, respectively. coded values of four operating variables (A, B, C and D) were set at five levels: -1.5 (minimum), -1, 0 (central), +1, and +1.5 (maximum). These five levels were assessed based on the full face-centered CCD experimental plan. Accordingly, 30 experiments (=2k +2k +6), where k is the number of factors, were conducted with 25 experiments organized in a factorial design (including 16 factorial points, 8 axial points and 1 center point) and the remaining 5 involving the replication of the central point to get a good estimate of the experimental error. Repetition experiments were carried out followed by order of runs designed by Design of Experiments (DOE) as shown in Table 2. Table 1. Experimental range and levels of the independent variables Variables Pumice dosage, g/L Initial ammonium concentration, mg/L Mixing rate Contact time (min) Range and levels - (-1.5) 20 20 100 60 -1 40 40 150 90 0 60 60 200 120 1 80 80 250 150 + (1.5) 100 100 300 180 Table 2. Experimental conditions and results of central composite design Variables Factor 1 Run A: Pumice dosage g/l 60 60 20 60 100 100 60 60 20 100 20 100 100 60 20 20 80 100 60 100 40 100 60 20 60 20 60 60 60 20 Factor 2 B: Initial NH4+ Conc. mg/l 60 60 100 60 100 20 60 60 100 100 20 20 100 60 20 20 60 100 60 20 60 20 60 20 60 100 40 80 60 100 Factor 3 C: Mixing speed rpm 200 200 300 200 100 300 200 200 300 300 300 100 300 200 100 100 200 100 250 100 200 300 200 300 200 100 200 200 150 100 Factor 4 D: Contact time min 120 120 180 120 60 60 120 150 60 60 60 180 180 120 180 60 120 180 120 60 120 180 90 180 120 180 120 120 120 60 Actual % 45.3 40.2 24.12 37.6 32.77 66.25 43.1 47.69 20.7 37.88 37.79 61.3 44.5 43.6 39.5 32.64 51.2 39.5 46.3 57.18 29.88 71.4 35.13 43.9 39.4 19.5 41.67 31.16 39.7 18.3 % 40.64 40.64 26.33 40.64 33.38 64.69 40.64 42.05 20.70 39.33 38.48 64.37 44.97 40.64 38.16 32.52 46.24 39.01 42.13 58.73 35.03 70.32 39.23 44.12 40.64 20.37 46.03 35.25 39.15 14.74 Responses qe (180 min), mg/g Actual mg/g 0.140 0.125 0.330 0.115 0.090 0.050 0.120 0.133 0.273 0.110 0.117 0.040 0.114 0.117 0.100 0.105 0.114 0.113 0.140 0.037 0.120 0.040 0.102 0.110 0.106 0.270 0.080 0.128 0.108 0.246 mg/g 0.13 0.13 0.28 0.13 0.091 0.046 0.13 0.13 0.28 0.11 0.11 0.026 0.11 0.13 0.089 0.089 0.097 0.091 0.13 0.026 0.16 0.046 0.13 0.11 0.13 0.26 0.097 0.13 0.12 0.26 (moisture) which was adsorbed by the sample from the outside environment (Grim 1968). The water adsorbed from the outside environment also produces a deformation vibration peak around 1620 cm-1 (Ersoy et al. 2010). SEM images detailing the micro-structure (i.e., inner pore structure) of the pumice are displayed in Fig. 2. According to SEM images, the pumice sample has irregular or oval shaped and fibrous cavities (or pores), in addition, it may be said that these pores are either in closed or in open forms. The pore diameter of these cavities varied between 1 m and 30 m and these cavities generally did not intersect each other. It is reported that the connectivity of the pore structure of pumice may range from completely closed to completely open (Lura et al. 2004). Moreover, the glassy phase structure of the pumice can be seen in Fig. 2. According to the element analysis of the XRD, the biggest peaks were recorded as related to Si and Al. Model fitting and statistical analysis The design matrix in the term of coded factors and experimental results of the NH4+ removal and equilibrium adsorption capacity are presented in Table 2. The obtained results were then, subjected to RSM to evaluate the relationship between the pumice dosage , NH4+ concentration , mixing rate (C) and contact time (D). The values obtained from model fitting technique using the design expert software for both responses were seen to be sufficiently correlated to the observed values (Table 2). The response values slightly deviated from the experimental data. Multiple regression coefficients of a second-order polynomial models describing the NH4+ removal and equilibrium adsorption capacity are summarized in Table 2. For determination of coefficients' significance and the suitability of proposed models, the ANOVA technique was performed. The significance or insignificant of each coefficient was determined by P-value. The P-values lower than 0.05 indicated that coefficient of a model term is significant. The resulting ANOVA are listed in Table 3 for each of the model terms. Based on the P-values, the first order effects of all variables (A, B, D and C) and two level interaction of A and B (AB) are significant for both process responses. Other model terms are not significant (with probability values larger than Mathematical modeling After conducting the experiments, the relationship between the dependent and independent variables was calculated using the following equation (Box and Draper 1987): Y = 0 + i X i + j X j + ii X i + jj X j + ij X i X j + ... (2) Where Y, i, j, ß, X represent process response, linear coefficient, quadratic coefficient, regression coefficient and coded independent variables respectively. Model terms are accepted or rejected based on the probability of error (P) value with 95% confidence level. The results obtained from CCD were entirely examined by means of analysis of variance (ANOVA) by Design Expert software. Three-dimensional (3D) plots and their respective contour plots were obtained due to the effect of the two factors while other factors can be changed by default. Therefore, the results of CCD can be presented in 3D presentations with contours. This will help to study the simultaneous interaction of the variables on the responses. Results and discussion Pumice characterization The chemical analysis of pumice granular samples was determined using XRD. The results showed the composition as SiO2 (67%), Al2O3 (15.6%), K2O (5.4%), Fe2O3 (1.1%), Na2O (3.7%), CaO (3.9%), MgO (0.2%) and etc. As seen from the values, the two most important components of the sample are SiO2 and Al2O3. FTIR spectrum is the feature of a particular compound that gives the information about its functional groups, molecular geometry and inter/intra molecular interactions. Fig. 1 shows the FTIR spectra of pumice granular. The peaks at 467and 798 cm-1 in the FTIR spectrum of the pumice (Fig. 1) may have resulted from the Si­O bending strength vibrations of the SiO2 which constituted the structure. The strong peak around 1070 cm-1 may have resulted from Si­O stretching vibrations (Ersoy et al. 2010). The other important peak was detected around 3406 cm-1. This peak pointed out the OH stretching vibrations of the water Fig. 1. FTIR spectrum of the pumice Fig. 2. SEM images of the pumice sample scaled up by ×15 times, scaled up by×300 times 0.05). In order to simplify the model, these insignificant model terms (A2, B2, C2, D2, AC, AD, BC and BD) were eliminated. The analysis of variance (ANOVA) was employed to examine the significance of the obtained models under an extensive range of operating conditions. The results are depicted in Table 4 (numerically) and Figs. 3­5 (graphically). The ANOVA results revealed that the quadratic model was significant at the 5% confidence level respectively for NH4+ removal and equilibrium adsorption capacity since P values were less than 0.05 for both models. Fisher F-test value explains the distribution of the actual data around the fitted model. According to the calculated F-test values (86.66 for NH4+ removal and 92.38 for equilibrium adsorption capacity) along with a very low probability value for both responses (p model < 0.0001), it can be concluded that both models are highly significant. The fitting of models to the empirical data was tested by calculating determination coefficient (R2). The high values of R2 (0.95 for NH4+ removal and 0.94 for equilibrium adsorption capacity) designate the capability of developed models to satisfactorily describe the system behavior within the investigated range of operating parameters. Adequate Precision (AP) compares the range of the values at the design points to the average prediction error. Ratios greater than 4 indicate adequate model discrimination and can be used to navigate the design space defined by the CCD (Mason et al. 2003). In this study, the ratios of AP value were found to be 39.32 and 33.6 for NH4+ removal and equilibrium adsorption capacity respectively, which are greater than 4 and considered satisfied results. Other important information about fitting, reliability, adequacy, the correlation between the observed and data, and homogeneous and heterogeneous error variances of the models performance can be attained in the diagnostic plots (Figs. 3­5). These figures give an obvious sight about any deficiency of the models fitting to the experimental data Figs. 3a and 3b display the normal probability of the residuals for both responses, to confirm whether the standard deviations between the actual and the response values follow a normal distribution (Lee and Gilmore 2005). Points and points' clusters in Figs. 3a and 3b indicate that experimental values are distributed relatively near to the straight line and show satisfactory correlation between these values. Therefore, there are no serious violations in the hypothesis that errors are normally distributed and independent of each other, and also that the error variances are homogeneous and residuals are independent (Kusic et al. 2011). The prediction of Amunium ion removal efficincy and qe (Eqs. 2 and 3) was compared with the experimental values given in Tab. 2 and shown in Fig. 4. It can be seen from Fig4 that in the both models equations predictions are statistically a satisfactory match with the experimental values. Therefore, this model can be used to navigate the design space. In fact, the figures 4a and 4b depict the relationship between the values of R2 and adjusted R2. The R2 values of 0.95 and 0.94 are in reasonable agreement with the Adj R2 values of 0.93 and 0.91 for the NH4+ removal and equilibrium adsorption capacity, which shows a very good agreement between the and actual data. The plots of residual against responses (Figs. 5a and 5b) check the assumption of constant variance. It can be seen that all points of experimental runs are randomly distributed and all values lie within the range of -3 and 3. These findings indicate that the models proposed by RSM are satisfactory and that the constant variance assumptions were confirmed. Interactive effect of the variables on the NH4+ removal efficiency RSM is applied to build up mathematical relationships between the process variables and the responses of interest, using the simulation results. RSM is one of the best known empirical modeling techniques, which helps in carrying out the analysis of experiments with the least experimental effort. From the ANOVA indicated in Table 4, the main order effects of all variables (A, B, C and D) and the two level interaction between pumice dosage and initial concentration of NH4+ (AB) are the significant model terms on the both process responses. ANOVA results also indicate that the effects of A and B do not depend Table 3. Estimated regression coefficients and corresponding to ANOVA results from the data of central composite design experiments before elimination of insignificant model terms Model variables C B2 -15.82 -0.038 8.24 0.038 41.64 2.8 *10-4 0.011 57.27 0.011 5.07 41.59 0.6 0.3 Not Not Not 0.04 < 0.0001 0.36 Not 0.89 10.10 18.43 1.63 0.2 1.18 0.26 1.3 *10-3 2.88 41.64 1.3 *10-2 3.69 1.01 0.06 0.33 Not Not 3.2 *10-4 1.4 *10-3 1.58 1.4 *10-3 0.14 1.28 0.71 0.3 Not Not 18.43 2.88 57.27 10.10 1.58 0.038 0.004 0.004 8.40 8.40 0.74 0.40 Not sig 8.24 8.24 0.84 0.84 0.84 0.84 0.84 0.004 1.14 1.3 *10-4 1.14 1.3 *10-3 0.1 1.22 0.75 0.28 Not Not 0.86 Not sig 0.042 -0.027 -0.005 -0.005 10.52 4.16 -1.89 0.79 0.31 -0.72 -0.27 0.15 0.84 0.36 0.36 0.03 C2 D2 AB AC AD BC BD CD 2.98 0.83 146.41 146.41 12.96 0.003 Not Not sig 0.03 0.72 0.0003 0.93 5.59 0.13 11.60 0.007 1.2 *10-2 1.3 *10-4 130.99 0.007 1.5 *10-3 3.4 *10-5 146.41 0.08 0.004 0.038 0.83 8.24 0.001 0.014 2.82 0.68 D A2 Coefficient Estimate -10.78 -0.058 Standard Error Sum of squares (SS) Mean square (MS) F-value Prob> F < 0.0001 < 0.0001 < 0.0001 < 0.0001 Significant or insignificant model termes Table 4. Analysis of variance (ANOVA) for fit of Ammonium ion removal (%) and qe(mg/g) from central composite design after elimination of insignificant model terms Model Type of model Significant model terms F value Prob > F SD R2 Adj. R2 Pred. R2 Mean Probability for Adeq. PRESS lack of fit precision 39.3 309.15 0.4520 Ammonium Quadratic A, B, C, D, AB ion removal, model % Quadratic model A, B, C, AB < 0.0001 3.16 0.95 < 0.0001 0.018 0.94 R2: determination coefficient, Adj. R2: adjusted R2, Adeq. precision: adequate precision, SD: standard deviation, PRESS: residual error sum of squares Normal % Probability Normal % Probability -1.67 -0.80 -1.97 -0.68 Fig. 3. Normal probability plots of residual for NH4+ removal (%), qe (mg/g) Actual Actual Fig. 4. vs. actual values plots for NH4+ removal (%), qe (mg/g) on the level of factors C and D, whereas factors A and B do affect each other. Hence A and B are said to be interacting. The following regression equation is the reduced quadratic models in terms of implied factors for the both responses: Based on the results obtained from ANVA (Table 4), the response surface models constructed in this study for predicting NH4+­N removal efficiency and equilibrium adsorption capacity were considered reasonable. The final regression models, in terms of their coded factors as determined by the Design Experts software for both responses, are expressed by the following equations. Amonium ion removal efficincy (%)= =+40.64+11.21A-10.78B+2.98C +2.82D-1.89 AB (2) -1.50 -1.45 -3.00 -3.00 Run Number Fig. 5. Plots of residual versus responses for NH4+removal (%) and qe (mg/g) B: Ammonium ion conc. mg/l40 D: Contact time, min Mixing speed, rpm Fig. 6. Response surface plots for NH4+ removal (%) with respect to: . pumice dosage and ammonium ion concentration, mixing speed (rpm) and contact time(min) =+0.13­0.058 A+0.059 B+0.1 C­0.027 AB (3) The relationship between the responses and variables was visualized by the three-dimensional surface plots to see the influence of the parameters. Fig. 6a shows threedimensional contour plots of the model for variation in NH4+ removal, as a function of pumice dosage and initial NH4+ concentration at constant values of mixing rate (200 rpm) and contact time (60 min). It is clear from this figure that one of the most significant variables on the response was pumice dosage. The adsorption of NH4+ increased with increase in pumice dosage, presumably from the greater availability of active binding sites and to the presence of a greater surface area for adsorption. But the amount of NH4+ ions adsorbed per unit mass of adsorbent decreased with increase in adsorbent dose (Fig 6a). The phenomenon that the reduction of the adsorption unit of NH4 -N on the pumice granular was mainly due to the fact that the total adsorption sites increased with the increase of the dosage of pumice granular, but the total amount of NH4 -N in the solution still remained the same. As shown in Fig. 6a, an increase in the initial ammonium concentration from 20 to 100 mg/L caused an increase in the response. Increasing the initial NH4+ ions concentration would increase the mass transfer driving force and therefore the rate at which NH4+ ions pass from the solution to the particle surface. Figure 6b shows the three-dimensional plot of the model for NH4+ removal with respect to mixing rate and contact time the design space, with pumice dosage at its high level (100 g/L) and initial ammonium concentration at its low level (20 mg/l). As shown in this figure, the contact time (range 60­120 min) did not show any strong effect on the process, as only a small difference is observed in the results obtained with different contact time. However, the improvement of NH4+ removal efficiency (about 6%) was observed with increase of contact time from 60 to 120 min. Therefore, contact time of 60 min was chosen as optimal parameter in the work. These results also show that NH4+ adsorption by pumice granular was very fast at the initial period of contact time and then to become slow with the increase of contact time. This was caused by strong attractive forces between the NH4+ ions and the sorbent; fast diffusion onto the external surface was followed by a fast pore diffusion into the intra particle matrix to attain rapid equilibrium. As shown in Fig. 6b, an increase in the mixing rate increases the NH4+ removal. It can be stated that increase the mixing rate cause increases the contact between adsorbent and adsorbate, this event improved the access of NH4+ ions to the active binding sites of pumice granular as a result the NH4 increases. Perturbation plot (Fig. 7) also shows the comparative effects of the variables on NH4+ removal efficiency. In Fig. 7, steep lines in pumice dosage and ammonium concentration show that the response of NH4+ removal efficiency was very sensitive to these factors. The maximum modeled NH4+ removal was 70% at pumice dosage, initial NH4+ concentration and mixing rate of 100 g/l, 20 mg/l and 300 rpm, respectively. Meanwhile, the minimum response (20.6%) was obtained at pumice dosage, initial NH4+ concentration and mixing rate of 20 g/l, 100 mg/l and 100 rpm, respectively. Fig. 8 shows response surface and contour plots of the model for variation in qe, as a function of pumice dosage and initial ammonium concentration with three different mixing rates (100, 200 and 300). As can be seen in the figures (8a­c), the same trends were found as the mixing rate changed from 100 to 300. In Fig. 8, the response increased upon increasing the NH4+ concentration and decreasing the pumice dosage. It is demonstrated by perturbation plot (Fig. 9). The highest qe was attained when ammonium ion concentration and mixing rate were high, while the pumice dosage was low as can be seen from Fig. 8. It must be noted that qe values were calculated based on contact time of 180 min. It was observed that the equilibrium adsorption capacity decreased with increasing the amount of pumice dosage from 20 to 100 g/L. The optimum for the maximum equilibrium adsorption capacity was found to be 0.28 mg/g. According to obtained results, increases the adsorbent dosage can be attributed to increase the pumice surface area and the availability of more adsorption sites, nevertheless, the values of equilibrium adsorption capacity decreased with an increase in the adsorbent dosage. The primary factor explaining this characteristic is that adsorption sites remain unsaturated DC CD -0.973 -0.466 (c) Fig. 8. Response surface plots for qe (mg/g) with respect topumice dosage and ammonium ion concentration at constant value of mixing rate of 100, 200 and (c) 300 rpm Dev iation f rom Ref erence Point Fig. 7. Perturbation plot for NH4+ removal Bekaroglu, S.K., Yigit, N., Karanfil, T. & Kitis, M. (2010) The adsorptive removal of disinfection by-product precursors in a high-SUVA water using iron oxide-coated pumice and volcanic slag particles, Journal of Hazardous Materials, 183(1), pp. 389­394. Bolen, W.P. (2003) Pumice and pumicite. US Geological Survey Minerals Yearbook. Box, G.E. & Draper, N.R. (1987) Empirical model-building and response surfaces, Wiley New York 1987. Demir, A., Gunay, A. & Debik, E. (2002) Ammonium removal from aqueous solution by ion-exchange using packed bed natural zeolite, Water SA, 28(3), pp. 329­336. Ersoy, B., Sariisik, A., Dikmen, S. & Sariisik, G. (2010) Characterization of acidic pumice and determination of its electrokinetic properties in water, Powder Technology, 197(1), pp. 129­135. Grim, R.E. (1968) Clay mineralogy. McGraw-Hill Book Company, New York, 596. Gunduz, L., Sariisik, A., Tozacan, B., Davraz, M., Ugur, I. & Cankiran, O. (1998) Pumice technology. Skin 1. Karapinar, N. (2009) Application of natural zeolite for phosphorus and ammonium removal from aqueous solutions, Journal of hazardous materials, 170(2), pp. 1186­1191. Khosravi, R., Fazlzadehdavil, M., Barikbin, B. & Taghizadeh, A.A. (2014) Removal of hexavalent chromium from aqueous solution by granular and powdered Peganum Harmala, Applied Surface Science, 292, pp. 670­677. Khuri, A.I. & Cornell, J.A. (1996) Response surfaces: designs and analyses, CRC press 1996. Kusic, H., Koprivanac, N. & Bozic, A.L. (2011) Treatment of chlorophenols in water matrix by UV/ferrioxalate system: Part I. Key process parameter evaluation by response surface methodology, Desalination, 279(1), pp. 258­268. Lee, K.-M. & Gilmore, D.F. (2005) Formulation and process modeling of biopolymer (polyhydroxyalkanoates: PHAs) production from industrial wastes by novel crossed experimental design, Process Biochemistry, 40(1), pp. 229­246. Lura, P., Bentz, D.P., Lange, D.A., Kovler, K. & Bentur, A. (2004) Pumice aggregates for internal water curing, pp. 22­24, RILEM Publications SARL Evanston. Mansouri, Y., Zinatizadeh, A.A., Mohammadi, P., Irandoust, M., Akhbari, A. & Davoodi, R. (2012) Hydraulic characteristics analysis of an anaerobic rotatory biological contactor (AnRBC) using tracer experiments and response surface methodology (RSM), Korean Journal of Chemical Engineering, 29(7), pp. 891­902. Mason, R.L., Gunst, R.F. & Hess, J.L. (2003) Statistical design and analysis of experiments: with applications to engineering and science, John Wiley & Sons 2003. Moraci, N. & Calabrò, P.S. (2010) Heavy metals removal and hydraulic performance in zero-valent iron/pumice permeable reactive barriers, Journal of Environmental Management, 91(11), pp. 2336­2341. Ozturk, B. & Yildirim, Y. (2008) Investigation of sorption capacity of pumice for SO 2 capture, Process Safety and Environmental Protection, 86(1), pp. 31­36. Panuccio, M.R., Sorgonà, A., Rizzo, M. & Cacco, G. (2009) Cadmium adsorption on vermiculite, zeolite and pumice: batch experimental studies, Journal of Environmental Management, 90(1), pp. 364­374. Pirsaheb, M., Dargahi, A., Hazrati, S. & Fazlzadehdavil, M. (2014) Removal of diazinon and 2, 4-dichlorophenoxyacetic acid (2, 4-D) from aqueous solutions by granular-activated carbon, Desalination and Water Treatment, 52(22­24), pp. 4350­4355. Vassileva, P. & Voikova, D. (2009) Investigation on natural and pretreated Bulgarian clinoptilolite for ammonium ions removal from aqueous solutions, Journal of hazardous materials, 170(2), pp. 948­953. during the adsorption reaction whereas the number of sites available for adsorption increased by increasing the adsorbent dosage. This means that the higher ratio of the ammonium ion concentration to the surface binding site on the sorbent could be obtained at lower pumice dosage and higher initial NH4+ ion concentration. DC CD -0.973 -0.466 Dev iation f rom Ref erence Point Fig. 9. Perturbation plot for qe ( mg/g) Conclusion This study has demonstrated the possibility of response surface methodology application for studying the effect of process factors on removal of NH4+ ions from aqueous solution by pumice granular. Based on the response surface plots, among the studied variables, pumice dosage and initial ammonium concentration significantly affected the NH4+ removal. The maximum NH4+ was modeled to be 70.3% under pumice dosage, initial NH4+ ion concentration, mixing rate and contact time of 100 g/l, 20 mg/l, 300 rpm and 180 min, respectively, while the observed value was 71.4%. Based on the equilibrium study, the optimum for the maximum equilibrium adsorption capacity was found to be 0.28 mg/g.

Journal

Archives of Environmental Protectionde Gruyter

Published: Jun 1, 2016

References