A Unified Equilibrated Flux Recovery Framework with Robust A Posteriori Error Estimation

We introduce the Equilibrated Averaging Residual Method (EARM), a unified equilibrated flux-recovery framework for elliptic interface problems that applies to a broad class of finite element discretizations. The method is applicable in both two and three dimensions and for arbitrary polynomial orders, and it enables the construction of computationally efficient recovered fluxes. We develop EARM for both discontinuous Galerkin (DG) and conforming finite element discretizations. For DG methods, EARM can be applied directly and yields an explicit recovered flux that coincides with state-of-the-art conservative flux reconstructions. For conforming discretizations, we further propose the Orthogonal Null-space--Eliminated EARM (ON-EARM), which ensures uniqueness by restricting the correction flux to the orthogonal complement of the divergence-free null space. We prove local conservation and establish a robust a~posteriori error estimator for the recovered flux in two dimensions, with robustness measured with respect to jumps in the diffusion coefficient. Numerical results in two and three dimensions confirm the theoretical findings.