$hp$-a posteriori error estimates for hybrid high-order methods applied to biharmonic problems

TL;DR

This paper develops residual-based $hp$-a posteriori error estimators for hybrid high-order methods solving biharmonic problems, providing bounds and demonstrating efficiency through numerical experiments.

Contribution

It introduces a novel residual-based $hp$-a posteriori error estimator for HHO methods on polytopal meshes for biharmonic equations, with new bounds and decomposition techniques.

Findings

The estimators effectively bound the error in biharmonic problem solutions.

Numerical experiments confirm the theoretical efficiency of the proposed estimators.

The approach handles both conforming and nonconforming error components.

Abstract

We derive a residual-based $hp$-a posteriori error estimator for hybrid high-order (HHO) methods on simplicial meshes applied to the biharmonic problem posed on two- and three-dimensional polytopal Lipschitz domains. The a posteriori error estimator hinges on an error decomposition into conforming and nonconforming components. To bound the nonconforming error, we use a $C^1$-partition of unity constructed via Alfeld splittings, combined with local Helmholtz decompositions on vertex stars. For the conforming error, we design two residual-based estimators, each associated with a specific interpolation operator. In the first setting, the upper bound for the conforming error involves only the stabilization term and the data oscillation. In the second setting, the bound additionally incorporates bulk residuals, normal flux jumps, and tangential jumps. Numerical experiments confirm the theoretical findings and demonstrate the efficiency of the proposed estimators.