Virtual element approximations for two species model of the Advection-Diffusion-Reaction in poroelastic media

This paper proposes virtual element methods for approximating the mathematical model consisting of coupled poroelastic and Advection-Diffusion-Reaction (ADR) equations. The space discretization relies on virtual element spaces containing piecewise linear polynomials as well as non-polynomials for displacement, pressure and concentrations, and piecewise constant for total pressure; a backwardEuler scheme is employed for the approximation of time derivative. Using standard techniques of explicit schemes, the well-posedness of the resultant fully discrete scheme is derived. Moreover, under certain regularity assumptions on the mesh, optimal apriori error estimates are established by introducing suitable projection operators. Several numerical experiments are presented to validate the theoretical convergence rate and exhibit the proposed scheme’s performance.