TY - JOUR
T1 - The local discontinuous galerkin method for the fourth-order euler-bernoulli partial differential equation in one space dimension. Part I
T2 - Superconvergence error analysis
AU - Baccouch, Mahboub
N1 - Funding Information:
Acknowledgments The author would also like to thank the anonymous referees for their constructive comments and remarks which helped improve the quality and readability of the paper. This research was partially supported by the NASA Nebraska Space Grant Program and UCRCA at the University of Nebraska at Omaha.
PY - 2014/6
Y1 - 2014/6
N2 - In this paper we develop and analyze a new superconvergent local discontinuous Galerkin (LDG) method for approximating solutions to the fourth-order Euler-Bernoulli beam equation in one space dimension. We prove the $$L 2$$ L 2 stability of the scheme and several optimal L 2 error estimates for the solution and for the three auxiliary variables that approximate derivatives of different orders. Our numerical experiments demonstrate optimal rates of convergence. We also prove superconvergence results towards particular projections of the exact solutions. More precisely, we prove that the LDG solution and its spatial derivatives (up to third order) are O (h k + 3 / 2) super close to particular projections of the exact solutions for k th-degree polynomial spaces while computational results show higher O (h k + 2) convergence rate. Our proofs are valid for arbitrary regular meshes and for P k polynomials with ≥ 1k 1, and for periodic, Dirichlet, and mixed boundary conditions. These superconvergence results will be used to construct asymptotically exact a posteriori error estimates by solving a local steady problem on each element. This will be reported in Part II of this work, where we will prove that the a posteriori LDG error estimates for the solution and its derivatives converge to the true errors in the L 2 -norm under mesh refinement.
AB - In this paper we develop and analyze a new superconvergent local discontinuous Galerkin (LDG) method for approximating solutions to the fourth-order Euler-Bernoulli beam equation in one space dimension. We prove the $$L 2$$ L 2 stability of the scheme and several optimal L 2 error estimates for the solution and for the three auxiliary variables that approximate derivatives of different orders. Our numerical experiments demonstrate optimal rates of convergence. We also prove superconvergence results towards particular projections of the exact solutions. More precisely, we prove that the LDG solution and its spatial derivatives (up to third order) are O (h k + 3 / 2) super close to particular projections of the exact solutions for k th-degree polynomial spaces while computational results show higher O (h k + 2) convergence rate. Our proofs are valid for arbitrary regular meshes and for P k polynomials with ≥ 1k 1, and for periodic, Dirichlet, and mixed boundary conditions. These superconvergence results will be used to construct asymptotically exact a posteriori error estimates by solving a local steady problem on each element. This will be reported in Part II of this work, where we will prove that the a posteriori LDG error estimates for the solution and its derivatives converge to the true errors in the L 2 -norm under mesh refinement.
KW - Fourth-order Euler-Bernoulli equation
KW - Local discontinuous Galerkin method
KW - Optimal error estimates
KW - Projection
KW - Stability
KW - Superconvergence
UR - http://www.scopus.com/inward/record.url?scp=84900807304&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84900807304&partnerID=8YFLogxK
U2 - 10.1007/s10915-013-9782-0
DO - 10.1007/s10915-013-9782-0
M3 - Article
AN - SCOPUS:84900807304
SN - 0885-7474
VL - 59
SP - 795
EP - 840
JO - Journal of Scientific Computing
JF - Journal of Scientific Computing
IS - 3
ER -