Design and Optimization of a Centrifugal Pump for Slurry Transport Using the Response Surface Method

: Since centrifugal pumps consume a mammoth amount of energy in various industrial applications, their design and optimization are highly relevant to saving maximum energy and increasing the system’s efﬁciency. In the current investigation, a centrifugal pump has been designed and optimized. The study has been carried out for the speciﬁc application of transportation of slurry at a ﬂow rate of 120 m 3 /hr to a head of 20 m. For the optimization process, a multi-objective genetic algorithm (MOGA) and response surface methodology (RSM) have been employed. The process is based on the mean line design of the pump. It utilizes six geometric parameters as design variables, i.e., number of vanes, inlet beta shroud, exit beta shroud, hub inlet blade draft, Rake angle, and the impeller’s rotational speed. The objective functions employed are pump power, hydraulic efﬁciency, volumetric efﬁciency, and pump efﬁciency. In this reference, ﬁve different software packages, i.e., ANSYS Vista, ANSYS DesignModeler, response surface optimization software, and ANSYS CFX, were coupled to achieve the optimized design of the pump geometry. Characteristic maps were generated using simulations conducted for 45 points. Additionally, erosion rate was predicted using 3-D numerical simulations under various conditions. Finally, the transient behavior of the pump, being the highlight of the study, was evaluated. Results suggest that the maximum ﬂuctuation in the local pressure and stresses on the cases correspond to a phase angle of 0 ◦ –30 ◦ of the casing that in turn corresponds to the maximum erosion rates in the region.


Introduction
Centrifugal pumps are being broadly used in coal, chemical, mining, and metallurgical industries mostly to transport solid-liquid mixtures referred to as slurries [1]. These pumps involve three-dimensional flows that are inherently turbulent and comprise secondary flows due to their complex internal passages [2]. Therefore, a good design for slurry transportation has always been a challenge, and several experimental and numerical studies to design and optimize centrifugal pumps can be found in the literature. George et al. performed a numerical study and provided a detailed analysis to minimize cavitation by improving the pump's performance characteristics [3]. In another study [4], a detailed methodology for the inverse design method for a centrifugal impeller using the software HELIOX has been provided. Afterward, the authors also used numerical simulation to validate their results. Shojaeefard et al. [5] performed numerical studies to show the effects of specific geometric parameters on the pump impeller's performance. They concluded that the passage width and blade outlet angle significantly affect the centrifugal pump's performance. Kim et al. [6] performed a design optimization study for centrifugal pumps using numerical calculations. The authors proposed optimized geometry for the pump impeller and pump casing. Nourbakhsh et al. [7] compared the multi-objective particle swarm optimization method with the NSGA II algorithm to design and optimize the centrifugal pump. They tried to optimize variables β hub 1 , β shr 1 , and β 2 . They consequently found important trade-offs in the optimum design of centrifugal pumps which they have proposed based on the Pareto front of two conflicting objective functions. Derakhshan et al. [8] optimized the impeller blade geometry using the artificial bee colony algorithm. They coupled three-dimensional Navier-Stokes optimization with the artificial bee colony algorithm to redesign the geometry and improve the performance of the barked 32-160 pump as a case study. They obtained 3.59% improvement in the efficiency and 6.89 m increase in the pressure head.
In the current study, a centrifugal pump impeller and casing geometry have been designed and optimized by coupling five different software packages, namely, ANSYS Vista CPD, response surface optimization software, DesignModeler, ANSYS TurboGrid, and ANSYS CFX. ANSYS Vista CPD obtains geometry based on required parameters using mean line calculations. The response surface optimization method was implemented to optimize the geometry obtained based on baseline calculations. More than one thousand design points were generated by varying geometrical parameters at the inlet and outlet of the impeller geometry. Three-dimensional geometries of impeller and casing were obtained using ANSYS DesignModeler software, and 3-D meshes of the geometries were generated in the software ANSYS TurboGrid. Finally, numerical simulations were carried out at different operating conditions to get the pump geometry's characteristic maps. The pump geometry has been specifically optimized for the slurry transport of 120 m 3 to a head of 20 m in this study. After achieving the optimized geometry, complete characteristic maps were generated. Moreover, the erosion rate density of the designed pump geometry has been predicted using the Finnie model. Lastly, the pump's transient behavior has been evaluated and reported in the current study, which has never been reported in the literature to the author's best knowledge.

Pump and Volute Casing Design
The pump was designed using Vista CPD, which employs a 1-D approach for the centrifugal pump's preliminary design in the current study. Assumptions and input parameters for the design are listed in Table 1. The pump's efficiency was calculated using built-in correlations in Vista CPD based on the past data for a given range of machinespecific speeds. Definitions of hydraulic efficiency, volumetric efficiency, mechanical efficiency, and pump efficiency are as follows.

Hydraulic Efficiency (η H )
Hydraulic efficiency is the most critical efficiency involved in the pump design that accounts for the head losses based on the hydrodynamic design, i.e., turning losses and friction losses, etc.
Here, H i is the ideal head, while losses based on hydrodynamic design are represented by H loss .

Volumetric Efficiency (η v )
A certain amount of fluid (Q leak ) from the main flow leaks from the tip clearance region to the eye of the pump. Volumetric efficiency in terms of the leaked flow (Q leak ) and actual flow through the outlet (Q) is thus given in the following equation.
where P sha f t is the power applied to drive the pump's shaft and P disk is the power lost due to viscous and mechanical friction.

Pump Efficiency η p
Pump efficiency can be defined as in the equation given below.
The shaft minimum diameter factor was calculated depending upon the maximum allowable shear stress of the shaft. The safety element is applied to the shaft minimum diameter, referred to here as the shaft minimum diameter factor, whereas the term D hub D sha f t is the ratio of hub diameter to shaft diameter and taken as 1.5 in the current study.

Calculation of Leading Blade Angles
Leading-edge blade angles were computed using the Contingent method [9]. These angles were calculated relative to the shroud leading-edge blade angle as given below. It should be noted here that superscripts hub and shr correspond to the parameters related to hub and shrould of the rotor. Tip diameter is the diameter of the impeller at the mean line of the trailing edge. In this study, the tip diameter was computed by specifying the head coefficient (ψ). Head coefficient could be defined as below.
Blade angle at the trailing edge is an essential factor that defines the impeller width at the trailing edge (b 2 ) because the flow rate at the outlet is a function of meridional velocity (C m ) and cross-sectional area at the outlet A out as defined in the equation below. Impeller geometry computed by ANSYS Vista CPD is shown in Figure 1.
coefficient could be defined as below.
Blade angle at the trailing edge is an essential factor that defines the impeller width at the trailing edge ( 2 ) because the flow rate at the outlet is a function of meridional velocity ( ) and cross-sectional area at the outlet as defined in the equation below. Impeller geometry computed by ANSYS Vista CPD is shown in Figure 1.

Specific Velocity
Pump efficiency will be maximum when the value of a specific speed (Ω ) approaches 1. Specific speed (Ω ) is defined in the equation below.
Suction specific speed (Ω ) is a non-dimensional parameter that can be useful in the evaluation of pump cavitation performance. NPSH is the net positive suction head.

Specific Velocity
Pump efficiency will be maximum when the value of a specific speed (Ω s ) approaches 1. Specific speed (Ω s ) is defined in the equation below.
Suction specific speed (Ω ss ) is a non-dimensional parameter that can be useful in the evaluation of pump cavitation performance. NPSH is the net positive suction head.

Pump Impeller Optimization
Pump geometry optimization was carried out using the RSM (response surface method) [10] coupled with a genetic algorithm (GA) [11]. Response surface methodology (RSM) is a tool used for optimization, regression analysis, and understanding the impact of design parameters on the objective function. It is well known and has been frequently adopted in the literature [12] based on its accuracy, trivial computational costs, and ease in implementation [13].
For the current study, the second-degree model of RSM has been adopted as given by Equation (11).
where x i and f (x) denote variables and the system's responses, respectively, while α 0 , β i,i , γ i,j , and γ i are regression coefficients. On the other hand, the term n corresponds to the total number of design variables (e.g., for this study; x 1 → x 6 ) as listed in Table 2. A list of the objective functions with respective constraints is given in Table 3.  Table 3. Response variables and their constraints.

Response Variables Constraints
Shaft power Minimize To develop the design of experiments (DOE) for the current work, a central composite design (CCD) option was utilized. The opted method of the design of the experiment, i.e., CCD, comprises an implanted factorial or fractional factorial design model with center points augmented with a group of "star points". This method permits approximation of curvature [14] if required. If the distance from the center of the design space to a factorial point is ±1 unit for each factor, the distance from the center of the design space to a star point is|α| > 1. The exact value of α rests on defined properties anticipated for the design and the number of factors involved in the system. To uphold rotatability, the assessment of α is influenced by the number of experimental runs in the factorial design part of the CCD: α = (2 n ) 1 4 Further details on the RSM and CCD can be found in the literature [10,[13][14][15]. The list of design variables considered for the pump optimization problem with their upper and lower bounds is provided in Table 2. In contrast, the list of response variables is given in Table 3.
The various design combinations based on the selected design variables ( Table 2) are based on the CCD are listed in Table 4. These designs were constructed and evaluated using Vista CPD. Computed results were later used to perform regression analysis by the response surface methodology. Details for input parameters, output parameters, and constraints are listed in Table 5. Five input parameters with defined ranges were selected, i.e., number of vanes, inlet beta shroud, exit beta shroud, hub inlet blade angle, and the impeller's rotational speed.    Table 5. Details for inlet and outlet parameters.

Range Parameters Constraints
Number of vanes (N)

Multi-Objective Genetic Algorithm
The commonly implemented optimization algorithms in the literature are genetic algorithms [13,[16][17][18][19][20][21] and particle swarm optimization (PSO) [22][23][24]. Both GA and PSO are iteration based and begin with a collection of initial data values and heuristic algorithms; however, the former utilizes discrete data and the latter works on the continuous data. It is registered in the published studies previously that both algorithms are extremely accurate; however, PSO presents comparatively better proficiency. At the same time, genetic algorithm (GA) is more valuable for constraint satisfaction challenges.
In light of the above discussion, the algorithm chosen for optimizing the current study is a multi-objective genetic algorithm (MOGA). Design variables mentioned above with their upper and lower limits have been displayed in Table 2. The solved regression model, i.e., solved Equation (11) for the RSM serves as an objective function for the MOGA. Previously conducted studies for the optimization of the pump have utilized complex and computationally expensive methods. Jiang et al. [25] utilized a CFD-DEM (computational fluid dynamics-discrete element method) approach for the fresh concrete pumping loss optimization. Ping et al. [26] utilized deep learning for the optimization of a multistage centrifugal pump. Similarly, Xiao and Tan [27] optimized pressure fluctuations for a multiphase flow in the pump. On the other hand, this is the first study, to the authors' best knowledge, that utilizes response surface methodology coupled with 3D-RANS (Reynoldsaveraged Navier-Stokes) simulations to compute the objective function that the MOGA would finally utilize to optimize the pump geometry.
Once the unknown coefficients of Equation (11) were found through RSM, it was used as an objective function for the multi-objective genetic algorithm (MOGA) [11]. For the current optimization process, Equation (11) was solved through RSM for shaft power, hydraulic efficiency, volumetric efficiency, and pump efficiency. For the optimization process, pump power was a constraint designated to be minimized, while hydraulic, volumetric, and pump efficiency were constraints to be maximized, as shown in Table 5. MOGA is based on the stochastic method and utilizes the principle of survival of the fittest [16]. The genetic algorithm is initialized through a random population of selected sizes based on the opted design variables that define the first generation. Every individual from the population is tested against the defined fitness function. In this process, each generation's population gets improved through an iterative process using mutation, crossover, and elite selection. A schematic of the process is shown in Figure 2. Population size for the current study was selected as 100 combinations of the defined variables. The crossover fraction and mutation fraction values for the current problem were chosen as 0.9 and 0.01. A 10% elite of the total population was carried out as the next generation unchanged. Figure 3 shows a variation of pump efficiency with rotational speed and hub draft angle. It can be observed that the effects of rotational speed on the pump efficiency are significant, but the hub draft angle has a minimum impact. It can also be noted that pump efficiency increases with the increase of the rotational speed, and then it starts decreasing after attaining a maximum value. A similar phenomenon was observed for all the hub inlet draft angles.

Governing Equations
On the other hand, Figure 4 shows a graph between shaft power, the blade angle beta, and the impeller's rotational speed. It could also be depicted from Figure 4 that the count of vanes has a minimal impact on the power required to run the pump, and again, the shaft power varies significantly with the pump's rotational speed. An increase in impeller rpm from 1400 to 1533 leads to a decrease in impeller power, but a further increase in rotational speed from 1533 to 1800 increases the pump's required power. Table 6 shows the optimized parameter based on response surface optimization. A geometric comparison between base geometry and optimized geometry is shown in Figure 5.   Figure 3 shows a variation of pump efficiency with rotational speed and hub draft angle. It can be observed that the effects of rotational speed on the pump efficiency are significant, but the hub draft angle has a minimum impact. It can also be noted that pump efficiency increases with the increase of the rotational speed, and then it starts decreasing after attaining a maximum value. A similar phenomenon was observed for all the hub inlet draft angles.
Design variables and upper/lower bounds (Table 2) Generate initial population

Governing Equations
On the other hand, Figure 4 shows a graph between shaft power, the blade angle beta, and the impeller's rotational speed. It could also be depicted from Figure 4 that the count of vanes has a minimal impact on the power required to run the pump, and again, the shaft power varies significantly with the pump's rotational speed. An increase in impeller rpm from 1400 to 1533 leads to a decrease in impeller power, but a further increase in rotational speed from 1533 to 1800 increases the pump's required power. Table 6 shows the optimized parameter based on response surface optimization. A geometric compari- beta, and the impeller's rotational speed. It could also be depicted from Figure 4 that the count of vanes has a minimal impact on the power required to run the pump, and again, the shaft power varies significantly with the pump's rotational speed. An increase in impeller rpm from 1400 to 1533 leads to a decrease in impeller power, but a further increase in rotational speed from 1533 to 1800 increases the pump's required power. Table 6 shows the optimized parameter based on response surface optimization. A geometric comparison between base geometry and optimized geometry is shown in Figure 5. In the previous section, an optimized model of the pump has been achieved based on the mean line 1-D computational model. Although 1-D calculations are robust, they are not as precise as 3-D calculations. Thus, this section provides a 3-D computational model to evaluate the performance of the optimized model accurately. Moreover, the design process has been conducted to design a pump specifically for slurry transport. Since the impact of slurry on the pump geometry, i.e., erosion and performance with slurry as a working fluid, could not be predicted using mean line design calculations, a 3-D computational model was necessary.   A steady form of the transport equations of continuity and momentum were solved for both phases (water and sand) using the Eulerian approach. Eulerian description of the equations is given in Equations (11)-(14) [28].
Base design Optimized Figure 5. Comparison between base geometry and optimized geometry.
In the previous section, an optimized model of the pump has been achieved based on the mean line 1-D computational model. Although 1-D calculations are robust, they are not as precise as 3-D calculations. Thus, this section provides a 3-D computational model to evaluate the performance of the optimized model accurately. Moreover, the design process has been conducted to design a pump specifically for slurry transport. Since the impact of slurry on the pump geometry, i.e., erosion and performance with slurry as a working fluid, could not be predicted using mean line design calculations, a 3-D computational model was necessary.
A steady form of the transport equations of continuity and momentum were solved for both phases (water and sand) using the Eulerian approach. Eulerian description of the equations is given in Equations (11)-(14) [28].
Continuity Equation for solid and fluid phase Momentum Equation for solid and fluid phase In the above equation α f and α s are the volume of fluid and solid while U, τ are velocity vector and shear stress tensor, respectively. Everywhere, subscript f and s are used for fluid and solid, respectively. The term p that appears in Equations (13) and (14) is an averaged pressure term and is the same for fluid and solid phases. It is presumed at the interface between the two phases [29]. F l and F vm,s is the lift force and the virtual mass force per unit mass for either phase, respectively. F l , s is the lift force that could be calculated using the formula given by Equation (16) such that F l , s = −F l,f .
Virtual mass effects come into play when the acceleration of the solid phase is observed in the fluid phase, and that results in a virtual mass force F vm,s .
Coupling of solid with fluid could be either one-way or two-way coupling based on the value of the variable β. If the value of β is less than 0.14, the one-way coupling is favored, while if it is more than or equal to 0.14, the two-way coupling is employed. Definition of β is defined in the equation below.
The current study's scope includes the wear prediction that occurs on the blade walls and the casing. The casing is likely to be more exposed to erosion than the blades, where higher tangential velocities tend to erode the casing.
Amongst available Reynolds-averaged Navier-Stokes (RANS) turbulence models, there exist two categories on the basis of how the boundary layer is treated in those models. In the first type of the RANS model, the boundary layer is not resolved, but a wall function is used to estimate the velocity profile, e.g., k − turbulence model. The model uses two transport equations for the solution of turbulence kinetic intensity and dissipation. Simultaneously, turbulent viscosity is modeled as a product of turbulent velocity and turbulent length scale. The model is extremely computationally economical and efficient in resolving the turbulence in the mainstream flow. However, it does not solve the boundary but models it using wall functions. Therefore, it is not considered suitable for flows with boundary layer separation, flows with sudden changes in the mean strain rate, flows in rotating fluids, and flows over curved surfaces. On the other hand, the shear stress turbulence model combines two models [21,[30][31][32], i.e., k − ω and k − turbulence model according to the distance from the wall. Wilcox k − ω is utilized to resolve the flow in the vicinity of the walls for accurate prognostication of the boundary layer, whereas the k − turbulence model is solving the flow in fully developed regions to gain advantage from its robustness, economy, and free stream independence [28,30]. Therefore, the shear stress turbulence model accounts for the transport of turbulent shear stress and avoids over-prediction of eddy viscosity. Numerous researchers [33][34][35][36][37][38] have adopted the shear stress transport (SST) turbulence model established by Menter [30] and achieved accurate aerodynamic predictions of rotating machines and complex flows. Based on the above discussion, the SST turbulence model is adopted for the current work as the flow is rotating. The accurate prediction of the boundary layer is critical for the precise evaluations of the erosion.

Erosion Model
Erosion initiates as particles collide against walls. This phenomenon is quite complicated and is a function of many parameters, such as multiphase effects, the effect of local cavities, properties of solid particles, shape, and size of particles, and local turbulent field. A simplified model suggested by Hutchings [39] shows that erosion rate density is modeled as a function of impact velocity and impact angle. E = kV n p f (γ) (19) In the relation as mentioned above, E characterizes a dimensionless mass, V p represents the velocity of the particle, and f (γ) represents the impact angle in radians. The term n is assigned with different values for different metal ranges from 2.3 to 2.5. It is challenging to develop an erosion model that applies to a wide range of flow conditions and various geometrical features because tuning of model parameters always requires experimental data. Finnie [40] anticipated a model for ductile materials with a value of n = 2 that relates with the kinetic energy of the particles. Finnie's model has been assumed to represent the wear in the current problem.

Computational Domain and Boundary Conditions
The current investigation, impeller, pump, and casing form the computational domain, as shown in Figure 6. Impeller geometry and the hub are modeled as a rotating domain that revolves at 1533 rpm, while the shroud is modeled as a stationary wall. On the other hand, the casing has been modeled as a stationary domain. The frozen rotor has been modeled as an interface between the impeller domain and casing domain for steady-state cases. In contrast, a transient stage interface type is implemented for unsteady simulations. At the inlet of the impeller, the total pressure boundary condition has been imposed, while the mass flow condition has been imposed on the outlet of the casing. Four different flow rates were used, i.e., Q Q η(max) = 0.7, 0.8, 0.9, and 1.0. Here, Q is the volume flow rate while Q η(max) is the flow rate equivalent to the pump's maximum efficiency, i.e., the point for which the rotor has been optimized. A volume concentration of 20% has been used for all the flow rates. Details of the simulations and particle distribution are listed in Table 7. On this account, four sets of steady simulations were conducted with different flow rates.

Mesh Generation
Structured mesh using hexahedral elements was generated for the impeller domain, while unstructured mesh was created for the stationary domain (volute casing). ANSYS TurboGrid was utilized to develop the structured mesh for the impeller domain with controlled topology. Meshes for the impeller domain and casing domain and the topology of the mesh are shown in Figure 7. The O-grid around the blade was constructed to capture the boundary layer accurately.
The value of y + was guaranteed to be less than 1 for both casing and impeller do-   For the unsteady simulation time step (∆t), it was chosen equivalent to the time through which the rotor revolves through 1 • . As mentioned above, the optimized value of the rotational speed of the rotor is 1533 rpm (ω = 25.55 rad/s); thus, time to complete 1 • = ∆t = 0.11 ms. Meanwhile, T can be referred to as time to complete one complete revolution (T = 2Π ω = 0.039 s). For the transient simulations, the total simulation time was chosen as 1.1 s, which corresponds to 10,000 time steps or 27.77 revolutions. It should be noted here that periodicity was achieved after ten revolutions; therefore, data of the first ten rotations was ignored.
It should be noted here that the Eulerian formulation has been used by providing individual densities of the solid phase and the liquid phase as opted by Adnan and Kim [41,42]. Although employing rheological models can improve the accuracy further, they are, however, not included in the current work to reduce the computational cost.

Mesh Generation
Structured mesh using hexahedral elements was generated for the impeller domain, while unstructured mesh was created for the stationary domain (volute casing). ANSYS TurboGrid was utilized to develop the structured mesh for the impeller domain with controlled topology. Meshes for the impeller domain and casing domain and the topology of the mesh are shown in Figure 7. The O-grid around the blade was constructed to capture the boundary layer accurately.  The value of y + was guaranteed to be less than 1 for both casing and impeller domains for all simulations to get maximum advantage to form the SST turbulence model [30].
The value y + is critical to get precise predictions for the SST turbulence model. Few turbulence models, such as the k − turbulence model, do not solve the boundary layer but model it using the wall function. On the contrary, the SST turbulence model resolved the boundary layer, and it is essential to maintain y + < 2 in the whole domain, and there must be fifteen nodes within the boundary layer, at least, for accurate predictions [30]. Consequently, it is ensured while generating the mesh that the value of y + max < 2 is everywhere in the domain. Further, the final mesh contains 27 elements within the boundary layer thickness.
The value of y + is dependent on the value of the ∆y (distance of the first node from the wall), D H (hydraulic diameter), and Reynolds number (Re). Equation (22) could be employed to estimate the value of the ∆y that in turn will be used to compute the value of the y + [43].
Four meshes were generated by varying different parameters listed in Table 1 to optimize the mesh for rotating and stationary domains.
For the mesh optimization, four different meshes were generated with parameters listed in Table 8. The average corrosion rate was set as a parameter for optimization. Computed results for Mesh 3 and 4 matches closely; however, based on less computational power and time requirements, M3 was selected as the final mesh for all the simulations conducted for the current study.

Validation of the Computational Model
Validation of the existing numerical model was completed by assessing the numerical model results utilized in the current study with the experimental research available in [44]. Experimental results introduced in [44] were acquired for a channel flow problem. Dimensions of the channel used were 6 m × 2 m × 0.05 m. The geometrical model of the channel and mesh distribution is shown in Figure 3. Particles used for the study have an average diameter of 165 µm and density of 2680 kg/m 3 . The volume concentration of the particles at the inlet of the domain was 8.41%, while average velocity at the inlet of the domain of the fluid was 1.66 m/s. The same geometry (Figure 8a) and other properties mentioned above used by [44] were implemented for the numerical model discussed above, and a comparison of the results is shown in Figure 8b,c. Figure 8b,c compare the concentration of particles and velocity profile for experimental results and results produced by the current numerical model. The comparison study shows that both works are in close agreement, and the present numerical model could predict erosion rate. The same validation procedure has been adopted in the literature as well by [29]. generated at the 50% span. It could be seen qualitatively in Figure 9 that erosion rate density is much higher at the casing than the erosion rates on the blades of the pump. It can be observed from velocity contours and vectors in Figure 10 that maximum velocity appears to be concentrated near the pump casing. In contrast, the velocities' magnitudes are comparatively low in the flow region near the blades of the impeller. This is why higher erosion rate density can be observed in the casing region and not on the blades of the impeller.
(a) Computational domain and its mesh of the channel geometry used for the validation of the numerical model.

Computational Results
A centrifugal pump for slurry transport has been designed and optimized in the current study using mean line design calculations. To evaluate the optimized pump performance under slurry transport and to develop the characteristic maps of the pump, a 3-D computational model was developed and solved numerically using commercial software ANSYS CFX. Moreover, the pump's transient behavior under slurry flow has been solved to predict the areas under the influence of the maximum erosion rate. Figures 9 and 10 show the qualitative results of pressure, erosion rate, and velocity contours. Figure 9 shows pressure and erosion rate density contours on the walls of both rotating and stationary domains.
In contrast, Figure 10 shows contours and velocity vectors plotted on a turbo surface generated at the 50% span. It could be seen qualitatively in Figure 9 that erosion rate density is much higher at the casing than the erosion rates on the blades of the pump. It can be observed from velocity contours and vectors in Figure 10 that maximum velocity appears to be concentrated near the pump casing. In contrast, the velocities' magnitudes are comparatively low in the flow region near the blades of the impeller. This is why higher erosion rate density can be observed in the casing region and not on the blades of the impeller.

Characteristic Maps of the Optimized Pump Geometry
Characteristic maps of the centrifugal pump were generated using 45 simulations at different values of Q N s and specific speed (Ω s ). Five curves were plotted to complete the characteristic maps for the pump using values of Q N s being 1.0, 0.1, 0.01, 0.001, and 0.001, where Q is the volume flow rate term and N s is the specific speed term. The relation between Ω s and N s is given in the equation below. Each curve was plotted using the results from nine simulations, as shown in Figure 11. It could be seen from Figure 11 that at specific speed values under 0.25, the efficiency of the pump is exceedingly small. Although, efficiency starts rising rapidly as specific speed rises above 0.25. As soon as specific speed approaches 0.5, efficiency increases faintly for a specific speed range from 0.5 to 1.0. Efficiency is maximum at a specific speed Ω s = 1 for all the curves. Moreover, for a specific speed range of 0.5 to 3.0, efficiency of the pump decreases for the curves at different values of Q N s .

Characteristic Maps of the Optimized Pump Geometry
Characteristic maps of the centrifugal pump were generated using 45 simulations at different values of and specific speed (Ω ). Five curves were plotted to complete the characteristic maps for the pump using values of being 1.0, 0.1, 0.01, 0.001, and 0.001, where is the volume flow rate term and is the specific speed term. The relation between Ω and is given in the equation below. Each curve was plotted using the results from nine simulations, as shown in Figure 11. It could be seen from Figure 11 that at specific speed values under 0.25, the efficiency of the pump is exceedingly small. Although, efficiency starts rising rapidly as specific speed rises above 0.25. As soon as specific speed approaches 0.5, efficiency increases faintly for a specific speed range from 0.5 to 1.0. Efficiency is maximum at a specific speed Ω = 1 for all the curves. Moreover, for a specific speed range of 0.5 to 3.0, efficiency of the pump decreases for the curves at different values of .    Figure 12 shows a variation of shear stress along the pump casing wall at various flow rates. Shear stress fluctuates twice at the throat of the casing near point C. It then keeps on increasing until point "D", where it achieves maximum value. After that, it starts to decrease along with the casing at a steady rate. Variation of this shear stress depends on the local velocities and concentration of the sand particles. It can be seen from Figure  11 that wall shear stresses also vary with flow rates. Maximum shear stress was observed for a flow rate corresponding to maximum efficiency, i.e.,   Figure 13 shows a variation of erosion rate density along the casing wall from point A to J. Similar to the shear stress graph shown in Figure 12, erosion rate curves fluctuated near point C, but contrary to the previous graph, it dips down from point C to D. That increases marginally as it moves further. This is because erosion is a complex phenomenon that depends on many parameters, i.e., particle impact angle, particle size, particle velocity, concentration, etc., other than wall shear stress.

Transient Behavior of the Pump under a Slurry Flow Condition
This part of the results elaborates the variations in shear stress and pressure data during one rotor pitch rotation. The impeller of the pump consists of seven blades; thus, the elapsed angle corresponding to one pitch is 360°7 = 51.42°≈ 50. Furthermore, it is mentioned above that 1°of rotation corresponds to one time step (Δ = 0.11 ms); therefore,  Figure 13 shows a variation of erosion rate density along the casing wall from point A to J. Similar to the shear stress graph shown in Figure 12, erosion rate curves fluctuated near point C, but contrary to the previous graph, it dips down from point C to D. That increases marginally as it moves further. This is because erosion is a complex phenomenon that depends on many parameters, i.e., particle impact angle, particle size, particle velocity, concentration, etc., other than wall shear stress.

Transient Behavior of the Pump under a Slurry Flow Condition
This part of the results elaborates the variations in shear stress and pressure data during one rotor pitch rotation. The impeller of the pump consists of seven blades; thus, the elapsed angle corresponding to one pitch is 360 • 7 = 51.42 • ≈ 50. Furthermore, it is mentioned above that 1 • of rotation corresponds to one time step (∆t = 0.11 ms); therefore, one pitch will lapse nearly 51 time steps.
In the current study, result variations for one pitch are reported. Figure 14 shows pressure contours on the turbo surface extracted at 50% of the blade span corresponding to different phases, i.e., 0 • −50 • . When the blade trailing edge is coincident with the line joining the throat of the casing to the axis of rotation, the blade phase angle is taken as θ = 0 • as shown in Figure 14a. It can be seen that for phase angle θ = 0 • − 30 • , maximum pressure at the casing exists then increases at some other locations, and it starts decreasing as the phase angle further increases from θ = 30 • − 50 • . On the other hand, the variation of pressure along the casing is shown in Figure 15. Pressure increases clearly and significantly everywhere in the casing as the phase angle increases from 0 • − 30 • and then it decreases from 30 • − 50 • . The pressure curve in Figure 15 from point C to D shows a quite different distribution compared to all other curves plotted for various phase angles. A maximum difference of pressure was observed as 49.91%, just ahead of point C, between the forces recorded at a phase angle of 0 • and 30 • . This difference was recorded at 9.46% at the exit of the pump casing, reflecting that a fluctuating pressure will be recorded at the exit of the pump with a fluctuating frequency equal to the number of the pump's impeller blades. The magnitude of flection for the current design was observed as 9.46%. Its dependence on other parameters, however, requires further investigations.  A variation of the wall shear stresses and the centrifugal pump's casing at different phases of rotation is displayed in Figure 16. It can be seen that clearly there are three peaks of the wall shear stress-Peek 1 just ahead of C, Peak 2 and Peak 3 just after D. Location of Peak 1 and 2 are same for all phase angles; however, the location of Peak 3 changes slightly with the phase angle.

Conclusions
In the current study, design and optimization of the centrifugal pump geometry has been carried out to transport slurry with a flow rate of 120 m 3 s −1 to a given head of 20 m.

Conclusions
In the current study, design and optimization of the centrifugal pump geometry has been carried out to transport slurry with a flow rate of 120 m 3 s −1 to a given head of 20 m. Five different software packages were coupled to develop the optimization strategy. Optimized geometry was further investigated for performance prediction and erosion prediction under steady and transient conditions. The following points were concluded from the current study.

•
Response surface optimization has proven to be an effective method for the optimization process in the ANSYS workbench environment, where it could be coupled with another meshing software and the Navier-Stokes solver. • Pressure varies slightly, but shear stresses on the wall vary significantly with flow rate variations.

•
Maximum pressure at the casing exits at phase angle θ = 0 • −30 • and also increases at some other locations, and then it starts decreasing as the phase angle further increases from θ = 30 • −50 • • Fluctuating pressure was recorded at the pump's exit with a fluctuating frequency equal to the number of impeller blades of the pump. The magnitude of fluctuation for the current design was observed as 9.46%. Its dependence on other parameters requires further investigation.