The phenomenological theory of continuous thickening of flocculated suspensions in an ideal cylindrical thickener is extended to vessels having varying cross-section, including divergent or convergent conical vessels. The purpose of this contribution is to draw attention to the corresponding mathematical model, whose key ingredient is a strongly degenerate parabolic partial differential equation. For ideal (non-flocculated) suspensions, which do not form compressible sediments, the mathematical model reduces to the kinematic approach by Anestis, who developed a method of construction of exact solutions by the method of characteristics. The difficulty lies in the fact that characteristics and iso-concentration lines, unlike the conventional Kynch model for cylindrical vessels, do not coincide, and one has to resort to numerical methods to simulate the thickening process. A numerical algorithm is presented and employed for simulations of continuous thickening. Implications of the mathematical model are also demonstrated by steady-state calculations, which lead to new possibilities in thickener design. (C) 2004 Elsevier B.V. All rights reserved.