A linearity-preserving vertex interpolation algorithm for cell-centered finite volume approximations of anisotropic diffusion problems

Purpose A finite volume scheme for diffusion equations on non-rectangular meshes is proposed in [Deyuan Li, Hongshou Shui, Minjun Tang, J. Numer. Meth. Comput. Appl., 1(4)(1980)217-224 (in Chinese)], which is the so-called nine point scheme on structured quadrilateral meshes. The scheme has both cell-centered unknowns and vertex unknowns which are usually expressed as a linear weighted interpolation of the cell-centered unknowns. The critical factor to obtain the optimal accuracy for the scheme is the reconstruction of vertex unknowns. However, when the mesh deformation is severe or the diffusion tensor is discontinuous, the accuracy of the scheme is not satisfactory, and the author hope to improve this scheme. Design/methodology/approach The authors propose an explicit weighted vertex interpolation algorithm which allows arbitrary diffusion tensors and does not depend on the location of discontinuity. Both the derivation of the scheme and that of vertex reconstruction algorithm satisfy the linearity preserving criterion which requires that a discretization scheme should be exact on linear solutions. The vertex interpolation algorithm can be easily extended to 3 D case. Findings Numerical results show that it maintain optimal convergence rates for the solution and flux on 2 D and 3 D meshes in case that the diffusion tensor is taken to be anisotropic, at times heterogeneous, and/or discontinuous. Originality/value This paper proposes a linearity preserving and explicit weighted vertex interpolation algorithm for cell-centered finite volume approximations of diffusion equations on general grids. The proposed finite volume scheme with the new interpolation algorithm allows arbitrary continuous or discontinuous diffusion tensors; the final scheme is applicable to arbitrary polygonal grids, which may have concave cells or degenerate ones with hanging nodes. The final scheme has second-order convergence rate for the approximate solution and higher than first-order accuracy for the flux on 2 D and 3 D meshes. The explicit weighted interpolation algorithm is easy to implement in three dimensions in case that the diffusion tensor is continuous or discontinuous.