Research Article  Open Access
Ran Tao, Ruofu Xiao, Wei Yang, Fujun Wang, "A Comparative Assessment of SpalartShur Rotation/Curvature Correction in RANS Simulations in a Centrifugal Pump Impeller", Mathematical Problems in Engineering, vol. 2014, Article ID 342905, 9 pages, 2014. https://doi.org/10.1155/2014/342905
A Comparative Assessment of SpalartShur Rotation/Curvature Correction in RANS Simulations in a Centrifugal Pump Impeller
Abstract
RANS simulation is widely used in the flow prediction of centrifugal pumps. Influenced by impeller rotation and streamline curvature, the eddy viscosity models with turbulence isotropy assumption are not accurate enough. In this study, SpalartShur rotation/curvature correction was applied on the SST  turbulence model. The comparative assessment of the correction was proceeded in the simulations of a centrifugal pump impeller. CFD results were compared with existing PIV and LDV data under the design and low flow rate offdesign conditions. Results show the improvements of the simulation especially in the situation that turbulence strongly produced due to undesirable flow structures. Under the design condition, more reasonable turbulence kinetic energy contour was captured after correction. Under the low flow rate offdesign condition, the prediction of turbulence kinetic energy and velocity distributions became much more accurate when using the corrected model. So, the rotation/curvature correction was proved effective in this study. And, it is also proved acceptable and recommended to use in the engineering simulations of centrifugal pump impellers.
1. Introduction
Reynoldsaveraged NavierStokes (RANS) simulation provides effective solutions for numerical simulations in engineering. By solving the timeaveraged NS equations, RANS simulation reduces the consumption of computing resources on the premise of keeping sufficient accuracy. Based on the Boussinesq assumption [1, 2], eddy viscosity models (EVM) directly established the relationship between Reynolds stress tensor and the traceless mean strain rate tensor. Because of the simplicity and stability, EVMs are widely used in engineering simulations. However, because of the coordinate invariance, EVMs are insensitive to streamline curvature and system rotation [3]. Moreover, with the turbulence isotropy assumption, the rotation effect on turbulent flow is difficult to describe. However, under the influences of the nonlinear wave induced by the Coriolis force, the flow structure will change [4]. So, a proper correction for EVM is necessary in the RANS simulation of rotating and swirling flow.
To describe the effects of rotation or streamline curvature on a turbulent flow, Bradshaw [5] first proposed the gradient Richardson number . Combined with the lowReynoldsnumber model, Khodak and Hirsch [6] introduced a new form of number which allows including the influence of curvature and rotation on the threedimensional turbulent flow. It improved the prediction accuracy of mean velocity and Reynolds stresses in verification cases. On the basis of Galilean invariance, Spalart and Shur [3] differentiated the property between streamline curvature and system rotation and also introduced a new number. By establishing intermediate variables, this new number was used to correct the turbulence production term of EVM. This correction method was applied to SpalartAllmaras model by Shur et al. [7] and proved to be much more accurate than the original SA model. Moreover, Dufour et al. [8, 9] corrected the model and compared it with the SpalartShur corrected SA model through a compressor case. Improvements were obtained after the rotation and curvature correction for the corrected model. Smirnov and Menter [10] also applied the SpalartShur correction to SST model that is known as the SSTRC model. The computational accuracy was proved to be significantly improved with just a little increase of the time costs. Then, Dhakal and Walters [11] corrected the SST model based on the correction method by York et al. [12]. The new model was compared with the SSTRC model and was proved accurate without compromising stability and efficiency. With the modified meanflow time scale, Hellsten [13] also introduced a new number to correct the equation of SST model. This new corrected model, which is called SSTRCHellsten model, was sensitized on the effects of system rotation and streamline curvature with slight improvement on the numerical behavior. Based on SpalartShur correction method and Hellsten’s time scale, Zhang and Yang [14] corrected the SA model and got reasonable results in a Uturn duct case.
In general, all the correction methods mentioned above were theoretically feasible for swirling turbulent flow. For centrifugal pumps, the internal flow regime varies with the operating condition [15]. Under offdesign conditions, the turbulent flow in the impeller is much more complicated due to the vortex, backflow, and other secondary flow structures. Hence, it is necessary to ensure the simulation accuracy of centrifugal pumps under the influence of system rotation and streamline curvature. However, in the vast majority of RANS simulations of centrifugal pumps, the rotation and curvature effects are not considered when using EVM. In order to evaluate the SpalartShur rotation/curvature correction in the centrifugal pump cases, RANS simulations were conducted with both the original and the corrected SST models. The simulation results were compared with the experimental data by particle image velocimetry and laser Doppler velocimetry [16, 17].
2. Turbulence Modeling
As mentioned above, the SSTRC turbulence model [10] was used in the turbulent flow simulation. Compared with the original SST model [18], in consideration of the turbulence anisotropy, the rotation/curvature correction coefficient was introduced as a multiplier of the turbulence production term in the turbulence kinetic energy equation and specific dissipation rate equation. So, the production term is defined as follows: where is given empirically with specific limiters as follows: where , , and are constants valued as 1.0, 2.0, and 1.0 respectively. The remaining functions are defined as follows: where is the rotation rate of the reference frame and the term represents the Lagrangian derivative of the strain rate tensor. is the strain rate tensor and is the rotation rate tensor by
The in (1) is the Reynolds stress tensor. In the Boussinesq assumption [2], is proportional to the strain rate tensor by where is the turbulent eddy viscosity. By applying the rotation/curvature correction, anisotropic effects were considered in the simulation when solving the turbulence equations.
3. Numerical Simulation
To assess the SpalartShur rotation/curvature correction in RANS simulation of centrifugal pump, the internal flow in the impeller was numerically simulated. Based on the numerical results and available experimental data [16], the flow details were comparatively investigated under both the design and low flow rate offdesign conditions.
3.1. Pump Impeller Model
The scheme of the investigated centrifugal pump impeller is shown in Figure 1. The specificspeed of this impeller is about 26.3 calculated by where is the rotational speed of 725 r/min, is the design flow rate of m^{3}/s, and is the design head of 1.75 m. The geometrical parameters of impeller are shown in Table 1 and illustrated in Figure 1. Two operating conditions including and were simulated and performed in this study.

3.2. Flow Domain Discretization
The flow domain consisted of the impeller only. For a better geometric adaptability, tetrahedral mesh elements were used to discretize the impeller domain. For the usage of wall functions, of the first element outside the walls should be set in the loglayer. So, prism boundary layers were used in the near wall region. Then, a meshsize independence check was conducted to compromise the accuracy and costs in the simulation. By modifying the mesh size, the residuals of head and hydraulic efficiency was ensured to be less than . By modifying the boundary layer height, the values were controlled within the range from 1.53 to 26.97 so that the nearwall region could be solved with wall functions. The final mesh scheme is shown in Table 2 and Figure 2.

3.3. Simulation Settings
In this study, transient numerical simulations were conducted. Three dimensional incompressible NS equations were solved in the simulation process. The fluid medium was set as water at 25 degree centigrade (density kg/m^{3} and dynamic viscosity kg/m·s). The reference frame was set as rotational with the speed of 725 r/min. The reference pressure was 1 Atm. Mass flow inlet was set at the impeller inflow with the velocity normal to the boundary. Pressure outlet was set at the impeller outflow with a value of 0 Pa. Noslip wall condition was given on the solid wall boundaries including hub, shroud and blades.
4. Results and Discussion
In order to fully assess the impacts of SpalartShur correction, as mentioned above, the simulations were proceeded under the design () and offdesign () conditions. The design condition is the most important operating condition of pump. A correct simulation of the flow details is obviously significant. The offdesign condition is also crucial. Undesirable flow structures make the flow hard to predict. For this reason, the improvements of simulation accuracy are necessary. Hence, based on the CFD results, the comparative assessment of the rotation/curvature correction and discussions are given as follows.
4.1. Flow under the Design Condition
Under the design condition (), the velocity field on the spanwise 50% surface was simulated and compared with the LDV data [16] as shown in Figure 3. It can be seen that the flow regime is smooth, stable, and uniform among all the impeller passages.
(a) SST
(b) SSTRC
(c) LDV
The correction coefficient contour on the spanwise 50% surface is shown in Figure 4. It indicates the enhancement or reduction of local turbulence production under the influence of rotation and curvature. As the multiplier of turbulence production term, the value of was almost 1.0 in the vast majority of impeller domain. However, small scale of reduction was detected at the blade leading edge (LE) and in the near suction surface (SS) region.
Figure 5 shows the comparison map of turbulence kinetic energy . In the impeller, fluid separated from the blade surface while flowing around the LE. High region occurred due to the small scale local separation. Also, in the near SS region, high turbulence occurred due to the local unattached flow. After the rotation/curvature correction, the range of high region became smaller than before because of the turbulence reduction. As shown in Figure 5, by comparing the value and pattern, the result by SSTRC model was more consistent with the PIV data [16].
(a) SST
(b) SSTRC
(c) PIV
Additionally, the relative velocity distributions on the spanwise 50% surface are analyzed and shown as follows. The radial component () and tangential component () in three different radius’ positions were plotted, respectively, in Figure 6 (). Considering the flow uniformity among all the passages, the velocity distributions were compared based on just a single passage. Seen from the velocity curves, the velocities magnitude by SSTRC model was just slightly bigger than that by SST model. In consideration of the deviations between the PIV and LDV data [16], both the original and the corrected turbulence models get the similar velocity distribution regularities.
(a) Radius
(b) Radius
(c) Radius
In general, the impact of rotation/curvature correction on the simulation accuracy was not obvious under . Nevertheless, there was no obvious extra time cost when using SSTRC model instead of SST model. So, it would be reasonable to use the SSTRC model in the simulation of centrifugal pump under the design condition.
4.2. Flow under the Low Flow Rate OffDesign Condition
Under the low flow rate offdesign condition (0.25), the velocity field on the spanwise 50% surface was also simulated and compared with the LDV data [16] as shown in Figure 7. With the flow rate decreasing, flow regime in the impeller became undesirable with secondary flow structures. As shown in the velocity field map, fluid did not flow along the direction of blade geometry. Back flow from outlet to inlet occurred in the passage. Lateral secondary flow from blade pressure surface (PS) to suction surface (SS) also occurred. Under the influence of all the disordered flow structures, the flow uniformity among all the passages disappeared. Some passages were blocked by secondary flow, but some other passages were smoother. As shown in the LDV experiment, the blocked passage and unblocked passage occurred alternately in the impeller [16]. This phenomenon was also captured by numerical simulations. In this situation, the streamline curvature under 0.25 became more obvious than that under 1.0.
(a) SST
(b) SSTRC
(c) LDV
The correction details under 0.25 are shown in Figure 8. As plotted in the contour, reductions of turbulence production were detected in the vast majority of impeller domain. But influenced by the flow regime, the rotation/curvature correction was different in each passage. The passages marked “A” and “B” in Figure 8 represented the unblocked and blocked passages, respectively. In passageA, the low region occurred at the whole blade SS. The high region occurred at the PS near LE and the midpassage near trailing edge (TE). In passageB, the low region occurred at SS near LE and midpassage near TE. The high region occurred at SS near TE.
Influenced by the differences of coefficient , the turbulence kinetic energy had also changed after the correction. Figure 9 shows the contour on the spanwise 50% surface under 0.25. As shown in the contours, the turbulence kinetic energy was low in the unblocked “A” passages and was high in the blocked “B” passages. Due to the flow separation at LE and backflow at TE, the high region occurred at LE and TE in “B” passages. On the contrary, because of the wellbehaved flow regime, the value was lower in “A” passages. Compared with the PIV data [16], the location and intensity of were more reasonable after correction. Moreover, the alternately blocked phenomena were more obvious by using the SSTRC model.
(a) SST
(b) SSTRC
(c) PIV
To assess in details, the relative velocity distributions on the spanwise 50% surface are also plotted in Figure 10. Considering the differences of flow regime among all the passages, two adjacent passages (“A” and “B”) were analyzed. The distributions of velocity components changed after correction. Compared with the experimental data [16], it can be seen that the results became more accurate. In particular, flow severely separated at blade LE with a stronger streamline curvature and higher turbulence. So, corrected by reducing the production of turbulence, radial and tangential velocities became more consistent with experiments especially at 0.5.
(a) Radius
(b) Radius
(c) Radius
All in all, the correction impact was obvious under 0.25. So, it is strongly recommended to use the corrected model under the low flow rate offdesign condition.
5. Conclusions
By comparatively assessing the SpalartShur correction in the RANS simulations in centrifugal pump impeller under different operating conditions, conclusions can be drawn as follows.(1)In the RANS simulations, the isotropous description of turbulence model is not perfect enough. The flow in a centrifugal pump impeller is strongly affected by the system rotation and streamline curvature. With the pump rotation, separation flow occurred at blade leading edge. Under different operating conditions, the scale of separation is also different. Particularly under low flow rate offdesign conditions, fluid does not flow along the blade geometry; secondary flow structures become more and more obvious and occurred everywhere in the pump impeller passages.(2)By the supplements of descriptions of turbulence anisotropy, turbulence production term is corrected. Verified by comparing the CFD results with experimental data, improvements are found after correction. Under the design condition, the impact of correction is not obvious but theoretically reasonable. Under the low flow rate offdesign condition, simulation accuracy is significantly improved especially in the strong separation region. Moreover, there is no obvious extra time cost when using the corrected model. Hence, in the RANS simulations of centrifugal pumps, if the flow regime is undesirable with strong secondary flow structures, it will be very appropriate to apply the SpalartShur rotation/curvature correction to the SST model and other eddy viscosity turbulence models.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors would like to acknowledge the financial support given by the National Natural Science Foundation of China (no. 51139007) and National “Twelfth FiveYear” Plan for Science & Technology Support (no. 2012BAD08B03).
References
 J. Boussinesq, Essai sur la Théorie des Eaux Courantes, Imprimerie Nationale, Paris, France, 1877.
 F. G. Schmitt, “About Boussinesq's turbulent viscosity hypothesis: historical remarks and a direct evaluation of its validity,” Comptes Rendus Mécanique, vol. 335, no. 910, pp. 617–627, 2007. View at: Publisher Site  Google Scholar
 P. R. Spalart and M. Shur, “On the sensitization of turbulence models to rotation and curvature,” Aerospace Science and Technology, vol. 1, no. 5, pp. 297–302, 1997. View at: Publisher Site  Google Scholar
 Y.N. Huang and H.Y. Ma, “Extended intrinsic mean spin tensor for turbulence modelling in noninertial frame of reference,” Applied Mathematics and Mechanics, vol. 29, no. 11, pp. 1463–1475, 2008. View at: Publisher Site  Google Scholar
 P. Bradshaw, “The analogy between streamline curvature and buoyancy in turbulent shear flow,” Journal of Fluid Mechanics, vol. 36, no. 1, pp. 177–191, 1969. View at: Publisher Site  Google Scholar
 A. Khodak and C. Hirsch, “Secondorder nonlinear kε models with explicit effect of curvature and rotation,” in Proceedings of the 3rd ECCOMAS Computational Fluid Dynamics Conference, pp. 690–696, 1996. View at: Google Scholar
 M. L. Shur, M. K. Strelets, A. K. Travin, and P. R. Spalart, “Turbulence modeling in rotating and curved channels: assessing the SpalartShur correction,” AIAA Journal, vol. 38, no. 5, pp. 784–792, 2000. View at: Publisher Site  Google Scholar
 G. Dufour, J.B. Cazalbou, X. Carbonneau, and P. Chassaing, “Assessing rotation/curvature corrections to eddyviscosity models in the calculations of centrifugalcompressor flows,” Journal of Fluids Engineering, Transactions of the ASME, vol. 130, no. 9, 2008. View at: Publisher Site  Google Scholar
 J.B. Cazalbou, P. Chassaing, G. Dufour, and X. Carbonneau, “Twoequation modeling of turbulent rotating flows,” Physics of Fluids, vol. 17, no. 5, Article ID 055110, 2005. View at: Publisher Site  Google Scholar
 P. E. Smirnov and F. R. Menter, “Sensitization of the SST turbulence model to rotation and curvature by applying the SpalartShur correction term,” Journal of Turbomachinery, vol. 131, no. 4, Article ID 041010, 2009. View at: Publisher Site  Google Scholar
 T. P. Dhakal and D. K. Walters, “A Threeequation variant of the SST kw model sensitized to rotation and curvature effects,” Journal of Fluids Engineering, vol. 133, no. 11, Article ID 111201, 2011. View at: Publisher Site  Google Scholar
 W. D. York, D. K. Walters, and J. H. Leylek, “A simple and robust linear eddyviscosity formulation for curved and rotating flows,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 19, no. 6, pp. 745–776, 2009. View at: Publisher Site  Google Scholar
 A. Hellsten, “Some improvements in Menter's SST kω turbulence model,” in Proceedings of the 29th AIAA Fluid Dynamics Conference, AIAA 982554, Albuquerque, NM, USA, June 1998. View at: Google Scholar
 Q. Zhang and Y. Yang, “A new simpler rotation/curvature correction method for SpalartAllmaras turbulence model,” Chinese Journal of Aeronautics, vol. 26, no. 2, pp. 326–333, 2013. View at: Publisher Site  Google Scholar
 C. E. Brennen, Hydrodynamics of Pumps, Cambridge University Press, Cambridge, UK, 2011.
 N. Pedersen, P. S. Larsen, and C. B. Jacobsen, “Flow in a centrifugal pump impeller at design and offdesign conditions—part I: particle image velocimetry (PIV) and laser Doppler velocimetry (LDV) measurements,” Journal of Fluids Engineering, vol. 125, no. 1, pp. 61–72, 2003. View at: Publisher Site  Google Scholar
 R. K. Byskov, C. B. Jacobsen, and N. Pedersen, “Flow in a centrifugal pump impeller at design and offdesign conditions—part II: large eddy simulations,” Journal of Fluids Engineering, Transactions of the ASME, vol. 125, no. 1, pp. 73–83, 2003. View at: Publisher Site  Google Scholar
 F. R. Menter, M. Kuntz, and R. Langtry, “Ten years of industrial experience with the SST turbulence model,” Turbulence, Heat and Mass Transfer, vol. 4, pp. 625–632, 2003. View at: Google Scholar
Copyright
Copyright © 2014 Ran Tao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.