We propose a cell-centered symmetric scheme which combines the advantages of MPFA (multipoint flux approximation) schemes such as the L or the O scheme and of hybrid schemes: it may be used on general non-conforming meshes, it yields a 9-point stencil on two-dimensional quadrangular meshes, it takes into account the heterogeneous diffusion matrix, it is coercive and it can be shown to converge. The scheme relies on the use of special points, called harmonic averaging points, located at the interfaces of heterogeneity.
Nous proposons un schéma ayant ses inconnues aux centres des mailles, combinant les avantages des schémas à flux multi-points et des schémas hybrides : il possède un stencil à 9 points en 2D, respecte les hétérogénéités de la matrice de diffusion, et il est coercif ; de plus, on peut montrer qu'il converge. Le schéma est basé sur l'utilisation de points situé aux interfaces d'hétérogénéité, en lesquels la formule de la moyenne harmonique est utilisable.
Accepted:
Published online:
@article{CRMATH_2009__347_11-12_673_0, author = {Agelas, L\'eo and Eymard, Robert and Herbin, Rapha\`ele}, title = {A nine-point finite volume scheme for the simulation of diffusion in heterogeneous media}, journal = {Comptes Rendus. Math\'ematique}, pages = {673--676}, publisher = {Elsevier}, volume = {347}, number = {11-12}, year = {2009}, doi = {10.1016/j.crma.2009.03.013}, language = {en}, url = {http://www.numdam.org/articles/10.1016/j.crma.2009.03.013/} }
TY - JOUR AU - Agelas, Léo AU - Eymard, Robert AU - Herbin, Raphaèle TI - A nine-point finite volume scheme for the simulation of diffusion in heterogeneous media JO - Comptes Rendus. Mathématique PY - 2009 SP - 673 EP - 676 VL - 347 IS - 11-12 PB - Elsevier UR - http://www.numdam.org/articles/10.1016/j.crma.2009.03.013/ DO - 10.1016/j.crma.2009.03.013 LA - en ID - CRMATH_2009__347_11-12_673_0 ER -
%0 Journal Article %A Agelas, Léo %A Eymard, Robert %A Herbin, Raphaèle %T A nine-point finite volume scheme for the simulation of diffusion in heterogeneous media %J Comptes Rendus. Mathématique %D 2009 %P 673-676 %V 347 %N 11-12 %I Elsevier %U http://www.numdam.org/articles/10.1016/j.crma.2009.03.013/ %R 10.1016/j.crma.2009.03.013 %G en %F CRMATH_2009__347_11-12_673_0
Agelas, Léo; Eymard, Robert; Herbin, Raphaèle. A nine-point finite volume scheme for the simulation of diffusion in heterogeneous media. Comptes Rendus. Mathématique, Volume 347 (2009) no. 11-12, pp. 673-676. doi : 10.1016/j.crma.2009.03.013. http://www.numdam.org/articles/10.1016/j.crma.2009.03.013/
[1] Discretization on non-orthogonal, quadrilateral grids for inhomogeneous, anisotropic media, J. Comput. Phys., Volume 127 (1996) no. 1, pp. 2-14
[2] Discretization on unstructured grids for inhomogeneous, anisotropic media. Part i: Derivation of the methods, SIAM J. Sci. Comput., Volume 19 (1998), pp. 1700-1716
[3] Direct Methods for Sparse Matrices, Monographs on Numerical Analysis, vol. XIII, Clarendon Press, Oxford, 1986 (341 pp)
[4] A new finite volume scheme for anisotropic diffusion problems on general grids: convergence analysis, C. R. Math. Acad. Sci. Paris, Volume 344 (2007) no. 6, pp. 403-406
[5] Benchmark on anisotropic problems, SUSHI: a scheme using stabilization and hybrid interfaces for anisotropic heterogeneous diffusion problems (Eymard, R.; Hérard, J.-M., eds.), Finite Volumes for Complex Applications, vol. V, Wiley, 2008, pp. 801-814
[6] R. Eymard, T. Gallouët, R. Herbin, Discretization of heterogeneous and anisotropic diffusion problems on general non-conforming meshes, SUSHI: a scheme using stabilization and hybrid interfaces, IMAJNA, 2009, in press; see also http://hal.archives-ouvertes.fr/
[7] Benchmark on discretization schemes for anisotropic diffusion problems on general grids for anisotropic heterogeneous diffusion problems (Eymard, R.; Hérard, J.-M., eds.), Finite Volumes for Complex Applications, vol. V, Wiley, 2008, pp. 659-692
Cited by Sources:
☆ Work supported by Groupement de Recherche MOMAS, PACEN/CNRS.