A parallel stabilized finite element method based on the lowest equal-order elements for incompressible flows

Based on a fully overlapping domain decomposition technique, a parallel stabilized equal-order finite element method for the steady Stokes equations is presented and studied. In this method, each processor computes a local stabilized finite element solution in its own subdomain by solving a global problem on a global mesh that is locally refined around its subdomain, where the lowest equal-order finite element pairs (continuous piecewise linear, bilinear or trilinear velocity and pressure) are used for the finite element discretization and a pressure-projection-based stabilization method is employed to circumvent the discrete inf-sup condition that is invalid for the used finite element pairs. The parallel stabilized method is unconditionally stable, free of parameter and calculation of derivatives, and is easy to implement based on an existing sequential solver. Optimal error estimates are obtained by the theoretical tool of local a priori error estimates for finite element solutions. Numerical results are also given to verify the theoretical predictions and illustrate the effectiveness of the method.