Abstract. A novel generalized Newtonian model has been implemented in a lattice-Boltzmann (LB) algorithm. Viscosity is given as a function of the principal invariants of the deformation rate tensor. The collision term of the LB method is linked to the second moment of the velocity distribution via a non-linear algebraic equation. It is shown that this equation can be efficiently solved for an arbitrary rheological law.