A locally conservative LDG method for the incompressible Navier-Stokes equations

Bernardo Cockburn, Guido Kanschat, Dominik Schötzau

Research output: Contribution to journalArticlepeer-review

247 Scopus citations


In this paper a new local discontinuous Galerkin method for the incompressible stationary Navier-Stokes equations is proposed and analyzed. Four important features render this method unique: its stability, its local conservativity, its high-order accuracy, and the exact satisfaction of the incompressibility constraint. Although the method uses completely discontinuous approximations, a globally divergence-free approximate velocity in H(div; Ω) is obtained by simple, element-by-element post-processing. Optimal error estimates are proven and an iterative procedure used to compute the approximate solution is shown to converge. This procedure is nothing but a discrete version of the classical fixed point iteration used to obtain existence and uniqueness of solutions to the incompressible Navier-Stokes equations by solving a sequence of Oseen problems. Numerical results are shown which verify the theoretical rates of convergence. They also confirm the independence of the number of fixed point iterations with respect to the discretization parameters. Finally, they show that the method works well for a wide range of Reynolds numbers.

Original languageEnglish (US)
Pages (from-to)1067-1095
Number of pages29
JournalMathematics of Computation
Issue number251
StatePublished - Jul 2005


  • Discontinuous Galerkin methods
  • Finite element methods
  • Incompressible Navier-Stokes equations


Dive into the research topics of 'A locally conservative LDG method for the incompressible Navier-Stokes equations'. Together they form a unique fingerprint.

Cite this