A posteriori error estimates for nonconforming discretizations of singularly perturbed biharmonic operators

Dietmar Gallistl,Shudan Tian
2024-10-17
Abstract:For the pure biharmonic equation and a biharmonic singular perturbation problem, a residual-based error estimator is introduced which applies to many existing nonconforming finite elements. The error estimator involves the local best-approximation error of the finite element function by piecewise polynomial functions of the degree determining the expected approximation order, which need not coincide with the maximal polynomial degree of the element, for example if bubble functions are used. The error estimator is shown to be reliable and locally efficient up to this polynomial best-approximation error and oscillations of the right-hand side.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to provide a reliable and locally effective a posteriori error estimation method for the non - conforming discretization of the singularly perturbed biharmonic operator. Specifically, the article mainly focuses on the following two aspects: 1. **Pure biharmonic equation**: \[ \Delta^2 u = f \quad \text{in} \quad \Omega \] where \( u=\frac{\partial u}{\partial n} = 0\) on the boundary \( \partial\Omega\). 2. **Fourth - order singularly perturbed problem**: \[ \varepsilon^2\Delta^2 u-\Delta u = f\quad \text{in} \quad \Omega \] also with the clamped boundary condition \( u = \frac{\partial u}{\partial n}=0\) on \( \partial\Omega\), where \( \varepsilon\) is a small positive parameter. ### Main challenges For singularly perturbed problems (i.e., when \( \alpha = 1\)), some non - conforming finite element methods (such as Morley elements) perform poorly and may even fail to converge. Therefore, this paper aims to improve the performance of these methods by introducing new residual - based a posteriori error estimators and ensure their reliability and local effectiveness. ### Method overview The authors propose a new error estimator that contains a new term to deal with the weak continuity problem caused by high - order polynomial shape functions. Specifically, for each simplex \( T\), the error estimator contains the following term: \[ |||u_h-\Pi_\ell u_h|||_T \] where \( \Pi_\ell\) is the \( L^2\) projection onto polynomials of degree not exceeding \( \ell\), and \( \ell\) is the desired approximation order. ### Key contributions - **Reliability**: It is proved that the proposed error estimator is reliable and can accurately estimate the error. - **Local effectiveness**: It is shown that the error estimator is locally effective up to the best piecewise polynomial approximation error and the oscillation of the right - hand - side data. - **Wide applicability**: This method is applicable to many existing non - conforming finite element methods, including Morley elements, Specht elements, etc. Through these improvements, this paper provides a more robust and accurate a posteriori error estimation framework for singularly perturbed biharmonic problems.