Modeling and Simulation of Hydroxyapatite Recovery in the Desliming Circuit of the Tapira Industrial Plant, Brazil

: The modeling and simulation of industrial mineral processing operations are traditionally used for cyclone sizing and optimizations of industrial operations. However, the main models used are based on the total population of particles in the pulp, thus not distinguishing the individual minerals. This article presents the results of an innovative method that investigated the optimization of the metallurgical recovery of P 2 O 5 in the desliming circuit of a phosphate ore processing plant in Brazil. A survey campaign was carried out in the existing industrial circuit, followed by determining the partition curves for the overall particles and specifically for the hydroxyapatite particles. The results were used to calibrate the Narasimha–Mainza cyclone model. From a Base Case determined with reference to the industrial survey, three optimization scenarios were simulated through cyclone geometries and respective operating conditions changes. Simulated scenarios indicated the possibility of P 2 O 5 metallurgical recovery increasing from 9.4% to 12.7% compared to the Base Case.


Introduction
Processing weathered apatitic phosphate ores includes separating the slimes' fraction to increase the selectivity in the downstream flotation concentration process [1][2][3][4].Although beneficial to the concentration, this condition also results in significant losses of phosphorusbearing minerals to the slimes.This combination results in the need to control each stage of an industrial desliming circuit to balance the benefits and losses of hydroxyapatite (HA) from the process.This last aspect is critical as slimes are diverted to tailings dams, which either reduces the life of these structures or increases the mass to be dewatered and piled up in high-cost operations, as recently occurring in Brazilian mining [5][6][7].
Cyclones are widely used in industrial classification operations mainly due to their versatility, high capacity and relative simplicity of operation.The forces acting on the particles determine downward trajectories for the relatively coarser ones and upward trajectories for the finer ones.The geometry of the cyclone, in combination with the operating conditions, modulates the partition of the solids in each particle size fraction.According to Delboni Jr. [8], the most widely used models for cyclone simulation were developed using hard ores, thus different from the relatively high amount of natural fines observed in most weathered Brazilian phosphate ores.Accordingly, the various mathematical models created in the last four decades [9][10][11][12][13] have been based on cyclones operating in comminution circuits, thus distinct from desliming circuits.Also, these models considered only one partition curve despite the individual minerals being classified.
Lynch and Rao [9] presented individual partition curves for the minerals limestone, silica and coal.Later in the same work, the authors showed the partition curve for the Warrego industrial operation and a theoretical partition curve based on the differences in specific weights of the minerals presented quartz and magnetite.In reference to the Warrego ore partition curve, the authors considered that the tail was because the d50c values for the ore were not the weighted average for the d50c values for the components but were biased toward the lower value.Further calculations showed that the shape of the curve was affected by the size distributions of the components and their relative proportions.The conclusion was that a composite reduced efficiency curve may change, but the curves for the components remain constant.
More recently, the Narasimha-Mainza [14] semi-mechanistic mode was developed based on a comprehensive database considering different operating and design conditions, showing greater accuracy than previously widely used models.Accordingly, Equations ( 1)-( 4) are used to predict the performance of the classification operation regarding, respectively, the recovery of water to the underflow (Rf), the corrected cut size (d50c), volumetric flow of slurry (Q) in the feed and sharpness of the separation (α).A fundamental point of the Narasimha-Mainza model is that although the density of solids is included in the equations, it refers to the bulk of particles, since the differences in performance for particles of different densities were not reported, or even validated.

Objectivea
The main objective of this work was to investigate the optimization of the metallurgical recovery of P 2 O 5 from the desliming stage of the phosphate ore beneficiation industrial circuit of Tapira-MG, Brazil, through a novel specific calibration and simulation of HA using the Narasimha-Mainza model.

Tapira Industrial Mineral Processing Plant
The Tapira phosphate ore processing plant, located in the state of Minas Gerais, Brazil, was implemented at the end of the 1970s and currently processes an average of 15.9 Mt per year, which results in a nominal feed flow rate of 1.950 t/h of solids.The annual production is 2 Mt of concentrate, thus ranking Tapira as the largest phosphate concentrator in Brazil.The plant feed grade of 8.5% P 2 O 5 is upgraded to 35% P 2 O 5 through initial multi-stage crushing steps, followed by a screening that splits the crushed feed into two different circuits.The hard ore circuit processes the coarse fraction, i.e., −19.1 + 7.0 mm, while the soft ore circuit is fed with the −7.0 mm fraction, representing nearly 80% of overall plant feed.
The soft ore circuit feed slurry is pumped and split into four parallel lines, each one including an initial classification stage carried out in a nest of 508 mm diameter cyclones, followed by processing the cyclone's coarse fraction in a dedicated grinding circuit, which also includes magnetite separation in magnetic drums, as shown in the flow sheet of Figure 1.
while the soft ore circuit is fed with the −7.0 mm fraction, representing nearly 80% of overall plant feed.
The soft ore circuit feed slurry is pumped and split into four parallel lines, each one including an initial classification stage carried out in a nest of 508 mm diameter cyclones, followed by processing the cyclone's coarse fraction in a dedicated grinding circuit, which also includes magnetite separation in magnetic drums, as shown in the flow sheet of Figure 1.The combined second desliming stage underflow is pumped to another classification stage, resulting in coarse and fine flotation feed circuits.In contrast, the combined second desliming stage overflow is transferred to the tailing dam.
The optimal operation of the second desliming stage consists of efficiently separating the slimes with the maximum possible recovery of HA to the combined underflow.While the former is essential for a high flotation performance, the latter modulates the P2O5 recovery in the industrial circuit.Based on the overall plant feed, the mass split to the two parallel second desliming stages is rated at 15%.Due to practical and operational aspects, the present study focused exclusively on one of the two secondary desliming lines.As the equipment and operating conditions are identical, the results were considered replicable for either of the two lines.

Method
A comprehensive survey was carried out at the second desliming stage at the Tapira industrial plant.It included sampling feed, underflow and overflow of one nest of cyclones, obtaining detailed data and information from the dedicated process information management system (PIMS), and measuring vortex and apex from selected equipment.The collected samples were processed in the Tapira Process Laboratory Unit (LPUT) to determine respective solids concentration, size distribution and chemical analysis of sizeby-size products.Further sample characterization was conducted at the Mineral Processing Laboratory of the Polytechnic School of the University of São Paulo (LTM-EPUSP) for assessing the size distribution of fines by Cyclosizer, as well as X-ray Diffraction (XRD) analysis at the Technological Ore Characterization Laboratory of the Polytechnic School  The combined second desliming stage underflow is pumped to another classification stage, resulting in coarse and fine flotation feed circuits.In contrast, the combined second desliming stage overflow is transferred to the tailing dam.
The optimal operation of the second desliming stage consists of efficiently separating the slimes with the maximum possible recovery of HA to the combined underflow.While the former is essential for a high flotation performance, the latter modulates the P 2 O 5 recovery in the industrial circuit.Based on the overall plant feed, the mass split to the two parallel second desliming stages is rated at 15%.Due to practical and operational aspects, the present study focused exclusively on one of the two secondary desliming lines.As the equipment and operating conditions are identical, the results were considered replicable for either of the two lines.

Method
A comprehensive survey was carried out at the second desliming stage at the Tapira industrial plant.It included sampling feed, underflow and overflow of one nest of cyclones, obtaining detailed data and information from the dedicated process information management system (PIMS), and measuring vortex and apex from selected equipment.The collected samples were processed in the Tapira Process Laboratory Unit (LPUT) to determine respective solids concentration, size distribution and chemical analysis of size-by-size products.Further sample characterization was conducted at the Mineral Processing Laboratory of the Polytechnic School of the University of São Paulo (LTM-EPUSP) for assessing the size distribution of fines by Cyclosizer, as well as X-ray Diffraction (XRD) analysis at the Technological Ore Characterization Laboratory of the Polytechnic School of the University of São Paulo (LCT-EPUSP).Mass balancing experimental data and mathematical model calibration were carried out using Excel 365 and JKSimMet 6.2 software.The flowsheet described in Figure 2 shows the sequence described above.

Survey Campaign
The industrial sampling campaign was carried out during a steady-state operation of the Tapira processing plant.Sampling points of the second desliming stage are listed in Table 1.During the stabilization and sampling periods, 13 cyclones were in full operation.
of the University of São Paulo (LCT-EPUSP).Mass balancing experimental data and mathematical model calibration were carried out using Excel 365 and JKSimMet 6.2 software.The flowsheet described in Figure 2 shows the sequence described above.

Survey Campaign
The industrial sampling campaign was carried out during a steady-state operation of the Tapira processing plant.Sampling points of the second desliming stage are listed in Table 1.During the stabilization and sampling periods, 13 cyclones were in full operation.1 consisted of combined flows, respectively, of the underflow and overflow of the 13 cyclones in operation.Sampling was carried out incrementally every 10 min for two hours to obtain sufficient mass for the sample processing for the three flows sampled.Circuit stabilization was assessed by three main indicators listed in the PIMS the plant fresh feed rate, cyclone nest feed pressure and pulp flow rate, and feed sump level.The cyclone nest feed flow rate of solids was calculated based on the pulp flow rate, as indicated by the PIMS, and the concentration of the solids obtained from the cyclone feed sample.Figure 3 shows a photograph of the secondary desliming stage installed at the Tapira industrial plant.Sampling points 2 and 3 in Table 1 consisted of combined flows, respectively, of the underflow and overflow of the 13 cyclones in operation.Sampling was carried out incrementally every 10 min for two hours to obtain sufficient mass for the sample processing for the three flows sampled.Circuit stabilization was assessed by three main indicators listed in the PIMS the plant fresh feed rate, cyclone nest feed pressure and pulp flow rate, and feed sump level.The cyclone nest feed flow rate of solids was calculated based on the pulp flow rate, as indicated by the PIMS, and the concentration of the solids obtained from the cyclone feed sample.Figure 3 shows a photograph of the secondary desliming stage installed at the Tapira industrial plant.

Sample Processing
The collected samples were initially weighed, filtered and then dried in an oven at the LPUT to determine the concentration of solids by weight (Cw).
Dry samples were split and screened at the following laboratory sieves: 0.500; 0.230; 0.152; 0.106; 0.074; 0.053; 0.044; and 0.037 mm.In each case, triplicate screening was carried out to obtain sufficient chemical analysis mass.Except for the fraction passing at 0.037 mm sieve, all other fractions were oven-dried, weighted and prepared for chemical analysis at the Mosaic laboratories.
Fractions passing at 0.037 mm were tested at a Cyclosizer installed at the LTM-EPUSP.Due to mass constrain as stipulated in Cyclosizer standard procedure, each sample from each stream, i.e., cyclone feed, underflow and overflow, was divided into six sub-samples for testing at the Cyclosizer, therefore totaling 18 Cyclosizer tests.Based on the results, only 0.025 mm (products of C1, C2 and C3 from the Cyclosizer) and 0.013 mm (products of C4 and C5 from the Cyclosizer) were adopted for reporting the Cyclosizer results.The obtained −0.037 + 0.025 mm, −0.025 + 0.013 mm and −0.013 mm fractions of cyclone feed, underflow and overflow were returned to the Tapira for chemical analysis.

Sample Processing
The collected samples were initially weighed, filtered and then dried in an oven at the LPUT to determine the concentration of solids by weight (Cw).
Dry samples were split and screened at the following laboratory sieves: 0.500; 0.230; 0.152; 0.106; 0.074; 0.053; 0.044; and 0.037 mm.In each case, triplicate screening was carried out to obtain sufficient chemical analysis mass.Except for the fraction passing at 0.037 mm sieve, all other fractions were oven-dried, weighted and prepared for chemical analysis at the Mosaic laboratories.
Fractions passing at 0.037 mm were tested at a Cyclosizer installed at the LTM-EPUSP.Due to mass constrain as stipulated in Cyclosizer standard procedure, each sample from each stream, i.e., cyclone feed, underflow and overflow, was divided into six subsamples for testing at the Cyclosizer, therefore totaling 18 Cyclosizer tests.Based on the results, only 0.025 mm (products of C1, C2 and C3 from the Cyclosizer) and 0.013 mm (products of C4 and C5 from the Cyclosizer) were adopted for reporting the Cyclosizer results.The obtained −0.037 + 0.025 mm, −0.025 + 0.013 mm and −0.013 mm fractions of cyclone feed, underflow and overflow were returned to the Tapira for chemical analysis.
The cyclone feed sample was further prepared for semi-quantitative mineralogical analysis via X-ray Diffraction (XRD) at LCT-EPUSP.

Mass Balance and Calibration of the Narasimnha-Mainza Model
The mass balance of experimental data was carried out based on the concentration of solids and size distributions of the cyclone feed, underflow and overflow, as well as on the cyclone feed solids' flow rate.Data from mass balance were further used to calibrate the Narasimha-Mainza model, thus resulting in the calibrated model, here referred to as Base Case.

Determination of the Partition Curve and Cyclone Model Calibration for HA
In addition to calibrating the Narasimnha-Mainza model for overall cyclone feed, underflow and overflow, the same model was calibrated specifically for the mineral HA, following the adopted procedures described in Sections 2.5.1-2.5.3.The cyclone feed sample was further prepared for semi-quantitative mineralogical analysis via X-ray Diffraction (XRD) at LCT-EPUSP.

Mass Balance and Calibration of the Narasimnha-Mainza Model
The mass balance of experimental data was carried out based on the concentration of solids and size distributions of the cyclone feed, underflow and overflow, as well as on the cyclone feed solids' flow rate.Data from mass balance were further used to calibrate the Narasimha-Mainza model, thus resulting in the calibrated model, here referred to as Base Case.

Determination of the Partition Curve and Cyclone Model Calibration for HA
In addition to calibrating the Narasimnha-Mainza model for overall cyclone feed, underflow and overflow, the same model was calibrated specifically for the mineral HA, following the adopted procedures described in Sections 2.5.1-2.5.3.

Determination of HA Distributions
The calculation of HA content in each size fraction of each stream around the cyclone was based on the XRD analysis results shown in Appendix A, which indicated that the only phosphorus-bearing mineral in the ore is HA.
The conversion of P 2 O 5 grades obtained from chemical analysis to HA grades was stoichiometric, i.e., by calculating the grade of P 2 O 5 in the HA formula.The result indicated a factor of 2.36.Equation ( 5) was used to calculate the experimental distributions of apatite in each particle size fraction for the feed, overflow and underflow. where: Minerals 2024, 14, 272 6 of 16 CAP i = HA content in size fraction i; x i = percent retained in size fraction i; a i = P 2 O 5 grade in size fraction i; i = number of size fractions; cf = conversion factor of P 2 O 5 to HA (2.36).
The results obtained from Equation ( 5) were balanced for consistency in the HA distributions in the cyclone feed, underflow and overflow.The consistency criterion adopted using a simple two-product formula to calculate the overall mass partition (Rm) around the cyclone streams.

Parameterization of HA Partition Curve
The partition of apatite by particle size fraction (Ui) was calculated.The obtained corrected partition distribution curve was parameterized by Whiten's formula, shown in Equation (6).
The parameter yi represents the ratio between particle size and d50c, the latter corresponding to the size associated with a 50% corrected partition.In the same Equation ( 6), the α parameter represents the dispersion of the corrected partition distribution around d50c.

Model Calibration
The parameters d50c and α, obtained from the parameterization of the HA partition curve, were used to calibrate the respective model equations in Narasimha-Mainza, i.e., Equations ( 2) and ( 4).The calibration of Equation ( 2), relative to parameter d50c, consisted of calculating the constant KD, while the calibration of Equation ( 4), relative to parameter α consisted of calculating the constant Kα.

Simulations
The simulations were carried out in two steps.In the first, the calibrated Narasimha-Mainza model was used to simulate overall mass and water flow rates around cyclone feed, underflow and overflow.The second simulation obtained the HA partitions based on specifically calibrated Equations ( 2) and ( 4), also based on the Narasimha-Mainza model.Accordingly, the obtained size distributions and HA distributions of underflow and overflow were used to calculate the P 2 O 5 grade in each size fraction of these two streams.The feed size distribution and HA distributions were kept constant for the cyclone feed stream.
Three complete simulations were conducted to increase the overall HA recovery to the deslimed flow (cyclone underflow), together with relatively small variations in the amount of fines in the same stream.

Results and Discussion
This section describes the results obtained from the industrial survey, sample characterization, mass balance and model calibration results.The simulation obtained for the second desliming stage results of the Tapira industrial plant is also presented and discussed.

Industrial Sampling Results
The main characteristics of the cyclones and the selected operating variables obtained during the survey campaign, are listed in Table 2.

Mass Balance and Model Calibration
This section presents the mass balance results and the calibration of the Narasimha-Mainza model.

Experimental and Mass Balance Data
Figure 4 shows experimental and mass balanced particle size distributions.Table A1 in Appendix B illustrates all the detailed data from the sampling.It can be seen that the experimental and balanced results are very close, illustrating the quality of the sampling.

Model Calibration for Overall Flows
The selected parameters for model calibration, as obtained from mass balance calculations, are listed in Table 3.The calibration constants obtained for the Narasimha-Mainza from Equations ( 1)-( 4) based on mass-balanced overall flows are listed in Table 4.

Model Calibration for Overall Flows
The selected parameters for model calibration, as obtained from mass balance calculations, are listed in Table 3.The calibration constants obtained for the Narasimha-Mainza from Equations ( 1)-( 4) based on mass-balanced overall flows are listed in Table 4.The corrected partition curve obtained from model calibration is shown in Figure 5.

HA Balancing and Partition Modelling
The results obtained from the HA partition curve parameterization and model calibration are presented in the following sections.

HA Distribution
Table 5 presents the experimental and balanced HA distributions for each of the three cyclone flows.Table A2 in Appendix B shows for each stream around the cyclones, the balanced particle size, together with experimental P2O5 grades and the resulted experimental HA distribution.Figure 6 shows the cyclone nest feed, overflow, and underflow distributions' balanced HA distributions.

HA Balancing and Partition Modelling
The results obtained from the HA partition curve parameterization and model calibration are presented in the following sections.

HA Distribution
Table 5 presents the experimental and balanced HA distributions for each of the three cyclone flows.Table A2 in Appendix B shows for each stream around the cyclones, the balanced particle size, together with experimental P 2 O 5 grades and the resulted experimental HA distribution.

Model Calibration for HA Partition
The balanced HA distributions were used to calculate the respective HA partitions, parameterized in terms of d50c and α, as shown in Table 6.Based on these two parameters, the constants KD and Kα, for Equations ( 2) and ( 4) of the Narasimha-Mainza model were calculated.The obtained values are shown in Table 7.

Selected HA Partition Curve Parameters
Corrected d50c (mm) 0.030 Sharpness of Efficiency Curve-α 1.75 HA Model Parameters d50c Constant-KD 0.032 Sharpness Constant-Kα 3.235 The obtained HA partition curve is shown in Figure 7.

Model Calibration for HA Partition
The balanced HA distributions were used to calculate the respective HA partitions, parameterized in terms of d50c and α, as shown in Table 6.Based on these two parameters, the constants KD and Kα, for Equations ( 2) and ( 4) of the Narasimha-Mainza model were calculated.The obtained values are shown in Table 7.The obtained HA partition curve is shown in Figure 7.

Model Calibration for HA Partition
The balanced HA distributions were used to calculate the respective HA partitions, parameterized in terms of d50c and α, as shown in Table 6.Based on these two parameters, the constants KD and Kα, for Equations ( 2) and ( 4) of the Narasimha-Mainza model were calculated.The obtained values are shown in Table 7.

Selected HA Partition Curve Parameters
Corrected d50c (mm) 0.030 Sharpness of Efficiency Curve-α 1.75 The obtained HA partition curve is shown in Figure 7.The first calibration exercise for the global sample indicated relatively small sharpness of the overall partition curve (α = 0.302), which may be misinterpreted as a low performance of the classification process.The second step, though, revealed a relatively high performance of the classification process specifically for the HA (α = 1.75).Thus, such a result is paramount, as an apparently low-performance classification is a high-performance one.This aspect was already masterfully assessed in the seminal work by Lynch and Rao [9].

Simulations
The simulations aimed at increasing the recovery of HA in the deslimed flow (cyclone underflow) while resulting in relatively small variations in the amount of fines in the same stream as compared with the Base Case.Simulated conditions included changes in the geometry of the cyclones, i.e., apex and vortex diameters, conical section inclination and changing the number of units in operation for adjusting the operating pressure.Table 8 shows the simulated changes as compared to the Base Case.The results obtained from the three simulations are presented in Table 9, together with the Base Case.
Compared with the Base Case, the three simulations significantly increased the flow rate of solids in the cyclone underflow.Accordingly, simulation 1 resulted in 146 t/h compared to 132 t/h of solids in the Base Case, an increase of 10.6%.Similarly, simulations 2 and 3 resulted in underflow solids flow rate increases of 12.9% and 16.7%, respectively.The same simulations also indicated moderate fines increases in the underflow, in this case, the fraction finer than 0.013 mm.

HA Simulations
The same conditions selected for the three simulations were used to simulate the HA distributions to cyclone' underflow and overflow in each case.Table 10 presents the results obtained for the HA distributions and the Base Case.
Table 10 indicates significantly higher HA partitions for the underflow relative to the Base Case.Considering the 66% HA partition figure as a reference in the Base Case, the three simulations resulted in HA partitions in the range of 72%-74%.The relative HA increases to underflow were 9.4%, 11.1% and 12.7%, respectively, for simulations 1, 2 and 3.
Table 11 shows the HA partition curve parameters for the Base Case and those calculated for the other three simulations.The respective HA partition curves were plotted based on these parameters, as shown in Figure 8.  Table 11 shows no significant differences among the d50c and α parameters of the HA partition curves calculated for the Base Case and the three additional simulations.However, the Graph in Figure 8 shows significant partition increases of HA in fractions of intermediate size, particularly in the −0.075 + 0.013 mm range.Conversely, simulated HA partitions for fractions coarser than 0.075 mm and finer than 0.013 mm virtually coincided with the Base Case.The combination of the simulated significant increase in HA in the −0.075 + 0.013 mm range and the high HA content in the same size range in the cyclone feed resulted in the simulated HA increase in the 9.4%-12.7%range.

Comparisons between Mass and HA Recoveries
Table 12 shows the absolute and relative global mass and P2O5 recoveries to the cyclone underflow for the Base Case and all three simulations.Table 12 shows more significant mass recoveries than P2O5 recoveries for all three simulations, thus indicating that the conventional approach of estimating cyclone performance based solely on mass recovery is not necessarily adequate for assessing metallurgical recovery.The HA simulated partitions based on calibrated equations to estimate partition curve parameters provided a suitable method for predicting the distribution of HA in each size fraction, as well as on the overall streams of the cyclone.
Further detailed simulations could have explored the reduction in the total fines of the cyclone underflow stream compared with the Base Case.Such exercises were not carried out as a comprehensive approach would include assessing the downstream flotation process, which was out of the scope of the present work.Table 11 shows no significant differences among the d50c and α parameters of the HA partition curves calculated for the Base Case and the three additional simulations.However, the Graph in Figure 8 shows significant partition increases of HA in fractions of intermediate size, particularly in the −0.075 + 0.013 mm range.Conversely, simulated HA partitions for fractions coarser than 0.075 mm and finer than 0.013 mm virtually coincided with the Base Case.The combination of the simulated significant increase in HA in the −0.075 + 0.013 mm range and the high HA content in the same size range in the cyclone feed resulted in the simulated HA increase in the 9.4%-12.7%range.

Comparisons between Mass and HA Recoveries
Table 12 shows the absolute and relative global mass and P 2 O 5 recoveries to the cyclone underflow for the Base Case and all three simulations.Table 12 shows more significant mass recoveries than P 2 O 5 recoveries for all three simulations, thus indicating that the conventional approach of estimating cyclone perfor-

Figure 1
Figure 1 also shows that the joined overflow from the initial classification stage, referred to as natural fines, is further classified in two parallel lines, each consisting of 16 cyclones of 254 mm in diameter, called the second desliming stage.The latter is the focus of the present work.The combined second desliming stage underflow is pumped to another classification stage, resulting in coarse and fine flotation feed circuits.In contrast, the combined second desliming stage overflow is transferred to the tailing dam.The optimal operation of the second desliming stage consists of efficiently separating the slimes with the maximum possible recovery of HA to the combined underflow.While the former is essential for a high flotation performance, the latter modulates the P2O5 recovery in the industrial circuit.Based on the overall plant feed, the mass split to the two parallel second desliming stages is rated at 15%.Due to practical and operational aspects, the present study focused exclusively on one of the two secondary desliming lines.As the equipment and operating conditions are identical, the results were considered replicable for either of the two lines.

Figure 1
Figure 1 also shows that the joined overflow from the initial classification stage, referred to as natural fines, is further classified in two parallel lines, each consisting of 16 cyclones of 254 mm in diameter, called the second desliming stage.The latter is the focus of the present work.The combined second desliming stage underflow is pumped to another classification stage, resulting in coarse and fine flotation feed circuits.In contrast, the combined second desliming stage overflow is transferred to the tailing dam.The optimal operation of the second desliming stage consists of efficiently separating the slimes with the maximum possible recovery of HA to the combined underflow.While the former is essential for a high flotation performance, the latter modulates the P 2 O 5 recovery in the industrial circuit.Based on the overall plant feed, the mass split to the two parallel second desliming stages is rated at 15%.Due to practical and operational aspects, the present study focused exclusively on one of the two secondary desliming lines.As the equipment and operating conditions are identical, the results were considered replicable for either of the two lines.

Figure 2 .
Figure 2. Method adopted in the present work.

Figure 2 .
Figure 2. Method adopted in the present work.

Figure 4 .
Figure 4. Size distributions of the cyclone feed, underflow and overflow.

Figure 4 .
Figure 4. Size distributions of the cyclone feed, underflow and overflow.

Minerals 2024 , 16 Figure 5 .
Figure 5. Corrected partition curve for overall flows around the cyclone nest.

Figure 5 .
Figure 5. Corrected partition curve for overall flows around the cyclone nest.

Figure 6 .
Figure 6.HA experimental and balanced data.

Figure 8 .
Figure 8. Base Case and simulated HA partition curves.

Figure 8 .
Figure 8. Base Case and simulated HA partition curves.

Table 2 .
Cyclone geometry and operating conditions.

Table 6 .
Selected parameters for the HA partition curve.

Table 6 .
Selected parameters for the HA partition curve.

Table 6 .
Selected parameters for the HA partition curve.

Table 12 .
Summary of overall mass and P2O5 recovery for the simulated scenarios.

Table 11 .
d50c e alfa parameters for the HA Base Case and simulations.

Table 12 .
Summary of overall mass and P 2 O 5 recovery for the simulated scenarios.