Saturation estimates for hp -finite element methods

In this paper we will prove saturation estimates for the adaptive $$hp$$ h p -finite element method for linear, second order partial differential equations. More specifically we will consider a sequence of nested finite element discretizations where we allow for both, local mesh refinement and locally increasing the polynomial order. We will prove that the energy norm of the error on the finer level can be estimated by the sum of a contraction of the old error and data oscillations. We will derive estimates of the contraction factor which are explicit with respect to the local mesh width and the local polynomial degree. In order to cover $$p$$ p -refinement of finite element spaces new polynomial projection operators will be introduced and new polynomial inverse estimates will be derive

On p-Robust Saturation for hp-AFEM

We consider the standard adaptive finite element loop SOLVE, ESTIMATE, MARK, REFINE, with ESTIMATE being implemented using the $p$-robust equilibrated flux estimator, and MARK being D\"orfler marking. As a refinement strategy we employ $p$-refinement. We investigate the question by which amount the local polynomial degree on any marked patch has to be increase in order to achieve a $p$-independent error reduction. The resulting adaptive method can be turned into an instance optimal $hp$-adaptive method by the addition of a coarsening routine

Adaptive boundary element methods with convergence rates

This paper presents adaptive boundary element methods for positive, negative, as well as zero order operator equations, together with proofs that they converge at certain rates. The convergence rates are quasi-optimal in a certain sense under mild assumptions that are analogous to what is typically assumed in the theory of adaptive finite element methods. In particular, no saturation-type assumption is used. The main ingredients of the proof that constitute new findings are some results on a posteriori error estimates for boundary element methods, and an inverse-type inequality involving boundary integral operators on locally refined finite element spaces.Comment: 48 pages. A journal version. The previous version (v3) is a bit lengthie

Space-Time Isogeometric Analysis of Parabolic Evolution Equations

We present and analyze a new stable space-time Isogeometric Analysis (IgA) method for the numerical solution of parabolic evolution equations in fixed and moving spatial computational domains. The discrete bilinear form is elliptic on the IgA space with respect to a discrete energy norm. This property together with a corresponding boundedness property, consistency and approximation results for the IgA spaces yields an a priori discretization error estimate with respect to the discrete norm. The theoretical results are confirmed by several numerical experiments with low- and high-order IgA spaces

A saturation property for the spectral-Galerkin approximation of a Dirichlet problem in a square

Both practice and analysis of adaptive $p$-FEMs and $hp$-FEMs raise the question what increment in the current polynomial degree $p$ guarantees a $p$-independent reduction of the Galerkin error. We answer this question for the $p$-FEM in the simplified context of homogeneous Dirichlet problems for the Poisson equation in the two dimensional unit square with polynomial data of degree $p$. We show that an increment proportional to $p$ yields a $p$-robust error reduction and provide computational evidence that a constant increment does not