摘要

In this article, we introduce and analyze a weak Galerkin finite element method for numerically solving the coupling of fluid flow with porous media flow. Flows are governed by the Stokes equations in primal velocity-pressure formulation and Darcy equation in the second order primary formulation, respectively, and the corresponding transmission conditions are given by mass conservation, balance of normal forces, and the Beavers-Joseph-Saffman law. By using the weak Galerkin approach, we consider the two-dimensional problem with the piecewise constant elements for approximations of the velocity, pressure, and hydraulic head. Stability and optimal error estimates are obtained. Finally, we provide several numerical results illustrating the good performance of the proposed scheme and confirming the optimal order of convergence provided by the weak Galerkin approximation.