Numerical solution of the Giesekus model for incompressible free surface flows without solvent viscosity