English

The local discontinuous Galerkin method for a singularly perturbed convection-diffusion problem with characteristic and exponential layers

Numerical Analysis 2022-09-22 v1 Numerical Analysis

Abstract

A singularly perturbed convection-diffusion problem,posed on the unit square in R2\mathbb{R}^2, is studied; its solution has both exponential and characteristic boundary layers. The problem is solved numerically using the local discontinuous Galerkin (LDG) method on Shishkin meshes. Using tensor-product piecewise polynomials of degree at most k>0k>0 in each variable, the error between the LDG solution and the true solution is proved to converge, uniformly in the singular perturbation parameter, at a rate of O((N1lnN)k+1/2)O((N^{-1}\ln N)^{k+1/2}) in an associated energy norm, where NN is the number of mesh intervals in each coordinate direction.(This is the first uniform convergence result proved for the LDG method applied to a problem with characteristic boundary layers.) Furthermore, we prove that this order of convergence increases to O((N1lnN)k+1)O((N^{-1}\ln N)^{k+1}) when one measures the energy-norm difference between the LDG solution and a local Gauss-Radau projection of the true solution into the finite element space.This uniform supercloseness property implies an optimal L2L^2 error estimate of order (N1lnN)k+1(N^{-1}\ln N)^{k+1} for our LDG method. Numerical experiments show the sharpness of our theoretical results.

Keywords

Cite

@article{arxiv.2209.10143,
  title  = {The local discontinuous Galerkin method for a singularly perturbed convection-diffusion problem with characteristic and exponential layers},
  author = {Yao Cheng and Martin Stynes},
  journal= {arXiv preprint arXiv:2209.10143},
  year   = {2022}
}

Comments

33pages, 7 figures, 7 tables