文摘
Mathematical modeling of pattern formation in developmental biology leads to non-linear reaction–diffusion systems which are usually highly stiff in both diffusion and reaction terms. In this paper, the Direct Meshless Local Petrov–Galerkin (DMLPG) procedure is applied to find the numerical solution of some non-linear time-dependent reaction–diffusion systems such as Schnakenberg model, Gierer–Meinhardt model, FitzHugh–Nagumo model and Gray–Scott model. As far as we know, it is the first time that DMLPG method is applied for solving non-linear partial differential equations (PDEs) and systems of PDEs. Computational efficiency is the most significant advantage of the DMLPG method in comparison with the classic Meshless Local Petrov–Galerkin (MLPG) method. This is due to the fact that DMLPG shifts the numerical integrations over low-degree polynomials instead of over complicated moving least squares (MLS) shape functions and this reduces the computational costs, significantly. The main aim of this paper is to show that the DMLPG method is also suitable for solving the non-linear time-dependent systems, especially reaction–diffusion systems. Numerical results support the good efficiency of the proposed method for solving non-linear reaction–diffusion systems. Also it is shown that DMLPG provides considerable savings in computational time in comparison with the classical MLPG method.