Plitt Modified Hydrocyclone Alpha

Summary

The Plitt Modified Hydrocyclone Alpha is intended to provide a modified Plitt-family option to represents hydrocyclone classification using a modified empirical Plitt-family formulation. The model calculates pressure, volumetric split, corrected cut size, sharpness of separation, water bypass and the size-by-size partition to underflow.

DPSIM model key: DPSIM.Classification.PlittModifiedHydrocycloneAlpha
Category: Classification
Subcategory: Hydrocyclones
Display name: Plitt Modified Hydrocyclone Alpha

Parameters

# Parameter Description
1 # in parallel Number of hydrocyclones operating in parallel. The feed volumetric flowrate is divided by this value before applying the hydraulic correlations.
2 Hydrocyclone diameter (inch) Hydrocyclone body diameter used in the pressure, split, sharpness and cut-size correlations.
3 Height (inch) Free vortex height or effective distance term used in the hydrocyclone correlations.
4 Apex diameter (inch) Underflow or apex diameter.
5 Vortex diameter (inch) Overflow or vortex finder diameter.
6 Inlet diameter (inch) Feed inlet diameter.
7 Coarse bypass (fraction) Fractional coarse bypass correction. It reduces the classificable part of the corrected partition curve. The implementation limits this value so that R_f+B_pc does not exceed 1.
8 Pressure adjustment factor Multiplying factor applied to the calculated pressure/head correlation.
9 By-pass adjustment factor Multiplying factor used in the liquid/feed bypass calculation.
10 Enable Pressure Control (0-disabled/1-enabled) DPSIM control extension. When enabled, the model can change the number of operating hydrocyclones based on calculated pressure.
11 High Pressure Threshold (kPa) Pressure above which the model opens one additional hydrocyclone, subject to the maximum number available.
12 Low Pressure Threshold (kPa) Pressure below which the model closes one hydrocyclone, subject to the minimum number available.
13 Minimum HC available Minimum number of hydrocyclones allowed by the pressure-control extension.
14 Maximum HC available Maximum number of hydrocyclones allowed by the pressure-control extension.
15 Open HC change delay (min) Minimum delay between consecutive hydrocyclone opening or closing actions.
16 [Component] d50 factor Component-specific multiplying factor applied to the calculated corrected cut size.
17 [Component] imperfection factor Component-specific multiplying factor applied to the sharpness parameter m. Although the interface label uses “imperfection factor”, the implementation multiplies m directly.
18 [Component] split factor Component-specific multiplying factor applied to the volumetric split S before calculating the component bypass and component partition.

Derived parameters

# Derived parameter Description Unit
1 Pressure Calculated pressure after applying the pressure adjustment factor. kPa

Model Description

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

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

Q_i=Q_F/N_parallel

Where:

Symbol Description Unit
Q_i Volumetric feed flowrate per operating hydrocyclone. m3/h
Q_F Total feed volumetric flowrate. m3/h
N_parallel Number of hydrocyclones in parallel. dimensionless

The model uses the feed volumetric solids fraction:

φ_v=V_S^F/V_P^F

Where:

Symbol Description Unit
φ_v Feed volumetric solids fraction. fraction
V_S^F Volumetric flow of solids in the feed. m3/h
V_P^F Volumetric flow of pulp in the feed. m3/h

The implemented head-like pressure correlation is:

H=A_1 Q_i^1.46 exp(-7.63φ_v+10.79φ_v^2)/(D_c^0.20 h^0.15 D_i^0.51 D_o^1.65 D_u^0.53)

The pressure is then calculated as:

P=H ρ_p 0.433 6.895 F_H

Where:

Symbol Description Unit
H Internal head-like pressure term calculated by the model. model unit
P Calculated pressure. kPa
A_1 Alpha pressure constant, equal to 9.67960521515 in the implementation. model constant
D_c Hydrocyclone diameter. inch
h Height parameter. inch
D_i Inlet diameter. inch
D_o Vortex finder diameter. inch
D_u Apex diameter. inch
ρ_p Feed pulp specific gravity. dimensionless
F_H Pressure adjustment factor. dimensionless

Difference from Plitt Original: Alpha does not use the original metric pressure equation directly. It uses an alternative fitted head/pressure correlation with geometry entered in inches and flowrate in m3/h per cyclone.

The base volumetric split ratio is calculated as:

S=A_3 h^0.19 (D_u/D_o)^2.64 exp(-4.33φ_v+8.77φ_v^2)/(H^0.54 D_c^0.38)

Where:

Symbol Description Unit
S Volumetric underflow-to-overflow split ratio. dimensionless
A_3 Alpha split constant, equal to 54.9640283243 in the implementation. model constant

Difference from Plitt Original: the original Plitt split equation uses a different coefficient, includes (D_u^2+D_o^2)^0.36, and uses an exponential term based on percent solids by volume. Alpha uses a different fitted split equation with a quadratic solids-concentration term.

For each component c, the component split is:

S_c=S F_(S,c)

Where:

Symbol Description Unit
S_c Component-adjusted split ratio. dimensionless
F_(S,c) Component split factor. dimensionless

The sharpness parameter is calculated as:

m=exp(A_4-1.58 S_c/(S_c+1))(D_c^2 h/Q_i)^0.15

The component-adjusted sharpness is:

m_c=m F_(m,c)

Where:

Symbol Description Unit
m Base sharpness of separation. dimensionless
m_c Component-adjusted sharpness of separation. dimensionless
A_4 Alpha sharpness constant, equal to 0.522977577202 in the implementation. model constant
F_(m,c) Component sharpness factor. dimensionless

Difference from Plitt Original: Plitt relates sharpness to the volumetric recovery to underflow and the geometric term D_c^2 h/Q. Alpha keeps the same general dependence but uses the modified Alpha coefficient A4 and the component-adjusted split.

The corrected cut size is calculated as:

d_(50,c)=A_2 D_c^0.44 D_i^0.58 D_o^1.91 exp(11.12φ_v)/(D_u^0.80 h^0.37 Q_i^0.44 (ρ_(s,c)-1)^0.5)

The calibrated component cut size is:

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

Where:

Symbol Description Unit
d_(50,c) Corrected cut size for component c before user factor. µm
d_(50,c)^* Calibrated corrected cut size for component c. µm
A_2 Alpha cut-size constant, equal to 1.40133094302 in the implementation. model constant
ρ_(s,c) Component solids specific gravity. dimensionless
F_(d50,c) Component d50 factor. dimensionless

Difference from Plitt Original: the Alpha cut-size equation does not use the original Plitt exponents and constant. The dependence on D_o is stronger, the density term is simplified as (ρ_s-1)^0.5, and the solids-concentration correction is exp(11.12φ_v) rather than the original metric expression.

The corrected Rosin-Rammler partition curve is:

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 size of size class i, calculated from the DPSIM size mesh. µm

The model calculates a hypothetical solids recovery from the feed PSD and the corrected curve:

R_sc=sum_i f_i Y'_i

Where:

Symbol Description Unit
R_sc Hypothetical solids recovery to underflow used in the bypass equation. fraction
f_i Feed retained mass fraction in size class i. fraction

The water bypass to underflow is calculated as:

B_pw=(S_c/(S_c+1)-φ_v R_sc)/(1-φ_v(1-λ(1-R_sc)))

The feed bypass term is:

R_f=λ B_pw

Both B_pw and R_f are limited to the range from 0 to 1.

Where:

Symbol Description Unit
B_pw Water bypass fraction to underflow. fraction
R_f Feed bypass or fines bypass term used in the solids partition curve. fraction
λ By-pass adjustment factor. dimensionless

Difference from Plitt Original: Plitt Original calculates liquid recovery to underflow through a relationship involving volumetric recovery and solids recovery. Alpha uses a modified lambda-adjusted bypass expression and does not solve a full iterative Plitt liquid-recovery calculation.

The coarse bypass is limited internally:

B_pc^*=min(B_pc,1-R_f)

The final partition to underflow is:

E_(c,i)=R_f+(1-R_f-B_pc^*)Y'_(c,i)

Where:

Symbol Description Unit
E_(c,i) Final partition to underflow for component c and size class i. fraction
B_pc^* Limited coarse bypass term. fraction

Difference from Plitt Original: the explicit coarse bypass term B_pc is a DPSIM/Alpha extension. It is not part of the original Plitt corrected partition formulation.

The underflow and overflow component-size masses are calculated as:

M_(U,c,i)=E_(c,i) M_(F,c,i)

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

The total retained size distributions are calculated by summing over components:

M_(U,i)=sum_c M_(U,c,i)

M_(O,i)=sum_c M_(O,c,i)

p_(U,i)=M_(U,i)/sum_j M_(U,j)

p_(O,i)=M_(O,i)/sum_j M_(O,j)

The component fractions in each size interval are:

z_(U,c,i)=M_(U,c,i)/M_(U,i)

z_(O,c,i)=M_(O,c,i)/M_(O,i)

Where:

Symbol Description Unit
M_(F,c,i) Feed mass flowrate of component c in size class i. tph
M_(U,c,i) Underflow mass flowrate of component c in size class i. tph
M_(O,c,i) Overflow mass flowrate of component c in size class i. tph
p_(U,i) Underflow retained size fraction in size class i. fraction
p_(O,i) Overflow retained size fraction in size class i. fraction
z_(U,c,i) Component fraction of component c in underflow size class i. fraction
z_(O,c,i) Component fraction of component c in overflow size class i. fraction

The water split is calculated from the base hydraulic water bypass:

W_U=B_pw,base W_F

W_O=W_F-W_U

Where:

Symbol Description Unit
W_U Water flowrate to underflow. tph
W_O Water flowrate to overflow. tph
W_F Feed water flowrate. tph
B_pw,base Base hydraulic water bypass before component-specific split factors. fraction

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

The model preserves total solids by splitting each component-size mass between underflow and overflow. Water is split using the base hydraulic water bypass.

The parameters d50 factor, imperfection factor and split factor are component-specific calibration factors.

The pressure-control extension is not part of the original Plitt model.

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 operating hydrocyclones is increased by one, limited by the maximum number available. If pressure is below the low threshold, the number of operating hydrocyclones is decreased by one, limited by the minimum number available.

The hydraulic constants A1, A2, A3 and A4 are empirical and must be calibrated before design use.

References

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