Evaluation of Reaction Parameters of the Phenol Steam Reforming over Ni / Co on ZrO 2 Using the Full Factorial Experimental Design

Full factorial experimental design with 32 runs was used to investigate the significant and interaction variable of the reaction parameters on phenol steam reforming toward hydrogen production. Effects of selected factors on the phenol conversion (Y1) and hydrogen yield (Y2) were evaluated. These factors were as follows: (A) temperature (500 and 800  ̋C); (B) feed flow rate (0.16–0.46 mL/min); (C) catalyst weight (0.1–0.3 g); (D) Ni-Co ratio (0–1); and (E) phenol concentration in the feed (2–10 wt %). Ni and Co over ZrO2 support for catalytic performance of phenol steam reforming (SRP) was prepared by the impregnation method. The result indicated that all the main independent variables had significant influence on the dependent variable of Y1 and Y2 with a range of 2.7%–96.8% and 21.4%–72.4%, respectively. Additionally, some interaction variables like AE, BE, CE, and DE have also influenced the Y1 and Y2 responses. This design showed that the best initial conditions that produced maximum Y1 and Y2 responses were at 800  ̋C, 0.16 mL/min feed flow rate, 0.3 g of catalyst, 0 ratio of Ni-Co (Co/ZrO2), and 10 wt % of phenol in the feed, where the phenol conversion was predicted to be 94.98% and the hydrogen yield was predicted to be 67.4%. Within the limits the variables were examined, a regression model which well-fitted the experimental data was proposed. The regression model were reduced to simplify and to get the significant regression coefficient with p-value less than 0.05.


Introduction
Due to fossil fuels exhaustion, the growth of concern about global warming, and the increase in major air toxics, hydrogen has emerged as one of the best clean energy vectors [1].H 2 can be produced from maintainable liquid fuels or phenols that consist of around 38 wt % [2] of the unwanted component of pyrolysis oil.However, this component can be transformed into H 2 rich gas, which makes phenol an ideal potential H 2 source.Based on previous research, there are two successive reactions that may take place in the steam reforming of phenol [3].One of the reactions for H 2 production from phenol is steam reforming of phenol (SRP), which can produce more H 2 yield based on stoichiometry in Equations ( 1) and ( 2 Ni-based catalysts with a high activity for the phenol steam reforming process have been investigated extensively in recent years [4][5][6] due to the ability to cleave C-C, O-H, and C-H bonds [7].In addition, cobalt can assist with C-C bond cleaving at temperatures as low as 400 ˝C; this shows high production of H 2 and CO 2 [8] due to the highly favourable water gas shift reaction [9][10][11].However, there are a few items that can be identified as the catalyst in the reforming of phenol: the conversion of phenol, and the yield of H 2 , CO or CO 2 .Many issues affect this performance, such as the steam to phenol (S/P) ratio, catalyst composition, flow rate of carrier gas, operating temperature, reactant feed rate and space time, preparation method, catalyst particle size, amount of catalyst, and ratio of active metals in multiple active metal catalyst use.Exothermic reactions, high phenol conversion, and the rate of gas production all occur at a high temperature.The same results were found in previous works regarding the temperature effect on phenol conversion.Some researchers have found that the conversion of phenol increased with the increasing of temperature [6,[12][13][14].Based on the works of Rioche et al. [13] and Shurong et al. [15], the increase in catalyst activity is based on the increasing temperature.They have found that a complete phenol conversion occurs at high temperatures.There has been no research studying the amount of catalyst, the feed reactant flow, or the metal ratio in the steam reforming of phenol.However, some previous works did use other feeds like methanol, acetic acid, glycerol, and alcohol.
Regarding previous works, increasing the space time [16,17] or expanding the contact time by increasing the amount of the catalyst [18,19] can result in an increase in the activity of the catalyst.However, the catalytic performance also depends on the kind of catalyst [17,20,21] and the technique of catalyst preparation [17].
DOE methods are more suitable than one-factor-at-a-time techniques because by studying multiple factors at the same time, a high efficiency can be obtained.These techniques have the following common benefits [22]: (1) more information per experiment than accidental methods; (2) a decrease in the number and fee of experiments; (3) making the calculation of the interaction amongst possible variables within the series studied, leading to better facts about the process; and (4) simplifying the description of the significant operation conditions for the scale-up process.Therefore, these techniques have been considered as powerful tools for process examination [23][24][25][26][27].
A numerical design of experiment (DOE) was applied in order to investigate the interaction amongst factors and the significance of factors on phenol conversion.The DOE capabilities improve the processes by screening the factors to determine which are important for explaining the process variation.After screening the factors, DOE helped to understand how these factors interact during the process.It can then find the factor settings that produced optimal process performance.Thus the DOE method is a great tool for process performance testing [28][29][30].
The aim of the current study is to investigate the relationship between the five control parameters, significant variables, and the interaction between variables via full factorial design.Therefore, a detailed study of DOE in the catalytic performance of steam reforming of phenol (SRP) for hydrogen production was carried out.The performance of the catalyst was stated according to the phenol conversion and hydrogen yield.

Catalyst Preparation
The impregnation method was applied to produce the Ni or Co over ZrO 2 catalyst.The advantage of the impregnation method is that is produces a high concentration of active materials on the catalyst surface [31].The procedure of this method was reported by Athanasios et al. [32].The techniques used to make the catalyst were as follows: firstly, 9 g of Zirconium oxide (ZrO 2 ) was dispersed in 250 mL of deionized water and stirred for 2 h at 90 ˝C.While stirring the slurry, 4.94 g of active metal (nickel nitrate hexahydrate (Ni(NO 3 ) 2 ¨6H 2 O) or cobalt nitrate hexahydrate (Co(NO 3 ) 2 ¨6H 2 O) was added to the slurry.The mixture was held until it became an extremely viscous paste.The paste was dried at 110 ˝C for 24 h and the catalyst was calcined at 750 ˝C for about three hours.Then, the catalyst was pelletized and sieved on two layers of 35 and 34 mesh to acquire particles sizes between 1.0 mm and 1.4 mm.

Catalyst Characterization
The crystalline structure of the reduced catalyst was determined by X-ray diffraction (High-Resolution X-ray Diffractometer brand Bruker D8 Advance, BRUKER AXS INC., Fitchburg, WI, USA) using a Cu Kα radiation at 40 kV and 30 mA.Diffraction angles were measured in steps of 0.02 ˝, at 1 s/step in the range of 20 to 80 (2θ).The size of the metal crystallites was calculated from line broadening with the help of the Scherrer equation [33].The total surface area of catalysts was measured using the multipoint BET-N 2 on a surface area analyzer (Micromeritics, Gemini 2360, Norcross, GA, USA).Prior to the analysis, the samples were degassed at 200 ˝C for 2 h to remove unwanted components on the surface before outgassing the samples.
A chemisorption analyser (Micromeritics, Chemisorb 2720, Norcross, GA, USA), was used to evaluate the reducibility of the active metals (nickel and cobalt) on the support by applying a temperature-programmed reduction of hydrogen (TPR-H 2 ).Before the reduction, 25 mg of the catalyst was treated at 300 ˝C under high purity 99.99% helium with the flow rate of 20 mL/min for an hour to remove moisture and other gas impurities.The TPR-H 2 profile was obtained by ramping up the temperature at 10 ˝C/min, 20 mL/min of 10% (vol.)H 2 /Ar between 50 ˝C and 800 ˝C.
Thermogravimetric analysis (TGA) profiles were carried out in a Perkin Elmer TGA instrument (Waltham, MA, USA) operated under nitrogen flowing in the heating rate of 10 ˝C/min.Exothermic weight loss was observed at the temperature range between 30 and 950 ˝C.

Catalyst Performance Test
A 0.2 g weight of catalyst was employed to investigate the performance of the catalyst.The catalyst bed was located inside a 6-mm ID quartz tubular reactor.A coaxially centred thermocouple was incorporated into the reactor.The nitrogen stream passed through the catalyst at 300 ˝C and subsequently reduced in situ with a 30 mL/min stream of pure hydrogen for 1 hour at 600 ˝C.The mixture of phenol and water was fed into the pre-heater by means of an HPLC pump (Bio-RadTM, Series 1350, Berkeley, CA, USA) before mixing with carrier N 2 .The samples in reduced forms were tested in the steam reforming of phenol (H 2 O/phenol(mol/mol) = 9:1) at an atmospheric pressure and a temperature range of 500 to 800 ˝C, maintaining the reaction for 6 h in order to check if there was any deactivation of the catalysts.The reaction products were analysed on-line by GC with TCD (Agilent 6890N, Agilent Technologies, Santa Clara, CA, USA ) equipped with Carboxen Plot 1010 capillary column (Fused Silica, 30 m ˆ0.53 mm, Sigma-Aldrich Corporation, St. Louis, MO, USA) connected in series, using argon as the carrier gas.The 95% level of confidence interval was chosen for the uncertainty in figure values for gas production and phenol conversion.The flow diagram of this work was the same as the study done by Abdullah et al. [34] and is shown in Figure 1.As a result of the products variation, the catalyst performance was characterized by two parameters; phenol conversion percentage and product yield.They were calculated according to Equations ( 3) and ( 4

Statistical Design of Experiment
A full experiment was carried out in order to study the relationship between the control variables, significant variables, and the interactions between the test parameters.A 2 5 -factorial design was applied in the present investigation.This required a selection of appropriate responses, factors, and levels.Phenol conversion and hydrogen yield were the responses selected.The factors and the experimental levels for each factor were selected based on literature values, available resources, and results from preliminary experiments.The selection of the factors was based on the operating conditions that had a significant influence on the SRP, according to previous experiments (with different feeds) [12][13][14][16][17][18].The considered factors chosen were temperature, A, Feed flow rate, B, Catalyst amount, C, the catalyst ratio of Ni to Co over ZrO2, D, and phenol concentration in the feed, E. The factors and their levels were assigned in Table 1.A Minitab statistical software (®16.2.4 , Minitab, Inc., Pennsylvania State University, University Park, PA, USA, 2013) was used to generate the design matrix and to conduct the data analysis.Also, the interaction among the factors was determined using the same software.The design matrix of un-coded values for the factors and the response for all the experimental runs are given in Table 2. Five factors were analyzed at two levels, using a full factorial 2 5 resulting in 32 runs, in a random order to minimize the effect of uncontrolled variables on responses (Y1 and Y2).
Since one of the aims of this study was to establish a relationship between the five control variables and the responses of the system (phenol conversion and hydrogen yield) a linear regression model was fitted to predict the responses as a function of independent variables and their interactions.In general, the response for the linear regression model is described in Equation ( 5):

Statistical Design of Experiment
A full experiment was carried out in order to study the relationship between the control variables, significant variables, and the interactions between the test parameters.A 2 5 -factorial design was applied in the present investigation.This required a selection of appropriate responses, factors, and levels.Phenol conversion and hydrogen yield were the responses selected.The factors and the experimental levels for each factor were selected based on literature values, available resources, and results from preliminary experiments.The selection of the factors was based on the operating conditions that had a significant influence on the SRP, according to previous experiments (with different feeds) [12][13][14][16][17][18].The considered factors chosen were temperature, A, Feed flow rate, B, Catalyst amount, C, the catalyst ratio of Ni to Co over ZrO 2 , D, and phenol concentration in the feed, E. The factors and their levels were assigned in Table 1.A Minitab statistical software ( ® 16.2.4, Minitab, Inc., Pennsylvania State University, University Park, PA, USA, 2013) was used to generate the design matrix and to conduct the data analysis.Also, the interaction among the factors was determined using the same software.The design matrix of un-coded values for the factors and the response for all the experimental runs are given in Table 2. Five factors were analyzed at two levels, using a full factorial 2 5 resulting in 32 runs, in a random order to minimize the effect of uncontrolled variables on responses (Y 1 and Y 2 ).
Since one of the aims of this study was to establish a relationship between the five control variables and the responses of the system (phenol conversion and hydrogen yield) a linear regression model was fitted to predict the responses as a function of independent variables and their interactions.In general, the response for the linear regression model is described in Equation ( 5): where Y i is the predicted response, β 0 is the offset term, β n is the regression coefficient that corresponds to the main effects and interaction effects, and ε is the error term.The level of significance of each coefficient was considered by comparing the smaller p-value (<0.05) or greater t-value.The t-test was used to determine the significance of the regression coefficients of the parameters.The p-values were used as tools to check the significance of each variable and the interaction among the variables.As the first step, a significance test was performed on all the possible terms in models for phenol conversion (Y 1 ) and hydrogen yield (Y 2 ).The results are illustrated in Pareto charts and Normal plots.Then a reduced model was obtained by removing the insignificant terms from the first model with a new reduced regression equation.The R 2 value was used to determine the fitting quality of the observed results to the proposed model.The regression analysis, linear regression plots (observed versus predicted value), and Pareto chart of the dependent variables, main effect and interaction plots, surface and contour plots, and response optimizer plots were plotted.The XRD pattern of the Ni/ZrO 2 and Co/ZrO 2 catalysts contains Bragg peaks at around 2θ = 44.52 ˝, which corresponds to metallic Ni.This can be attributed to the presence of the cubic [JCPDS 45-1027] structure of Ni and is in agreement with previous research [6,[35][36][37].The diffraction peaks at 2θ = 44.37 correspond to metallic Co [JCPDS 01-1254], similar to reports from previous works [38,39].It can be inferred that nickel oxide (NiO) and cobalt oxide (CoO) were reduced at 600 ˝C to metallic Ni [40,41] and metallic Co [42], respectively.The relative intensity of XRD peaks of the supported catalysts depends on the nature of the metal oxide precursors formed during the calcination step [43].The XRD patterns of the ZrO 2 support resulted in a number of peaks located at 2θ of 24.47  211), (300), (320), (400), (410), (411), (331), (332), and (422) crystal planes, which can be attributed to the monoclinic ZrO 2 (m-ZrO 2 ) and tetragonal ZrO 2 (t-ZrO 2 ) [JCPDS 83-0943] structure of the ZrO 2 .This conforms with previous research [44][45][46][47].
To have better insights into the nature of the cobalt and nickel phases, TPR-H2 of the calcined catalysts was carried out.This helped to explore the level of reduction of the Co and Ni oxides' species.Figure 3 shows that only single reduction zones were observed.This corresponds to the area below 400 °C (~280-450 °C), corresponding to the reduction of NiO [48,49] and CoO [50][51][52] particles with weak metal-support interaction.Thus, the reduction of the NiO and CoO to metallic Ni and Co agrees with the XRD analysis.In particular, the reductions were peaks shifted to higher temperatures as the Co/ZrO2 catalyst.This suggests that higher temperatures are needed in order to reduce the cobalt oxides (>452 °C) completely, as compared to nickel oxides [53].Therefore, these catalysts were reduced at 600 °C.To have better insights into the nature of the cobalt and nickel phases, TPR-H 2 of the calcined catalysts was carried out.This helped to explore the level of reduction of the Co and Ni oxides' species.Figure 3 shows that only single reduction zones were observed.This corresponds to the area below 400 ˝C (~280-450 ˝C), corresponding to the reduction of NiO [48,49] and CoO [50][51][52] particles with weak metal-support interaction.Thus, the reduction of the NiO and CoO to metallic Ni and Co agrees with the XRD analysis.In particular, the reductions were peaks shifted to higher temperatures as the Co/ZrO 2 catalyst.This suggests that higher temperatures are needed in order to reduce the cobalt oxides (>452 ˝C) completely, as compared to nickel oxides [53].Therefore, these catalysts were reduced at 600 ˝C.

Factors Screening in a Full 2 5 Factorial Design for the Initial Regression Model
Table 2 gives the experimental values of the two main responses in each experiment.The 2-level factorial was employed to improve the correlation between the five process variables, temperature, feed flow rate, amount of catalyst and the catalyst of Ni to Co ratio, and phenol concentration in the feed (coded as A, B, C, D, and E, respectively) to responses of phenol conversion (Y1) and hydrogen yield (Y2). Thirty-two sets of experiments were performed to explore the effects of the reaction parameters on the phenol conversion and the generation of hydrogen in the phenol steam reforming process.It was found that the conversion of the phenol varied from 2.7% to 96.8% and H2 yield was in the range of 21.4%-72.4%.The maximum catalyst activity for Y1 of 96.8 and Y2 of 72.4% was found to be at 800 °C, 0.16 mL/min feed flow rate, 0.3 g of catalyst, 0 ratio of Ni-Co (Co/ZrO2), and 10 wt % of phenol in the feed (run number 32).
To identify whether the results are "real" or "chance," a Normal plot (Figure 4) was used to display the estimates of the main effects among the factors as well as their interaction.To obtain this Normal plot, effects were ranked in ascending order and plotted against their corresponding cumulative probability [54].All the effects that tended to fall along a straight line in the plot (distributed with mean zero and marked in black) do not demonstrate any significant effect on the responses.In comparison, the significant effects were located far from the straight line and had non-zero means.They appear in red.The significant effects get more distanced from the straight line as they become larger.From Figure 4, it can be seen that all of the main parameters (A, B, C, D, and E) and their interactions (AE, BE, CE, and DE) are found to be "real" parameters for Y1 and Y2.The red marks, which lie to the left of the fitted line, had a negative effect on the responses [55].

Factors Screening in a Full 2 5 Factorial Design for the Initial Regression Model
Table 2 gives the experimental values of the two main responses in each experiment.The 2-level factorial was employed to improve the correlation between the five process variables, temperature, feed flow rate, amount of catalyst and the catalyst of Ni to Co ratio, and phenol concentration in the feed (coded as A, B, C, D, and E, respectively) to responses of phenol conversion (Y 1 ) and hydrogen yield (Y 2 ).Thirty-two sets of experiments were performed to explore the effects of the reaction parameters on the phenol conversion and the generation of hydrogen in the phenol steam reforming process.It was found that the conversion of the phenol varied from 2.7% to 96.8% and H 2 yield was in the range of 21.4%-72.4%.The maximum catalyst activity for Y 1 of 96.8 and Y 2 of 72.4% was found to be at 800 ˝C, 0.16 mL/min feed flow rate, 0.3 g of catalyst, 0 ratio of Ni-Co (Co/ZrO 2 ), and 10 wt % of phenol in the feed (run number 32).
To identify whether the results are "real" or "chance," a Normal plot (Figure 4) was used to display the estimates of the main effects among the factors as well as their interaction.To obtain this Normal plot, effects were ranked in ascending order and plotted against their corresponding cumulative probability [54].All the effects that tended to fall along a straight line in the plot (distributed with mean zero and marked in black) do not demonstrate any significant effect on the responses.In comparison, the significant effects were located far from the straight line and had non-zero means.They appear in red.The significant effects get more distanced from the straight line as they become larger.From Figure 4, it can be seen that all of the main parameters (A, B, C, D, and E) and their interactions (AE, BE, CE, and DE) are found to be "real" parameters for Y 1 and Y 2 .The red marks, which lie to the left of the fitted line, had a negative effect on the responses [55].
(distributed with mean zero and marked in black) do not demonstrate any significant effect on the responses.In comparison, the significant effects were located far from the straight line and had non-zero means.They appear in red.The significant effects get more distanced from the straight line as they become larger.From Figure 4, it can be seen that all of the main parameters (A, B, C, D, and E) and their interactions (AE, BE, CE, and DE) are found to be "real" parameters for Y1 and Y2.The red marks, which lie to the left of the fitted line, had a negative effect on the responses [55].The evaluation of the effect of each variable on the Y 1 and Y 2 may also be verified from the Pareto chart.A Pareto chart (Figure 5) visually represents the absolute values of the effects of main factors and the effects of n-way interaction of factors.The chart includes a vertical reference line at the critical t-value for an alpha of 0.05 (magnitude for a 95% confidence level) to indicate that the factors that extend past this line are potentially important.For a 95% confidence level and 32 degrees of freedom (DF), t-value is equal to 4.89 for Y 1 and 3.04 for Y 2 .It can again be seen from the Pareto chart that all effects are statistically significant.The significance of the n-way interactions could also be compared from the Pareto charts.These results were in agreement with those obtained from the analysis of the Normal plot.These plots clearly indicate that the significance of the variables and their n-way interaction towards Y 1 and Y 2 decreased in this order: The evaluation of the effect of each variable on the Y1 and Y2 may also be verified from the Pareto chart.A Pareto chart (Figure 5) visually represents the absolute values of the effects of main factors and the effects of n-way interaction of factors.The chart includes a vertical reference line at the critical t-value for an alpha of 0.05 (magnitude for a 95% confidence level) to indicate that the factors that extend past this line are potentially important.For a 95% confidence level and 32 degrees of freedom (DF), t-value is equal to 4.89 for Y1 and 3.04 for Y2.It can again be seen from the Pareto chart that all effects are statistically significant.The significance of the n-way interactions could also be compared from the Pareto charts.These results were in agreement with those obtained from the analysis of the Normal plot.These plots clearly indicate that the significance of the variables and their n-way interaction towards Y1 and Y2 decreased in this order: Based on the experimental data, regression models were fitted for Y1 and Y2, as shown in Equations ( 6) and (7), respectively.The adequacy of the initial model was tested a via parity plot for observed versus predicted values, as demonstrated in Figure 6.As seen in Figure 6, the high value of the correlation coefficient (R 2 = 1) demonstrates good correlation between the observed and the predicted responses by the initial models.The analysis of the regression coefficient for initial models Y1 and Y2 is shown in Tables 3 and 4. In such analysis, a modal term is considered significant if its p-value is less than 0.05 (bold).It was found that among the variables and their interaction, A, B, C,  D, E, AE, BE, CE  Based on the experimental data, regression models were fitted for Y 1 and Y 2 , as shown in Equations ( 6) and (7), respectively.The adequacy of the initial model was tested a via parity plot for observed versus predicted values, as demonstrated in Figure 6.As seen in Figure 6, the high value of the correlation coefficient (R 2 = 1) demonstrates good correlation between the observed and the predicted responses by the initial models.The analysis of the regression coefficient for initial models Y 1 and Y 2 is shown in Tables 3 and 4. In such analysis, a modal term is considered significant if its p-value is less than 0.05 (bold).It was found that among the variables and their interaction, A, B, C, D, E, AE      In the reduced regression model, all of the insignificant terms were removed and only the variables with a p-value of lower than 0.05 were considered.The Pareto chart of the model after reduction at p = 0.05 with 10 degrees of freedom (DF) for Y 1 and nine for Y 2 is presented in Figure 7.All the values were greater than 2.07 (p = 0.05).Moreover, the model is statistically significant as indicated by a p-value less than 0.05 (95% confidence).
The reduced linear regression model for production of hydrogen from phenol was fitted for Y 1 and Y 2 through the application of multiple regression analysis on the experimental data and subsequent removal of insignificant terms as shown in Equations ( 8) and (9).From the analysis of the regression coefficient, the value of the model helps in making inferences about whether there is a need for a complex model for better fit.The significances of the models were determined by the experimental F-value derived from the regression coefficient analysis.On the other hand, the tabulated F-values are calculated based on the degree of freedom for regression and residual at a specific confidence level.If the F test for the model is significant at the 5% level (p < 0.05), then the model is fit and can adequately account for the variation observed.The ratio of the mean square due to regression to that due to error is depicted by the F value.Generally, a good model should have its calculated F value several times greater than the tabulated.Tables 5 and 6 depict the reduced model coefficients of the factors and interactions, and the p-values of the effects in the model are shown for Y1 and Y2, respectively.According to the analysis of the regression coefficient of the reduced regression models, the models are highly significant.This was evident from the Fisher (F) test.The ratio of the mean square regression to mean square residual was 121.1 for Y1 and 60.13 for Y2 responses.Also, the probability value was very low (p-value = 0.000).
The coefficient of determination (R 2 ) and adjusted R 2 are among other criteria that are commonly used to show the best fitness of a model.Tables 5 and 6 show the calculated R 2 and adjusted R 2 values for the reduced models developed with respect to phenol conversion (Y1) and the hydrogen yield (Y2) responses, respectively.From these values, it was clear that the regression model was adequate.The ratio of the mean square due to regression to that due to error is depicted by the F value.Generally, a good model should have its calculated F value several times greater than the tabulated.Tables 5 and 6 depict the reduced model coefficients of the factors and interactions, and the p-values of the effects in the model are shown for Y 1 and Y 2 , respectively.According to the analysis of the regression coefficient of the reduced regression models, the models are highly significant.This was evident from the Fisher (F) test.The ratio of the mean square regression to mean square residual was 121.1 for Y 1 and 60.13 for Y 2 responses.Also, the probability value was very low (p-value = 0.000).
Y 1 " ´1.87 `0.01A ´3.13B `4.99C `0.24D `1.18E `0.01AE ´6.3BE `6.23CE ´1.94D (8) The coefficient of determination (R 2 ) and adjusted R 2 are among other criteria that are commonly used to show the best fitness of a model.Tables 5 and 6 show the calculated R 2 and adjusted R 2 values for the reduced models developed with respect to phenol conversion (Y 1 ) and the hydrogen yield (Y 2 ) responses, respectively.From these values, it was clear that the regression model was adequate.The model reduction has little influence on the predicted values of Y 1 and Y 2 , as indicated by the small change of the R 2 values from the analysis of the regression coefficient results (Tables 5 and 6).The validity of the regression model developed was further tested by drawing the parity charts.The parity plot helps in detecting any outliers and provides a prompt indication of the accuracy of the correlation.The observed values and predicted values of the responses are scattered close to the 95% prediction limits with 45 ˝, indicating an almost perfect fit of the developed reduced model.Figure 8 shows the parity plot of the experimental and predicted values for the reduced model.In all the cases, one or two data points out of 32 fall outside the 95% prediction limits (as shown with blue dotted lines).The model reduction has little influence on the predicted values of Y1 and Y2, as indicated by the small change of the R 2 values from the analysis of the regression coefficient results (Tables 5 and 6).The validity of the regression model developed was further tested by drawing the parity charts.The parity plot helps in detecting any outliers and provides a prompt indication of the accuracy of the correlation.The observed values and predicted values of the responses are scattered close to the 95% prediction limits with 45°, indicating an almost perfect fit of the developed reduced model.Figure 8 shows the parity plot of the experimental and predicted values for the reduced model.In all the cases, one or two data points out of 32 fall outside the 95% prediction limits (as shown with blue dotted lines).

Effects of Operating Conditions on Reaction Responses
3.4.1.Main Effect Plots for Y1 and Y2 Graphs of the main effects were plotted in order to identify the factor level providing the minimum response value, as shown in Figure 9. Figure 9 presents the mean changes that occurred in

Main Effect Plots for Y 1 and Y 2
Graphs of the main effects were plotted in order to identify the factor level providing the minimum response value, as shown in Figure 9. Figure 9 presents the mean changes that occurred in both responses (Y 1 and Y 2 ) when the levels of the main factors were changed from the lower level, passing through the central point, to a higher level.The reference line represents the overall mean.A main effect is present when the mean response changes across the levels of a factor.This plot shows that the Y 1 and Y 2 changes depend on the levels of factor E (phenol concentration) with the steepest line.The steeper the slope of the line, the greater the magnitude of the main effect.Therefore, the concentration of phenol should be studied in detail.It can be seen that for the temperature there were great changes in the mean phenol conversion (Y 1 ) and hydrogen yield (Y 2 ), characterized by a great degree of departure from the horizontal (but less than the phenol concentration main effect).Therefore, a higher temperature is desirable to get higher phenol conversion and hydrogen yield.This is because of the favorable hydrocarbon conversion at high temperatures [56].Feed flow rate and Ni/Co ratio showed slight negative changes in mean Y 1 and Y 2 .Y 1 and Y 2 decreased from 39.8% and 40.1% to 27.6% and 33.5% with the increase of feed flow rate from 0.16 to 0.46 mL/min, respectively.On the other hand, Y 1 and Y 2 also decreased from 39.5% and 40.6% to 28.02% and 32.97% with the increase of the Ni to Co ratio from 0 to 1. Nianjun et al.
[52] used Ni x Co y /Al 2 O 3 catalysts to investigate the activity of glycerin liquid reforming.They noticed that an increase in Co content results in maximum H 2 and minimum CH 4 production.They mentioned that low content of Ni is unfavorable for carbon-carbon breaking.Consequently, the performance of the reforming process was lowered.Hence, it can be concluded that temperature and phenol concentration variables highly influenced Y 1 and Y 2 , whereas feed flow rate, Ni/Co ratio, and catalyst weight had minor effects due to their lines being less sharp against the x-axis.These patterns were previously identified by statistical significance (Pareto charts and analysis of the regression coefficient data).This is because of the favorable hydrocarbon conversion at high temperatures [56].Feed flow rate and Ni/Co ratio showed slight negative changes in mean Y1 and Y2.Y1 and Y2 decreased from 39.8% and 40.1% to 27.6% and 33.5% with the increase of feed flow rate from 0.16 to 0.46 mL/min, respectively.On the other hand, Y1 and Y2 also decreased from 39.5% and 40.6% to 28.02% and 32.97% with the increase of the Ni to Co ratio from 0 to 1. Nianjun et al.
[52] used NixCoy/Al2O3 catalysts to investigate the activity of glycerin liquid reforming.They noticed that an increase in Co content results in maximum H2 and minimum CH4 production.They mentioned that low content of Ni is unfavorable for carbon-carbon breaking.Consequently, the performance of the reforming process was lowered.Hence, it can be concluded that temperature and phenol concentration variables highly influenced Y1 and Y2, whereas feed flow rate, Ni/Co ratio, and catalyst weight had minor effects due to their lines being less sharp against the x-axis.These patterns were previously identified by statistical significance (Pareto charts and analysis of the regression coefficient data).

Interaction Plots for Y1 and Y2
To visualize the impact of each factor combination and identify which factors were most influential, interaction plots were constructed.These also consider interactions among variables and can be used to optimize the operating parameters in multivariable systems [57].The mean response of the two factors is depicted by the plot for all the possible combinations of their settings.When there is interaction between the two factors, the lines appear as non-parallel [58].Figure 10 presents the full interactions effect plot matrix for five factors, temperature (A), liquid feed rate (B), catalyst weight (C), Ni to Co ratio (D), and phenol concentration (E) at two collocation levels on phenol conversion (Y1) and hydrogen yield (Y2).An interaction exists between variables when the effect of one variable on the response depends on the level of another variable [59].Also, the presence of interaction results in the combined effect being more or less than that expected when the effects are individually added together [60].Figure 10 shows that phenol concentration factor (E) is a predominant factor in Y1 and Y2.This means that E interacted strongly with all other factors.There

Interaction Plots for Y 1 and Y 2
To visualize the impact of each factor combination and identify which factors were most influential, interaction plots were constructed.These also consider interactions among variables and can be used to optimize the operating parameters in multivariable systems [57].The mean response of the two factors is depicted by the plot for all the possible combinations of their settings.When there is interaction between the two factors, the lines appear as non-parallel [58].Figure 10 presents the full interactions effect plot matrix for five factors, temperature (A), liquid feed rate (B), catalyst weight (C), Ni to Co ratio (D), and phenol concentration (E) at two collocation levels on phenol conversion (Y 1 ) and hydrogen yield (Y 2 ).An interaction exists between variables when the effect of one variable on the response depends on the level of another variable [59].Also, the presence of interaction results in the combined effect being more or less than that expected when the effects are individually added together [60].Figure 10 shows that phenol concentration factor (E) is a predominant factor in Y 1 and Y 2 .This means that E interacted strongly with all other factors.There was a strong two-way interaction between the main effects of A and E, B and E, C and E, D and E, and C and E, as evident from the non-parallel nature of the effect lines for the interaction on Y 1 and Y 2 responses.However, it can be proved by regression coefficient analysis.The graph also shows that there is a possible interaction between C and B factors with respect to Y 1 and Y 2 due to their un-parallel lines.The rest of the graphs are parallel, suggesting a general lack of interaction or weak interaction between the different variables with respect to their influence on Y 1 and Y 2 .

Surface and Contour Plots for Phenol Conversion (Y1) Response
As is seen from the analysis of the regression coefficient results in Table 5, all of the variables (temperature (A), feed flow rate (B), catalyst weight (C), Ni-Co ratio (D), and phenol concentration (E)) and their interactions (AE, BE, CE, and DE) are statistically significant and affect the phenol conversion (Y1) response.The response surface plots are useful in understanding the relationship between the main and interaction effects of the variables [61].The contour plots are curved lines because the model contains the interactions of the factors (AE, BE, CE, and DE).The surface and contour plots given in Figures 11 and 12 show the relative effects of any two variables when the remaining variables are kept constant.In general, the results of surface and contour plots reveal that the phenol concentration plays an important role in catalyst activity towards phenol conversion.The surface and contour plots for AE in Figures 11a and 12a show the interactive effect of phenol concentration (E) and temperature (A) by keeping B, C, and D constant.The phenol conversion is relatively high for high levels of A and E. In the same way, the surface and contour plot show the interactive effect of B and E through keeping A, C, and D constant.Figures 11b and 12b show the Y1 response of the mutual interactions between the variables.It was found to be high at a low flow rate and high phenol concentration.Furthermore, when the catalyst weight ranged from 0.1 to 0.3 g the effect on the Y1 response was significant (Figures 11c and 12c).On the other hand, Figures 11d and  12d show that the NiCo ratio and phenol concentration (D versus E) have a negative effect on Y1 response.In conclusion, surface and contour plots prove the results of the regression coefficient data and the Pareto charts.

Surface and Contour Plots for Phenol Conversion (Y 1 ) Response
As is seen from the analysis of the regression coefficient results in Table 5, all of the variables (temperature (A), feed flow rate (B), catalyst weight (C), Ni-Co ratio (D), and phenol concentration (E)) and their interactions (AE, BE, CE, and DE) are statistically significant and affect the phenol conversion (Y 1 ) response.The response surface plots are useful in understanding the relationship between the main and interaction effects of the variables [61].The contour plots are curved lines because the model contains the interactions of the factors (AE, BE, CE, and DE).The surface and contour plots given in Figures 11 and 12 show the relative effects of any two variables when the remaining variables are kept constant.In general, the results of surface and contour plots reveal that the phenol concentration plays an important role in catalyst activity towards phenol conversion.The surface and contour plots for AE in Figures 11a and 12a show the interactive effect of phenol concentration (E) and temperature (A) by keeping B, C, and D constant.The phenol conversion is relatively high for high levels of A and E. In the same way, the surface and contour plot show the interactive effect of B and E through keeping A, C, and D constant.Figures 11b and 12b show the Y 1 response of the mutual interactions between the variables.It was found to be high at a low flow rate and high phenol concentration.Furthermore, when the catalyst weight ranged from 0.1 to 0.3 g the effect on the Y 1 response was significant (Figures 11c and 12c).On the other hand, Figures 11d and 12d show that the NiCo ratio and phenol concentration (D versus E) have a negative effect on Y 1 response.In conclusion, surface and contour plots prove the results of the regression coefficient data and the Pareto charts.
and high phenol concentration.Furthermore, when the catalyst weight ranged from 0.1 to 0.3 g the effect on the Y1 response was significant (Figures 11c and 12c).On the other hand, Figures 11d and  12d show that the NiCo ratio and phenol concentration (D versus E) have a negative effect on Y1 response.In conclusion, surface and contour plots prove the results of the regression coefficient data and the Pareto charts.13 and 14 illustrate that E is more significant with a lower level of B and D, and a higher level of A and C factors.The effect of two-way interaction on Y2 response is almost the same as the Y1 response and proves the analysis of the regression coefficient data in Table 6.13 and 14 illustrate that E is more significant with a lower level of B and D, and a higher level of A and C factors.The effect of two-way interaction on Y 2 response is almost the same as the Y 1 response and proves the analysis of the regression coefficient data in Table 6.
temperature (A), feed flow rate (B), catalyst weight (C), Ni-Co ratio (D), and phenol concentration (E) upon the hydrogen yield (Y2) response.As in the previous section, the highest Y2 values are observed at the highest phenol concentration.Figures 13 and 14 illustrate that E is more significant with a lower level of B and D, and a higher level of A and C factors.The effect of two-way interaction on Y2 response is almost the same as the Y1 response and proves the analysis of the regression coefficient data in Table 6.

Response Optimization
Finally, the Minitab Response Optimizer tool helps to obtain the best possible result for phenol conversion (Y1) and hydrogen yield (Y2).Moreover, it can help solve problems involving several design variables and one or more responses simultaneously.It also has an interactive graphic mode that helps in investigating local behaviors around an optimized solution.By using a response optimizer, the optimal conditions for the suggested reduced model (Equations ( 8) and ( 9)) were found to be 800 °C, 0.16 mL/min feed flow rate, 0.

Response Optimization
Finally, the Minitab Response Optimizer tool helps to obtain the best possible result for phenol conversion (Y 1 ) and hydrogen yield (Y 2 ).Moreover, it can help solve problems involving several design variables and one or more responses simultaneously.It also has an interactive graphic mode that helps in investigating local behaviors around an optimized solution.By using a response optimizer, the optimal conditions for the suggested reduced model (Equations ( 8) and ( 9)) were found to be 800 ˝C, 0.16 mL/min feed flow rate, 0.3 g catalyst weight, 0 Ni/Co ratio, and 10 wt % of phenol concentration in order to obtain 94.98% of Y 1 and 67.4% of Y 2 responses (Figure 15a,b, respectively).
design variables and one or more responses simultaneously.It also has an interactive graphic mode that helps in investigating local behaviors around an optimized solution.By using a response optimizer, the optimal conditions for the suggested reduced model (Equations ( 8) and ( 9)) were found to be 800 °C, 0.16 mL/min feed flow rate, 0.3 g catalyst weight, 0 Ni/Co ratio, and 10 wt % of phenol concentration in order to obtain 94.98% of Y1 and 67.4% of Y2 responses (Figure 15a,b, respectively).

Coke Analysis
The changes occurring during the steam reforming of Ni/ZrO2 and Co/ZrO2 catalysts were studied through thermo-gravimetric analyses (TGA) of the used catalysts for different phenol concentrations.As is shown in Figure 16, they had different amounts corresponding to the weight

Coke Analysis
The changes occurring during the steam reforming of Ni/ZrO 2 and Co/ZrO 2 catalysts were studied through thermo-gravimetric analyses (TGA) of the used catalysts for different phenol concentrations.As is shown in Figure 16, they had different amounts corresponding to the weight loss associated with decomposition/oxidation of graphite from the catalyst.The carbon deposit types depend carbon oxidation temperature.The curves of weight loss as a function of temperature can be divided into three parts: below 200 ˝C (WL 1 ), 200-600 ˝C (WL 2 ), and above 600 ˝C (WL 3 ).In the literature, the region below 200 ˝C has been attributed to water and moisture evaporation.The weight loss between 200 ˝C and 600 ˝C has been ascribed to less stable deposits [62], whereas the weight loss above 600 ˝C is generally attributed to oxidation of coke deposits with a more stable [62] and/or different degree of graphitization [63].In general, Ni/ZrO 2 has less coke deposited on its surface than Co/ZrO 2 .The TGA result shows that the Ni/ZrO 2 catalyst for 2% phenol concentration in the feed did not suffer from any coke formation on its surface in all WL 1 , WL 2 , and WL 3 regions.This catalyst profile exhibited additional weight gained; this might be due to the oxidation of the metallic active sites [64,65].A significant mass loss was found for both Ni/ZrO 2 and Co/ZrO 2 catalysts when the phenol concentration increased from 2% to 10% and the temperature rose above 600 ˝C due to the increase in carbon molecules in the feed.Therefore, a higher phenol concentration leads to higher catalyst activity in terms of phenol conversion and hydrogen production but causes the carbon deposition to increase.loss associated with decomposition/oxidation of graphite from the catalyst.The carbon deposit types depend carbon oxidation temperature.The curves of weight loss as a function of temperature can be divided into three parts: below 200 °C (WL1), 200-600 °C (WL2), and above 600 °C (WL3).In the literature, the region below 200 °C has been attributed to water and moisture evaporation.The weight loss between 200 °C and 600 °C has been ascribed to less stable deposits [62], whereas the weight loss above 600 °C is generally attributed to oxidation of coke deposits with a more stable [62] and/or different degree of graphitization [63].In general, Ni/ZrO2 has less coke deposited on its surface than Co/ZrO2.The TGA result shows that the Ni/ZrO2 catalyst for 2% phenol concentration in the feed did not suffer from any coke formation on its surface in all WL1, WL2, and WL3 regions.This catalyst profile exhibited additional weight gained; this might be due to the oxidation of the metallic active sites [64,65].A significant mass loss was found for both Ni/ZrO2 and Co/ZrO2 catalysts when the phenol concentration increased from 2% to 10% and the temperature rose above 600 °C due to the increase in carbon molecules in the feed.Therefore, a higher phenol concentration leads to higher catalyst activity in terms of phenol conversion and hydrogen production but causes the carbon deposition to increase.

Conclusions
In this study, hydrogen production from catalytic steam reforming of phenol was performed over Ni/ZrO2 (Ni-Co ratio = 1) and Co/ZrO2 (Ni-Co ratio = 0) catalysts.In order to investigate the promising catalyst activity towards phenol conversion (Y1) and hydrogen yield (Y2), as well as reduce the total number of experiments, a statistical design of experiment (DOE) using the 2 5 full

Conclusions
In this study, hydrogen production from catalytic steam reforming of phenol was performed over Ni/ZrO 2 (Ni-Co ratio = 1) and Co/ZrO 2 (Ni-Co ratio = 0) catalysts.In order to investigate the promising catalyst activity towards phenol conversion (Y 1 ) and hydrogen yield (Y 2 ), as well as reduce the total number of experiments, a statistical design of experiment (DOE) using the 2 5 full factorial method (five variables, each variable at two levels) was performed.The five variables were temperature (A), liquid feed flow rate (B), catalyst weight (C), Ni to Co ratio (D), and phenol concentration (E).Thirty-two sets of experiments were performed to explore the effects of the reaction variables and their interaction on the phenol conversion and the hydrogen yield.The catalyst performance testing shows that 800 ˝C, 0.16 mL/min feed flow rate, 0.3 g of catalyst, 0 ratio of Ni-Co (Co/ZrO 2 ), and 10 wt % of phenol in the feed gave a phenol conversion of 96.8% and a hydrogen yield of 72.4%.The statistical DOE indicated that within the limits that were being examined, all the variables had a significance level of 5% (p < 0.05).In addition, the interaction among A versus E, B versus E, C versus E, and D versus E were found to have a significant effect on Y 1 and Y 2 responses.This was proven by Pareto charts and Normal plots.The analysis of the regression coefficient results for the reduced regression model indicates that the model can fit the experimental data for both Y 1 and Y 2 responses with correlation coefficients (R 2 ) of 0.968 and 0.936, respectively.The parity charts of predicted versus observed responses also confirm the fit of the experimental data with the regression model.The p-value of the reduced model was zero and that model is considered to be statistically significant.The surface and contour plots show that the E variable is more significant with higher levels of A and C and lower levels of B and D. The response optimizer tool indicated that the optimum conditions for the reduced model were 800 ˝C, 0.16 mL/min, 0.3 g catalyst, 0 Ni-Co ratio, and 10 wt % of phenol concentration to obtain 94.98% of Y 1 and 67.4% of Y 2 , respectively.

Figure 1 .
Figure 1.Schematic flow diagram of experimental setup.

Figure 1 .
Figure 1.Schematic flow diagram of experimental setup.

Figure 6 .
Figure 6.Parity charts of predicted versus observed responses for the initial models of (a) phenol conversion (Y1) and (b) hydrogen yield (Y2), including all main and interaction independent variables.

Figure 6 .
Figure 6.Parity charts of predicted versus observed responses for the initial models of (a) phenol conversion (Y 1 ) and (b) hydrogen yield (Y 2 ), including all main and interaction independent variables.

Figure 7 .
Figure 7. Pareto chart of the reduced model for (a) Y1 and (b) Y2.

Figure 7 .
Figure 7. Pareto chart of the reduced model for (a) Y 1 and (b) Y 2 .

Figure 8 .
Figure 8. Parity charts of predicted versus observed responses for the reduced models of phenol conversion (Y1) (a) and hydrogen yield (Y2) (b) with a 95% prediction band.

Figure 8 .
Figure 8. Parity charts of predicted versus observed responses for the reduced models of phenol conversion (Y 1 ) (a) and hydrogen yield (Y 2 ) (b) with a 95% prediction band.

3. 4 . 4 .
Figures 13 and 14 give surface and contour plots indicating the effect of interactive effect temperature (A), feed flow rate (B), catalyst weight (C), Ni-Co ratio (D), and phenol concentration (E) upon the hydrogen yield (Y 2 ) response.As in the previous section, the highest Y 2 values are observed at the highest phenol concentration.Figures13 and 14illustrate that E is more significant with a lower level of B and D, and a higher level of A and C factors.The effect of two-way interaction on Y 2 response is almost the same as the Y 1 response and proves the analysis of the regression coefficient data in Table6.

Figure 15 .
Figure 15.Results from the Response Optimizer tool for Y1 (a) and Y2 (b).

Figure 15 .
Figure 15.Results from the Response Optimizer tool for Y 1 (a) and Y 2 (b).

Table 1 .
Selected factors and levels for the experimental study.

Table 2 .
Experiment matrix for the 25factorial designs and experimental results.

Table 2 .
Cont.XRD analysis was used to identify the crystal phases of the catalysts.The XRD patterns of the Ni/ZrO 2 and Co/ZrO 2 catalysts after reduction at 600 ˝C are displayed in Figure2.The XRD spectrum was characterized by several intense peaks between the diffraction angles of 20 ˝and 80 ˝.

Table 3 .
Analysis of the regression coefficient of the initial model for phenol conversion (Y1).

Table 3 .
Analysis of the regression coefficient of the initial model for phenol conversion (Y 1 ).

Table 4 .
Analysis of the regression coefficient of the initial model for hydrogen yield (Y 2 ).

Table 5 .
The regression coefficient value of the reduced model for Y1.

Table 5 .
The regression coefficient value of the reduced model for Y 1 .

Table 6 .
The regression coefficient value of the reduced model for Y 2 .

Table 6 .
The regression coefficient value of the reduced model for Y2.