Hi,
I have found out that one of the best treatments to stabilizing a reaction-dominated diffusion-advection-reaction (DAR) PDE is using bubble stabilization, see here for example. I also know that Fenics offers bubble enrichment for a DAR system in 2D in the following form:
P1 = FiniteElement('P', triangle,1)
PB = FiniteElement('B', triangle,3)
NEE = NodalEnrichedElement(P1, PB)
element = MixedElement([NEE,NEE,NEE])
Mixed_Space = FunctionSpace(mesh, element)
I was wondering if this is what the paper is talking about. Is it really as easy as this?