In general, particle simulation methods such as the MPS(Moving Particle Simulation) or SPH(Smoothed Particle Hydrodynamics) methods have some serious drawbacks for pressure solutions. The pressure field shows spurious high fluctuations both temporally and spatially. It is well known that pressure fluctuation primarily occurs because of the numerical approximation of the partial differential operators. The MPS and SPH methods employ a pre-defined kernel function in the approximation of the gradient and Laplacian operators. Because this kernel function is constructed artificially, an accurate solution cannot be guaranteed, especially when the distribution of particles is irregular. In this paper, we propose a particle simulation method based on the moving least-square technique for solving the partial differential operators using a Taylor-series expansion. The developed method was applied to the hydro-static pressure and dam-broken problems to validate it.