A fast solver for the Stokes equations with distributed forces in complex geometries
Abstract
We present a new method for the solution of the Stokes equations. The main features of our method are: (1) it can be applied to arbitrary geometries in a blackbox fashion; (2) it is secondorder accurate; and (3) it has optimal algorithmic complexity. Our approach, to which we refer as the embedded boundary integral method (EBI), is based on Anita Mayo's work for the Poisson's equation: "The Fast Solution of Poisson's and the Biharmonic Equations on Irregular Regions", SIAM Journal on Numerical Analysis, 21 (1984) 285299. We embed the domain in a rectangular domain, for which fast solvers are available, and we impose the boundary conditions as interface (jump) conditions on the velocities and tractions. We use an indirect boundary integral formulation for the homogeneous Stokes equations to compute the jumps. The resulting equations are discretized by Nyström's method. The rectangular domain problem is discretized by finite elements for a velocitypressure formulation with equal order interpolation bilinear elements ( Q1 Q1). Stabilization is used to circumvent the inf sup condition for the pressure space. For the integral equations, fast matrixvector multiplications are achieved via an Nlog N algorithm based on a block representation of the discrete integral operator, combined with (kernel independent) singular value decomposition to sparsify lowrank blocks. The regular grid solver is a Krylov method (conjugate residuals) combined with an optimal twolevel Schwartzpreconditioner. For the integral equation we use GMRES. We have tested our algorithm on several numerical examples and we have observed optimal convergence rates.
 Publication:

Journal of Computational Physics
 Pub Date:
 January 2004
 DOI:
 10.1016/j.jcp.2003.08.011
 Bibcode:
 2004JCoPh.193..317B