Equivalence between lowest-order mixed finite element and multi-point finite volume methods on simplicial meshes
Department of Mathematics, Faculty of Nuclear
Sciences and Physical Engineering, Czech Technical University in
Prague, Trojanova 13, 12000 Prague 2, Czech Republic.
2 Laboratoire de Mathématiques, Analyse Numérique et EDP, Université de Paris-Sud, Bât. 425, 91405 Orsay, France. email@example.com
Revised: 16 December 2005
We consider the lowest-order Raviart–Thomas mixed finite element method for second-order elliptic problems on simplicial meshes in two and three space dimensions. This method produces saddle-point problems for scalar and flux unknowns. We show how to easily and locally eliminate the flux unknowns, which implies the equivalence between this method and a particular multi-point finite volume scheme, without any approximate numerical integration. The matrix of the final linear system is sparse, positive definite for a large class of problems, but in general nonsymmetric. We next show that these ideas also apply to mixed and upwind-mixed finite element discretizations of nonlinear parabolic convection–diffusion–reaction problems. Besides the theoretical relationship between the two methods, the results allow for important computational savings in the mixed finite element method, which we finally illustrate on a set of numerical experiments.
Mathematics Subject Classification: 76M10 / 76M12 / 76S05
Key words: Mixed finite element method / saddle-point problem / finite volume method / second-order elliptic equation / nonlinear parabolic convection–diffusion–reaction equation.
© EDP Sciences, SMAI, 2006