Plitt Hydrocyclone Bravo

Summary

The Plitt Hydrocyclone Bravo model represents hydrocyclone classification using an alternative Plitt-type empirical formulation. It calculates pressure, hydraulic head, volumetric split, corrected cut size, sharpness of separation, water bypass and size-by-size partition to underflow.

This model is intended to provide an alternative empirical hydrocyclone calculation structure while retaining the general classification framework of Plitt-type hydrocyclone models.

DPSIM model key: DPSIM.Classification.PlittHydrocycloneBravo
Category: Classification
Subcategory: Hydrocyclones
Display name: Plitt Hydrocyclone Bravo

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 Cut size mode (0-exp/1-alt/2-alt2) Selects the cut-size calculation branch. Mode 0 uses an exponential solids-concentration correction. Modes 1 and 2 use the alternative cut-size branch.
8 Alternative cut-size parameter Parameter used only by the alternative cut-size branches.
9 Flow split correction Multiplying factor applied to the calculated flow split S.
10 Sharpness correction divisor Divisor applied to the calculated sharpness m. Higher values reduce the sharpness of the corrected partition curve.
11 Cut size correction Multiplying factor applied directly to the calculated cut size.
12 Minimum flow split S Lower bound imposed on the calculated flow split S.
13 Maximum flow split S Upper bound imposed on the calculated flow split S.
14 Maximum sharpness m Upper bound imposed on the calculated sharpness parameter m.
15 Include water bypass in solids partition (0/1) If set to 1, the solids partition includes the water bypass term. If set to 0, the solids partition uses only the corrected classification curve.
16 Minimum class partition Lower bound imposed on the final class partition to underflow.
17 Maximum class partition Upper bound imposed on the final class partition to underflow.
18 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.
19 High Pressure Threshold (kPa) Pressure above which the model opens one additional hydrocyclone, subject to the maximum number available.
20 Low Pressure Threshold (kPa) Pressure below which the model closes one hydrocyclone, subject to the minimum number available.
21 Minimum HC available Minimum number of hydrocyclones allowed by the pressure-control extension.
22 Maximum HC available Maximum number of hydrocyclones allowed by the pressure-control extension.
23 Open HC change delay (min) Minimum delay between consecutive hydrocyclone opening or closing actions.
24 [Component] d50 factor Component-specific multiplying factor applied to the calculated corrected cut size.
25 [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.
26 [Component] split factor Component-specific multiplying factor applied to the volumetric split S before calculating component bypass and component partition.

Derived parameters

# Derived parameter Description Unit
1 Pressure Calculated cyclone pressure drop. kPa

Model Description

The Plitt Hydrocyclone Bravo 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 feed volumetric solids fraction is:

φ_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 pressure correlation is calculated from the implemented pressure term:

P_int=129.72875 Q_i^1.78 exp(0.55φ_v)/(D_c^0.37 D_i^0.94 h^0.28 (D_u^2+D_o^2)^0.87)

The reported pressure is:

P=0.0980665 P_int

The hydraulic head used by the split equation is:

H=P_int/(9.8100004196 ρ_p)

Where:

Symbol Description Unit
P_int Internal pressure/head term used by the model. model unit
P Reported pressure drop. kPa
H Hydraulic head term used in the split equation. m-equivalent model basis
D_c Hydrocyclone diameter. inch
D_i Inlet diameter. inch
D_o Vortex finder diameter. inch
D_u Apex diameter. inch
h Height parameter. inch
ρ_p Feed pulp specific gravity. dimensionless

Difference from Plitt Original: the pressure equation retains the general Plitt-type structure but uses constants transformed for the implemented unit basis. The model uses inch geometry and m3/h per cyclone rather than directly using the original metric units.

The base volumetric split ratio is:

S=3.3411661493 (D_u/D_o)^3.31 h^0.54 (D_u^2+D_o^2)^0.36 exp(-0.8884871132φ_v) F_S/(H^0.24 D_c^1.11)

The calculated value is limited between the minimum and maximum flow split parameters:

S^*=min(S_max,max(S_min,S))

Where:

Symbol Description Unit
S Calculated volumetric underflow-to-overflow split ratio before limiting. dimensionless
S^* Limited volumetric split ratio. dimensionless
F_S Flow split correction. dimensionless
S_min Minimum flow split S. dimensionless
S_max Maximum flow split S. dimensionless

For each component c, the component-adjusted 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

Difference from Plitt Original: Bravo uses an alternative split coefficient and solids-concentration term, and it imposes explicit user-defined limits on S. These limits are not part of the original Plitt equation.

The water bypass to underflow is calculated as the volumetric recovery to underflow:

B_pw=S_c/(S_c+1)

The feed bypass term used in the partition curve is:

R_f=B_pw

Where:

Symbol Description Unit
B_pw Water bypass fraction to underflow. fraction
R_f Bypass term used in the solids partition curve. fraction

Difference from Plitt Original: the original Plitt formulation relates liquid recovery to volumetric recovery and solids recovery. Bravo does not perform an iterative Plitt liquid-recovery calculation; it uses R_f=B_pw=S/(S+1).

The sharpness of separation is calculated as:

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

The result is limited by the maximum sharpness parameter:

m^*=min(m,m_max)

The component-adjusted sharpness is:

m_c=m^* F_(m,c)

Where:

Symbol Description Unit
m Calculated sharpness before maximum limit. dimensionless
m^* Sharpness after maximum limit. dimensionless
m_c Component-adjusted sharpness. dimensionless
F_m Sharpness correction divisor. dimensionless
m_max Maximum sharpness m. dimensionless
F_(m,c) Component sharpness factor. dimensionless

Difference from Plitt Original: Bravo uses an alternative coefficient for the sharpness equation, applies a user correction divisor, and limits the maximum value of m.

The corrected cut size has two calculation branches.

For cut size mode 0, the exponential branch is:

d_50=F_cut 4418.82577186 D_c^0.46 D_i^0.60 D_o^1.21 exp(3.9300207955φ_v)/(D_u^0.71 h^0.38 Q_i^0.45 31.6227766017)

For cut size modes 1 and 2, the alternative branch is:

d_50=F_cut 4418.82577186 D_c^0.46 D_i^0.60 D_o^1.21 φ_v^0.41 A_alt^0.35/(D_u^0.71 h^0.38 Q_i^0.45 31.6227766017)

The model then applies a component density correction, a solids-load correction and the component d50 factor:

d_(50,c)^*=d_50 F_(ρ,c) F_load F_(d50,c)

with:

F_(ρ,c)=sqrt((ρ_p-1)/(ρ_(s,c)-1))

and:

F_load=(M_S^F/(100 N_parallel))^(-0.0465008346)

Where:

Symbol Description Unit
d_50 Base corrected cut size. µm
d_(50,c)^* Component-adjusted corrected cut size. µm
F_cut Cut size correction. dimensionless
A_alt Alternative cut-size parameter. dimensionless
F_(ρ,c) Component density correction factor. dimensionless
F_load Solids-load correction factor. dimensionless
F_(d50,c) Component d50 factor. dimensionless
M_S^F Feed dry solids flowrate. tph
ρ_(s,c) Component solids specific gravity. dimensionless

Difference from Plitt Original: Bravo includes alternative cut-size modes, a solids-load correction, a component density correction based on pulp density and component SG, and a global cut-size correction factor. These are not part of the original Plitt equation.

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

If “Include water bypass in solids partition” is enabled, the final partition is:

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

If the option is disabled, the final partition is:

E_(c,i)=Y'_(c,i)

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

E_(c,i)^*=min(E_max,max(E_min,E_(c,i)))

Where:

Symbol Description Unit
E_(c,i) Final partition to underflow before limiting. fraction
E_(c,i)^* Final partition to underflow after limiting. fraction
E_min Minimum class partition. fraction
E_max Maximum class partition. fraction

Difference from Plitt Original: Bravo provides an explicit switch to include or exclude bypass in the solids partition curve and applies minimum and maximum class partition limits. These are DPSIM compatibility and numerical robustness extensions.

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 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 pressure-control action is applied after the current model calculation. Therefore, the new number of hydrocyclones affects the following calculation cycle.

The alternative cut-size modes, load correction, density correction, split limits, maximum sharpness and class-partition limits are DPSIM calibration and robustness extensions.

The use of R_f=B_pw=S/(S+1) is a simplification relative to the original Plitt liquid-recovery treatment.

The model should be calibrated against plant surveys, pilot data or vendor data before being used for design conclusions.

References

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