Four-field finite element solver and sensitivities for quasi-Newtonian flows