Remediation of contaminated groundwater by meshless local weak forms

In this paper, four meshless local weak form methods such as direct meshless local Petrov-Galerkin method (DMLPG), meshless local Petrov-Galerkin method (MLPG), local weak radial point interpolation method (LWRPIM) and local weak moving kriging method (LWMK) are applied to find the numerical solution of coupled non-linear advection-diffusion-reaction system arising in the prevention of groundwater contamination problem. A comparison between these methods is done from the perspective of accuracy and computational efficiency. An efficient fourth-order exponential time differencing Runge-Kutta method is utilized for the time discretization. The computational efficiency is the most significant advantage of the DMLPG method in comparison with the other meshless local weak form methods, because DMLPG reduces the computational costs, significantly. This is due to the fact that DMLPG shifts the numerical integrations over low-degree polynomials rather than over complicated shape functions. The main aim of this paper is to show that the meshless local weak form methods can be used for solving the system of non-linear partial differential equations especially coupled non-linear advection-diffusion-reaction system. The numerical results confirm the good efficiency of the proposed methods for solving our model.