Convergence analysis of mixed numerical schemes for reactive flow in a porous medium

K. Kumar, I.S. Pop, F.A. Radu

    Research output: Contribution to journalArticleAcademicpeer-review

    23 Citations (Scopus)
    100 Downloads (Pure)


    This paper deals with the numerical analysis of an upscaled model describing the reactive flow in a porous medium. The solutes are transported by advection and diffusion and undergo precipitation and dissolution. The reaction term and, in particular, the dissolution term have a particular, multivalued character, which leads to stiff dissolution fronts. We consider the Euler implicit method for the temporal discretization and the mixed finite element for the discretization in space. More precisely, we use the lowest order Raviart--Thomas elements. As an intermediate step we consider also a semidiscrete mixed variational formulation (continuous in space). We analyze the numerical schemes and prove the convergence to the continuous formulation. Apart from the proof for the convergence, this also yields an existence proof for the solution of the model in mixed variational formulation. Numerical experiments are performed to study the convergence behavior. Keywords: numerical analysis, reactive flows, weak formulation, implicit scheme, mixed finite element discretization
    Original languageEnglish
    Pages (from-to)2283-2308
    Number of pages26
    JournalSIAM Journal on Numerical Analysis
    Issue number4
    Publication statusPublished - 2013


    Dive into the research topics of 'Convergence analysis of mixed numerical schemes for reactive flow in a porous medium'. Together they form a unique fingerprint.

    Cite this