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.