A new approach of superconvergence analysis of a low order nonconforming MFEM for reaction–diffusion equation

Houchao Zhang,An Wang

DOI: https://doi.org/10.1186/s13661-018-1091-y

2018-11-14

Boundary Value Problems

Abstract:In this paper, a low order nonconforming mixed finite element method (MFEM) is studied with EQ1rotdocumentclass[12pt]{minimal}usepackage{amsmath}usepackage{wasysym}usepackage{amsfonts}usepackage{amssymb}usepackage{amsbsy}usepackage{mathrsfs}usepackage{upgreek}setlength{oddsidemargin}{-69pt}egin{document}$EQ_{1}^{mathrm{rot}}$end{document} element and zero order Raviart–Thomas (R–T) element for a class of nonlinear reaction–diffusion equations. On the one hand, a priori bound is proved using Lyapunov functional, which leads to the global existence and uniqueness of the approximation solutions. Further, the superclose estimates of order O(h2)documentclass[12pt]{minimal}usepackage{amsmath}usepackage{wasysym}usepackage{amsfonts}usepackage{amssymb}usepackage{amsbsy}usepackage{mathrsfs}usepackage{upgreek}setlength{oddsidemargin}{-69pt}egin{document}$O(h^{2})$end{document} for original variable u in broken H1documentclass[12pt]{minimal}usepackage{amsmath}usepackage{wasysym}usepackage{amsfonts}usepackage{amssymb}usepackage{amsbsy}usepackage{mathrsfs}usepackage{upgreek}setlength{oddsidemargin}{-69pt}egin{document}$H^{1}$end{document} norm and the flux p=∇udocumentclass[12pt]{minimal}usepackage{amsmath}usepackage{wasysym}usepackage{amsfonts}usepackage{amssymb}usepackage{amsbsy}usepackage{mathrsfs}usepackage{upgreek}setlength{oddsidemargin}{-69pt}egin{document}$oldsymbol{p}= abla u$end{document} in (L2)2documentclass[12pt]{minimal}usepackage{amsmath}usepackage{wasysym}usepackage{amsfonts}usepackage{amssymb}usepackage{amsbsy}usepackage{mathrsfs}usepackage{upgreek}setlength{oddsidemargin}{-69pt}egin{document}$(L^{2})^{2}$end{document} norm are obtained respectively for a semi-discrete scheme. On the other hand, a linearized Crank–Nicolson fully-discrete scheme is established and the superclose estimates of order O(h2+τ2)documentclass[12pt]{minimal}usepackage{amsmath}usepackage{wasysym}usepackage{amsfonts}usepackage{amssymb}usepackage{amsbsy}usepackage{mathrsfs}usepackage{upgreek}setlength{oddsidemargin}{-69pt}egin{document}$O(h^{2}+ au^{2})$end{document} are also obtained unconditionally by making full use of the special characters of the above mixed finite elements (MFEs) and two different splitting techniques, which are used to deal with the consistency errors and nonlinear terms, respectively. These approaches circumvent the restrictive condition on a time step size arising as an inverse inequality used to prove a posteriori bounds in L∞documentclass[12pt]{minimal}usepackage{amsmath}usepackage{wasysym}usepackage{amsfonts}usepackage{amssymb}usepackage{amsbsy}usepackage{mathrsfs}usepackage{upgreek}setlength{oddsidemargin}{-69pt}egin{document}$L^{infty}$end{document} norm, which is necessary for nonlinear problems for the conventional finite element analysis. What is more, the corresponding global superconvergent results are derived through interpolated postprocessing techniques. Finally, numerical results are provided to confirm the theoretical analysis. Here h is the subdivision parameter and τ is the time step.

mathematics, applied