A perturbed twofold saddle point-based mixed finite element method for the Navier-Stokes equations with variable viscosity Artículo académico