The advection penalty does not work well unless the problem is advection only! Do not use for advection and diffusion problems. For this, we need a real slope limiter.