Volume 32, pp. 1-16, 2008.

On the role of boundary conditions for CIP stabilization of higher order finite elements

Friedhelm Schieweck

Abstract

We investigate the Continuous Interior Penalty (CIP) stabilization method for higher order finite elements applied to a convection diffusion equation with a small diffusion parameter $\epsilon$. Performing numerical experiments, it turns out that strongly imposed Dirichlet boundary conditions lead to relatively bad numerical solutions. However, if the Dirichlet boundary conditions are imposed on the inflow part of the boundary in a weak sense and additionally on the whole boundary in an $\epsilon$-weighted weak sense due to Nitsche then one obtains reasonable numerical results. In many cases, this holds even in the limit case where the parameter of the CIP stabilization is zero, i.e., where the standard Galerkin discretization is applied. We present an analysis which explains this effect.

Full Text (PDF) [322 KB], BibTeX

Key words

diffusion-convection-reaction equation, finite elements, Nitsche type boundary conditions, error estimates

AMS subject classifications

65N15, 65N30, 65N50