Abstract

For functions uin H^1(Omega ) in a bounded polytope Omega subset {mathbb {R}}^dd=1,2,3 with plane sides for d=2,3 which are Gevrey regular in overline{Omega }backslash {mathscr {S}} with point singularities concentrated at a set {mathscr {S}}subset overline{Omega } consisting of a finite number of points in overline{Omega }, we prove exponential rates of convergence of hp-version continuous Galerkin finite element methods on affine families of regular, simplicial meshes in Omega . The simplicial meshes are geometrically refined towards {mathscr {S}} but are otherwise unstructured.

Highlights

  • Many nonlinear PDEs admit solutions which are smooth in a bounded, polytopal domain ⊂ R, but exhibit isolated point singularities at a set S ⊂

  • Research performed in part while the authors were visiting the Erwin Schrödinger Institute (ESI) in Vienna, Austria, during the ESI thematic term “Numerical Analysis of Complex PDE Models in the Sciences” from June-August 2018

  • We show that the output of Algorithm 1 can, for each fixed 0 < σ < 1, be identified with sequences Mκ,σ of regular, simplicial meshes which are σ -geometrically refined towards S in

Read more

Summary

Introduction

Many nonlinear PDEs admit solutions which are smooth in a bounded, polytopal domain ⊂ R, but exhibit isolated point singularities at a set S ⊂. We prove an exponential convergence result for C0-conforming hp-FEM on regular, simplicial mesh families with isotropic, geometric refinement towards the singular point(s) c ∈ S. These meshes are in addition required to be shape-regular. The rate (1) coincides, in the cases d = 1, 2 and for analytic solutions, i.e. when δ = 1, with the exponential convergence rate bounds obtained in [18,19] for corner singularities on structured geometric meshes (consisting of axiparallel quadrilaterals with inserted triangles to remove irregular nodes).

Analytic regularity
Weighted Sobolev spaces
Examples of boundary value problems with Gevrey-regular solutions
Linear elliptic boundary value problems in polygons
Three-dimensional problems
Local polynomial spaces
Mesh layers
Statement of the exponential convergence result
Univariate hp-projectors and hp error bounds
Tensor projector on the unit cube
Transformation formula
Element interpolants
Exponential convergence in broken Sobolev norms
Enforcement of homogeneous Dirichlet boundary conditions
Concluding remarks
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call