A new method to generate a quasi-equal distribution of particles on arbitrary domains is provided. Based on an advancing front mesh generation method, triangles are formed at the boundary. These triangle center-points are used for the boundary condition of the Delaunay mesh generation method which generates inward triangles. All vertices of triangles are the particles in meshfree methods. An adaptive method is given to readjust distribution of particles, according to the gradient of fluid density.