A unified stabilized method for Stokes’ and Darcy's equations

We use the lowest possible approximation order, piecewise linear, continuous velocities and piecewise constant pressures to compute solutions to Stokes equation and Darcy's equation, applying an edge stabilization term to avoid locking. We prove that the formulation satisfies the discrete inf-sup condition, we prove an optimal a priori error estimate for both problems. The formulation is then extended to the coupled case using a Nitsche-type weak formulation allowing for different meshes in the two subdomains. Finally, we present some numerical examples verifying the theoretical predictions and showing the flexibility of the coupled approach.