Physically realistic methods are a necessity to close the gap between reality and numerical result in modelling of tumbling mills. A problem is that tumbling mills often operate in a metastable state because of the difficulty to balance the rate of replenishment of large ore particles from the feed with the consumption in the charge. Understanding of the charge motion within the mill is of significance in mill optimisation. Both the breakage of ore particles and the wear of liners/ball media are closely linked to the charge motion. In this work, a ball charge and its interaction with the mill structure is modelled with the smoothed particle hydrodynamic (SPH) method. The mesh free formulation and the adaptive nature of the SPH method result...