A dual PetrovGalerkin finite element method for the convectiondiffusion equation
We present a minimum-residual finite element method (based on a dual PetrovGalerkin formulation) for convectiondiffusion problems in a higher order, adaptive, continuous Galerkin setting. The method borrows concepts from both the Discontinuous PetrovGalerkin (DPG) method by Demkowicz and Gopalakrishnan (2011) and the method of variational stabilization by Cohen, Dahmen, and Welper (2012), and it can also be interpreted as a variational multiscale method in which the fine-scales are defined through a dual-orthogonality condition. A key ingredient in the method is the proper choice of dual norm used to measure the residual, and we present two choices which are observed to be robust in both convection and diffusion-dominated regimes, as well as a proof of stability for quasi-uniform meshes and a method for the weak imposition of boundary conditions. Numerically obtained convergence rates in 2D are reported, and benchmark numerical examples are given to illustrate the behavior of the method.
petrovgalerkin convectiondiffusion high order stabilized finite elements