Data from ball mills of a wide range of sizes are analysed using detailed mathematical models and non-linear least squares fitting techniques. In this way model parameters are obtained describing the behaviour of ball mills over a wide range of mill sizes and operating conditions. These parameters define the flow and breakage of each size fraction of ore in the mill. Flow of coarse particles through the mill is found to depend on the proportion of coarse particles present. The breakage rates of coarse particles are found to decrease markedly with mill diameter for mills over 3 m in diameter.Simulation of grinding circuit throughput at constant productsizing for different mill diameters shows that throughput is proportional to theoretical mill power consumption for mills up to about 3 m in diameter. Thereafter the throughput gradually falls below that indicated by power consumption, the decrease becoming quite significant for large diameter mills. Simple qualitative physical arguments are presented to explain this behaviour. The models developed have applications to ball mill scale up. Work on their further development is continuing.