Mathematical Methods for Heat Transfer and Thermodynamic Analysis of Conductive, Convective, and Radiative Media
View this Special IssueResearch Article  Open Access
Bin Zheng, Yongqi Liu, Lichen Zou, Ruiyang Li, "Heat Transfer Characteristics of Calcined Petroleum Coke in Waste Heat Recovery Process", Mathematical Problems in Engineering, vol. 2016, Article ID 2649383, 9 pages, 2016. https://doi.org/10.1155/2016/2649383
Heat Transfer Characteristics of Calcined Petroleum Coke in Waste Heat Recovery Process
Abstract
This paper reports the results of heat transfer characteristics of calcined petroleum coke in waste heat recovery process. The model of heat exchanger was set up. The model has been used to investigate the effects of porosity (0.58 to 0.79), equivalent heat conductivity coefficient (0.9 to 1.1), and equivalent specific heat (0.9 to 1.1). The calculated values of calcined petroleum coke temperature showed good agreement with the corresponding available experimental data. The temperature distribution of calcined petroleum coke, the calcined petroleum coke temperature at heat exchanger outlet, the average heat transfer coefficient, and the heat recovery efficiency were studied. It can also be used in deriving much needed data for heat exchanger designs when employed in industry.
1. Introduction
With the rapid development of economy in China, energy consumption has been increasing on a vast scale in recent years. The equivalent coal consumption of China was 3.75 billion tons in 2013 [1]. To solve this problem, it is widely believed that energy recycle is an effective measure to decrease energy consumption. In China’s industrial production, waste heat resource is rich and various, such as waste heat from high temperature gas, coke slag, and other solid particles [2]. The researches on waste heat utilization of high temperature gas and liquid are relatively wide, and technology of waste heat recovery is relatively mature. However, the researches and technology of waste heat utilization of high temperature solid are not inadequate.
Calcined petroleum coke is one of the important basic raw materials. It is widely used in producing anode of aluminum electrolytic, graphite electrode, recarburizer, industrial silicon, and other carbon products. The production capacity of calcined petroleum coke in China is the largest in the world, over 70% of which is produced in tank calcined furnaces. When calcined petroleum coke is discharged from tank calcined furnaces, the temperature is 1000°C. The heat of high temperature calcined petroleum coke is about 33.5% of the total heat of calcination process [3]. The waste heat of it has high possibility to be utilized.
Waste heat reutilization of the solid particle has been studied in some field. Cai and Ding [4, 5] summarized the current energy consumption condition of steel industry and distribution of waste heat resource of steel production process. Wang et al. [6] analyzed the thermal equilibrium of molten slag heat recovery system and discussed the relationship between waste heat utilization efficiency and the amount of cooling air. The cooling process of steel slag particles in the moving bed was analyzed. Chen [7] simulated the waste heat recovery system of high temperature slag with fluent software, analyzed the flow field and temperature field of steel slag, and researched the change of temperature of air at outlet and heat recovery efficiency when the wind speed and air temperature at inlet are varied. Liu [8] proposed a new technology of fluidized bed. The waste heat of slag can be efficiently used with this technology. The flow and heat transfer process were simulated in fluidized bed with CFD software. Herz et al. [9] researched the heat transfer coefficient between moving material layer and the wall of bed through experiment. The effects of operational parameters, moving speed of material layer, and depth of buried material on heat transfer coefficient were studied. Atmakidis and Kenig [10] simulated the movement of irregularly arranged particle near the wall with CFD. Sundaresan and Kolar [11] researched the axial heat transfer characteristics in circulating fluidized bed and analyzed surface heat transfer coefficient of several vertical pipes. The results showed that the length and location of the tube affected the heat transfer coefficient. Guo et al. [12] simulated the heat transfer in gassolid fluidized bed with EulerEuler method and obtained the effects of Reynolds number and porosity on heat transfer.
The investigations on the reutilization of calcined petroleum coke waste heat are lacking. In this paper, heat transfer characteristics of calcined petroleum coke in waste heat recovery process were studied.
2. Physical and Mathematical Description of the Problem
The waste heat of calcined petroleum coke is recovered by a heat exchanger. The heat exchanger of calcined petroleum coke includes an internal heat exchanger and an external heat exchanger. Its structure is shown in Figure 1. The external heat exchanger consists of tubes and a diaphragm wall, and the ends of tubes are held together by two ring pipes. Similarly the internal heat exchanger also consists of tubes, with the ends of tubes held together by a straight pipe. The calcined petroleum coke flows from top to bottom in the heat exchanger surrounded by tubes and diaphragm walls, and the water flows from bottom to top in the tubes. The temperature of calcined petroleum coke is high; the water temperature is low; so the heat is transferred to water from calcined petroleum coke.
In this paper, the calculated model of calcined petroleum coke heat exchanger is equal to the experimental heat exchanger model in scale. Considering the structure characteristic of the heat exchanger, the calculation model is appropriately simplified.(1)The lower ring pipe and the upper ring pipe of the internal heat exchanger and the external heat exchanger will not be calculated.(2)The calculated model includes two external heat exchanger tubes, one internal heat exchanger tube, and the region between the tubes and thermal barriers.(3)The tube pitch is equal from top to bottom. Additionally the symmetry plane of calculation model is the plane where the short axis of the oval heat exchanger exists.(4)The calcined petroleum coke particles of high temperature and the gas among them are regarded as homogenous continuum. The physical parameters of calculation model are calculated as the porosity of calcined petroleum coke.(5)Considering the thermal insulation layer outside the heat exchanger, it is assumed that there is no heat loss from the heat exchanger lateral.
According to the simplification above, the model has been set up by Gambit. The lengths of and axis are 78 mm and 262 mm. The height of the model is 1550 mm. It is shown as Figure 2.
When heat exchanges between the high temperature calcined petroleum coke and water, the mass, momentum, and energy conservation exist. In this paper, the unsteady model is used for calculation. The mass, momentum, and energy conservation equations are shown as follows.
The mass conservation equation is as follows:where is density, kg/m^{3}, is time, s, and , , and are the velocity vectors in , , and coordinate directions, m/s.
The momentum conservation equation is as follows:where is pressure on the infinitesimal body, Pa, , , and are stick stress component of on the infinitesimal body caused by molecular viscosity, N/m^{2}, , , and are body forces, N/m^{3}, and is velocity vector quantity, m/s.
The energy conservation equation is as follows:where is specific heat of water, J/(kg·K), is effective heat transfer coefficient, W/(m·K), is temperature of fluid, K, and is viscous dissipation term.
The equations of the flow model in the pipe are as follows.
Turbulence intensity is
Turbulent kinetic energy is
Turbulence kinetic energy dissipation rate iswhere is turbulence intensity, is Reynolds number of fluid, is average fluid turbulent velocity, m/s, is empirical constant, 0.09, and is characteristic length, m.
According to the simplified model, the heat transfer of the calcined petroleum coke in the heat exchanger mainly includes the heat transfer among calcined petroleum coke particles and the heat transfer between calcined petroleum coke and the wall of the heat exchanger. The heat transfer process in this paper is steady with no inner heat source, so the heat conduction differential equation is as follows:where is density, kg/m^{3}, is specific heat, J/(kg·K), is temperature, K, and is coefficient of heat conductivity, W/(m·K).
The discrete points within the domain are solved. The implicit format of the unsteady heat conduction equation is used to solve problems in this paper:where is time step, s, is temperature, K, is the serial number of node; is the serial number of temporal sequence; is thermal diffusivity, m^{2}/s, and is distance between two adjacent nodes, m.
The heat transfer equation of calcined petroleum coke and the wall of heat exchanger is as follows:where is heat flux, W, is heat transmission coefficient, W/(m^{2}·K), is heat exchange area, m^{2}, and is average temperature difference, K.
The calcined petroleum coke velocity, water velocity in the tube, and calcined petroleum coke temperature in the heat exchange inlet are constant. The calcined petroleum coke particles and the gas among them are regarded as homogenous continuum. The physical parameters of calculation model are calculated as the porosity of calcined petroleum coke. The effective heat conductivity coefficient of homogenous continuum is calculated using (10) [13]:where is effective thermal conductivity coefficient, W/(m·K), is porosity, is thermal conductivity coefficient of gas, W/(m·K), and is thermal conductivity coefficient of solid, W/(m·K).
The effective specific heat of homogenous continuum is calculated usingwhere is effective specific heat, J/(kg·K), is porosity, is specific heat of gas, J/(kg·K), and is specific heat of solid, J/(kg·K).
When the heat exchanger model is meshed, the internal boundary is refined locally. The quantity and quality of the grid have a great influence on the calculation results. If the grids are much less, it cannot meet the precision request. If the grids are too dense, it will require sufficient computer memory. And the calculating time will increase. So it is critical to confirm the quantity of grids. The test of grid independence has been widely used to solve this problem. There are four schemes for meshing. The quantity of grids is 130442, 215817, 313721, and 467604. When calculation is completed, the temperature of calcined petroleum coke at outlet is used as the evaluation standard. The calculation results are shown as in Table 1.

Table 1 showed that the relative change of Scheme is within 2% compared to Scheme . There will be little influence on calculation results if grids are much dense. Therefore Scheme not only reaches the calculation accuracy but also saves much time. So Scheme was used for calculation.
The velocity of water in tube is 1 m/s, the inlet temperature of water is 300 K, the inlet temperature of calcined petroleum coke is 1073 K, and the velocity of calcined petroleum coke in the heat exchanger is 6 × 10^{−5 }m/s.
The heat recovery efficiency is calculated bywhere is heat recovery efficiency, %, is water flow mass of 1st external heat exchange tube, kg·s^{−1}, is enthalpy value of water in 1st tube outlet, kJ·kg^{−1}, is enthalpy value of water in 1st tube inlet, kJ·kg^{−1}, is water flow mass of 2nd external heat exchange tube, kg·s^{−1}, is enthalpy value of water in 2nd tube outlet, kJ·kg^{−1}, is enthalpy value of water in 2nd tube inlet, kJ·kg^{−1}, is water flow mass of internal heat exchange tube, kg·s^{−1}, is enthalpy value of water in internal heat exchange tube outlet, kJ·kg^{−1}, is enthalpy value of water in internal heat exchange tube inlet, kJ·kg^{−1}, is specific heat capacity, kJ·kg^{−1}·K^{−1}, is calcined petroleum coke flow mass, kg·s^{−1}, and is temperature of calcined petroleum coke when it flows into the heat exchanger, K.
The definition of the equivalent heat conductivity coefficient is the ratio of the actual heat conductivity coefficient to the standard heat conductivity coefficient. The standard heat conductivity coefficient is the actual heat conductivity coefficient of the calcined petroleum coke porosity which is 0.7.
The definition of the equivalent specific heat is the ratio of the actual specific heat to the standard specific heat. The standard specific heat is the actual specific heat of the calcined petroleum coke porosity which is 0.7.
3. Experimental System
The experimental system of the waste heat utilization exchanger was built in Weifang Lianxing New Materials Technology Co., Ltd. The heat transfer experiments were carried out by the experimental system (Figure 3). The experimental system is composed of the waste heat exchanger, the calcined petroleum coke supply system, the water circulation system, and the measurement system.
The experimental system is installed in the tank calcined furnace. The high temperature calcined petroleum coke is supplied directly from the tank calcined furnace. The water circulation system is composed of cooling pond, pumps, valves, sight glass, a drum, and so on. There is a WSMD twophase flowmeter in the heat exchanger outlet, which is used to measure steam dryness and flow. There is a FLUXUS F601 ultrasonic flowmeter in the heat exchanger inlet, which is used to measure the water flow. Its measuring accuracy is ±0.5%. Its measuring velocity range is from 0.01 m/s to 25 m/s. Its repeatability is 0.15%.
The measurement system is composed of a temperature sensor, a pressure sensor, and a piece of data acquisition instrument. The temperature sensor is used to measure the temperature of water in the heat exchanger inlet, the temperature of calcined petroleum coke in heat exchanger, the temperature of steam in the heat exchanger outlet, and so on. The pressure sensor is used to measure the water pressure and the steam pressure. The data acquisition instrument is used to record the data. The temperature measuring points of calcined petroleum coke in the heat exchanger are shown in Figure 4.
4. Results and Discussions
4.1. Model Validation
The experimental data were used to validate the developed model. Simulations were conducted for the same operational conditions as those employed in the experimental investigation. Figure 5 shows the temperatures of numerical and experimental results in different planes when the porosity is 0.7. The temperature difference of numerical and experimental results is 16 K and the relative error is 1.78%. The results are in good agreement with the experimental value obtained in the same conditions.
4.2. Temperature Distribution Characteristics in the Heat Exchanger
Figure 6 shows the temperature distribution of calcined petroleum coke in the heat exchanger for the porosity which is 0.7. Figure 6 indicates that the temperature of calcined petroleum coke is low near the wall of the heat exchanger tube and the temperature differences are small in different planes. The temperature difference of the first plane and the fifth plane is only 19 K, while, in the center of the internal and external heat exchangers, the temperature of calcined petroleum coke is high and the temperature differences of different planes are great. The temperature difference of the first plane and the fifth plane is 359 K. The main reason is that the coefficient of the heat conductivity of calcined petroleum coke is very low. The longer the distance of the heat transfer is, the bigger the heat resistance of calcined petroleum coke is and accordingly the higher the temperature is.
4.3. Effect of Calcined Petroleum Coke Porosity
Figure 7 shows the effects of porosity on calcined petroleum coke temperature at the outlet of the heat exchanger. Figure 7 shows that when the porosity increases from 0.58 to 0.79, the average temperature of calcined petroleum coke decreases from 541 K to 489 K and the maximum temperature decreases from 650 K to 582 K. The main reason is that, with the increase of porosity, the effective specific heat and the effective density decrease and the initial total heat quantity reduces, which make the temperature of calcined petroleum coke decrease at the heat exchanger outlet. However, with the increase of porosity, the effective heat conductivity coefficient decreases. So the ability of the heat conductivity decreases, which makes the temperature of calcined petroleum coke increase at the heat exchanger outlet. But the effect of the initial total heat quantity reduction is a major influence factor, so the calcined petroleum coke temperature decreases with the increase of porosity.
Figure 8 shows the effects of porosity on average heat transfer coefficient. As can be seen from Figure 8, the average heat transfer coefficient of the internal heat exchanger is higher than the external heat exchanger. When the porosity increases from 0.58 to 0.79, the average heat transfer coefficients of the internal and the external heat exchanger decrease by 6.4 and 5.6 W/(m^{2}·K), respectively. The reason for this is that, with the increase of porosity, the effective specific heat and effective heat conductivity coefficient decrease and the average heat flux of the heat exchanger tube wall decreases. At the same time, with the decrease of the effective specific heat, the initial total heat quantity reduces while the cooling rate of calcined petroleum coke increases, so the average heat transfer temperature difference decreases. However, the effect of heat exchanger tube average heat flux is a major influence factor, so the average heat transfer coefficient decreases with the increase of porosity.
Figure 9 shows the effects of porosity on heat recovery efficiency. Figure 9 shows that the heat recovery efficiency increases by 5.2% when the porosity increases from 0.58 to 0.79. The main reason is that when the porosity increases, the initial total heat quantity and the heat recovery quantity reduce. But the reduce rate of the initial total heat quantity is higher than that of the reduce rate of heat recovery quantity, and thus the heat recovery efficiency increases.
4.4. Effect of Equivalent Heat Conductivity Coefficient
The normal value of heat conductivity coefficient is the effective heat conductivity coefficient when the calcined petroleum coke porosity is 0.7. The definition of the equivalent heat conductivity coefficient is a ratio of the objective heat conductivity coefficient and the normal value.
Figures 10–12 show the effects of equivalent heat conductivity coefficient at the calcined petroleum coke temperature at the heat exchanger outlet, the average heat transfer coefficient, and the heat recovery efficiency. As can be seen from Figures 10–12, when the equivalent heat conductivity coefficient increases from 0.9 to 1.1, the average temperature reduces by 35 K and the maximum temperature decreases by 50 K, the average heat transfer coefficient of the internal and external heat exchanger increases by 2.8 and 2.6 W/(m^{2}·K), respectively, and the heat recovery efficiency increases by 4.3%. The reason for this is that, with the increase of the equivalent heat conductivity coefficient, the heat resistance of calcined petroleum coke decreases, the heat transfer quantity in unit time increases, the temperatures at the heat exchanger outlet reduce, and the average heat transfer coefficient and the heat recovery efficiency increase.
4.5. Effect of Equivalent Specific Heat
The normal value of specific heat is the effective specific heat when the calcined petroleum coke porosity is 0.7. The definition of equivalent specific heat is a ratio of the objective specific heat and the normal value.
Figure 13 shows the effects of equivalent specific heat on calcined petroleum coke temperature at the heat exchanger outlet. Figure 13 shows that when the equivalent specific heat increases from 0.9 to 1.1, the average temperature of calcined petroleum coke increases from 501 K to 535 K and the maximum temperature increases from 595 K to 644 K. The main reason is that when the equivalent specific heat increases, the initial total heat quantity heat increases and when the heat conductivity coefficient of calcined petroleum coke is a constant, the calcined petroleum coke temperature at the heat exchanger outlet increases.
Figure 14 shows the effects of equivalent specific heat on average heat transfer coefficient. As can be seen from Figure 14, when the equivalent specific heat increases from 0.90 to 1.10, the average heat transfer coefficients of the internal and external heat exchanger decrease by 0.16 and 0.28 W/(m^{2}·K), respectively. In other words, the effect is very small.
Figure 15 shows the effects of the equivalent specific heat on heat recovery efficiency. Figure 15 shows that the heat recovery efficiency decreases by 4.1% when the equivalent specific heat increases from 0.90 to 1.10. The main reason is that when the heat transfer time is a constant, with the increase of equivalent specific heat, the initial total heat quantity heat increases, the heat transfer sufficiency decreases, and the heat recovery efficiency reduces.
5. Conclusions
The main conclusions of the present study are as follows.(1)The temperature of calcined petroleum coke is low near the wall of the heat exchanger tube, and the temperature differences are small in different planes, while, in the center of the internal and external heat exchangers, the temperature of calcined petroleum coke is high and the temperature differences of different planes are great.(2)When the porosity increases from 0.58 to 0.79, the average temperature of calcined petroleum coke at the outlet of the heat exchanger decreases from 541 K to 489 K and the maximum temperature decreases from 650 K to 582 K. The average heat transfer coefficients of the internal and external heat exchanger decrease by 6.4 and 5.6 W/(m^{2}·K), respectively. The heat recovery efficiency increases by 5.2%.(3)With the increase of the equivalent heat conductivity coefficient (from 0.9 to 1.1), the average temperature reduces by 35 K and the maximum temperature decreases by 50 K, the average heat transfer coefficients of the internal and external heat exchanger increase by 2.8 and 2.6 W/(m^{2}·K), respectively, and the heat recovery efficiency increases by 4.3%.(4)If the other conditions remain the same, when the equivalent specific heat increases from 0.9 to 1.1, the average temperature of calcined petroleum coke increases from 501 K to 535 K and the maximum temperature increases from 595 K to 644 K. The average heat transfer coefficient of the internal and external heat exchanger decreases by 0.16 and 0.28 W/(m^{2}·K), respectively. The heat recovery efficiency decreases by 4.1%.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by Shandong Provincial Natural Science Foundation, China (ZR2013EEQ005), and Shandong Provincial Science and Technology Development Program, China (2013GGX10404).
References
 D. Bob, BP Statistical Review of World Energy, BP, London, UK, 2014.
 Q. Zhao, Y. Wang, and X. Wang, “Technical advances status of China's waste heat utilization,” Industrial Boiler, vol. 25, no. 5, pp. 8–15, 2009. View at: Google Scholar
 B. Zheng, Y. Liu, and Z. Wang, “Study on thermophysical properties of calcined petroleum coke,” Carbon Techniques, vol. 32, no. 2, pp. 33–35, 2013. View at: Google Scholar
 J.J. Cai, J.J. Wang, and C.X. Chen, “Recovery of residualheat integrated steelworks,” Iron & Steel, vol. 42, no. 6, pp. 1–7, 2007. View at: Google Scholar
 Y. Ding and D.M. Shi, “Highefficiency utilization of waste heat at fully integrated steel plant,” Iron & Steel, vol. 46, no. 10, pp. 88–98, 2011. View at: Google Scholar
 B. Wang, X.C. Wang, Y.C. Yuan, and Q.P. Zhou, “Advances in the study of the blast furnace slag waste heat recovery technologies,” Journal of Engineering for Thermal Energy and Power, vol. 29, no. 2, pp. 113–121, 2014. View at: Google Scholar
 Y. Chen, Numerical simulation study on waste heat recovery system of high temperature steel slag [M.S. thesis], Shandong University, Shandong, China, 2014.
 F. Liu, Research on a novel waste heat recovery equipment for blast furnace slag and heat transfer numerical simulation [M.S. thesis], Qingdao Technological University, Shandong, China, 2010.
 F. Herz, I. Mitov, E. Specht, and R. Stanev, “Experimental study of the contact heat transfer coefficient between the covered wall and solid bed in rotary drums,” Chemical Engineering Science, vol. 82, pp. 312–318, 2012. View at: Publisher Site  Google Scholar
 T. Atmakidis and E. Y. Kenig, “Numerical analysis of mass transfer in packedbed reactors with irregular particle arrangements,” Chemical Engineering Science, vol. 81, pp. 77–83, 2012. View at: Publisher Site  Google Scholar
 R. Sundaresan and A. K. Kolar, “Axial heat transfer correlations in a circulating fluidized bed riser,” Applied Thermal Engineering, vol. 50, no. 1, pp. 985–996, 2013. View at: Publisher Site  Google Scholar
 X.Y. Guo, H.S. Chai, and D.H. Chao, “Numerical simulation of large particle fluidized bed and comparison of gasparticle heat transfer models,” Journal of University of Shanghai for Science and Technology, vol. 34, no. 1, pp. 81–87, 2012. View at: Google Scholar
 W. Woodside and J. H. Messmer, “Thermal conductivity of porous media. I. Unconsolidated sands,” Journal of Applied Physics, vol. 32, no. 9, pp. 1688–1699, 1961. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Bin Zheng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.