Abstract. For Laplace's equation, we discuss whether it is possible to construct a linear positive finite volume (FV) scheme on arbitrary unstructured grids. Dealing with the arbitrary grids, we state a control volume which guarantees a positive FV scheme with linear reconstruction of the solution. The control volume is defined by a property of the analytical solution to the equation and does not depend on the grid geometry. For those problems where the choice of the control volume is prescribed a priori, we demonstrate how to improve positivity of the linear FV scheme by using corrected reconstruction stencils. The difficulties arising when grids with no geometric restrictions are used for the discretization are discussed. Numerical examples illustrating the developed approach to the stencil correction are given.Key words. Laplace's equation, finite volume scheme, positivity, stencil correction
AMS subject classifications. 65N12, 74S10, 35J05PII. S1064827500368925
Introduction.A discrete Laplace operator is often considered to be a good model for investigating a discretization of partial differential equations which contain diffusion operator ∇·(D∇). Two important examples are given by convection-diffusion equations and Navier-Stokes equations with possible applications that include the problems of fluid dynamics, chemical engineering, and environmental pollution. For the numerical solution of these equations, what is desirable are discretization schemes which satisfy a discrete maximum principle (monotone schemes); otherwise one can expect strong oscillations or even divergency of the solution.There are two possible approaches for development of monotone schemes on unstructured grids. The first approach is to use the grids with some geometric constraints on the triangulation. It is well known that triangulations with no obtuse triangles allow us to construct monotone schemes. The relevant examples are given in [6,11]. However, nonobtuse triangulations can be used on very few practical problems. In the process of the grid generation it is often required to resolve some complicated features of the problem geometry that makes strict angle control to be difficult.Looking for a wider class of triangulations, in the two-dimensional case a Delaunay triangulation [10] is very attractive. For Laplace's equation, in the two-dimensional case the Delaunay triangulation provides positivity (that guarantees the discrete maximum principle) of the linear finite element/finite volume scheme (Barth [3]). This important property may be applied for the solution of a wide range of problems, even more general than discretization of the Laplace equation. For instance, Xu and Zikatanov [17] developed a linear monotone finite element scheme for convectiondiffusion equations in any spatial dimension. To obtain a positive discretization of the diffusion operator they assumed the restrictive geometrical conditions which in the two-dimensional case mean that the triangulation is a Delaunay triangulation.