Primal dual mixed finite element methods for indefinite advection--diffusion equations

We consider primal-dual mixed finite element methods for the advection--diffusion equation. For the primal variable we use standard continuous finite element space and for the flux we use the Raviart-Thomas space. We prove optimal a priori error estimates in the energy- and the $L^2$-norms for the primal variable in the low Peclet regime. In the high Peclet regime we also prove optimal error estimates for the primal variable in the $H(div)$ norm for smooth solutions. Numerically we observe that the method eliminates the spurious oscillations close to interior layers that pollute the solution of the standard Galerkin method when the local Peclet number is high. This method, however, does produce spurious solutions when outflow boundary layer presents. In the last section we propose two simple strategies to remove such numerical artefacts caused by the outflow boundary layer and validate them numerically.Comment: 25 pages, 6 figures, 5 table

Reliable a-posteriori error estimators for hphp-adaptive finite element approximations of eigenvalue/eigenvector problems

We present reliable a-posteriori error estimates for $hp$-adaptive finite element approximations of eigenvalue/eigenvector problems. Starting from our earlier work on $h$ adaptive finite element approximations we show a way to obtain reliable and efficient a-posteriori estimates in the $hp$-setting. At the core of our analysis is the reduction of the problem on the analysis of the associated boundary value problem. We start from the analysis of Wohlmuth and Melenk and combine this with our a-posteriori estimation framework to obtain eigenvalue/eigenvector approximation bounds.Comment: submitte

Analytic Regularity for Linear Elliptic Systems in Polygons and Polyhedra

We prove weighted anisotropic analytic estimates for solutions of second order elliptic boundary value problems in polyhedra. The weighted analytic classes which we use are the same as those introduced by Guo in 1993 in view of establishing exponential convergence for hp finite element methods in polyhedra. We first give a simple proof of the known weighted analytic regularity in a polygon, relying on a new formulation of elliptic a priori estimates in smooth domains with analytic control of derivatives. The technique is based on dyadic partitions near the corners. This technique can successfully be extended to polyhedra, providing isotropic analytic regularity. This is not optimal, because it does not take advantage of the full regularity along the edges. We combine it with a nested open set technique to obtain the desired three-dimensional anisotropic analytic regularity result. Our proofs are global and do not require the analysis of singular functions.Comment: 54 page

hp-version time domain boundary elements for the wave equation on quasi-uniform meshes

Solutions to the wave equation in the exterior of a polyhedral domain or a screen in $\mathbb{R}^3$ exhibit singular behavior from the edges and corners. We present quasi-optimal $hp$-explicit estimates for the approximation of the Dirichlet and Neumann traces of these solutions for uniform time steps and (globally) quasi-uniform meshes on the boundary. The results are applied to an $hp$-version of the time domain boundary element method. Numerical examples confirm the theoretical results for the Dirichlet problem both for screens and polyhedral domains.Comment: 41 pages, 11 figure