A FAST SCHUR COMPLEMENT METHOD FOR THE SPECTRAL ELEMENT DISCRETIZATION OF THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

The weak formulation of the incompressible Navier-Stokes equations in three space dimensions is discretized with spectral element approximat ions and Gauss-Lobatto-Legendre quadratures. The Uzawa algorithm is ap plied to decouple the Velocities from the pressure. The equation that results for the pressure is solved by an iterative method. Within each pressure iteration, a Helmholtz operator has to be inverted. This can efficiently be done by separating the equations for the interior node s from the equations at the interfaces, according to the Schur method. Fast diagonalization techniques are applied to the interior variables of the spectral elements. Several ways to deal with the resulting int erface problem are discussed. Finally, a comparison is made with a mor e classical method. (C) 1995 Academic Press, Inc.