Thermal Analysis and Energy Efﬁciency Improvements in Tunnel Kiln for Sustainable Environment

: Kiln is a prime need in the ceramics industry, where energy loss is a major part which consumes about 60% production cost through thermal energy for different applications. Higher density of ﬁred and tunnel kiln refractory material lowers the thermal diffusivity and the proper selection of ﬁred material minimizes the energy loss along the kiln. In particular, this research analysed the results of a heat recovery system comprised of a metallic recuperator which gives around 8% energy savings in natural gas consumption. In this work, detailed power quality analysis of low-power factor motors of a tunnel kiln was carried out and a power factor improvement solution was suggested to save electrical energy with payback period of 0.8 y. The motor operating at a low-power factor consumes more reactive power which does not produce beneﬁcial work. A low-power factor around 0.4 causes network power loss, increases in transformer loss and voltage drops. The solution with accumulative capacitance power of 148.05 uF was installed to achieve the power factor to 0.9. Flu gas analyzer was installed to monitor the range of O 2 in pre-heating, oxidation, and ﬁring zones of the kiln which should be ≥ 8% and 3%, respectively. Regression analysis for thermal energy consumption of a tunnel kiln is done to ﬁnd the forecast thermal energy consumption. This analysis can be used to ﬁnd operational efﬁciency, supporting decisions regarding dependent variable of thermal energy consumption and independent variable of production. This research is very helpful for the ceramics industry to mitigate the energy loss at SMEs as well as in mass production level.


Introduction
Firing and drying process of ceramics manufacturing requires intensive thermal energy consumption. These processes can be improved by improving the product quality, By using the suitable thermal insulation, it was revealed that approximately 35% of energy gain and reduction of temperature inside the wall from 249.34 • C to 79.49 • C can be achieved [22]. Recovered air mass flow optimization study revealed that reduction of natural gas by 4.6% daily consumption can be achieved by passing the air to pre heat the sanitary ware from cooling to preheating zone [23]. In ceramics manufacturing, the efficiency of a tunnel kiln is not more than 30% due to its convective and radiation losses. The mass ratio between tunnel kiln car and fired material should be more than 50%. Different thermal properties of fired material and constructive material can influence the specific energy losses [24]. The simulation of a dimension-less model in preheating zone predicting gas flow, evaporation of bound water, and normative heat transfer between gas was developed. The temperature profile results in entering the ambient air in preheating zone and not entering the preheating zone to monitor the quality of the fired bricks examined [25].
AC induction motors usually work with a load much less as compared to their maximum rating. Many surveys show that the average load factor for these motors is 57-68% of their ratings, and machine often operate with more less load, even in between 3-16% [20]. It was observed that 44% of motors installed in industrial sector are loaded at less than 40% of their rated power [26]. It is observed that AC motors have the highest efficiency at 75% of load [27]. The main purpose is to save energy by improving the power factor of the motor of the tunnel kiln [28]. This can be achieved by decreasing the distribution losses and improving the electrical load operation for a better energy efficiency level. The power factor improvement of an industrial AC motor can be achieved by installing the power factors. For industrial operation, AC induction motors are widely used. It is estimated that 70-80% is the average efficiency of the electricity consumed by these motors in the whole world. Power factor correction is gained by adding the capacitive load to offset the inductive load in the power network [29,30]. By reducing the tunnel channel height, aerodynamic characteristics such as pressure and velocity can be changed by modeling. The uniform velocity distribution along the tunnel length was visualized in a simulation. [31]. The experimental analysis and thermodynamics calculation of a kiln for ceramic porcelain tiles was carried out [32]. Exergy analysis of ceramics kilns expressed that 83% of accumulative energy input of kiln was demolished and only 10% of exergy input can be controlled in the physio-chemical transformation of the entire ceramic configuration. Over 61% of heat recovery can be achieved by increasing the entire energy performance of the kiln [33]. Specific energy consumption of table ware ceramics kiln was analyzed and determined to be 19,363 MJ/t. It was observed that there is a linear relation between kiln loading capacity and thermal energy consumption LPG, an increase in one kg of product of ceramics per m 3 loading can reduce the energy consumption by 0.67 kg LPG per kg. Therefore, it was observed that increasing the loading capacity of the kiln, reducing the actual heat capacity, and repairing the kiln insulator which reduce the specific thermal energy consumption of the kiln [34].
The investigation of the sanitary ware ceramics industry tunnel kiln was done, and it was observed that 76% of thermal energy consumed in this particular kiln. The rest of the energy was consumed in the drying process of the ceramics industry [35]. Best available techniques (BAT) for energy efficiency for ceramics production such as installation of superlative burners to reduce the thermal energy of the entire kiln, a large number of well-controlled burners which can optimize the temperature profile in entire tunnel kiln, installation of metallic recuperator heat exchangers to save thermal energy consumption, introduction of high loading rate and reduction in kiln energy consumption, and introduction of environmental management structure to monitor the flu gas regularly were mentioned [36]. The convective heat transfer analysis of tunnel kiln was investigated, and three turbulence models were accounted for in this study: k-e, k-w and RNG k-w. It was observed that 17% of convective heat transfer enhanced by the longitudinal bricks arrangement and 27% was enhanced by transverse lattice setting arrangement of bricks in the tunnel kiln [37]. The regression analysis depends on the dependent and independent variable and their correlation, interception and slope [38].
The thermogravimetric analysis of gaseous fuels with high volatile matter was discussed to enhance the system efficiency [39]. A parametric study of inside pressure and its study was well defined [40]. Turbulence flow of gaseous fluid and its characteristic with solid particles in the enclosed chamber were introduced [41]. The detailed analysis of nano particles and their energy distributive with respect to solar thermal power was defined [42]. Development of solid-based absorbent and effective carbon capture applications with respect to a sustainable environment was discussed [43]. Evaporator cooling and its effect on environmental sustainability in Pakistan was introduced for form house [44]. Environmental sustainability is very important because greenhouse gases such as CO and CO 2 continuously damage the ozone layer. The stack analysis of kilns is very important to monitor CO and CO 2 [45]. The detail of heat transfer mechanism in heat exchangers was discussed. In the kiln, it is very important to maintain the pressure draught for air and fuel chemistry. The main factor is that by increasing the heat transfer, there will be a high range of pressure drop [46].
Techno-economic analysis of thermal energy consumption gives a better idea to predict future energy consumption demand. There is almost 8% savings, which has been achieved after installing the heat recovery system. a detailed analysis of thermal energy consumption of the tunnel kiln has been carried out, and the results of pre, post and forecast are mentioned in this paper. Thermal energy consumption from July 2019 to January 2020 is expressed in a linear relation. Thermal energy consumption in July 19 was 4190 MMBtu, whereas in August 2019, the thermal energy consumption was 4158 MMBtu. January 2020, the thermal energy consumption is slighter higher i.e., 4199 MMBtu. The production values from July 2019 to January 2020 are in the range of 178t to 195t. The techno-economic analysis predicts the future thermal energy consumption, which gives the statistic view for a tunnel kiln. The forecast future consumption from Febuary-2020 to September 2020 is in the range of 4150 to 4189 MMBtu.
In the present research, thermal energy consumption analysis for pre-and post-heat recovery system at the tunnel kiln were conducted. Power quality analysis of tunnel kiln motors were also conducted to suggest a power factor improvement solution with a payback period of 0.8 months. Regression analysis for thermal energy consumption of the tunnel kiln was also done to find the forecast thermal energy consumption.

Materials and Methods
The chemical composition of ceramics body comprises aluminum oxide Al2O3 (15%), P-I (14%), potash S-2 gulf (17%), potash feldspar S-3 (4%), TA clay (22%), China clay D.W (13%), Black clay HK 3 4 (10%), and A A-7 (5%). These composition materials with described percentages are mixed in ball mills for 13 h to prepare the ceramics body for firing in the tunnel kiln. Results of shrinkage for six batches showed good compromise between fired maturity and shrinkage. Porcelain material has the highest level of fired shrinkage compared to stoneware. The set standard for the test of shrinkage was 10.5 ± 1.5% for diameter wise, and 8 ± 1% for length wise. The mechanical strength of green, glazed, and unglazed batches is mentioned below. Table 1 shows the properties of ceramics bodies which were fired in tunnel kiln. The set testing standard for the green body, un-glazed, and glazed batches were 50 ± 10, 850 ± 100. Experimentally, it is observed that when the ceramics body is heated up, its size and volume is increased in small increments, which is called thermal expansion. The values of thermal expansion for six batches are mentioned below. The set testing standard for thermal expansion was 0.38 ± 0.02. tunnel kiln car heating capacity is taken Cp 980 J/kgK, density ρ car 1250 kg/m 3 , thermal conductivity λ car 1 W/m K, total cycle time of tunnel kiln 70 h, height of the tunnel kiln car h car 1.1 m, initial temperature 303 K, length 2.25 m, useful length 1.08 m, velocity of the kiln car 3.0 m/h.  Table 2 shows the properties of natural gas and air in the tunnel kiln, in which detail of gas fuel and air contents is discussed. This table helps to identify the air flow mass rate in the tunnel kiln.

Experimental Setup
The recovered air optimization from cooling zone to firing zone resulted in reduction of natural gas consumption. This mechanism provides the waste heat recovery using the air of 167 • C from tunnel kiln cooling zone to preheating zone burners. Detail of the installed system is mentioned in the Figure 1. The experimental setup is of a tunnel kiln of the ceramics industry under steady-state condition. Some assumptions are made to conduct the analysis, i.e., gas forms are ideal, combustion during the analysis is complete, and ambient temperature is 30 °C. For preheating zone of kiln, the combustion air, the tunnel kiln is equipped with a metallic radiation recuperator with separate design for two groups which is arranged above the cooling zone. The recuperator is shielded against the combustion ducts by the ceiling with bay. It consists of heat-resistant steel tubes and has a collecting box on either end. The recuperator is attached by a tunnel iron structure by two bracket plates. The plates near the burning zone rest loosely on two angles at the furnace roof so as to allow free expan- The experimental setup is of a tunnel kiln of the ceramics industry under steadystate condition. Some assumptions are made to conduct the analysis, i.e., gas forms are ideal, combustion during the analysis is complete, and ambient temperature is 30 • C. For preheating zone of kiln, the combustion air, the tunnel kiln is equipped with a metallic Processes 2021, 9, 1629 6 of 21 radiation recuperator with separate design for two groups which is arranged above the cooling zone. The recuperator is shielded against the combustion ducts by the ceiling with bay. It consists of heat-resistant steel tubes and has a collecting box on either end. The recuperator is attached by a tunnel iron structure by two bracket plates. The plates near the burning zone rest loosely on two angles at the furnace roof so as to allow free expansion, while the rear end is rigidly mounted. The Preheating Zone Mass balance of air is given as follows: .
Typical heat transfer equation in the tunnel kiln preheating zone is written below The above equation gives the proper idea of heat loss by air in the preheating zone which is always equal to the sanitary ware heat transfer and inner area of the kiln by the method of convection method, plus the heat loss due to evaporation of water from sanitary ware and heat absorbed by the evaporated water.
The water evaporation from sanitary ware is expressed by the diffusion of the moisture that is equal to: Diffusion coefficient is governed as: D is the coefficient of diffusion for water vapours and tortuosity of the sanitary ware is equal to The Firing Zone Firing zone mass balance of the tunnel kiln is expressed as: .
where in the location of group burners: .
Mass balance between two group of burners: In the firing zone, the heat loss by the air and gained by the air is equal to the accumulative heat transfer between the sanitary ware and air and in between the wall and air, in addition to heat released by natural gas combustion: Processes 2021, 9, 1629 7 of 21 The heat equation in the burner group location: And in between the group of two burners: The Cooling Zone Air mass balance in the cooling zone is expressed below whereas, in the equation, the last term is the air mass flux that is injected in the cooling zone. Air mass flux . m a i is greater than zero, and the air mass extracted . m a e in the cooling zone is less than zero.
The heat mass balance in the cooling zone of tunnel kiln is expressed as In the above equation, the heat attained by the air is equal to convective heat transfer between the sanitary ware, air and wall, plus a source at the right side of equation.
At cooling zone injection position of the tunnel kiln At the extraction position of cooling zone of the tunnel kiln The rates of all contents of material in physical energy is stated below, 'h' is expressed as enthalpy and 'T' is the required temperature and 'cp' is expressed as specific heat capacity, where 'i' means ith components in equation mentioned above, Highest activity which can be obtained by mass and heat transfer is known as exergy, and by this, the system reaches equilibrium with the environmental conditions. Due to the effect of irreversibility the exergy can be destroyed or consumed that can arise in the system as divergent to energy conservation.
The exergy equation of concerned material is expressed in the general way mentioned below,Ė x =Ėx ph +Ėx ch (18) By this analysis we know that 'Ėxch' is expressed as the chemical exergy rate and 'Ėxph' is known as physical exergy rate and The physical exergy equation of formed gas can be stated aṡ where si, p, and R are the gas constant and 'o' is expressed as the environmental state in the equation above.
The overall gas content specific chemical exergy is calculated in the equation mentioned belowē Molar gas constant is denoted by R and 'x' is molar rate. Specific chemical exergy of hydrocarbon fuel such as 'Ca Hb' is formulated. In Figure 2, the conductivity of the tunnel kiln car is examined with respect to energy loss. The tunnel kiln furniture is made up of the refractory material and constructive material and the selected kiln is also made of refractory blocks and steel structure. Experimentally, it is observed that by increasing the conductivity of the tunnel kiln, car energy loss increased. It is observed that the thermal conductivity of insulated material is less than the constructed material. Thermal conductivity of the material is ranged from 0.02 W/m K to 2 W/m K. The thermal conductivity of the material of the tunnel kiln is 0.5 W/m K at 1300 • C, and the energy loss will be 0.43 MJ/kg. Therefore, by this graph, it is clear that if the thermal conductivity is low, there is low energy loss in the tunnel kiln. In general, the most important part of this research is to conserve energy. Figure 3, represents the heat transfer along the cycle time for a tunnel kiln. After 10 h, heat transfer value was 8 W/m 2 K, whereas after 20 h, the heat transfer was increased only 20% i.e., 10 W/m 2 K. After 30 h in firing zone, the heat transfer value was quite similar and only 1 W/m 2 K increment was observed. Convective heat transfer will increase by increasing the temperature of the tunnel kiln.

Effect of Heat Transfer for Different Cycle Times of Kiln
than the constructed material. Thermal conductivity of the material is ranged from 0.02 W/m K to 2 W/m K. The thermal conductivity of the material of the tunnel kiln is 0.5 W/m K at 1300 °C, and the energy loss will be 0.43 MJ/kg. Therefore, by this graph, it is clear that if the thermal conductivity is low, there is low energy loss in the tunnel kiln. In general, the most important part of this research is to conserve energy.  Figure 3, represents the heat transfer along the cycle time for a tunnel kiln. After 10 h, heat transfer value was 8 W/m 2 K, whereas after 20 h, the heat transfer was increased only 20% i.e., 10 W/m 2 K. After 30 h in firing zone, the heat transfer value was quite similar and only 1 W/m 2 K increment was observed. Convective heat transfer will increase by increasing the temperature of the tunnel kiln.     The recommended standard value for kiln draft pressure inside is 1.0 mm H 2 O, so to maintain the pressure up to the standard value, installation of the heat recovery system some experiments were performed. After installation of the heat recovery system, it was observed that the kiln pressure range was under range. Before installing the heat recovery system, the values of kiln pressure were 0.85, 0.88, 0.87, 0.84, 0.85, 0.89 mm H 2 O which was not appropriate for standard operating of the tunnel kiln. However, after installing the heat recovery system, the pressure inside the tunnel kiln was 1.06, 1.02, 1.01, 1.05, 1.06 mmH 2 O which was under range as per standard for the tunnel kiln. It is mandatory to maintain the kiln pressure up to 1.0 mmH 2 O for complete draught inside the kiln. If the range of the kiln pressure is not up to standard limit, the natural gas flow will be disturbed.

Effect of Heat Transfer for Different Cycle Times of Kiln
was not appropriate for standard operating of the tunnel kiln. However, after installing the heat recovery system, the pressure inside the tunnel kiln was 1.06, 1.02, 1.01, 1.05, 1.06 mmH2O which was under range as per standard for the tunnel kiln. It is mandatory to maintain the kiln pressure up to 1.0 mmH2O for complete draught inside the kiln. If the range of the kiln pressure is not up to standard limit, the natural gas flow will be disturbed.

Effect of Density of Ceramics Material of Tunnel Kiln on Energy Loss
In Figure 5, it is observed that by increasing the density of refractory material, the energy loss decreased. Industrial data showed that the density of the material fired in the tunnel kiln was higher than 2000 kg/m 3 . By increasing the density of the fired material, higher conductive resistance occurs. If the solid density increased, the thermal diffusivity of the car will be decreased. As a result, by increasing the density of material, the heat transfer along the kiln car will decrease and energy loss will minimize. In Figure 5, it is observed that an increase in density of ceramics material will result in the energy loss being decreased. One other recommendation is that if we reduce the slab size of the tunnel kiln car and increase the density of the tunnel kiln car, then thermal energy consumption along the kiln will reduce. If there is a proper selection of material having a density of 2000 kg/m 3 , the energy loss will be minimized, i.e., 0.39 MJ/kg.

Effect of Density of Ceramics Material of Tunnel Kiln on Energy Loss
In Figure 5, it is observed that by increasing the density of refractory material, the energy loss decreased. Industrial data showed that the density of the material fired in the tunnel kiln was higher than 2000 kg/m 3 . By increasing the density of the fired material, higher conductive resistance occurs. If the solid density increased, the thermal diffusivity of the car will be decreased. As a result, by increasing the density of material, the heat transfer along the kiln car will decrease and energy loss will minimize. In Figure 5, it is observed that an increase in density of ceramics material will result in the energy loss being decreased. One other recommendation is that if we reduce the slab size of the tunnel kiln car and increase the density of the tunnel kiln car, then thermal energy consumption along the kiln will reduce. If there is a proper selection of material having a density of 2000 kg/m 3 , the energy loss will be minimized, i.e., 0.39 MJ/kg. The thermal energy consumption along the temperature profile of the tunnel kiln is mentioned in the Figure 6. By increasing the temperature of kiln, the thermal energy consumption of tunnel kiln is in an increasing trend. The thermal capacity of the burner installed in the tunnel kiln is 2000 kcal/h. The total length of the ceramics kiln is taken 100 m, which includes 30 m pre-heating zone, 20 m firing zone, 30 m cooling zone. Total 44 tunnel kiln cars were traveling in the kiln. The process is totally sequential because the process of the tunnel kiln is always continuous. Different products are loaded on the tun-

Effect of Temperature Profile and Natural Gas Consumption of Tunnel Kiln
The thermal energy consumption along the temperature profile of the tunnel kiln is mentioned in the Figure 6. By increasing the temperature of kiln, the thermal energy consumption of tunnel kiln is in an increasing trend. The thermal capacity of the burner installed in the tunnel kiln is 2000 kcal/h. The total length of the ceramics kiln is taken 100 m, which includes 30 m pre-heating zone, 20 m firing zone, 30 m cooling zone. Total 44 tunnel kiln cars were traveling in the kiln. The process is totally sequential because the process of the tunnel kiln is always continuous. Different products are loaded on the tunnel kiln cars, which is heated up at the maximum temperature and cooled down at ambient temperature. Natural gas consumption along the tunnel kiln temperature is explained below. In Figure 6, thermal energy consumption along the temperature profile of the tunnel kiln is explained. There are three zone of the tunnel kiln i.e., pre-heating, oxidation and firing zone. Thermal energy consumption at the temperature of 202 • C to 560 • C in preheating is 0.591 MMBtu/h. Thermal energy consumption at a temperature of 742 • C to 939 • C in oxidation zone is (0.739~0.750) MMBtu/h. Thermal energy consumption at a temperature of 1076 • C to 1183 • C in firing is 1.846~1.898 MMBtu/h.

Effect of Temperature Profile and Natural Gas Consumption of Tunnel Kiln
The thermal energy consumption along the temperature profile of the tunnel kiln is mentioned in the Figure 6. By increasing the temperature of kiln, the thermal energy consumption of tunnel kiln is in an increasing trend. The thermal capacity of the burner installed in the tunnel kiln is 2000 kcal/h. The total length of the ceramics kiln is taken 100 m, which includes 30 m pre-heating zone, 20 m firing zone, 30 m cooling zone. Total 44 tunnel kiln cars were traveling in the kiln. The process is totally sequential because the process of the tunnel kiln is always continuous. Different products are loaded on the tunnel kiln cars, which is heated up at the maximum temperature and cooled down at ambient temperature. Natural gas consumption along the tunnel kiln temperature is explained below. In Figure 6, thermal energy consumption along the temperature profile of the tunnel kiln is explained. There are three zone of the tunnel kiln i.e., pre-heating, oxidation and firing zone. Thermal energy consumption at the temperature of 202 °C to 560 °C in preheating is 0.591 MMBtu/h. Thermal energy consumption at a temperature of 742 °C to 939 °C in oxidation zone is (0.739~0.750) MMBtu/h. Thermal energy consumption at a temperature of 1076 °C to 1183 °C in firing is 1.846~1.898 MMBtu/h.

Flu Gas Analysis in Pre-Heating and Firing Zone
Experimentally, oxygen analysis was done by Flu-gas analyser, which has a basic O 2 and CO sensing probe. The Measuring range of the instrument for O 2 is up to 29 percent volumetrically and for CO, the measuring range is up.to 2000 ppm with a resolution of one ppm [38]. Table 3 describes the experimental values of different zones in the tunnel kiln which have been compared with the international standard. The level of O 2 was analysed before installation of heat recovery system in oxidation, preheating, reduction, and firing zones of the tunnel kiln. It was observed that the level of O 2 before installation of the heat recovery system was out of range. The standard values for O 2 in the pre-heating and oxidation zone should not be greater than ≥8%. O 2 must be zero in this zone. The standard values for O 2 in firing zone should not be greater than 3%. It has been shown in the below table that the O 2 in firing zone was a little bit out of range before installing the heat recovery system. It was observed that the level of O 2 , after installation of the heat recovery system are in range. The standard values for O 2 in the pre-heating and oxidation zone should not be greater than or equal to 8%. The standard values for O 2 must be zero in this zone. The standard values for O 2 in the firing zone should not greater than 3%. It has been shown in the below table that the O 2 in natural zone is in range after installing the heat recovery system. The detailed analysis of exergy and energy approach is performed on the system of the selected tunnel kiln. The waste heat recovery is achieved by adding the recovered air mass from cooling zone to preheating zone of tunnel kiln. The energy consumption required for the firing zone of the tunnel kiln is determined under normal conditions. Before installing the heat recovery system, the natural gas consumption and combustion air are 5.80 m 3 /h and 0.33 m 3 /h. The total consumption of natural gas in the tunnel kiln was 4171 MMBtu per month before installing the heat recovery system.
Heat recovery system gives the decrease in natural gas consumption to 0.29 m 3 /h and 8% natural gas saving has been achieved. The heat recovery system was installed in February 2020. For the preheating of the combustion air tunnel, the kiln is equipped with a metallic radiation recuperator with separate design for two groups which is arranged above the cooling zone. The recuperator is shielded against the combustion ducts from the ceiling with bay. It consists of heat-resistant steel tubes and has a collecting box on either end. The recuperator is attached with a tunnel iron structure by two bracket plates. The plates near the burning zone rest loosely on two angles at the furnace roof so as to allow free expansion, while the rear end is rigidly mounted. Figure 7 shows the comparison of thermal energy consumption before and after installing the heat recovery system. The detailed analysis of thermal consumption is carried out, and thermal energy consumption in June 2019 was 4190 MMBtu and thermal consumption in July 2019 was 4158 MMBtu, which is 7.63% lower. Thermal energy consumption in August 2019 was 4165 MMBtu, which is 5.96% lower. The average value of thermal energy consumption after installing the heat recovery system is 8.5% lower, which is a significant value. The values of February-2020 to September-2020 are at low trends in terms of thermal energy consumption.

Power Quality Analysis of Tunnel Kiln Motors
In Figure 8, detailed power quality analysis of all the motors was carried out at field. Motors of entry fan, combustion air fan 1, 2, and 3 have a low power factor and 8%, 17.27%, 20%, and 42.86% loaded. These motors were found oversized and underloaded which were replaced with more energy efficient, properly sized motors. To save energy, a power factor improvement was suggested, and there was a potential to save 8% energy by investing in a payback period of 0.8 y.
Efficiency of induction motors varies with the variation of load. Some motors utilize power efficiently at rated conditions. However, other induction motors waste a huge portion of energy when they run normally. Efficiency of the motor can be described by power factor (PF) at specific load scaling from 0 to 1. The motor operating at a low power factor consumes high reactive power, which does not give useful work. A low power factor increases the network loss, transformer losses and voltage drop in the system. Power supply companies impose penalties to industrial consumers on low power factor running loads. Industrial consumers encourage electricity supplier to improve their power factor for different reasons. The power factor is the ratio of real power P(kW) and S (kVA) total drawn power of the motor. The closer the ratio of P/S to one, the higher will be the efficiency. Thermal Energy Consumption (MMBtu) Before heat recovery system After heat recovery system Figure 7. Thermal energy consumption before and after heat recovery system.

Power Quality Analysis of Tunnel Kiln Motors
In Figure 8, detailed power quality analysis of all the motors was carried out at field. Motors of entry fan, combustion air fan 1, 2, and 3 have a low power factor and 8%, 17.27%, 20%, and 42.86% loaded. These motors were found oversized and underloaded which were replaced with more energy efficient, properly sized motors. To save energy, a power factor improvement was suggested, and there was a potential to save 8% energy by investing in a payback period of 0.8 y. Efficiency of induction motors varies with the variation of load. Some motors utilize power efficiently at rated conditions. However, other induction motors waste a huge portion of energy when they run normally. Efficiency of the motor can be described by power Power factor is corrected by adding opposite reactive power by adding capacitors in parallel to load. The circuit impedance with a power factor compensation capacitor is required to improve the power factor of motor system. Capacitance C(F) is the value of capacitor needed to add to raise the power factor to 0.95, which gives the significant savings.
In Table 4, reactive power of capacitor is calculated, Capacitor of 49.47, 20.53, 37.32, and 40.73 uF value is required for 7.5, 1.5, 1.5, and 2.2 kW to improve the power factor of required motors. The current power factor of these low power factors is 0.20, 0.23, 0.26, and 0.39 which needs to improve by using these capacitor banks. In Figure 9, the detailed quality analysis has been done to evaluate the real, apparent, and reactive power of the tunnel kiln motors. The rated power and the percentage-loaded condition is analysed to find the power factor of the tunnel kiln motors which are under normal conditions. Thirteen motors of different ratings are analysed to find the results of percentage loading. Motor ratings of 7.5 kW are around 40-48% loaded at the tunnel kiln. The detailed power quality analysis is performed to check the quality index of the tunnel kiln motors. Motor ratings of 2.2 kW are 22-58% loaded at the tunnel kiln. The electrical load of the motors covers around 30% load of the total energy consumption of the tunnel kiln. Major barriers which affects the implementation of motor efficiency measures are operational issues and allocation of resources in the industry. Power factor improvement systems increases the system efficiency as well as the electrical energy savings. Capacitors of 49.47, 20.53, 37.32, and 40.73 uF value are required for 7.5, 1.5, 1.5, and 2.2 kW to improve the power factor upto 0.95 of tunnel kiln motors. The cost of power factor improvement for these low-power factor motors is 110,000 pkr including all taxes, and the typical annual savings achieved after installation is about 162,153 pkr per annum. These motors are running continuously 24 h/d and 350 d/y. The kiln remains closed for annual maintenance and cleaning for 2 w/y. Simple payback can be achieved by dividing the investment cost of the power factor improvement system and savings per annum. This solution is viable to achieve the power factor improvement up to 0.95 with a suitable payback period of 0.8 y.

Regression Analysis of Thermal Consumption of Tunnel Kiln
In Figure 10, techno-economic analysis of thermal energy consumption and its pro- Major barriers which affects the implementation of motor efficiency measures are operational issues and allocation of resources in the industry. Power factor improvement systems increases the system efficiency as well as the electrical energy savings. Capacitors of 49. 47, 20.53, 37.32, and 40.73 uF value are required for 7.5, 1.5, 1.5, and 2.2 kW to improve the power factor upto 0.95 of tunnel kiln motors. The cost of power factor improvement for these low-power factor motors is 110,000 pkr including all taxes, and the typical annual savings achieved after installation is about 162,153 pkr per annum. These motors are running continuously 24 h/d and 350 d/y. The kiln remains closed for annual maintenance and cleaning for 2 w/y. Simple payback can be achieved by dividing the investment cost of the power factor improvement system and savings per annum. This solution is viable to achieve the power factor improvement up to 0.95 with a suitable payback period of 0.8 y.

Regression Analysis of Thermal Consumption of Tunnel Kiln
In Figure 10, techno-economic analysis of thermal energy consumption and its production data of tunnel kiln ceramics was analysed. It is observed that the regression during the set period of June 2019 to January 2020 was examined critically. The regression line during this period is linear and the values of the dependent variable (thermal energy consumption) are denoted with Y-axis, and value of production independent variable is denoted with X-axis. The regression expression is mentioned below is y = 3.2799x + 3567.4. The correlation between X-variable and Y-variable is R = 3.2733, Interception between X-variable and Y-variable is 3567.4 [47]. Regression analysis for thermal energy consumption of the tunnel kiln is useful to finding the forecast consumption. This analysis can be used to find predict analytics, operational efficiency, and supporting decisions regarding two dependent and independent variables. In Figure 10, the relationship of production (Tons) and thermal energy consumption (MMBtu) of studied tunnel kiln is linear. Thermal energy consumption from July 2019 to January 2020 is expressed in linear relation. Thermal energy consumption in July 19 was 4190 MMBtu, whereas in August 19, the thermal energy consumption was 4158 MMBtu. While in January 20 the thermal energy consumption is only slighter higher, i.e., 4199 MMBtu. The production values from July 2019 to January 2020 are in the range of 178 t to 195 t. The techno-economic analysis predicts the future thermal energy consumption, which gives the statistic view for the tunnel kiln.

Regression of Production Data
In Figure 11, detailed regression analysis of production data for the months of June 2019 to January 2020 is done. The figure shows the linear trend across the thermal consumption and production of tunnel kiln. The slope value of the analysis is 3.273 ± 0.7403, intercept value is 3567 ± 136.6, and the regression value is 0.8747. In Figure 10, the relationship of production (Tons) and thermal energy consumption (MMBtu) of studied tunnel kiln is linear. Thermal energy consumption from July 2019 to January 2020 is expressed in linear relation. Thermal energy consumption in July 19 was 4190 MMBtu, whereas in August 19, the thermal energy consumption was 4158 MMBtu. While in January 20 the thermal energy consumption is only slighter higher, i.e., 4199 MMBtu. The production values from July 2019 to January 2020 are in the range of 178 t to 195 t. The techno-economic analysis predicts the future thermal energy consumption, which gives the statistic view for the tunnel kiln.

Regression of Production Data
In Figure 11, detailed regression analysis of production data for the months of June 2019 to January 2020 is done. The figure shows the linear trend across the thermal consumption and production of tunnel kiln. The slope value of the analysis is 3.273 ± 0.7403, intercept value is 3567 ± 136.6, and the regression value is 0.8747.

Residual Analysis of Regression
In Figure 12, residual analysis of regression is carried out, and value of percentile 0 to 90 and residual −25 to +25 is selected for better outcome of results. The independent value of production and dependent value of thermal energy consumption is selected for residual analysis. which gives the statistic view for the tunnel kiln.

Regression of Production Data
In Figure 11, detailed regression analysis of production data for the months of June 2019 to January 2020 is done. The figure shows the linear trend across the thermal consumption and production of tunnel kiln. The slope value of the analysis is 3.273 ± 0.7403, intercept value is 3567 ± 136.6, and the regression value is 0.8747. Figure 11. Regression of production data. Figure 11. Regression of production data. In Figure 12, residual analysis of regression is carried out, and value of percentile 0 to 90 and residual −25 to +25 is selected for better outcome of results. The independent value of production and dependent value of thermal energy consumption is selected for residual analysis.     Figure 13 shows the comparison of forecast thermal energy consumption and post heat recovery energy consumption. This comparison gives the clear idea of significant savings in thermal energy consumption after installing the heat recovery system on tunnel kiln.
In Table 5, the detailed analysis of thermal energy consumption of the tunnel kiln has been carried out, the results of pre-, post-and forecast are mentioned. Techno-economic analysis of thermal energy consumption gives a better idea to predict future energy consumption demand. There is almost 8% savings, which has been achieved after installing the heat recovery system. Thermal energy consumption from July 2019 to January 2020 is expressed in linear relation. Thermal energy consumption in July 2019 was 4190 MMBtu, whereas in August 2019, the thermal energy consumption was 4158 MMBtu. In January 2020, the thermal energy consumption is only slighter higher i.e., 4199 MMBtu. The production values from July 2019 to January 2020 are in the range of 178 t to 195 ton. The techno-economic analysis predicts future thermal energy consumption, which gives a statistic view for the tunnel kiln. The forecast future consumption from February 2020 to September 2020 is in the range of 4150 to 4189 MMBtu.  Figure 13 shows the comparison of forecast thermal energy consumption and post heat recovery energy consumption. This comparison gives the clear idea of significant savings in thermal energy consumption after installing the heat recovery system on tunnel kiln. In Table 5, the detailed analysis of thermal energy consumption of the tunnel kiln has been carried out, the results of pre-, post-and forecast are mentioned. Techno-economic analysis of thermal energy consumption gives a better idea to predict future energy consumption demand. There is almost 8% savings, which has been achieved after installing the heat recovery system. Thermal energy consumption from July 2019 to January 2020 is expressed in linear relation. Thermal energy consumption in July 2019 was 4190 MMBtu, whereas in August 2019, the thermal energy consumption was 4158 MMBtu. In January  In Table 6, the detailed summary output of first regression is analyzed. The multiple regression is a significant relation of two set independent and dependent variable which is 0.121283. Variation analysis of 16 Nos observations has been analysed, and significant value of regression is achieved i.e, 0.6545557.
In Table 7, the residual output and probability output of thermal energy consumption have been analyzed. This analysis helps to predict the future energy consumption in the tunnel kiln. The values of predicted percentile and residuals are mentioned below.

Conclusions
Installing the heat recovery system on the cooling zone to preheat the ware of preheating zone of the tunnel kiln gives 8% natural gas consumption. O 2 analysis after installing the heat recovery system is in range which clearly depicts that the combustion efficiency is improved after installing the heat recovery system. Mechanical strength, shrinkage (diameter, length wise), and thermal expansion of proper selected ceramics body material are in range of standard values. Density of selected material i.e., (2.1~2.3 g/cm 3 ) gives better results to improve the thermal energy consumption along the tunnel kiln. Detailed power quality analysis of low power factor motors of the tunnel kiln was carried out and a power factor improvement solution was suggested to save electrical energy with a payback period of 0.8 y. Regression analysis for thermal energy consumption of the tunnel kiln is done to find the forecast thermal energy consumption of the tunnel kiln. This analysis gives the predicted analytics, regarding dependent variable thermal energy consumption and independent production variables. The model will be validated with a numerical simulation. It is necessary to re-investigate the experimental setup of the kiln and continue this study for a variable air velocity and inlet gas at different zones. The results obtained with simulation should be compared with experimental results.