Size-reduction systems have been extensively used in industry for many years. Nevertheless, reliable engineering tools to be used to predict the comminution of particles are scarce. Computational fluid dynamics(CFD)–discrete element model(DEM) numerical simulation may be used to predict such a complex phenomenon and therefore establish a proper design and optimization model for comminution systems.They may also be used to predict attrition in systems where particle attrition is significant. Therefore,empirical comminution functions(which are applicable for any attrition/comminution process), such as:strength distribution, selection, equivalence, breakage, and fatigue, have been integrated into the threedimensional CFD–DEM simulation tool. The main drawback of such a design tool is the long computational time required owing to the large number of particles and the minute time-step required to maintain a steady solution while simulating the flow of particulate materials with very fine particles.The present study developed several methods to accelerate CFD–DEM simulations: reducing the number of operations carried out at the single-particle level, constructing a DEM grid detached from the CFD grid enabling a no binary search, generating a sub-grid within the DEM grid to enable a no binary search for fine particles, and increasing the computational time-step and eliminating the finest particles in the simulation while still tracking their contribution to the process.The total speedup of the simulation process without the elimination of the finest particles was a factor of about 17. The elimination of the finest particles gave additional speedup of a factor of at least 18.Therefore, the simulation of a grinding process can run at least 300 times faster than the conventional method in which a standard no binary search is employed and the smallest particles are tracked. Size-reduction systems have been extensively used in industry for many years. Nevertheless, reliable engineering tools to be used to predict the comminution of particles are scarce. Computational fluid dynamics (CFD) -discrete element model (DEM) numerical simulation may be used to predict such a complex phenomenon and therefore establish a proper design and optimization model for comminution systems. They may also be used to predict attrition in systems where particle attrition is significant. Therefore, empirical comminution functions (which are applicable for any attrition / comminution process) , such as: strength distribution, selection, equivalence, breakage, and fatigue, have been integrated into the threedimensional CFD-DEM simulation tool. The main drawback of such a design tool is the long computational time required due to the large number of particles and the minute time-step required to maintain a steady solution while simulating the flow of particulate materials with very fin e particles. The present study developed several methods to accelerate CFD-DEM simulations: reducing the number of operations carried out at the single-particle level, constructing a DEM grid detached from the CFD grid enabling a no binary search, generating a sub-grid within the DEM grid to enable a no binary search for fine particles, and increasing the computational time-step and eliminating the fine particles in the simulation while still tracking their contribution to the process. the total speedup of the simulation process without the elimination of the finest particles was a factor of about 17. The elimination of the finest particles gave additional speedup of a factor of at least 18.Therefore, the simulation of a grinding process can run at least 300 times faster than the conventional method in which a standard no binary search is employed and the smallest particles are tracked.