Abstract

We propose and analyze a high order unfitted mixed finite element method for the pseudostress-velocity formulation of the Stokes problem with Dirichlet boundary condition on a fluid domain Omega with curved boundary Gamma. The method consists of approximating Omega by a polyhedral subdomain D-h, with boundary Gamma(h), where a Galerkin method is applied to approximate the solution, and on a transferring technique, based on integrating the extrapolated discrete gradient of the velocity, to approximate the Dirichlet boundary data on the computational boundary Gamma(h). The associated Galerkin scheme is defined by Raviart-Thomas elements of order k >= 0 for the pseudostress and discontinuous polynomials of degree k for the velocity. Provided suitable hypotheses on the mesh near the boundary Gamma, we prove well-posedness of the Galerkin scheme by means of the BabuSka-Brezzi theory and establish the corresponding optimal convergence of O(h(k+1)). Moreover, for the case when Gamma(h) is constructed through a piecewise linear interpolation of Gamma, we propose a reliable and quasi-efficient residual-based a posteriori error estimator. Its definition makes use of a postprocessed velocity with enhanced accuracy to achieve the same rate of convergence of the method when the solution is smooth enough. Numerical experiments illustrate the performance of the scheme, show the behavior of the associated adaptive algorithm and validate the theory. (C) 2019 Elsevier B.V. All rights reserved.