[A linear scheme satisfying a maximum principle for anisotropic diffusion operators on distorted grids]
We introduce a generalized finite-difference method for anisotropic diffusion operators on distorted grids. We calculate the second-order derivatives in space using a Taylor expansion. The resulting global matrix associated to the scheme is an M-matrix. Thanks to a certain assumption on the grid properties, we show the convergence of the scheme. We show the robustness of the method in comparison with analytical solutions and results obtained by other numerical schemes.
Nous décrivons une méthode de type différences finies généralisées pour des opérateurs de diffusion anisotropes sur des maillages déformés en dimension 2. Les dérivées d'ordre 2 sont calculées à l'aide de développements en séries de Taylor. La matrice associée au schéma est une M-matrice. Une légère hypothèse sur le maillage permet de montrer la convergence du schéma. L'efficacité de la méthode est mise en évidence par comparaison à des solutions analytiques et aux résultats obtenus par d'autres schémas numériques.
Accepted:
Published online:
@article{CRMATH_2009__347_1-2_105_0, author = {Le Potier, Christophe}, title = {Un sch\'ema lin\'eaire v\'erifiant le principe du maximum pour des op\'erateurs de diffusion tr\`es anisotropes sur des maillages d\'eform\'es}, journal = {Comptes Rendus. Math\'ematique}, pages = {105--110}, publisher = {Elsevier}, volume = {347}, number = {1-2}, year = {2009}, doi = {10.1016/j.crma.2008.11.006}, language = {fr}, url = {http://www.numdam.org/articles/10.1016/j.crma.2008.11.006/} }
TY - JOUR AU - Le Potier, Christophe TI - Un schéma linéaire vérifiant le principe du maximum pour des opérateurs de diffusion très anisotropes sur des maillages déformés JO - Comptes Rendus. Mathématique PY - 2009 SP - 105 EP - 110 VL - 347 IS - 1-2 PB - Elsevier UR - http://www.numdam.org/articles/10.1016/j.crma.2008.11.006/ DO - 10.1016/j.crma.2008.11.006 LA - fr ID - CRMATH_2009__347_1-2_105_0 ER -
%0 Journal Article %A Le Potier, Christophe %T Un schéma linéaire vérifiant le principe du maximum pour des opérateurs de diffusion très anisotropes sur des maillages déformés %J Comptes Rendus. Mathématique %D 2009 %P 105-110 %V 347 %N 1-2 %I Elsevier %U http://www.numdam.org/articles/10.1016/j.crma.2008.11.006/ %R 10.1016/j.crma.2008.11.006 %G fr %F CRMATH_2009__347_1-2_105_0
Le Potier, Christophe. Un schéma linéaire vérifiant le principe du maximum pour des opérateurs de diffusion très anisotropes sur des maillages déformés. Comptes Rendus. Mathématique, Volume 347 (2009) no. 1-2, pp. 105-110. doi : 10.1016/j.crma.2008.11.006. http://www.numdam.org/articles/10.1016/j.crma.2008.11.006/
[1] The Andra Couplex1 test: comparisons between finite element, mixed hybrid finite element and finite volume element discretizations: simulation of transport around a nuclear waste disposal site, Comput. Geosci., Volume 8 (2004) no. 2, pp. 187-201
[2] Discrete maximum principle for Galerkin approximations of the Laplace operator on arbitrary meshes, C. R. Acad. Sci. Paris, Ser. I, Volume 338 (2004) no. 8, pp. 641-646
[3] F. Dabbène, Mixed hybrid finite elements for transport of polluants by undergrounds water, in: Proceeding of the 10th International Conference on Finite Elements in Fluids, Tucson, USA, 1998
[4] Maximization of a linear function of variables subject to linear inequalities (Koopmans, T.C., ed.), Activity Analysis of Production and Allocation, Wiley, New York, 1951, pp. 339-347 (Chapter 21)
[5] Finite volume method (Ciarlet, P.G.; Lions, J.L., eds.), Handbook of Numerical Analysis, vol. VII, 2000
[6] I. Faille, Modélisation bidimensonnelle de la génèse et de la migration des hydrocarbures dans un bassin sédimentaire, Thèse de l'université Joseph Fourier-Grenoble 1, 1992
[7] Improvements of generalized finite difference method and comparison with other meshless method, Appl. Math. Modelling, Volume 27 (2003) no. 10, pp. 831-847
[8] A. Genty, C. Le Potier, Avoiding negative concentrations for transport calculations of radionuclides in high-level radioactive waste disposal system, in: IAHR, International Groundwater Symposium, Istanbul, 2008, pp. 286–292
[9] Benchmark on discretization schemes for anisotropic diffusion problems on general grids, 2008 http://www.latp.univ-mrs.fr/fvca5 (in: 5th International Symposium on Finite Volumes for Complex Applications, June 08–13)
[10] Finite volume scheme for highly anisotropic diffusion operators on unstructured meshes, C. R. Acad. Sci. Paris, Ser. I, Volume 340 (2005), pp. 921-926
[11] Finite volume monotone scheme for highly anisotropic diffusion operators on unstructured triangular meshes, C. R. Acad. Sci. Paris, Ser. I, Volume 341 (2005), pp. 787-792
[12] C. Le Potier, Finite volume scheme satisfying maximum and minimum principles for anisotropic diffusion operators, in: Finite Volumes for Complex Applications V, Aussois, 2008, pp. 103–118
[13] A monotone finite element scheme for convection–diffusion equations, Math. Comp., Volume 66 (1999) no. 228, pp. 1429-1446
Cited by Sources: