Analysis and Computation of a Weak Galerkin Scheme for Solving the 2D/3D Stationary Stokes Interface Problems with High-Order Elements.

In this paper, we present a high-order weak Galerkin finite element method (WG-FEM) for solving the stationary Stokes interface problems with discontinuous velocity and pressure in ℝ<italic>d</italic> (<italic>d</italic> = 2, 3). This WG method is equipped with stable finite elements consisting of usual polynomials of degree <italic>k</italic> ≥ 1 for the velocity and polynomials of degree <italic>k</italic> – 1 for the pressure, both are discontinuous. Optimal convergence rates of order <italic>k</italic> + 1 for the velocity and order <italic>k</italic> for the pressure are established in <italic>L</italic>2-norm on hybrid meshes. Numerical experiments verify the expected order of accuracy for both two-dimensional and three-dimensional examples. Moreover, numerically it is shown that the proposed WG algorithm is able to accommodate geometrically complicated and very irregular interfaces having sharp edges, cusps, and tips. [ABSTRACT FROM AUTHOR]