The Scientific World Journal

The Scientific World Journal / 2014 / Article

Research Article | Open Access

Volume 2014 |Article ID 869120 |

Anam Asghar, Abdul Aziz Abdul Raman, Wan Mohd Ashri Wan Daud, "A Comparison of Central Composite Design and Taguchi Method for Optimizing Fenton Process", The Scientific World Journal, vol. 2014, Article ID 869120, 14 pages, 2014.

A Comparison of Central Composite Design and Taguchi Method for Optimizing Fenton Process

Academic Editor: Seval Kutlu Akal Solmaz
Received21 Apr 2014
Revised22 Jul 2014
Accepted06 Aug 2014
Published27 Aug 2014


In the present study, a comparison of central composite design (CCD) and Taguchi method was established for Fenton oxidation. , Dye : Fe+2, H2O2 : Fe+2, and pH were identified control variables while COD and decolorization efficiency were selected responses. orthogonal array and face-centered CCD were used for the experimental design. Maximum 99% decolorization and 80% COD removal efficiency were obtained under optimum conditions. squared values of 0.97 and 0.95 for CCD and Taguchi method, respectively, indicate that both models are statistically significant and are in well agreement with each other. Furthermore, Prob > less than 0.0500 and ANOVA results indicate the good fitting of selected model with experimental results. Nevertheless, possibility of ranking of input variables in terms of percent contribution to the response value has made Taguchi method a suitable approach for scrutinizing the operating parameters. For present case, pH with percent contribution of 87.62% and 66.2% was ranked as the most contributing and significant factor. This finding of Taguchi method was also verified by 3D contour plots of CCD. Therefore, from this comparative study, it is concluded that Taguchi method with 9 experimental runs and simple interaction plots is a suitable alternative to CCD for several chemical engineering applications.

1. Introduction

Fenton oxidation is an efficient and widely applied AOPs, which utilizes ferrous iron (Fe+2) and hydrogen peroxide (H2O2) under acidic conditions to produce hydroxyl radical (HO) (reaction 1). Formation of HO is advantageous because it is highly oxidative () and nonselective in nature [1, 2]. However, the performance of Fenton reaction is influenced by several parameters such as pH, reaction time, initial concentrations of Dye, H2O2, and Fe+2 catalysts [3, 4]. Moreover, excess consumption of chemicals (Fe+2 salt, H2O2 acid/base) and high cost of H2O2 have made this process economically nonviable [5, 6]. Therefore, a systematic way of planning, execution, and statistical evaluation of process is required which is only possible through optimization process [7]. Consider In context to any engineering problem, optimization refers to improving the performance of system or process by applying several variables in different combinations to get the best possible result [8]. Optimization techniques are broadly classified into two categories: univariate and multivariate approach. Univariate approach, also known as one-factor-at-time (OFAT) approach, involves variation of one parameter at a time. Nevertheless, multivariate approach has several advantages over OFAT. For example [9], (1) it provides global knowledge in its whole experimental domain while OFAT gives local knowledge where experiment is performed; (2) it is possible to study the interaction between the factors and nonlinear relationships with the responses, (3) the number of experiments required to optimize the process is considerably lesser than that of OFAT approach; (4) within the experimental domain at each point, quality of the information is higher and can be known through the leverage. The stepwise execution of optimization is presented in Figure 1.

Driven by the need of reducing the number of experiments, cost, time, and physical efforts, design of experiment (DOE) is an important statistical and mathematical tool for solving complex and multifactor engineering problems. It includes response surface methodology (RSM), factorial design, and in some special cases artificial neural network (ANN) [10]. However, central composite design (CCD), D-optimal, and Box-Behnken are found to be widely used optimization techniques for Fenton oxidation [1113] because of the advantage of optimizing multifactor problems with optimum number of experimental runs. Besides, these methods have limitation of increased number of experiments if several factors were selected for process optimization. For example, Box-Behnken suggests 54 experiments if six factors are required to be studied and with CCD this number increases up to 80. Thus, with multiple factors these techniques are not appropriate because it increases the cost of chemicals, time, and physical efforts. Therefore, a simplified design strategy is required that can be used to overcome these problems.

Taguchi method is a robust statistical tool that allows the independent evaluation of the responses with minimum number of experiments. It employs orthogonal arrays for experimental design and ratio instead of responses itself to determine the optimum settings of control factors and thus neglects the variations caused by uncontrollable factors [14]. With this method, experimental results can be analyzed through ratio and ANOVA with simultaneously evaluating the significance of the factors in terms of their contribution to the response values.

Available works on optimization of Fenton oxidation mainly focus on CCD under RSM [1517]. Nevertheless, fewer studies based on Taguchi method are found in the literature [1820] but which technique is better for Fenton oxidation is still not inclusive. Therefore, this study was planned to compare two experimental design techniques and Fenton oxidation with four parameters such as , Dye : Fe+2, H2O2 : Fe+2, and pH was selected as a case study. For comparison, frequently used CCD and robust Taguchi experimental design technique were chosen with a special focus on comparing the optimization and detailed statistical analysis of the experimental results.

2. Background

2.1. Response Surface Methodology and Central Composite Design

RSM, a multivariate statistical tool, consists of a group of mathematical and statistical techniques that are based on the fit of empirical models to the experimental data obtained in relation to experimental design [21]. It employs lower order polynomial [22] and it has already been proved to be a reliable statistical method for chemical process applications [23, 24].

In RSM category, CCD which is appropriate for fitting second order polynomial equations has been frequently discussed for optimizing several research problems. A CCD has three groups of design points [25]:(a)two-level factorial or fractional factorial design points (), consisting of possible combinations of +1 and −1 levels of factor;(b) axial points (sometimes called star points) fixed axially at a distance say from the center to generate quadratic terms;(c)center points which represent replicate terms; center points provide a good and independent estimate of the experimental error.

Considering these points, the number of experiments designed by CCD will be where is the total number of experiments, is the number of factors studied, and is the number of replicates. Central composite design under RSM is normally performed either by using design expert software or Minitab. In this study, design expert (Version: is used as optimization software. The steps that will be followed for the central composite design (CCD) are presented in Figure 2.

In CCD, value of alpha is important to calculate as it could determine the location of axial points in experimental domain. Depending on alpha value, design is spherical, orthogonal, rotatable, or face centered. Practically, it is in between face centered and spherical and is calculated as Value of alpha equals 1 is desirable because it ensures the position of axial point within factorial portion region. It is called face centered design and offers three levels for the factors to be put in the experimental design matrix.

Experimental results obtained are analyzed using response surface regression procedure of statistical analysis system. Corelation between responses and independent variables is obtained by fitting them into second order polynomial equation [10] Here, represents the responses, is the total number independent factors, is an intercept, , , and with represent the coefficient values for linear, quadratic, and interaction effects, respectively, and and in the above equation show the coded levels for independent variables [26].

2.2. Taguchi Method and Signal to Noise Ratio

Taguchi method is a robust and multiparameter optimization statistical technique which employs fewer numbers of experiments to identify and optimize parameters to achieve desired response [27, 28]. It is a simple and systematic way of determining the effect of factors on responses and optimum condition of the factors. Taguchi method utilizes full fractional design called orthogonal arrays and ANOVA as a tool for analysis [29]. Orthogonal arrays are the minimum set of experiments which represents the various combinations of factors. Output of the orthogonal arrays is optimized with respect to signal to noise ratio () of the responses instead of the responses itself and thus it reduces the process variability [30]. This feature marks the difference between the conventional statistical technique and Taguchi method. The stepwise procedure for Taguchi method is represented in Figure 3.

In Taguchi method, ratio is the measure of the deviation of the response from the desired value. Here, “signal” implies the mean value while “noise” shows the standard deviation term. It means that lower variability in the process is ensured through maximizing the ratio. However, depending on the type of response desired, Taguchi classified ratio into three categories: smaller-the-better, larger-the-better, and nominal-the-better [31]:Smaller is better: Larger is better: Nominal is better: where ;  .

And represents response variables and “” donates the number of experiments.

One of the distinct features of Taguchi method is that it determines optimum value in the form of ratio. The predicted ratio at optimal process conditions can be computed by the following mathematical equation [32]: where shows the mean of all ratios, is the ratio at the optimal level for each parameters, and is the number of process parameters that significantly affect the process.

3. Case Study

In this study, Fenton oxidation of synthetic dye acid blue 113 has been taken as a case study for comparing two optimization techniques, that is, CCD and Taguchi method.

3.1. Experimental Details

For performing Fenton oxidation, all reagents used were of analytical grade and procured from Merck Sdn Bhd, Malaysia. Synthetic dye acid blue 113 was purchased from Sigma-Aldrich (M) Sdn Bhd, Malaysia.

For performing Fenton oxidation, stock solutions of 33.3 g/L of H2O2 and 100 g/L of FeSO4·7H2O were prepared. All solutions were prepared by using distilled water. Fenton oxidation experiments were performed in batch laboratory scale Erlenmeyer flask (500 mL capacity) equipped with a magnetic stirrer. For each test, 150 mL of acid blue 113 dye solution was placed into flask and pH was adjusted by using 0.5 M H2SO4 or 1 M NaOH solution according to the designed experiment. Desired amounts of FeSO4·7H2O and H2O2 solutions were added to the dye solution, respectively, while stirring at 150 rpm. The start time was recorded when H2O2 solution was added. After 90 minutes, pH of the solution was measured and sufficient amount of 1 M NaOH solution was added to increase pH above 11 to stop the reaction. This is because in alkaline conditions Fe+2 precipitates out as Fe+3 and H2O2 decomposes into oxygen. Then, after 2 hrs of settling time, the supernatant was filtered through 0.45 μm Millipore filter and subjected to various analyses. All experiments were performed at room temperature and atmospheric pressure.

3.2. Analysis

pH measurements of the solution were carried out with the aid of pH meter (CyberScan pH 300, Eutectic instruments). Quantitative measurement of the residual H2O2 in treated wastewater was carried out by using peroxide test strips (Merck). The raw and treated samples were scanned using UV-spectrophotometer (Spectroquant Pharo 300, Merck). And decolorization efficiency of the treated sample was calculated as follows: COD measurements were made according to the standard method (APHA, AWWA, and WFE, 1998). For this, COD test cell supplied by Merck was heated in thermoreactor (Spectroquant TR 420) after adding required amount of the sample followed by the subsequent measurement in UV-spectrophotometer according to the standard method. And COD removal efficiency is calculated as

4. Result and Discussion

4.1. Central Composite Design

All experiments for Fenton oxidation were designed according to RSM using central composite design (CCD) with the aid of design expert (Version: According to the literature review, the most important parameters, which affect the efficiency of Fenton process, are pH, initial concentrations of Dye, Fe+2, and H2O2 [4, 33, 34]. Therefore, , pH, H2O2 : Fe+2 (wt/wt), and Dye : Fe+2 (wt/wt) ratios were chosen as control variables to be optimized using CCD as given in Table 1.

Independent numerical variablesCodedLow actual value High actual value

Dye (mg/L)100300
H2O2 : Fe+2 (wt/wt)525
Dye : Fe+2 (wt/wt)1050

Total 30 runs with 16 factorial, 8 Axial, and 6 center points were suggested by design expert to optimize the responses, that is, COD and decolorization efficiency. Based on the proposed model, the following quadratic equation was developed to predict the dependent variables (responses) in terms of independent variables and their interactions: where is the response variable (COD removal or decolorization efficiency), is constant, , , , and are coefficient for linear effects, , , , and are quadratic coefficient, and , , , , , and are interaction coefficients [36], respectively.

4.2. Model Results for Fenton Oxidation of Acid Blue 113
4.2.1. Model Fitting and Analysis of Variance (ANOVA)

Different concentrations of C.I. acid blue dye 113 were prepared and Fenton oxidation was performed under different conditions according to the experimental runs as suggested by CCD model. Based on experimental results, the following empirical second order polynomial equations were developed showing the interactions between the proposed independent variables to obtain decolorization and COD removal efficiencies: The objective of this empirical model was to adequately describe the interaction of factors influencing the process efficiency at the concentration ranges investigated. Experimental and predicted values on COD removal and decolorization efficiencies are provided in Table 2. The observed COD removal and decolorization values vary between 30.2–87.7% and 56.1–99.23% which are in good agreement with the predicted values as shown in Figure 4.

Independent variables ()Dependent variables (; (%))
Actual valuesPredicted values
RunDye (mg/L)H2O2 : Fe+2 (wt/wt)Dye : Fe+2 (wt/wt)pHCOD (%)Decolorization (%)COD (%)Decolorization (%)


Mathematical equation developed after fitting the function to the data may give misleading results and cannot describe the domain of the model adequately [37]. That is why ANOVA is an integral part of the data analysis and is the more reliable way to evaluate the quality of the model fitted [21]. Table 3 shows the analysis of variance for COD and color removal efficiency.

SourceSum of squaresDegree of freedomMean square valueProb >  

COD removal
 H2O2 : Fe+2188.481188.4813.940.0012
 Dye : Fe+2168.361168.3614.950.002
Color Removal
 H2O2 : Fe+2631.551631.55102<0.0001
 Dye : Fe+238.90138.906.280.0242
 pH1331.9711331.97 215.12<0.0001

Values of Prob > less than 0.0500 imply that the model terms are significant while values greater than 0.1000 are demonstrated as insignificant for the regression model. In the current study, values 32.41 and 49.73 for COD and color removal denote the model as significant. The goodness of the fit of the model is also checked by value. In both cases, the values for this regression coefficient were 0.9789 and 0.9581 which implies that this model is statistically significant and is in reasonable agreement with the adjusted . Moreover, “Adeq precision” is used to determine the signal to noise () ratio to determine the validity of the model. A ratio greater than 4 is recommended. In our case, values of 25.9 and 22.3 for color and COD removal indicate an adequate signal. They also show that this model can be used to navigate the design space. Another checkpoint for validating the experimental data is the analysis of normal probability plots. The normal probability plot indicates whether the residuals follow a normal distribution, in which case the points will follow a straight line, as it is in the current study (Figures A.1 and A.2 in Supplementary Material available online at Considering the above explained ANOVA results, it can be concluded that this model explained the Fenton reaction and can be employed to navigate the design space in terms of decolorization and COD removal efficiencies.

4.2.2. Response Surface Plotting and Optimization of Operating Parameters

Graphical interpretation of interactions by using three and two dimensional plots of regression model is highly recommended and is used to assess the interactive effects between the process variables and treatment efficiencies of Fenton process [3840]. The surface response plots of COD and color removal efficiencies of Fenton reaction regarding interaction effects between dye concentration and other selected independent variables are presented in supplementary data in Figures A.3–A.6.

The interaction for assessing decolorization efficiency implies that Dye/Fe+2 and H2O2/Fe+2 have proportional effect on color removal efficiency. It can be seen from the plots that there is an increase in decolorization efficiency with increase in control variables and dye concentration. The reason for this increasing trend is that HO radically produced is reactive enough for the cleavage of the Azo bonds (–N=N–) and at high values of H2O2/Fe+2 and Dye/Fe+2, higher concentrations of radical are available to react with the high concentration of dye [33]. In this study, 95% decolorization of acid blue dye was observed within 10 minutes of reaction time under optimum conditions. And over 90% of the color removal was observed for most of the treated samples. However, from data available in Table 2 and contour plots (Figure A.3), it was found that pH exhibits inverse relationship with dye decolorization efficiency. This is because at higher pH values, decomposition of H2O2 to oxygen and conversion of ferrous ions to ferric hydroxocomplex take place which will in turn make radical unavailable for the reaction to take place [41].

From aforementioned discussion, it is confirmed that Fenton’s reagents are active for decolorization. However, efficient degradation of dye is difficult to achieve and it requires optimization of the Fenton oxidation. Unlike decolorization, concentration of Fenton’s reagents is not directly linked with COD removal efficiency. In order to make things understandable and comparable, a detailed analysis of the effect of Fenton’s reagent on COD removal efficiency is presented in Table 4. Moreover, three- and two-dimensional contour plots obtained by CCD model are also provided in supplementary data (Figures A.4–A.6).

COD (%)
Dye : Fe+2
H2O2 : Fe+2
pHCOD (%)Reason

Effect of Dye : Fe+2 ratio (Figure A.4)10010–2315 365–72.5 (increased)At low dye concentrations, increase in Dye : Fe+2 (wt/wt) results in a decrease in Fe+2 concentration and increase in H2O2 addition (H2O2 : Fe+2) which increases the production of HO radical for dye degradation
35–5015 372–65 (decreased) Scavenging of HO radical by H2O2 at low concentrations of Fe+2 and high concentrations of H2O2 [3]
250–30010–2515 380 (increased)Optimum amounts of H2O2 and Fe+2 result in the production of HO radical adequate enough for maximum dye degradation

Effect of H2O2 : Fe+2 ratio (Figure A.5)100305–10370–72
100–150305373–68 (decreased) Less availability of HO
1003020–25373–68 (decreased) Scavenging of HO radical by excess amount of H2O2
3003010–25374–80 (increased)COD removal efficiency increases from 74% to 80% because of the proportionate amount of HO radical production

pH100–30010–505–2538080% COD removal efficiency because of the availability of Fe+2 and H2O2 in aqueous medium (optimum conditions of other variables)

(Figure A.6) 100–30010–505–25960Decomposition of H2O2 to H2O and O2 at pH above 4 [3]
Deactivation of ferrous catalyst with the formation of ferric hydroxocomplex [35]

Optimization study of the experimental results was performed by keeping all responses within desired ranges by using response surface methodology. In present studies, dye concentration was targeted to the maximum and other variables were kept in range. Based on the suggested values as given in Table 5, experiment was conducted to validate the optimized results. According to the results obtained approximately, 79.1% of the COD removal was achieved with 98.401% decolorization efficiency. This indicates a good agreement of the experimental and predicted results under optimized conditions. Comparison of the results with previously conducted studies indicates that efficiency of the process is improved in terms of both chemical consumptions and COD removal efficiency. For example, Meriç et al. [42] obtained 71% COD and 99% color removal at optimum conditions of 100 mg/L of FeSO4 and 400 mg/L of H2O2 for 100 mg/L of Reactive Black 5. In another study, 2000 mg/L of H2O2 and 200 mg/L of FeSO4 were used for 88% COD removal of 200 mg/L of Reactive Black 5 [43]. It implies that CCD under RSM is a suitable tool for optimizing Fenton treatment of the recalcitrant wastewater with improved efficiency and less consumption of chemicals.

Dye (mg/L)Dye : Fe+2 (wt/wt)H2O2 : Fe+2 (wt/wt)pH Predicted responses
COD (%)Decolorization (%)Desirability


5. Taguchi Method

5.1. Experimental Design

Taguchi method was used to determine the optimum conditions for Fenton oxidation and to make comparison with the results obtained with CCD. Minitab 16 was used for the orthogonal experimental design. Same as before, , Dye/Fe+2 (wt/wt), H2O2/Fe+2 (wt/wt), and pH were chosen as control factors for optimization through Taguchi orthogonal arrays experimental design. Each factor was varied at three levels (Table 6) and orthogonal array was selected to determine the optimal conditions with minimum number of experiments. The number of experiments required was reduced to 9 which was considerably lesser than that of CCD. It implies that only 9 experiments with different combination of parameters are required to study Fenton oxidation, which in conventional full factorial design would be experimental runs.

ParametersLevel 1Level 2Level 3

Dye : Fe+2103050
H2O2 : Fe+251525

In second stage, experiments were performed and response values were obtained. For result analysis, Taguchi method follows entirely different steps in contrast to CCD. As followed by previous studies, response values were converted into ratio which was used to analyze the results [18, 19, 4446]. For Fenton oxidation, maximum COD and color removal percentages are desired that is why “larger is better” ratio formula, as given by (3), was used to determine the value for each response. The experimental results and ratio for each run are given in Table 7.

Dye (mg/L)Dye : Fe+2 (wt/wt)H2O2 : Fe+2 (wt/wt)pH Actual values/ ratio
COD (%)Decolorization (%)CODDecolorization


In comparison to CCD, ranking of operating parameters in terms of contribution to response values is possible with Taguchi method. It also suggests the use of Taguchi method for screening the input variables during the initial stages of process investigation. In the current study, four parameters were selected and from Table 8, it is clear that for both types of responses pH was the most contributing factor. This observation was also confirmed from the literature that Fenton oxidation shows maximum efficiency at pH 3 and deviation from this value decreases the efficiency of the process [41].

LevelDyeDye : Fe+2 (wt/wt)H2O2 : Fe+2 (wt/wt)pH

COD removal


5.2. Statistical Analysis

In CCD, analysis of variance (ANOVA) and regression models are used as evaluation criteria for the statistical analysis of the results. However, Taguchi method uses signal to noise ratio () as a main approach for the analysis. Nevertheless, ANOVA can be employed for the evaluation of experimental results with a main objective of determining the contribution of each factor to variance of result. It is similar to regression analysis which is used to study and model the relationship between response and one or more independent variables [47].

Table 9 shows ANOVA results for both COD and color removal efficiency of acid blue 113 dye. From the table, it is clear that pH exhibits maximum contribution with percent contribution of 66.2% and 87.62% for both color and COD removal, respectively. It is also in agreement with the values reported in Table 8. Moreover, Sohrabi et al. [19] also reported pH to be the most significant factor (percent contribution: 84.90%) while demonstrating Fenton and photo-Fenton process. This observation can be supported by the fact that, at high pH values, decomposition of H2O2 to oxygen and formation of hydroxocomplex take place which results in low efficiency of Fenton oxidation even at optimum conditions of Fenton’s reagent [41]. Thus, with confidence level of 95%, it has been confirmed that Taguchi method with minimum number of experiments can be used as an alternative to CCD for Fenton oxidation. In order to confirm it, optimized values for the responses were computed to countercheck it with that of CCD.

FactorsDegree of freedomSum of squaresMean Square ratio valuePercent contribution (%)

COD (%)
 Dye : Fe+221671670.150.8604.89
 H2O2 : Fe+222112110.200.8266.17%
Decolorization (%)
 Dye : Fe+2218.
 H2O2 : Fe+

5.3. Confirmation Test and Optimization of Operating Parameters by Taguchi Method

In Taguchi’s method, confirmation test is important to verify the experimental results. However, this method has limitations for multiple response processes like Fenton oxidation. This is because unlike CCD it suggests individual sets of optimized parameters for individual responses. In the present case study, decolorization and COD removal efficiencies were selected as responses for determining the optimized set of the input variables. Figure 5 shows the best conditions for the Fenton oxidation. The highest ratios corresponding to the COD and color removal suggest the best levels for each of the parameters. It can be seen from the figure that pH shows the highest ratio for both COD and color removal efficiencies. Thus, it confirms the maximum contribution of pH to COD and color removal efficiencies.

Figure 5 shows the effect of operating parameters on COD and color removal efficiency in terms of ratio. The highest ratios corresponding to the COD and color removal is desirable and suggests the best levels for each of the parameters. In the present example, Dye (), Dye : Fe+2 (), H2O2 : Fe+2 (), and pH () were optimized conditions for COD removal efficiency while Dye (), Dye : Fe+2 (), H2O2 : Fe+2 (), and pH () were the optimized set for decolorization. Optimized values corresponding to these levels are given in Table 10.

DyeDye : Fe+2H2O2 : Fe+2pHPredicted Actual

COD removal 3001025379.06 81.2
Decolorization 3001525399.3199.04

Evaluation of operating parameters by using interaction plots indicates that dye concentration of 200 mg/L in different combinations with other parameters shows maximum COD and decolorization efficiency. The interaction plots were produced in Minitab and are provided in supplementary data in Figures A.7 and A.8. From the figures, it can be seen that pH at shows maximum efficiency for all combinations of operating parameters. However, dye concentration at in combination with Dye : Fe+2 of 30 showed 80% COD removal efficiency. Same results were obtained when H2O2 : Fe+2 ratio was set at 25. Similar trend was observed while studying the interactive effects of operating parameters for determining decolorization efficiency and maximum 99% decolorization efficiency was obtained. Comparison of the findings of the interactive plots of both types of statistical techniques showed that Taguchi method with simple graphical presentation can well explain the interaction of operating parameters. The results obtained were in good agreement with each other. Thus, Taguchi method can be used as a substitute for CCD for assessing the experimental results.

Predicted values of the responses for these optimized values can be computed by using (5). This equation considers only significant parameters. Although only pH is significant parameter according to ANOVA analysis, Dye, Dye : Fe+2, and H2O2 : Fe+2 have to be optimized to maximize the COD removal and decolorization efficiency. The predicted ratios for optimized conditions were computed and obtained values were 37.96 and 39.94 for COD and color removal, respectively. The predicted and actual values obtained as a result of confirmation test are also given in Table 9. The predicted values obtained through Taguchi method were almost close to those obtained through CCD except H2O2 : Fe+2 ratios. The H2O2 : Fe+2 ratio is higher when compared to that obtained with CCD; the main reason for this is that with design expert, it is possible to select the desired range for experimental parameters. With CCD, H2O2 : Fe+2 ratio was selected as the lowest one. However, experimental results obtained under optimized conditions (Taguchi method) are close to predicted values as well as those obtained in CCD. This shows that for optimizing operating parameters for Fenton oxidation, Taguchi method is suitable and economical approach and can be used as a substitute to central composite design.

6. Conclusion

This study was planned to compare two optimization techniques such as central composite design and Taguchi orthogonal array and Fenton oxidation with four parameters; that is, , Dye : Fe+2, H2O2 : Fe+2, and pH were selected as a case study. From this study, the following points can be drawn as conclusion.(i)Taguchi method offered nine experiments for analyzing the Fenton oxidation process while CCD suggested 30 experiments.(ii)At optimized conditions, Fenton process would be able to achieve 81% and 79% COD removal efficiency and 99% decolorization efficiency for both CCD and Taguchi method, respectively.(iii) ratios and ANOVA under Taguchi method showed pH to be the most contributed factor with percent contribution of 87.62% and 66.2% for COD and decolorization efficiency, respectively.(iv)Interactive study of operating parameters by 3D contour plots produced by CCD also exhibits pH and H2O2 : Fe+2 the most significant factors. However, quantification of contribution is not possible with CCD.

Thus, it can be concluded that Taguchi method is a robust statistical tool for experimental design and process optimization. Data analysis and optimization of operating parameters are possible with fewest numbers of experiments, less computational experience, and graphs obtained which are easy to read and understand. The optimized values obtained for both cases were in good agreement with each other which shows the potential of Taguchi method to be used in chemical engineering applications. Therefore, it can be concluded that it can be used as a substitute for central composite design for Fenton oxidation process.


AOPs:Advance oxidation process
Fe+2:Ferrous iron
Fe+3:Ferric iron
H2O2:Hydrogen peroxide
HO:Hydroxyl radical
OFAT:One factor at time
DOE:Design of experiments
CCD:Central composite design
RSM:Response surface methodology
ANN:Artificial neural network
ANOVA:Analysis of variance
S/N:Signal to noise ratio.

Conflict of Interests

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


The authors are grateful to the Faculty of Engineering, University of Malaya, University of Malaya High Impact Research Grant (HIR-MOHE-D000037-16001) from the Ministry of Higher Education Malaysia, and University of Malaya Bright Sparks Unit which financially supported this work.

Supplementary Materials

The supplementary data for supporting the information provided in Table 4, Section 4.1. and interaction of parameters for Taguchi method has been provided in Supplementary Material

  1. Supplementary Material


  1. R. Andreozzi, V. Caprio, A. Insola, and R. Marotta, “Advanced oxidation processes (AOP) for water purification and recovery,” Catalysis Today, vol. 53, no. 1, pp. 51–59, 1999. View at: Publisher Site | Google Scholar
  2. L. Gomathi Devi, S. Girish Kumar, K. Mohan Reddy, and C. Munikrishnappa, “Photo degradation of Methyl Orange an azo dye by Advanced Fenton Process using zero valent metallic iron: Influence of various reaction parameters and its degradation mechanism,” Journal of Hazardous Materials, vol. 164, no. 2-3, pp. 459–467, 2009. View at: Publisher Site | Google Scholar
  3. K. Ntampegliotis, A. Riga, V. Karayannis, V. Bontozoglou, and G. Papapolymerou, “Decolorization kinetics of Procion H-exl dyes from textile dyeing using Fenton-like reactions,” Journal of Hazardous Materials, vol. 136, no. 1, pp. 75–84, 2006. View at: Publisher Site | Google Scholar
  4. B. H. Hameed and T. W. Lee, “Degradation of malachite green in aqueous solution by Fenton process,” Journal of Hazardous Materials, vol. 164, no. 2-3, pp. 468–472, 2009. View at: Publisher Site | Google Scholar
  5. X. Zhu and B. E. Logan, “Using single-chamber microbial fuel cells as renewable power sources of electro-Fenton reactors for organic pollutant treatment,” Journal of Hazardous Materials, vol. 252-253, pp. 198–203, 2013. View at: Publisher Site | Google Scholar
  6. B. H. Diya'uddeen, A. R. Abdul Aziz, and W. M. A. W. Daud, “On the limitation of Fenton oxidation operational parameters: a review,” International Journal of Chemical Reactor Engineering, vol. 10, no. 1, article R2, 2012. View at: Publisher Site | Google Scholar
  7. A. Ellert and A. Grebe, “Process optimization made easy: design of experiments with multi-bioreactor system BIOSTAT Qplus,” Nature Methods, vol. 8, no. 4, 2011. View at: Google Scholar
  8. P. W. Araujo and R. G. Brereton, “Experimental design II. Optimization,” TrAC Trends in Analytical Chemistry, vol. 15, no. 2, pp. 63–70, 1996. View at: Publisher Site | Google Scholar
  9. R. Leardi, “Experimental design in chemistry: a tutorial,” Analytica Chimica Acta, vol. 652, no. 1-2, pp. 161–172, 2009. View at: Publisher Site | Google Scholar
  10. M. B. Kasiri, H. Aleboyeh, and A. Aleboyeh, “Modeling and optimization of heterogeneous photo-fenton process with response surface methodology and artificial neural networks,” Environmental Science and Technology, vol. 42, no. 21, pp. 7970–7975, 2008. View at: Publisher Site | Google Scholar
  11. M. Azami, M. Bahram, and S. Nouri, “Central composite design for the optimization of removal of the azo dye, Methyl Red, from waste water using Fenton reaction,” Current Chemistry Letters, vol. 2, no. 2, pp. 57–68, 2013. View at: Google Scholar
  12. D. B. Hasan, A. R. A. Aziz, and W. M. A. Wan Daud, “Using D-optimal experimental design to optimise remazol black B mineralisation by Fenton-like peroxidation,” Environmental Technology, vol. 33, no. 10, pp. 1111–1121, 2012. View at: Publisher Site | Google Scholar
  13. E. C. Catalkaya and F. Kargi, “Response surface analysis of photo-fenton oxidation of simazine,” Water Environment Research, vol. 81, no. 7, pp. 735–742, 2009. View at: Publisher Site | Google Scholar
  14. T. J. S. Anand, “DOE based statistical approaches in modeling of laser processing—review and suggestion,” International Journal of Engineering & Technology IJET-IJENS, vol. 10, no. 4, pp. 1–7, 2010. View at: Google Scholar
  15. L. A. Lu, Y. S. Ma, A. Daverey, and J. G. Lin, “Optimization of photo-Fenton process parameters on carbofuran degradation using central composite design,” Journal of Environmental Science and Health - Part B Pesticides, Food Contaminants, and Agricultural Wastes, vol. 47, no. 6, pp. 553–561, 2012. View at: Publisher Site | Google Scholar
  16. F. Torrades, S. Saiz, and J. A. García-Hortal, “Using central composite experimental design to optimize the degradation of black liquor by Fenton reagent,” Desalination, vol. 268, no. 1–3, pp. 97–102, 2011. View at: Publisher Site | Google Scholar
  17. A. Zuorro, M. Fidaleo, and R. Lavecchia, “Response surface methodology (RSM) analysis of photodegradation of sulfonated diazo dye Reactive Green 19 by UV/H2O2 process,” Journal of Environmental Management, vol. 127, pp. 28–35, 2013. View at: Publisher Site | Google Scholar
  18. T. C. Cheng, K. S. Yao, Y. H. Hsieh, L. L. Hsieh, and C. Y. Chang, “Optimizing preparation of the TiO2 thin film reactor using the Taguchi method,” Materials and Design, vol. 31, no. 4, pp. 1749–1751, 2010. View at: Publisher Site | Google Scholar
  19. M. R. Sohrabi, A. Kavaran, S. Shariati, and S. Shariati, “Removal of Carmoisine edible dye by Fenton and photo Fenton processes using Taguchi orthogonal array design,” Arabian Journal of Chemistry, 2014. View at: Publisher Site | Google Scholar
  20. K. C. Namkung, A. E. Burgess, D. H. Bremner, and H. Staines, “Advanced Fenton processing of aqueous phenol solutions: a continuous system study including sonication effects,” Ultrasonics Sonochemistry, vol. 15, no. 3, pp. 171–176, 2008. View at: Publisher Site | Google Scholar
  21. M. A. Bezerra, R. E. Santelli, E. P. Oliveira, L. S. Villar, and L. A. Escaleira, “Response surface methodology (RSM) as a tool for optimization in analytical chemistry,” Talanta, vol. 76, no. 5, pp. 965–977, 2008. View at: Publisher Site | Google Scholar
  22. L. Ye, M. Ying, L. Xu, C. Guo, L. Li, and D. Wang, “Optimization of inductive angle sensor using response surface methodology and finite element method,” Measurement, vol. 48, pp. 252–262, 2014. View at: Google Scholar
  23. T. Ölmez, “The optimization of Cr(VI) reduction and removal by electrocoagulation using response surface methodology,” Journal of Hazardous Materials, vol. 162, no. 2-3, pp. 1371–1378, 2009. View at: Publisher Site | Google Scholar
  24. S. Baroutian, M. K. Aroua, A. A. A. Raman, and N. M. N. Sulaiman, “A packed bed membrane reactor for production of biodiesel using activated carbon supported catalyst,” Bioresource Technology, vol. 102, no. 2, pp. 1095–1102, 2011. View at: Publisher Site | Google Scholar
  25. D. Montgomery, Design and Analysis of Experiments, John Wiley & Sons, New York, NY, USA, 2001.
  26. A. El-Ghenymy, S. Garcia-Segura, R. M. Rodríguez, E. Brillas, M. S. El Begrani, and B. A. Abdelouahid, “Optimization of the electro-Fenton and solar photoelectro-Fenton treatments of sulfanilic acid solutions using a pre-pilot flow plant by response surface methodology,” Journal of Hazardous Materials, vol. 221-222, pp. 288–297, 2012. View at: Publisher Site | Google Scholar
  27. N. Ali, V. F. Neto, S. Mei et al., “Optimisation of the new time-modulated CVD process using the Taguchi method,” Thin Solid Films, vol. 469-470, pp. 154–160, 2004. View at: Publisher Site | Google Scholar
  28. S. Maghsoodloo, G. Ozdemir, V. Jordan, and C. Huang, “Strengths and limitations of taguchi's contributions to quality, manufacturing, and process engineering,” Journal of Manufacturing Systems, vol. 23, no. 2, pp. 73–126, 2004. View at: Publisher Site | Google Scholar
  29. G. Barman, A. Kumar, and P. Khare, “Removal of congo red by carbonized low-cost adsorbents: process parameter optimization using a Taguchi experimental design,” Journal of Chemical and Engineering Data, vol. 56, no. 11, pp. 4102–4108, 2011. View at: Publisher Site | Google Scholar
  30. S. K. Gauri and S. Chakraborty, “Multi-response optimisation of WEDM process using principal component analysis,” The International Journal of Advanced Manufacturing Technology, vol. 41, no. 7-8, pp. 741–748, 2009. View at: Publisher Site | Google Scholar
  31. P. J. Ross, Taguchi Techniques for Quality Engineering, McGraw-Hill, New York, NY, USA, 1998.
  32. M. Kaladhar, K. V. Subbaiah, C. Srinivasa Rao, and K. Narayana Rao, “Application of Taguchi approach and utility concept in solving the multi-objective problem when turning AISI 202 austenitic stainless steel,” Journal of Engineering Science and Technology Review, vol. 4, no. 1, pp. 55–61, 2011. View at: Google Scholar
  33. C. L. Hsueh, Y. H. Huang, C. C. Wang, and C. Y. Chen, “Degradation of azo dyes using low iron concentration of Fenton and Fenton-like system,” Chemosphere, vol. 58, no. 10, pp. 1409–1414, 2005. View at: Publisher Site | Google Scholar
  34. S. Wang, “A Comparative study of Fenton and Fenton-like reaction kinetics in decolourisation of wastewater,” Dyes and Pigments, vol. 76, no. 3, pp. 714–720, 2008. View at: Publisher Site | Google Scholar
  35. M. Neamtu, A. Yediler, I. Siminiceanu, and A. Kettrup, “Oxidation of commercial reactive azo dye aqueous solutions by the photo-Fenton and Fenton-like processes,” Journal of Photochemistry and Photobiology A: Chemistry, vol. 161, no. 1, pp. 87–93, 2003. View at: Publisher Site | Google Scholar
  36. M. Y. Can, Y. Kaya, and O. F. Algur, “Response surface optimization of the removal of nickel from aqueous solution by cone biomass of Pinus sylvestris,” Bioresource Technology, vol. 97, no. 14, pp. 1761–1765, 2006. View at: Publisher Site | Google Scholar
  37. B. K. Körbahti, “Response surface optimization of electrochemical treatment of textile dye wastewater,” Journal of Hazardous Materials, vol. 145, no. 1-2, pp. 227–286, 2007. View at: Publisher Site | Google Scholar
  38. B. Bianco, I. de Michelis, and F. Vegliò, “Fenton treatment of complex industrial wastewater: optimization of process conditions by surface response method,” Journal of Hazardous Materials, vol. 186, no. 2-3, pp. 1733–1738, 2011. View at: Publisher Site | Google Scholar
  39. C. T. Benatti, C. R. G. Tavares, and T. A. Guedes, “Optimization of Fenton's oxidation of chemical laboratory wastewaters using the response surface methodology,” Journal of Environmental Management, vol. 80, no. 1, pp. 66–74, 2006. View at: Publisher Site | Google Scholar
  40. M. Ahmadi, F. Vahabzadeh, B. Bonakdarpour, E. Mofarrah, and M. Mehranian, “Application of the central composite design and response surface methodology to the advanced treatment of olive oil processing wastewater using Fenton's peroxidation,” Journal of Hazardous Materials, vol. 123, no. 1–3, pp. 187–195, 2005. View at: Publisher Site | Google Scholar
  41. Y. W. Kang and K. Hwang, “Effects of reaction conditions on the oxidation efficiency in the Fenton process,” Water Research, vol. 34, no. 10, pp. 2786–2790, 2000. View at: Publisher Site | Google Scholar
  42. S. Meriç, D. Kaptan, and T. Ölmez, “Color and COD removal from wastewater containing Reactive Black 5 using Fenton's oxidation process,” Chemosphere, vol. 54, no. 3, pp. 435–441, 2004. View at: Publisher Site | Google Scholar
  43. M. E. Argun and M. Karatas, “Application of Fenton process for decolorization of reactive black 5 from synthetic wastewater: kinetics and thermodynamics,” Environmental Progress & Sustainable Energy, vol. 30, no. 4, pp. 540–548, 2011. View at: Publisher Site | Google Scholar
  44. A. Chowdhury, R. Chakraborty, D. Mitra, and D. Biswas, “Optimization of the production parameters of octyl ester biolubricant using Taguchi's design method and physico-chemical characterization of the product,” Industrial Crops and Products, vol. 52, pp. 783–789, 2014. View at: Google Scholar
  45. Ö. Gökkuş, Y. Ş. Yildiz, and B. Yavuz, “Optimization of chemical coagulation of real textile wastewater using Taguchi experimental design method,” Desalination and Water Treatment, vol. 49, no. 1–3, pp. 263–271, 2012. View at: Publisher Site | Google Scholar
  46. O. Gökkuş and Y. Ş. Yıldız, “Investigation of the effect of process parameters on the coagulationflocculation treatment of textile wastewater using the taguchi experimental method,” Fresenius Environmental Bulletin, vol. 23, no. 2, pp. 1–8, 2014. View at: Google Scholar
  47. L. Hsieh, H. Kang, H. Shyu, and C. Chang, “Optimal degradation of dye wastewater by ultrasound/Fenton method in the presence of nanoscale iron,” Water Science and Technology, vol. 60, no. 5, pp. 1295–1301, 2009. View at: Publisher Site | Google Scholar

Copyright © 2014 Anam Asghar et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.