A new development of Finite Volumes (FV, for short) and its theoretical analysis are the purpose of this work. Recall that FV are known as powerful tools to address equations of conservation laws (mass, energy, momentum,...). Over the last two decades investigators have succeeded in putting in place a mathematical framework for the theoretical analysis of FV. A perfect illustration of this progress is the design and mathematical analysis of Discrete Duality Finite Volumes (DDFV, for short). We propose now a new class of DDFV for 2nd order elliptic equations involving discontinuous diffusion coefficients or nonlinearities. A one-dimensional linear elliptic equation is addressed here for illustrating the ideas behind our numerical strategy. The algebraic structure of the discrete system we have got is different from that of standard DDFV. The main novelty is that the so-called diamond mesh elements are confined in homogeneous zones for flow problems governed by piecewise constant coefficients. This is got from our judicious definition of the primal mesh. The gain is that there is no need to compute homogenized coefficients to be allocated to the so-called diamond cells as required to conventional DDFV. Notice that poor homogenized permeability allocated to diamond elements leads to poor approximations of fluxes across grid-block interfaces. Moreover for 1-D flow problems in a porous medium involving permeability discontinuities (piecewise constant permeability for instance) the proposed FV scheme leads to a symmetric positive-definite discrete system that meets the discrete maximum principle; we have shown its second order convergence under relevant assumptions.
Published in | American Journal of Applied Mathematics (Volume 12, Issue 4) |
DOI | 10.11648/j.ajam.20241204.12 |
Page(s) | 91-110 |
Creative Commons |
This is an Open Access article, distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium or format, provided the original work is properly cited. |
Copyright |
Copyright © The Author(s), 2024. Published by Science Publishing Group |
Diffusion-reaction Problems, New Finite Volume Scheme, Diffusion Coefficient Discontinuities, Second Order Convergence
[1] | C. M. Marle, Cours de production: Les écoulements polyphasiques en milieu poreux, Tome IV, Ed.Technip, 1984, 300 pages. |
[2] | P. Lemonnier, and B. Bourbiaux, Simulation of Naturally Fractured Reservoirs. State of the Art, Oil / Gas Science and Technology, Review of ”Institut Francais du Pétrole”, March 2010. |
[3] | S. C. Brenner, and L. R. Scott, The Mathematical Theory of Finite Element Methods, Book series: Texts in Applied Mathematics Springer (Editor), 2008. |
[4] | A. Njifenjou, Overview on conventional finite volumes for elliptic problems involving discontinuous diffusion coefficients. Part I: Focus on the one dimension space models, Research Gate Preprint, September 2022. |
[5] | H. Brézis, Functional Analysis, Sobolev spaces and Partial Differential Equations, Springer (2011). |
[6] | A. Njifenjou, Introduction to Finite Element Methods, Research Gate Preprint, September 2021. |
[7] | A. Njifenjou, Geometrical technique for proving the discrete maximum principle for conventional finite volumes, Research Gate preprint (version1), April 2023. |
[8] | K. Domelevo, and P. Omnes, A finite volume method for the Laplace equation on almost arbitrary two-dimensional grids. ESAIM, Math. Model. Numer. Anal. 39, No6, 1203-1249, 2005. |
[9] | F. Hermeline, A finite volume method for the approximation of diffusion operators on distorted meshes. J. Comput. Phys., 160, No2, pp 481-499, 2000. |
[10] | A. Njifenjou, and I. Moukouop-Nguena, Traitement des anisotropies de perméabilité en simulation d’écoulement en milieu poreux par les volumes finis. Proceedings of an International Conference on ”Systèmes Informatiques pourla Gestiondel’ Environnement”, (MauriceTchuente, Ed.), Douala (Cameroon), pp 245-259. |
[11] | A. Njifenjou, H. Donfack, and I. Moukouop-Nguena, Analysis on general meshes of a discrete duality finite volume method for subsurface flow problems, Computational Geosciences, 17, (2013), pp 391-415. |
[12] | F. Hermeline, A finite volume method for approximating 3D diffusion operators on general meshes, J. Comput. Phys., Vol. 228, No17, (2009), pp 5763-5786. |
[13] | A. Kinfack Jeutsa, A. Njifenjou, and J. Nganhou, Convergence Analysis on Unstructured Meshes of a DDFV Method for Flow Problems with Full Neumann Boundary Conditions, Journal of Applied Mathematics, Hindawi Ed., Volume 2016, September 2016. |
[14] | A. Kinfack Jeutsa, A. Njifenjou, and J. Nganhou, Convergence Analysis on Quadrilateral Grids of a DDFV Method for Subsurface Flow Problems in Anisotropic Heterogeneous Porous Media with Full Neumann Boundary Conditions, African Diaspora Journal of Mathematics, Volume 19, Number 2, pp.128 (2016). |
[15] | A. Njifenjou, H. Donfack, and A. Kinfack Jeutsa, A discrete duality finite volume method for flow problems with prescribed periodic boundary conditions. HAL open science. 2019. hal-02266479. |
[16] | F. Boyer, and F. Nabet, A DDFV method for a Cahn-Hilliard/Stokes phase field model with dynamic boundary conditions, ESAIM: Mathematical Modelling and Numerical Analysis, EDP Sciences, 2017, 51 (5), pp. 1691-1731. |
[17] | S. Delcourte, K. Domelevo and P. Omnes, A Discrete Duality Finite Volume Approach to Hodge Decomposition and div-curl Problems on Almost Arbitrary Two-Dimensional Meshes, SIAM Journal on Numerical Analysis, Vol.45, 3 (2007). |
[18] | P. Omnes, A second-order convergence of a function reconstructed from finite volume approximation of the Laplace equation on Delaunay-Voronoi meshes, ESAIM:M2AN, Vol. 45, No4(2011), pp 627-650. |
[19] | G. Galindez, Darlan Karlo Elisiário de Carvalho, and P. R. M. Lyra, Numerical simulation of 1-D oil and water displacements in petroleum reservoirs using the correction procedure via reconstruction (CPR) method, Computational Geosciences, February 2020, 24(2). |
[20] | P G. Ciarlet, The finite element method for elliptic problems, North-Holland, Amsterdam (1978). |
[21] | R. Eymard, T. Gallouet, and R. Herbin, Finite volume methods. In P. G. Ciarlet and J.-L. Lions, editors, Techniques of Scientific Computing, Part III, Handbook of Numerical Analysis. North-Holland, Amsterdam, 2000. |
[22] | F. Brezzi, K. Lipnikov, and M. Shashkov, Convergence of the mimetic finite difference method for diffusion problems on polyhedral meshes, SIAM J. Numer. Anal., 43(5):1872 1896, 2005. |
[23] | R. Herbin and F. Hubert, Benchmark on discretization schemes for anisotropic diffusion problems on general grids. In Finite volumes for complex applications V. ISTE, London, 2008. |
[24] | F. Brezzi, K. Lipnikov, and M. Shashkov, Convergence of the mimetic finite difference method for diffusion problems on polyhedral meshes with curved faces, Mathematical Models and Methods in Applied Sciences, 16(02): 275-297, 2006. |
[25] | Ph. G. Ciarlet, Introduction a l’analyse numérique matricielle et a l’optimisation, Masson (Editor), 1982. |
APA Style
Njifenjou, A., Mansou, A. T., Sali, M. (2024). A New Second-order Maximum-principle-preserving Finite-volume Method for Flow Problems Involving Discontinuous Coefficients. American Journal of Applied Mathematics, 12(4), 91-110. https://doi.org/10.11648/j.ajam.20241204.12
ACS Style
Njifenjou, A.; Mansou, A. T.; Sali, M. A New Second-order Maximum-principle-preserving Finite-volume Method for Flow Problems Involving Discontinuous Coefficients. Am. J. Appl. Math. 2024, 12(4), 91-110. doi: 10.11648/j.ajam.20241204.12
AMA Style
Njifenjou A, Mansou AT, Sali M. A New Second-order Maximum-principle-preserving Finite-volume Method for Flow Problems Involving Discontinuous Coefficients. Am J Appl Math. 2024;12(4):91-110. doi: 10.11648/j.ajam.20241204.12
@article{10.11648/j.ajam.20241204.12, author = {Abdou Njifenjou and Abel Toudna Mansou and Moussa Sali}, title = {A New Second-order Maximum-principle-preserving Finite-volume Method for Flow Problems Involving Discontinuous Coefficients}, journal = {American Journal of Applied Mathematics}, volume = {12}, number = {4}, pages = {91-110}, doi = {10.11648/j.ajam.20241204.12}, url = {https://doi.org/10.11648/j.ajam.20241204.12}, eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.ajam.20241204.12}, abstract = {A new development of Finite Volumes (FV, for short) and its theoretical analysis are the purpose of this work. Recall that FV are known as powerful tools to address equations of conservation laws (mass, energy, momentum,...). Over the last two decades investigators have succeeded in putting in place a mathematical framework for the theoretical analysis of FV. A perfect illustration of this progress is the design and mathematical analysis of Discrete Duality Finite Volumes (DDFV, for short). We propose now a new class of DDFV for 2nd order elliptic equations involving discontinuous diffusion coefficients or nonlinearities. A one-dimensional linear elliptic equation is addressed here for illustrating the ideas behind our numerical strategy. The algebraic structure of the discrete system we have got is different from that of standard DDFV. The main novelty is that the so-called diamond mesh elements are confined in homogeneous zones for flow problems governed by piecewise constant coefficients. This is got from our judicious definition of the primal mesh. The gain is that there is no need to compute homogenized coefficients to be allocated to the so-called diamond cells as required to conventional DDFV. Notice that poor homogenized permeability allocated to diamond elements leads to poor approximations of fluxes across grid-block interfaces. Moreover for 1-D flow problems in a porous medium involving permeability discontinuities (piecewise constant permeability for instance) the proposed FV scheme leads to a symmetric positive-definite discrete system that meets the discrete maximum principle; we have shown its second order convergence under relevant assumptions.}, year = {2024} }
TY - JOUR T1 - A New Second-order Maximum-principle-preserving Finite-volume Method for Flow Problems Involving Discontinuous Coefficients AU - Abdou Njifenjou AU - Abel Toudna Mansou AU - Moussa Sali Y1 - 2024/08/26 PY - 2024 N1 - https://doi.org/10.11648/j.ajam.20241204.12 DO - 10.11648/j.ajam.20241204.12 T2 - American Journal of Applied Mathematics JF - American Journal of Applied Mathematics JO - American Journal of Applied Mathematics SP - 91 EP - 110 PB - Science Publishing Group SN - 2330-006X UR - https://doi.org/10.11648/j.ajam.20241204.12 AB - A new development of Finite Volumes (FV, for short) and its theoretical analysis are the purpose of this work. Recall that FV are known as powerful tools to address equations of conservation laws (mass, energy, momentum,...). Over the last two decades investigators have succeeded in putting in place a mathematical framework for the theoretical analysis of FV. A perfect illustration of this progress is the design and mathematical analysis of Discrete Duality Finite Volumes (DDFV, for short). We propose now a new class of DDFV for 2nd order elliptic equations involving discontinuous diffusion coefficients or nonlinearities. A one-dimensional linear elliptic equation is addressed here for illustrating the ideas behind our numerical strategy. The algebraic structure of the discrete system we have got is different from that of standard DDFV. The main novelty is that the so-called diamond mesh elements are confined in homogeneous zones for flow problems governed by piecewise constant coefficients. This is got from our judicious definition of the primal mesh. The gain is that there is no need to compute homogenized coefficients to be allocated to the so-called diamond cells as required to conventional DDFV. Notice that poor homogenized permeability allocated to diamond elements leads to poor approximations of fluxes across grid-block interfaces. Moreover for 1-D flow problems in a porous medium involving permeability discontinuities (piecewise constant permeability for instance) the proposed FV scheme leads to a symmetric positive-definite discrete system that meets the discrete maximum principle; we have shown its second order convergence under relevant assumptions. VL - 12 IS - 4 ER -