Convergence of the MAC Scheme for the Stokes/Darcy Coupling Problem

Abstract

In this paper, we extend the MAC scheme for Stokes problem to the Stokes/Darcy coupling problem. The interface conditions between two separate regions are discretized and well-incorporated into the MAC grid setting. We first perform the stability analysis of the scheme for the velocity in both Stokes and Darcy regions and establish the stability for the pressure in both regions by considering an analogue of discrete divergence problem. Following the similar analysis on stability, we perform the error estimates for the velocity and the pressure in both regions. The theoretical results show the first-order convergence of the scheme in discrete L-2 norms for both velocity and the pressure in both regions. Moreover, in fluid region, the first-order convergence for the x-derivative of velocity component u and the y-derivative of velocity component v is also obtained in discrete L-2 norms. However, numerical tests show one order better for the velocity in Stokes region and the pressure in Darcy region.