In this article, we analyze an Euler implicit-mixed finite element scheme for a porous media solute transport model. The transporting flux is not assumed given, but obtained by solving numerically the Richards equation, a model for subsurface fluid flow. We prove the convergence of the scheme by estimating the error in terms of the discretization parameters. In doing so we take into account the numerical error occurring in the approximation of the fluid flow. The article, is concluded by numerical experiments, which are in good agreement with the theoretical estimates.
|Journal||Numerical Methods for Partial Differential Equations|
|Publication status||Published - 2010|