Plitt Original Hydrocyclone

Summary

The Plitt Original Hydrocyclone model represents hydrocyclone classification using the empirical model proposed by L. R. Plitt in “A Mathematical Model of the Hydrocyclone Classifier”, CIM Bulletin, December 1976. The model estimates corrected cut size, flow split, pressure drop, sharpness of separation and the particle-size partition to underflow.

The model should be used when hydrocyclone performance must be estimated from cyclone geometry, feed flowrate, feed solids concentration, and calibration factors.

In this DPSIM implementation, the Plitt model is implemented using the metric equations from the original paper. The model includes calibration factors and an optional pressure-control extension that adjusts the number of hydrocyclones in operation based on calculated pressure.

DPSIM model key: DPSIM.Classification.PlittOriginal
Category: Classification
Subcategory: Hydrocyclones
Display name: Plitt Original Hydrocyclone

Parameters

# Parameter Description
1 # in parallel Number of hydrocyclones operating in parallel. The feed flowrate is divided by this value before applying the Plitt equations.
2 Hydrocyclone diameter (inch) Cyclone body diameter. The value is converted internally from inches to centimetres.
3 Height (inch) Free vortex height used in the Plitt equations. The value is converted internally from inches to centimetres.
4 Apex diameter (inch) Underflow or apex diameter. The value is converted internally from inches to centimetres.
5 Vortex diameter (inch) Overflow or vortex finder diameter. The value is converted internally from inches to centimetres.
6 Inlet diameter (inch) Feed inlet diameter. The value is converted internally from inches to centimetres.
7 Carrier liquid specific gravity Specific gravity of the liquid phase. This value is used in the density-difference term of the cut-size equation.
8 Global d50 calibration factor Global multiplying factor applied to the calculated Plitt cut size.
9 Global bypass calibration factor Global multiplying factor applied to the calculated liquid bypass to underflow.
10 Global partition efficiency factor Global multiplying factor applied to the calculated sharpness of separation.
11 Minimum class partition Lower limit imposed on the calculated class partition to underflow.
12 Maximum class partition Upper limit imposed on the calculated class partition to underflow.
13 Enable Pressure Control (0-disabled/1-enabled) DPSIM control extension. When enabled, the model can change the number of operating hydrocyclones according to calculated pressure.
14 High Pressure Threshold (kPa) Pressure above which the model opens one additional hydrocyclone, subject to the maximum number available.
15 Low Pressure Threshold (kPa) Pressure below which the model closes one hydrocyclone, subject to the minimum number available.
16 Minimum HC available Minimum number of hydrocyclones allowed by the pressure-control extension.
17 Maximum HC available Maximum number of hydrocyclones allowed by the pressure-control extension.
18 Open HC change delay (min) Minimum delay between consecutive hydrocyclone opening or closing actions.
19 [Component] d50 calibration Component-specific multiplying factor applied to the calculated Plitt cut size.
20 [Component] bypass calibration Component-specific multiplying factor applied to the liquid bypass term used in the partition curve.
21 [Component] partition efficiency calibration Component-specific multiplying factor applied to the sharpness of separation.

Derived parameters

Derived parameters

# Derived parameter Description Unit
1 Pressure Calculated cyclone pressure drop. kPa
2 Hydraulic head Hydraulic head calculated from pressure and pulp density. m
3 Flow split S Volumetric underflow-to-overflow split calculated by the Plitt flow split equation. dimensionless
4 Liquid bypass Rf to underflow Calculated liquid recovery to underflow, obtained by iteration from the volumetric split and solids recovery. fraction
5 Solids recovery Rs to underflow Calculated total solids recovery to underflow. fraction
6 Reference component d50 Calculated cut size for the first component, after global and component calibration. µm
7 Reference sharpness Calculated sharpness of separation before component-specific diagnostics. dimensionless

Model Description

The Plitt Original Hydrocyclone model receives one feed stream and generates two product streams: underflow and overflow. In DPSIM, the product port represents the underflow and the tail port represents the overflow.

The feed volumetric flowrate is divided by the number of hydrocyclones in parallel:

Q_cyc=Q_F/N_parallel

The individual cyclone flowrate is converted to litres per minute:

Q=1000 Q_cyc/60

The geometric parameters are converted from inches to centimetres:

D_cm=2.54 D_in

The feed solids volume fraction is converted to percent:

φ=100 φ_v

Where:

Symbol Description Unit
Q_cyc Volumetric feed flowrate per cyclone. m3/h
Q Volumetric feed flowrate per cyclone used in the Plitt equations. L/min
N_parallel Number of hydrocyclones in parallel. dimensionless
φ_v Feed solids volume fraction. fraction
φ Feed solids volume percent. %

The pressure drop is calculated as:

P=1.88 Q^1.78 exp(0.0055φ)/(D_c^0.37 D_i^0.94 h^0.28 (D_u^2+D_o^2)^0.87)

Where:

Symbol Description Unit
P Pressure drop. kPa
D_c Hydrocyclone diameter. cm
D_i Inlet diameter. cm
D_u Apex diameter. cm
D_o Vortex finder diameter. cm
h Free vortex height. cm

The hydraulic head is calculated from pressure and pulp density:

H=P/(9.81 ρ_p)

Where:

Symbol Description Unit
H Hydraulic head. m
ρ_p Feed pulp specific gravity. dimensionless

The flow split is calculated as:

S=1.9 (D_u/D_o)^3.31 h^0.54 (D_u^2+D_o^2)^0.36 exp(0.0054φ)/(H^0.24 D_c^1.11)

The volumetric recovery to underflow is:

R_v=S/(S+1)

Where:

Symbol Description Unit
S Volumetric underflow-to-overflow flow split. dimensionless
R_v Volumetric recovery to underflow. fraction

For each component c, the corrected cut size is calculated as:

d_(50,c)=52.5 D_c^0.46 D_i^0.60 D_o^1.21 exp(0.063φ)/(D_u^0.71 h^0.38 Q^0.45 (ρ_(s,c)-ρ_l)^0.5)

Then calibration factors are applied:

d_(50,c)^*=d_(50,c) F_d50 F_(d50,c)

Where:

Symbol Description Unit
d_(50,c) Plitt corrected cut size for component c. µm
d_(50,c)^* Calibrated cut size for component c. µm
ρ_(s,c) Component solids specific gravity. g/cm3
ρ_l Carrier liquid specific gravity. g/cm3
F_d50 Global d50 calibration factor. dimensionless
F_(d50,c) Component d50 calibration factor. dimensionless

The sharpness of separation is calculated as:

m=1.94 exp(-1.58 R_v) (D_c^2 h/Q)^0.15

The calibrated sharpness for component c is:

m_c=m F_m F_(m,c)

Where:

Symbol Description Unit
m Plitt sharpness of separation. dimensionless
m_c Calibrated sharpness for component c. dimensionless
F_m Global partition efficiency factor. dimensionless
F_(m,c) Component partition efficiency calibration factor. dimensionless

The corrected partition curve is calculated as:

Y'(c,i)=1-exp(-0.693 (d_i/d(50,c)^*)^m_c)

Where:

Symbol Description Unit
Y'_(c,i) Corrected partition to underflow for component c and size class i. fraction
d_i Representative particle size of size class i. µm

The liquid bypass to underflow is calculated iteratively. Starting from R_f=R_v, the model calculates solids recovery:

R_s=sum_c sum_i M_(c,i)^F (R_f+(1-R_f)Y'_(c,i))/M_S^F

The liquid recovery to underflow is updated as:

R_f=(R_v-R_s φ_v)/(1-φ_v)

The iteration is repeated until convergence or until the maximum iteration count is reached.

After convergence, the liquid bypass is calibrated:

R_f^*=R_f F_b

For each component:

R_(f,c)=R_f^* F_(b,c)

Where:

Symbol Description Unit
R_f Liquid recovery to underflow calculated by iteration. fraction
R_f^* Calibrated global liquid recovery to underflow. fraction
R_(f,c) Component-adjusted bypass term. fraction
R_s Solids recovery to underflow. fraction
F_b Global bypass calibration factor. dimensionless
F_(b,c) Component bypass calibration factor. dimensionless

The final class partition to underflow is:

Y_(c,i)=R_(f,c)+(1-R_(f,c))Y'_(c,i)

The partition is limited by the minimum and maximum class partition parameters:

Y_(c,i)^*=min(Y_max,max(Y_min,Y_(c,i)))

The underflow and overflow component masses are calculated as:

M_(c,i)^U=Y_(c,i)^* M_(c,i)^F

M_(c,i)^O=M_(c,i)^F-M_(c,i)^U

The stream size distributions and component-by-size matrices are rebuilt from these component masses.

The water split is calculated from the calibrated global liquid bypass:

M_W^U=R_f^* M_W^F

M_W^O=M_W^F-M_W^U

Pressure Control Extension

The pressure-control extension is not part of the original Plitt model. It is a DPSIM operational logic for adjusting the number of hydrocyclones in operation.

When pressure control is enabled, the model compares the calculated pressure with the high and low pressure thresholds. If pressure is above the high threshold, the number of hydrocyclones in parallel is increased by one, limited by the maximum available hydrocyclones. If pressure is below the low threshold, the number of hydrocyclones in parallel is decreased by one, limited by the minimum available hydrocyclones.

The model uses the metric Plitt equations internally. The user enters hydrocyclone geometry in inches, and the model converts these dimensions to centimetres before applying the equations.

The product stream corresponds to the underflow. The tail stream corresponds to the overflow.

The model preserves total solids and water mass by splitting the feed into underflow and overflow. Solids are split size-by-size and component-by-component according to the calculated partition curve. Water is split according to the calculated liquid bypass to underflow.

The d50, bypass and partition efficiency calibration parameters are DPSIM extensions. They allow calibration of the original Plitt equations against plant data, pilot data or vendor information.

The optional feed-size correction to d50 presented by Plitt is not implemented in this DPSIM version.

The model should be calibrated before being used for design conclusions, especially outside the operating range of the original Plitt data set.

References

Plitt, L. R. (1976). A Mathematical Model of the Hydrocyclone Classifier. CIM Bulletin, December 1976, 114–123.