Simulation research of a counter-flow rotary kiln hazardous waste incineration system

As industrialization accelerates and the amount of hazardous waste generated gradually increases, the means of disposal of hazardous waste is of increasing concern. In this paper, a 40 t/d counter-flow rotary kiln incineration system owned by a Jiangsu environmental protection company was researched. The software Aspen Plus was used to build the mixed pyrolysis model and the software Fluent was used to build the computational fluid dynamics model of the incineration system. The influence of the calorific value of the hazardous waste, the operating temperature and the air supply on the operational effectiveness of the incineration system were analyzed by varying the simulation conditions. The results show that the SOx and NOx content of the product is lower when the operating temperature is above 800 °C. The incineration system could only operate above 800 °C when the calorific value of the hazardous waste is not less than 1500 kcal/kg. The incineration system operated best at a primary air velocity of 1.5 m/s. The simulation results in this paper serve as a guide for the operation of counter-flow rotary kiln incineration systems.


Introduction
Since the reform and opening-up, all industries in China have been developing rapidly and the economic level has greatly improved. However, the issue of environmental pollution is also increasingly serious (Nanda and Berruti 2021). Regarding pollutants, hazardous waste is considered as one of the most difficult wastes to address due to its reactivity, flammability, toxicity and corrosiveness. In addition, hazardous waste is extremely harmful to both people and the environment, whereby it can endanger human health (Tekade 2021) and cause serious damage to soil (Wang and Geng 2015) and water bodies (Hadas et al. 2021;Visvanathan 1996). Therefore, an efficient and safe disposal of hazardous waste is essential.
Hazardous waste treatment methods include physical and chemical methods, solidification, safe landfill, and incineration (Makarichi et al. 2018;Paula et al. 1998;Singhabhandhu and Tezuka 2010). Among these, physical and chemical methods cannot completely treat hazardous waste and still require a conjunction with other treatment methods (Huang et al. 1993;Murphy and McAllister 2001). Solidification and safe landfill methods are simple to operate, however, they have long treatment cycles and cannot achieve a rapid and effective treatment of large quantities of hazardous waste (Conner and Hoeffner 1998;Malviya and Chaudhary 2006). Incineration is the most widely used method due to its high efficiency and thorough treatment. The choice of a suitable incinerator is a key factor in determining the effectiveness of incineration. In order to adapt to the complex and variable composition of hazardous waste, rotary kiln incinerators, with a wide range of applications, are often used. Liu et al. (2005) used experimental and simulation methods to research the main factors affecting the movement characteristics of hazardous waste in a rotary kiln. The results showed that the rotational speed and inclination angle of the rotary kiln affect the residence time and filling rate of hazardous waste in the rotary kiln, which, in turn, affects the actual operating effect of the rotary kiln. Andréa et al. (1998) researched the practical effect of coprocessing hazardous waste treatment in cement kilns and showed that hazardous waste can act as a secondary fuel required for cement kiln incineration, where the incineration efficiency of hazardous waste was greater than 99%. The durability of the refractory in the rotary kiln is the key to the operational life of the rotary kiln. Adrian et al. (2016) conducted an experimental study on the causes of refractory brick breakage in rotary kilns, which showed that the breakage of refractory bricks was not due to thermochemical corrosion, but rather to the penetration of molten material into the refractory bricks in the rotary kiln. Then, they developed a mullite-zirconia composite as a new refractory material, which had excellent durability (Adrian et al. 2021). Zhang and Wang (2022) experimentally investigated the operating effect of rotary kilns under different hazardous waste dispensing conditions and showed that the lower the ignition temperature of the hazardous waste, the more stable the operation of the rotary kiln. Wajda et al. (2022) provided an algorithmic procedure for solving the optimal operating conditions of rotary kilns to achieve the highest energy use efficiency. A comparative analysis of the evolution of hydrocarbon composition in pilot rotary kiln incinerators and rotary kiln incinerators was carried out by Lester et al. (1991). The results showed that the experimental results in the pilot reactor were consistent with the actual operation of the rotary kiln incinerator and that the pilot system exhibited a good response to the actual operation of the rotary kiln. The incineration of hazardous waste in the rotary kiln produced dioxins and other harmful substances.
The above studies were very informative, though they did not provide a good explanation for the incomplete burning of hazardous waste and the uneven temperature distribution in rotary kilns, which often occurs in rotary kilns. The matching of the primary air volume to the calorific value of the hazardous waste was the key to the uniformity of the rotary kiln temperature, which affected the incineration of hazardous waste. The software Aspen Plus was widely used to simulate pyrolysis and gasification processes of coal and solid waste (Ismail et al. 2020). In addition, the operating temperature of the rotary kiln and the secondary combustion chamber affected pollutant emissions. Hence, a 40 t/d counter-flow rotary kiln incineration system owned by a Jiangsu environmental protection company was researched. The software Aspen Plus was used to model the mixed pyrolysis of hazardous waste in a rotary kiln and to analyze the content of pollutants such as SO x and NO x in the gaseous products at different temperatures. The software Fluent was used to simulate the temperature distribution and the concentration distribution of various substances in the rotary kiln and second combustion chamber. The effect of primary air flow and the calorific value of the hazardous waste on the operation of the rotary kiln was also researched. The simulation results were of certain significance for the actual operation of rotary kilns. Figure 1 shows a schematic diagram of a counter-flow rotary kiln with a horizontal inclination of 2.5° and a rotational speed of 0.2-0.5 r/min. Hazardous waste was used as the incineration material. The component content of the hazardous waste is given in Table 1. The hazardous waste entered the rotary kiln from the kiln head and moved from the kiln head to the kiln tail under the rotation of the rotary kiln, during which the hazardous waste underwent drying, pyrolysis and combustion processes in turn, and the ash was discharged from the kiln tail. The primary air entered the rotary kiln at the end of the kiln, where it was first preheated by the ash, and then participated in the combustion, pyrolysis and drying processes of the hazardous waste in turn, before forming a mixture that entered the second combustion chamber.

Aspen plus modeling
In this paper, a hybrid pyrolysis model for the counter-flow rotary kiln incineration system was simulated using the Aspen Plus software. Three processes of drying, pyrolysis and combustion occurred sequentially in the rotary kiln, so the hybrid pyrolysis model was built by first creating a single model for each process and then connecting the single models based on the direction of material and energy flow between the models. As shown in Fig. 2, the model consisted of 13 unit operation modules, 20 material flow units (black solid line) and 6 energy flow units (red dashed line), which enabled a counter-flow reaction process between the hazardous waste and primary air in the kiln. The energy was supplied from the combustion model to the drying and pyrolysis models through the counter flow of the energy. The model contained a total of five reactors and the operating conditions for the five reactors are given in Table 2. The components of hazardous waste were complex and the chemical reactions involved in the pyrolysis and incineration of hazardous waste were numerous. Therefore, a Gibbs reactor was used to simulate the generation of various gaseous and solid products, which followed the Gibbs free energy minimum principle.

Fluent modeling
In this paper, a CFD model of the counter-flow rotary kiln incineration system was developed using the software Fluent. The physical properties of the hazardous waste were defined in the UDF user-defined mode in Fluent, so that it existed as a "pseudo-fluid". By setting up the pyrolysis gasification reaction of the material in the flow field area of the kiln, the pyrolysis and combustion reactions of the hazardous waste occurred on a time scale appropriate to the kiln. The pyrolysis and combustion reactions of hazardous waste were adapted to the time scale of the kiln. The component transport model was used to simulate the flow of each component and pyrolysis gas in the kiln, resulting in the simultaneous simulation of the gas-solid two-phase flow, the pyrolysis process of solid materials and the combustion process of materials and pyrolysis gas throughout the kiln.

Geometric model and boundary conditions
The geometry of the counter-flow rotary kiln is shown in Fig. 3. It can be seen that the number of inlets and outlets of the rotary kiln was four, divided into the following four categories: material inlet, primary air inlet, mixed gas outlet and slag outlet. The boundary conditions for each inlet and outlet are given in Table 3.

Mathematical-physical model
In the CFD model, the standard k−ε model was chosen as the turbulent flow model and a component transport model was used to simulate the chemical reaction process. Only volumetric reactions were considered in the model and the effects of energy diffusion and thermal diffusion on the reaction were considered.
(1) Turbulence model. The equation of turbulent kinetic energy k was as follows: The equation of turbulent energy dissipation rate ε was as follows: where t is the time, ρ d was is the density of the mixture, v d is the speed of the mixture, t d is the turbulent viscosity of the mixture, G kd is the turbulent energy, and C 1e , C 2e , σ k and σ ε are turbulence model coefficients.
(1)  The reaction equations for material and fluid dynamics in a rotary kiln generally included four equations for the conservation of mass, the conservation of momentum, the conservation of energy and the conservation of component transport.
where ρ is the density, ρ' is the average density, u' is the instantaneous speed, p is the individual gas partial pressure, μ is the kinetic viscosity, H is the specific enthalpy of fluids, T is the temperature, λ is the thermal conductivity of fluids, Y f is the mass fraction, D is the diffusion coefficient, R f is the production rate of components per unit volume, and S p , S N , S H and S Y are source items.
(3) Chemical reaction model. Seven chemical reactions were initially set up in the simulation. The chemical reactions were calculated using the Arrhenius rate model set up with the equation: where K is the Arrhenius constant, A is the pre-exponential factor, Ea is the activation energy, and R is the molar gas constant.
The main reaction equations and kinetic parameters are given in Table 4.

Grid independence test and model validation
The grid independence of the CFD model was checked using the temperature distribution in the rotary kiln at 1.5 m/s as a test indicator. The grid shape was rectangular and the grid quality is given in Table 5. Figure 4 shows the results of the independence test for the grid. It can be observed that the temperature distribution was too far from the expected error when using 900,000 grids, while the simulation with 950,000 grids exhibited the most consistent temperature distribution, as expected. When the number of grids was increased to one million, the temperature distribution curve was also very close to the expected distribution, but at this point, the error from the expected distribution was already small. The number of grids should be as small as possible for the simulation to meet our expectations, so from this grid independence check, it can be determined that a grid of 950,000 should be chosen.

Fig. 4
Results of the grid independence test for the temperature field the CFD model of the counter-flow rotary kiln incineration system was accurate.

Calculation of mass balance and energy balance
The PR-BM method was chosen in Aspen Plus to simulate the pyrolysis incineration process of hazardous waste, and mass-energy accounting was carried out using the following method.
where, M loss represents the mass loss in rotary kilns; M hw and M air represent the mass flow of hazardous waste and air, respectively; M fg and M ash represent the mass flow of flue gas and ash, respectively; E represents the energy of the various components. To make the model more realistic for engineering purposes, a mass loss of 1% and a heat loss of 2% were considered in the aspen modelling. A mass and energy balance using the simulation results showed that M loss and E loss were 16.67 kg/h and 123,333 kcal/h, respectively. The above results satisfied the mass and energy balance.

Effect of operating temperature
The yields of reducing gases, sulfides and nitrides obtained from mixed pyrolysis model simulations at 500-900 °C are given in Fig. 5. It can be seen that with the increase of operating temperature, the contents of CO and H 2 in the mixed gas gradually increased, but the contents of CH 4 gradually decreased. With the increase of the operating temperature, the proportion of H 2 S in gaseous sulfides increased and the proportion of SO 2 decreased. The proportion of HCN and NH 3 increased and the proportion of NO decreased. SO 2 and NO were mainly produced by the combustion of fuel. The kiln tail had sufficient air and high temperatures to facilitate the production of sulfur SO x and NO x .
The drying and pyrolysis of hazardous waste occurred in the front and middle sections of the rotary kiln, while air was fed from the tail section of the kiln. Therefore, there was insufficient air during the pyrolysis of hazardous waste. Therefore, it could be assumed that as the reaction temperature increased, more N and S elements in the hazardous waste participated in the pyrolysis reaction to form H 2 S, NH 3 and HCN. The result led to a reduction in the proportion of SO x and NO x in the mixture. In addition, as can be seen in Fig. 5, the yield of various reducing gases tended to be stable, the contents of SO x and NO x were lower, and the emission of pollutants reached a stable level when the operating temperature was above 800 °C. Therefore, the

Effect of the hazardous waste calorific value
When the calorific value of the hazardous waste was low, there was a risk that the counter-flow rotary kiln incineration system could not reach 800 °C. The high level of gaseous pollutants generated by the operation of the system indicated that the hazardous could not be incinerated in the counter-flow rotary kiln. Hence, the actual operation of the counter-flow rotary kiln incineration system was simulated using a mixed pyrolysis model with hazardous waste of different calorific values. Figure 6 shows the calorific values (residual heat) of the reducing gases in the gas mixture obtained from the simulated operation of the mixed pyrolysis model for hazardous waste with calorific values in the range of 1000 kcal/kg to 5000 kcal/kg. Residual heat is the amount of heat that can be released by complete combustion of the combustible gas in the exit gas of the rotary kiln. The simulation results obtained from the mixed pyrolysis model showed that the main flammable gases in the flue gas were CO, H 2 and CH 4 . Hydrocarbon gases with C > 2 were less than five per cent of the total flammable gas. Thus, only H 2 , CO and CH 4 were considered when calculating the residual heat of the mixture. The calculation formula was as follows: where Q is the residual heat, V H 2 , V CO and V CH 4 are the gas volumes, and q H 2 , q CO and q H 4 are the calorific values.
The total energy consumption, the calorific value of hazardous waste and the total residual heat satisfied the following equation: As can be seen from Fig. 6, when the hazardous waste calorific value was 1500 kcal/kg, the residual heat in the mixture entering the second combustion chamber was 3.65 kcal/kg, which indicated that the counter-flow rotary kiln incineration system could satisfy its own energy balance at an operating temperature of 850 °C. The result indicated that the calorific value of 1500 kcal/kg was the critical value at 850 °C.

Effect of primary air volume and primary air speed
Most of the existing rotary kiln incineration systems use the negative pressure suction air supply method. In the actual operation, the primary air volume and velocity are difficult to control, making the rotary kiln incineration system less stable. Therefore, the CFD model was used to adjust the primary air supply method to active air supply, and the primary air volume and speed were adjusted to investigate the effect of primary air volume and speed on the counter-flow rotary kiln incineration system. The cross-sectional area of the air inlet was constant. When the primary air velocity was adjusted to 1.2 m/s, 1.5 m/s and 2.0 m/s, the corresponding air flows were 4560 Nm 3 /s, 5700 Nm 3 /s and 7600 Nm 3 /s, respectively. Figure 7 shows the temperature and CO concentration distributions inside the rotary kiln for the three air velocities mentioned above, respectively. H 2 , CO and CH 4 were all mainly produced in the pyrolysis zone of the rotary kiln, with rotary kiln pyrolysis occurring mainly in the middle section of the rotary kiln. Therefore, the effect of primary air volume and primary air velocity on the distribution of H 2 and CH 4 concentrations was similar to CO. Therefore, their concentration distribution can be referred to that of CO.
As can be seen from Fig. 7a, the variation in air speed had a significant effect on the overall temperature distribution in the kiln. The high temperature region of the rotary kiln gradually moved towards the rotary kiln inlet, when the primary air speed gradual increased. In addition, with the gradual increase in primary air volume, the kiln tail temperature gradually decreased due to the increased heat absorption by cold air. When the temperature of the rotary kiln was too low, the ash was not burned completely and, thus, the burning effect of the rotary kiln was reduced. Therefore, it was not reasonable to supply the rotary kiln with air at a speed of 2.0 m/s. As can be seen from Fig. 7b, the high concentration region of CO showed Residual heat = Calorific value of hazardous waste − Energy consumption Fig. 6 Residual heat in the gas mixture the same characteristics as the high temperature region, both moving to the left of the rotary kiln as the primary air velocity increased. In addition, the increase in primary air volume was beneficial for the incineration of hazardous waste, resulting in a gradual reduction in the CO concentration in the rotary kiln. Therefore, hazardous waste was burned better at a wind speed of 1.5 m/s compared to 1.2 m/s. In general, 1.5 m/s was the optimum air supply speed for a counter-flow rotary kiln incineration system, which corresponded to a primary air flow of 5700 Nm 3 /s.

Simulation results for the second combustion chamber
The simulation results of the rotary kiln showed that the combustible gases entering the second combustion chamber were mainly H 2 , CO and CH 4 . In addition, CH 4 was used to supplement the combustion in order to achieve an exhaust gas temperature of more than 1150 °C. As a result, there was no significant production of SO x and NO x in the secondary combustion chamber. Therefore, the distribution of temperature and CO concentration in the second combustion chamber were given. Figure 8 shows the temperature distribution in the second combustion chamber. The temperature at the entrance to the second combustion chamber was low, which was related to the lower flue gas temperature at the exit of the rotary kiln and the addition of supplementary fuel. However, the combustion of combustible gases in the flue gas and the make-up fuel further rose the temperature in the secondary combustion chamber, resulting in a greater temperature at the outlet of the secondary combustion chamber compared to the industry emission standard of 1100 °C. Figure 9 shows the distribution of CO concentrations in the second combustion chamber. The concentration of CO in the secondary combustion chamber gradually decreased along the outlet direction. The rotary kiln outlet gas contained a certain amount of CO, which resulted in the highest CO concentration at the inlet of the second combustion chamber. Under the conditions of adding supplementary fuel to fuel the combustion, the CO in the second combustion chamber gradually decreased in the direction of the flue gas flow and was reduced to a minimum when the flue gas left the second combustion chamber.

Conclusions
The software Aspen Plus and Fluent were used to establish a mixed pyrolysis model and CFD model of the counterflow rotary kiln incineration system, respectively. The effects of the operating temperature, hazardous waste calorific value, primary air velocity and primary air volume on the system operation effect were investigated. The main conclusions obtained are as follows.
(1) As the operating temperature increased, the content of SO x and NO x gradually decreased, and the content of these two gas products tended to remain constant when the operating temperature of the counter-flow rotary kiln incineration system was greater than 800 °C. The result indicated that the incineration system should meet the condition of an operating temperature greater than 800 °C. (2) At an operating temperature of 850 °C, the incineration system was unable to meet its own heat balance through the combustion of hazardous waste when the calorific value of the hazardous waste was below 1500 kcal/kg. Therefore, when using this incineration system for the incineration of hazardous waste materials, the hazardous waste should be simply screened according to the different calorific values. (3) The active air supply method helped to stabilize the operating conditions of the incineration system, but the primary air speed and the size of the primary air volume affected the actual operating effect of the incineration system. When the primary air volume and speed were too low, the mixed pyrolysis of the hazardous waste was not complete and the thermal reduction rate of the hazardous waste was not up to standard. When the primary air volume and velocity were too high, the temperature field in the kiln was not uniform and the temperature at the kiln head was too high. The simulation showed that the incineration system operated best at a wind speed of 1.5 m/s, which corresponded to a primary air flow of 5700 Nm 3 /s. (4) The combustible gases generated in the rotary kiln were completely combusted in the secondary combustion chamber by means of CH 4 make-up combustion. The temperature in the second combustion chamber could be made greater than 1100 °C by supplementary combustion, which met the industry requirements for controlling pollutant emissions.