El documento trata sobre el flujo de Stokes sujeto a las condiciones de contorno de fuga de umbral en dos y tres dimensiones espaciales. La formulación velocidad-presión conduce al problema de tipo desigualdad que es aproximado por los elementos finitos mixtos P1-burbuja/P1. El sistema algebraico resultante es no suave. Se resuelve mediante la variante de seguimiento de trayectoria del método del punto interior, y mediante la implementación de conjunto activo del método de Newton semi-suave. Los sistemas lineales internos se resuelven mediante el método de gradiente conjugado precondicionado. Experimentos numéricos ilustran la escalabilidad de los algoritmos. La novedad de este trabajo consiste en aplicar estrategias duales para resolver el problema.