In this paper we propose stabilized finite element methods for both Stokes’ and Darcy’s problems that accommodate any interpolation of velocities and pressures. Apart from the interest of this fact, the important issue
is that we are able to deal with both problems at the same time, in a completely unified manner, in spite of the fact
that the functional setting is different. Concerning the stabilization formulation, we discuss the effect of the choice of the length scale appearing in the expression of the stabilization parameters, both in what refers to stability and to accuracy. This choice is shown to be crucial in the case of Darcy’s problem. As an additional feature of this work, we treat two types of stabilized formulations, showing that they have a very similar behavior.