Hydrocyclone Cut-Size Estimation Using Artificial Neural Networks
Hydrocyclone Cut-Size Estimation Using Artificial Neural Networks
including Biosystems
June 6-8, 2016. NTNU, Trondheim, Norway
Abstract: The hydrocyclone is widely used throughout the mineral processing industry when working with
slurries. It is either used for classifying, desliming or dewatering. Hydrocyclones are inexpensive,
application-efficient and relatively small to employ. In order to quantify its separation efficiency, models
are utilised to estimate the cut-size and sharpness of classification coefficient, usually in the form of a
partition curve. Most models are based on experimentally obtained data and are therefore not always
universally applicable. Over the last decade researchers have started employing Artificial Neural Networks
(ANNs) in order to obtain a dynamic model. This study endeavoured to use experimentally acquired data
to develop models that predict the cut-size. The models are discussed and evaluated in detail and the best
predicting model was compared to a conventional model from literature.
Keywords: Artificial Neural Network, Modelling, Hydrocyclone, Cut-size, Partition curve, Plitt-Flintoff.
quantitatively depict the separation efficiency of the
1. INTRODUCTION
hydrocyclone. Ideally a hydrocyclone is operated at
Ever since the hydrocyclone became popular in the mineral conditions where a specific cut-size and sharpness of
processing industry, there have been researchers that worked classification is achieved. These parameters cannot however
on developing a model describing its separation efficiency. be monitored in real-time (Frachon & Cilliers, 1999). This
Bradley (1965) published a book in which the known paper thus aimed at developing ANNs, based on
hydrocyclone fundamentals and research of that time were experimentally obtained data, which could predict the cut-size
detailed. He also compiled an extensive list of theoretical and sharpness of classification parameters at certain operating
equations estimating the cut-size and sharpness of conditions. An interesting and relatively new operating
classification, among others. These equations however were variable, called the angle of discharge, is included in some of
not always relevant to industrial hydrocyclones. The next the models.
important contributions made regarding the modelling of
hydrocyclones were by Lynch & Rao (1975). Their research A brief overview of a hydrocyclone is given in Section 2 in
was mainly focused around developing empirical equations for terms of what it is, where it is used, variables associated with
industrial hydrocyclones. In 1976, Plitt published a paper on it and how the performance is described. Section 3 focusses
his mathematical model of a hydrocyclone. His model is one on the ANN models that were developed, discussing their
of the most popular and most extensively referenced models. inputs, architecture and sample division. In order to check the
Plitt developed the model by utilising empirical data along adequacy of the developed ANN models, Analysis of Variance
with hydrocyclone variables that were deemed important in (ANOVA) studies were done and is shown in Section 4.1. To
describing the hydrocyclone’s operation. Flintoff et al. (1987) determine which one of the models performed the best, three
later reviewed the mathematical model and published a revised popular error metrics were utilised. Finally the best ANN
article on it where useful calibration factors were added to the model was compared to Plitt-Flintoff’s mathematical model
mathematical model to further improve estimations. estimations, in order to determine whether the ANN would
perform better than the mathematical model, given in Section
With the advance in computational capability, models now 4.2. The paper is concluded by outlining the work done and
include Computational Fluid Dynamics (CFD), System discussing the most important aspects of the study’s findings.
Identification and expanded statistical correlations. In 1997 H.
Eren et al. were some of the first researchers that employed
Artificial Neural Networks (ANNs) to predict the Particle Size 2. HYDROCYCLONE OVERVIEW
Distributions (PSD) and cut-size of various hydrocyclones (H
Eren, Fung, & Wong, 1997; Halit Eren, Fung, & Wong, 1997). 2.1 A general description
A hydrocyclone’s separation efficiency is mainly described by A hydrocyclone is a static, conical apparatus that is generally
the cut-size (𝑑50 ) and the sharpness of classification (𝑚). used within the mineral processing industry to separate solids
These two parameters are used within a partition curve to from water, better known as slurries. The separation of the
slurries is based on sedimentation, where the necessary swirl Table 1. Hydrocyclone variables
motion is generated by the slurry being fed into the
hydrocyclone by means of a pump (Bradley, 1965). Two Design variables
vortices form within the hydrocyclone referred to as the Hydrocyclone diameter 𝐷𝑐
primary vortex and the secondary vortex as depicted in Fig. 1 Feed inlet diameter 𝐷𝑖
(a). The primary vortex moves downwards and carries the Vortex finder diameter 𝐷𝑜
coarse particles to an opening called the underflow. The Spigot opening diameter 𝐷𝑢
secondary vortex carries the lighter particles, along with most Cone angle 𝜃
of the water, upwards to an opening called the overflow Free vortex height ℎ
(Frachon & Cilliers, 1999).
Operating variables
2.2 Hydrocyclone variables Inlet flow rate 𝑄𝑖
Overflow flow rate 𝑄𝑜
When working with hydrocyclones, two groups of variables Underflow flow rate 𝑄𝑢
are observed. The design variables (shown in Fig. 1 (b)) Pressure 𝑃
include variables that are dependent on the hydrocyclone’s Volumetric solid concentration 𝜙
size and design proportions. The operating variables are
Solid density 𝜌𝑆
independent of the hydrocyclone’s design and solely relate to
the operating conditions of the hydrocyclone. It should Overflow density 𝜌𝑜
however be noted that these two groups of variables cannot be Angle of discharge 𝜔
considered separately because of interactions that occur
between them. Table 1 summarises the main variables
allocated to the two groups.
Tangential inlet
Overflow
Overflow
Feed Cylindrical
Feed
section
Vortex finder
Secondary
vortex
Spigot
Underflow Underflow Fig. 2. A partition curve displaying the cut-size and sharpness
(a) (b) of classification coefficient.
Fig. 1. A graphical representation of a hydrocyclone
2.4 Experimental setup
depicting the (a) vortices and (b) relevant design variables
(adapted) (Frachon & Cilliers, 1999). A refurbished hydrocyclone test-rig, shown in Fig. 3, was
instrumented and used in order to acquire the necessary
2.3 Performance of a hydrocyclone
experimental data. The hydrocyclone’s dimensions are
A hydrocyclone’s performance is mainly described by its tabulated in Table 2. An analog pressure gauge was utilised to
separation efficiency which is directly quantifiable by a measure the inlet pressure and a Doppler flow meter to
partition curve. The partition curve describes the weight (or measure the inlet flow rate.
percentage) fraction of each particle size in the feed that might Table 2. The hydrocyclone dimensions
report to the underflow on the y-axis, to the specific particle
size, on the x-axis. An example of a partition curve is given in Design variable dimensions
Fig. 2 showing how the two efficiency indicating parameters, Hydrocyclone diameter 100 mm
cut-size and sharpness of classification coefficient, are related Feed inlet diameter 33.4 mm
to it. The cut-size, indicated as 𝑑50 , is defined as the size of Overflow diameter 34.0 mm
the particle in the Particle Size Distribution (PSD) that has a Free vortex height 531 mm
50% probability of reporting to either the underflow or the
overflow of the hydrocyclone (Bradley, 1965). The sharpness
of classification coefficient (𝑚) is a parameter that supplies a 3. ARTIFICIAL NEURAL NETWORKS
measure for the gradient of the partition curve. Ideally an 𝑚
> 3 is required in order to obtain sufficiently sharp separation 3.1 Artificial Neural Network overview
(Plitt, 1976).
997
IFAC DYCOPS-CAB, 2016
June 6-8, 2016. NTNU, Trondheim, Norway
998
IFAC DYCOPS-CAB, 2016
June 6-8, 2016. NTNU, Trondheim, Norway
to be significant and thereby considered adequate (Devore & of all three models are satisfactorily higher than the identified
Farnum, 2005). The ANOVA of the models are tabulated in critical F-value and it can thus be concluded that all three
Table 4. When evaluating the table it is seen that the F-values models are adequate.
The process was repeated for all three models and the results 30 30 30
30 32 34 36 38 40 30 32 34 36 38 40 30 32 34 36 38 40
are shown in Fig. 5. The dashed line represents the best fit Actual d50 (7m) Actual d50 (7m) Actual d50 (7m)
linear regression line between the actual and predicted cut- Fig. 5. The actual cut-size versus the predicted cut-size for (a)
sizes. The higher the coefficient-values the stronger the Model 1, (b) Model 2 and (c) Model 3.
relationship is expected to be. By comparing the coefficients
of the three models, one can observe that Model 3 performed The predicted cut-size per sample was then plotted on the same
the best. This was however the first assessment and some graph. The models’ predicted cut-size were expected to at
additional evaluation is needed, as these coefficients should least follow the trend of the actual cut-size. The graphs
never be the only measures examined (Devore & Farnum, obtained for the three models are given in Fig. 6. showing all
2005). the samples as used for training, validation and testing
(unknown samples). Very small differences are observed
The second measure used to evaluate the models was to between the three models but it is noted that the Model 2 shows
visually compare the predicted cut-sizes. The actual cut-size the worst correlation of the three. In order to better view the
was plotted per sample depicting the expected experimental testing samples’ actual and predicted cut-sizes, they were
error, calculated as 2.95%, as error bars.
999
IFAC DYCOPS-CAB, 2016
June 6-8, 2016. NTNU, Trondheim, Norway
separately plotted and are shown in Fig. 7. It was found once and examined. After computing the calibration factor for
more that the correlation of Model 2 was the worst. known samples, (5) was employed to estimate the cut-size for
the same unknown testing samples used to evaluate Model 3.
To further investigate the models’ performance, three popular
error metrics were calculated, assessed and compared. The 42
(a) Model 1 Actual vs Predicted d50
Actual
metrics used in this study was the Mean Square Error (MSE), 40
38
Model 1
d50 (7m)
the Root Mean Square Error (RMSE) and the Mean Absolute 36
34
Error (MAE). When comparing the three models in terms of 32
30
the three metrics, one would like to see the same ranking order 28
10 12 21 25 28 29 31 37
of the models for each one of the metrics. The metric error (b) Model 2 Actual vs Predicted d50
42
equations are given in (2), (3) and (4), respectively. 40 Actual
Model 2
38
d 50 (7m)
36
1 34
𝑀𝑆𝐸 = ∑𝑛𝑖=1(𝑦𝑖 − 𝑦̂𝑖 )2 , (2) 32
𝑛 30
28
1 6 7 12 19 25 26 31 39
𝑅𝑀𝑆𝐸 = √ ∑𝑛𝑖=1(𝑦𝑖 − 𝑦̂𝑖 )2 , (3) (c) Model 3 Actual vs Predicted d50
𝑛 42
Actual
1 40
𝑀𝐴𝐸 = ∑𝑛𝑖=1|𝑦𝑖 − ŷ𝑖 |, (4) 38
Model 3
d50 (7m)
𝑛 36
34
32
35
Model 3 0.5508 0.7422 0.5757
30
25
The conventional model that was investigated was the
1 2 4 5 6 7 8 9 111315 181922 232427 303233 343539 4041 3 141617 202636 381012 212528 293137
(b) Model 2 Actual vs Predicted d50
mathematical model Plitt and Flintoff revised in 1987.
45
Actual Training Validation Testing Equation (5) shows the mathematical model that was used to
40 calculate the cut-size (𝑑50 ) when the hydrocyclone design
d50 (7m)
1000
IFAC DYCOPS-CAB, 2016
June 6-8, 2016. NTNU, Trondheim, Norway
It is therefore concluded that not only can an ANN model be to see which variable(s) and/or ANN attributes used with
developed and used to predict the cut-size at specific operating Model 2 were responsible for the poor predictions.
conditions, but that an ANN model could substitute the
ACKNOWLEDGEMENTS
conventional mathematical model.
Actual vs Predicted d50 This work is based on the research supported in part by the
38
Actual National Research Foundation (NRF) of South Africa (Grant
Model 3
37 Plitt ! Flinto, Number 91093). Any opinion, finding and conclusion or
36 recommendation expressed in this material is that of the
author(s) and the NRF does not accept any liability in this
d50 (7m)
35
regard.
34
33
REFERENCES
32 Bradley, D. (1965). The Hydrocyclone - Bradley. London:
31
Pergamon Press Ltd.
1 2 11 16 23 24 31 36
Sample number
Cilliers, J. J., Austin, R. C., & Tucker, J. P. (1992). An
Fig. 8. Comparing the actual and predicted cut-size of evaluation of formal experimental design procedures. In
unknown samples employing Model 3 and Plitt-Flintoff’s International Conference on Hydrocyclones (pp. 31–49).
mathematical model.
Table 6. Summary of error functions when evaluating the Devore, J., & Farnum, N. (2005). Applied Statistics for
unknown samples for Model 3 and Plitt-Flintoff model Engineers and Scientists (2nd ed.). Brooks/Cole, Cengage.
Error functions Average Eren, H., Fung, C. C., & Gupta, A. (1996). Application of
Name of artificial neural network in estimation of hydrocyclone
MSE RMSE MAE errors parameters. In Diversity , the key to prosperity : the AusIMM
Model 3 0.5508 0.7422 0.5757 0.6229 1996 Annual Conference (pp. 225–229). Perth, Western
Plitt-Flintoff 0.8292 0.9106 0.7073 0.8157 Australia.
1001