We present and analyse in this paper a novel colocated Finite Volume scheme for the solution of the Stokes problem.
It has been developed following two main ideas.
On one hand, the discretization of the pressure gradient term is built as the discrete transposed of the velocity divergence term, the latter being evaluated using a natural finite volume approximation; this leads to a non-standard interpolation formula for the expression of the pressure on the edges of the control volumes.
On the other hand, the scheme is stabilized using a finite volume analogue to the Brezzi-Pitkäranta technique.
We prove that, under usual regularity assumptions for the solution (each component of the velocity in H2(∞) and pressure in H1(∞)), the scheme is first order convergent in the usual finite volume discrete H1 norm and the L2 norm for respectively the velocity and the pressure, provided, in particular, that the approximation of the mass balance flux is of second order.
With the above-mentioned interpolation formulae, this latter condition is satisfied only for particular meshes: acute angles triangulations or rectangular structured discretizations in two dimensions, and rectangular parallelepipedic structured discretizations in three dimensions.
Numerical experiments confirm this analysis and show, in addition, a second order convergence for the velocity in a discrete L2 norm.