On incompressible flow simulations using Scott-Vogelius finite elements

Please always quote using this URN: urn:nbn:de:0296-matheon-11661
  • Recent research has shown that in some practically relevant situations like multi-physics flows[11] divergence-free mixed finite elements may have a significantly smaller discretization error than standard non-divergence-free mixed finite elements. In order to judge the overall performance of divergence-free mixed finite elements, we investigate linear solvers for the saddle point linear systems arising in $((P_k)^d,P_{k-1}^{disc})$ Scott-Vogelius finite element implementations of the incompressible Navier-Stokes equations. We investigate both direct and iterative solver methods. Due to discontinuous pressure elements in the case of Scott-Vogelius elements, considerably more solver strategies seem to deliver promising results than in the case of standard mixed finite elements like Taylor-Hood elements. For direct methods, we extend recent preliminary work using sparse banded solvers on the penalty method formulation to finer meshes, and discuss extensions. For iterative methods, we test augmented Lagrangian and H-LU preconditioners with GMRES, on both full and statically condensed systems. Several numerical experiments are provided that show these classes of solvers are well suited for use with Scott-Vogelius elements, and could deliver an interesting overall performance in several applications.

Export metadata

  • Export Bibtex
  • Export RIS

Additional Services

Share in Twitter Search Google Scholar
Author:Ben Cousins, Sabine Le Borne, Alexander Linke, Leo G. Rebholz, Zhen Wang
Referee:Alexander Mielke
Date of first Publication:2012/09/01
Institute:Freie Universit├Ąt Berlin
MSC-Classification:35-XX PARTIAL DIFFERENTIAL EQUATIONS / 35Qxx Equations of mathematical physics and other areas of application [See also 35J05, 35J10, 35K05, 35L05] / 35Q30 Navier-Stokes equations [See also 76D05, 76D07, 76N10]
65-XX NUMERICAL ANALYSIS / 65Fxx Numerical linear algebra / 65F05 Direct methods for linear systems and matrix inversion
65-XX NUMERICAL ANALYSIS / 65Fxx Numerical linear algebra / 65F10 Iterative methods for linear systems [See also 65N22]
65-XX NUMERICAL ANALYSIS / 65Nxx Partial differential equations, boundary value problems / 65N30 Finite elements, Rayleigh-Ritz and Galerkin methods, finite methods
Preprint Number:983

$Rev: 13581 $