Abstract
The paper presents a new analytical method called the local fractional Laplace variational iteration method (LFLVIM), which is a combination of the local fractional Laplace transform (LFLT) and the local fractional variational iteration method (LFVIM), for solving the two-dimensional Helmholtz and coupled Helmholtz equations with local fractional derivative operators (LFDOs). The operators are taken in the local fractional sense. Two test problems are presented to demonstrate the efficiency and the accuracy of the proposed method. The approximate solutions obtained are compared with the results obtained by the local fractional Laplace decomposition method (LFLDM). The results reveal that the LFLVIM is very effective and convenient to solve linear and nonlinear PDEs.