A parallel two-level Grad-div stabilized finite element algorithm for steady incompressible Stokes equations is proposed. Basic idea of the algorithm is to solve global Grad-div stabilized solution in a coarse mesh firstly, and then correct it in parallel on overlapping fine mesh subdomains. With reasonable selection of stabilization parameters and mesh sizes, an optimal convergence rate can be obtained. Numerical results verify efficiency of the algorithm.
