  In this paper we will prove that suitable weak solutions of three dimensional
Navier-Stokes equations in bounded domain can be constructed by a particular
type of artificial compressibility approximation.
