A projection method for constructing a mass conservative velocity field