Enforcing Neumann Boundary Conditions with Polynomial Extension Operators to Acheive Optimal Convergence Rates on Polytopial Meshes in the Finite Element Method

TL;DR

This paper develops new analysis techniques to prove optimal error estimates for Neumann boundary condition approximations in finite element methods on polytopial meshes, improving upon prior work that lacked $L^2$ optimality.

Contribution

It introduces novel analysis that establishes optimal error estimates for Neumann boundary conditions in both $W^1_ abla$ and $L^2$ norms, completing previous results.

Findings

Proves optimal $W^1_ abla$ error estimates for Neumann problems.

Establishes optimal $L^2$ error estimates for Neumann problems.

Completes the theoretical framework for boundary condition enforcement in finite element methods.

Abstract

In \cite{cheung2019optimally}, the authors presented two finite element methods for approximating second order boundary value problems on polytopial meshes with optimal accuracy without having to utilize curvilinear mappings. This was done by enforcing the boundary conditions through judiciously chosen polynomial extension operators. The $H^1$ error estimates were proven to be optimal for the solutions of both the Dirichlet and Neumann boundary value problems. It was also proven that the Dirichlet problem approximation converges optimally in $L^2$. However, optimality of the Neumann approximation in the $L^2$ norm was left as an open problem. In this work, we seek to close this problem by presenting new analysis that proves optimal error estimates for the Neumann approximation in the $W^1_\infty$ and $L^2$ norms.