The main parameters that affect the heat transfer effect of vertical sinter cooling furnace include operating parameters and structural parameters. The operating parameters mainly include cooling air inlet temperature, sinter inlet air temperature equivalent diameter, and the structural parameters include cooling section height and cooling section diameter.
3.1. Effect of Cooling Air Inlet Temperature on Gas–Solid Heat Transfer Characteristics
The heat transfer process in the furnace was numerically simulated by changing only the inlet air temperature. The inlet air temperature parameters are shown in
Table 2, in which the inlet temperature is 750 °C, the height of the cooling section is 7 m, the diameter of the cooling section is 11.25 m, and the equivalent diameter of the sinter is 14.45 mm.
In order to obtain the effect of sinter equivalent diameter on the temperature field in the furnace, it is necessary to analyse the axial and radial temperature distribution of the cooling air in the furnace with the equivalent diameter.
The central section of sinter vertical cooling furnace is selected as the characteristic section to reflect the change in axial cooling air temperature in the furnace. The z = 3.5 m section is the characteristic section, which reflects the change in radial cooling air temperature in the furnace.
As can be seen from
Figure 5, the temperature of the cooling air in the central section of the furnace gradually becomes larger when the inlet temperature of the cooling air becomes larger. This is due to the fact that with the increase in the inlet air temperature, when other parameters in the furnace remain unchanged, the heat exchange effect between the cooling air and the hot sintered ore is enhanced, so that the temperature of the cooling air also increases. In addition, the temperature of the cooling air is lower in the furnace near the wall and at the hood. This is due to the high flow rate of the cooling air at the hood and the shorter heat exchange time with the hot sinter, so the temperature of the cooling air at this place is lower compared with that in the centre of the furnace. As the sintered ore particles near the wall are larger, and there are more gaps between the particles, which is conducive to the flow of cooling air, the contact area between the cooling air and the sintered ore is larger, and the heat exchange effect is stronger, so the temperature of the cooling air in this area is also lower.
Figure 6 shows that the cooling air temperatures show a trend of higher temperatures in the centre part than in the surroundings in all this cross-section, and with the increase in the inlet air temperature, the temperature of the cooling air in this cross-section also increases. The main reason is that the sintered ore in the centre part of the furnace is more closely stacked, and the sintered ore particles in this part are smaller and at higher temperatures, while the smaller gaps between the particles of this part of the sintered ore lead to a higher cooling air velocity in this part, and the heat exchange between the cooling air and the sintered ore is stronger in this region, and the temperature of the cooling air is also higher. In addition, with the increase in inlet air temperature and the constant flow rate of cooling air inlet, the apparent flow rate of gas in the furnace increases. This phenomenon makes the time required for the sinter at the same temperature inlet to be cooled to a certain temperature decrease continuously. When the sinter moves at a certain speed, the contact time between the cooling wind and the sinter increases, so the temperature of the cooling wind also increases successively.
A radial straight line is taken for a certain height in the furnace to reflect the variation law of each parameter of the section at the same height in the furnace. In order to clearly observe the changes of the section at the same height and select the section with drastic parameter changes for analysis, the line graph with z = 1 (the radial direction is taken at 1 m above the hood) is selected for data processing.
Figure 7 shows the variation in cooling air temperature along the axial direction at different cooling air inlet temperatures. As can be seen from
Figure 5, the temperature of the cooling air in the pre-storage section under different cooling air inlet temperature is basically unchanged, because in the pre-storage section of the cooling air flow is extremely slow, this part of the basic is not involved in the heat transfer, so the temperature of the cooling air in the pre-storage section of the basic remains unchanged in each operating condition. In the main cooling section from −2 m to 2 m, the heat exchange between the sinter and the cooling air is the most intense, and the heat exchange effect of this part is greatly affected by the inlet air temperature. As the inlet air temperature increases, the temperature difference between the cooling air and the sinter decreases, and the heat exchange efficiency decreases. Therefore, when the inlet air temperature of the cooling air rises, the temperature of the cooling air in the furnace is also higher. The section from 2 m to 7 m is the secondary heat exchange area, and the temperature of sintered ore in this area is higher. With the increase in inlet wind temperature, the heat exchange effect between sintered ore and cooling wind becomes weaker, the heat exchange coefficient also decreases, and the inlet wind temperature has little influence on the cooling wind temperature in this part. From the figure, it can be seen that the temperature of the cooling air in this part also basically does not change. Starting from −2 m is the exit section of the sintered ore; the cooling wind is basically not involved in heat exchange in this part, so the cooling wind temperature in this area is not affected by the inlet wind temperature, but this area is close to the inlet wind cap. So, with the increase in the inlet wind temperature, the temperature of the cooling wind in this area is also higher, but basically unchanged within the same working condition.
Figure 8 shows the radial variation in cooling air temperature under different cooling air inlet temperatures. In z = 1 m cross-section of the furnace cooling air temperature along the radial first increased and then decreased, showing the shape of the “M” type distribution. The distribution of cooling air temperature along the radial direction at the same cross-section is not much affected by the working condition parameters, basically showing a low temperature in the central part of the chamber. The temperature on both sides of the air temperature rises, and the temperature of the cooling air immediately adjacent to the wall is decreasing. The reason is that the inlet air cap is cone-shaped, so the cooling wind enters the furnace through the cross duct, resulting in the centre part of the cooling wind temperature is small. And, with the cooling wind from the hood to the surrounding flow process and sintered ore cooling heat exchange, in the flow process of cooling, wind temperature increases and sintered ore is cooled. Near the wall, due to the wall effect, the cooling air temperature is slightly reduced at the wall.
Table 3 shows the values of sintered ore and cooling air outlet parameters under different operating conditions, and the results are analysed below.
Figure 9 shows the variation in cooling air and sinter outlet temperature for different sinter inlet air temperatures. When the inlet temperature of the sintered ore is 750 °C, the inlet velocity of the cooling air is 7.8 m/s, the height of the cooling section is 7 m, and the diameter of the cooling section is 10 m, the outlet temperatures of both the sintered ore and the cooling air increase with the increase in the inlet air temperature, and the increase in the outlet temperature of the sintered ore has a linear distribution, while the increase in the outlet temperature of the cooling air decreases gradually. The reason is that when the cooling air inlet speed is unchanged, the influence of gas–solid heat transfer in the furnace is strengthened as the cooling air inlet temperature increases, so that the temperature at the cooling air outlet increases. Meanwhile, the sinter outlet temperature increases at the same rate as the cooling air inlet temperature increases.
Figure 10 shows the increase in exergy and the variation in exergy efficiency of the cooling air in the cooling furnace for different cooling air inlet air temperatures. As the cooling wind inlet wind temperature increases, the temperature at the cooling wind outlet also increases, which results in a larger heat exchange between the wind and the mine. At the same time, the speed of the inlet cooling wind is constant, and the heat exchange effect in the furnace increases, resulting in the exergy carried by the cooling wind at the exit also increases. The exergy efficiency is affected by the increase in cooling wind exergy and sinter ore inlet exergy as shown in the formula. When the inlet temperature of the sintered ore is stable, the increase in the inlet air temperature leads to the increase in the exergy of the cooling air at the outlet, and then the exergy efficiency also increases.
3.2. Influence of Inlet Air Temperature on the Heat Exchange Effect of Vertical Sinter Ore Cooling Furnace
Numerical simulation of the heat transfer process in the furnace is carried out by changing the inlet ore temperature. The parameters of the inlet mine temperature are shown in
Table 4, where the inlet air temperature is 20 °C, the height of the cooling section is 7 m, the diameter of the cooling section is 11.25 m, and the equivalent diameter of the sintered ore is 14.45 mm.
Figure 11 shows that with the increase in the sinter ore inlet temperature, the temperature of the sinter ore in the central section of the furnace gradually becomes larger. The reason is that with the increase in the inlet ore temperature, the cooling air and the hot sintered ore also undergo intense heat exchange, and the heat exchange effect is enhanced. However, the contact time between the sintered ore and the cooling air is constant due to the other parameters in the furnace remaining unchanged, which makes the sintered ore discharged before it is completely cooled. As a result, the temperature of the sinter in the central section of the furnace increases as the sinter inlet temperature increases. At the same time, as the sinter inlet temperature increases, the temperature of the sintered ore near the wall is lower and the temperature of the sintered ore in the furnace tends to be the same. This is because the higher temperature inlet sintered ore is not cooled to the same extent in the same time, so the temperature distribution in the furnace is not consistent.
Figure 12 shows the variation in cooling air temperature with sinter ore inlet temperature at z = 3.5 m section in the furnace, and the following analyses are made. The temperature of the sintered ore at the same height in the furnace increases gradually with the increase in the inlet temperature of the sintered ore. This is also because the high-temperature sintered ore is not sufficiently cooled in the furnace, and the higher the inlet temperature is, the heat exchange within the same cooling time is basically the same, so the overall temperature of the sintered ore in the furnace also increases. In addition, the temperature of the sintered ore near the wall is lower than the temperature of the sintered ore in the centre of the furnace.
Figure 13 shows the variation in sinter ore temperature along the axial direction in the furnace for different sinter ore inlet temperatures. With the increase in the sinter ore inlet temperature, the overall trend of the sinter ore temperature in the furnace increases, but the inlet ore temperature does not have much influence on the heat exchange process in the furnace. The temperatures of cooling air in the pre-storage section and sinter outlet section are basically unchanged under the same operating conditions. The reason is that the simulation process treats the high-temperature sinter as an isotropic porous medium, so when the inlet boundary value of the sinter temperature is changed, the overall heat transfer in the furnace is also enhanced with the increase in the inlet mine temperature, which makes the temperature of the cooling air at the same position also increase, but it remains unchanged under the same working condition.
Figure 14 shows the radial variation in sinter temperature at z = 1 m with different sinter inlet temperatures, and the following analyses are made. With the increase in sinter inlet temperature, the sinter in the furnace at the same height along the radial distribution of double “hump” type, basically flat at the cap, in the middle of the furnace diameter to reach the maximum, and near the wall and gradually decreased.
Table 5 shows the values of sinter and cooling air outlet parameters under different schemes.
Figure 15 shows the trend of cooling wind and sinter outlet temperature with the sinter inlet temperature. When the inlet temperature of the cooling air is 20 °C, the inlet air velocity of the cooling air is 7.8 m/s, the height of the cooling section and the diameter of the cooling section are 7 m and 10.0 m, respectively, the outlet temperatures of the cooling air and the sintered ore both increase with the increase in the sintered ore inlet temperature. The reason is that when the cooling wind flow rate is constant, the cooling time of the cooling wind on the sintered ore is also relatively constant. The sinter that is not completely cooled within the same time is discharged, which increases the temperature of both the sinter and the cooling air at the outlet. In addition, the rate of increase in the outlet mine temperature becomes faster with the increase in the inlet mine temperature, while the rate of increase in the outlet temperature of the cooling air and the change in the inlet temperature of the sintered ore are basically the same, and both of them show a linear relationship.
Figure 16 shows the amount of exergy change and the change in energy efficiency of the cooling air in the furnace at different sinter ore inlet temperatures. Both the exergy increase and the energy efficiency of the cooling air in the furnace increase with the increase in the sinter inlet temperature, and basically show a linear distribution. The reason is that the cooling air outlet temperature increases linearly with the increase in sinter inlet temperature, and the exergy increases with the increase in sinter inlet temperature after the cooling air absorbs heat, and the increase is basically the same, as shown in the equation of exergy calculation.
3.3. Effect of Particle Equivalent Diameter on Gas–Solid Heat Transfer Characteristics
The heat transfer process in the furnace is numerically simulated by changing the sinter equivalent diameter. The parameters of the sinter equivalent diameter were selected as shown in
Table 6, in which the inlet air temperature was 20 °C, the inlet ore temperature was 750 °C, the height of the cooling section was 7 m, and the diameter of the cooling section was 11.25 m.
Figure 17 illustrates the cooling air temperature distribution with sinter equivalent diameter in the centre profile of the shaft furnace. The sinter equivalent diameter has little influence on the overall heat transfer process in the furnace. The high-temperature sintered ore enters the furnace from the roof inlet and moves slowly from top to bottom, while the cooling air enters the furnace from the air cap and moves from bottom to top. In this process, the sintered ore moves to the cooling section and starts to be cooled, and the cooling wind absorbs heat and is gradually heated, and changing the equivalent diameter of the sintered ore in the process of gas–solid heat transfer does not have much effect on the heat transfer process in the furnace, so the trend of change in the figure is not obvious. The temperature of sintered ore in the centre part of the furnace tends to be the same, while the temperature near the wall decreases, because the sintered ore particles contact the wall of the cooling furnace, the hot sintered ore not only exchanges heat with the cooling air, but also exchanges heat with the wall at lower temperature, so the temperature of the sintered ore near the wall is lower.
Figure 18 demonstrates the variation in cooling air temperature with equivalent diameter at the 3.5 m cross-section in the furnace. The cooling air temperature in the furnace does not change significantly at this cross-section, and the diameter of the sinter particles does not have much influence on the heat exchange process in the furnace. The temperature of the cooling air is the highest in the centre of the furnace, and gradually decreases along the radial direction to the surrounding area. The reason is that the cooling air along the radial flow to the surrounding process has been and hot sinter ore heat transfer, and basically in thermal equilibrium, the temperature and the sinter ore temperature is basically the same, the temperature is the lowest near the wall.
Figure 19 shows the variation in cooling air temperature along the axial direction in the furnace for different sinter equivalent diameters, and the trend in the figure is analysed. With the increase in sinter equivalent diameter, the overall change in cooling air temperature in the furnace is not obvious, and the temperature of cooling air in the pre-storage section and sinter outlet section is basically unchanged under different working conditions. The reason is that the heat exchange effect in the furnace remains unchanged for a certain sinter inlet temperature and cooling air inlet temperature, which makes the effect of changing the sinter equivalent diameter on the temperature field in the furnace insignificant.
Figure 20 shows the radial variation in cooling air temperature with different equivalent diameters, and the trend is analysed in the figure. With the increase in the equivalent diameter of the sinter ore, the cooling air temperature in the cross-section along the radial direction of the middle of the low sides of the high class “M” type distribution. However, with the equivalent diameter from 18.80 mm to 27.30 mm, the temperature of the cooling air in the cooling furnace at the very centre of the cooling air is higher, because with the increase in the equivalent diameter of the sintered ore, the sintered ore and the cooling air at the part of the air cap heat transfer is the most intense, and at the same time, in the area nearer to the centre of the furnace, the lower the porosity of the sintered ore, the higher the flow rate of the cooling air, the convection in the region of the best heat transfer, so there is a higher temperature.
Table 7 shows the values of sinter and cooling air outlet parameters under different schemes.
Figure 21 demonstrates the variation in outlet cooling wind and sintered ore temperature with equivalent diameter. As can be seen from
Figure 19, when the cooling air inlet temperature is 20 °C, the cooling air inlet wind speed is 7.8 m/s, the height of the cooling section and the diameter of the cooling section are 7 m and 10.0 m, respectively, the outlet wind temperature decreases with the increase in the equivalent diameter of the sintered ore and the cooling air temperature decreases at a faster rate; the outlet mine temperature also increases with the increase in the equivalent diameter of the sintered ore and the magnitude of the increase in temperature becomes greater with the increase in the equivalent diameter of the sintered ore. Increases in the rate of temperature increase with the equivalent diameter of the sintered ore. This is because the expansion of the diameter of the sintered ore particles, the cooling air and sintered ore contact area becomes larger, resulting in a convective heat transfer coefficient with the expansion of the equivalent diameter of the particles, and as the gas–solid heat transfer becomes smaller, the export temperature of the ore rises and the export air temperature decreases.
Figure 22 demonstrates the amount of exergy change and the change in exergy efficiency of the cooling air in the cooling furnace for different equivalent diameters. The increase in cooling air exergy decreases as the equivalent diameter increases. The reason is that the change in equivalent diameter of sintered ore makes the stacking density of particles change, and the larger the equivalent diameter is, the smaller the stacking density becomes, which makes the effect of heat exchange between the sintered ore and the cooling air increase. When the particle size decreases when the cooling air and sintered ore heat transfer is more gentle, the sintered ore heat transfer time is relatively longer, so is the equivalent diameter of the cooling air when the change in exergy is smaller, and expanding the equivalent diameter increases the cooling air exergy changes.
3.4. Establishment of Gas–Solid Two-Phase Heat Transfer Correlation in a Furnace
In order to study the gas–solid heat transfer process in vertical cooling furnace, the dimensionality of heat transfer coefficient is analysed by the Π theorem.
According to the theoretical analysis, the factors affecting the heat transfer coefficient
h in the furnace are air density, air flow rate
u, diameter of sinter particles
d, air viscosity
, specific pressure heat capacity
Cp, sinter thermal conductivity
, etc., then there is the following relationship between the formula:
The seven variables mentioned above are all composed of the basic measures mass M, length L, time T and temperature Θ, so d, , , u are chosen as the basic measures, and the rest of the variables can be expressed by the four basic measures.
According to the Π theorem, the number of variables minus the number of basic measures is the number of Π-numbers, and three Π-numbers, Π
1, Π
2 and Π
3, can be obtained, then Equation (11) can be expressed as:
The quantitative representation of all variables is shown in
Table 8.
(1) Π1 = =
Then, the dimensionless form of Π
1 is:
(2) Π2 = =
Then, the dimensionless form of Π
2 is:
(3) Π2 =
Then, the dimensionless form of Π
3 is:
It can be known from fluid mechanics that
Nu =
,
Re =
,
Pr =
, so Π
1 =
Nu, Π
2 =
Re, Π
3 =
Pr; therefore, the following dimensionless equation is obtained:
The bed porosity is also a major factor affecting the heat exchange process between the sinter particles and the cooling air in the furnace as analysed in the previous section, so the dimensionless number is added to the equation of magnitude. In addition, the temperature difference between sinter and cooling air obtained from the simulation calculation exists in the cooling section in the region of more than 50 °C. At the same time, the cooling section is the most important gas–solid exchange process. At the same time, the cooling section is the most important gas–solid heat transfer area, so it is necessary to consider the change in the relevant physical properties in the heat transfer correlation equation, and it is not possible to distinguish the influence of the change in different physical properties on the heat transfer process only by the difference of the
Pr number index. Therefore, a correction coefficient
Ct [
25] on temperature is introduced with the following expression:
Combining the above analyses, the following dimensionless equations can be obtained:
By a process similar to the analysis of the bed pressure drop correlation equation, it can be seen that Equation (15) is also written in the form of a power function:
It can be shown that Equation (16) is the initial form of the correlation equation for the heat transfer of the flow within the bed of sintered ore particles obtained by using the method of magnitude analysis, where k, a, b, c, and d are coefficients to be determined.
The pending coefficients of the dimensionless equations were obtained from the data obtained from the design conditions using multiple regression analysis, and fitted based on Matlab 2022. The pending coefficients of the fitted equation are shown in
Table 9.
The empirical coefficients were substituted into the heat transfer correlation equation, which was fitted to obtain the correlation equation for heat transfer in the furnace:
The conditions of use for this bed pressure drop correlation are as follows: , the range of cooling air inlet temperature is 20 °C~70 °C, the range of sinter inlet temperature is 600 °C~750 °C, the range of sinter equivalent diameter is 6.05 mm~27.30 mm, the range of bed porosity is 0.35~0.6, the range of height of cooling section is 6 m~8.5 m, and the range of bed diameter is 9 m~11.5 m.
Figure 23 shows the error analysis of the heat transfer correlation formula and the calculation results. As can be seen from
Figure 23, the calculated results by the heat transfer correlation formula basically fall within the range of 30%, indicating that the calculation accuracy of this correlation formula is high, and the heat transfer correlation formula obtained by fitting is suitable for the prediction of the heat transfer effect in the furnace under this condition, which can provide an important reference for the subsequent research.