The purpose of this work is to introduce and validate a new staggered control volume method for the simulation of 2Daaxisymmetric incompressible ¯ows.
The present study introduces a numerical procedure for solving the Navier±Stokes equations using the primitive variable formulation. The proposed method is an extension of the staggered grid methodology to unstructured triangular meshes for a control volume approach which features ease of handling of irregularly shaped domains. Two alternative elements are studied: transported scalars are stored either at the sides of an element or at its vertices, while the pressure is always stored at the centre of an element. Two interpolation functions were investigated for the integration of the momentum equations: a skewed mass-weighted upwind function and a ¯ow-oriented exponential shape function. The momentum equations are solved over the covolume of a side or of a vertex and the pressure±velocity coupling makes use of a localized linear reconstruction of the discontinuous pressure ®eld surrounding an element in order to obtain the pressure gradient terms. The pressure equation is obtained through a discretization of the continuity equation which uses the triangular element itself as the control volume.
The method is applied to the simulation of the following test cases: backward-facing step ¯ow, ¯ow over a two-dimensional obstacle and ¯ow in a pipe with sudden contraction of cross-sectional area. All numerical investigations are compared with experimental data from the literature. A grid convergence and error analysis study is also carried out for ¯ow in a driven cavity.
Results compared favourably with experimental data and so the new control volume scheme is deemed well suited for the prediction of incompressible ¯ows in complex geometries.