Newton-Multigrid FEM Solver for the Simulation of Quasi-Newtonian Modeling of Thixotropic Flows