A superconvergent local discontinuous Galerkin method for nonlinear two-point boundary-value problems

Research output: Contribution to journalArticle

4 Citations (Scopus)

Abstract

In this paper, we present and analyze a superconvergent and high order accurate local discontinuous Galerkin (LDG) method for nonlinear two-point boundary-value problems (BVPs) of the form u = f (t, u), which arise in a wide variety of engineering applications. We prove the L2 stability of the LDG scheme and optimal L2 error estimates for the solution and for the auxiliary variable that approximates the first-order derivative. The order of convergence is proved to be p + 1, when piecewise polynomials of degree at most p are used. Our numerical experiments demonstrate optimal rates of convergence. Moreover, we show that the derivatives of the LDG solutions are superconvergent with order p + 1 toward the derivatives of Gausss-Radau projections of the exact solutions. Finally, we prove that the LDG solutions are superconvergent with order p + 3/2 toward Gauss-Radau projections of the exact solutions. Our computational results indicate that the observed numerical superconvergence rate is p + 2. Our proofs are valid for arbitrary regular meshes using piecewise polynomials of degree p ≥ 1 and for the periodic, Dirichlet, and mixed boundary conditions. All proofs are valid under the hypotheses of the existence and uniqueness theorem for BVPs. Several numerical results are presented to validate the theoretical results.

Original languageEnglish (US)
Pages (from-to)697-718
Number of pages22
JournalNumerical Algorithms
Volume79
Issue number3
DOIs
StatePublished - Nov 1 2018

Fingerprint

Local Discontinuous Galerkin Method
Discontinuous Galerkin
Nonlinear Boundary Value Problems
Galerkin methods
Two-point Boundary Value Problem
Boundary value problems
Piecewise Polynomials
Derivatives
Derivative
Exact Solution
Polynomials
Projection
Valid
Optimal Rate of Convergence
Optimal Error Estimates
Superconvergence
Auxiliary Variables
Existence and Uniqueness Theorem
Mixed Boundary Conditions
Order of Convergence

Keywords

  • A priori error estimates
  • Gauss-Radau projections
  • L stability
  • Local discontinuous Galerkin method
  • Nonlinear two-point boundary-value problems
  • Superconvergence

ASJC Scopus subject areas

  • Applied Mathematics

Cite this

A superconvergent local discontinuous Galerkin method for nonlinear two-point boundary-value problems. / Baccouch, Mahboub.

In: Numerical Algorithms, Vol. 79, No. 3, 01.11.2018, p. 697-718.

Research output: Contribution to journalArticle

@article{a2cb7762adcb4fa19e62f58896cf6585,
title = "A superconvergent local discontinuous Galerkin method for nonlinear two-point boundary-value problems",
abstract = "In this paper, we present and analyze a superconvergent and high order accurate local discontinuous Galerkin (LDG) method for nonlinear two-point boundary-value problems (BVPs) of the form u″ = f (t, u), which arise in a wide variety of engineering applications. We prove the L2 stability of the LDG scheme and optimal L2 error estimates for the solution and for the auxiliary variable that approximates the first-order derivative. The order of convergence is proved to be p + 1, when piecewise polynomials of degree at most p are used. Our numerical experiments demonstrate optimal rates of convergence. Moreover, we show that the derivatives of the LDG solutions are superconvergent with order p + 1 toward the derivatives of Gausss-Radau projections of the exact solutions. Finally, we prove that the LDG solutions are superconvergent with order p + 3/2 toward Gauss-Radau projections of the exact solutions. Our computational results indicate that the observed numerical superconvergence rate is p + 2. Our proofs are valid for arbitrary regular meshes using piecewise polynomials of degree p ≥ 1 and for the periodic, Dirichlet, and mixed boundary conditions. All proofs are valid under the hypotheses of the existence and uniqueness theorem for BVPs. Several numerical results are presented to validate the theoretical results.",
keywords = "A priori error estimates, Gauss-Radau projections, L stability, Local discontinuous Galerkin method, Nonlinear two-point boundary-value problems, Superconvergence",
author = "Mahboub Baccouch",
year = "2018",
month = "11",
day = "1",
doi = "10.1007/s11075-017-0456-0",
language = "English (US)",
volume = "79",
pages = "697--718",
journal = "Numerical Algorithms",
issn = "1017-1398",
publisher = "Springer Netherlands",
number = "3",

}

TY - JOUR

T1 - A superconvergent local discontinuous Galerkin method for nonlinear two-point boundary-value problems

AU - Baccouch, Mahboub

PY - 2018/11/1

Y1 - 2018/11/1

N2 - In this paper, we present and analyze a superconvergent and high order accurate local discontinuous Galerkin (LDG) method for nonlinear two-point boundary-value problems (BVPs) of the form u″ = f (t, u), which arise in a wide variety of engineering applications. We prove the L2 stability of the LDG scheme and optimal L2 error estimates for the solution and for the auxiliary variable that approximates the first-order derivative. The order of convergence is proved to be p + 1, when piecewise polynomials of degree at most p are used. Our numerical experiments demonstrate optimal rates of convergence. Moreover, we show that the derivatives of the LDG solutions are superconvergent with order p + 1 toward the derivatives of Gausss-Radau projections of the exact solutions. Finally, we prove that the LDG solutions are superconvergent with order p + 3/2 toward Gauss-Radau projections of the exact solutions. Our computational results indicate that the observed numerical superconvergence rate is p + 2. Our proofs are valid for arbitrary regular meshes using piecewise polynomials of degree p ≥ 1 and for the periodic, Dirichlet, and mixed boundary conditions. All proofs are valid under the hypotheses of the existence and uniqueness theorem for BVPs. Several numerical results are presented to validate the theoretical results.

AB - In this paper, we present and analyze a superconvergent and high order accurate local discontinuous Galerkin (LDG) method for nonlinear two-point boundary-value problems (BVPs) of the form u″ = f (t, u), which arise in a wide variety of engineering applications. We prove the L2 stability of the LDG scheme and optimal L2 error estimates for the solution and for the auxiliary variable that approximates the first-order derivative. The order of convergence is proved to be p + 1, when piecewise polynomials of degree at most p are used. Our numerical experiments demonstrate optimal rates of convergence. Moreover, we show that the derivatives of the LDG solutions are superconvergent with order p + 1 toward the derivatives of Gausss-Radau projections of the exact solutions. Finally, we prove that the LDG solutions are superconvergent with order p + 3/2 toward Gauss-Radau projections of the exact solutions. Our computational results indicate that the observed numerical superconvergence rate is p + 2. Our proofs are valid for arbitrary regular meshes using piecewise polynomials of degree p ≥ 1 and for the periodic, Dirichlet, and mixed boundary conditions. All proofs are valid under the hypotheses of the existence and uniqueness theorem for BVPs. Several numerical results are presented to validate the theoretical results.

KW - A priori error estimates

KW - Gauss-Radau projections

KW - L stability

KW - Local discontinuous Galerkin method

KW - Nonlinear two-point boundary-value problems

KW - Superconvergence

UR - http://www.scopus.com/inward/record.url?scp=85037632607&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85037632607&partnerID=8YFLogxK

U2 - 10.1007/s11075-017-0456-0

DO - 10.1007/s11075-017-0456-0

M3 - Article

AN - SCOPUS:85037632607

VL - 79

SP - 697

EP - 718

JO - Numerical Algorithms

JF - Numerical Algorithms

SN - 1017-1398

IS - 3

ER -