A Posteriori Error Estimation for the Discontinuous Galerkin Method Applied to Nonlinear Volterra Integro-differential Equations

Abstract

In this paper, we present and rigorously analyze a novel residual-based a posteriori error estimator for the discontinuous Galerkin (DG) method applied to nonlinear Volterra integro-differential equations. We first prove that the DG solution using piecewise polynomials of degree

p⩾1

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p \geqslant 1$$\end{document}

achieves the optimal convergence rate of

\mathcal{O}(h^{p+1})

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathcal {O}(h^{p+1})$$\end{document}

in the

L^2

\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. Moreover, we show that the DG solution exhibits

\mathcal{O}(h^{2p+1})

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathcal {O}(h^{2p+1})$$\end{document}

superconvergence at downwind points. These results are used to show that the DG solution is

\mathcal{O}(h^{p+2})

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathcal {O}(h^{p+2})$$\end{document}

superconvergent towards a specific Gauss-Radau projection of the exact solution. This superconvergence result enables a decomposition of the discretization error into two components: a dominant term proportional to a right Radau polynomial of degree

(p+1)

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$(p+1)$$\end{document}

and a remainder that converges at a rate of

\mathcal{O}(h^{p+2})

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathcal {O}(h^{p+2})$$\end{document}

in the

L^2

\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. The global superconvergence properties are used to construct a residual-based a posteriori error estimator. The proposed estimator is computationally efficient, easy to implement, and asymptotically exact. It is constructed by solving local residual problems without imposing initial conditions on each element. We establish convergence of the a posteriori error estimator to the true error at the optimal rate of

\mathcal{O}(h^{p+2})

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathcal {O}(h^{p+2})$$\end{document}

under mesh refinement. Finally, we prove that the global effectivity index of the estimator approaches unity at a rate of

\mathcal{O}(h)

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\mathcal {O}(h)$$\end{document}

. Several numerical examples are provided to validate the theoretical results and to illustrate the accuracy and reliability of the proposed estimator.

Keywords

Discontinuous Galerkin (DG) method / Nonlinear Volterra integro-differential equations / Superconvergence / Right Radau polynomial / A posteriori error estimation / 45D05 / 65R20 / 65L60

References

[1] Abramowitz M, Stegun IA. Handbook of Mathematical Functions, 1965, New York, Dover

[2] Adjerid S, Baccouch M. Asymptotically exact a posteriori error estimates for a one-dimensional linear hyperbolic problem. Appl. Numer. Math., 2010, 60: 903-914

[3] Ainsworth M, Oden JT. A Posteriori Error Estimation in Finite Element Analysis, 2000, New York, Wiley

[4] Baccouch M. Analysis of a posteriori error estimates of the discontinuous Galerkin method for nonlinear ordinary differential equations. Appl. Numer. Math., 2016, 106: 129-153

[5] Brunner H. Collocation Methods for Volterra Integral and Related Functional Differential Equations, 2004, Cambridge, Cambridge University Press 15

[6] Brunner, H., Schötzau, D.: hp-discontinuous Galerkin time-stepping for Volterra integrodifferential equations. SIAM J. Numer. Anal. 44(1), 224–245 (2006)

[7] Ciarlet PG. The Finite Element Method for Elliptic Problems, 1978, Amsterdam, North-Holland Pub. Co.

[8] Delfour M, Dubeau F. Discontinuous polynomial approximations in the theory of one-step, hybrid and multistep methods for nonlinear ordinary differential equations. Math. Comput., 1986, 47: 169-189

[9] Delfour M, Hager W, Trochu F. Discontinuous Galerkin methods for ordinary differential equation. Math. Comput., 1981, 154: 455-473

[10] Estep D. A posteriori error bounds and global error control for approximation of ordinary differential equations. SIAM J. Numer. Anal., 1995, 32: 1-48

[11] Jiang Y-J. On spectral methods for Volterra-type integro-differential equations. J. Comput. Appl. Math., 2009, 230(2): 333-340

[12] Johnson C. Error estimates and adaptive time-step control for a class of one-step methods for stiff ordinary differential equations. SIAM J. Numer. Anal., 1988, 25: 908-926

[13] Lakshmikantham V. Theory of Integro-differential Equations, 1995, Boca Raton, CRC Press1

[14] LeSaint, P., Raviart, P.A.: Mathematical aspects of finite elements in partial differential equations. In: de Boor, C. (ed) On a Finite Element Method for Solving the Neutron Transport Equation, pp. 89–123. Academic Press, New York (1974)

[15] Mustapha K. A superconvergent discontinuous Galerkin method for Volterra integro-differential equations, smooth and non-smooth kernels. Math. Comput., 2013, 82(284): 1987-2005

[16] Mustapha K, Ryan JK. Post-processing discontinuous Galerkin solutions to Volterra integro-differential equations: analysis and simulations. J. Comput. Appl. Math., 2013, 253: 89-103

[17] Reed, W.H., Hill, T.R.: Triangular mesh methods for the neutron transport equation. Technical Report LA-UR-73-479, Los Alamos Scientific Laboratory, Los Alamos (1973)

[18] Tang, T., Xu, X., Cheng, J.: On spectral methods for Volterra integral equations and the convergence analysis. J. Comput. Math. 825–837 (2008)

[19] Verfürth R. A Review of A Posteriori Error Estimation and Adaptive Mesh-Refinement Techniques, 1996, New York, Wiley-Teubner

[20] Wang Z-Q, Guo Y-L, Yi L-J. An hp\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$hp$$\end{document}-version Legendre-Jacobi spectral collocation method for Volterra integro-differential equations with smooth and weakly singular kernels. Math. Comput., 2017, 86(307): 2285-2324

[21] Wei Y, Chen Y. Convergence analysis of the spectral methods for weakly singular Volterra integro-differential equations with smooth solutions. Adv. Appl. Math. Mech., 2012, 4(1): 1-20

[22] Wen J, Huang C. Multistep Runge-Kutta methods for Volterra integro-differential equations. J. Comput. Appl. Math., 2024, 436 ArticleID: 115384

[23] Wen J, Huang C, Li M. Stability analysis of Runge-Kutta methods for Volterra integro-differential equations. Appl. Numer. Math., 2019, 146: 73-88

[24] Yi, L.: An h-p version of the continuous Petrov-Galerkin finite element method for nonlinear Volterra integro-differential equations. J. Sci. Comput. 65(2), 715–734 (2015)

[25] Yi, L., Guo, B.: An h-p version of the continuous Petrov-Galerkin finite element method for Volterra integro-differential equations with smooth and nonsmooth kernels. SIAM J. Numer. Anal. 53(6), 2677–2704 (2015)

[26] Yi L, Zhang M, Mao X. Superconvergent postprocessing of the discontinuous Galerkin time stepping method for nonlinear Volterra integro-differential equations. J. Comput. Appl. Math., 2023, 427 ArticleID: 115140

[27] Zhang S, Lin T, Lin Y, Rao M. Defect correction and a posteriori error estimation of Petrov-Galerkin methods for nonlinear Volterra integro-differential equations. Appl. Math., 2000, 45: 241-263

[28] Zhang, S.-H., Lin, T., Lin, Y.-P., Rao, M.: Extrapolation and a-posteriori error estimators of Petrov-Galerkin methods for non-linear Volterra integro-differential equations. J. Comput. Math. 19(4), 407–422 (2001)

RIGHTS & PERMISSIONS

The Author(s)

Just Accepted

This article has successfully passed peer review and final editorial review, and will soon enter typesetting, proofreading and other publishing processes. The currently displayed version is the accepted final manuscript. The officially published version will be updated with format, DOI and citation information upon launch. We recommend that you pay attention to subsequent journal notifications and preferentially cite the officially published version. Thank you for your support and cooperation.