Generalization of robustness test procedure for error estimators. Part II: test results for error estimators using SPR and REP

In this part of the paper we shall use the formulation given in the first part to assess the quality of recovery-based error estimators using two recovery methods, i.e. superconvergent patch recovery (SPR) and recovery by equilibrium in patches (REP). The recovery methods have been shown to be asymptotically robust and superconvergent when applied to two-dimensional problems. In this study we shall examine the behaviour of the recovery methods on several three-dimensional mesh patterns for patches located either inside or at boundaries. This is performed by first finding an asymptotic finite element solution, irrespective of boundary conditions at far ends of the domain, and then applying the recovery methods. The test procedure near kinked boundaries is explained in a step-by-step manner. The results are given in a series of tables and figures for various cases of three-dimensional mesh patterns. It has been experienced that the full superconvergent property is generally lost due to presence of boundary layer solution and the definition of the recoveries near boundaries though the results of the robustness test is still within an acceptable range