Numerical Simulation and Optimization of Perforated Tube Trolley in Circular Cooler

Show more

1. Introduction

In the iron and steel production process, the circular cooler is the most widely used equipment for the sintering process [1] . The sintering furnace cooling process directly determines the development of iron and steel enterprises. After decades of development and research, the equipment manufacturing process has become more and more matured [2] . FENG Junxiao has established two- dimensional heat transfer model of sintered porous media. After simulation, he got the best operating conditions of the parameters, the thickness of the material layer and the particle size [3] . The parameters such as the thickness and air volume of the sintering process were studied by Luo Yuanqiu. The results showed that the air volume and the thickness have an optimal value [4] . The gas-solid heat transfer coefficient of sintered ore in different working conditions was studied experimentally by ZHAO Bin. The porosity and the Nusselt number were sorted according to the experimental data. Through the ring cooling machine sintered ore channel gas-solid heat test, he has got the relationship between the product and the Reynolds number Curve [5] . The three-dimensional physical model of the two-energy equation of local non-thermodynamically balanced heat transfer is established by ZHANG Xiaohui. The numerical simulation of the two-phase heat transfer process and its influence parameters is carried out by using Fluent simulation software. And the orthogonal test method was used to combine the optimal parameters, and finally the conclusion of the experiment was used [6] .

However, the cooling process of the circular cooler trolley which layer by layer is still relatively simple, resulting in poor effect of the top area of the sinter and the longer cooling time. This paper starts from the internal structure of the trolley to make the design and optimization of the structure so that the cooling time of the sinter is greatly reduced and the recovery of waste heat is greatly improved.

2. Research Model

450 m^{2} sintering ring cooling machine as the research object, the size of the trolley 4 m, 2.5 m, 1.6 m. During the working process, the circular cooler trolley can be regarded as fixed, because the running speed of the tank is much lower than that of the air passing through the trolley. This process can be treated as a fixed bed of gas-solid heat exchange in the numerical simulation [7] . The built-in perforated tube trolley consists of a trolley wall, a grate plate and three perforated tubes. Below the grate plate, there are the main inlet cooling inlet, and the three sides of the three tubes are used as auxiliary cooling inlet. In order to facilitate the simulation calculation, the three-dimensional model of the trolley is simplified as a two-dimensional model. Built-in perforated tube trolley structure diagram and vertical sectional view as Figure 1 and Figure 2 shown.

Figure 1. Schematic diagram of the built-in perforated tube structure of the circular cooler.

Figure 2. Vertical sectional view of the ring trolley structure model.

3. Theoretical Basis and Basic Assumptions

In the process of cooling, there are heat conduction, convection and radiation. [8] Therefore, the heat exchange in the circular cooler can be solved by the equations of continuity equation, momentum conservation equation and energy conservation equation to describe.

3.1. Continuity Equation

According to the law of conservation of mass, the continuity equation can be obtained:

$\frac{\partial \rho}{\partial t}+\frac{\partial \left(\rho u\right)}{\partial x}+\frac{\partial \left(\rho v\right)}{\partial y}+\frac{\partial \left(\rho w\right)}{\partial z}=0$ (1)

In the above equation, $\rho $ is the density; $t$ is the time; $u$ , $v$ and $w$ is the component of the velocity vector, respectively x, y and z in the direction.

3.2. Momentum Conservation Equation

$\begin{array}{l}\frac{\partial \left(\rho u\right)}{\partial t}+\frac{\partial \left(\rho uu\right)}{\partial x}+\frac{\partial \left(\rho uv\right)}{\partial y}+\frac{\partial \left(\rho uw\right)}{\partial z}\\ =\frac{\partial}{\partial x}\left(\mu \frac{\partial u}{\partial x}\right)+\frac{\partial}{\partial y}\left(\mu \frac{\partial u}{\partial y}\right)+\frac{\partial}{\partial z}\left(\mu \frac{\partial u}{\partial z}\right)-\frac{\partial p}{\partial x}+{S}_{u}\end{array}$ (2a)

$\begin{array}{l}\frac{\partial \left(\rho v\right)}{\partial t}+\frac{\partial \left(\rho vu\right)}{\partial x}+\frac{\partial \left(\rho vv\right)}{\partial y}+\frac{\partial \left(\rho vw\right)}{\partial z}\\ =\frac{\partial}{\partial x}\left(\mu \frac{\partial v}{\partial x}\right)+\frac{\partial}{\partial y}\left(\mu \frac{\partial v}{\partial y}\right)+\frac{\partial}{\partial z}\left(\mu \frac{\partial v}{\partial z}\right)-\frac{\partial p}{\partial y}+{S}_{v}\end{array}$ (2b)

$\begin{array}{l}\frac{\partial \left(\rho w\right)}{\partial t}+\frac{\partial \left(\rho wu\right)}{\partial x}+\frac{\partial \left(\rho wv\right)}{\partial y}+\frac{\partial \left(\rho ww\right)}{\partial z}\\ =\frac{\partial}{\partial x}\left(\mu \frac{\partial w}{\partial x}\right)+\frac{\partial}{\partial y}\left(\mu \frac{\partial w}{\partial y}\right)+\frac{\partial}{\partial z}\left(\mu \frac{\partial w}{\partial z}\right)-\frac{\partial p}{\partial z}+{S}_{w}\end{array}$ (2c)

In the above equation, $\rho $ is the density; $t$ is the time; $u$ , $v$ and $w$ is the velocity vector in the x, y and z direction of the component; $p$ is the pressure on the fluid micro-body; $\mu $ is the dynamic viscosity; ${S}_{u}$ , ${S}_{v}$ and ${S}_{w}$ is the momentum conservation equation of the generalized source.

3.3. Energy Conservation Equation

Numerical simulations for built-in porous trolleys can be considered as porous media models. Without considering convection acceleration terms and diffusion terms, according to Darcy’s law [9] , the energy equation for porous media is as follows：

$\begin{array}{l}\frac{\partial}{\partial t}\left[\epsilon {\rho}_{f}{E}_{f}+\left(1-\epsilon \right){\rho}_{s}{E}_{s}\right]+\nabla \left[u\cdot \left({\rho}_{f}{E}_{f}+p\right)\right]\\ ={S}_{f}^{h}+\nabla \left[{k}_{eff}\nabla T-\left({\displaystyle \underset{i}{\sum}{h}_{i}{J}_{i}}\right)+\left(\stackrel{\xaf}{\stackrel{\xaf}{\tau}}\cdot u\right)\right]\end{array}$ (3)

${E}_{f}$ is the total energy of the fluid; ${\rho}_{f}$ is the density of the fluid; ${E}_{s}$ is the total energy of the solid medium; ${\rho}_{s}$ is the density of the solid medium; $\epsilon $ is the porosity; $T$ is the temperature; ${k}_{eff}$ is the effective thermal conductivity of the porous medium; ${S}_{f}^{h}$ is the source of the enthalpy of the fluid.

3.4. Basic Assumptions

Because the circulation, heat transfer and mass transfer process are very complicated [10] , in order to facilitate the numerical simulation of the sinter in the trolley, the following assumptions are put forward:

1) The sinter particles are idealized particles and are evenly distributed in the trolley, there is no temperature difference between the particles;

2) The wall in the trolley does not participate in heat transfer;

3) The operation parameters of the cooling trolley are not changed with time.

After making the assumptions, the traditional trolley is simulated. The results are compared with the relevant experimental papers, the results found that the error within 10%, so that the assumption is valid.

4. Optimization Design and Performance Analysis

4.1. Optimized Design

Put forward four factors: 1) Three pipe diameter d; 2) both sides of the tube from the grate plate height ${h}_{1}$ ; 3) the middle tube from the grate plate height ${h}_{2}$ ; 4) the distance between the two sides of the tube from the side $l$ . Those factors constitute a four factors and three levels orthogonal experiment. Combined with the trolley structure and perforated tube for the seamless steel pipe and other factors, perforated tube diameter ranges 200 mm from 400 mm, set d = 270 mm, 310 mm, 350 mm; taking the sinter stacking height and center layout principle into account, ${h}_{1}$ 250 - 500 mm and ${h}_{2}$ 550 - 950 mm for the best, take ${h}_{1}$ were 320 mm, 370 mm, 420 mm, ${h}_{2}$ were 640 mm, 740 mm, 840 mm; taking into account the sinter no clogging problems, $l$ should range 400 mm from 800 mm, take 500 mm, 600 mm, 700 mm. Structure parameters of calculation area shown in Figure 3.

The inlet air velocity is 1.1 m/s, the inlet temperature is 298 K, the initial temperature of sinter is 1073 K, and 373 K is set as cooling node. According to the principle of orthogonal test, select the nine groups of numerical combinations. Nine groups of perforated tube trolley were simulated by Fluent respectively. The value selection of the nine groups of study objects and the cooling time of the sinter from 1073 K to 373 K are shown in Table 1.

According to the Fluent simulation calculation result, the cooling time of the traditional trolley sinter from 1073 K to 373 K is 2130 s. It can be seen that the trolley built-in three porous tube reduced the average cooling time by 45% compared with the traditional one.

Figure 3. Structure parameters of calculation area.

Table 1. 9 sets of study object numerical parameters and cooling time.

4.2. Analysis of Waste Heat Recovery Utilization Amount

During the sintering production, great amount of sinter cooled air temperature ranges from 423 K to 623 K has high recovery value [11] . According to the calculation formula of waste heat recovery [12] :

$Q={c}_{p}m\Delta T=\left({c}_{p}\rho \Delta T\right)\cdot \left({q}_{v}\cdot t\right)$ (4)

In the formula,
$Q$ is the heat capacity of the waste heat recovery, (kJ);
${c}_{p}$ is the specific heat capacity of the air, 0.725 J/(kg∙˚C);
$\rho $ is the density of the air, 1.29 kg/m^{3};
$\Delta T$ is the temperature difference of the recovered value of the exhaust gas;
${q}_{v}$ is the volumetric flow rate of the inlet velocity, 4.4 m^{3}/s;
$t$ is waste heat recovery stage of the time, (s).

Table 1 in the nine sets of fixed value of the import wind speed numerical simulation calculation, get the corresponding waste heat utilization section as shown in Table 2.

The four factors and three levels orthogonal test is used to analyze the waste heat utilization as shown in Table 3.

The greater amount of waste heat recovery, indicating that the more heat can be used, and the relative energy consumption in the cooling process less. The greater the difference is, indicating the impact of this level on the results of the greater, the level has more important influence on the result. So the importance order of four parameters is l > d > h_{1} > h_{2}. The higher the K value is, the better level is in the factor. Select the largest set of data from K to get the preferred group d (3) h_{1} (3) h_{2} (3) l (1).

The simulation of the preferred group structure trolley, get up to 373 K when the time is 1155 s, the waste heat recovery and utilization time is 490 s, and the waste heat recovery is 403.28 kJ. In the nine groups of structures, the waste heat recovery of Group 9 is the highest. Group 9 and the Preferred group cooling time, waste heat recovery compared to Figure 4 and Figure 5.

Through the comparison between the Preferred group and the Group 9, the waste heat recovery ratio of the Group 9 structure trolley is much higher than that of the Preferred group, and the cooling time is very close.

In the cooling time, the fastest cooling rate of 373 K is the Group 8 trolley, but the Group 8 can use the waste heat recovery is only 57% to Group 9.

In the comprehensive comparison, the Group 9 is selected as the optimal structure combination.

Table 2. Perforated tube trolley and traditional trolley waste heat utilization.

Table 3. Analysis of orthogonal test results.

Figure 4. The cooling process of Group 9 and the preferred group.

Figure 5. The waste heat recovery of group 9 and the preferred group.

4.3. Comparison of Temperature Fields

In order to obtain the real-time effect of the built-in perforated tube cooling process, the temperature field analysis were carried out for the Group 9 and the traditional group structure trolley. The cloud field of the Group 9 structure trolley and traditional structure trolley at 800 s, 1200 s and 1600 s are shown in Figures 6-11 by Fluent numerical simulation.

In the traditional structure of the trolley, the cooling air can only enter the sintered area from beneath the grate plate. During the process of sintering the area, the air is continuously heated, so the high temperature sinter is cooled from bottom to top. In the built-in perforated tube structure trolley, the cooling air is equivalent to four outlets-the grate plate and the full side of the perforated tube 1, 2, 3. The upper of perforated tube part mainly by the perforated tube side of the cooling air to cool, the lower mainly from the grate plate into the cooling air to cool. Through the comparison of the temperature field of the Group 9 and the traditional trolley, it is found that the Group 9 trolley is more well-distri- buted than the traditional structure trolley in the sintering of the cooling.

5. Conclusions

1) Perforated tube structure trolley has greatly changed the temperature field of the cooling air in the sinter area, and accelerated the air flow velocity at the wall surface and the tapered area. It achieved the cooling effect better than that of the traditional structure trolley. It is of great significance to improve the cooling efficiency of the sinter and reduce the energy consumption during cooling.

2) The orthogonal test design method was used to construct the four factors and three levels test; the results show that the order of influence on the heat recovery is l > d > h_{1} > h_{2}. The optimal combination of three kinds of perforated tube trolley is the 9th group; the structural parameters are d = 350 mm, l = 500 mm, h_{1} = 420 mm, h_{2} = 740 mm; Compared with the traditional trolley, the best

Figure 6. Group 9 structure trolley 800 s air temperature distribution cloud.

Figure 7. Traditional structure trolley 800 s air temperature distribution cloud.

Figure 8. Group 9 structure trolley 1200 s air temperature distribution cloud.

Figure 9. Traditional structure trolley 1200 s air temperature distribution cloud.

Figure 10. Group 9 structure trolley 1600 s air temperature distribution cloud.

Figure 11. Traditional structure trolley 1600 s air temperature distribution cloud.

structure of the porous tube trolley has reduced 41% cooling time and increased 50% waste heat recovery

3) In this paper, when the model is numerically simulated, the sintering process is simplified; it only considered the convective heat transfer and ignored other forms of heat transfer, such as the radiation heat transfer of the sinter and the thermal conductivity of the trolley wall. Therefore, follow-up work proposed to take the radiation heat transfer model and wall heat transfer into consideration, to make the numerical simulation more consistent with reality.

Acknowledgements

The authors would like to acknowledge the Key Subject of Hunan Provincial Department of Education (17A183) support and encouragement to accomplish this work.

References

[1] Tang, X.J. (2008) China’s Sintering Industry, Technological Progress. Sintering Pellets, 33, 1-4.

[2] Pan, H. (2007) Sintering Circular Cooler Machine System Test Technology and Optimization Research. Liaoning University of Science and Technology.

[3] Feng, J.X., Liang, L.L., Zhang, M., et al. (2010 Mathematical Model of Thermal Process of Pellets in Ring Cooler. Reported, 32, 1596-1600.

[4] Dong, H., Li, J., Luo, Y.Q., et al. Experimental study on the Cooling Process of Sinter. Journal of Northeast University (Natural Science), No. 5, 689-692.

[5] Zhao, B., Zhao, L.J., Qu, T.T., et al. (2013) Relevant Correlation Field of Gas-Solid Heat Transfer Experiment in Sinter Channel of Annular Cooling Machine. Thermal Science and Technology, 12, 302-306.

[6] Zhang, X.H., Zhang, J.Y., Dai, C.D., et al. (2011) Numerical Simulation and Optimization of the Cooling Process of Sinter. Acta Chimica Sinica, 62, 3081-3087.

[7] Wang, F.J. (2013) Computational Fluid Dynamics Analysis. Tsinghua University Press, Beijing, 11.

[8] Shang, X.D. (2013) Optimization of Multiphase Flow of Blast Furnace Based on Computational Fluid Dynamics. Beijing University of Chemical Technology.

[9] Fluent User’s Guide. Fluent Inc. 14.5.

[10] Zhang, J.Y., Tian, W.Y., Dai, C.D., et al. (2012) Simulation and Optimization of Stratified Fabric of Circular Cooler. Acta Chimica Sinica, 63, 1385-1390.

[11] Zhao, H.G., Hu, S.Z. and Mao, Y.Z. (2009) Lean Iron Sintering Machine Exhaust Waste Heat Recovery Practice. Shandong Metallurgy, 31, 12-14.

[12] Zhu, M.K., Qing, D.F., Liu, S.X., et al. (2016) Research on Built-In Horizontal Pipe Structure and numerical Simulation of Annular Cooling Machine. Journal of South of China University (Natural Science Edition), 30, 22-27.