An asymptotically optimal schur complement reduction for the Stokes equation

Boris N. Khoromskij, Gabriel Wittum*

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

Abstract

In this paper we develop an efficient Schur complement method for solving the 2D Stokes equation. As a basic algorithm, we apply a decomposition approach with respect to the trace of the pressure. The alternative stream function-vorticity reduction is also discussed. The original problem is reduced to solving the equivalent boundary (interface) equation with symmetric and positive definite operator in the appropriate trace space. We apply a mixed finite element approximation to the interface operator by P1 iso P2/P1 triangular elements and prove the optimal error estimates in the presence of stabilizing bubble functions. The norm equivalences for the corresponding discrete operators are established. Then we propose an asymptotically optimal compression technique for the related stiffness matrix (in the absence of bubble functions) providing a sparse factorized approximation to the Schur complement. In this case, the algorithm is shown to have an optimal complexity of the order O(N logq N), q = 2 or q = 3, depending on the geometry, where N is the number of degrees of freedom on the interface. In the presence of bubble functions, our method has the complexity O(N2 log N) arithmetical operations. The Schur complement interface equation is resolved by the PCG iterations with an optimal preconditioner.

Original languageEnglish (US)
Pages (from-to)345-375
Number of pages31
JournalNumerische Mathematik
Volume81
Issue number3
DOIs
StatePublished - Jan 1 1999

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics

Fingerprint Dive into the research topics of 'An asymptotically optimal schur complement reduction for the Stokes equation'. Together they form a unique fingerprint.

Cite this