In this study, the volumetric power draw P/V was determined as a factor in designing and identifying the optimal condition a successful aeration for stirred wastewater biological treatment vessels. The study was performed to characterize the volumetric power draw in the aerated stirred vessels by optimizing the operation variables. The concerning factors were improved by conjugating stirring and aeration with efficient and economic volumetric power draw condition. The drawn volumetric power was tested and analyzed for three independent parameters; impellers rotation speed (100–200 rpm), turbine blades submergence ratio S/W (0.33–1.67) and wastewater height level ratio H/D (1.37–1.58). A mathematical model was developed in the form of a nonlinear polynomial mathematical model to predict the P/V. The optimal values of the P/V and of relevant parameters were computed through the application of the Box–Wilson technique by application of the central composite rotatable design (CCRD) model. The volumetric power draw P/V and the relevant independent parameters are presented in optimal conditions surface plots that obtained from the nonlinear mathematical model. Optimum analysis result for the independent parameters showed low levels of impellers rotation speed and turbine submergence ratio draw lower P/V while wastewater height did not have a clear effect on P/V.

## LIST OF SYMBOLS

- C
propeller clearance, (m)

- D
rotated cone turbine or impeller diameter, (m)

- d
pitch blade propeller diameter (PBP, U), (m)

- g
acceleration of gravity (9.80665), (m

^{2}s^{−1}) - H
water level in the vessel, (m)

- N
rotation speed, (s

^{−1}) - N
_{p}power number, (P/ ρN

^{3}D^{5}) - S
turbine submergence, (m)

- t
time, (s)

- T
vessel diameter, (m)

- T
_{o}torque measured with filled vessel, (N m)

- T
_{oe}torque measured with empty vessel, (N m)

- V
water volume, (m

^{3}) - W
turbine blade width, (m)

## IMPELLERS TYPES ABBREVIATIONS

## INTRODUCTION

Aeration with stirring (mixing) are considered the most significant processes in gas–liquid dispersion and homogenization in wastewater treatment and in different biotechnological fields (Gimbun *et al.* 2009; Montante *et al.* 2010). Traditionally, stirred vessels are used in unit operations in wastewater biological treatment to increase flocculation rates and to prevent particle settling (Stenstorm & Rosso 2008). In wastewater treatment processes, aeration is needed for large gas handling capacity and an effective gas dispersion for generating as large an interfacial area as possible (Vasconcelos *et al.* 2000). The basic form of the process is a batch operation involves a single vessel which is filled, aerated, and then completely emptied (Kumar & Rao 2010).

Recently, multiple impellers are used in wastewater treatment (Shewale & Pandit 2006). Multiple impellers are preferred over a single impeller as they provide better gas utilization due to the higher gas phase residence time, narrower spread in the residence time distribution in the flow systems and higher surface area per unit liquid volume (Shewale & Pandit 2006). Different multiple impellers types and combinations are used in biological wastewater treatment, the choice depends usually on the way to gas dispersion, top to bottom liquid mixing condition and mass transfer potential created by the system (Stenstorm & Rosso 2008). The flow pattern in multiple impellers stirred tank is remarkably altered than that one impeller is used due to the interactions in circulation loops generated by impellers (Montante *et al.* 2010; Bao *et al.* 2012). In the biological treatment of wastewater, the aeration which accompanied with stirring can be accomplished either via air injection or by some impellers like turbines, which consume energy supplied in the form of mechanical work. Many types of the turbine can be used to achieve the aeration like rotating cones (Adachi 2015). With multiple impellers, an upper rotating cone is usually located close to the wastewater surface to achieve the aeration. While the lower impellers are placed inside the wastewater to assist the air bubbles dispersion, particle suspension, and mixing processes (Li *et al.* 2009).

The power draw of aeration or stirring in wastewater treatment had already been investigated individually by many studies to reduce the consumed energy (Davoody *et al.* 2015; Sun *et al.* 2016). Few of those studies have focused on the power draw practical application of aeration and stirring. Various trends were followed to characterize the drawn power by stirring and aeration in condition applicable for wastewater biological treatment. (Zlokarnik 1979) had investigated the power consumption for four different types of the surface aerator in wastewater treatment. He related the drawn power with the geometrical configuration of the implemented aerators and the wastewater surface condition. Rao *et al.* (2010) reached the same results for the power draw of subsurface aerators, they found that the geometrical configuration such as liquid height, tank diameter, and impeller clearance have most influences. Many studies (Wu 1995; Taghavi *et al.* 2011; Scargiali *et al.* 2013) have obtained practical results on drawn power in aerated and stirred vessels that agree with Zlokarnik (Zlokarnik 1979). They found for systems with multiple impellers, the power draw is usually related with impeller geometry and water surface condition. The performed studies were focused on characterizing and analysis the volumetric power draw in aeration vessels for different geometries, and modes of operations. The power draw has a great dependency on the flow pattern and geometry in aerated vessels of the wastewater (Karimi *et al.* 2013; Wan *et al.* 2016). To breakthrough this relation, many studies were performed to correlate power draw with flow characteristics, the geometrical and operational parameters in aerated stirred vessels (Hiraoka *et al.* 2001; Issa 2016).

Different works have been made for determining optimal process parameters in aeration processes or stirring separately. Generally, the optimum parameter determination may vary as a minimum or a maximum of a parameter due to design consideration. Process parameters optimization is performed by applying multivariate statistical analysis for efficient operation and design purposes. Among the applied techniques for this purpose, the response surface methodology (RSM) or it called factorial design is now widely used. This method is used instead of typical and conventional methods for multifactor experimental design. Various statistical analysis methods failed to detect the true optimum conditions because some of them overlap between parameters effects and other methods need so long time to have results (Tanyildizi *et al.* 2005). To optimize the main response surface in engineering cases, the RSM which mixes mathematical and statistical methods can accomplish the analysis for many affecting parameters of the work (Aslan 2007). Many methods are used to design the experimental run in order modeling and analysis. These methods commonly are the full or partial factorial methods or the central composite rotatable design (CCRD) method. At minimum three levels for each variable are required to estimate the coefficients of quadratic terms in the response model (Box & Wilson 1951). Among the different RSM experimental designs, the CCRD is a more efficient method in order to find the needed information with the minimum required a number of experiments of the dependent variable (Fakheri *et al.* 2012). Considerable works have been conducted the CCRD method in diverse process design (Alalayah *et al.* 2010; Parilti 2010; Kamal *et al.* 2014). By applying CCRD method on aerated and stirred vessel power requirement for different impeller types (Afshar Ghotli *et al.* 2013) figured out this method helps to identify which impeller is more economical and efficient.

In this work the volumetric power draw in a laboratory model, which is geometrically similar to an actual stirred and aerated vessel for biological wastewater treatment is determined in term of three important relevant parameters. These parameters are impellers rotation speed N, turbine submergence ratio S/W, and wastewater level ratio H/D. The aim of this work is to apply the CCRD method for analyzing and optimizing the effects of the impeller rotational speed, wastewater level, and turbine blade submergence on the volumetric power drawn in an aerated and stirred with a dual impeller in wastewater biological treatment vessel. The importance of the chosen relevant parameters comes from their evident effect on the volumetric power drawn that figured out experimentally. This step is necessary to characterize the process drawn energy that dissipated in the wastewater biological treatment vessel during the aerobic process which in turn decides the cost and efficiency of stirring and aeration processes by determining the most affecting parameters.

## METHODS

### RSM

*y*represents the system response variable that needs to be optimized, and the

*x*,

_{1}*x*are action variables, which also they are called factors.

_{2},…., x_{k}An important assumption of this method is that the conducted experiment has insignificant errors when it deals with independent variables in which they are continuous and controllable. The task then is the variable of response surface should have an accurate practical relation with independent variables. The design consists of factorial design points with 2^{k} runs, 2^{k} axial or star points, and n replications at the center of the design (Kwak 2005; Fakheri *et al.* 2012). Usually; RMS method implements a second-order polynomial to represent the mathematical relationships for response variables. Multiple linear regression analysis is applied to estimate the developed mathematical models for each response surface variable. Various terms of a quadratic model are resulting like linear, squared, and interaction (Kwak 2005).

### CCRD

The CCRD was first found by Box and Wilson and then developed by Box and Hunter (Box & Wilson 1951; Box & Hunter 1957). CCRD method requires a fewer number of experiments than the full factorial method, provides much information, and has described the steady-state responses of various operations (Aslan 2008). According to these merits, in this study, it was decided to use CCRD to design and optimize the volumetric power draw P/V in the aerated and stirred wastewater vessel. The affecting parameters under a steady state condition are impeller rotational speed N (rpm), turbine blades submergence ratio S/W and water height level ratio H/D. The required number of experiments in CCRD method covers the standard 2^{3} factorials with its origin at the center, 2^{3} points fixed axially at a distance from the center to generate the quadratic terms, and repeat runs at the center, the axial points are chosen such that they allow rotatability (Aslan 2007). Repeats of the experimental runs at the center level are essential to obtain an independent estimate of the experimental error (Box & Hunter 1957). The optimization of interested parameters may require many tests in ordinary ways to measure the impact of investigated factor on the oxygen transfer rate and to study their interactions; one parameter is varied and the other parameters are kept constant.

### Experimental setup

The types of impellers were used in the study are: a turbine of the type of an inverted pitched rotated cone (IBRC) (used particularly for surface aeration), and a pitch blade propeller, up pumping flow (PBPU). This multiple impeller combination is used typically for mixing and aeration purposes with its up-pumping in the biological wastewater treatment operation, where the turbine (IBRC) is placed at the wastewater surface. in all these cases an acceptable mixing will provide good contact interfacial area between the contents, where the function of the lower propeller (PBPU) is re-directing the liquid flow toward the upper turbine (IBRC) for ensuring continuous feeding intake and to achieve the well-mixing and distribution of the contents of the treatment vessel. The geometrical ratio for the propeller (d/T) = 0.15. A number of blades of the turbine (IBRC) of is 12 (blades width 0.24 m), with diameter ratio D/T = 0.24. Propeller clearances ratio were kept constant with the ratio C/T = 0.2. It is necessary to mention that impellers choice were made on the basis of the requirements found at the real plant and result in a sufficient oxygen supply for the operation conditions in an aerobic wastewater vessel.

### Power consumption calculation

### Experimental design

*X*

_{i}) are coded as a dimensionless values

*x*i while

*X*

_{0}corresponds to the central value of the independent variables and

*ΔX*is the step change (Tanyildizi

_{i}*et al.*2005) as follows.

The coefficients of the yielded second order polynomial were calculated and analyzed by using mathematical model equations that derived by computer simulation programming with the applying least squares method using Matlab R2014a software. Depending on the preliminary experimental results, the chosen levels for the independent variables, impellers rotational speed (*x*_{1}), turbine blades submergence ratio (*x*_{2}), and wastewater height level ratio (*x*_{3}) are listed in Table 1. CCRD method of the experiments design was used, where the values of independent variables were coded as the variables, *x*, in the range of +1 and –1 levels. The dependence of drawn volumetric power P/V on the parameters has been determined with the first and second-grade polynomials like linear and non-linear models.

. | Coded Variables . | Real Variables . | ||||
---|---|---|---|---|---|---|

Run No. . | x_{1} N, (-)
. | x_{2} S/W, (-)
. | x_{3} H/D,(-)
. | X_{1} N, (rpm)
. | X_{2} S/W, (-)
. | X_{3} H/D, (-)
. |

1 | −1 | −1 | −1 | 122 | 0.61 | 1.43 |

2 | −1 | −1 | +1 | 122 | 0.61 | 1.58 |

3 | −1 | +1 | −1 | 122 | 1.39 | 1.43 |

4 | −1 | +1 | +1 | 122 | 1.39 | 1.58 |

5 | +1 | −1 | −1 | 178 | 0.61 | 1.43 |

6 | +1 | −1 | +1 | 178 | 0.61 | 1.58 |

7 | +1 | +1 | −1 | 178 | 1.39 | 1.43 |

8 | +1 | +1 | +1 | 178 | 1.39 | 1.58 |

9 | 0 | 0 | 0 | 150 | 1.0 | 1.47 |

10 | 0 | 0 | 0 | 150 | 1.0 | 1.47 |

11 | 0 | 0 | 0 | 150 | 1.0 | 1.47 |

12 | 0 | 0 | 0 | 150 | 1.0 | 1.47 |

13 | 0 | 0 | 0 | 150 | 1.0 | 1.47 |

14 | −1.73 | 0 | 0 | 100 | 1.0 | 1.47 |

15 | +1.73 | 0 | 0 | 200 | 1.0 | 1.47 |

16 | 0 | −1.73 | 0 | 150 | 0.33 | 1.47 |

17 | 0 | +1.73 | 0 | 150 | 1.67 | 1.47 |

18 | 0 | 0 | −1.73 | 150 | 1.0 | 1.37 |

19 | 0 | 0 | +1.73 | 150 | 1.0 | 1.58 |

. | Coded Variables . | Real Variables . | ||||
---|---|---|---|---|---|---|

Run No. . | x_{1} N, (-)
. | x_{2} S/W, (-)
. | x_{3} H/D,(-)
. | X_{1} N, (rpm)
. | X_{2} S/W, (-)
. | X_{3} H/D, (-)
. |

1 | −1 | −1 | −1 | 122 | 0.61 | 1.43 |

2 | −1 | −1 | +1 | 122 | 0.61 | 1.58 |

3 | −1 | +1 | −1 | 122 | 1.39 | 1.43 |

4 | −1 | +1 | +1 | 122 | 1.39 | 1.58 |

5 | +1 | −1 | −1 | 178 | 0.61 | 1.43 |

6 | +1 | −1 | +1 | 178 | 0.61 | 1.58 |

7 | +1 | +1 | −1 | 178 | 1.39 | 1.43 |

8 | +1 | +1 | +1 | 178 | 1.39 | 1.58 |

9 | 0 | 0 | 0 | 150 | 1.0 | 1.47 |

10 | 0 | 0 | 0 | 150 | 1.0 | 1.47 |

11 | 0 | 0 | 0 | 150 | 1.0 | 1.47 |

12 | 0 | 0 | 0 | 150 | 1.0 | 1.47 |

13 | 0 | 0 | 0 | 150 | 1.0 | 1.47 |

14 | −1.73 | 0 | 0 | 100 | 1.0 | 1.47 |

15 | +1.73 | 0 | 0 | 200 | 1.0 | 1.47 |

16 | 0 | −1.73 | 0 | 150 | 0.33 | 1.47 |

17 | 0 | +1.73 | 0 | 150 | 1.67 | 1.47 |

18 | 0 | 0 | −1.73 | 150 | 1.0 | 1.37 |

19 | 0 | 0 | +1.73 | 150 | 1.0 | 1.58 |

To design the experiments, the operating ranges for three variables is first identified as; impellers rotation speed (100–200 rpm), turbine blades submergence ratio S/W (0.33–1.67) and wastewater height level ratio H/D (1.37- 1.58).

*x*

_{j}) and the corresponding real process variables (

*X*) according to Equation (3) for three variables are as follows:

_{j}For the general form of the quadratic second-degree polynomial model (Equation (5)), the surface contains linear terms in *x*_{n}, square terms in *x*_{nn}^{2}, and cross-product term in *x*_{n}*x*_{j}.

For 2^{3} factorial designs, the total number of the treatment combinations is (2^{n} + 2n + 1). However for three variables, the required experiments are 15 units. The rotatable most likely to be useful in practice belong to a series that is also central. The design is subdivided into three parts; (1) the eight points constitute a 2^{3} factorial, (2) the six points are the extra points included to form a central composite design, and (3) five points are added at the center. According to Equation (5) and upon the completion of each test a total 19 runs were planned to be carried out due to the applied statistical method (see Table 1).

## RESULTS AND DISCUSSION

### Development of a three-variable model depending upon the experimental design

The coded values for the first eight experiments were at two levels −1 and +1 only which are 2^{3} factorial design representing the linear components. Seven extra points were added to the first eight points. These seven points include; one point added the center of the design which was repeated 5 times to estimate the external errors and to facilitate the interpretation of the functional relationship. The other six points were added in pairs along the three coordinate axes at a distance of a = 1.73 from the center of the *x*_{1}, *x*_{2,} and *x*_{3}. This arrangement was found to develop a suitable second order equation fitted to any observed response with studied variables.

The experimental results for the various combinations of the tested variables for the 19 experimental runs are shown in Table 2. Depending on these obtained results the second order polynomial coefficients in Equation (5) were calculated using a computer simulation programming by applying least squares method using Matlab R2014a software. Statistical analysis of the model was performed to evaluate for the predicted P/V by Equation (5) (see Table 2).

. | Coded Variables . | P/V, (watt m^{−3}). | |||
---|---|---|---|---|---|

Run No. . | x_{1} N, (-)
. | x_{2} S/W, (-)
. | x_{3} H/D,(-)
. | Experimental . | Predicted . |

1 | −1 | −1 | −1 | 40.051 | 39.814 |

2 | −1 | −1 | +1 | 44.297 | 45.189 |

3 | −1 | +1 | −1 | 57.112 | 54.993 |

4 | −1 | +1 | +1 | 62.318 | 62.222 |

5 | +1 | −1 | −1 | 70.596 | 70.538 |

6 | +1 | −1 | +1 | 76.567 | 78.533 |

7 | +1 | +1 | −1 | 108.386 | 107.338 |

8 | +1 | +1 | +1 | 117.106 | 117.187 |

9 | 0 | 0 | 0 | 75.534 | 75.752 |

10 | 0 | 0 | 0 | 76.003 | 75.752 |

11 | 0 | 0 | 0 | 75.632 | 75.752 |

12 | 0 | 0 | 0 | 76.043 | 75.752 |

13 | 0 | 0 | 0 | 75.544 | 75.752 |

14 | −1.73 | 0 | 0 | 32.267 | 33.093 |

15 | +1.73 | 0 | 0 | 107.833 | 107.214 |

16 | 0 | −1.73 | 0 | 49.784 | 48.226 |

17 | 0 | +1.73 | 0 | 93.028 | 94.792 |

18 | 0 | 0 | −1.73 | 65.789 | 67.711 |

19 | 0 | 0 | +1.73 | 82.602 | 80.879 |

. | Coded Variables . | P/V, (watt m^{−3}). | |||
---|---|---|---|---|---|

Run No. . | x_{1} N, (-)
. | x_{2} S/W, (-)
. | x_{3} H/D,(-)
. | Experimental . | Predicted . |

1 | −1 | −1 | −1 | 40.051 | 39.814 |

2 | −1 | −1 | +1 | 44.297 | 45.189 |

3 | −1 | +1 | −1 | 57.112 | 54.993 |

4 | −1 | +1 | +1 | 62.318 | 62.222 |

5 | +1 | −1 | −1 | 70.596 | 70.538 |

6 | +1 | −1 | +1 | 76.567 | 78.533 |

7 | +1 | +1 | −1 | 108.386 | 107.338 |

8 | +1 | +1 | +1 | 117.106 | 117.187 |

9 | 0 | 0 | 0 | 75.534 | 75.752 |

10 | 0 | 0 | 0 | 76.003 | 75.752 |

11 | 0 | 0 | 0 | 75.632 | 75.752 |

12 | 0 | 0 | 0 | 76.043 | 75.752 |

13 | 0 | 0 | 0 | 75.544 | 75.752 |

14 | −1.73 | 0 | 0 | 32.267 | 33.093 |

15 | +1.73 | 0 | 0 | 107.833 | 107.214 |

16 | 0 | −1.73 | 0 | 49.784 | 48.226 |

17 | 0 | +1.73 | 0 | 93.028 | 94.792 |

18 | 0 | 0 | −1.73 | 65.789 | 67.711 |

19 | 0 | 0 | +1.73 | 82.602 | 80.879 |

*x*i) and two-factor interactions (

*x*ij) were estimated from the experimental results using a computer simulation applying the method of least squares in Matlab R2014a software. From the experimental results in Table 2, the second-order response functions representing P/V can be expressed as a function of the three operating parameters of P/V, namely the N, S/W and H/D. The relationship between responses drawn volumetric power (P/V) and operating parameters were obtained for coded unit as follows:

^{2}value of 0.9723 for P/V) and F-value equals 148.05.

### Effect of variables on volumetric power draw (P/V)

In Figure 3, the effect of impellers rotation speed and turbine blades submergence on the P/V at the center level of water height can be seen. A higher volumetric power draw can be consumed with increasing the turbine submergence and impellers rotation speed. Figure 4 shows the effect of water height and turbine blades submergence on P/V at center level of rotation speed. The general form of the three-dimensional relationship shows high P/V occurs with a high level of turbine blades submergence for all wastewater height levels.

Figure 4 shows that higher levels of turbine blade submergence have a remarkable effect on P/V. This effect comes from that the immersed blades draw more power to reach the required rotation compared with other conditions. It can be seen also that a higher wastewater height cannot reduce the consumed P/V. Figure 4 shows that the center level of impellers rotation speed is in good condition to obtain a lower P/V.

Figure 5 shows the effect of impellers rotation speed and wastewater height on P/V at center level of turbine blades submergence. It can be seen that as the rotation speed is increased, the P/V is increased progressively. It can be noticed that levels of wastewater height are in acceptable condition to obtain a moderate P/V but the extreme levels of impellers rotation speed are highly power consumable.

From the 3-D response surface figures that impellers rotation speed has a considerable impact on P/V. The minimum level of impellers rotation speed is determined as an optimum to achieve minimum P/V of (0–20 watt/m^{3}). The second important operating parameter is the turbine blades submergence while wastewater height has a minor effect on P/V.

## CONCLUSION

A derived mathematical model was derived by using the experimental data and by application of CCRD for optimizing the surface aerator performance. CCRD was used to design an experimental program for modeling the effects of rotation speed, wastewater level, and turbine blades submergence on the performance of surface aerator volumetric power draw in an agitated tank.

Predicted values of P/V from the developed mathematical model (Equation (6)) were found to be in good agreement with the observed values of P/V (R^{2} value of 0.9545). The results show that impellers rotation speed has a significant effect on the obtained P/V, whilst the turbine blades submergence was less significant and wastewater height has a small effect on P/V. The optimum value of impellers rotation speed is at its minimum level 100 rpm, at this level lower volumetric power was drawn (0–20 watt/m^{3}). Generally, Low turbine submergence ratio keeps the P/V at low levels (0–20 watt/m^{3}). For the third independent variable (H/T) it is hard to distinguish an optimum level for it in term of the other variables. The described analysis tool here will basically help to find the best performance for the used impellers.

## ACKNOWLEDGEMENTS

The author wishes to express his gratitude to the contribution of the LGC-INP Toulouse and the contribution of Biotrade-Toulouse is also gratefully acknowledged by providing the necessary facilities to perform this study.