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.