Error estimation for non-linear singularly perturbed reaction-diffusion parabolic problems via element-free Galerkin method

The current study aims to develop an error analysis for non-linear parabolic singularly perturbed reaction-diffusion problems using element-free Galerkin method. A robust numerical methodology is introduced based on combining the implicit Crank-Nicolson scheme for temporal derivatives and the element-free Galerkin (EFG) method for spatial derivatives. The moving least-squares (MLS) approximation has been employed to generate the shape functions. Essential boundary conditions have been enforced by the incorporation of the Lagrange multiplier method. Due to the presence of steep boundary layers in the solution of the considered problem, a piecewise-uniform layer-adapted Shishkin's technique has been used to generate nodal points at the transition point. The stability and error analysis of the present method on a discrete L2-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L<^>{2}-$$\end{document}norm is analyzed in an innovative theoretical framework. The & varepsilon;\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon$$\end{document}-uniform convergency of the fully-discrete EFG method is shown to be O(tau 2+dsm)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {O}(\tau <^>{2}+d_{s}<^>{m})$$\end{document}, where tau\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau$$\end{document} and dsm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_{s}<^>{m}$$\end{document} are the time step size and size of the influence domain, respectively. The Lagrange multiplier method has been incorporated to deal with the implementation of essential boundary conditions. Lastly, a few numerical experiments are performed to validate the theoretical results and verify the computational consistency and robustness of the proposed scheme. The L infinity\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L_{\infty }$$\end{document} errors and the convergence rate have been presented.