Minerals Beneficiation - Simulation of Nonlinear Grinding Systems: Rod-Mill Grinding

Simulation of nonlinear grinding systems is discussed in the context of the size-discretized batch-grinding model. A linear approximation of environment-dependent (or nonlinear) selection functions is postulated for systems which exhibit feed grinding orders less than unity. For this, a convex linear combination of zero and first-order selection function forms is used to obtain an analytical solution that can be used for simulation. An example in which marked nonlinear behavior occurs is rod-mill grinding. Accordingly, a series of grinding experiments was carried out in a batch laboratory rod mill to test these techniques; agreement between experiment and the model was found to be very good. Recently, significant progress has been made in the simulation of linear grinding systems by using a population balance mathematical model. Herbst and Fuer-stenaul presented a parameter estimation scheme that can be used with the size-discretized batch-grinding model, and Gumtz and Fuerstenau3 have extended these techniques to the simulation of locked-cycle grinding. Grandy et aL." summarized the simulation of linear systems and noted important alternate solutions to the batch-grinding equation. Kelsall et a1.'-" discussed the dynamic simulation of linear batch grinding, as well as linear open and closed-circuit grinding. In all of these, the key word is linear, which certainly restricts the engineer trying to simulate a variety of grinding systems. This paper discusses an approach for the simulation of non-linear grinding systems, particularly those which exhibit grinding orders less than unity. Experimentally, dry rod-mill grinding is used to evaluate the technique presented. Batch Grinding Model The size-discretized batch-grinding model has been discussed in detail by Mika et al.'.' and by Herbst and Fuerstenau,' and therefore, will only be briefly developed here. The model is formulated by making a differential mass balance of each size fraction in such a manner that the equation parameters are precisely defined in terms of selection and breakage functions,' which can be either calculated'.' or estimated1 from experimental data. The model is phenomenological because it averages all stress-application events in the system and, therefore, does not distinguish the manner in which individual particles are stressed.' For this average environment, the resistance of all particles in a size fraction to breakage is defined by that size-frac-tion's selection function; and when fracture occurs from an average event, the breakage product size distribution is defined again by that size-fraction's breakage function. The model is deterministic and is derived in terms of a continuous time variable and a discretized size variable,' i.e., the n mass fractions associated with the n defined size intervals (see Table 1 for nomenclature).