AbstractThe paper presents an efficient finite volume method for unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements for both single‐ and two‐phase flows. Mathematical model used in computations is based on the ensemble averaged conservation equations. These equations are solved for each phase and in case of single‐phase flow reduce to the transient Reynolds‐averaged Navier–Stokes (TRANS) equations. Transient flow induced by rotating impellers is thus resolved in time. The use of unstructured grids allows an easy and flexible meshing for the entire flow domain. Polyhedral cell volumes are created on the arbitrary mesh interface placed between rotating and static parts. Cells within the rotating parts move each time step and the new faces are created on the arbitrary interfaces only, while the rest of the domain remain ‘topologically’ unchanged. Implicit discretization scheme allows a wide range of time‐step sizes, which further reduce the computational effort. Special attention is given to the interpolation practices used for the reconstruction of the face quantities. Mass fluxes are recalculated at the beginning of each time step by using an interpolation scheme, which enhances the coupling between the pressure and velocity fields. The model has been implemented into the commercially available CFD code AVL SWIFT (AVL AST, SWIFT Manual 3.1, AVL List GmbH, Graz, Austria, 2002). Single‐phase flow in a mixing vessel stirred by a six‐bladed Rushton‐type turbine and two‐phase flow in aerated stirred vessel with the four‐blade Rushton impeller are simulated. The results are compared with the available experimental data, and good agreement is observed. The proposed algorithm is proved to be both stable and accurate for single‐phase as well as for the two‐phase flows calculations. Copyright 2004 John Wiley & Sons, Ltd.