Numerical simulation of incompressible fluid flow using locally solenoidal elements

In this paper, we present a numerical study of the performance of a discontinuous Galerkin formulation for the Navier-Stokes equations. This method is characterized by the fact that the velocity field is approximated using piecewise polynomial functions that are totally discontinuous across interelement boundaries and which are pointwise divergence-free on each element (locally solenoidal In particular, numerical results are presented for two well-known benchmark problems.