Efficient implementation of the exact artificial boundary condition for the exterior problem of the Stokes system in three dimensions

In this paper the Stokes system in an unbounded domain is solved by the artificial boundary method. The novelty lies in an operator form of the exact Dirichlet-to-Neumann (DtN) mapping. With the help of the Chebyshev rational approximation of the square root function, we derive a highly accurate approximate DtN mapping, which can be numerically implemented without resorting to the eigen-decomposition in terms of the vectorial spherical harmonics. In addition, we develop an efficient block preconditioner for the augmented truncated saddle point problem. Numerical experiments demonstrate the effectiveness of the proposed method.