<p>We consider a discontinuous Galerkin method for the numerical solution of boundary value problems in two-dimensional domains with curved boundaries. A key challenge in this setting is the potential loss of convergence order due to approximating the physical domain by a polygonal mesh. Unless boundary conditions can be accurately transferred from the true boundary to the computational one, such geometric approximation errors generally lead to suboptimal convergence. To overcome this limitation, a higher-order strategy based on polynomial reconstruction of boundary data was introduced for classical finite element methods in [<CitationRef CitationID="CR31">31</CitationRef>, <CitationRef CitationID="CR32">32</CitationRef>] and in the finite volume context in [<CitationRef CitationID="CR8">8</CitationRef>, <CitationRef CitationID="CR14">14</CitationRef>]. More recently, this approach was extended to discontinuous Galerkin methods in [<CitationRef CitationID="CR35">35</CitationRef>], leading to the DG–ROD method, which restores optimal convergence rates on polygonal approximations of domains with curved boundaries. In this work, we provide a rigorous theoretical analysis of the DG–ROD method, establishing existence and uniqueness of the discrete solution and deriving error estimates for a two-dimensional linear advection-diffusion-reaction problem with homogeneous Dirichlet boundary conditions on both convex and non-convex domains. Following and extending techniques from classical finite element methods [<CitationRef CitationID="CR32">32</CitationRef>], we prove that, under suitable regularity assumptions on the exact solution, the DG–ROD method achieves optimal convergence despite polygonal approximations. Finally, we illustrate and confirm the theoretical results with a numerical benchmark considering triangular meshes.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Optimal Error Estimates for a Discontinuous Galerkin Method on Curved Boundaries with Polygonal Meshes

  • Adérito Araújo,
  • Milene Santos

摘要

We consider a discontinuous Galerkin method for the numerical solution of boundary value problems in two-dimensional domains with curved boundaries. A key challenge in this setting is the potential loss of convergence order due to approximating the physical domain by a polygonal mesh. Unless boundary conditions can be accurately transferred from the true boundary to the computational one, such geometric approximation errors generally lead to suboptimal convergence. To overcome this limitation, a higher-order strategy based on polynomial reconstruction of boundary data was introduced for classical finite element methods in [31, 32] and in the finite volume context in [8, 14]. More recently, this approach was extended to discontinuous Galerkin methods in [35], leading to the DG–ROD method, which restores optimal convergence rates on polygonal approximations of domains with curved boundaries. In this work, we provide a rigorous theoretical analysis of the DG–ROD method, establishing existence and uniqueness of the discrete solution and deriving error estimates for a two-dimensional linear advection-diffusion-reaction problem with homogeneous Dirichlet boundary conditions on both convex and non-convex domains. Following and extending techniques from classical finite element methods [32], we prove that, under suitable regularity assumptions on the exact solution, the DG–ROD method achieves optimal convergence despite polygonal approximations. Finally, we illustrate and confirm the theoretical results with a numerical benchmark considering triangular meshes.