The local discontinuous galerkin method for the fourth-order euler-bernoulli partial differential equation in one space dimension. Part II: A posteriori error estimation

Mahboub Baccouch

Research output: Contribution to journalArticlepeer-review

21 Scopus citations

Abstract

In this paper new a posteriori error estimates for the local discontinuous Galerkin (LDG) method for one-dimensional fourth-order Euler-Bernoulli partial differential equation are presented and analyzed. These error estimates are computationally simple and are obtained by solving a local steady problem with no boundary condition on each element. We use the optimal error estimates and the superconvergence results proved in Part I to show that the significant parts of the discretization errors for the LDG solution and its spatial derivatives (up to third order) are proportional to (k+1) -degree Radau polynomials, when polynomials of total degree not exceeding k are used. These results allow us to prove that the k -degree LDG solution and its derivatives are O(h k+3/2) superconvergent at the roots of (k+1) -degree Radau polynomials. We use these results to construct asymptotically exact a posteriori error estimates. We further apply the results proved in Part I to prove that, for smooth solutions, these a posteriori LDG error estimates for the solution and its spatial derivatives at a fixed time tconverge to the true errors at O(hk+5/4 rate. We also prove that the global effectivity indices, for the solution and its derivatives up to third order, in the L2 -norm converge to unity at O(h1/2) rate. Our proofs are valid for arbitrary regular meshes and for Pk polynomials with k ≥ 1, and for periodic and other classical mixed boundary conditions. Our computational results indicate that the observed numerical convergence rates are higher than the theoretical rates. Finally, we present a local adaptive procedure that makes use of our local a posteriori error estimate.

Original languageEnglish (US)
Pages (from-to)1-34
Number of pages34
JournalJournal of Scientific Computing
Volume60
Issue number1
DOIs
StatePublished - Jul 2014

Keywords

  • A posteriori error estimates
  • Adaptive mesh method
  • Fourth-order Euler-Bernoulli equation
  • Local discontinuous Galerkin method
  • Radau points
  • Superconvergence

ASJC Scopus subject areas

  • Software
  • Theoretical Computer Science
  • Numerical Analysis
  • General Engineering
  • Computational Theory and Mathematics
  • Computational Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'The local discontinuous galerkin method for the fourth-order euler-bernoulli partial differential equation in one space dimension. Part II: A posteriori error estimation'. Together they form a unique fingerprint.

Cite this