Rotary mills aim to effectively reduce the size of particles through a process called comminution. Modelling comminution in rotary mills is a challenging task due to substantial material deformation and the intricate interplay of particle kinematics of segregation, mixing, crushing, and abrasion. Existing particle-based simulations tend to provide predictions that cannot cope with the large number of particles within rotary mills, their wide range of sizes, and the physics dictating the crushing of individual particles. Similarly, there is currently no deterministic modelling means to determine the evolving population of particle sizes at any point in time and space within the mill. The aim of this two-part contribution is to address these gaps by advancing a framework for a novel stochastic comminution model for rotary mills, which has a well-defined deterministic continuum limit and can cope with arbitrarily large numbers of particles. This work describes the basic physics and structure of the new model within a heterarchical framework for ball and autogenous grinding mills. The primary focus of this Part I paper is to develop a computational model for the integration of motion of material along streamlines inside a mill. Coupled to this process is the kinetic physics dictating particle crushing. In a subsequent work, Part II, segregation and mixing will be added to this model such that realistic behaviour from the mill can be observed.Graphical