Here, we describe a numerical method which combines a mixed nite volume//nite element formulation on unstructured triangulations and a linearized implicit method for time advancing. A particular attention has been paid to the source term implicitation. This method is applied to three typical low-Reynolds number k ? " models, namely those of Lam-Bremhorst, Nagano-Tagawa and Speziale-Abid-Anderso...