This presentation discusses an anisotropic adaptive strategy in the context of the 3rd AIAA high-lift prediction workshop. If anisotropic mesh adaptation has proven its reliability for inviscid flows [1, 2], additional challenges remain to be solved to have the full gain of adaptivity, including early asymptotic (spatial second) order convergence, early capturing of the scales of then physical phenomena,. .. Several (fundamental) modifications are needed in the classical adaptive loop to address complex viscous effects. This includes the way error estimates are evaluated, how viscous solutions are interpolated between (anisotropic) meshes, and finally the generation of the boundary layer mesh to comply with the metric size prescription. Addressing fully each of this component remains a challenge on itself. In this presentation, we propose two error estimates to address the RANS equations, this implies in particular to treat appropriately the considered turbulence model. Here, we only consider the one equation Spalart-Allmaras turbulence model.