Acta Polytechnica https://doi.org/10.14311/AP.2021.61.0199 Acta Polytechnica 61(1):199–218, 2021 © 2021 The Author(s). Licensed under a CC-BY 4.0 licence Published by the Czech Technical University in Prague NEW METHOD EVALUATION OF DETAIL MATERIAL AND HEAT FLOWS FOR SINGLE STRING CEMENT CLINKER PLANT Prihadi Setyo Darmantoa, ∗, I Made Astinaa, Alfian Kusuma Wardhanaa, Alfi Amaliab, Arief Syahlanb a Institut Teknologi Bandung, Faculty of Mechanical and Aerospace Engineering, Jalan Ganesha 10, 40132 Bandung Indonesia b Indonesia Cement and Concrete Institute, Jalan Ciangsana Raya, 16968 Bogor, Indonesia ∗ corresponding author: prihadisetyo@gmail.com Abstract. Material flow in each main equipment of a cement clinker plant, which is very useful for controlling the process, is impossible to be measured during operation due to very high temperatures. This paper intends to overcome the difficulties associated with the measurement of these material flow values. This study presents a new method of calculating material flow (gas and solid) in each main equipment of a single string conventional suspension preheater type of a cement clinker plant. Using the proposed method, mass flow rate at a clinker cooler, kiln, suspension preheater (SP) and even each cyclone separator can be calculated with a heat conservation error less than 1 %. With the application of the least square method for solving the overdetermined system of mass and heat conservation equations obtained in the cyclones of SP, the flow of gas and solid materials entering and exiting each cyclone that cannot be measured directly in the operating plant can be approached. Based on the operation temperature data of gas and solid flows monitored in the control room of an Indonesian cement plant as a case study, the mass flow rate of gas and solid entering and exiting as well as separation efficiency of each cyclone can be calculated. The results show that the separation efficiencies of cyclones 1, 2, 3 and 4 are 95 %, 91.89 %, 84.09 % and 79.51 % respectively. Finally, this study will be very useful by providing data that are impossible to gather by a direct measurement in an operating plant, due to a very high process temperature constraint, for operational control needs, new equipment design, process simulation using computational fluid dynamics (CFD) software and even modification of existing equipment. The proposed method can be applied to all types of modern cement clinker plant configurations, either with or without a calciner including the double strings. Keywords: Mass & heat conservations, suspension preheater, cyclone, separation efficiency. 1. Introduction Due to the fact that processes in a cement manufacturing plant are highly exothermic, it is necessary to conduct a precise heat consumption analysis in an effort to optimize heat conservation and efficiency [1]. Understanding the actual heat consumption and its conservation efforts have been main issues in cement industries [2]. Therefore, a heat audit at a cement plant is conducted to identify opportunities for decreasing heat consumption, increasing the productivity and improving the production has been carried out [3]. The heat consumption assessment of Portland cement production in a Thailand’s cement plant was also reported and found 3.29 GJ per ton of cement [4]. Based on the audit results, there are several methods to improve heat conservation in cement plants: improving the quality of kiln feed in the kiln, modifying the preheater cyclone, repairing the control system, utilizing waste heat and improving the combustion process at the kiln [5]. The realization of the utilization of waste heat discharged into the environment through the preheater and cooler to 5.26 MW of electrical heat in a cement plant in Ethiopia has been reported [6]. Previous studies using heat and exergy analyses in kilns [7], preheater and calciner [8], as well as a thermodynamic analysis of processes in a raw mill [9] and a heat transfer analysis of preheater hot gas [10] found that there are many opportunities left out unexplored for heat conservation. One of the notable efforts in heat conservation is by using renewable materials as an alternative of using fossil fuel, such as: industrial waste [11], biomass (e.g. rice husk and bamboo) [12, 13] and tyres [14]. No less important is the use of alternative materials as a clinker substitution to reduce the production cost and maintain their environment sustainability (e.g. steelmaking slag [15], fly ash [16], pozzolanic material [17], limestone and micro-silica [18]). Modifying processes and equipment was done in an effort to increase the heat efficiency. A top cyclone preheater modification was conducted to increase its separation efficiency and to reduce the return dust, which can contaminate the fine coal in the coal mill was explored in [19, 20]. A modification of the process using co-fuel combustion and modelling of calcination process in the calciner was explored in [21, 22]. A modification 199 https://doi.org/10.14311/AP.2021.61.0199 https://creativecommons.org/licenses/by/4.0/ https://www.cvut.cz/en P. S. Darmanto, I M. Astina, A. K. Wardhana et al. Acta Polytechnica of the tertiary duct of kiln hood and modernization of the kiln using pyroclone calciner was explored in [23, 24]. A development of the design concept as well as the evolution of the main burner was explored in [25]. These studies used computational fluid dynamics (CFD) to understand the mechanism of any equipment modification. CFD simulations required precise and accurate data such as mass flowrate of the gas and material entering into the equipment [26, 27], but it was rarely mentioned where the data came from. As simulation results, heat consumption, life cycle inventory, and waste generated from the cement manufacturing plant can be carried out profoundly [28, 29]. However, obtaining the precise input data from a currently active cement manufacturing plan remain as prevailing issue. The issues are due to limited measurement instruments and methods and very high temperature conditions, which sometimes make the measurements impossible to do if the factory is running. The input data measurement for the simulations can be obtained from a direct possible measurement in a cement plant. These data consist of operational parameters: feed rate, clinker production, fuel mass rate, air consumption rate for transporting fuel, clinker cooling air rate and equipment surface temperature. Meanwhile, the material mass conservation between two interconnected devices cannot be measured directly when the plant is in operation. The examples are between the cyclones in the suspension preheater (SP) during the preheating process of materials; between the cyclone and the calciner during a calcination process; and the between kiln entrance and the clinker cooler. These mass conservation calculations such as the input and output of a material flow rate in each cyclone stage and calciner are necessary for the design and modification of a cyclone and calciner. The difficulties in achieving direct measurements in the cement manufacturing equipment are still critical due to technical limitations such as high temperature, dusty environment and others [30]. Furthermore, the separation efficiency of each cyclone is closely related to the mass flow rate and heat in the whole SP. For example, to obtain cyclone separation efficiency values, the design results are normally used in the calculation of mass and heat conservations [31]. However, this certainly does not match the conditions of the real equipment when the factory is running. Thus, parametrization based on a real operating condition value is necessary to obtain a realistic approach for the preheater system design modification. A single string preheater is currently still used in a cement manufacturing plant. This technology consumes high thermal heat while having low production capacity. In order to increase the efficiency of the fuel used and increase the production, a plant modification using a newly designed calciner shall be conducted [32]. Thus, the aim of this research is to propose a new method for estimating the value of cyclone efficiency by utilizing detail mass and heat conservation equations on each cyclone in the preheater system, without making direct mass flow rate measurements on the related equipment. To simplify the results of modelling and calculation, the equation of mass and heat conservations at cyclone preheaters are formed as a matrix and solved using the least squares method. The equation of this study is an overdetermined equation system that could be solved by using an approximation method. From the estimated cyclone separation efficiency value, the mass flow rate in each component of the cyclone can be approximated. The literature survey indicates that studies of detail materials and heat flow evaluation of the cyclone separator in the SP are limited in number and scope. Generally, only limited to kiln [33–35], clinker cooler [36] and all the main equipment of a cement plant as a whole [37]. The additional objective of this study is to contribute providing detailed data required for an equipment design, modification, control of plant operation, and detail heat audit. 2. Materials and method 2.1. Materials The SP configuration of cement plants depends on the plant manufacturers. According to the classification proposed by Schmidt [38], the configuration of the studied dry process cement plant suspension preheater is a conventional suspension preheater (SP) type as shown in Fig. 1, while a schematic diagram of the whole cement plant main equipment is presented in Fig. 2. SP consists of four cyclone separators and its riser duct is installed in series. Fig. 1 shows raw-meal being inserted into SP through the riser duct of the top cyclone, while the material coming out of the preheater is transported towards the kiln for further processing to become clinker. The fuel used is coal, and is introduced into the kiln using the main burner. Table 1 shows the chemical composition of kiln feed and coal that are normally used in an Indonesian cement plant, including in this study. 2.2. Mass and heat conservation equations From Fig. 2, with the assumption that the plant condition is steady and the chemical composition of the kiln feed and fuel (coal) is as presented in Table 1 and the ash contained in fuel forms clinker, the material flow conservation of the plant can be written as Eq. 1. The production of the clinker mcli can be evaluated by Eqs. (2) to (8). mkf + mcoal + mtr−air + mcoal−air = mcli + mH2O−kf + musnep_kf _1 + mgas−kf (1) 200 vol. 61 no. 1/2021 New method evaluation of detail material and heat flows. . . Kiln feed composition % of mass Coal composition % of mass CaO 43.29 H2 4.44 MgO 1.43 C 73.40 SiO2 13.82 N2 1.30 Al2O3 3.37 O2 11.10 Fe2O3 1.77 S 0.50 H2O 0.32 H2O 4.70 LOI (CO2, Na2O, K2O & SO3) 36.00 Ash + dust 13.00 Table 1. Chemical composition of kiln-feed and coal. Figure 1. Suspension preheater (SP). Figure 2. Material flow in cement clinker’s plant main equipment. 201 P. S. Darmanto, I M. Astina, A. K. Wardhana et al. Acta Polytechnica mcli = msep−kf c + mash (2) msep−kf c = (1 −LOIkf ) ·msep−kf (3) mgas−kf = (LOIkf −H2Okf ) ·msep−kf (4) msep−kf = ηC1.A&1.B ·mkf (5) munsep_kf _1 = (1 −ηC1.A&1.B ) ·mkf (6) mash = ashcoal ·mcoal (7) mH2O−kf = H2Okf ·msep−kf (8) If it is assumed as complete combustion, the mass flow rate of the hot gas (mhg) is equal to the flue gas resulted from coal combustion process. The value of mhg can be estimated by Eq. (9): mhg = (1 −ashcoal) ·mcoal + mtr−air + mcoal−air (9) Based on the Fig. 2, the heat conservation equation of the plant can be written in Eq. (10):∑ Enin = ∑ Enout (10) where ∑ Enin is the sum of the heat flow entering the plant equipment, while ∑ Enout is the sum of the heat flow exiting the plant. The heat flow or heat entering the plant consists of: (1.) Heat flow of kiln feed entering the SP (Eq. 11): Enkf = mkf ·hkf (Tkf ) (11) (2.) Heat flow of coal entering the kiln (Eq. 12): Encoal = mcoal ·hcoal(Tcoal) (12) (3.) Heat flow of air entering the plant through clinker cooler and main burner (Eq. 13): Enair = (mtr−air + mcool−air ) ·hair (Tair ) (13) (4.) Heat flow resulted by coal combustion process (Eq.14): Encoal−comb = mcoal ·NHVcoal (14) Meanwhile, the flow of heat exiting plant’s main equipment, with the assumption that the temperature of dust and gas exiting from SP is equal Thg1, contains of: (1.) Heat flow of clinker product at temperature Tcli (Eq. 15): Encli = mcli ·hcli(Tcli) (15) (2.) Heat flow of coal combustion process gas exiting from the top of the cyclone of SP (Eq. 16): Enhg = mhg ·hhg (Thg1) (16) (3.) Heat flow of vapour created by the evaporation of water content in kiln feed exits from SP (Eq. 17): Envapor−kf = mH2O−kf ·hvapor (Thg1) (17) It is noted that the mass of vapour created by the coal combustion process is included in the mass flow of flue gas. (4.) Heat flow of unseparated kiln feed exiting from the top of the cyclone of SP (Eq. 18): Enunsep−kf −1 = munsep−kf −1 ·hkf (Thg1) (18) (5.) Heat flow of kiln feed gas resulting from the calcination process, mostly CO2, exits from SP (Eq. 19): Engas−kf = mgas−kf ·hCO2 (Thg1) (19) 202 vol. 61 no. 1/2021 New method evaluation of detail material and heat flows. . . Material A B C Raw meal 0.206 101 -37 Clinker 0.186 54 0 Coal 0.262 390 0 Air 0.237 23 0 CO2 0.196 118 -43 O2 0.218 30 0 Vapor (H2O) 0.443 39 28 Table 2. Values of A, B, and C of Eq. (24) [41]. (6.) Heat of the clinker formation where the value can be calculated by Bogue’s equation multiplied by the rate of the clinker production [39] (Eq. 20): Enf orm = mcli{7.646(CaCO3)kf + 6.48(MgCO3)kf + 4.11(Al2O3)kf− − 5.176(SiO2)kf − 0.59(Fe2O3)kf} (20) where (CaCO3)kf , (MgCO3)kf , (Al2O3)kf , (SiO2)kf , and (Fe2O3)kf are mass based percentage of each substance contained in the kiln feed. It should be noted that for each kg of clinker produced, the heat of the clinker formation consists of calcination heat (Encalc), which is endothermic and the heat of clinkerization or sintering (Enclink), which is exothermic. Hence, the heat of the clinker formation is the difference between the calcination and the sintering heat. Calcination is the decomposition reaction of CaCO3 to form CaO and CO2 while sintering is the reaction process of the formation of oxides contained in the clinker. Sintering occurs in the temperature range of 1250 to 1450 °C, thus it only occurs in the kiln. The heat of this calcination reaction which can occur in the SP and kiln is considered constant at 425 kcal/kg of CaCO3 [40]. By determining the calcination heat per kg of CaCO3, the calcination heat per kg of clinker produced can be calculated. Finally, we can summarize the calculation in Eq. (21). Enclink = Encalc −Enf orm (21) (7.) Evaporation heat of the water content in the kiln feed and coal exits from SP (Eq. 22): Enevap = (mH2O−kf + mH2O−coal)hf g (22) where hf g is the enthalpy of the water evaporation process. (8.) Heat loss by radiation and convection (Qloss) through the overall surface area of the main equipment Atot that can be calculated by the proposed formula [41] (Eq. 23): Qloss = Atot(4 · 10−8(T 4surf −T 4 amb) + 80.33 Å Tsurf −Tamb 2 ã−0.724 (Tsurf −Tamb)1.333) (23) Enthalpy of the substances can be calculated, if their temperatures are known from a measurement, using Eq. (24) h(T) = ∫ T2 T1 CpdT = A + BT 2 · 10−6 + CT 3 · 10−9 (24) where Cp is the specific heat of the substance and T1 is the reference temperature. The value of constants A, B, and C for limited substances and gas are shown in Table 2. 2.3. Methodology In general, a cement plant is equipped with measuring devices and instrumentations to monitor operating parameters such as temperatures and draft pressures. A weighing feeder is applied for measuring the mass flow rate of the kiln feed and coal. An oxygen content monitoring system is put in the kiln and SP gas outlets. An obstruction meters for measuring the air flow are generally available at the control system of the plant. In addition, the flow of dust returning from the unseparated kiln feed leaving for the SP is also measured. This allows an estimation of the separation efficiency of the top cyclone by using Equation (6). The first methodology of the study is the evaluation of mass conservation of the whole plant based on the measurement results of these instrumentation devices using Equations (1) to (9). The second methodology is to solve the mass and heat rates in each main equipment in more detail, (e.g. clinker cooler, kiln, SP and each cyclone separator). Additional 203 P. S. Darmanto, I M. Astina, A. K. Wardhana et al. Acta Polytechnica Figure 3. Materials flow in planetary cooler. mass and heat conservation equations are required for each equipment. These equations are interdependent and can be formed into linear algebraic equations and solved simultaneously. However, especially for cyclones, because the number of equations is greater than the number of parameters to be searched for, the approximation solution taken is by using the least square method with the minimum square of a residual error as the criteria. 2.4. Mass and heat conservation equations of planetary clinker cooler For the planetary clinker cooler, a schematic diagram of mass flow entering and leaving the planetary cooler is presented in Fig. 3. The mass conservation equation can be written as Eq. (25): mcli_k + mcool−air = mcli + mclid_to_k + mcool_air_k (25) where mcli, mcli_k , mcool−air, mcool_air_k and mclid_to_k are the mass flow rates of the clinker product, clinker entering the cooler, cooling air entering the cooler, cooling air entering the kiln (combustion air) and clinker dust returning to kiln, respectively. With the assumption that there is no leakage in this cooler, the mass flow rate of the cooling air is equal to the combustion air minus the transporting air. The rate of the clinker mass from leaving the kiln and entering the cooler must be equal to the sum of the clinker production and the returning clinker dust back to the kiln. This clinker dust is carried by the combustion air back from the cooler to the kiln. The percentage of clinker dust returning to the clinker production will be assumed in the calculation of the mass conservation in the cooler. The heat conservation calculation error in the cooler and kiln is not more than 1 %. With limits determined like this and with the measurement data from the plant, the calculation of the mass and heat conservations in the cooler can be carried out. The planetary cooler heat conservation equation is written in Eq. (26): Encli_k + Encool−air = Encli + Enclid_to_k + Encool_air_k + Qloss (26) The flow of heat entering the planetary cooler consists of: (1.) Heat flow of clinker from kiln at a temperature Tcli_k: Encli_k = mcli_k ·hcli(Tcli_k) (27) (2.) Heat of cooling air: Encool_air = mcool_air ·hair (Tair ) (28) Meanwhile, the flow of heat leaving for the cooler, assuming that the temperature of the dust returning to kiln and the gas exiting from the cooler is equal to Thg_c consists of: (1.) Heat flow of clinker product Encli at a temperature Tcli as mentioned in Eq. (15). (2.) Heat of combustion air entering the kiln Ensec_air_c at temperture Thg_c is: Ensec_air_c = mcool_air_k ·hhg_c(Thg_c) (29) (3.) Heat flow of clinker dust returning to kiln Enclid_k at temperture Thg_c: Enclid_k = mclid_to_k ·hcli(Thg_c) (30) Heat loss by radiation and convection Qloss_c through the surface area of the clinker cooler (Acooler) can be calculated by Eq. (23). In Eqs. (25) and (26), the value of the mass rate of the clinker dust returning to the kiln (mclid_to_k) and combustion air temperature cannot be measured directly at the plant. However, the mass and heat conservation in the kiln can be evaluated with the minimum acquisition limit between the heat entering and exiting the kiln. Using this condition, the Eqs. (27) to (30) can be solved by completing the heat conservation of the cooler. 204 vol. 61 no. 1/2021 New method evaluation of detail material and heat flows. . . Figure 4. Materials flow in the kiln. 2.5. Mass and heat conservation equations of the kiln For kiln system, the mass flow in and from the kiln is described in Fig. 4. The equation of mass conservation is written in Eq. (31): msep−kf −4 + mclid_to_k + mcoal + mcomb_air = mcli_k + mclid_to_SP + mhg_to_SP (31) where msep−kf −4 = mass flow rate of kiln feed from SP. mcoal = mass flow rate of coal. mcomb_air = mass flow rate of combustion air (cooling air mcool-air and transporting air mtr-air). mclid_to_SP = mass flow of clinker dust entering the SP. mhg_to_SP = mass flow rate hot gas entering the SP. The hot gas entering the SP consists of flue gas, as mentioned in Eq. (9), and CO2, including other gasses from the calcination and burning process. The portion of the calcinated kiln feed can be evaluated using equation (4) multiplied by the percentage of the kiln feed calcinated in the kiln (%Calc_k). Thus, the value of mhg_to_SP can be described as Eq. (32), mhg_to_SP = mhg + (%Calc_k ·mgas_kf ) (32) The other notations were already mentioned in the mass conservation of the cooler. While the heat conservation equation of the kiln is as follows (Eq. 33), Enkf _to_k + Encoal + Encoal−comb + Ensec_air_c + Entr_air + Enclid_k = Encli_k + Enclid_to_SP + + Enf orm_k + Enhg_to_SP + Qloss_k (33) The flow of heat entering the kiln consists of: (1.) Heat flow of kiln feed from the SP at a temperature Tin_k (Eq. 34): Enkf _to_k = msep−kf −4 ·hkf (Tin_k) (34) (2.) Heat of coal entering the kiln as defined in Equation (12). (3.) Heat from coal combustion process (Eq. 14). (4.) Heat flow of cooling and transporting air (Eq. 13). (5.) Heat of clinker dust returning from cooler as written in Eq. (30). While the heat flows leaving for the kiln consist of: (1.) Heat flow of clinker product Encli_k to the cooler (Eq. 27). (2.) Heat flow of clinker dust returning to the SP Enclid_to_SP in which its temperature is assumed equal to the temperature of the hot gas leaving for kiln Tin_k: Enclid_to_SP = mclid_to_SP ·hcli(Tin_k) (35) (3.) Net heat of calcination and clinkerization process in the kiln Enf orm_k that can be assumed equal to the sum of calcination heat and clinkerization heat. Hence, the value of Enf orm_k can be written as follows, Eq. (36). Enf orm_k = Encalc_k + Enclink (36) where Encalc_k = %Calc_k ·Encalc. Here, %Calc_k can be estimated indirectly when the evaluation of the heat conservation in the kiln is carried out. The value of %Calc_k also depends on the percentage of the kiln feed that was calcined in the SP (%Calc_SP). Therefore, %Calc_k = 1 − %Calc_SP. In this study, we assume that the calcination process in the SP occurs at the cyclone 3 (%Calc_SP_C3) and the lowest cyclone (%Calc_SP_C4) and their riser ducts. Thus %Calc_SP = %Calc_SP_C3 + %Calc_SP_C4. 205 P. S. Darmanto, I M. Astina, A. K. Wardhana et al. Acta Polytechnica (4.) The temperature of the heat flow of hot gas from the kiln to the SP is Tin_k. The total mass of this hot gas is the sum of the mass flow rate of the combustion gas mhg (Eq. 9) and CO2 produced by the calcination process of the kiln feed in the kiln equal to (%Calc_k ·mgas−kf ). Therefore, the heat flow of the hot gas can be calculated as follows (Eq. 37), Enhg_to_SP = mhg ·hhg (Tink ) + (%Calc_k ·mgas−kf ) ·hCO2(Tink ) (37) The heat loss by radiation and convection Qloss_k, and the heat loss through the surface area of the kiln (Akiln) can be calculated by Equation (23). 2.6. Detail mass and heat conservation equations of SP system and each cyclone separator The flow of solid materials and gas of the SP is presented in Fig. 1 and its related mass conservation equation can be written as follows (Equation 38), mkf + mclid_to_SP + mhg_to_SP = msep−kf −4 + mH2O−kf + munsep_kf _1 + mhg + mgas−kf (38) It should be noted that the values of mkf _to_k and mclid_to_SP are known by solving the kiln mass conservation equation. Heat flow into the SP consists of: (1.) Heat of kiln feed entering the SP, Enkf . (2.) Heat of clinker dust returning from the kiln, Enclid_to_SP , which is represented by Eq. (35). (3.) Heat flow of hot gas from kiln, which was presented by Eq. (37) and consisted of the sum of the heat flow rate of the combustion gas mhg and LOI produced by the calcination process of the kiln feed. While the heat flows leaving the SP consist of: (1.) Heat flow of unseparated kiln feed by top cyclones (Eq. 18). (2.) Heat flow of coal combustion hot gas exiting from the top cyclone as formulated by Eq. (16). (3.) Heat of evaporated water content in the kiln feed as written by Eq. (17). (4.) Heat of the gas product of the kiln feed calcination and burning process, mainly CO2. Its value is approximated by Eq. (19). (5.) Heat flow of the kiln feed entering the kiln at a temperature Tink as presented in Equation (34). (6.) Heat of the clinker formation in SP Enf orm_SP that can be assumed equal to the heat of the clinker formation Enf orm (Eq. 20) minus the heat formation absorbed in the kiln Enf orm_k: Enf orm_SP = Enf orm −Enf orm_k (39) It should be noted that the heat formation in the SP only consists of the calcination heat. Its value can be estimated based on the percentage of the calcined kiln feed in the cyclones 3 and 4. We can write this term as Enf orm_SP = (%Calc_SP_C3 + %Calc_SP_C4) ·Encalc. (7.) Heat loss by radiation and convection Qloss_SP through the SP’s surface area (ASP ) that can be calculated by Eq. (23). Finally, the heat conservation equation for the SP is written in Eq. (40): Enkf + Enclid_to_SP + Enhg_to_SP = Enunsep−kf −1 + Enclid_to_SP + Enkf _to_k+ + Enf orm_SP + Engas−kf + Envapor−kf + Qloss_SP (40) For understanding the flows of mass and heat at each cyclone separator that is very useful for the plant operation and design as well as equipment modification, its formulation will be derived one by one. The design of top cyclone 1A and 1B consist of 2 cyclones of a diameter smaller than their lower cyclones. The reason of this design is to achieve a higher separation efficiency. The schematics of these top cyclones are presented in Fig. 5a, while the lower cyclones are depicted in Fig. 5b. For the top cyclones, assuming no flow leakage, the incoming mass flow consists of: (1.) kiln feed, mkf at temperature Tkf . (2.) flow of unseparated kiln feed from the cyclone 2, munsep_kf _2 where its temperature is Thg2. (3.) gas flow from the cyclone below at a temperature Thg2, which consists of combustion gas (mhg )in and kiln feed gas (CO2 and others), (mgas−kf )in. 206 vol. 61 no. 1/2021 New method evaluation of detail material and heat flows. . . Figure 5. Mass flow of top (left) and other (right) cyclones. While the outflow of mass consists of: (1.) separated kiln feed to the cyclone 2, msep−kf −1. (2.) unseparated kiln feed at a temperature Tkf 1, (munsep_kf _1). (3.) combustion gas, kiln feed gas and water vapour (mH2O−kf ) at temperature Tkf 1. The mass conservation equation for the top cyclone can be written as Eq. (41). mkf + munsep_kf _2 + (mhg )in + (mgas−kf )in = msep−kf −1 + mH2O−kf + munsep_kf _1+ + (mhg )out + (mgas−kf )out (41) For a general SP operation, the gas temperature that heats the kiln feed is slightly higher than the kiln feed temperature when exiting the cyclone (Tkf 1). These temperatures are measured and can be read from the control room. The rate of the kiln feed and unseparated dust filtered by baghouse can also be measured. Assuming that the dust from cyclone 2 will be filtered back by the top cyclone. Thus, the material that has been filtered by the top cyclone is a part of the circulating material in the the SP. Furthermore, the ratio of (mkf −munsep_kf _1) and mkf can be used for the calculation of the value of its separation efficiency, ηC1.A&1.B. In this case, (mgas−kf )in = (mgas−kf )out = mgas−kf and (mhg )in = (mhg )out = mhg and it will be valid for the following derived equations. Based on this operating conditions, the top cyclone separation efficiency can be calculated by Eq. (5) and its heat conservation can be written as follows (Eq. 42), Enkf + munsepkf2 ·hkf (Thg2) + mhg·hhg (Thg2) + mgas−kf·hCO2 (Thg2) = Enhg + Engas−kf + + Envapor−kf + Enunsep−kf −1 + msep−kf −1·hkf (Tkf 1) + Qloss_C1 (42) where Qloss_C1 is the heat loss by radiation and convection through the top cyclone surface area (AC1) that can be calculated by Eq. (23). Using the same notation of materials and gas flows as denoted in Fig. 5right, the mass conservation equation of the cyclone 2 can be written as Eq. (43). msep−kf −1 + munsep_kf _3 + (mhg )in + (mgas−kf )in = msep−kf −2 + munsep_kf _2+ + (mhg )out + (mgas−kf )out (43) While its heat conservation equation is (Eq. 44): msep−kf −1 ·hkf (Tkf 1) + munsep_kf _3 ·hkf (Thg3) + mhg ·hhg (Thg3) + mgas−kf ·hCO2 (Thg3) = = msep−kf −2 ·hkf (Tkf 2) + munsep_kf _2 ·hkf (Thg2) + mhg ·hhg (Thg2) + mgas−kf ·hCO2 (Thg2) + + Qloss_C2 (44) For the cyclone 3, the temperature of the exiting gas and solid materials are Thg3 and Tkf 3, respectively. Meanwhile, the temperature of the gas entering this cyclone from the lowest cyclone is Thg4. The calcination reaction starts at a temperature of around 600 °C [42–44] and the kiln feed temperature from cyclone 3 (Tkf 3) is 207 P. S. Darmanto, I M. Astina, A. K. Wardhana et al. Acta Polytechnica higher than 600 °C. Hence, a part of the kiln feed is calcined. By assuming that the percentage of calcined kiln feed in cyclone 3 is %Calc_SP_C3, then the CO2 formed from this reaction and the associated reaction heat must be taken into account in the mass and heat conservation equation. Hence, the mass conservation equation for cyclone 3 can be written as follows (Eq. 45), msep−kf −2 + munsep_kf _4 + (mhg )in + (mgas−kf ) in = msep−kf −3 + munsep_kf _3+ + (mhg )out + (mgas−kf )out (45) where (mgas−kf ) out− (mgas−kf ) in = %Calc_SP_C3 ·mgas−kf = (msep−kf −2 + munsep_kf _4)− − (msep−kf −3 + munsep_kf _3) (46) and its heat conservation equation is: msep−kf −2 ·hkf (Tkf 2) + munsep_kf _4 ·hkf (Thg4) + mhgmhg ·hhg (Thg4) + + mgas−kf ·hCO2 (Thg4) = msep−kf −3 ·hkf (Tkf 3) + munsep_kf _3 ·hkf (Thg3) + + Enf orm_SP _C3 + mhg ·hhg (Thg3) + mgas−kf ·hCO2 (Thg3) + Qloss_C3 (47) It should be noted that the kiln feed calcination process in cyclone 3 (msep−kf c−3) continues in the lowest cyclone due to the higher temperature of the hot gas coming from the kiln. Therefore, a part of the kiln feed mass transforms to CO2 in this lowest cyclone. The mass conservation equation of the lowest cyclone is finally written as (Eq. 48), msep−kf −3 + (mhg )in + mclid_to_SP + (%Calc_k ·mgas−kf ) = = msep−kf −4 + munsep_kf _4 + (mhg )out + (mgas−kf ) out (48) The gas produced by the calcination process in the lowest cyclone is equal to (Eq. 49), (mgas−kf ) out− (%Calc_kmgas−kf ) = %Calc_SP_C4mgas−kf = = (msep−kf −3 + mclid_to_SP ) − (msep−kf −4 + munsep_kf _4) (49) The related heat conservation equation can be given as Equation (50), msep−kf −3 ·hkf (Tkf 3) + Enclid_to_SP + Enhg_to_SP = munsep_kf _4 ·hkf (Thg4) + + mhg ·hhg (Thg4) + mgas−kf ·hCO2 (Thg4) + Enkf _to_k + Enf orm_SP _C4 + Qloss_C4 (50) The mass and heat flows in cyclones 2, 3 and 4 (lowest cyclone) cannot be calculated sequentially. This is due to the unknown value of the mass of the kiln feed separated by cyclones 2 and 3 and unseparated kiln feed from cyclones 3 and 4. These four unknown parameters are interdependent and must be solved simultaneously. To evaluate all four parameters simultaneously, equations (43) to (45), (47), (48) and (50) are used and written in the form of a matrix (Eq. 51): [A] × [X] = [B] (51) The matrices [A], [X] and [B] are represented by Equations (51), (52) and (53), respectively. Matrix [A] 6 × 4 is a multiplier constant of unknown variables, while the matrix [B] 6 × 1 is a constant resulting from the known constant value of the equation (51). [A] =   1 0 −1 0 hkf (Tkf 2) 0 −hkf (Thg3) 0 −1 1 1 −1 −hkf (Tkf 2) hkf (Tkf 3) hkf (Thg3) −hkf (Thg4) 0 −1 0 1 0 −hkf (Tkf 3) 0 hkf (Thg4)   (52) [X] =   msep−kf −2 msep−kf −3 munsep_kf _3 munsep_kf _4   (53) 208 vol. 61 no. 1/2021 New method evaluation of detail material and heat flows. . . [B] =   msep−kf −1 −munsep_kf _2 msep−kf −1 ·hkf (Tkf 1) + mhg · {hhg (Thg3) −hhg (Thg2)}+ mgas−kf · {hCO2 (Thg3) −hCO2 (Thg2)}− munsep−kf −2 ·hkf (Thg2) −Qloss_C2 0 mhg · {hhg (Thg4) −hhg (Thg3)}+ mgas−kf · {hCO2 (Thg4) −hCO2 (Thg3)}− Enf orm_SP _C3 −Qloss_C3 mclid_to_SP −msep−kf −4− {(1 − %Calck)·mgas−kf} Enclid_to_SP + Enhg_to_SP −mhg ·hhg (Thg4)− mgas−kf ·hCO2 (Thg4) −Enkf _to_k− Enf orm_SP _C4 −Qloss_C4   (54) Matrix [X] 6×1 is the variable which shall be solved. From the equation (51), it is known that the number of unknown variables is less than the number of equations (overdetermined equations). A system that has a number of equations greater than the number of variables usually do not have a solution or cannot be solved directly [45]. The conservation of mass and heat equations in the form of matrix as written in Eqs. (52) to (54) are stand-alone equations so that the rank matrices [A] and [A|B] cannot be reduced from one another because the matrix compiler equation is not a combination of other equations. Thus, the overdetermined equation systems usually have inconsistent variables. However, the solution can still be found using the Least Squares method, which is to find the value to fulfil the equation of results with the minimum square of the residual error [46]. The method used is by least squaring the matrix [A] by multiplying the transpose matrix itself [A]T so that the matrix [A]T [A] has the same number of rows and columns. The matrix [A]T [B] will change the row of the matrix [B] so the matrix [B] has the same line as the number of rows of [A]T as shown in Eq. (55). The variable values in the matrix [X] can be found by multiplying the inverse matrix [AT A]−1 on both sides of the equation as described in Eq. (56), that can be solved using mathematical software. [A]T [A] × [X] = [A]T [B] (55)[ AT A ]−1 [ AT A ] × [X] = [ AT A ]−1 [A]T [B] (56) From equation (56), the values of X ≈ HLS are called as unique least squares approximate solution that comply the Least Square method, and are defined by Eq. (57). XLS = [ AT A ]−1 [ AT B ] (57) Knowing the approximating values of X, the approximate value of the separation efficiency of the cyclones 2, 3, and 4 can be calculated using the following Eqs. (58)-(60) [47, 48]. ηC2 = msep−kf −2 msep−kf −1 + munsep_kf _3 (58) ηC3 = msep−kf −3 msep−kf −2 + munsep_kf _4 − (%Calc_SP_C3 ·mgas−kf ) (59) ηC4 = msep−kf −4 msep−kf −3 + mclid_to_SP − (%Calc_SP_C4·mgas−kf ) (60) The approximate value of the cyclone separation efficiency can be substituted into the mass conservation equation of each cyclone to recalculate the value of X that meets the mass conservation law. The final results of X can be used to evaluate the flow of heat in each cyclone with a minimum error. Hence, it can be accepted from an engineering point of view (i.e. error < 1 %). The flowchart of the calculation steps is shown in Figure 6. All above mentioned equations and methods of calculations have been adopted as base in developing the substantial engineering equation solver (i.e., in the form of executable programs) and used to calculate detail materials and heat flows and conservations. 209 P. S. Darmanto, I M. Astina, A. K. Wardhana et al. Acta Polytechnica No Data input: - Kiln feed and fuel chemical composition and heating value - operational parameters(feeding, temperatures of gas, materials, surface, clinker, fuel flow, primary air flow, return dust flow, cooling air) - first estimating of top cyclone separation efficiency - percentage of return dust from cooler to kiln and from kiln to SP. - Calculation of plant mass & heat balance and klinker product using Eqs. 1-9 - Calculation of SP, kiln and clinker cooler's mass & heat balance using Eqs. 10-40 - Using data of material's flow exiting from C4 and top cyclones, calculation detail mass balance of each cyclones using least square method (Eqs. 41-57) - Calculation of separating efficiency of each cyclone using Eqs. 58-60 - Subtitution to mass & heat balances equations of detail each cyclone using Eqs. 41-50 ΔECyclone < 1% best approximation results report end ΔmCyclone = 0 No Using new values of cyclones separation efficiency slight tuning of gas and/or material's temperatures of the cyclone Yes start Yes Figure 6. Flowchart of the calculation steps. 3. Results and discussion 3.1. Operational data parameters A cement factory plant in Indonesia, Tonasa 2, was used as the object of this study. Daily average data at a normal operation capacity were used for the current analysis. These data were taken from the control room and direct field measurement. The data for input values are classified into several groups (i.e. kiln feed, coal, SP, kiln, and cooler data). The coal consumption is 271 tons per day (TPD) or 3.137 kg/s with 6,100 kcal/kg of coal nett heating value. The rate of the transporting air is 2.08 m3/s. The rate of the inlet kiln feed is 3,195 TPD. The rate of the return dust is 158.8 TPD (1.85 kg/s). The remaining oxygen percentage in the kiln and top cyclone outlets are 2.5 % and 2.7 %, respectively, and can be used to estimate the kiln excess air and SP false air [41]. The ambient temperature is 33 °C. Based on these measurement data, the production of clinker, flue gas and kiln feed gas can be evaluated using Eqs. (2) to (9). Meanwhile, the data for kiln, planetary cooler, and SP are based on measurement results and general assumptions of normal cement plant operations and are given in Tables 3 and 4. Another parameter that can’t be measured directly in the plant is the mass flow rate of the re-circulated dust from the cooler to the kiln and from the kiln to the SP, in this study, their values are assumed as 15 % and 17 % of the clinker production rate, respectively. 210 vol. 61 no. 1/2021 New method evaluation of detail material and heat flows. . . Parameters Unit Kiln Cooler Average surface temperature (°C) 330.7 189.2 Cylinder diameter (m) 4.5 2.3 Length of cylinder (m) 75 16 Number of cylinder unit 1 11 Clinker product temperature (°C) - 190 Entering air temperature (°C) - 33 Estimated return dust from cooler (%) 15 - Table 3. Input data for kiln and planetary cooler. Parameters Unit Cyclone 1 Cyclone 2 Cyclone 3 Cyclone 4 Average surface temperature (°C) 164.1 194.2 198.8 226 Total surface area (m2) 450.7 376.9 376.5 455.3 Hot gas inlet temperature (°C) 560 694 867 1190 Hot gas exit temperature (°C) 360 560 694 867 Separated kiln feed temperature (°C) 345 545 680 837 Estimated return dust from kiln (%) - - - 17 Estimated kiln feed calcinated exit SP (%) - - - 25 Table 4. Input data for each cyclone of the SP. 3.2. Plant mass and heat conservations Using the available input data, the mass and heat conservations for the case study plant can be evaluated using Eqs. (10) to (23). The results of the mass and heat conservations calculation is presented in the form of per unit mass of clinker products (i.e., per kg of clinker produced). With the kiln feed mass fed to the plant of 3,195 TPD and the measured top cyclone separation efficiency of around 95 %, the clinker production is approximately 1,978 TPD (22.9 kg/s). For the whole clinker plant, the overall mass and heat conservations per unit mass of clinker produced are given in Table 5. Generally, the heat consumption at a higher capacity of the cement plant is lower. From the results, it can clearly be seen that the formation of heat per kg of clinker is around 425 kcal. The heat of this clinker formation generally consists of the calcination heat of CaCO3, and clinkerization heat to form clinker-contained oxides. The calcination process is endothermic, while the clinkerization process is an exothermic process. Because both processes occur at a high temperature range (600 °C - 1,450 °C), an additional heat is needed to maintain the temperature of the kiln and Suspension Preheater (SP). This is the reason why the plant heat consumption is far above its clinker formation heat. This excess heat is finally carried away by the flue gas, CO2 and other gasses out of the SP. This hot gas, with a low oxygen content, is used as a heat source for a drying of coal and raw materials. In modern cement factories, the flue gas is even used as for electricity generation. In addition, the excess heat is also carried out by clinker products. The clinker heat is recovered in the cooler to increase the combustion air temperature and reduce the coal requirements. In effort to reduce the production cost and supporting environmental sustainability, part of the fuel used is low rank coal and alternative fuels. However, by shifting to lower rank coal with a high water content and alternative fuels, the specific heat consumption could increase to ±850 kcal/kg of the produced clinker [49]. 3.3. Cooler and kiln mass and heat conservations The calculation results of the mass and heat conservation at the cooler and kiln, using Eqs. (24) to (36), are given in Tables 6 and 7, respectively. The absolute value of the incoming and outgoing difference heat of the cooler is less than 1 %. The temperature of the combustion air entering the kiln is around 960 °C and its recovery heat efficiency is 74.55 %. This combustion air temperature is in the range of the value measured by Alfi et.al [49]. Meanwhile, the percentage of the heat recovery in the cooler is relatively high when compared to the grate cooler, which is widely used in large capacity factories. This is caused by the conditions on the planetary cooler with the entirety of the cooling air being used as a combustion air. Whereas, on the grate cooler, not all of the cooling air is used as combustion air. The consequence is that the clinker temperature coming out of the cooler is generally higher on the planetary cooler than in the case of the grate cooler. The heat loss from the cooler is large enough that it can be further utilized (i.e., for drying coal in Tonasa 2 plant [50]). In Table 7, the mass and heat conservation per kg of clinker produced in the kiln can be seen. The difference between the rate of heat entering and exiting the kiln is insignificant, with an approximate value of 0.46 %. This value is obtained by assuming that 17 % of clinker from 211 P. S. Darmanto, I M. Astina, A. K. Wardhana et al. Acta Polytechnica Parameters Inlet Outletmass (kg) heat (kcal) mass (kg) heat(kcal) Kiln feed in 1.6155 20.541 - - Coal to main burner 0.1370 837.102 - - Total combustion air 1.5850 11.442 - - Clinker production - - 1.0000 37.289 Combustion gas including excess air - - 1.7042 154.620 Kiln feed gas - - 0.5528 46.744 Return dust - - 0.0805 6.887 Clinker heat formation - - - 425.162 Convection & radiation loss - - - 193.285 Water evaporation process - - - 6.741 Total of mass & heat flow 3.3375 869.085 3.3375 870.728 Table 5. Plant materials and heat flows per kg of clinker produced. Parameters Inlet Outletmass (kg) heat (kcal) mass (kg) heat(kcal) Clinker from kiln 1.1500 428.488 - - Cooling air 1.3239 9.440 - - Clinker product - - 1.0000 37.289 Return clinker dust to kiln - - 0.1500 34.201 Combustion air - - 1.3239 328.863 Convection & radiation loss - - - 37.375 Total of mass & heat flow 2.4739 437.928 2.4739 437.728 Table 6. Materials and heat flows of cooler per kg of clinker produced. Parameters Inlet Outletmass (kg) heat (kcal) mass (kg) heat(kcal) Kiln feed from SP 1.5636 346.300 - - Coal to main burner 0.1370 837.102 - - Combustion air 1.5189 330.393 - - Return clinker dust from cooler 0.1500 34.201 - - Produced clinker - - 1.1500 428.488 Recirculating clinker dust to SP - - 0.1700 50.630 Flue gas to SP - - 1.6381 537.000 Kiln feed gas - - 0.4114 134.900 Net of calcination and clinkerization heat - - - 293.310 Evaporating of coal moisture content - - - 3.846 Convection & radiation loss - - - 92.730 Total of mass & heat flow 3.3695 1547.996 3.3695 1540.904 Table 7. Materials and heat flows of kiln per kg of clinker produced. 212 vol. 61 no. 1/2021 New method evaluation of detail material and heat flows. . . Parameters Inlet Outletmass (kg) heat (kcal) mass (kg) heat(kcal) Kiln feed 1.6155 20.442 - - Recirculating clinker dust from kiln 0.1700 50.630 - - Flue gas from kiln 1.7042 541.660 - - Kiln feed gas from kiln 0.4114 134.900 - - Kiln feed to the kiln - - 1.5636 346.300 Combustion gas including excess air - - 1.7042 154.620 Kiln feed gas out - - 0.5528 46.744 Return dust flow - - 0.0805 6.887 Heat of calcination process - - - 131.800 Evaporating of kiln feed water content - - - 2.894 Convection & radiation loss - - - 63.180 Total of mass & heat flow 3.9011 747.632 3.9011 752.425 Table 8. Materials and heat flows of SP per kg of produced clinker. kiln is returned back to the SP. By limiting the difference between the heat flows in and out of the kiln to ≤ 1 %, we can estimate the percentage of the kiln feed calcined in the SP and in the kiln. The calculation results of this heat conservation show that about 25 % of the kiln feed is calcined in the SP while the remaining 75 % are calcinated in the kiln. In addition, using Equations (20) and (21), we can evaluate that the heat released during the sintering process in the kiln is about 102.2 kcal/kg of clinker produced. We can also analyse that the heat carried by the clinker to the cooler, and by the gas from the fuel combustion and calcination of the kiln feed out from the kiln to the SP is very large. In the cooler, the heat carried by the clinker is mostly recovered by the combustion air. In the SP, it is used for heating and the initial calcination process of the kiln feed. Moreover, it also appears that the heat loss to the environment (i.e. 92.73 kcal/kg of produced clinker) is also significant (i.e., greater than in the cooler). The value obtained is still in the range of the results of Radwan (2012) (50 − 140 kcal/kg of produced clinker) [5]. In general, the specific heat loss will decrease with an increasing plant capacity. However, there are still opportunities to recover this heat loss for the purpose of increasing the plant efficiency. 3.4. SP and detail each cyclone mass and heat conservations The focus of this research is to estimate the separation efficiency of cyclones in the SP with the smallest possible error of the difference between the flows of heat in and out. This certainly cannot be obtained by a direct measurement in the field. Hence, by using Equations (55) to (60) based on the data from gas and material temperature measurements at the inlet and outlet of each cyclone, the value of the separation efficiency approximation can be obtained. The results of calculating the mass and heat conservation in the SP are given in Table 8. From the results, we can conclude that the difference between the outgoing and incoming heat rates is insignificant (0.64 %). The percentage of calcined kiln feed in this SP is around 25 % with around 8 % in cyclone 3 and the remaining 17 % in the lowest cyclone. Using Equations (52) to (54) and completion methods as described in equations (55) to (57), the mass and heat conservations in each cyclone can be obtained with the difference between the heat in and out being less than 1 %. The mass and heat conservations in cyclones 1, 2 and 3, 4 are given in Tables 9 and 10, respectively. From the results of these calculations using Equations (58) to (60), the estimated separation efficiency for cyclone 2 is 91.89 %, 84.09 % for cyclone 3 and 79.51 % for the lowest cyclone. The separation efficiency of the top cyclone has been evaluated by Eq. (6). The results of the mass and heat conservation calculation for each cyclone is highly useful for designing equipment, controlling the process and modifying equipment (e.g., estimating the rate of kiln feed flow to the calciner). Some studies that simulate the process in a calciner [22, 51] require data of the entering kiln feed flow, which generally come out of cyclone 3. Unfortunately, these studies did not explain the origin of this data, which of course cannot be measured directly at the plant while the plant is operating. With the results of the calculations as shown in Table 9, the process simulation that occurs within the calciner is expected to be closer to the real conditions, because the data can be obtained in a more accurate way. The results of calculating the separation efficiency of these cyclones depend on the measured temperature data. For example, in cyclone 3, the variation in temperature of material coming out of the cyclone 2 as well as the gas entering the cyclone 3 is ±2 °C, which practically does not affect the value of its separation efficiency. The lower value of the separation efficiency for the lower cyclone position is a common design referred to in cement plants. The top cyclone is designed with a high efficiency since the separation function is more dominant. 213 P. S. Darmanto, I M. Astina, A. K. Wardhana et al. Acta Polytechnica Parameters Cyclone 1 Cyclone 2mass (kg) heat (kcal) mass (kg) heat(kcal) Kiln feed 1.6155 20.670 - - Kiln feed from upper cyclone - - 1.6966 138.396 Kiln feed dust from lower cyclone 0.1665 23.406 0.3572 64.028 Flue gas from lower cyclone 1.7042 247.467 1.7042 311.981 Kiln feed gas from lower cyclone 0.5479 76.000 0.5479 97.440 Flue gas from kiln - - - - Kiln feed gas from kiln - - - - Total of mass & heat flow in - 367.543 4.3059 611.845 Separated kiln feed to the lower cyclone 1.6966 138.396 1. 8873 257.199 Unseparated dust to upper cyclone/duct 0.0805 6.887 0.1665 23.406 Kiln feed gas 0.5528 47.648 0.5479 76.000 Flue gas including excess air 1.7042 153.543 1.7042 247.467 Kiln feed moisture evaporating process - 2.894 - - Heat of calcination process - - - - Convection & radiation heat loss - 20.550 - 11.630 Total of mass & heat flow out 4.0341 367.543 4.3059 615.702 Table 9. Materials and heat flows per kg of produced clinker of Cyclones 1 and 2. Parameters Cyclone 3 Cyclone 4mass (kg) heat (kcal) mass (kg) heat(kcal) Kiln feed from upper cyclone 1.8873 257.199 1.8875 330.592 Kiln feed dust from lower cyclone 0.4029 92.839 0.1700 50.630 Flue gas from lower cyclone 1.7042 397.520 - - Kiln feed gas from lower cyclone 0.5024 115.400 - - Flue gas from kiln - - 1.7042 541.526 Kiln feed gas from kiln - - 0.4114 134.900 Total of mass & heat flow in 4.4968 862. 958 4.1731 1057.648 Separated kiln feed to the lower cyclone 1.8875 330.592 1.5636 346.306 Unseparated dust to upper cyclone/duct 0.3572 64.028 0.4029 92.839 Kiln feed gas 0.5479 97.440 0.5024 115.400 Flue gas including excess air 1.7042 311.981 1.7042 397.520 Kiln feed moisture evaporating process - - - - Heat of calcination process - 43.950 - 87.890 Convection & radiation heat loss - 12.330 - 18.670 Total of mass & heat flow out 4.4968 860.321 4.1731 1058.625 Table 10. Materials and heat flows per kg of produced clinker of Cyclones 3 and 4. 214 vol. 61 no. 1/2021 New method evaluation of detail material and heat flows. . . In addition to controlling the process, a high efficiency of the top cyclone is also intended to reduce the amount of pollutant emissions released into the environment. Whereas for the 3 lower cyclones, their heat transfer function is more important. From Tables 8 and 9, it can be seen that the heat carried by the gas coming out of the top-cyclone is significantly high (more than 200 kcal/kg of clinker produced). This heat can be recovered for drying fine coal, lime stone and clay as reported by Alfi et.al [52]. 4. Conclusions Equations and a new approach to precisely calculate the mass and heat conservations in a single string cement plant are presented in this paper. Using the proposed new method, the detail mass and heat conservation in each main equipment (e.g., cooler, kiln, SP and cyclones) can be evaluated based on common data available in the control room and measuring the surface temperature of the main equipment. As a case study, operational data of Tonasa 2 cement plant were used for the detailed mass and heat conservation evaluation. The conclusions that can be drawn from the results of this study are as follows: (1.) The advantage of this proposed method over the conventional method of calculating the mass and energy flow of a clinker plant is that it can estimate the mass and energy flow in each cyclone so that the separation efficiency of each cyclone, which is impossible to measure, can be evaluated. This cannot be done with conventional methods, which only evaluate the SP as a whole. The estimated separation efficiency value of each cyclone can be obtained with less than 1 % error in heat conservation. (2.) This proposed method can also be used to estimate several parameters that cannot be obtained by a direct measurement of a running plant (e.g. combustion air temperature at the exit from the cooler, and clinkerization heat). (3.) The results obtained are expected to be used as additional data in controlling operations, designing new equipment and modifying the processes and dimensions of an existing equipment. (4.) The proposed calculation method can be applied to all types of modern cement clinker plant configurations, either with or without a calciner including the double strings. It is also expected to contribute to design process improvement, operation simulation of a calciner and controlling the daily operation of a clinker plant. List of symbols A surface area [m2] ash Mass percentage of ash in coal [%] En Heat flow rate [kcal/s, kW] m Mass flow rate [kg/s] NHV Nett heating value [kcal/kg, kJ/s] η Materials separation efficiency of the cyclone [%] LOI Loss of ignition [%] H2O Mass flow rate of vapor [kg/s] hx(T x) Enthalpy of substance x at temperature Tx [kcal/kg] hf g Enthalpy of evaporation [kcal/kg] hkf Enthalpy of kiln feed [kcal/kg] (CaCO3)kf Mass percentage of CaCO3 [%] (MgCO3)kf Mass percentage of MgCO3 [%] (SiO2)kf Mass percentage of SiO2 [%] (Al2O3)kf Mass percentage of Al2O3 [%] (Fe2O3)kf Mass percentage of Fe2O3 [%] T Temperature [řC] Qloss Radiation and convection heat loss [kcal/s, kW] %Calc_k Percentage of kiln feed calcined in the kiln [%] %Calc_SP Calcined kiln feed percentage in SP [%] ∆Ecyclone Error of energy conservation law (the difference of heat in-flow and heat out-flow of the cyclone) ∆mcyclone Error of mass conservation law (the difference of mass in and mass out of the cyclone) Subscripts amb ambient ash ash in coal c clinker cooler calc calcination process 215 P. S. Darmanto, I M. Astina, A. K. Wardhana et al. Acta Polytechnica C Cyclone Separator C1 Cyclone 1 (top) C2 Cyclone 2 C3 Cyclone 3 C4 Cyclone 4 (lowest) cli clinker clid_to_k Clinker dust return to kiln coal coal coal − comb coal combustion process cool − air clinker cooling air cool_air_k cooling air to kiln evap evaporation process form clinker formation process gas − kf kiln feed gas hg hot gas hgi hot gas exit from Cyclone i H2O − kf kiln feed evaporated water k Kiln kf kiln feed kfi separated kiln feed from Cyclone i SP Suspension Preheater sec_air_c Combustion air from cooler to kiln sep − kf separated kiln feed sep − kfc separated kiln feed to form clinker sep − kf − i separated kiln feed by cyclone i (i = 1A and 1B, 2, 3 and 4) surf surface area tr − air transporting air tot whole plant unsep_kf_i unseparated kiln-feed by cyclone i, ( i = 1A and 1B, 2, 3 and 4) vapor water vapor vapor − kf vapor from kiln feed References [1] R. Virendra, B. Sudheer Prem Kumar, J. Suresh Babu, D. Rajani Kant. Detailed energy audit and conservation in a cement plant. International Research Journal of Engineering and Technology 2(1):248 – 256, 2015. [2] N. Anantharaman. Energy audit in cement industry (1500 tpd). International Journal of Science Technology & Engineering 3(10):12 – 18, 2017. [3] S. Sattari, A. Avami. Assessment of energy-saving opportunities of cement industry in Iran. In Proceeding of the 3rd IASME/WSEAS International Conference on Energy, Environment, Ecosystems and Sustainable Development, pp. 585 – 593. Citeseer, Agios Nikolaos, Greece, 2007. [4] P. Khongprom, U. Suwanmanee. Environmental benefits of the integrated alternative technologies of the Portland cement production: A case study in Thailand. Engineering Journal 21(7):15 – 27, 2017. doi:10.4186/ej.2017.21.7.15. [5] A. M. Radwan. Different possible ways for saving energy in the cement production. Advances in Applied Science Research 3(2):1162 – 1174, 2012. [6] T. T. Ayu, M. H. Hailu, F. Y. Hagos, S. M. Atnaw. Energy audit and waste heat recovery system design for a cement rotary kiln in Ethiopia: A case study. International Journal of Automotive and Mechanical Engineering 12:2983 – 3002, 2015. doi:10.15282/ijame.12.2015.14.0249. [7] M. Parmar, D. Solanki, B. Vegada. Energy and Exergy analysis of cement rotary kiln. International Journal of Advanced Engineering and Research Development 3(4):284 – 293, 2016. [8] A. I. Okoji, D. E. Babatunde, A. N. Anozie, J. A. Omoleye. Thermodynamic analysis of raw mill in cement industry using aspen plus simulator. In IOP Conference Series: Materials Science and Engineering, vol. 413, p. 012048. 2018. doi:10.1088/1757-899X/413/1/012048. [9] A. Kolip, A. F. Savas. Energy and exergy analyses of a parallel flow, four-stage cyclone precalciner type cement plant. International Journal of Physical Sciences 5(7):1147 – 1163, 2010. doi:10.5897/IJPS.9000219. [10] S. Jonnalagadda, S. Reddy. Heat transfer analysis of recuperative air preheater. International Journal of Engineering Research and Management 4(1):105 – 110, 2017. 216 http://dx.doi.org/10.4186/ej.2017.21.7.15 http://dx.doi.org/10.15282/ijame.12.2015.14.0249 http://dx.doi.org/10.1088/1757-899X/413/1/012048 http://dx.doi.org/10.5897/IJPS.9000219 vol. 61 no. 1/2021 New method evaluation of detail material and heat flows. . . [11] L. K. Nørskov. Combustion of solid alternative fuels in the cement kiln burner. Ph.D. thesis, Technical University of Denmark, 2012. [12] H. Mikulčić, E. von Berg, M. Vujanović, N. Duić. Numerical study of co-firing pulverized coal and biomass inside a cement calciner. Waste Management & Research 32(7):661 – 669, 2014. doi:10.1177/0734242X14538309. [13] C. Y. H. Chao, P. C. W. Kwong, J. H. Wang, et al. Co-firing coal with rice husk and bamboo and the impact on particulate matters and associated polycyclic aromatic hydrocarbon emissions. Bioresource Technology 99(1):83 – 93, 2008. doi:10.1016/j.biortech.2006.11.051. [14] Verein Deutscher Zementwerke. Activity Report: Process Technology of Cement manufacture. Utilisation of used tyres in cement works. Tech. rep., Verein Deutscher Zementwerke, 2003 - 2005. [15] S. Kourounis, S. Tsivilis, P. E. Tsakiridis, et al. Properties and hydration of blended cements with steelmaking slag. Cement and Concrete Research 37(6):815 – 822, 2007. doi:10.1016/j.cemconres.2007.03.008. [16] M. Varma, P. Gadling. Additive to cement, A pozzolanic material-fly ash. International Journal of Engineering Research 3(5):558 – 564, 2016. doi:10.17950/ijer/v5i3/010. [17] D. A. Y. Ghassan K. Al-Chaar, L. A. Kallemeyn. The Use of Natural Pozzolan in Concrete as an Additive or Substitute for Cement. ERDC/CERL TR-11-46. Tech. rep., Federal University of Technology, Minna, Niger. [18] A. Allahverdi, S. Salem. Studies on main properties of ternary blended cement with limestone powder and microsilica. Iranian Journal of Chemical Engineering 4(1):3 – 13, 2007. [19] D. Paa, P. Darmanto. Studi Numerik Pengaruh Laju Umpan Kiln terhadap Rugi Tekanan dan Efisiensi Pemisahan Top Siklon Suatu Pabrik Semen. In Proceeding Seminar Nasional Tahunan Teknik Mesin XIII (SNTTM XIII), pp. 625 – 630. 2014. [20] P. S. Darmanto, I. M. Astina, A. Syahlan. Design and Implementation of Deduster for Improving Fine Coal Quality in a Cement Plant. In The International Conference on Fluid and Thermal Heat Conversion (FTEC). Tongyeong, South Korea, 2009. [21] H. Mikulčić, E. Von Berg, M. Vujanović, N. Duić. Numerical study of co-firing pulverized coal and biomass inside a cement calciner. Waste management & research 32(7):661 – 669, 2014. doi:10.1177/0734242X14538309. [22] A. C. Kahawalage, M. C. Melaaen, L.-A. Tokheim. Numerical modeling of the calcination process in a cement kiln system. In Linköping Electronic Conference Proceedings, vol. 138, pp. 83 – 89. Reykjavik, Iceland, 2017. doi:10.3384/ecp1713883. [23] G. Borsuk, J. Wydrych, B. Dobrowolski. Modification of the inlet to the tertiary air duct in the cement kiln installation. Chemical and Process Engineering 37(4):517 – 527, 2016. doi:10.1515/cpe-2016-0042. [24] Claus Bauer. Modernization and production increase with cement kilns. Humbolt report. Tech. rep., KHD Humboldt Wedag AG, 2000. [25] X. D’Hubert. Latest burner profiles. Global Cement Magazine 03, 2017. [26] N. Gopani, A. Bhargava. Design of high efficiency cyclone for tiny cement industry. International Journal of Environmental Science and Development 2(5):350 – 254, 2011. doi:10.7763/IJESD.2011. V2.150. [27] H. Mikulčić, E. Von Berg, M. Vujanović, et al. Numerical analysis of cement calciner fuel efficiency and pollutant emissions. Clean Technologies and Environmental Policy 15(3):489 – 499, 2013. doi:10.1007/s10098-013-0607-5. [28] A. T. Gebremariam. CFD Modelling and Experimental Testing of Thermal Calcination of Kaolinite Rich Clay Particles - An Effort towards Green Concrete. Ph.D. thesis, Aalborg University, Denmark, 2015. [29] Y. Sonavane, E. Specht. Numerical analysis of the heat transfer in the wall of rotary kiln using finite element method ANSYS. In Seventh International Conference on CFD in the Minerals and Process Industrie CSIRO, pp. 1 – 5. Melbourne, Australia, 2009. [30] E. Copertaro, A. A. Estupinan Donoso, B. Peters. A discrete-continuous method for predicting thermochemical phenomena in a cement kiln and supporting indirect monitoring. Engineering Journal 22(6):165 – 183, 2018. doi:10.4186/ej.2018.22.6.165. [31] Z. Dulaimi, H. Hameed, A. Ali, M. Alfahham. Investigation the effect of calcinations degree and rotary kiln gases bypass opining in the preheating system for dry cement industries. International Journal of Latest Trends in Engineering and Technology 10(3):119 – 127, 2018. doi:10.21172/1.103.21. [32] F. L. Smidth. Preheater calciner systems. www.flsmidth-prod-cdn.azureedge.net/-/media/brochures/brochures-products/pyro/2000-2017/ preheater-calciner-systems.pdf?rev=3519484b-6b80-47ee-811a-b69f889f353f, 2011. [33] A. Atmaca, R. Yumrutaş. Analysis of the parameters affecting energy consumption of a rotary kiln in cement industry. Applied Thermal Engineering 66(1-2):435 – 444, 2014. doi:10.1016/j.applthermaleng.2014.02.038. [34] S. B. Nithyananth, H. Rahul. Thermal heat audit of kiln system in a cement plant. International Journal of Modern Engineering Research 5(12):73 – 79, 2015. 217 http://dx.doi.org/10.1177/0734242X14538309 http://dx.doi.org/10.1016/j.biortech.2006.11.051 http://dx.doi.org/10.1016/j.cemconres.2007.03.008 http://dx.doi.org/10.17950/ijer/v5i3/010 http://dx.doi.org/10.1177/0734242X14538309 http://dx.doi.org/10.3384/ecp1713883 http://dx.doi.org/10.1515/cpe-2016-0042 http://dx.doi.org/10.7763/IJESD.2011. V2.150 http://dx.doi.org/10.1007/s10098-013-0607-5 http://dx.doi.org/10.4186/ej.2018.22.6.165 http://dx.doi.org/10.21172/1.103.21 www.flsmidth-prod-cdn.azureedge.net/-/media/brochures/brochures-products/pyro/2000-2017/preheater-calciner-systems.pdf?rev=3519484b-6b80-47ee-811a-b69f889f353f www.flsmidth-prod-cdn.azureedge.net/-/media/brochures/brochures-products/pyro/2000-2017/preheater-calciner-systems.pdf?rev=3519484b-6b80-47ee-811a-b69f889f353f http://dx.doi.org/10.1016/j.applthermaleng.2014.02.038 P. S. Darmanto, I M. Astina, A. K. Wardhana et al. Acta Polytechnica [35] P. Prasanth, G. Sudhakar. Analysis of heat loss in kiln in cement industry - A review. In International Conference on Explorations and Innovations in Engineering & Technology (ICEIET –2016). [36] D. Touil, H. Belabed, C. Frances, S. Belaadi. Heat exchange modeling of a grate clinker cooler and entropy production analysis. International Journal of Heat and Technology 23(1):61 – 68, 2005. [37] L. Farag. Energy and exergy analyses of egyptian cement kiln plant with complete kiln gas diversion through by pass. International Journal of Advances in Applied Sciences 1(1):35 – 44, 2012. doi:10.11591/ijaas.v1i1.757. [38] F. L. Smidth. Dry Process kiln. FL Smidth Inc, USA, 2004. [39] W. H. Duda. Cement Data Book I: International Process in the Cement Industry. French & European Pubs, 3rd edn., 1985. [40] D. K. Fidaros, C. A. Baxevanou, C. D. Dritselis, N. S. Vlachos. Numerical modelling of flow and transport processes in a calciner for cement production. Powder Technology 171(2):81 – 95, 2007. doi:10.1016/j.powtec.2006.09.011. [41] F. L. Smidth. Plant Services Devision. Heat Conservations. International cement production seminar, Lecture 5.13A. Tech. rep., FL Smidth Inc, 1990. [42] T. Chatterjee. Burnability and clinkerization of cement raw mixes. In Advances in Cement Technology, pp. 69 – 113. Pergamon Press, 1983. doi:10.1016/B978-0-08-028670-9.50009-0. [43] P. A. Aisop, H. Chen, H. Tseng. The cement plant operations handbook. Tradeship Publications Ltd, Surrey, UK, 5th edn., 2007. [44] Holcim Group Support Ltd. Cement Manufacturing Services. Reference Guide for Process Performance Engineers. Thermal Process and Materials Technology, Edition 3.0. Tech. rep., Holcim Group Support Ltd., 2006. [45] A. Howard, R. Chris. Elementary Linear Algebra. John Wiley and Sons Inc., United States, 9th edn., 2005. [46] I. Markovsky. Low rank approximation - Algorithms, implementation, applications. Springer-Verlag London, 2012. doi:10.1007/978-1-4471-2227-2. [47] J. Chen, M. Shi. Analysis on cyclone collection efficiencies at high temperatures. China Particuology 1(1):20 – 26, 2003. doi:10.1016/S1672-2515(07)60095-5. [48] A.-N. Huang, N. Maeda, D. Shibata, et al. Influence of a laminarizer at the inlet on the classification performance of a cyclone separator. Separation and Purification Technology 174:408 – 416, 2017. doi:10.1016/j.seppur.2016.09.053. [49] A. Amalia, A. Syahlan, P. S. Darmanto. Heat auditing of Gresik and Tonasa plants. Internal Project Report of Indonesian Cement and Concrete Institute. Tech. rep., Indonesian Cement and Concrete Institute, 2006. [50] A. Amalia, A. Syahlan, P. S. Darmanto. Design and implementation of hot gas system for raw coal drying in Tonasa 3 Plant. Internal Project Report of Indonesian Cement and Concrete Institute. Tech. rep., Indonesian Cement and Concrete Institute, 2006. [51] H. Mikulčić, E. Von Berg, M. Vujanović, et al. CFD analysis of a cement calciner for a cleaner cement production. Chemical Engineering Transactions 29:1513 – 1518, 2012. doi:10.3303/CET1229253. [52] A. Amalia, A. Syahlan, P. S. Darmanto. Modification of hot gas utilization for drying lime stone and clay at Tonasa 2 and 3 plants. Final Project Report of Indonesian Cement and Concrete Institute. Tech. rep., Indonesian Cement and Concrete Institute, 2007. 218 http://dx.doi.org/10.11591/ijaas.v1i1.757 http://dx.doi.org/10.1016/j.powtec.2006.09.011 http://dx.doi.org/10.1016/B978-0-08-028670-9.50009-0 http://dx.doi.org/10.1007/978-1-4471-2227-2 http://dx.doi.org/10.1016/S1672-2515(07)60095-5 http://dx.doi.org/10.1016/j.seppur.2016.09.053 http://dx.doi.org/10.3303/CET1229253 Acta Polytechnica 61(1):199–218, 2021 1 Introduction 2 Materials and method 2.1 Materials 2.2 Mass and heat conservation equations 2.3 Methodology 2.4 Mass and heat conservation equations of planetary clinker cooler 2.5 Mass and heat conservation equations of the kiln 2.6 Detail mass and heat conservation equations of SP system and each cyclone separator 3 Results and discussion 3.1 Operational data parameters 3.2 Plant mass and heat conservations 3.3 Cooler and kiln mass and heat conservations 3.4 SP and detail each cyclone mass and heat conservations 4 Conclusions List of symbols References