A Nyström-based finite element method on polygonal elements

Akash Anand,Jeffrey S. Ovall,Steffen Weisser

DOI: https://doi.org/10.1016/j.camwa.2018.03.007

2017-08-24

Abstract:We consider families of finite elements on polygonal meshes, that are defined implicitly on each mesh cell as solutions of local Poisson problems with polynomial data. Functions in the local space on each mesh cell are evaluated via Nyström discretizations of associated integral equations, allowing for curvilinear polygons and non-polynomial boundary data. Several experiments demonstrate the approximation quality of interpolated functions in these spaces.

Numerical Analysis

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to develop a finite element method based on Nyström discretization to solve second - order linear elliptic boundary value problems on polygonal meshes. Specifically, the author aims to evaluate basis functions and their derivatives by using the Nyström method, thereby forming a finite - element linear system. This method allows the handling of polygonal elements with curved edges and can better resolve singular behavior in the local space.

### Problem Background

Traditional finite element methods usually rely on elements of simple geometric shapes such as triangles, quadrilaterals or hexahedra. However, in high - precision calculations, these simple elements may limit the flexibility and adaptability of the mesh. Therefore, in recent years, more and more research has focused on developing flexible finite - element discretization schemes that can be used on polygonal and polyhedral meshes. These schemes not only simplify the mesh refinement and coarsening processes but also can handle complex geometric shapes more naturally.

### Core Contributions of the Paper

1. **Application of the Nyström Method**: The paper proposes a method based on Nyström discretization to solve the second - kind integral equations associated with each mesh element. This is different from the traditional boundary element method (BEM - FEM), which uses the first - kind integral equations. The Nyström method is simpler when dealing with high - order discretization, can better resolve singular behavior in the local space, and allows the handling of polygonal elements with curved edges.

2. **Verification by Numerical Experiments**: The author verifies the effectiveness of this method through a series of numerical experiments and shows the approximation quality of the interpolation function on different types of polygonal meshes.

3. **Theoretical Analysis**: The paper discusses in detail the definitions of local and global discrete approximation spaces and provides a theoretical framework on how to implicitly express functions on each mesh element.

### Mathematical Formulas

The key formulas involved in the paper include:

- **Definition of the Poisson Problem**:
\[
v\in V_{m}(K)\text{ if and only if }\Delta v\in P_{m - 2}(K)\text{ inside }K\text{, and }v\in P_{m}(\partial K)\text{ on }\partial K
\]

- **Integral Equation of Nyström Discretization**:
\[
w(x)=\int_{\partial K}\frac{\partial G(x,y)}{\partial n(y)}\varphi(y)ds(y)\quad\text{for }x\in K
\]
where \(\varphi\) satisfies the second - kind integral equation:
\[
\frac{\varphi(x)}{2}-\int_{\partial K}\frac{\partial G(x,y)}{\partial n(y)}\varphi(y)ds(y)= - g(x)\quad\text{for }x\in\partial K
\]

- **Treatment of Dirichlet Boundary Conditions**:
\[
\left(1-\frac{\alpha}{2}\right)\varphi(z)-\int_{\partial K}\frac{\partial G(z,y)}{\partial n(y)}\varphi(y)ds(y)= - g(z)
\]
where \(z\) is a vertex on the boundary with an internal angle of \(\alpha\pi\).

### Summary

By introducing the Nyström method, this paper provides a new finite - element method suitable for second - order linear elliptic boundary value problems on polygonal meshes. This method not only improves computational efficiency but also can better handle complex geometric shapes and singular behavior.