The aim of this method is focus on developments of the innovative numerical methods, which based on the algorithm of particle swarm optimization (PSO) to combine with the method of fundamental solutions (MFS) for construction the proposed novel numerical models of potential flow field. The proposed model is a kind boundary-type meshless method for simulation the flow field with different scenarios then to calculate the forcing. The so-called method of fundamental solutions applied in this article is considered a kind standard boundary-type meshless numerical method, which only needs to set the points on the boundary for simulations directly and which can obtain the high accurate results. Therefore, the problems of mesh generation and related procedures of numerical quadrature can be easy avoided. However, the users were usually troubled by determiniding the suiteble space position of the source points when MFS was implemented. On the other words, the numerical results are sensitive with the space position of source points for MFS implements. In order to search suiteble space position of the source points, the algorithm of particle swarm optimization is applied in this article. The so-called PSO can be considered a kind of metaheuristic algorithm without training procedures, which based on the optimization theory for particles group motion. Thus the proposed model based on the PSO-MFS, can be considered the simple, innovative and high accurate numerical method with efficiency. After the developments of the potential flow model is completed, the benchmark problems of are also tested directly and shows the ability of the proposed mothod. The proposed model in this artical is based on potential flow theory and considered mass conservation and momentum conservation at same time. It is convinced that the achievements based on this PSO-MFS can become the effective numerical tool for science and engineering.