Document Type : Research Paper
Department of Applied Mathematics, Imam Khomeini International University, Qazvin, 34149-16818, Iran
In this article, the meshless local radial point interpolation (MLRPI) methods are applied to simulate two dimensional wave equation subject to given appropriate initial and Neumann’s boundary conditions. In MLRPI method, all integrations are carried out locally over small quadrature domains of regular shapes such as square or circle. The radial point interpolation method is proposed to construct shape functions for MLRPI. A weak formulation with a Heaviside step function transforms the set of governing equations into local integral equations on local sub domains where Neumann’s boundary condition is imposed naturally. A two-step time discretization method with the help of Crank-Nicolson technique is employed to approximate the time derivatives. Convergence studies in the numerical example show that MLRPI method possesses excellent rates of convergence.