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.
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.
|
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.
|
|
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
|
|
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
|