CFD-Based Structural Optimization of Rotor Cage for High-Efﬁciency Rotor Classiﬁer

: Due to the uneven materials dispersion and high dust concentration in industrial applications of turbo air classiﬁers, a high-efﬁciency rotor classiﬁer was designed. Numerical simulations by ANSYS-Fluent 19.0, the effects of rotor cage shape, the number of blades, and the blade proﬁle on the inner ﬂow ﬁeld, as well as classiﬁcation performance, were investigated. The simulation results indicated a signiﬁcant improvement in ﬂow ﬁeld distribution near the classiﬁcation surface with the conical rotor cage. Furthermore, there was an average reduction of 10.1% in cut size, as well as a 23.6% increase in classiﬁcation accuracy. When the number of blades was 36, the ﬂow ﬁeld distribution between the blades was relatively uniform and a smaller cut size was obtained at a higher classiﬁcation accuracy. A streamline blade with 52 ◦ as the inlet installation angle effectively reduced the impact of the airﬂow on the blade and eliminated the inertia anti-vortex between blades. The cut size reduction was 4.7–6.3%, with a basically unchanged classiﬁcation accuracy. The material classiﬁcation experimental results were in agreement with the simulated results. The discrete phase model (DPM) could well-predict the cut sizes and classiﬁcation accuracy, but it could not present the ﬁshhook effect. The present study provides theoretical guidance for the structural optimization of an air classiﬁer with a rotor cage.


Introduction
An air classifier is a key piece of equipment in the powder preparation process that is widely used in the mineral processing, food processing, and fine chemical industries [1][2][3][4]. The performance requirements for air classifiers are increasing with applications of ultra-fine powders. High-efficiency and energy-saving air classifiers are receiving widespread attention.
The turbo air classifier is a mainstream dynamic classifier for industrial applications that has a simple structure and controllable product granularity. However, with the surge in demand for powder production capacity, many problems have arisen in the industrial application of turbo air classifiers due to structural limitations [5,6]. On issue is a high dust concentration in the annular area; material classification is mainly concentrated in the annular area, so if a material cannot be sorted out in time, this would be a high dust concentration in the annular area, which would increase the chance of particle collision and reduce classification accuracy. A second problem is the uneven dispersion of materials; the number of feed inlets cannot be indefinitely increased with the increase in the demand for capacity, which has led to the uneven distribution of materials on the spreading plate and resulting in the formation of a thick material curtain, which is not conducive for material classification. A third problem is that the guide vanes can become severely worn and cannot be easily replaced. The materials in the annular area mainly rely on the impact guide blade to settle and be collected as coarse powder. The guide vanes wear out as a result and cannot be easily replaced due to structural limitations. In recent years, with impact guide blade to settle and be collected as coarse powder. The guide vanes wear o as a result and cannot be easily replaced due to structural limitations. In recent years, wi the rapid development of computational fluid dynamics (CFD), researchers have focuse more attention on the study of powder classification technology and the development classification equipment, and they have made progress in many respects. In order to solv these problems, many scholars at home and abroad have become devoted to the researc of optimizing the inner flow field distribution [7][8][9] and structurally improving the turb air classifier [10][11][12]. However, it should be noted that most studies have focused on loc optimization. The problems that have occurred in the industrial applications of turbo a classifiers have still not been positively solved. Since classification equipment is gradual moving towards high production levels and large scales, problems such as uneven mat rial spreading and dust concentration are becoming more obvious. Therefore, in order address these problems, a high-efficiency rotor classifier was designed. When the hig efficiency rotor classifier was used in combination with the static classifier, the abov problems were effectively avoided, thereby replacing the turbo air classifier. The structu is shown in Figure 1. A static classifier for the pre-classification of materials can effective reduce the content of coarse particles and decrease dust concentration in the classificatio chamber. Meanwhile, the wind feeding method allows for the material to be fully di persed. As the core component of the classifier, the structural parameters of the rotor cag directly affect the flow field distribution and classification performance of the classifie Hence, many researchers have focused on the rotor cage. Since material classification mainly concentrated near the classification surface, the flow field distribution at the e trance of the rotor cage has been investigated. According to the results, maintaining th stability of the flow field on the classification surface is conducive to improving the cla sification performance [13][14][15]. The present research mainly focused on turbo air classi ers. Hence, the flow field distribution inside the classifier with a downward air inlet r mains unclear. For a high-efficiency rotor classifier with a downward air inlet, furth study is needed to determine how to ensure the uniformity of airflow distribution ne the classifying surface. The rotor cage blade is the core part of the rotor cage, and it pla a crucial role in material classification. The blade pitch has a great influence on the pre sure drop and classification efficiency of the classifier. For instance, Yu et al. [16] analyze As the core component of the classifier, the structural parameters of the rotor cage directly affect the flow field distribution and classification performance of the classifier. Hence, many researchers have focused on the rotor cage. Since material classification is mainly concentrated near the classification surface, the flow field distribution at the entrance of the rotor cage has been investigated. According to the results, maintaining the stability of the flow field on the classification surface is conducive to improving the classification performance [13][14][15]. The present research mainly focused on turbo air classifiers. Hence, the flow field distribution inside the classifier with a downward air inlet remains unclear. For a high-efficiency rotor classifier with a downward air inlet, further study is needed to determine how to ensure the uniformity of airflow distribution near the classifying surface. The rotor cage blade is the core part of the rotor cage, and it plays a crucial role in material classification. The blade pitch has a great influence on the pressure drop and classification efficiency of the classifier. For instance, Yu et al. [16] analyzed the effect of the inner and outer rotor cage diameters on the flow field of a turbo air classifier. The results indicated that the inner and outer diameters of the rotor cage had a significant effect on the velocity distribution in the channel. The inner and Processes 2021, 9, 1148 3 of 22 outer diameters of the rotor cage should not be too large or too small. Jiang et al. [17] performed numerical simulations on rotor cages with different blade pitches in order to observe the flow field changes and energy losses in the classifier. It was found that there was an optimal blade pitch in the rotor cage, where the flow field in the classifier was the most stable and the classifying performance was the best. Liu et al. [18] combined their experiments with actual production experience, and five different blade numbers were selected for the numerical analysis. According to the research, the rough running phenomenon improved as the number of blades increased, and the coarse particle content gradually decreased. However, an excessive amount of blades would increase the chances of particle and blade collision, preventing the particles from being sorted out in time. Hence, a reasonable number of blades is conducive to improving the classification performance of air classifiers. As a necessary path in the material sorting process, the flow field distribution in the rotor cage channel directly determines classification efficiency [19,20]. In order to improve the flow field distribution in the rotor cage channel, some researchers have made improvements to the blade profile. The results have indicated that a positively bended arc blade can significantly improve the flow field distribution in turbo air classifiers, reducing the inertial rotating vortex and radial velocity fluctuation between the blades [21][22][23][24]. The Tianjin Cement Industry Design and Research Institute carried out numerical simulations of pressure loss and classification efficiency for blades widely used in the industry, and the researchers proposed a new type of U-shaped blade classifier with low pressure loss and high efficiency [25]. However, most previous studies have been limited to specific equipment, and the designed blades have not been universally applicable. On the contrary, when using the flow field distribution to determine the inlet and outlet installation angle, using the single arc method to design the blade can be suitable for different equipment.
In the present study, the numerical simulation of a high-efficiency rotor classifier for industrial applications was conducted with the aid of the ANSYS-Fluent 19.0 software. The effects of rotor cage shape, the number of blades, and the blade profile on the inner flow field, as well as the classification performance of the classifier, were investigated. According to the results, the conical rotor cage was more suitable for the classifier in the form of a downward air inlet. An optimum number of blades existed for the conical rotor cage. The designed streamlined blade could effectively improve the flow field distribution between the blades. Furthermore, the simulation results were verified by powder classification experiments. The present study not only designed a high-efficiency rotor classifier but also provided a new strategy for the design and optimization of rotor cages.

Structural Characteristics
The performance of an air classifier mainly depends on the two stages of material dispersion and classification, and the organic combination of these two is the key to optimize the design of such a classifier [26]. Hence, a high-efficiency rotor classifier was designed with the full consideration of material dispersion and classification. The structural characteristics are as follows: (1) The adequate dispersion of the material is a prerequisite to ensure high-efficiency classification. Therefore, the traditional approach of spreading the material was abandoned and replaced by wind feeding, and the material was carried into the highefficiency rotor classifier by airflow. The wind feeding method made the material fully mixed with the airflow, thus ensuring the dispersion of the material. (2) The wind feeding method decreased the coarse particle content and fully dispersed the material. This is conducive to the timely sorting out of material particles and the effective reduction of the dust concentration in the classifying chamber. (3) The material classification effect directly determines classification efficiency and energy consumption. Due to its stability and controllability, the centrifugal field was chosen as the classification force field.  (4) In the high-efficiency rotor classifier, the guide vanes were removed, and the conical rotor cage was chosen to strengthen and stabilize the centrifugal classification force field. The swirling flow and tangential shear force generated by the rotation of the conical rotor cage increased the material dispersion performance and improved the classification efficiency. (5) The spindle in a conventional turbo air classifier not only needed to drive the rotation of the rotor cage but also dispersed the material. As a comparison, the spindle in the high-efficiency rotor classifier only needed to drive the rotation of the rotor cage, which effectively reduced energy consumption.

Working Principle
The working principle for the high-efficiency rotor classifier is presented in Figure 2. After the milled material is dispersed and pre-classified by the static classifier, the coarse material is returned to the mill, while the fine material is carried into the high-efficiency rotor classifier with the airflow. As the airflow moves to the cone part, the flow velocity decreases due to the sudden increase in the cross-section of the channel. Coarse particles are collected as coarse powder, since these bounce back from the collision with the guide cone under inertia. The remaining material continues to be carried by the airflow through the guide cone rectification into the classifying chamber. In the classifying chamber, the material particles are mainly subjected to the air drag force and the centrifugal force generated by the rotor cage. The coarse particles are thrown out of the rotor cage due to the large centrifugal force, hit the wall, lose their kinetic energy, fall down, and are collected as coarse powder. The fine particles are carried by the airflow through the classifying surface into the rotor cage, discharged from the fine powder outlet, and collected as finished products.
(3) The material classification effect directly determines classification efficiency and energy consumption. Due to its stability and controllability, the centrifugal field was chosen as the classification force field.
(4) In the high-efficiency rotor classifier, the guide vanes were removed, and the conical rotor cage was chosen to strengthen and stabilize the centrifugal classification force field. The swirling flow and tangential shear force generated by the rotation of the conical rotor cage increased the material dispersion performance and improved the classification efficiency.
(5) The spindle in a conventional turbo air classifier not only needed to drive the rotation of the rotor cage but also dispersed the material. As a comparison, the spindle in the high-efficiency rotor classifier only needed to drive the rotation of the rotor cage, which effectively reduced energy consumption.

Working Principle
The working principle for the high-efficiency rotor classifier is presented in Figure 2. After the milled material is dispersed and pre-classified by the static classifier, the coarse material is returned to the mill, while the fine material is carried into the high-efficiency rotor classifier with the airflow. As the airflow moves to the cone part, the flow velocity decreases due to the sudden increase in the cross-section of the channel. Coarse particles are collected as coarse powder, since these bounce back from the collision with the guide cone under inertia. The remaining material continues to be carried by the airflow through the guide cone rectification into the classifying chamber. In the classifying chamber, the material particles are mainly subjected to the air drag force and the centrifugal force generated by the rotor cage. The coarse particles are thrown out of the rotor cage due to the large centrifugal force, hit the wall, lose their kinetic energy, fall down, and are collected as coarse powder. The fine particles are carried by the airflow through the classifying surface into the rotor cage, discharged from the fine powder outlet, and collected as finished products.

Geometry Modeling and Meshing
The modeling for the computational model of the high-efficiency rotor classifier was performed using SolidWorks. The main structural dimensions of the classifier are presented in Figure 3. The dimensions for the blades were 20 mm in length, 2 mm in thickness, and 150 mm in height, and the blades were evenly distributed along the outer edge of the rotor cage. In order to avoid backflow during the calculation, the central outlet pipe was extended up to three times of the diameter of the pipe. The computational model was divided into six regions, and each region was meshed using ICEM CFD in order to ensure the quality of the mesh. Except for the coarse powder outlet area, in which an unstructured mesh was used, a structured mesh was used for all other areas, and these areas were connected through the interface.

Geometry Modeling and Meshing
The modeling for the computational model of the high-efficiency rotor classifie performed using SolidWorks. The main structural dimensions of the classifier are sented in Figure 3. The dimensions for the blades were 20 mm in length, 2 mm in t ness, and 150 mm in height, and the blades were evenly distributed along the outer of the rotor cage. In order to avoid backflow during the calculation, the central outle was extended up to three times of the diameter of the pipe. The computational mode divided into six regions, and each region was meshed using ICEM CFD in order to e the quality of the mesh. Except for the coarse powder outlet area, in which an unstruc mesh was used, a structured mesh was used for all other areas, and these areas were nected through the interface. In order to verify the mesh independence, the average radial velocity of the ai at the entrance of the rotor cage with different mesh numbers was analyzed and pared. The results are presented in Table 1. When the number of meshes exceeded million, there was no significant difference in the average radial velocity. Hence, in sideration of computation time and accuracy, the 2.02 million mesh model was chose the simulation. The number of grids in each region is shown in Table 2. In order to verify the mesh independence, the average radial velocity of the airflow at the entrance of the rotor cage with different mesh numbers was analyzed and compared. The results are presented in Table 1. When the number of meshes exceeded 2.02 million, there was no significant difference in the average radial velocity. Hence, in consideration of computation time and accuracy, the 2.02 million mesh model was chosen for the simulation. The number of grids in each region is shown in Table 2.

Boundary Conditions and Parameter Settings
The boundary conditions were set by combining the structural characteristics and flow field characteristics of the high-efficiency rotor classifier. The SIMPLEC algorithm was adopted for pressure-velocity coupling, and the QUICK difference scheme was used for convective diffusion. The no-slip boundary conditions were used for the wall boundary, and the near-wall surface was treated with standard wall functions. The simulation working conditions were set as follows: an air volume of 1200 m 3 /h, rotation speed of 500 rpm, and feeding speed of 650 kg/h. The inlet was defined as a velocity inlet with a flow velocity of 13.8 m/s, assuming that the airflow was uniformly distributed in the inlet cross-section and direction perpendicular to the cross-section. The outlet was defined as a fully developed pipe flow (outflow). The multiple reference coordinates system model (MRF) was adopted to simulate the rotational motion of the rotor cage in the classification chamber, with the direction of rotation was defined as clockwise. The rotor cage blades were set to 'moving wall,' with a rotational speed of zero relative to the rotating area.
The discrete phase model (DPM) was adopted to simulate the trajectory of the particles. In order to simplify the model, the coarse powder outlet and sloping surface on the right side were set to 'trap.' This means that the calculation would be aborted if the material particles touched this surface, and they were collected as coarse powder. Considering the interaction between particles and the discrete vortex, the discrete random walk model (DRW) was used to simulate the dispersion effect of the turbulent flow on particles. Since material particles above 120 µm were already sorted out in the pre-classification stage, merely the classification efficiency of material particles within 0-120 µm needed to be discussed. In the velocity inlet for establishing 60 groups of surface injection source, the density of the material particles was set at 2750 kg/m 3 of calcium carbonate, and the particle size distribution was set to 0-120 µm. The tromp curves were obtained by tracking the particle trajectory, and the ratio of the number of particles trapped (trap) to the total number of particles tracked (track) was taken as the classification efficiency.

Turbulence Model
The three-dimensional steady simulation was performed using ANSYS-Fluent 19.0 (Canonsburg, PA, USA). For the case of incompressible flow, the mass and momentum equations are as follows: where ρ refers to the fluid density, u i refers to the fluid velocity, x i refers to the position, Reynolds stress, according to the law of time homogenization, in order to better predict the turbulent flow in the air classifier. The prediction of the partial classification efficiency and pressure drop also agreed well with the experimental data [2,27,28]. The transport equations for the transport of the Reynolds stresses may be written as follows: where ρ refers to the fluid density, u i refers to the fluctuating velocity to the direction i, x k refers to the positional length, D ij refers to the diffusion term, P ij refers to the stress production term, G ij refers to the buoyancy production term, ϕ ij refers to the pressure strain term, ε ij refers to the dissipation term, F ij refers to the production by system rotation term, and S user refers to the user-defined source term. Since the fluid was considered incompressible, Equation (3) with G ij = 0. The RSM model required the following empirical constants:

Discrete Phase Model
Since the volumetric amount of the inlet particles was below the traditional dilute limit of 0.1%, a one-way coupling was assumed [29]. Hence, the interaction between particles and the effect of particles on the gas phase were neglected. In the Lagrangian coordinate system, the DPM predicted the trajectory of the particle based on the forces on the particle, and the equation for motion was derived as follows: where m p refers to the particle mass, t refers to the time, refers to the drag force, and τ r refers to the relaxation time of the particle, which can be calculated using the following equation: where µ refers to the molecular viscosity of the fluid, d p refers to the particle diameter, and Re refers to the relative Reynolds number, which was derived as follows:

Experimental Setup
The scheme for the experimental setup is presented in Figure 4. The experimental setup mainly comprised the high-efficiency rotor classifier, cyclone collector, pulse bag filter, high-pressure induced draft fan, and control system. The analysis of the particle size distribution of the samples was performed using the Mastersizer 2000 laser particle sizer. The classifier pressure drop (∆P) was calculated by the average inlet static pressure minus the average outlet static pressure, which could be measured using a U-type manometer.
drop (ΔP), which were obtained from simulations and experiments. The characteri particle sizes, such as d50, d25, and d75, could be obtained by partial classification efficie curves. Then, K was calculated according to the formula K = d25/d75. A lower value for and a higher value for K indicate better classification performance.

Effect of the Rotor Cage Shape on Velocity Distribution near the Classification Surface
In order to clarify the distribution characteristics of the flow field in the classif multiple cross-sections were selected as the research object in the classifier, and veloc cloud diagrams and local vector diagrams of different cross-sections were obtained shown in Figure 5. It can be observed from the figure that the airflow entered from the inlet, entered the classification chamber through the guide cone, entered the rotor c under the action of the central negative pressure, and finally discharged from the f powder outlet. The velocity reached the maximum at the outlet of the fine powder, wh was favorable to the output of fine powder, while the velocity was smaller in the ne wall area of the classification chamber, which was favorable to the fall of the coarse po der. According to the velocity vector diagram, it can be observed that after the air fl was rectified by the guide cone, most of the air flow smoothly entered the rotor cage a a small part of the air flow continued to move in the axial direction. Inside the rotor ca the airflow moved upward at an angle of approximately 60°, causing part of the airfl to collide with the inner wall at the bottom of the rotor cage and deflect to form a vor This seriously affected the conveyance and classification of the material, which in t compressed the effective classification area. While outside the rotor cage, the airflow c tinued to move in the axial direction, colliding with the wall at the top of the classificat chamber and turning to form a vortex. Then, the vortex formed the secondary scour The evaluation of the performance of the high-efficiency rotor classifier was performed using parameters, such as cut size (d 50 ), classification accuracy (K), and pressure drop (∆P), which were obtained from simulations and experiments. The characteristic particle sizes, such as d 50 , d 25 , and d 75 , could be obtained by partial classification efficiency curves. Then, K was calculated according to the formula K = d 25 /d 75 . A lower value for d 50 and a higher value for K indicate better classification performance.

Effect of the Rotor Cage Shape on Velocity Distribution near the Classification Surface
In order to clarify the distribution characteristics of the flow field in the classifier, multiple cross-sections were selected as the research object in the classifier, and velocity cloud diagrams and local vector diagrams of different cross-sections were obtained, as shown in Figure 5. It can be observed from the figure that the airflow entered from the air inlet, entered the classification chamber through the guide cone, entered the rotor cage under the action of the central negative pressure, and finally discharged from the fine powder outlet. The velocity reached the maximum at the outlet of the fine powder, which was favorable to the output of fine powder, while the velocity was smaller in the near-wall area of the classification chamber, which was favorable to the fall of the coarse powder. According to the velocity vector diagram, it can be observed that after the air flow was rectified by the guide cone, most of the air flow smoothly entered the rotor cage and a small part of the air flow continued to move in the axial direction. Inside the rotor cage, the airflow moved upward at an angle of approximately 60 • , causing part of the airflow to collide with the inner wall at the bottom of the rotor cage and deflect to form a vortex. This seriously affected the conveyance and classification of the material, which in turn compressed the effective classification area. While outside the rotor cage, the airflow continued to move in the axial direction, colliding with the wall at the top of the classification chamber and turning to form a vortex. Then, the vortex formed the secondary scouring for the material after classification, reducing the chance for fine particles to be collected as coarse powder.
Processes 2021, 9, x FOR PEER REVIEW 9 for the material after classification, reducing the chance for fine particles to be collect coarse powder. In order to obtain a higher classification accuracy and efficiency, the "aerodyn screen" must maintain the same mesh size, that is, to ensure a stable and evenly di uted airflow near the classification surface [6]. According to the flow field distrib characteristics, since the high-efficiency rotor classifier was combined with the down air inlet form, it remained difficult to ensure that the velocity distribution near the c fication surface remained consistent along the axial direction. In order to determin suitable rotor cage shape for the high-efficiency rotor classifier, the velocity distribu near the classifying surface of the column rotor cage and conical rotor cage were ana and compared.
The outer surface of the rotor cage was taken as the object of study. The radia tangential velocity distributions of the airflow were obtained by taking points alon axial direction at the entrance of the rotor cage, as shown in Figure 6. The cloud dia for the radial and tangential velocity distributions corresponded to different rotor shapes, as shown in Figure 7. From the figure, it can be observed that the radial and gential velocities of the conical rotor cage decreased to different degrees, when comp to the column rotor cage, and that the velocity distribution was more uniform with sm gradients. The decreases in radial and tangential velocities were mainly correlated t rotor cage shape. Since the column rotor cage was combined with the downward air form, most of the airflow moved to the upper part of the rotor cage. Meanwhile, a vortex was formed at the bottom of the rotor cage, which made it difficult for the ai to evenly enter the rotor cage. Hence, the radial velocity of the airflow presented a d bution of upper large and lower small. Influenced by the vortex, the tangential veloc airflow also had a large velocity gradient. However, the conical rotor cage was comb In order to obtain a higher classification accuracy and efficiency, the "aerodynamic screen" must maintain the same mesh size, that is, to ensure a stable and evenly distributed airflow near the classification surface [6]. According to the flow field distribution characteristics, since the high-efficiency rotor classifier was combined with the downward air inlet form, it remained difficult to ensure that the velocity distribution near the classification surface remained consistent along the axial direction. In order to determine the suitable rotor cage shape for the high-efficiency rotor classifier, the velocity distributions near the classifying surface of the column rotor cage and conical rotor cage were analyzed and compared.
The outer surface of the rotor cage was taken as the object of study. The radial and tangential velocity distributions of the airflow were obtained by taking points along the axial direction at the entrance of the rotor cage, as shown in Figure 6. The cloud diagram for the radial and tangential velocity distributions corresponded to different rotor cage shapes, as shown in Figure 7. From the figure, it can be observed that the radial and tangential velocities of the conical rotor cage decreased to different degrees, when compared to the column rotor cage, and that the velocity distribution was more uniform with smaller gradients. The decreases in radial and tangential velocities were mainly correlated to the rotor cage shape. Since the column rotor cage was combined with the downward air inlet form, most of the airflow moved to the upper part of the rotor cage. Meanwhile, a large vortex was formed at the bottom of the rotor cage, which made it difficult for the airflow to evenly enter the rotor cage. Hence, the radial velocity of the airflow presented a distribution of upper large and lower small. Influenced by the vortex, the tangential velocity of airflow also had a large velocity gradient. However, the conical rotor cage was combined with the downward air inlet form, which increased the cross-sectional area at the entrance of the rotor cage due to the structural characteristics of large upper and small lower, resulting in a decrease in radial velocity. The vortex at the bottom of the rotor cage was weakened by the structural limitation, thus decreasing the tangential velocity at the bottom of the rotor cage. Since the conical rotor cage was more catered to the airflow direction, the airflow could evenly enter the rotor cage, making the velocity distribution more uniform. The tangential velocity of the airflow at the entrance of the rotor cage was small, primarily because the tangential velocity of the airflow was mainly generated by the rotation of the rotor cage. However, as the airflow entered the rotor cage channel, the tangential velocity gradually increased until it remained stable. Hence, compared with the column rotor cage, the conical rotor cage made the flow field near the classifying surface more evenly distributed and the classifying force field more stable. Therefore, material particles were basically subjected to the same force at any position near the classification surface, which was conducive to improving the classification performance.
with the downward air inlet form, which increased the cross-sectional area at the entrance of the rotor cage due to the structural characteristics of large upper and small lower, resulting in a decrease in radial velocity. The vortex at the bottom of the rotor cage was weakened by the structural limitation, thus decreasing the tangential velocity at the bottom of the rotor cage. Since the conical rotor cage was more catered to the airflow direction, the airflow could evenly enter the rotor cage, making the velocity distribution more uniform. The tangential velocity of the airflow at the entrance of the rotor cage was small, primarily because the tangential velocity of the airflow was mainly generated by the rotation of the rotor cage. However, as the airflow entered the rotor cage channel, the tangential velocity gradually increased until it remained stable. Hence, compared with the column rotor cage, the conical rotor cage made the flow field near the classifying surface more evenly distributed and the classifying force field more stable. Therefore, material particles were basically subjected to the same force at any position near the classification surface, which was conducive to improving the classification performance.  The results for the discrete phase simulation using the column and conical rotor cage at different operating conditions are presented in Table 3. For each operation condition, the cut size using the column rotor cage was smaller than that using the conical rotor cage, with a decrease of 9.2-11.3%. In addition, K increased by 13.0-25.9%. According to the simulation results, it could be observed that the conical rotor cage could effectively improve the classification performance and had a better adaptability to different working The results for the discrete phase simulation using the column and conical rotor cage at different operating conditions are presented in Table 3. For each operation condition, the cut size using the column rotor cage was smaller than that using the conical rotor cage, with a decrease of 9.2-11.3%. In addition, K increased by 13.0-25.9%. According to the simulation results, it could be observed that the conical rotor cage could effectively improve the classification performance and had a better adaptability to different working conditions.

Effect of Blade Number on Classification Performance
The number of rotor cage blades is the key factor that affects the flow field between blades. When the number of blades increased and made the rotor cage channel narrower, the vortex intensity between blades was weakened and the flow field became more stable. However, an excessive number of blades could lead to a very narrow channel in the rotor cage, which would affect the material passage and thereby reduce the processing capacity of the classifier. Therefore, it is essential to select a reasonable number of rotor cage blades in real applications.
In order to investigate the effect of the number of rotor cage blades on classification performance, five different rotor cages with N = 24, N = 36, N = 48, N = 60, and N = 72 blade numbers were simulated and the radial velocity distribution cloud was obtained, as shown in Figure 8. Due to the existence of vortices at the bottom of the rotor cage, only the cross-section within the effective classification area needed to be analyzed. Therefore, the velocity distributions at an XOY-plane with Z = −37.5mm and Z = −75 mm were discussed as follows. The results for the discrete phase simulation using the column and conical rotor cage at different operating conditions are presented in Table 3. For each operation condition, the cut size using the column rotor cage was smaller than that using the conical rotor cage, with a decrease of 9.2-11.3%. In addition, K increased by 13.0-25.9%. According to the simulation results, it could be observed that the conical rotor cage could effectively improve the classification performance and had a better adaptability to different working conditions.

Effect of Blade Number on Classification Performance
The number of rotor cage blades is the key factor that affects the flow field between blades. When the number of blades increased and made the rotor cage channel narrower, the vortex intensity between blades was weakened and the flow field became more stable. However, an excessive number of blades could lead to a very narrow channel in the rotor cage, which would affect the material passage and thereby reduce the processing capacity of the classifier. Therefore, it is essential to select a reasonable number of rotor cage blades in real applications.
In order to investigate the effect of the number of rotor cage blades on classification performance, five different rotor cages with N = 24, N = 36, N = 48, N = 60, and N = 72 blade numbers were simulated and the radial velocity distribution cloud was obtained, as shown in Figure 8. Due to the existence of vortices at the bottom of the rotor cage, only the cross-section within the effective classification area needed to be analyzed. Therefore, the velocity distributions at an XOY-plane with Z = −37.5mm and Z = −75 mm were discussed as follows. As shown in the figure, the radial velocity between blades was distributed in a fusiform stepped manner, and a large velocity gradient generally existed. Mainly due to the limited number of rotor cage blades, the airflow in the channel was not completely in accordance with the movement of the blade profile. There was a relative axial rotational motion in addition to the relative motion between blades in the radial direction, and these two movements compounded the relative motion of the airflow in the channel. On the side of the blade proceeding surface, two kinds of motion in the same direction superimposed on each other, making the relative velocity increase. On the side of the blade receding surface, these two motions were in opposite directions, and these canceled each other out, thus causing the relative velocity to decrease. When N < 36, the blade spacing was larger. Furthermore, in the receding surface and the rotor cage inlet, there were different degrees of anti-vortex, and the flow field between blades was unstable. When N > 36, the blade spacing gradually decreased as the number of blades increased. Furthermore, the radial velocity increased and the velocity gradient in the channel became larger. Meanwhile, the increased pressure difference between the blades made the airflow on the receding side reverse, forming an inertial anti-vortex in the channel. When the number of blades was increased to 72, the rotor cage channel narrowed, the space for airflow movement was reduced, and the inertial anti-vortex between the blades was basically eliminated. However, excessively narrow channels would lead to an increased chance of collision between materials and blades, which would not be conducive for classification. When N = 36, the velocity gradient between blades was small, the radial velocity direction was basically the same, and the flow field was relatively stable.
According to the flow field distribution, it was observed that the flow field between the blades was relatively stable and uniform when N = 36. For further validation, the partial classification efficiency of the classifier was calculated by the DPM model. At different operating conditions, the partial classification efficiency curves that corresponded to different numbers of blades were obtained, as shown in Figure 9. The cut size and classification accuracy were calculated according to the partial classification efficiency curve, as shown in Table 4. As shown in the figure, the radial velocity between blades was distributed in a fusiform stepped manner, and a large velocity gradient generally existed. Mainly due to the limited number of rotor cage blades, the airflow in the channel was not completely in accordance with the movement of the blade profile. There was a relative axial rotational motion in addition to the relative motion between blades in the radial direction, and these two movements compounded the relative motion of the airflow in the channel. On the side of the blade proceeding surface, two kinds of motion in the same direction superimposed on each other, making the relative velocity increase. On the side of the blade receding surface, these two motions were in opposite directions, and these canceled each other out, thus causing the relative velocity to decrease. When N < 36, the blade spacing was larger. Furthermore, in the receding surface and the rotor cage inlet, there were different degrees of anti-vortex, and the flow field between blades was unstable. When N > 36, the blade spacing gradually decreased as the number of blades increased. Furthermore, the radial velocity increased and the velocity gradient in the channel became larger. Meanwhile, the increased pressure difference between the blades made the airflow on the receding side reverse, forming an inertial anti-vortex in the channel. When the number of blades was increased to 72, the rotor cage channel narrowed, the space for airflow movement was reduced, and the inertial anti-vortex between the blades was basically eliminated. However, excessively narrow channels would lead to an increased chance of collision between materials and blades, which would not be conducive for classification. When N = 36, the velocity gradient between blades was small, the radial velocity direction was basically the same, and the flow field was relatively stable.
According to the flow field distribution, it was observed that the flow field between the blades was relatively stable and uniform when N = 36. For further validation, the partial classification efficiency of the classifier was calculated by the DPM model. At different operating conditions, the partial classification efficiency curves that corresponded to different numbers of blades were obtained, as shown in Figure 9. The cut size and classification accuracy were calculated according to the partial classification efficiency curve, as shown in Table 4.
the blades was relatively stable and uniform when N = 36. For further validation, the partial classification efficiency of the classifier was calculated by the DPM model. At different operating conditions, the partial classification efficiency curves that corresponded to different numbers of blades were obtained, as shown in Figure 9. The cut size and classification accuracy were calculated according to the partial classification efficiency curve, as shown in Table 4.    As shown in the chart, the cut size and classification accuracy had some fluctuations as the number of blades increased, indicating the existence of a better number of blades. When N < 36, the Tromp curve was flat, the classification accuracy was low, and the classifying effect of the classifier on material particles was poor. This was mainly due to the existence of an inertia anti-vortex between the blades, resulting in the phenomenon of back-mixing in the side of the receding surface. Hence, the material could not be classified in time. Meanwhile, due to the large blade spacing, the radial velocity in the channel was small. As the particle size increased, the effect of centrifugal force on the particles increased. Hence, the probability for material particles to enter the rotor cage and complete the classification decreased. When N > 36, with the increase in the number of blades, the classification efficiency of materials with a smaller particle size improved, but the classification efficiency of materials with a larger particle size rapidly decreased. Moreover, the greater the number of blades, the faster the decrease rate. Therefore, the classification accuracy decreased with the increase in the number of blades at different working conditions. This was mainly because as the number of blades increased, the blade spacing decreased, the radial velocity between blades increased, and the inertial anti-vortex strength became weakened. The fine particles were classified according to the better airflow carrying properties. Hence, the classification efficiency for the fine particles improved. However, a narrow rotor cage channel could lead to an increased chance of collision between materials and blades, which would increase the randomness of the material classification and decrease the classification efficiency of materials with a larger particle size. When N = 36, no polarization of the Tromp curve was observed, and a smaller cut size was obtained at higher classification accuracy, which was better for improving the classification performance and optimizing the particle grade.

Effect of the Blade Profile on the Flow Field Distribution between Blades
The movement of airflow in a rotor cage channel presents a typical forced vortex motion. On one hand, the airflow is driven by the rotor cage blade to do rotational movement; on the other hand, this moves inward under the action of the central negative pressure [30]. In this case, the classification of material particles mainly relies on the air drag force and centrifugal force, and both of these forces come from the airflow. Therefore, achieving a high classification accuracy is necessary to ensure that the flow field between blades is uniform and stable. By analyzing the flow characteristics of the airflow in the rotor cage channel, the influence of the blade profile on the flow field distribution was investigated, and a blade profile suitable for the high-efficiency rotor classifier was designed.
Due to the existence of a vortex at the bottom of the rotor cage, two cross-sections at the axial positions of Z = −37.5 and −75 of the rotor cage were chosen as the object of study. Figure 10 presents the cloud diagram for the radial velocity distribution and relative velocity streamline of airflow corresponding to the straight blade profile. The negative values of velocity in the figure indicate that the radial velocity direction pointed toward the center of the rotor cage, while the positive values of velocity indicate that the radial velocity direction pointed outward. As shown in Figure 10a, the flow field in the channel was seriously polarized, and a large velocity gradient was present and could reach up to −6-0 m/s. Hence, the coarse particles could easily enter the rotor cage through the side of the proceeding surface, and these were collected as fine powder. The fine particles were easily carried by the airflow to deflect and collide with the blades, making it difficult to pass through the rotor cage. This was mainly due to the airflow in the channel induced by the rotor cage blade forced action, which formed a negative pressure area in the receding surface. This made it easy to form an inertia anti-vortex between the blades. As shown in Figure 10b, the relative velocity of the airflow on the proceeding side was greater than that on the receding side at different axial positions of the rotor cage, and a partial inertial antivortex was formed on the side of the receding surface. Affected by the inertial anti-vortex and deflection, the airflow impacted the blade at a larger angle. This increased the chance for the material particles to hit the blades, which was not conducive for classification.
By analyzing the flow characteristics of the airflow in the rotor cage channel, the blade profile of the streamlined blade was determined. The movement of the airflow in the rotor cage channel exhibited a composite movement that could be decomposed into implicate and relative motion. The implicated velocity is represented by → u , the relative velocity is represented by → w, and the absolute velocity is represented by → v . According to these three velocity vectors, the velocity triangle could be drawn at the blade inlet and outlet, as shown in Figure 11. In the figure, α in refers to the inlet installation angle, i.e., the angle between the tangent line at the airflow inlet and the tangent line of the blade; β 1 refers to the relative velocity angle, which is the angle between the relative velocity vector and the tangent line at the airflow inlet; and i 1 refers to the impact angle, which is the angle between the tangent line of the blade and the relative velocity vector. chance for the material particles to hit the blades, which was not conducive for classification. By analyzing the flow characteristics of the airflow in the rotor cage channel, the blade profile of the streamlined blade was determined. The movement of the airflow in the rotor cage channel exhibited a composite movement that could be decomposed into implicate and relative motion. The implicated velocity is represented by ⃗ , the relative velocity is represented by ⃗⃗ , and the absolute velocity is represented by . According to these three velocity vectors, the velocity triangle could be drawn at the blade inlet and outlet, as shown in Figure 11. In the figure, refers to the inlet installation angle, i.e., the angle between the tangent line at the airflow inlet and the tangent line of the blade; 1 refers to the relative velocity angle, which is the angle between the relative velocity vector and the tangent line at the airflow inlet; and 1 refers to the impact angle, which is the angle between the tangent line of the blade and the relative velocity vector. Since a conventional rotor cage has a straight blade profile, the inlet installation angle and outlet installation angle were both 90°. The relative velocity angle was calculated using Fluent 19.0, in which the relative velocity angle 1 at the inlet and 2 at the Since a conventional rotor cage has a straight blade profile, the inlet installation angle α in and outlet installation angle α out were both 90 • . The relative velocity angle was calculated using Fluent 19.0, in which the relative velocity angle β 1 at the inlet and β 2 at the outlet were 38 • and 100 • , respectively. Thus, the airflow impact angle i 1 at the inlet and i 2 at the outlet were determined to be 52 • and −10 • , respectively. According to the energy equation, it is clear that the greater the impact angle of airflow, the greater the resistance and energy loss. Therefore, the inlet angle of the streamlined blade was set to 38 • in order to reduce the impact of airflow on the blade. In addition, the outlet installation angle was set to 90 • , because the airflow in the blade end of the impact angle was small and the impact on the blade could be negligible. Meanwhile, the deflection of the blade profile could be reduced to facilitate the processing and manufacturing. According to the blade inlet and outlet installation angle, the streamlined blade was drawn using the single arc method, as shown in Figure 12. Since a conventional rotor cage has a straight blade profile, the inlet installation angle and outlet installation angle were both 90°. The relative velocity angle was calculated using Fluent 19.0, in which the relative velocity angle 1 at the inlet and 2 at the outlet were 38° and 100°, respectively. Thus, the airflow impact angle 1 at the inlet and  The cloud diagram for the radial velocity distribution and relative velocity streamlines corresponding to the streamlined blade is shown in Figure 13. In Figure 13a, it can be observed that the velocity gradient in the flow channel was significantly reduced and that the radial velocity basically stabilized from −5 to −1 m/s. Furthermore, the flow field distribution was more uniform and the inertial anti-vortex disappeared. In Figure 13b, it can be observed that the streamlines in the different axial positions of the rotor cage were basically coincident with the blade profile. The airflow impact angle at the inlet was basically 0 • , which reduced the impact of the airflow on the blade. The airflow at the outlet was also in good agreement with the blades, and there were no large deflection occurrences. The improvement of the inner flow field was mainly because the streamlined blade profile was more compatible with the flow characteristics of the airflow, reducing the aerodynamic resistance. This would be beneficial for the conveyance and classification of material particles, reducing the chance of particle back-mixing.
According to the discrete phase simulation results, d 50 and K are illustrated in Table 5. For each operating condition, the simulation results indicated that the streamlined blade rotor cage reduced the cut size by 4.7-6.3%, while the classification accuracy remained basically unchanged. The cut size reduction was mainly due to the flow field in the rotor cage channel, which became more uniform, and the velocity gradient decreased. As a result, the drag force of the particles decreased and materials with a larger particle size had difficulties in passing through the rotor cage. However, a range of application exists for the inlet installation angle at simulated operating conditions. When the inlet air velocity and rotational speed greatly differed from the simulated working condition, the inlet installation angle and airflow incidence angle were no longer equal, which was not conducive for classification. Therefore, the application range of the streamlined blade rotor cage needs to be further studied.
basically coincident with the blade profile. The airflow impact angle at the inlet was basically 0°, which reduced the impact of the airflow on the blade. The airflow at the outlet was also in good agreement with the blades, and there were no large deflection occurrences. The improvement of the inner flow field was mainly because the streamlined blade profile was more compatible with the flow characteristics of the airflow, reducing the aerodynamic resistance. This would be beneficial for the conveyance and classification of material particles, reducing the chance of particle back-mixing. Figure 13. The cloud diagram for the radial velocity distribution and relative velocity streamline corresponding to the streamlined blade.
According to the discrete phase simulation results, d50 and K are illustrated in Table  5. For each operating condition, the simulation results indicated that the streamlined blade rotor cage reduced the cut size by 4.7-6.3%, while the classification accuracy remained basically unchanged. The cut size reduction was mainly due to the flow field in the rotor Figure 13. The cloud diagram for the radial velocity distribution and relative velocity streamline corresponding to the streamlined blade.

Comparison of the Simulation and Experiment Results
In order to verify the reliability of the simulation results, powder classification experiments were conducted and the simulation results were compared with the experimental data. An optimized rotor cage (a conical rotor cage with 36 straight blades) was selected for the classification experiments. The experiment and simulation were performed at different inlet velocities and rotary speeds, including 13.8-500, 13.8-600, and 16.2-500. The particle size distribution when using calcium carbonate as the raw material is presented in Figure 14. The blue line represents the cumulative distribution curve, while the red line represents the distribution frequency of the particles. After each test, the coarse and fine powders were weighed for mass balance. The samples were analyzed using a laser particle sizer, and then the partial classification efficiency of the high-efficiency rotor classifier was calculated.

Comparison of the Simulation and Experiment Results
In order to verify the reliability of the simulation results, powder classification experiments were conducted and the simulation results were compared with the experimental data. An optimized rotor cage (a conical rotor cage with 36 straight blades) was selected for the classification experiments. The experiment and simulation were performed at different inlet velocities and rotary speeds, including 13.8-500, 13.8-600, and 16.2-500. The particle size distribution when using calcium carbonate as the raw material is presented in Figure 14. The blue line represents the cumulative distribution curve, while the red line represents the distribution frequency of the particles. After each test, the coarse and fine powders were weighed for mass balance. The samples were analyzed using a laser particle sizer, and then the partial classification efficiency of the high-efficiency rotor classifier was calculated. The partial classification efficiency (Tromp curve) for the simulation and experiment is presented in Figure 15. A certain difference between the simulated and experimental Tromp curves can be observed in the figure, mainly in the fishhook effect, but the overall trend was the same. The experimentally obtained Tromp curves exhibited the same S- Figure 14. Particle size analysis report for the calcium carbonate raw material.
The partial classification efficiency (Tromp curve) for the simulation and experiment is presented in Figure 15. A certain difference between the simulated and experimental Tromp curves can be observed in the figure, mainly in the fishhook effect, but the overall trend was the same. The experimentally obtained Tromp curves exhibited the same S-shaped trend, presenting a fishhook effect. As the particle size increased, the partial classification efficiency curve presented a trend of initially decreasing and subsequently increasing, while the fishhook dip appeared at approximately 6 µm. This phenomenon could not be predicted by the DPM model, which was due to the neglect of particle interactions during the simulation, resulting in a higher classification accuracy for the simulation compared to that of the experiment. A comparison between the simulated results and the experimental data is presented in Table 6. It can be observed that the cut size and classification accuracy obtained from the simulation were in better agreement with the experiment, with maximum errors of 9.4% and 9.8%, respectively. In addition, the variation of the pressure drop was well-predicted, with a maximum error of 6.0%. Since the relative errors were all within 10%, the computational model was correct and reasonable for predicting the performance of the classifier. However, the effect of the streamlined blades on the classification performance of the high-efficiency rotor classifier needs to be further verified. The application range of streamlined blades, as well as the effect on the fishhook effect and classification accuracy, are worthy of further studies. for predicting the performance of the classifier. However, the effect of the streamlined blades on the classification performance of the high-efficiency rotor classifier needs to be further verified. The application range of streamlined blades, as well as the effect on the fishhook effect and classification accuracy, are worthy of further studies.

Conclusion
In the present study, a high-efficiency rotor classifier was designed, and the flow behavior in the classifier was numerically simulated by the ANSYS-Fluent 19.0 software. The effects of rotor cage shape, the number of blades, and the blade profile on the inner

Conclusions
In the present study, a high-efficiency rotor classifier was designed, and the flow behavior in the classifier was numerically simulated by the ANSYS-Fluent 19.0 software. The effects of rotor cage shape, the number of blades, and the blade profile on the inner flow field distribution and the classification performance of the classifier were investigated. The conclusions are as follows: (1) The high-efficiency rotor classifier employed wind feed, which reduced the content of coarse particles and increased the dispersion of materials. This effectively solved the problems of uneven materials spreading and high dust concentration. Meanwhile, since no guide vanes were equipped in the high-efficiency rotor classifier, there was no need to consider the wear of the guide vanes. (2) According to the analysis of the flow field characteristics, it could be seen that the vortices existed at the bottom of the rotor cage and the column rotor cage was not conducive to maintaining the uniformity of airflow distribution near the classification surface. On the contrary, the conical rotor cage not only reduced the vortex at the bottom of the rotor cage but also improved the uniformity of the flow field. (3) As the number of blades increased, the cut size and classification accuracy exhibited some fluctuations, indicating the existence of a better number of blades. A reasonable number of blades resulted in a more uniform flow field distribution between the blades and no polarization of the Tromp curve. This contributed to the improvement of classification performance and the optimization of particle gradation. (4) Compared with a traditional straight blade, a streamlined blade was found to eliminate the inertial anti-vortex between blades, reducing the velocity gradient and making the flow field distribution more uniform. Meanwhile, the streamline of airflow in the channel basically coincided with the blade profile, which reduced the impact of the airflow on the blades. This is beneficial for the conveyance and classification of materials. (5) The experimental results indicated that the DPM could well-predict the d 50 and K values of the high-efficiency rotor classifier. However, the experimental partial classification efficiency presented an S-shape, known as 'fishhook,' which is a phenomenon that cannot be predicted.