We develop a new method for the efficient calculation of solenoidal vector fields on general regions. The method takes advantage of fast direct methods and uses boundary integral equations to satisfy boundary conditions. For the latter we give an effective scheme for computing far-field boundary influences (based on discrete charges). Examples and numerical results are given. The method is applicable to incompressible Navier-Stokes calculations.