- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Table of Contents

Advances in Mechanical Engineering

Volume 2013 (2013), Article ID 710769, 11 pages

http://dx.doi.org/10.1155/2013/710769

## Instability Analysis of a Model Pump-Turbine with MGV Based on Nonlinear Partially Averaged Navier-Stokes Methods

^{1}Beijing Institute of Control Engineering, China Aerospace and Technology Corporation, Beijing 100190, China^{2}State Key Laboratory of Hydro Science and Engineering, Tsinghua University, Beijing 100084, China^{3}Institute of Chemical Process Machinery, Zhejiang University, Hangzhou 310027, China

Received 2 July 2013; Revised 15 September 2013; Accepted 16 September 2013

Academic Editor: Jiin Yuh Jang

Copyright © 2013 Jintao Liu 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.

#### Abstract

Pump-turbines were always running at partial condition with the power grid changing. Flow separations and stall phenomena were obvious in the pump-turbine. Most of the RANS turbulence models solved the shear stress by linear difference scheme and isotropic models, so they could not capture all kinds of vortexes in the pump-turbine well. At present, partially-averaged Navier-Stokes (PANS) model has been found to be better than LES in simulating flow regions especially those with less discretized grid. In this paper, a new nonlinear PANS turbulence model was proposed, which was modified from RNG - turbulence model and the shear stresses were solved by Ehrhardt’s nonlinear methods. The nonlinear PANS model was used to study the instability of “S” region of a model pump-turbine with misaligned guide vanes (MGV). The opening of preopened guide vanes had great influence on the “S” characteristics. The optimal relative opening of the preopened guide vanes was 50% for the improvement of the “S” characteristics. Pressure fluctuations in the vaneless space were analyzed. It is found that the dominant frequency at the vaneless space was twice the blade passing frequency, while the second dominant frequency decreased as the preopening increased.

#### 1. Introduction

Pump-turbine, which was one of the key components in a pumped storage power station, played an important role in the optimization of power network. The starting period and stopping period of a pump-turbine ran at a high frequency, leading to a lot of instability problems due to the existence of “S” characteristics. When a pump-turbine ran in the “S” region, pressure fluctuation was very dangerous, even leading to some serious accidents. At present, misaligned guide vanes (MGV) were used in some pumped storage power stations to improve the “S” characteristics, but the influence of the relative opening of preopened guide vanes on the “S” characteristics was unclear.

Pump-turbines were always running at partial condition with the power grid changing, leading to flow separations and stall phenomena in the pump-turbine. The study of instability of a pump-turbine should overcome the difficulties in the simulation of strong vortex and separation flow. Large-scale coherent structures play a crucial role in the internal instabilities. For predicting the fluctuation of the dominant scale of vortexes, the traditional Reynolds-averaged Navier-Stokes (RANS) method suffers from inherent physical limitations [1, 2], which could not capture the flow with all kinds of scales. Large eddy simulation (LES) may not be computationally viable due to the large computational cost. Instead, the mixed computational methods, which combined the desirable aspects of RANS and LES, were used in engineering computation [3, 4]. Speziale [5] proposed a new turbulence model that combines the advantages of RANS method with those of LES [6]. Girimaji et al. [7] developed a bridging method inspired by the modeling paradigm proposed by Speziale [5]. The method was given the name partially averaged Navier-Stokes (PANS) model and was purported for any filter width [8].

At present, partially averaged Navier-Stokes (PANS) model has been found to be better than that of the LES in flow regions where simulations suffered from poor near-wall resolution [9]. Ji et al. [10, 11] studied the unsteady cavitating turbulence flow around a highly skewed model marine propeller based on the PANS method. Ma et al. [12] and Huang and Wang [13] used PANS model to investigate the cavitating flow around a hydrofoil. The PANS models mentioned above were all modified from standard - turbulence model, while the standard - turbulence model was poor in the simulation of strong swirling flows [14]. Most of the RANS turbulence models solved the shear stress by linear difference scheme and they were isotropic models [15], so they could not capture all kinds of vortexes in the pump-turbine well.

In this paper, a new nonlinear PANS turbulence model was proposed. The results based on the new nonlinear PANS turbulence model were compared with experimental results. Then the instability of a pump-turbine was studied by the nonlinear PANS model.

#### 2. Nonlinear PANS Model

For incompressible flow, the continuity equation and Reynolds averaged Navier-Stokes equations are as follows:

is partitioned into resolved and unresolved parts in the instantaneous velocity field, using an arbitrary homogeneous filter: where is the resolved velocity field and is the unresolved field. It is used by (3) instead of the resolved field The additional nonlinear term , which is the generalized central second moment, is defined as:

In RNG - turbulence mode, the equations for turbulent kinetic energy () and turbulent dissipation () were shown as follows: where where , , , , , , and denotes the production terms of turbulence kinetic energy.

For PANS methods,

Filter control parameters are and . These two parameters are the partial average-quantification. The natural to quantify the extent of PANS filtering using and . Considering large scales turbulence flows contain most of the kinetic energy, while the smallest scales flows produce much of the dissipation. The parameters dictated that . The smaller the , the finer the filter: represents RANS and indicates DNS. Unity value of implies that the RANS and PANS unresolved small scales are identical. Smaller values of would require the resolution of dissipative scales of motion.

Equating the source terms, in RNG - turbulence mode has a relationship with in PANS model shown in (7):

Then we can obtain PANS model by the modification of RNG - turbulence model, and the equations are Thus,

Considering the nonlinear turbulence flow in the pump-turbine, the shear stress was solved by nonlinear turbulence model which was proposed by Liu et al. [15]: where , , , , , , , is turbulence time scale, and is the turbulence velocity scale.

#### 3. Pump-Turbine Geometry

Parameters of the model pump-turbine are shown in Table 1. denotes the runner inlet diameter in pump mode; , , and are the numbers of stay vanes, guide vanes, and runner blades, respectively; denotes the rated head; denotes the rotational speed of the runner; denotes the rated discharge.

The pump-turbine’s structure is shown in Figure 1. The structure of MGV is shown in Figure 2.

#### 4. Simulation Conditions and Verification of Grid Independence

The model’s grids, which were composed of an unstructured hexahedron and tetrahedron, were developed using ICEM, which is a commercial software package used for CFD discretization. Hexahedral grids were used for the runner and draft tube, and mixed grids were used for the other components. The mesh of the runner is shown in Figure 3.

The commercial CFD code FLUENT was used to perform the simulations. The boundary conditions of the pump-turbine were different in the reverse pump mode and in the turbine mode. The velocity at the spiral casing inlet was given for the turbine mode, whereas the pressure was given for the reverse pump mode. The runner’s hydraulic region was set as a moving mesh model. The SIMPLEC algorithm was used to enforce mass conservation. For the quad/hex grids and complex flows in turbine mode, a second-order centered difference was used for the pressure interpolation to obtain better results. Taking the higher pressure gradient into consideration for reversed pump mode, PRESTO! (FLUENT Inc., 2005) discretization was used in the pressure interpolation. For present unsteady flow calculation, the time step was 0.00019 s. The nonlinear PANS model was used to calculate the “S” characteristics. For all calculations, simulations were run until convergence, which was determined by a reduction in the residual error to less than 0.0001.

The nonlinear PANS method was performed by the commercial software FLUENT by user defined function (UDF). Couple method was used to correct the pressure during the calculation. Second-order upwind scheme was used in resolving the Navies-Stocks equations. The rate of unsolved unresolved-to-total kinetic energy and the rate of unsolved unresolved-to-total kinetic energy dissipation were used in the calculation.

To verify the calculation accuracy in terms of the grid size, several grid scales were selected. The cell number varied from 4 million to 10 million. The efficiency of the pump-turbine with different cell numbers when the pump-turbine was run in turbine mode at a no-load opening is shown in Figure 4. The pressure distribution on the pressure side of the runner blade is shown in Figure 5. The pressure distribution performed by 8.2 million is almost the same as the result performed by 9.7 million. It can be seen that the degree of computational precision is satisfactory when the element number is more than 7 million. A mesh with about 9 million cells in total was chosen for the simulations. The number of nodes and elements in each part is shown in Table 2.

#### 5. Results and Discussion

##### 5.1. “S” Characteristics with MGV

“S” characteristics of the pump-turbine with MGV at no-load opening are shown in Figure 6. The simulation results of the “S” characteristics with synchronized guide vane (SyV) based on PANS model agree well with experimental data. The result of “S” curve based on RNG - model has large error compared with experimental result. The RNG - model could not capture the detailed flow phenomena in the pump-turbine with small opening of guide vanes, causing the error of external characteristics. The capability comparisons of RNG - model and PANS model are shown in the Appendix. MGV can improve the “S” curve, without changing the unit speed and unit flow rate at the runaway point. The flow rate of the pump-turbine with MGV is larger than the result of SyV in turbine mode, while it is smaller when it runs at reverse pump mode.

The head and flow rate curve is shown in Figure 7 and the torque of the runner and the flow rate is shown in Figure 8. MGV increases the discharge capacity of the pump-turbine within the same head. Flow rates of the two cases are very close at the runaway point. When the pump-turbine runs at reverse pump mode, the flow rate of the pump-turbine with SyV is larger than the unit with MGV. Downward displacement of the torque curve can be found in Figure 8 by the use of MGV.

##### 5.2. Flow in the Pump-Turbine with MGV

During the starting period, the “S” region with determines the instability of the pump-turbine, and the relative opening of preopened guide vanes has a relationship with the “S” characteristics, so the “S” characteristics with in different relative openings are calculated. The results of the “S” characteristics with different relative openings are shown in Figure 9. The openings of preopened guide vanes are 10°, 15°, and 21°, and their corresponding relative openings are 33%, 50%, and 70%, respectively. The existence of the “S” characteristic can be seen when the relative opening is less than 50%, but the increase of relative opening of guide vane improves the “S” region. The “S” characteristic disappears when the relative opening is 50%, and that is to say there is a one-to-one relationship between the and . The increase of relative opening of preopened guide vane will not make further improvement on the “S” region when the relative opening is larger than 50%. The increase of relative opening will cause flow instability in the pump-turbine, even leading to serious pressure fluctuation. The chosen of optimal relative opening should consider both the “S” characteristics and the flow instability.

##### 5.3. Flow in the Pump-Turbine with MGV

The runaway point is the point for starting load during the start-up process; it is very important for the instability of a pump-turbine. Streamline in the runner on blade-to-blade surface at no-load condition is shown in Figure 10. Flow in the runner with SyV has periodic distribution. MGV changes the flow in the runner. Different flow rates in the 20 passages of guide vane lead to reverse flow in the runner. The flow rate in the passage of the runner against the preopened guide vanes is larger than others, which may be the reason for the change of the performance of the pump-turbine.

Flow in the casing and the region of guide vanes is shown in Figure 11. No reverse flow can be seen with SyV. Because the relative opening of guide vanes is small, the resistance is large enough for the fluid flowing through the space between two adjacent guide vanes so that a water ring swirling at a high speed is formed in the vaneless space. The use of MGV destroys the water ring, as can be found in Figure 10(b).

Vortex rope in the draft tube at the runaway point is shown in Figure 12. The vortex rope is calculated by the pressure contour. Vortex rope in the draft tube has large volume with SyV compared with MGV. The use of MGV changes the shape of vortex rope into foam-like structure.

##### 5.4. Pressure Fluctuation in the Vaneless Space

Pressure fluctuations in the vaneless space at runaway point are shown in Figure 13. The fluctuation of the pressure in the pump-turbine with SyV is consistent, and the pressure fluctuation only contains high frequency component. The use of MGV increases the amplitude of the pressure fluctuation. With the increase of the relative opening of preopened guide vane, the low frequency of pressure fluctuation is more and more obvious. Low frequency is very dangerous to the operation of the pump-turbine.

The pressure fluctuations at frequency domain are calculated by fast Fourier transform (FFT), which are shown in Figure 14. The spectrum characteristics of the pressure fluctuation are shown in Table 3. Results of pressure fluctuation in the vaneless space agree well with the experimental data. The amplitude of pressure fluctuation increases as the opening of preopened guide vane increases. The amplitude is twice the result with SyV when the pre-opening is 15°. At the runaway point, the dominant frequency (DF) of pressure fluctuation is the blade passing frequency. The second dominant frequency (SDF) of pressure fluctuation is two times the blade passing frequency when the preopening is small. The increase of the preopening makes the second dominant frequency lower and lower.

#### 6. Conclusions

A PANS model was proposed based on RNG - turbulence model. The shear stress was resolved by nonlinear model. The nonlinear PANS model was used to simulate the flow field in a model pump-turbine and was proved to be accurate in the simulation of “S” characteristics. The “S” characteristics of a pump-turbine with MGV were investigated and compared with the results of SyV. The MGV improved the “S” characteristics by the changing of the flow field in the pump-turbine. It destroyed the water ring in the vaneless space and changed the periodic flow in the runner. The use of MGV aggravated the pressure fluctuation in the vaneless space. The amplitude of pressure fluctuation increased with the preopening increase. The dominant frequency of pressure fluctuation in the vaneless space was the blade passing frequency, but the second dominant frequency gradually decreased as the preopening increased.

#### Appendix

To check the capability of the nonlinear PANS model for the prediction of complex flow in turbo machineries, a closed-loop test rig of Pedersen et al. [16] shown in Figure 15 was investigated. Fluid was cycling in the closed-loop test rig. Pedersen et al. [16] investigated the flow field of the centrifugal pump by particle image velocimetry (PIV) method. The velocity distributions on the measuring surface were calculated and compared with experimental results.

When the pump runs at small flow rate condition , velocities at different diameters in the impeller on the measuring surface at the small flow rate condition are shown in Figure 16. is the rated flow rate. is the flow rate of the pump. Streamlines on the measuring surface in the impeller performed by three models are shown in Figure 17. Calculation result of the nonlinear PASN model shows that flow in the impeller has a periodic distribution. No periodic phenomenon can be seen from the result of RNG - turbulence model. The velocity distributions performed by the nonlinear PANS model agree well with the experimental result. The nonlinear PANS model can capture the large-scale vortex in the impeller accurately. The velocity at the impeller inlet is larger than the result of the nonlinear PANS model. The result of RNG -03B5 turbulence model varies greatly from the experimental result. The closed-loop test rig has a periodic distribution in the hydraulic region. The flow in the cycling system should have periodic distribution, so that the result of the nonlinear PANS model is more realistic.

#### 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 thank projects 51076077, 51176168, and 51306018 supported by National Natural Science Foundation of China and project supported by National Science and Technology Ministry (ID: 2011BAF03B01).

#### References

- M. R. Khorrami, B. A. Singer, and M. E. Berkman, “Time-accurate simulations and acoustic analysis of slat free shear layer,”
*AIAA Journal*, vol. 40, no. 7, pp. 1284–1291, 2002. View at Scopus - M. R. Khorrami, B. Singer, and D. Lockard, “Time-accurate simulations and acoustic analysis of slat free shear layer: part 2,” in
*Proceedings of the AIAA/CEAS Aeroacoustics Conference and Exhibit*, AIAA 2002-2579, Breckenridge, Colo, USA. - P. R. Spalart, “Trends in turbulence treatments,”
*AIAA*2000-2306, Fluids 2000, Denver, Colo, USA, 2000. - P. Batten, U. Goldberg, and S. Chakravarthy, “LNS-an approach towards embedded LES,”
*AIAA*2002-0427, 2002. - C. G. Speziale, “Computing non-equilibrium flows with time-dependent RANS and VLES,” in
*Proceedings of the 15th International Conference on Numerical Methods in Fluid Dynamics (ICNMFD ’96)*, Monterey, Calif, USA, 1996. - W. T. Su, X. B. Li, and F. C. Li, “Comparisons of LES and RANS computations with PIV experiments on a cylindrical cavity flow,”
*Advances in Mechanical Engineering*, vol. 2013, Article ID 592940, 10 pages, 2013. View at Publisher · View at Google Scholar - S. S. Girimaji, E. Jeong, and R. Srinivasa, “PANS method for turbulence: fixed point analysis and comparison with URANS,”
*Journal of Applied Mechanics*, vol. 73, pp. 422–429, 2006. - S. S. Girimaji, “Partially-averaged Navier-Stokes model for turbulence: a reynolds-averaged Navier-Stokes to direct numerical simulation bridging method,”
*Journal of Applied Mechanics*, vol. 73, no. 3, pp. 413–421, 2006. View at Publisher · View at Google Scholar · View at Scopus - S. Lakshmipathy,
*Partially averaged Navier-Stokes method for turbulence closures: characterization of fluctuations and extension to wall bounded flows [Ph.D. thesis]*, Texas A&M University, 2009. - B. Ji, X. W. Luo, Y. L. Wu, X. X. Peng, and H. Y. Xu, “Partially-Averaged Navier-Stokes method with modified
*k*-*ε*model for cavitating flow around a marine propeller in a non-uniform wake,”*International Journal of Heat and Mass Transfer*, vol. 55, pp. 6582–6588, 2012. - B. Ji, X. W. Luo, Y. L. Wu, and H. Y. Xu, “Unsteady cavitating flow around a hydrofoil simulated using partially-averaged Navier-Stokes model,”
*Chinese Physics Letters*, vol. 29, no. 7, Article ID 076401, 2012. - J. M. Ma, F. J. Wang, X. Yu, and Z. Q. Liu, “Partially averaged Navier-Stokes method for time-dependent turbulent cavitating flows,”
*Journal of Hydrodynamics*, vol. 23, no. 4, pp. 466–475, 2011. - B. Huang and G.-Y. Wang, “Partially averaged Navier-Stokes method for time-dependent turbulent cavitating flows,”
*Journal of Hydrodynamics*, vol. 23, no. 1, pp. 26–33, 2011. View at Publisher · View at Google Scholar · View at Scopus - V. Yakhot, S. A. Orszag, S. Thangam, T. B. Gatski, and C. G. Speziale, “Development of turbulence models for shear flows by a double expansion technique,”
*Physics of Fluids A*, vol. 4, no. 7, pp. 1510–1520, 1992. View at Scopus - J. T. Liu, S. H. Liu, Y. L. Wu, et al., “Numerical investigation of the hump characteristic of a pump-turbine based on an improved cavitation model,”
*Computers and Fluids*, vol. 68, pp. 105–111, 2012. - N. Pedersen, P. S. Larsen, and C. B. Jacobsen, “Flow in a centrifugal pump impeller at design and off-design 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 · View at Google Scholar · View at Scopus