Research Article | | Peer-Reviewed

A New Second-order Maximum-principle-preserving Finite-volume Method for Flow Problems Involving Discontinuous Coefficients

Received: 24 June 2024     Accepted: 30 July 2024     Published: 26 August 2024
Views:       Downloads:
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.

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

Keywords

Diffusion-reaction Problems, New Finite Volume Scheme, Diffusion Coefficient Discontinuities, Second Order Convergence

References
[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.
Cite This Article
  • 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

    Copy | Download

    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

    Copy | Download

    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

    Copy | Download

  • @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}
    }
    

    Copy | Download

  • 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  - 

    Copy | Download

Author Information
  • Department of Mathematics, National Advanced School of Engineering, University of Yaounde 1, Yaounde, Cameroon

  • Department of Mathematics, Faculty of Science, University of Douala, Douala, Cameroon

  • Department of Civil Engineering and Architecture, Ecole Nationale Superieure Polytechnique de Maroua, Maroua, Cameroon

  • Sections