In this paper, a numerical model for the flow of foams in flotation cells in up to three dimensions is described. By using finite elements to solve Laplace's equation for a potential scalar, the flow pattern and velocity of the foam can be obtained. This has allowed, for the first time, the carrying out of numerical investigations of foam flow patterns regardless of the symmetry of the system. The numerical model has been implemented in Fluidity, a general purpose finite element method code that allows simulations to be carried out on anisotropic unstructured meshes. Such a model is of particular interest since other phenomena occurring in flotation depend on the foam velocity. Convergence tests that verify the numerical model are presented, and simulation results that showcase important features of the present implementation are discussed. (C) 2011 Elsevier Ltd. All rights reserved.