The vortex particle method, which includes viscous effects, consists of diffusion of boundary vorticity and creation of the vortex particles, convection, particle strength exchange, and particle redistribution. Accuracy of the boundary element method is very important since it creates the particles around the body at every time step. A boundary element method based on source panel was investigated as part of computation of unsteady separated flows by rising the vortex particle method. The potential flows were computed around a circular cylinder and a square cylinder. The results around the circular cylinder were compared with the exact solution, and the distribution of vorticity, in particular near the sharp comers of the square cylinder, is scrutinized for different number of panels.