Pump performance prediction using potential flow method



Introduction

The design of hydraulic turbomachines, such as pumps and water turbines, has reached the stage were improvements can only be achieved through a detailed understanding of the internal flow. The prediction of the flow in such equipment is very complicated due to the three-dimensionality of the flow and the highly curved passages in rotating impellers. Furthermore, turbomachines show unsteady flow behaviour, especially under off-design conditions, as a result of the interaction between impeller and volute or stator. Considering these complexities, computer simulations of the flow are becoming increasingly important. Over the past few years, there is a tendency towards the development of numerical methods based on the Navier-Stokes equations, in order to account for viscous effects like wakes, boundary layers and separation bubbles. An open problem in such computations is the choice of an appropriate turbulence model, since the standard turbulence models appear not to be adequate in rotating systems. Furthermore, the computer time needed to compute the flow through a single impeller channel is enormeous, even when using supercomputers. As part of a design tool, these methods are of limited suitability.

Assumptions

Fortunately, the physical model can be simplified considerably without losing its overall validity. As a first assumption the fluid can be considered incompressible. Secondly, the bulk of the fluid can be regarded inviscid, as viscous forces are negligible when compared to inertia forces. Viscous effects are restricted to boundary layers and wakes behind the blades. The third assumption is that the flow enters the impeller free of vorticity. For an inviscid bulk flow this means that the flow remains irrotational. Combining these assumptions leads to the the incompressible potential flow. Many investigators have presented numerical methods for computing the potential flow inside a rotor channel in two or three dimensions, and in complete rotor-stator configurations in two dimensions. In order to impose the Kutta condition of smooth flow at the blades' trailing edges they either superimpose a number of subpotentials (reflecting unit flowrate, unit rotation and line vortices of unit strength, shed from points at the trailing edges) or determine the (varying) circulations iteratively.

Multi-block method

In order to reduce computing time and memory requirements as much as possible, an entirely new potential flow solver has been developed. It is capable of computing the unsteady three-dimensional potential flow inside a rotor-stator configuration, taking into account the varying circulation along the blades' span. The method is based on the Finite Element Method. Two techniques are responsible for its high efficiency: the method of substructuring (multi-block method) and the implicit imposition of Kutta conditions at the blades' trailing edges.

Without going to much into detail, the method is as follows. The flow region of interest is divided into subdomains or blocks, all having a topologically cubic shape (see the above figure). Blocks in the rotor of the pump are separated from blocks in the stator by a cylindrical or comical surface. Thus the rotation of the impeller with respect to the stator can be simulated by sliding the rotorblocks along this surface. A structured grid is generated for each of the blocks (figure below). The solution method consists of two steps. In the first one, internal degrees of freedom will be eliminated from the system of equations, for all blocks separately. The remaining degrees of freedom at internal block boundaries are then coupled in the second step. Kutta equations which prescribe flow conditions at trailing edges are added. After solving the resulting system of equations, the previously eliminated degrees of freedom are finally obtained.





PumpFlow

The above described flow solver has been integrated into the information system PumpFlow. This graphic computer system offers a variety of tools, necessary to analyse the performance of pumps or parts of pumps. The information system PumpFlow is being developed to offer the pump developer a powerful and yet easy-to-use tool for analyzing pump performance. It consists of the following sub-systems: These subsystems are embedded into a graphic user interface, which takes care of all necessary file manipulation.

The features of this system include:





Examples

The information system PumpSystem consists of 3 subsystems: PumpGrid, PumpFlow, and PumpPost. Below are a number of frames available as examples of what the different subsystems look like.

PumpGrid
The subsystem PumpGrid takes care of the design of the pump impeller and housing, and the subsequent generation of computational meshes for PumpFlow.

Frame 1 Camber surfaces of impeller blades, showing construction lines and points, in 3D view.
Frame 2 Camber construction lines and points in 3D view, showing values of blade angles in color.
Frame 3 Camber surfaces of impeller blades in 3D view, showing values of blade angles in color.
Frame 4 Hidden surface view of camber surfaces of impeller blades, showing values of blade angles in color.
Frame 5 Shaded view of camber surfaces of impeller blades and hub surface.
Frame 6 Plane view of impeller blade camber surface.
Frame 7 Meridional view of impeller blade camber surface.
Frame 8 2D view of blade angle along meridional coordinate of blade construction lines, radial coordinate displayed in color.
Frame 9 Contour lines of blade angle values in plane view of impeller blade.
Frame 10 Changing the shape of the camber surface by editing construction lines and points.
Frame 11 Adding a construction line or a construction point to a blade camber surface.
Frame 12 Reshaping a construction line of a blade camber surface.
Frame 13 Reshaping the camber surface.


PumpFlow
Frame 1 Selecting a project
Frame 2 Selecting meshblock for topology information
Frame 3 Inspecting topology of a blade surface
Frame 4 Modifying topology of a blade surface
Frame 5 Inspecting boundary condition for surface
Frame 6 Modifying boundary condition for surface
Frame 7 Block connections
Frame 8 Inspecting block connection
Frame 9 Modifying block connection
Frame 10 Inspecting sliding connection between rotor and stator
Frame 11 Modifying sliding connection between rotor and stator
Frame 12 Inspecting wake connection
Frame 13 Setting slitsurface
Frame 14 View object from different positions
Frame 15 Visualise pressure distribution and surface grid
Frame 16 Set view options
Frame 17 Shrink blocks to see block boundaries
Frame 18 View and print shaded surface
Frame 19 Set color palette
Frame 20 Inspect vector quantity




PUBLICATIONS

1. " Experimental investigation of hydrodynamic forces due to non-uniform entrance flow to a mixed-flow pump ," Nyirenda, P.J., van Esch, B.P.M., 2003, Conference on Modelling Fluid Flow (CMFF'03), Budapest, Hungary, pp. 984-990.
2. " Calculation of radial forces due to non-uniform entrance flow in a mixed-flow waterjet pump ," Bulten, N., van Esch, B.P.M., 2003, Conference on Modelling Fluid Flow (CMFF'03), Budapest, Hungary, pp. 977-983.
3. " Experimental and theoretical study of the flow in the volute of a low specific-speed pump ," Kelder, J.D.H., Dijkers, R.J.H., van Esch, B.P.M., Kruyt, N.P., 2001, Fluid Dynamics Research, vol. 28, pp. 267-280.
4. " Hydraulic performance of a mixed-flow pump: unsteady inviscid computations and loss models ," van Esch, B.P.M., Kruyt, N.P., 2001, ASME Journal of Fluids Engineering, vol. 123, pp. 257-264.
5. "A superelement-based method for computing unsteady three-dimensional potential flows in hydraulic turbomachines," N.P. Kruyt, B.P.M. van Esch, J.B. Jonker, 1999, Communications in Numerical Methods in Engineering, 15, issue 6, pp. 381-397. samenvatting & complete tekst
6. "Simulation of three-dimensional unsteady flow in hydraulic pumps," B.P.M. van Esch, 1997, PhD thesis, University of Twente, The Netherlands. pdf-format
7. "An inviscid-viscous coupling method for computing flows in entire pump configurations, " B.P.M. van Esch, N.P. Kruyt, J.B. Jonker, 1997, Proc. third international symposium on pumping machinery, ASME FED Summer Meeting, Vancouver, Canada. pdf-format
8. "A Tool for the Analysis of Unsteady Potential Flows in Centrifugal and Mixed-flow Pumps", N.P. Kruyt, B.P.M. van Esch, J.B. Jonker, 1996, Pumpentagung, 30 Sept - 02 Oct. 1996, Karlsruhe, Germany. abstract
9. "Analysis of the flow in a centrifugal pump using a multi-block finite element method for computing three-dimensional potential flows," B.P.M. van Esch, N.P. Kruyt, 1995, ERCOFTAC Seminar and Workshop i on 3D Turbomachinery Flow Prediction III, January 1995, Val d'Isere, France abstract
10. "An efficient method for computing three-dimensional potential flows in hydraulic turbomachines," B.P.M. van Esch, N.P. Kruyt, J.B. Jonker, 1995, Finite Elements in Fluids - New Trends and Applications, October 1995, Venice, Italy. abstract
11. "Analysis of three-dimensional potential flows in centrifugal and mixed-flow pumps by a finite element method," N.P. Kruyt, B.P.M. van Esch, J.B. Jonker, 1995, 10th Conference on Fluid Machinery, September 1995, Budapest, Hungary. abstract

RELATED PUBLICATIONS

10. "A parametric study of the cavitation inception behavior of a mixed-flow impeller using a three-dimensional potential flow model," M.J. van Os, J.G.H. Op de Woerd, J.B. Jonker, 1997, Proc. third international symposium on pumping machinery, ASME FED Summer Meeting, Vancouver, Canada. pdf-format
11. "The computation of fluid-induced forces on centrifugal impellers rotating and whirling in a volute casing," R.G.K.M. Aarts, J.B. Jonker, 1997, Proc. third international symposium on pumping machinery, ASME FED Summer Meeting, Vancouver, Canada. pdf-format


To view pdf-files, download the latest version of Acrobat Reader from Adobe

Download Acroread



For comments and questions regarding this page email to b.p.m.vanesch@wtb.tue.nl

Last update: June 22, 1999