Research Article - (2014) Volume 4, Issue 2

Application of Box Behnken Experimental Design to Analysis Reaction Conversion of n-Butyl Acetate to n-Butanol

Firuz Fakheri*, Ali Sharif Alam, Shokrollah Mohseni and Ramin Ramezani Kalhor
Shifa Pharmed Industrial Group, Karaj, Iran, E-mail: alamsharifali@gmail.com
*Corresponding Author: Firuz Fakheri, Shifa Pharmed Industrial Group, Karaj, Iran, Tel: 9367341446 Email:

Abstract

In this paper Box Behnken was used to study the effect and optimization of reaction temperature, percentage of n -butyl acetate, value of sodium hydroxide on the percentage of reaction conversion of n -butyl acetate to n -butanol. The mathematical relationship of percentage of reaction conversion on three significant independent variables can be approximated by a nonlinear polynomial model. The analysis of experimental design showed that the percentage of reaction conversion was increased by increasing the temperature and value of n -butyl acetate. Results show that the excess value of sodium hydroxide has less effect on the percentage reaction conversion than the other variables and the percentage reaction conversion was found to be strongly affected by the temperature of reaction.

Introduction

Organic solvents are widely used in the pharmaceutical industry. Solvent recovery in pharmaceutical industry is an important plant to reduce operating cost and environmental protection. Purification and recovery of solvents is essential for impurity-free sensitivity standards. During the recovery process chemical reaction is used to remove n-butyl acetate from waste n-butanol. n-butyl acetate cannot be completely removed from n-butanol-water mixtures by distillation because of the presence of the minimum ternary azeotrope [1-3].

In this study the parameters affecting the percentage reaction conversion in constant conditions of the experiment are: temperature of reaction (C°), percentage of n-butyl acetate (v/v), value of sodium hydroxide (ml).

The effect of these parameters may not be independent of each other and it is necessary to consider their interactions. Considerable research has been conducted for determining optimal process parameters. The optimization of the reaction variable has been carried out by using multivariate statistic techniques [4]. Design and methods such as Box Behnken, factorial design and Taguchi methods are now widely used in place of the one factor at a time experimental approach, which is time consuming and exorbitant in cost.

Box Behnken summarized as a collection of statistical tools and techniques for constructing and exploring an approximate functional relationship between an objective function (percentage of reaction conversion) and a set of independent variables. In this technique, the main purpose is to optimize the objective function that is influenced by various operating variables [5].

In this work the objective was to use the Box Behnken Design (BBD) for analyzing and optimizing the effects of the temperature of reaction (C°), percentage of n-butyl acetate (v/v), value of sodium hydroxide(ml) on the percentage of reaction conversion by using twophase stirred.

Materials and Methods

Box behnken

Response surface is as a function of independent variables where the response surface can be expressed as follows:

Y = f (x1,x2,....,xn) (1)

The goal is to optimize the response variable Y (percentage of reaction conversion) and xi, the variables of action called factors. An important assumption is that the operating variables are continuous and controllable by experiments with negligible errors. The design consists of factorial design points with k (k +1) runs and n replications at the center point of the design [4-6].

Among the many types of experimental design models, Box Behnken are commonly used to analysis of factors in three levels three factor experiments. It is necessary to find a suitable approximation for the true functional relationship between independent variables and the response surface. Usually a second-order polynomial Equation 2 is utilized in response surface methodology and is expressed as [7-8]:

image (2)

Where Y is the percentage reaction conversion, b0 is the value of fitted response at the center point of design, b1…3, b11…33 and b12,13,23 are the linear, quadratic and interaction terms respectively, and x1…3 the dimensionless coded variables. The second-order model includes all the terms in the first-order model, plus all quadratic terms like b11x12 and all cross product terms like B13x1x3 [6,9].

Experimental design

According to the Table 1 the hydrolysis of n-butyl acetate to n-butanol with sodium hydroxide solution was investigated. In this study hydrolysis reaction with 20% sodium hydroxide aqueous solution was used to convert n-butyl acetate to n-butanol. In constant state of reaction the time duration of hydrolysis reaction was 120 minutes. Gas chromatography applied to the quantitative determination of volume percent n-butyl acetate in waste water-butanol.

Response surfaces were generated by MINITAB (Demo version 14) statistical analysis software. Statistical analysis of the model was performed to evaluate the analysis of variance (ANOVA). Based on the preliminary experimental results, the levels chosen for the independent variables, Temperature of Reaction (x1), Volume of NaOH (x2) and Percent of BA (x3) were shown in Table 1. Box Behnken design (BBD) of the experiments was used, where the values of independent variables were coded as the variables, x, in the range of +1 and – 1 levels. This design was composed of a 12 factorial design (runs 1–12), 3 center points (runs 13–15) thus 15 experiments were needed in total. The mean value of the response (percentage of reaction conversion) obtained under the different experimental conditions are summarized in Table 2. The experimental data (Table 2) were analyzed using statistical methods appropriate to the experimental design used. Replications of factor combinations were necessary to estimate the experimental error [7].

Temperature of Reaction (C°) Volume of NaOH (ml) Percent of BA% (v/v)
Real value X1 Coded value x1 Real value X2 Coded value x2 Real value X3 Coded value x3
80 1 60 1 10 1
50 0 40 0 7 0
20 -1 20 -1 4 -1

Table 1: Experimental range and levels of the independent variable.

Run order Block x1 x2 x3 percentage of reaction conversion (%)
Experimental Predicted
1 1 0 -1 -1 14 17
2 1 1 1 0 100 95.2
3 1 -1 0 -1 9 1.25
4 1 -1 1 0 50 59.7
5 1 1 0 1 91 98.7
6 1 -1 -1 0 44 48.7
7 1 -1 0 1 95 88.2
8 1 0 -1 1 74 76
9 1 0 1 -1 27 25
10 1 1 -1 0 96 86.2
11 1 1 0 -1 57 63.7
12 1 0 1 1 91 88
13 1 0 0 0 84 84
14 1 0 0 0 84 84
15 1 0 0 0 84 84

Table 2: Parameter levels of Box Behnken (coded value) and the results of experimental and predicted values for percentage of reaction conversion.

Codification of the levels of the variable consists of transforming each studied real value into coordinates inside a scale with dimensionless values. The variable levels Xi were coded as xi according to the following equation such that X0 corresponded to the central value:

image (3)

Where xi is the dimensionless value of an independent variable; Xi is the real value of an independent variable; X0 is the real value of an independent variable at the center point and ΔXi is the step change [9].

According to Equation. 3, the relations between the coded and the real value were as given below:

image (4)

Results and Discussion

Estimated regression coefficient for response and analysis of variance (ANOVA)

The correlation coefficient may not be sufficient by itself to check the fitted. Thus, ANOVA analysis was also done in this work.

Importance of each coefficient was specified by T-values and P-values are listed in Table 3. The larger the magnitude of T-value and the smaller the P-value indicates the considerable significance of the corresponding coefficient [10].

Term Coef SE Coef T P
Constant 84.00   5.612 14.967 0.000
x1 30.00 3.437  8.874 0.000
x2 5.00 3.437  1.455 0.205
x3 21.75   7.719  2.818 0.037
x1*x1 -21.00 5.059 -4.151 0.009
x2*x2 -11.50  5.059 -2.273 0.072
x3*x3 0.00   5.059 0.000 1.000
x1*x2 1.00 4.861  0.206 0.845
x1*x3 -13.00 4.861 -2.675 0.044
x2*x3 -0.50   4.861 -0.103 0.922

Table 3: Results of the regression analysis of the BBD.

The results of the quadratic and linearity model in the form of ANOVA showed small probability value (P<0.1) indicating that the all of terms in the model have a significant effect (Table 3).

It can be seen that the variables with largest effect were the linear terms (x1 and x3) and the quadratic effect of value of sodium hydroxide (ml). Furthermore, the linear effect of the temperature of reaction (C°), percentage of n-butyl acetate (v/v), the quadratic effect of value of sodium hydroxide (ml) and temperature of reaction (C°) and interaction terms (x12 and x13) are more significant than the other factors, which can be neglected and eliminated from the model. Multiple regression analysis of the experimental data gave the following second order polynomial equation with a 95% confidence:

image (5)

Table 2 shows the results of Box Behnken design experiments for studying the effect of three independent variables, along with the predicted mean and observed responses.

The fit of the regression model attained was checked by the coefficient of determination R2. The following Figure 1 demonstrates a high degree of correlation between the experimental and predicted values. The value of the determination coefficient R2 (965%) suggests that only about 3.5% of the total variation is not explained by the model.

advanced-chemical-engineering-Experimental-values

Figure 1: Experimental values and values predicted by Equation (5).ectively.

The normal probability plotting of the residuals can be used to checking normality of the data [11]. Figure 2 shows the normal probability plot for the empirical model for percent of conversion. The data shows a straight line in the plot of the residuals, which represents a normal distribution and thus supports the adequacy of the leastsquares fit.

advanced-chemical-engineering-Normal-probability

Figure 2: Normal probability plot of the residuals.

Response surface plots

The purpose application of the Response Surface Method is to increase efficiency of hydrolysis reaction (Percent of conversion) and to understand how the Percent of conversion changes in a given direction by adjusting the design variables. In general, the response surface can be visualized graphically.

The function percent of conversion can be plotted versus the levels of temperature of reaction (C°), percentage of n-butyl acetate (v/v), and value of sodium hydroxide as shown as Figures 3a-3c. The curvatures in the 3D plots arise due to the quadratic dependency on response and parameters.

advanced-chemical-engineering-Response-Surface

Figure 3a: Effect of Percent of BA and Volume of NAOH on the Percent of Conversion (Y) in Response Surface (3D) Plot.

advanced-chemical-engineering-Temperature-Reaction

Figure 3b: Effect of Percent of BA and Temperature of Reaction on the Percent of Conversion (Y) in Response Surface (3D) Plot.

advanced-chemical-engineering-Temperature-Reaction

Figure 3c: Effect of Volume of NAOH and Temperature of Reaction on the Percent of Conversion (Y) in Response Surface (3D) Plot.

Figure 3a shows the effect of value of sodium hydroxide and percentage of n-butyl acetate on efficiency of hydrolysis reaction (Percent of conversion).

The rate of conversion reaction can be increased by increasing the value of sodium hydroxide, except value of sodium hydroxide of about X2=50 ml and higher.

The percent of conversion reached the maximum value in the middle level of value of sodium hydroxide. Addition of excess sodium hydroxide wills insignificant decrease the rate of conversion reaction and it caused formation of large amounts of the sodium acetate sediment in stirred reactor. The rate of conversion reaction can be increased by increasing the percentage of n-butyl acetate (Figure 3a).

The maximum value of conversion reaction was detected in a high level of temperature of reaction and percentage of n-butyl acetate. Results show that the temperature of reaction has more effect on the rate of conversion reaction than the other variables of the hydrolysis reaction (Figures 3b-3c).

Sometimes it is less complicated to view the response surface in two-dimensional graphs. The contour plots can show the contour lines of the independent variables that have the same response value Y. From the response surface contour plots, it is very easy and convenient to understand the interactions between two independent variables and also to locate their optimum levels.

Figures 4a-4c shows the contour plots of percent of conversion for each pair of parameters. An elliptical contour plot indicates that the interactions between them are significant. In Figure 4b indicating that there were significant interactions between the temperature of reaction and percentage of Butyl Acetate.

advanced-chemical-engineering-Percent-Volume

Figure 4a: Effect of Percent of BA and Volume of NAOH on the Percent of Conversion (Y) in Contour (2D) Plot.

advanced-chemical-engineering-Conversion-Contour

Figure 4b: Effect of Percent of BA and Temperature of Reaction on the Percent of Conversion (Y) in Contour (2D) Plot.

advanced-chemical-engineering-Temperature-Reaction

Figure 4c: Effect of Volume of NAOH and Temperature of Reaction on the Percent of Conversion (Y) in Contour (2D) Plot.

Conclusions

In this work, the effects of three experimental parameters on the percent of conversion were investigated statistically. BBD was used to design an experimental program for to analysis reaction conversion of n-butyl acetate to n-butanol. In this model the parameters affecting the reaction conversion in constant conditions of the experiment were: the temperature of reaction (x1), value of sodium hydroxide (x2) and percentage of n-butyl acetate (x3) . Mathematical model equations were derived for the percent of conversion by using the experimental data and the mathematical software package MATLAB (Demo version 14). The adequacy of the predictive model was verified effectively by the validation data. Predicted values were found to be in good agreement with the experimental values. In order to gain a better understanding of the three variables for optimal percent of conversion performance, the models were presented as 3-D response surface graphs. The optimum process parameters for percent of conversion were obtained as: temperature of reaction, 65°C; value of sodium hydroxide, 40 ml and percentage of n-butyl acetate, 8% v/v. The values of R-squared (R2=0.965) is close to 1, indicating a high degree of correlation between the response and the independents variables in two responses (experimental and predicted values). The value of the determination coefficient R2 (965%) suggests that only about 3.5% of the total variation are not explained by the model.

References

Citation: Fakheri F, Alam AS, Mohseni S, Kalhor RR (2014) Application of Box Behnken Experimental Design to Analysis Reaction Conversion of n-Butyl Acetate to n-Butanol. J Adv Chem Eng 4:110.

Copyright: 2014 Fakheri F, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.