In the last few years there has been a growing interest in the development of numerical techniques appropriate for the approximation of differential model problems presenting multiscale solutions. This is the case, for instance, with functions displaying a smooth behavior, except in certain regions where sudden and sharp variations are localized. Typical examples are internal or boundary layers. When the number of degrees of freedom in the discretization process is not sufficient to ensure a fine resolution of the layers, some stabilization procedures are needed to avoid unpleasant...