Study of Heat Transfer and the Hydrodynamic Performance of Gas–Solid Heat Transfer in a Vertical Sinter Cooling Bed Using the CFD-Taguchi-Grey Relational Analysis Method

: The vertical sinter cooling bed (VSCB) is a high-e ﬃ ciency energy-saving and environmentally friendly waste heat recovery equipment. In this study, a computational ﬂuid dynamics (CFD) convection model was established to reveal the typical factors on the thermodynamic performance in VSCB. Indeed, a multiple performance optimal algorithm based on the Taguchi-grey relational analysis (GRA) method was ﬁrst applied to investigate the e ﬀ ects of geometric and operational factors, including the diameter of the bed, height of the bed, air mass ﬂow rate, air inlet temperature, and sinter mass ﬂow rate, on improving the heat transfer ( Ex ) and hydrodynamic performance ( P drop ) and obtain the optimum combination of each factor in VSCB. The results found that the diameter of the bed was the most inﬂuential factor contributing the multiple types of performance with a contribution rate of 70.51%, followed by the air mass ﬂow rate (15.84%), while the height of the bed (0.27%) exerted a limited e ﬀ ect on the performance of multiple processes. The optimal combination of factors (A 1 B 5 C 5 D 5 E 1 ) was compared with the initially selected parameters by performing a conﬁrmation test. The performances of heat transfer and hydrodynamics were improved by the Taguchi with the GRA method. sinter ﬂux unit quantity of sinter granules mass ﬂow of of the sinter of in velocity of of thermodynamics, in the mass ﬂow rate in the air outlet temperature; of sinter outlet temperature in the air outlet temperature. pressure drop airﬂow sinter air mass ﬂow drop of airﬂow 5-fold the air mass ﬂow rate to 160 kg · s − 1 superﬁcial velocity of air an increase in air mass ﬂow rate correlated with the inertial and viscous terms within the pressure drop of airﬂow through the VSCB. in superﬁcial velocity of in the Abstract: The vertical sinter cooling bed (VSCB) is a high-efficiency energy-saving and environmentally friendly waste heat recovery equipment. In this study, a computational fluid dynamics (CFD) convection model was established to reveal the typical factors on the thermodynamic performance in VSCB. Indeed, a multiple performance optimal algorithm based The and the , and exergy Ex ( T the beginning of the process the maximum at 80 kg · s 1 , and with an increased air mass ﬂow rate. The maximum values of Ex ( T ) and Ex were 26.61 MJ · s − 1 and 26.36 MJ · s − 1 , respectively. The di ﬀ erence between Ex ( T ) and Ex increased with an increase in the air mass ﬂow rate, while the pressure Ex ( p ) decreased with an increase in the air mass ﬂow rate. The increase in the air mass ﬂow rate resulted in an increase in the superﬁcial velocity of the air and a loss of airﬂow resistance. Therefore, the pressure exergy decreased with an increase in the air mass ﬂow rate and subsequently increased the di ﬀ erence between Ex and Ex ( T ).


Introduction
The vertical cooling bed is widely used in various industries, such as the metallurgical [1], energy [2], chemical [3], and food [4] sectors. In particular, the steel industry accounts for 23.5% of the energy consumption by the entire industrial sector, which is deemed a process with extensive energy consumption and emission production [5]. Several types of vertical cooling packed beds have been employed in different processes during steel production, such as a slag furnace, vertical pellet furnace, coke dry quench, and sinter vertical packed bed. According to the published literature [6], the energy consumed by the sintering process accounts for 15% of the total energy consumption of the steel industry. Therefore, studies designed to improve the efficient utilization and recovery of waste heat resources from the sintering process are very important. Furthermore, researchers have focused on developing several methods to reduce energy consumption and increase the efficiency of utilization of available recovered energy waste [6][7][8]. In general, the circular and linear moving sinter cooling beds use gas leakage with low efficiency. The novelty vertical sinter cooling packed bed was invented for the efficient recovery and utilization of waste heat from the sinter cooling process. Although previous studies have focused on the design parameters, the most effective factor and the contribution rate of the factor to the thermodynamic characteristics have not been investigated using a systematic vertical sinter cooling bed was the objective function, and the effects of five design factors were investigated with the help of a three-dimensional computational fluid dynamics (CFD)-simulated model. In addition, The Taguchi-grey relational analysis method was used to determine the optimum combination of design factors, simultaneously impacting multiple target functions. Therefore, the aim of this study was to determine the optimal combination of multiple performance characteristics for simultaneously increasing the exergy of the air outlet and reducing the pressure drop in the vertical bed.

Vertical Sinter Cooling Bed (VSCB)
This technical innovation creates a variation in the heat transfer mechanism that transfers crossflow heat transmission in a circular cooler to counter contact current heat transfer in the VSCB. As shown in Figure 1, the sinter cooling process is operating in the cooling section, where heat transfer is initiated between hot sinter particles and ambient temperature air with direct counter-current contact modes. The blazing sinter block is initially filled in the pre-store chamber that guarantees continuous sinter granular charging and then packed into the cooling section; the air is simultaneously introduced from the center and circular channels and flows through center air cowl into the cooling section. The cooled sinter granule is unloaded using a discharge device at the bottom of the vertical bed and heated air flows through a ramp channel and is recovered at the outlet of the VSCB. The discharge device is composed of an electromagnetic vibrator and rotary sealing valve, ensuring the continuous and smooth discharge of sinter granules at a constant mass flow rate. As described above, the heat transfer process between sinter granular and air occurs in a confined space (cooling section), which contributes to improving the heat transfer efficiency and reducing air leakage.

Studies of the Heat Transfer Performance of the VSCB
In the past decade, many investigations aiming to improve heat transfer performance in terms of the exergy and energy have been performed. In general, theoretical, experimental, and simulated methods have been applied to improve the heat transfer performance of the VSCB. The published studies have mainly examined the effect of geometric and operating parameters on the heat transfer process in the VSCB. Feng et al. [9] presented the effects of the air inlet temperature and mass flow rate on the recovered exergy performance based on a three-dimensional simulated model. Zheng et al. [10,11] presented the effects of several design parameters on the volumetric exergy transfer

Studies of the Heat Transfer Performance of the VSCB
In the past decade, many investigations aiming to improve heat transfer performance in terms of the exergy and energy have been performed. In general, theoretical, experimental, and simulated methods have been applied to improve the heat transfer performance of the VSCB. The published studies have mainly examined the effect of geometric and operating parameters on the heat transfer process in the VSCB. Feng et al. [9] presented the effects of the air inlet temperature and mass flow rate on the recovered exergy performance based on a three-dimensional simulated model. Zheng et al. [10,11] presented the effects of several design parameters on the volumetric exergy transfer coefficient and volumetric transfer coefficient in terms of the second law of thermodynamics and pilot experimental method. Pan et al. [12] established a theoretical model of gas-solid heat exchange to determine the effects of the operational parameters on the heat transfer and pressure drop. Feng et al. [13] determined the overall heat transfer coefficient using a pilot experimental method. In addition, experimental research has been conducted, and the reliability of the Nusselt number has been obtained from regression analysis and verification. Zhang et al. [14] investigated the effects of structural and operating parameters on the exergy performance of waste heat recovery using a 3D calculated model. Four main parameters were studied in terms of exergy and energy consumption. Feng et al. [15] described the sinter bed system used to obtain the optimal combination of design parameters with the help of an orthogonal experiment to improve the heat transfer performance based on a 3D porous medium and non-equilibrium heat exchange model. Pan et al. [16] studied the volumetric heat transfer coefficient and air pressure drop of a sinter layer based on the experimental method. The research on the vertical sinter cooling bed also examined the hydrodynamic performance, such as the pressure drop and energy loss, which are deemed as functions with a negative effect on the effectiveness of the heat resource recovery process. Feng et al. [17] investigated the airflow pattern and pressure drop characteristics using an experimental approach. The flow regimes in the sinter bed were also partitioned into several parts in terms of the critical Reynolds number. The relationship between the pressure drop and the parameters of interest were established simultaneously. Tian et al. [18] studied the air pressure drop through a random distribution sinter packed bed by varying the irregular particle sizes, the porosity of the sinter layer, and sphericity for different flow regimes.
The aforementioned studies, examining the heat transfer performance in a vertical sinter cooling bed, have mainly determined the effects of structural and operational parameters on the energy and exergy performance, which is beneficial for providing comprehensive insights into the heat transfer characteristics in a vertical sinter bed. However, few studies have employed a systematic analytical approach with a statistical analysis of factors of interest that affect the heat transfer performance of the vertical sinter cooling bed along with a simultaneous assessment of multiple purpose functions. To our knowledge, the Taguchi with grey relational analysis (GRA) was established to analyze the effective parameters that affect multiple purpose functions and is necessary to obtain the relationship of the simulated and ideal results in the different design combinations.

Studies of the Statistical Analysis Method
The Taguchi grey relational analysis is a versatile technique that has been applied in various fields, such as iron and steel [19,20], heat exchanger [21], energy systems [22], refrigeration [23], powder metallurgy [24], building material [25], and chemical mechanical planarization [26]. In the present study, the main focus was the application of the Taguchi grey relational analysis method to determine the effects of the design parameters on the heat transfer recovery process. In the literature, many researchers have paid their attention to this topic. N. Celik et al. [27] simultaneously examined the effects of design parameters on the performance of a tabulated heat exchanger and the effects of the optimum design parameters on heat and pressure loss. Kumbhar et al. [28] used the CFD method to investigate the performance of a simple tube with twisted tapes. The twisted tape pitch and dimple pitch are the most important parameters of the system, and the best parameter combination in terms of the performance of heat transfer and friction was identified using the Taguchi-grey relational analysis. Acir et al. [29] investigated the optimal parameter, impacting a novel solar air heater system for efficient energy and exergy transfer using the Taguchi grey relational analysis method. Kuo et al. [30] determined the effects of selected factors on the performance of a flat-plate collector. The optimum operating factors for this system were determined using the Taguchi grey relational analysis. Gunes et al. [31] determined the effects of four optimal design parameters on the two response functions (heat transfer and pressure drop) using the Taguchi grey relational analysis.
The factors and selected range of the vertical sinter bed used in previous studies are listed in Table 1. The related factors and levels are presented in Table 1 as a reference for selecting the factors of interest. Feng et al. [32] Energy and exergy analysis 293-373 5-9 160-200 6-10 -Based on previous research, the present study mainly focused on the parameters of interest by determining the order of importance of effective factors and contribution rates of factors for simple or multiple target performance in a heat transfer process using a VSCB. The Taguchi grey relational analysis method was introduced to simultaneously improve heat transfer and hydrodynamic performance. First, the CFD convection model was established to determine the variations in heat transfer and hydrodynamic performance characteristics for various factors with different ranges. Then, the Taguchi L 25 (5 5 ) experimental array comprising five concerned factors, with five levels for each factor, was separately applied to investigate the effect of the factor on simple response performance (Ex or P drop ) using the Taguchi method. Furthermore, the Taguchi-grey relational analysis method was used to determine the order of effects of the factors of interest on the heat transfer and hydrodynamic performance (Ex and P drop ) and to simultaneously improve the performance of both responses.

Mathematical Description within the CFD
The airflow through the void space between the sinter ore and sinter particles potentially functions as a sieve, allowing the air to penetrate. The RANS (Reynolds Averaged Navier Stokes) k-ε turbulence model is used to describe the conservation of air hydrodynamic Continuity: Momentum equation: where ρ is the air density; u is air the velocity vector; p is the air pressure; I is the unit direction vector; µ is the dynamic viscosity; µ T is the eddy viscosity, which is defined as a function of the turbulent kinetic energy (k) and turbulent energy dissipation rate (ε); F is a volumetric force vector; S is the source term, which includes the viscosity loss and inertial loss term and controlled by the porous media of Navier-Stokes macro-equations. The source term is defined as follows: where 1/α is coefficient of viscous resistance; C 2 is the coefficient of inertial resistance, and u is the air velocity. In this study, the coefficients of viscous and inertial resistance were used to determine the air pressure drop through vertical sinter packed bed based on the modified Ergun correlation [15]. Then, 1/α and C 2 are calculated as follows: where ϕ is the porosity of vertical sinter cooling bed, and d pSM is the Sauter mean diameter (SMD) applied to describe the mean diameter of multi-dispersed sinter particles. The Sauter mean diameter (SMD) is calculated as follows: where x i is the size fraction of granular in the size interval i . In this study, the average diameter was 0.035 m. The turbulent transportation equation (RANS k-ε) is comprised of turbulent energy transportation and turbulent dissipation transportation equations. The RANS k-ε model determines the behavior of airflow through the vertical sinter packed bed in the turbulent flow range (Re > 300), which is calculated as follows: where Cµ = 0.09, Cε1 = 1.44, and Cε2 = 1.92 are turbulent constants; and σ k = 1.00 and σ ε = 1.3 are Prandtl numbers. The temperature differs between the air and sinter particles in the heat transfer process; thus, in this study, the local thermal non-equilibrium (LTNE) model was applied to describe the behavior of air-sinter heat transfer. The following conservation equations of air and sinter energy are used: Energy: Fluid (air): Solid (sinter): where ρ f and ρ s are the density of air and sinter, respectively; c f and c s are the specific heat of air and sinter, respectively; u f and u s are the velocity of the air and sinter, respectively; λ f and λ s are the thermal conductivity of air and sinter, respectively; T f and T s are the temperature of air and sinter, respectively; h v is volume coefficient of air-sinter heat transfer that connects air and sinter heat transfer process. The energy source term is added to the last position of the fluid and solid conservation energy equation, which controls the convective heat transfer rate between air and sinter. The Achenbach criterion relation is used to express the volume coefficient of heat transfer h v as follows: where h is the area coefficient of convective heat transfer, according to the results reported by Feng et al. [15], who described air and sinter heat transfer behaviors through a vertical sinter cooling packed bed. The modified correlation with the Nusselt number is determined as the coefficient of convective heat transfer h based on the area of heat transfer between air and sinter. The modified correlation, defined as the Nusselt number, is calculated as follows: where Re and Pr are the sinter Reynolds number and air Prandtl number, respectively. Re and Pr are calculated as follows: In this study, one of the target physical phenomena was the exergy of outlet air, which was significantly affected by the geometry and operating parameters as a response to the variations in several considered factors. Exergy represents the maximum amount of useful work that can be obtained from a vertical sinter packed bed and is deemed as a progressive tool for eliminating energy consumption. Thus, we chose it as a target value in the present study. The exergy of outlet air is calculated as follows: where m is the mass flow rate of air, kg·s −1 ; p 0 is the ambient pressure; ∆p is the pressure drop when air flows through the sinter bed, ∆p = p inlet − p outlet , p inlet , and p outlet are air pressure at the initial inlet and final outlet, respectively. T 0 is the ambient temperature and defined as 298 K. The physical properties of the sinter granular material and air in the calculation are listed in Table 2. Sinter ore density [14] kg·m −3 2800 Notably, the equation used to calculate the exergy of outlet air is comprised of two terms. The first term on the right is the temperature exergy that provides the impetus for obtaining useful work in the air and sinter heat transfer process. The second term on the right is pressure exergy related to the resistance of the vertical sinter bed. The pressure exergy describes the air pressure drop through the sinter bed related to energy consumption, and its role in eliminating useful work is obtained. In this regard, pressure exergy undermines the heat transfer process between air and sinter. The temperature exergy and pressure exergy are calculated as follows: The various assumptions are presented based on the criteria of a homogenous distribution and statistical analysis to ensure simple and smooth calculations. The cooling section is deemed as the system comprised of a homogenous distribution of sinter granular material and continuous airflow in the heat transfer process. Therefore, several assumptions related to the cooling section for the heat transfer calculation in VSCB are listed below.

1.
The sinter granules maintain the mass flow pattern in the cooling section, which is the precondition to ensure steady operation.

2.
The sinter granules are deemed a homogenous and isotropic bulk solid material.

3.
The heat dissipates at the outer wall of the VSCB is insulated by the insulation material.

Determination of the Grid Generation and Boundary Condition
The ICEM-CFD (Integrated Computer Engineering and Manufacturing code for Computational Fluid Dynamics) tool is usually applied to generate a mesh in a physical model with complex geometry. According to a previous study, the target cooling section is defined as a calculated domain [14,15,32]. Smaller calculated workloads, a more rapid convergence process, and more boundary conditions are adopted to obtain a higher mesh quality. Thus, the 3D structured mesh was established to calculate the study target, which was comprised of a quadrilateral and hexahedral grid. The hexahedral grid was applied to generate the fluid domain mesh, and the quadrilateral grid was employed to create the mesh of the faces. The ICEM was applied to create the surface and volume grid and exported to the finite element software ANSYS Fluent 16.0 for calculation. The pressure-velocity coupling method was applied to the SIMPLE algorithm to discretize the control element, and a second upstream scheme was used for pressure and energy spatial discretization. The air and sinter heat transfer cooling system is very complex. Several suitable UDF (user-defined functions) are available that utilize the C programming language to define the parameters, such as the properties of air and sinter, volume heat transmission coefficient (h v ) in energy conservation equations, flux and diffusivity coefficients, and source term (S i ), in user-defined scalar (UDS) transportation equations.
The boundary condition settings are described in Figure 2. The air inlet section is defined as the mass flow inlet boundary condition (the circular flow channel at the bottom of the bed, and center flow channel at the central air cowl) to describe the mass flow rate in compressible airflow. The air outlet section is defined as the pressure outlet boundary condition that defines the outlet with static pressure and is beneficial for better convergence during the simulation in which backflow occurred. The designed specific values for turbulent intensity and hydraulic diameter were used in the specification method. In addition, the relative pressure value at the air outlet was set to zero because the effective pressure drop between the air and wall is not considered for the resistance loss in the computation domain. The other parts of this domain were set as the adiabatic wall condition boundary.  Sinter ore density [14] kg·m −3 2800

Determination of the Grid Generation and Boundary Condition
The ICEM-CFD (Integrated Computer Engineering and Manufacturing code for Computational Fluid Dynamics) tool is usually applied to generate a mesh in a physical model with complex geometry. According to a previous study, the target cooling section is defined as a calculated domain [14,15,32]. Smaller calculated workloads, a more rapid convergence process, and more boundary conditions are adopted to obtain a higher mesh quality. Thus, the 3D structured mesh was established to calculate the study target, which was comprised of a quadrilateral and hexahedral grid. The hexahedral grid was applied to generate the fluid domain mesh, and the quadrilateral grid was employed to create the mesh of the faces. The ICEM was applied to create the surface and volume grid and exported to the finite element software ANSYS Fluent 16.0 for calculation. The pressurevelocity coupling method was applied to the SIMPLE algorithm to discretize the control element, and a second upstream scheme was used for pressure and energy spatial discretization. The air and sinter heat transfer cooling system is very complex. Several suitable UDF (user-defined functions) are available that utilize the C programming language to define the parameters, such as the properties of air and sinter, volume heat transmission coefficient (hv) in energy conservation equations, flux and diffusivity coefficients, and source term (Si), in user-defined scalar (UDS) transportation equations.
The boundary condition settings are described in Figure 2. The air inlet section is defined as the mass flow inlet boundary condition (the circular flow channel at the bottom of the bed, and center flow channel at the central air cowl) to describe the mass flow rate in compressible airflow. The air outlet section is defined as the pressure outlet boundary condition that defines the outlet with static pressure and is beneficial for better convergence during the simulation in which backflow occurred. The designed specific values for turbulent intensity and hydraulic diameter were used in the specification method. In addition, the relative pressure value at the air outlet was set to zero because the effective pressure drop between the air and wall is not considered for the resistance loss in the computation domain. The other parts of this domain were set as the adiabatic wall condition boundary.

Grid Independence Test and Model Validation
The computational domain is presented for this study in Figure 2. The function of the mesh independence test was to examine the effect of the number of elements on the simulation data and to obtain the applicable grid for subsequent investigations. The structured grids with quadrilateral elements were employed to display the solid and fluid zones for the non-equilibrium thermal model, where the inlet and outlet boundary conditions were increased for valid computation. The criteria used to judge a high quality of mesh were sufficiently fine to present an accurate and detailed description of flow modes and sufficiently coarse to simulate the process in a suitable time. Therefore, five mesh elements ranging from 19,232 to 399,722 were employed to simulate the heat transfer process. The simulated results are described in Table 3, which presented a gradual decrease in the deviation of previous and subsequent mesh elements overall. The air and sinter outlet temperatures and pressure drop were set as target functions to determine the reliable mesh element. As shown in Table 3, the relative deviation of the results between case 4 and case 5 was −0.003, −0.018, and −0.001, respectively. As described above, the case 4 (172,872) was the suitable mesh element used in subsequent calculations. A vertical sinter cooling device that includes heat transfer and pressure drop characteristics was established in this study to verify the reliability of the results from the simulation. The dimensions and components of the experimental device are presented in Figure 3. The calculated diameter and height of the sinter cooling bed were 0.6 m and 1.6 m, respectively. A U-tube water manometer was applied to measure the pressure drop of airflow through the sinter layer. The pressure measurement ranged from 0 Pa to 3000 Pa. Six thermocouples were used to measure the variation in air temperature and were located from the bottom to the top of the sinter layer bed at intervals of 200 mm. The outer wall of the device was wrapped with insulated Rockwool to reduce the amount of heat transmitted to the external environment. The air was continuously introduced by a centrifugal fan located at the bottom of the sinter bed, while a throttle valve was used to regulate the airflow rate through the vertical sinter bed. The heated sinter granules randomly fell onto the sinter bed until the appropriate height was attained. For the initial heat transfer between air and sinter granular material, the sinter granular material flew out with the assistance of a discharge device with a constant mass flow rate. The air temperature and pressure drop were measured using a pilot design of an experimental sinter vertical bed (Figure 3), and the results are shown in Table 4. The air outlet temperature and air pressure drop were selected as target validation objects to compare the experimental data and results of the simulation. Five operating cases were used to compare the results from the simulated and measured devices, as illustrated in Table 4. The mean deviation of the air outlet temperature and pressure drop was 4.25% and 3.53%, respectively, indicating that the calculated values were similar to the experimental measurements. For a further analysis of the reliability of data obtained from the simulation, the dataset obtained from the vertical sinter beds in published reports is introduced in this section. The comparison of measured and calculated air outlet temperatures of the investigated vertical sinter bed is shown in Figure 4. The specific data obtained from the measured and simulated air outlet temperature are described in Table 5. The deviation of all cases used for validation was less than 5% and was consistent with the reported data. Therefore, the reliability of the simulated model analyzed in this study was verified and used to accurately describe the heat transfer performance of the vertical sinter cooling bed. The air temperature and pressure drop were measured using a pilot design of an experimental sinter vertical bed (Figure 3), and the results are shown in Table 4. The air outlet temperature and air pressure drop were selected as target validation objects to compare the experimental data and results of the simulation. Five operating cases were used to compare the results from the simulated and measured devices, as illustrated in Table 4. The mean deviation of the air outlet temperature and pressure drop was 4.25% and 3.53%, respectively, indicating that the calculated values were similar to the experimental measurements. Table 4. Comparison of the results from the simulation and experiment. For a further analysis of the reliability of data obtained from the simulation, the dataset obtained from the vertical sinter beds in published reports is introduced in this section. The comparison of measured and calculated air outlet temperatures of the investigated vertical sinter bed is shown in Figure 4. The specific data obtained from the measured and simulated air outlet temperature are described in Table 5. The deviation of all cases used for validation was less than 5% and was consistent with the reported data. Therefore, the reliability of the simulated model analyzed in this study was verified and used to accurately describe the heat transfer performance of the vertical sinter cooling bed.

Results and Discussion
Many factors affect the fluid flow and thermodynamic performance of a sinter bed for recovering heat wasted during sinter cooling. The simulation studied the geometric and operational parameters for the heat transfer process in a VSCB and compared them to the thermodynamic characteristics. According to the results from previous studies and actual production of a VSCB for a sinter cooling procedure, the operational and geometric factors play important roles in exergy transfer and fluid flow in a vertical sinter cooling process. The operational parameters, namely, the air mass flow rate, the cooling air inlet temperature, and sinter inlet mass flow rate, and the structural parameters, namely, the height and diameter of the sinter bed, are selected as important factors that significantly alter the fluid flow and heat transfer behaviors. Therefore, the investigation of the effect of geometric and operating parameters on the sinter cooling procedure in the VSCB is necessary.

Effect of the Mass Flow Rate of Air
The mass flow rate of air is an important parameter in the vertical sinter cooling process. After setting the five parameters to specific values, the variations of the outlet temperature and pressure drop in the cooling bed were obtained for airflow mass rates of 60, 80, 100, 120, and 140 kg·s −1 , as shown in Figure 5. In this experiment, the sinter mass flow rate (100 kg·s −1 ), the effective height of the sinter bed (6 m), the diameter of the sinter bed (10 m), air inlet temperature (293 K), and sinter inlet temperature (973 K) were fixed. As shown in Figure 5, the air and sinter outlet temperatures peaked at 968 K and 512 K, respectively, when the air mass flow rate was 60 kg·s −1 . Moreover, the air outlet

Results and Discussion
Many factors affect the fluid flow and thermodynamic performance of a sinter bed for recovering heat wasted during sinter cooling. The simulation studied the geometric and operational parameters for the heat transfer process in a VSCB and compared them to the thermodynamic characteristics. According to the results from previous studies and actual production of a VSCB for a sinter cooling procedure, the operational and geometric factors play important roles in exergy transfer and fluid flow in a vertical sinter cooling process. The operational parameters, namely, the air mass flow rate, the cooling air inlet temperature, and sinter inlet mass flow rate, and the structural parameters, namely, the height and diameter of the sinter bed, are selected as important factors that significantly alter the fluid flow and heat transfer behaviors. Therefore, the investigation of the effect of geometric and operating parameters on the sinter cooling procedure in the VSCB is necessary.

Effect of the Mass Flow Rate of Air
The mass flow rate of air is an important parameter in the vertical sinter cooling process. After setting the five parameters to specific values, the variations of the outlet temperature and pressure drop in the cooling bed were obtained for airflow mass rates of 60, 80, 100, 120, and 140 kg·s −1 , as shown in Figure 5. In this experiment, the sinter mass flow rate (100 kg·s −1 ), the effective height of the sinter bed (6 m), the diameter of the sinter bed (10 m), air inlet temperature (293 K), and sinter inlet temperature (973 K) were fixed. As shown in Figure 5, the air and sinter outlet temperatures peaked at 968 K and 512 K, respectively, when the air mass flow rate was 60 kg·s −1 . Moreover, the air outlet temperature was 935 K for an air mass flow rate of 80 kg·s −1 , which exhibited a substantial decrease and an accelerated rate over the next three points to reach a temperature slightly higher than 648 K at 160 kg·s −1 . However, the decreasing trend of the sinter outlet temperature was gradually decelerated. A potential explanation for this finding is that the increase in the air mass flow rate leads to an increase in the superficial velocity of air, and the volumetric heat transfer coefficient increases with an increase in air superficial velocity. When the sinter inlet heat flux was constant, the unit quantity of heat transfer between air and sinter granules decreased as the air mass flow rate increased, resulting in the deceleration of the decreasing trend of the sinter outlet temperature. When the cross-sectional area of VSCB was constant, the air mass flow rate increased with an increase in the superficial velocity of air. In accordance with the first law of thermodynamics, an increase in the mass flow rate led to a decrease in the air outlet temperature; however, the decreasing trend of the sinter outlet temperature decelerated and led to a substantial decrease in the air outlet temperature. Furthermore, the pressure drop of airflow through the sinter bed layer increased with an increasing air mass flow rate, as shown in Figure 5. The pressure drop of airflow increased by approximately 5-fold when the air mass flow rate increased from 60 to 160 kg·s −1 . The superficial velocity of air increased with an increase in air mass flow rate and positively correlated with the inertial and viscous terms within the pressure drop of airflow through the VSCB. Therefore, the increase in the superficial velocity of air led to an increase in the loss of airflow resistance. The air mass flow rate also affects the heat transfer and hydrodynamic performances, and thus the variations in the enthalpy exergy Ex, temperature exergy Ex(T), and pressure exergy Ex(p) were measured for various air mass flow rates, as shown in Figure 6. The enthalpy exergy Ex and temperature exergy Ex(T) simultaneously increased at the beginning of the process until they reached the maximum value at 80 kg·s −1 , and then decreased with an increased air mass flow rate. The maximum values of Ex(T) and Ex were 26.61 MJ·s −1 and 26.36 MJ·s −1 , respectively. The difference between Ex(T) and Ex increased with an increase in the air mass flow rate, while the pressure Ex(p) decreased with an increase in the air mass flow rate. The increase in the air mass flow rate resulted in an increase in the superficial velocity of the air and a loss of airflow resistance. Therefore, the pressure exergy decreased with an increase in the air mass flow rate and subsequently increased the difference between Ex and Ex(T).

Effect of the Air Inlet Temperature
The variations in the air and sinter outlet temperatures and pressure drop of airflow in the sinter bed with changes in the air inlet temperature are depicted in Figure 7. Under this condition, a fixed sinter mass flow rate (100 kg/s), the effective height of the sinter bed (6 m), the diameter of the sinter bed (10 m), air mass flow rate (140 kg/s), and sinter inlet temperature (973 K) were used. The air and sinter outlet temperatures and pressure drop increased as the air inlet temperature increased because the increase in air inlet temperature increased the superficial velocity of air. However, a decreased cooling time between the gas and solid led to the inefficient transfer of heat between air and sinter, resulting in a linear increase in the sinter outlet temperature. The thermal energy transported in the sinter bed from the air inlet increased much faster than the energy transported out of bed by the sinter outlet, resulting in an increase in the air outlet temperature. Furthermore, when the increasing trend of the sinter outlet temperature presented linear growth, the increasing trend of air outlet temperature was gradually decelerated. In accordance with the first law of thermodynamics, the air and sinter outlet temperatures were increased by 5.3 K and 17 K, respectively, if the air inlet temperature increased by 20 K.

Effect of the Air Inlet Temperature
The variations in the air and sinter outlet temperatures and pressure drop of airflow in the sinter bed with changes in the air inlet temperature are depicted in Figure 7. Under this condition, a fixed sinter mass flow rate (100 kg/s), the effective height of the sinter bed (6 m), the diameter of the sinter bed (10 m), air mass flow rate (140 kg/s), and sinter inlet temperature (973 K) were used. The air and sinter outlet temperatures and pressure drop increased as the air inlet temperature increased because the increase in air inlet temperature increased the superficial velocity of air. However, a decreased cooling time between the gas and solid led to the inefficient transfer of heat between air and sinter, resulting in a linear increase in the sinter outlet temperature. The thermal energy transported in the sinter bed from the air inlet increased much faster than the energy transported out of bed by the sinter outlet, resulting in an increase in the air outlet temperature. Furthermore, when the increasing trend of the sinter outlet temperature presented linear growth, the increasing trend of air outlet temperature was gradually decelerated. In accordance with the first law of thermodynamics, the air and sinter outlet temperatures were increased by 5.3 K and 17 K, respectively, if the air inlet temperature increased by 20 K.

Effect of the Air Inlet Temperature
The variations in the air and sinter outlet temperatures and pressure drop of airflow in the sinter bed with changes in the air inlet temperature are depicted in Figure 7. Under this condition, a fixed sinter mass flow rate (100 kg/s), the effective height of the sinter bed (6 m), the diameter of the sinter bed (10 m), air mass flow rate (140 kg/s), and sinter inlet temperature (973 K) were used. The air and sinter outlet temperatures and pressure drop increased as the air inlet temperature increased because the increase in air inlet temperature increased the superficial velocity of air. However, a decreased cooling time between the gas and solid led to the inefficient transfer of heat between air and sinter, resulting in a linear increase in the sinter outlet temperature. The thermal energy transported in the sinter bed from the air inlet increased much faster than the energy transported out of bed by the sinter outlet, resulting in an increase in the air outlet temperature. Furthermore, when the increasing trend of the sinter outlet temperature presented linear growth, the increasing trend of air outlet temperature was gradually decelerated. In accordance with the first law of thermodynamics, the air and sinter outlet temperatures were increased by 5.3 K and 17 K, respectively, if the air inlet temperature increased by 20 K.  The variations in the enthalpy exergy Ex, temperature exergy Ex(T), and pressure exergy Ex(p) were calculated for various air inlet temperatures, as shown in Figure 8. The enthalpy exergy Ex and temperature exergy Ex(T) monotonically decreased with an increase in the air inlet temperature. The potential explanation for this result is that the increase in the air inlet temperature results in an increase in the air outlet temperature and, subsequently, an increase in air energy levels. As the air inlet temperature increased, the temperature difference between air and sinter granules decreased, leading to a decrease in the quality of the heat transfer process. Furthermore, as the air inlet temperature increased, the increase in the amplitude of the air energy level was smaller than gas-solid heat transfer recovery, leading to decreased Ex. For every 20 K increase in the air inlet temperature, the enthalpy exergy Ex decreased by a mean of 2.37 MJ·s −1 . The pressure exergy Ex(p) gradually decreased with an increase in air inlet temperature because the increase in the air inlet temperature increased the air volume flow rate and resulted in a larger gas resistance loss, leading to a gradual decrease in pressure exergy Ex(p). The variations in the enthalpy exergy Ex, temperature exergy Ex(T), and pressure exergy Ex(p) were calculated for various air inlet temperatures, as shown in Figure 8. The enthalpy exergy Ex and temperature exergy Ex(T) monotonically decreased with an increase in the air inlet temperature. The potential explanation for this result is that the increase in the air inlet temperature results in an increase in the air outlet temperature and, subsequently, an increase in air energy levels. As the air inlet temperature increased, the temperature difference between air and sinter granules decreased, leading to a decrease in the quality of the heat transfer process. Furthermore, as the air inlet temperature increased, the increase in the amplitude of the air energy level was smaller than gassolid heat transfer recovery, leading to decreased Ex. For every 20 K increase in the air inlet temperature, the enthalpy exergy Ex decreased by a mean of 2.37 MJ·s −1 . The pressure exergy Ex(p) gradually decreased with an increase in air inlet temperature because the increase in the air inlet temperature increased the air volume flow rate and resulted in a larger gas resistance loss, leading to a gradual decrease in pressure exergy Ex(p).

Effect of the Sinter Mass Flow Rate
In this section, the effect of the sinter mass flow rate on heat transfer and gas fluid performances was investigated. The variations in the air and sinter outlet temperatures and pressure drop with varying sinter import mass flow rate are depicted in Figure 9. The sinter import mass flow was set to 300 t/h to imitate the standard sinter production yield, whereas the other import mass flow rate was used to imitate the uncertainty and fluctuations during sinter production. The behavior of air and sinter outlet temperatures and pressure drop in the bed layer increased with an increase in the sinter import mass flow rate, as shown in Figure 9. The outlet temperature and pressure drop increased linearly with the increase in the sinter import mass flow rate. For each 5 t/h increase in the sinter import mass flow rate, the air and sinter outlet temperatures increased by a mean of 5.43 K and 4.08 K, respectively, and the loss of air resistance increased by a mean of 21 Pa. As the sinter import mass flow rate increased, a greater amount of thermal energy was carried in the sinter bed and increased the air and sinter temperatures, leading to increases in the superficial velocity of air in the sinter bed and the loss of air resistance.

Effect of the Sinter Mass Flow Rate
In this section, the effect of the sinter mass flow rate on heat transfer and gas fluid performances was investigated. The variations in the air and sinter outlet temperatures and pressure drop with varying sinter import mass flow rate are depicted in Figure 9. The sinter import mass flow was set to 300 t/h to imitate the standard sinter production yield, whereas the other import mass flow rate was used to imitate the uncertainty and fluctuations during sinter production. The behavior of air and sinter outlet temperatures and pressure drop in the bed layer increased with an increase in the sinter import mass flow rate, as shown in Figure 9. The outlet temperature and pressure drop increased linearly with the increase in the sinter import mass flow rate. For each 5 t/h increase in the sinter import mass flow rate, the air and sinter outlet temperatures increased by a mean of 5.43 K and 4.08 K, respectively, and the loss of air resistance increased by a mean of 21 Pa. As the sinter import mass flow rate increased, a greater amount of thermal energy was carried in the sinter bed and increased the air and sinter temperatures, leading to increases in the superficial velocity of air in the sinter bed and the loss of air resistance. The increases in enthalpy exergy Ex and temperature exergy Ex(T) with the increase in the sinter inlet mass flow rate are summarized in Figure 10. An increase in the sinter import mass flow rate might lead to an increase in the air outlet temperature and superficial velocity. The difference between the temperature exergy and enthalpy exergy remained constant, approximately from 288 t/h to 360 t/h. For each 10 t/h increase in the sinter mass flow rate, the temperature and enthalpy exergy increased by a mean value of 1.12 MJ·s −1 and 1.13 MJ·s −1 , respectively. Figure 10 also presents that the pressure exergy Ex(p) decreased as the sinter import mass flow rate increased. It was negatively correlated with the pressure drop, as the pressure drop of airflow increased with the sinter import mass flow rate, as reported in previous studies. Therefore, the increase in gas resistance loss led to a decrease in pressure exergy Ex(p).

Effect of the Diameter of the Sinter Bed
The variations in air and sinter outlet temperatures and pressure drop of airflow with the diameter of the sinter bed are presented in Figure 11. As the diameter of the sinter bed increased, the air outlet temperature gradually increased and showed a decelerating rate. While the sinter outlet The increases in enthalpy exergy Ex and temperature exergy Ex(T) with the increase in the sinter inlet mass flow rate are summarized in Figure 10. An increase in the sinter import mass flow rate might lead to an increase in the air outlet temperature and superficial velocity. The difference between the temperature exergy and enthalpy exergy remained constant, approximately from 288 t/h to 360 t/h. For each 10 t/h increase in the sinter mass flow rate, the temperature and enthalpy exergy increased by a mean value of 1.12 MJ·s −1 and 1.13 MJ·s −1 , respectively. Figure 10 also presents that the pressure exergy Ex(p) decreased as the sinter import mass flow rate increased. It was negatively correlated with the pressure drop, as the pressure drop of airflow increased with the sinter import mass flow rate, as reported in previous studies. Therefore, the increase in gas resistance loss led to a decrease in pressure exergy Ex(p). The increases in enthalpy exergy Ex and temperature exergy Ex(T) with the increase in the sinter inlet mass flow rate are summarized in Figure 10. An increase in the sinter import mass flow rate might lead to an increase in the air outlet temperature and superficial velocity. The difference between the temperature exergy and enthalpy exergy remained constant, approximately from 288 t/h to 360 t/h. For each 10 t/h increase in the sinter mass flow rate, the temperature and enthalpy exergy increased by a mean value of 1.12 MJ·s −1 and 1.13 MJ·s −1 , respectively. Figure 10 also presents that the pressure exergy Ex(p) decreased as the sinter import mass flow rate increased. It was negatively correlated with the pressure drop, as the pressure drop of airflow increased with the sinter import mass flow rate, as reported in previous studies. Therefore, the increase in gas resistance loss led to a decrease in pressure exergy Ex(p).

Effect of the Diameter of the Sinter Bed
The variations in air and sinter outlet temperatures and pressure drop of airflow with the diameter of the sinter bed are presented in Figure 11. As the diameter of the sinter bed increased, the air outlet temperature gradually increased and showed a decelerating rate. While the sinter outlet

Effect of the Diameter of the Sinter Bed
The variations in air and sinter outlet temperatures and pressure drop of airflow with the diameter of the sinter bed are presented in Figure 11. As the diameter of the sinter bed increased, the air outlet temperature gradually increased and showed a decelerating rate. While the sinter outlet temperature steadily decreased, a flatter decreasing trend was observed. With the increase in the diameter of the sinter bed, an obvious increase in the decreasing velocity of the sinter flow rate was detected at the constant sinter production yield. The cooling time between air and sinter became longer, and the heat transfer was complete in the sinter bed. Thus, the air outlet temperature increased, and the sinter outlet temperature decreased. On the other hand, with the increase in the diameter of the sinter bed, a significant decrease in the superficial velocity of air was observed, which reduced the efficiency of heat transfer between the air and sinter. Thus, the cooling rate of hot sinter decreased, leading to a decreasing trend and deceleration of the sinter outlet temperature. As the diameter of the sinter bed continued to increase, the decreasing trend of sinter outlet temperature was slowing, and the increasing trend of the air outlet temperature decreased as well. Notably, as the diameter increased from 6 to 11 m, the air outlet temperature increased by 14.2 K per meter, while the sinter outlet temperature decreased by 4 K per meter. As further shown in Figure 11, the pressure drop of airflow was substantially reduced as the diameter of the sinter bed increased. Because the superficial velocity of air decreased as the diameter of the sinter bed increased, as mentioned above, it further hindered the pressure drop of airflow. For each 1 m increase in the diameter, the average pressure drop increased by 4427 Pa.
Energies 2020, 13, 2225 15 of 30 temperature steadily decreased, a flatter decreasing trend was observed. With the increase in the diameter of the sinter bed, an obvious increase in the decreasing velocity of the sinter flow rate was detected at the constant sinter production yield. The cooling time between air and sinter became longer, and the heat transfer was complete in the sinter bed. Thus, the air outlet temperature increased, and the sinter outlet temperature decreased. On the other hand, with the increase in the diameter of the sinter bed, a significant decrease in the superficial velocity of air was observed, which reduced the efficiency of heat transfer between the air and sinter. Thus, the cooling rate of hot sinter decreased, leading to a decreasing trend and deceleration of the sinter outlet temperature. As the diameter of the sinter bed continued to increase, the decreasing trend of sinter outlet temperature was slowing, and the increasing trend of the air outlet temperature decreased as well. Notably, as the diameter increased from 6 to 11 m, the air outlet temperature increased by 14.2 K per meter, while the sinter outlet temperature decreased by 4 K per meter. As further shown in Figure 11, the pressure drop of airflow was substantially reduced as the diameter of the sinter bed increased. Because the superficial velocity of air decreased as the diameter of the sinter bed increased, as mentioned above, it further hindered the pressure drop of airflow. For each 1 m increase in the diameter, the average pressure drop increased by 4427 Pa. Under this condition, the sinter mass flow rate (100 kg·s −1 ), effective height of the sinter bed (6 m), the air inlet temperature (293 K), air mass flow rate (140 kg/s), and sinter inlet temperature (973 K) were fixed. The variations in the enthalpy exergy Ex, temperature exergy Ex(T), and pressure exergy Ex(p) were obtained for various diameters of the sinter bed, as shown in Figure 12. The Ex, Ex(T), and Ex(p) increased as the diameter of the sinter bed increased. An increase in the diameter likely increased the air outlet temperature and dramatically reduced the superficial velocity of air. As the pressure drop of airflow decreased, the pressure exergy Ex(p) in the sinter bed increased substantially. At the same time, an increase in the air outlet temperature led to a simultaneous increase in the temperature exergy Ex(T) and enthalpy exergy Ex. For each 1 m increase in diameter, the temperature exergy Ex(T) and enthalpy exergy Ex increased by an average of 0.8 MJ·s −1 and 1.4 MJ·s −1 , respectively. Furthermore, the increasing trend of pressure exergy Ex(p) was accelerated, which further limited the difference between the enthalpy exergy and temperature exergy as the diameter of the sinter bed steadily increased. Under this condition, the sinter mass flow rate (100 kg·s −1 ), effective height of the sinter bed (6 m), the air inlet temperature (293 K), air mass flow rate (140 kg/s), and sinter inlet temperature (973 K) were fixed. The variations in the enthalpy exergy Ex, temperature exergy Ex(T), and pressure exergy Ex(p) were obtained for various diameters of the sinter bed, as shown in Figure 12. The Ex, Ex(T), and Ex(p) increased as the diameter of the sinter bed increased. An increase in the diameter likely increased the air outlet temperature and dramatically reduced the superficial velocity of air. As the pressure drop of airflow decreased, the pressure exergy Ex(p) in the sinter bed increased substantially. At the same time, an increase in the air outlet temperature led to a simultaneous increase in the temperature exergy Ex(T) and enthalpy exergy Ex. For each 1 m increase in diameter, the temperature exergy Ex(T) and enthalpy exergy Ex increased by an average of 0.8 MJ·s −1 and 1.4 MJ·s −1 , respectively. Furthermore, the increasing trend of pressure exergy Ex(p) was accelerated, which further limited the difference between the enthalpy exergy and temperature exergy as the diameter of the sinter bed steadily increased.

Effect of the Height of the Sinter Bed
The effects of the height of the cooling sinter bed on the heat transfer and fluid flow of the sinter bed were analyzed. The variations in air and sinter outlet temperatures and pressure drop at various heights of the sinter bed are illustrated in Figure 13. When the height of the sinter bed increased, the air outlet temperature steadily increased, and a steady increasing trend was observed: the air outlet temperature increased by 67.8 K per meter. Furthermore, the sinter outlet temperature smoothly decreased, and a gradually decreasing trend was observed: the sinter outlet temperature decreased by 12.3 K per meter. As the height increased, the cooling time of heat transfer between the air and sinter would be prolonged, resulting in more sufficient heat transfer in the sinter bed, an increasing trend in the air outlet temperature, and a decreasing trend in the sinter outlet temperature. Moreover, the local temperature difference between the air and the sinter increased as the height of the sinter bed increased. A weaker heat transfer between air and sinter was observed, leading to a deceleration of the decreasing trend in the sinter outlet temperature. As shown in Figure 13, the pressure drop of the airflow stably increased when the height of the sinter bed increased. For each 1 m increase in height, the pressure drop of airflow increased by 1488.8 Pa.

Effect of the Height of the Sinter Bed
The effects of the height of the cooling sinter bed on the heat transfer and fluid flow of the sinter bed were analyzed. The variations in air and sinter outlet temperatures and pressure drop at various heights of the sinter bed are illustrated in Figure 13. When the height of the sinter bed increased, the air outlet temperature steadily increased, and a steady increasing trend was observed: the air outlet temperature increased by 67.8 K per meter. Furthermore, the sinter outlet temperature smoothly decreased, and a gradually decreasing trend was observed: the sinter outlet temperature decreased by 12.3 K per meter. As the height increased, the cooling time of heat transfer between the air and sinter would be prolonged, resulting in more sufficient heat transfer in the sinter bed, an increasing trend in the air outlet temperature, and a decreasing trend in the sinter outlet temperature. Moreover, the local temperature difference between the air and the sinter increased as the height of the sinter bed increased. A weaker heat transfer between air and sinter was observed, leading to a deceleration of the decreasing trend in the sinter outlet temperature. As shown in Figure 13, the pressure drop of the airflow stably increased when the height of the sinter bed increased. For each 1 m increase in height, the pressure drop of airflow increased by 1488.8 Pa.

Effect of the Height of the Sinter Bed
The effects of the height of the cooling sinter bed on the heat transfer and fluid flow of the sinter bed were analyzed. The variations in air and sinter outlet temperatures and pressure drop at various heights of the sinter bed are illustrated in Figure 13. When the height of the sinter bed increased, the air outlet temperature steadily increased, and a steady increasing trend was observed: the air outlet temperature increased by 67.8 K per meter. Furthermore, the sinter outlet temperature smoothly decreased, and a gradually decreasing trend was observed: the sinter outlet temperature decreased by 12.3 K per meter. As the height increased, the cooling time of heat transfer between the air and sinter would be prolonged, resulting in more sufficient heat transfer in the sinter bed, an increasing trend in the air outlet temperature, and a decreasing trend in the sinter outlet temperature. Moreover, the local temperature difference between the air and the sinter increased as the height of the sinter bed increased. A weaker heat transfer between air and sinter was observed, leading to a deceleration of the decreasing trend in the sinter outlet temperature. As shown in Figure 13, the pressure drop of the airflow stably increased when the height of the sinter bed increased. For each 1 m increase in height, the pressure drop of airflow increased by 1488.8 Pa.  The enthalpy exergy Ex and temperature exergy Ex(T) increased when the height of the sinter bed increased, and pressure exergy Ex(p) decreased as the height of the sinter bed increased, as shown in Figure 14. As the height of the sinter bed increased, the increasing rate of enthalpy exergy Ex slowed, but the pressure exergy Ex(p) decreased as the pressure drop of airflow increased. Furthermore, the increase in the rate of temperature exergy Ex(T) decelerated, which led to a greater difference between Ex and Ex(T). As a result, the increasing trend of the enthalpy exergy Ex of air decreased further. The enthalpy exergy Ex and temperature exergy Ex(T) increased when the height of the sinter bed increased, and pressure exergy Ex(p) decreased as the height of the sinter bed increased, as shown in Figure 14. As the height of the sinter bed increased, the increasing rate of enthalpy exergy Ex slowed, but the pressure exergy Ex(p) decreased as the pressure drop of airflow increased. Furthermore, the increase in the rate of temperature exergy Ex(T) decelerated, which led to a greater difference between Ex and Ex(T). As a result, the increasing trend of the enthalpy exergy Ex of air decreased further.

Taguchi Method and ANOVA
The Taguchi method is deemed as one of the popularly applied technical tools developed by Taguchi, which focuses on optimizing design experiments to determine the best factors and levels affecting the target response by reducing tedious and uneconomical computations and to obtain the order of importance of different factors for an objective response. The Taguchi method is the most effective optimization method to minimize the number of design experiments and contains fullfactorial information from the fewest experiments with the help of orthogonal arrays. In this study, orthogonal arrays described the particular arrangement of factorial design experiments in terms of the factors being analyzed. Several parameters and their combinations were present in the vertical sinter bed that modulated the efficiency of the heat transfer process and must be analyzed to comprehensively understand the effects of factors of interest on the heat transfer and hydrodynamic performances in the sinter cooling process. Therefore, the factors of interest and the range of levels are presented in Table 6.
The choice of an appropriate orthogonal array is an essential step of a statistical analysis, which is determined by the degree of freedom (DOF) theory [33]. The calculations from the L25 (5 5 ) orthogonal experiments were applied to identify the effects of factors and their levels on the target value. The corresponding results and signal to noise ratios are described in Table 7. According to previous studies, the air inlet temperature, sinter mass flow rate, height of the sinter bed, diameter of the sinter bed, and air mass flow rate were used for optimization.
According to the Taguchi method, three criteria determine optimization: higher is better, normal is better, and lower is better. In the present study, the higher is better criterion was applied to describe the heat transfer performance, and the lower is better criterion was employed to present the hydrodynamic performance in the sinter cooling process [34].

Taguchi Method and ANOVA
The Taguchi method is deemed as one of the popularly applied technical tools developed by Taguchi, which focuses on optimizing design experiments to determine the best factors and levels affecting the target response by reducing tedious and uneconomical computations and to obtain the order of importance of different factors for an objective response. The Taguchi method is the most effective optimization method to minimize the number of design experiments and contains full-factorial information from the fewest experiments with the help of orthogonal arrays. In this study, orthogonal arrays described the particular arrangement of factorial design experiments in terms of the factors being analyzed. Several parameters and their combinations were present in the vertical sinter bed that modulated the efficiency of the heat transfer process and must be analyzed to comprehensively understand the effects of factors of interest on the heat transfer and hydrodynamic performances in the sinter cooling process. Therefore, the factors of interest and the range of levels are presented in Table 6. The choice of an appropriate orthogonal array is an essential step of a statistical analysis, which is determined by the degree of freedom (DOF) theory [33]. The calculations from the L 25 (5 5 ) orthogonal experiments were applied to identify the effects of factors and their levels on the target value. The corresponding results and signal to noise ratios are described in Table 7. According to previous studies, the air inlet temperature, sinter mass flow rate, height of the sinter bed, diameter of the sinter bed, and air mass flow rate were used for optimization. According to the Taguchi method, three criteria determine optimization: higher is better, normal is better, and lower is better. In the present study, the higher is better criterion was applied to describe the heat transfer performance, and the lower is better criterion was employed to present the hydrodynamic performance in the sinter cooling process [34]. S/N = 10 log 1 n y 2 (17) S/N = −10 log 1 n y 2 (18) where y is the mean value of the response function, and n is the number of factors analyzed in the experiment.
A response table was constructed to describe the mean S/N (signal to noise) ratio and the significance level of factors and to obtain a comprehensive understanding of the effects of design factors on the response function, as shown in Table 8, for air outlet exergy. As described above, the criteria used in the statistical method to analyze the heat transfer performance was higher is better. Therefore, the greatest variation in the range of the S/N ratio (delta value) represented the most important factor in the corresponding response function. Thus, the best performance level of each factor had the largest value of the S/N ratio. As shown in Table 8, the most effective factor was the diameter of the bed, followed by the airflow rate and height of the bed. The optimal combination of factors for the air outlet exergy response was A 1 B 5 C 5 D 5 E 2 . The Taguchi method was used to analyze the effects of design factors on the pressure drop as a response. In Table 9, the diameter of the bed (D b ) was the main factor contributing to the response compared with other factors, and a better performance level was attributed to the higher S/N ratio among the five levels for each factor. Thus, the optimized combination of factors and levels were determined as A 1 B 5 C 1 E 5 D 1 for the pressure drop as a response function. ANOVA is a statistical analysis method developed by Ronald Fisher that divides aggregate variations in the observed data (response measurement) into two portions: systematic portions and random portions. The systematic portions and random portions are determined by varying factors being tested in the experiments. ANOVA is applied to determine the effects of independent variables (the various factors of interest) on the dependent variable (the response of physical phenomenon) in statistical analysis. Tables 10 and 11 present the results of the ANOVA for air outlet exergy and pressure drop, respectively. The results of ANOVA are presented in tables, and the effect of design factors on the air outlet exergy was reported as F-values and 95% confidence intervals. In this situation, if F 0.05 (4.4) > 6.388, according to the deciding value for the 95% confidence interval and 5% significance level, the amplitude of the variation in the tested factors would affect the air outlet exergy and pressure drop. The criteria for determining the effect of a design factor on response function mainly depends on the magnitude of the F-value. The design factor with the highest F-value represents the factor with the most significant effect on the response function.  The mathematical description of the ANOVA, such as a sum of squares (SS), mean of squares (MS), degree of freedom (DOF), and F-values, are calculated using Equations (19)- (26): (26) where CF is the correction factor, y i is the factors tested in the i th case, m is the total number of cases, SS factor is the total of squares of a single factor, y factor-i is the i th level of the factor of interest, N is the total number of levels of the factor of interest, SS T is the sum of squares of variation, SS e is the total of squares due to error, SS A is the sum of squares of factor A, f total is the total DOF, n is the number of standard experiments, f factor is the DOF for a single factor, k is the total number of levels of a single factor, F factor is the F-value of factor, and V factor and V error are the sum of squares of a single factor and error function, respectively. In this section, ANOVA was applied to analyze the effects of various design factors on a single response function. The method introduced in this study had two advantages: (1) the contribution rate and most significant parameters were obtained by analyzing the variations in design factors, and (2) the reliability of the results of the Taguchi analysis method was verified. As shown in Table 10, the diameter of the bed was the most significant factor affecting the heat transfer performance, with a contribution rate of 56.79%, followed by the height of the bed and air inlet temperature, with contribution rates of 16.32% and 9.37%, respectively. The air mass flow rate exerted little effect on the exergy recovered. According to the magnitude of the contribution rate for each factor, the significance of the design factors was determined as D b > H b > T g,in > M s > M a , consistent with the results acquired from the Taguchi method. The quantitative description of the effect of each factor on the air outlet exergy and contribution rates are presented in Figure 15. In addition, the contribution rates of design factors to the response of pressure drop determined using ANOVA are presented in Table 11. As shown in Table 11, the diameter of the bed exerted the greatest effect on the pressure drop, with a contribution rate of 57.08%, followed by the airflow rate (29.22%) and air inlet temperature (5.66%). The sinter flow rate and height of the bed exerted a limited effect on the hydrodynamic performance of the vertical sinter cooling bed. The contribution rate of error was 0.95%, which comprised of three error parts, such as calculation error, external sources of error, and uncontrollable error. The contribution of each factor to air pressure is presented in Figure 16. error parts, such as calculation error, external sources of error, and uncontrollable error. The contribution of each factor to air pressure is presented in Figure 16.

Grey Relational Analysis (GRA)
In the actual industry, production is usually required to meet multi-objective target functions, and thus the Taguchi method maybe not suitable to optimize the multi-objective performance. The Taguchi grey relational method is applied to determine the order of importance of design factors in modulating the performance of multiple objective characteristics, reveal the effect of geometric and operating factors on the heat transfer performance characteristics of a vertical sinter cooling bed, and overcome the limitation of the Taguchi analysis that is only suitable for a single response function. The advance of the GRA method is the ability to convert the bi-objective or multi-objective target optimization into a single target optimization problem. In this section, the GRA method was employed for optimization to simultaneously improve the heat transfer and hydrodynamic performance in the vertical sinter cooling process.
Based on the published literature, five factors-the airflow mass flow rate, air inlet temperature, sinter mass flow rate, height of the bed, and diameter of the bed-were chosen to optimize the performance characteristics. Two criteria had been established for the simultaneous optimization analysis: higher is better for describing heat transfer performance (air outlet exergy), and the lowest is better for the hydrodynamic performance (air pressure drop). The mathematical description of criteria for higher is better and the lowest is better is shown in Equation (27) and Equation (28), respectively.
The value of the grey relational grade calculated from the grade relational coefficient for each response must be calculated. The grade relational coefficient determines the relation between comparative and referenced normalized value. Equations (28) In the actual industry, production is usually required to meet multi-objective target functions, and thus the Taguchi method maybe not suitable to optimize the multi-objective performance. The Taguchi grey relational method is applied to determine the order of importance of design factors in modulating the performance of multiple objective characteristics, reveal the effect of geometric and operating factors on the heat transfer performance characteristics of a vertical sinter cooling bed, and overcome the limitation of the Taguchi analysis that is only suitable for a single response function. The advance of the GRA method is the ability to convert the bi-objective or multi-objective target optimization into a single target optimization problem. In this section, the GRA method was employed for optimization to simultaneously improve the heat transfer and hydrodynamic performance in the vertical sinter cooling process.
Based on the published literature, five factors-the airflow mass flow rate, air inlet temperature, sinter mass flow rate, height of the bed, and diameter of the bed-were chosen to optimize the performance characteristics. Two criteria had been established for the simultaneous optimization analysis: higher is better for describing heat transfer performance (air outlet exergy), and the lowest is better for the hydrodynamic performance (air pressure drop). The mathematical description of criteria for higher is better and the lowest is better is shown in Equation (27) and Equation (28), respectively.
The value of the grey relational grade calculated from the grade relational coefficient for each response must be calculated. The grade relational coefficient determines the relation between comparative and referenced normalized value. Equations (28)-(31) are used to describe the grade relational coefficient: ∆ min = min ∀j∈i min ∀k y 0 (k) − y i (k) (32) where ∆ 0i is the value of the deviation between comparative and reference sequences, y i (k) and y 0 (k) are the comparative and references sequences, respectively. ∆ min and ∆ min are the minimum and maximum values of ∆ 0i , respectively. ρ is the identification coefficient that ranges from 0 to 1. According to previous studies using GRA, the use of an arbitrary value for the identification coefficient ranging from 0 to 1 does not affect the sequence of the significance of the factors. Consistent with the conclusions of many research papers, ρ was equal to 0.5 in this study. Notably, the applicable method for determining the weight factor for each response function in GRA that is significant obtains the reliability of the results of the grey relational grade. In a previous study, factors were assigned equal weight factors in the multiple response functions to describe the relation of the grey relational coefficient for different responses, but it may not be a reliable method. In actual industrial production, the weight of each factor contributing to multiple physical phenomena differs. Therefore, a reliable method must be introduced to determine the weight factors for different responses based on the effects of various factors of interest on the performance in the vertical sinter cooling process.
where γ i is the grey relational grade for each combination orthogonal experiment, w k is the weight factor for kth response function, and n is the number of response functions. The mean S/N ratio scale of each factor represents the order of significance in the statistical analysis. The calculated delta value is applied to determine the most effective factor for each response function that is equal to the largest value minus the smallest value of the S/N ratio for each factor. The delta value for each factor is calculated to obtain the sum of delta values for each response function. The weight factor for each response is determined by calculating the ratio of the total of delta value for each response to the sum of delta values for all responses. The mathematical process is described in Equation (34): where m is the total number of response functions (in this study, the total number of responses was two), n is the number of the factors for each response function, and Delta describes the range from largest to smallest values of the S/N ratio for each factor.
According to Equation (34), the weight factor is calculated for each response function. The weight factors of air outlet exergy and pressure drop were 48.71% and 51.29%, respectively. The weight factor for air outlet exergy and pressure drop are listed in Tables 8 and 9, respectively. The grey relational grade is described as follows by considering weight factors: where GRC P and GRC Ex represent the grey relational coefficient of the pressure drop and air outlet exergy, respectively. As mentioned above, a significant advantage of the grey relational grade method is that it converts the two response functions (increasing air outlet exergy and reducing pressure drop) to a single target function (multi-objective performance). The GRA method is used for the optimization process to simultaneously improve multiple performance characteristics of a vertical sinter cooling bed. Table 12 shows the calculated values for normalized results, grey relational coefficients, grey relational grades, and the order of significance for GRA. The grey relational grade is calculated using Equation (34). The variations in the grey relational grade and corresponding S/N ratios are described in Figure 17. The case combining all considered factors, which exhibited the maximum grey relational grade, performed the best in terms of multiple performance characteristics. As shown in Table 7, case 5 (A 1 B 5 C 5 D 5 E 5 ) exhibited the best performance for multiple characteristics. Under the condition of the combination of factors included in case 5, the air outlet exergy and pressure drop of the vertical sinter cooling bed were 32.2 MJ·s −1 and 8328.9 Pa, respectively. In addition, the smallest grey relational grade was obtained for the combination of factors in case 23 (A 5 B 3 C 2 D 1 E 5 ). In this situation, the air outlet exergy and pressure drop of the cooling bed were 4.08 MJ·s -1 and 27668 Pa, respectively. Equation (34). The variations in the grey relational grade and corresponding S/N ratios are described in Figure 17. The case combining all considered factors, which exhibited the maximum grey relational grade, performed the best in terms of multiple performance characteristics. As shown in Table 7, case 5 (A1B5C5D5E5) exhibited the best performance for multiple characteristics. Under the condition of the combination of factors included in case 5, the air outlet exergy and pressure drop of the vertical sinter cooling bed were 32.2 MJ·s −1 and 8328.9 Pa, respectively. In addition, the smallest grey relational grade was obtained for the combination of factors in case 23 (A5B3C2D1E5). In this situation, the air outlet exergy and pressure drop of the cooling bed were 4.08 MJ·s -1 and 27668 Pa, respectively. A response table was generated for the grey relational grade to describe the order of significance for each factor. As shown in Table 13, the order of importance of each factor was Db > Ma> Tg,in> Ms> Hb. Each factor and the levels resulting in the best and worst performance of multiple processes in the vertical cooling bed are shown in Figure 18. The diameter of the bed (Db) was the most important factor, followed by the air mass flow rate (Ma) in terms of the effect of the tested factors on the performance of multiple processes in the cooling bed, as shown in Figure 18.
ANOVA was applied to reveal the contribution rates of the tested factors and verify the result of the response table for the grey relational grade. The results of the AVONA are shown in Table 14, and the diameter of the bed was the most important factor, followed by air mass flow rate and air inlet temperature for determining the performance of multiple processes. The height of the cooling bed and sinter mass flow exerted limited effects on multiple levels of performance of the cooling bed. The contributions of each factor to the grey relational grade are described in Figure 19. The results of the ANOVA were consistent with the response table (Table 13) for the grey relational grade.
In addition, the regression analysis of the equation, including the design factors for GRG (grey relational grade), is expressed in Equation (36). The R 2 value of 94.67% was obtained from the regression equation for the performance of multiple processes, and the standard deviation was 0.03687.  Table 12 and their S/N (signal to noise) ratio.  Table 12 and their S/N (signal to noise) ratio.
A response table was generated for the grey relational grade to describe the order of significance for each factor. As shown in Table 13, the order of importance of each factor was D b > M a > T g,in > M s > H b . Each factor and the levels resulting in the best and worst performance of multiple processes in the vertical cooling bed are shown in Figure 18. The diameter of the bed (D b ) was the most important factor, followed by the air mass flow rate (M a ) in terms of the effect of the tested factors on the performance of multiple processes in the cooling bed, as shown in Figure 18.    ANOVA was applied to reveal the contribution rates of the tested factors and verify the result of the response table for the grey relational grade. The results of the AVONA are shown in Table 14, and the diameter of the bed was the most important factor, followed by air mass flow rate and air inlet temperature for determining the performance of multiple processes. The height of the cooling bed and sinter mass flow exerted limited effects on multiple levels of performance of the cooling bed. The contributions of each factor to the grey relational grade are described in Figure 19. The results of the ANOVA were consistent with the response table (Table 13) for the grey relational grade.  Figure 18. The effect of each factor on the performance of the two responses.  Figure 19. Contributions of each factor to the grey relational grade.

Confirmation Test
In this section, the confirmation test was an important step to verify the optimization of multiple responses. Table 15 shows a comparison of the effects of the optimal factor combination (A 1 B 5 C 5 D 5 E 1 ) and initial factor combination (A 1 B 5 C 5 D 5 E 5 ) for multiple performance characteristics. The confirmation experiment, including the optimal combination of design parameters, resulted in an exergy of outlet air and air pressure drop of 29.89 MJ·s −1 and 2412.5 Pa, respectively. The increases in the GRG value and corresponding percentage compared with the initial settings were 0.0727 and 8.51%, respectively. Thus, the heat transfer and hydrodynamic performances of the sinter cooling process in the VSCB were improved with the aid of the Taguchi GRA method. In addition, the utilization of waste heat was increased, and energy loss was reduced.

Conclusions
The CFD numerical simulation, coupled with the Taguchi-grey relational analysis method was applied in the present study to investigate the effects of geometric and operational factors on heat transfer and hydrodynamic performance. The exergy of air outlet recovery and the air pressure for the heat transfer process in the VSCB were selected as the response performance functions. The main conclusions are listed below.
The effect of design factors on the exergy of air outlet recovery was studied. The results revealed that the diameter and height of bed were the most important factors, with a contribution rate of 56.79% and 16.32%, respectively. The optimal combination of factors for the exergy of air outlet recovery was A 1 B 5 C 5 D 5 E 2 . In addition, the maximum exergy of recovery value was 29.68 MJ·s −1 , based on the optimal situation from the thermodynamics analysis.
The impact of design factors on air pressure was evaluated. The optimal condition was A 1 B 5 C 1 D 5 E 1 , and the diameter of the bed and air mass flow rate were the most important factors affecting the hydrodynamic performance, with contribution rates of 57.08% and 29.22%, respectively. Furthermore, the air pressure was simulated as 3703.8 Pa, based on the combination of optimal factors from the hydrodynamics analysis.
The influence of geometric and structural factors on both of the maximum exergy of air outlet recovery and minimal air pressure drop was investigated. The optimal factor combination was A 1 B 5 C 5 D 5 E 1 , while the diameter of the bed and air mass flow rate exerted the dominant effects on multiple performance measures, with contributions of 70.51% and 15.84%, respectively. The optimum values for multiple performance characteristics based on the best combination were 29.89 MJ·s −1 and 2412.5 Pa, respectively.
The confirmation test was performed to compare the optimum condition and initial condition based on the Taguchi with a grey relational analysis method to improve the heat transfer and hydrodynamic performances. The optimum results revealed increases in the GRG value and corresponding ratio by 0.0727 and 8.51%, respectively, compared with the initial design.   Parameter of interest