Abstract
Approximations to solutions of the inhomogeneous boundary value problem for the Navier-Stokes equations are constructed via the discontinuous Galerkin method. The velocity held is approximated using piecewise polynomial functions that are totally discontinuous across interelement boundaries and which are pointwise divergence-free on each dement (locally solenoidal). The pressure is approximated by standard continuous piecewise polynomial functions.