Weak Galerkin method for the coupled Darcy-Stokes flow

A family of weak Galerkin finite element discretization is developed for solving the coupled Darcy-Stokes equation. The equation in consideration admits the Beaver-Joseph-Saffman condition on the interface. By using the weak Galerkin approach, in the discrete space we are able to impose the normal continuity of velocity explicitly. Or in other words, strong coupling is achieved in the discrete space. Different choices of weak Galerkin finite element spaces are discussed, and error estimates are given.