Finite volume schemes for the biharmonic problem on general meshes

Abstract

During the development of a convergence theory for Nicolaides’ extension of the classical MAC scheme for the incompressible Navier-Stokes equations to unstructured triangle meshes, it became clear that a convergence theory for a new kind of finite volume discretizations for the biharmonic problem would be a very useful tool in the convergence analysis of the generalized MAC scheme. Therefore, we present and analyze new finite volume schemes for the approximation of a biharmonic problem with Dirichlet boundary conditions on grids which satisfy an orthogonality condition. We prove that a piecewise constant approximate solution of the biharmonic problem converges in L 2 ( Ω ) L^2(\Omega ) to the exact solution. Similar results are shown for the discrete approximate of the gradient and the discrete approximate of the Laplacian of the exact solution. Error estimates are also derived. This part of the paper is a first, significant step towards a convergence theory of Nicolaides’ extension of the classical MAC scheme. Further, we show that finite volume discretizations for the biharmonic problem can also be defined on very general, nonconforming meshes, such that the same convergence results hold. The possibility to construct a converging lowest order finite volume method for the H 2 H^2 -regular biharmonic problem on general meshes seems to be an interesting result for itself and clarifies the necessary ingredients for converging discretizations of the biharmonic problem. All these results are confirmed by numerical results.