World Academy of Science, Engineering and Technology International Journal of Mechanical and Mechatronics Engineering Vol:5, No:1, 2011 Simulation of Heat Transfer in the Multi-Layer Door of the Furnace U. Prasopchingchana workers [7] examined the combined conduction and mixed Abstract—The temperature distribution and the heat transfer convection flow field and temperature distribution in an rates through a multi-layer door of a furnace were investigated. The enclosure in the presence of vertical partition with finite inside of the door was in contact with hot air and the other side of the thickness and thermal conductivity. Nouanegue et al. [8] door was in contact with room air. Radiation heat transfer from the investigated conjugate heat transfer by natural convection, walls of the furnace to the door and the door to the surrounding area conduction and radiation in open cavities in which a uniform was included in the problem. This work is a two dimensional steady heat flux is applied to the inside surface of the solid wall state problem. The Churchill and Chu correlation was used to find local convection heat transfer coefficients at the surfaces of the facing the opening. Chiu and Yan [9] carried out a numerical furnace door. The thermophysical properties of air were the functions study in the radiation effect on the characteristics of the mixed of the temperatures. Polynomial curve fitting for the fluid properties convection fluid flow and heat transfer in inclined ducts. The were carried out. Finite difference method was used to discretize for effect of conduction of horizontal walls on natural convection conduction heat transfer within the furnace door. The Gauss-Seidel heat transfer in a square cavity is numerically investigated by Iteration was employed to compute the temperature distribution in Mobedi [10]. Molla and Hossain [11] examined the effect of the door. thermal radiation on a study two-dimensional mixed The temperature distribution in the horizontal mid plane of the convection laminar flow of viscous incompressible optically furnace door in a two dimensional problem agrees with the one thick fluid along a vertical wavy surface. dimensional problem. The local convection heat transfer coefficients This research aims to study the coupling between multi- at the inside and outside surfaces of the furnace door are exhibited. layer conduction, natural convection and surface radiation on Keywords—Conduction, heat transfer, multi-layer door, natural the multi-layer door of the furnace. The work used the convection Churchill and Chu [12] empirical correlation in the simulation to compute the average convection heat transfer coefficients. I. INTRODUCTION Later, the local convection heat transfer coefficients were EAT loss from a multi-layer door of a furnace is a determined. Then, radiation heat transfer coefficients were Hcombination problem of multi-layer conduction, natural computed. Finally, the temperature distribution and total heat convection and radiation. At present, several research works loss from the inside of the furnace to room air through the [1-11] have been devoted to investigate the conduction, multi-layer door of the furnace were investigated. natural convection and radiation phenomena. Prasopchingchana and Laipradit [1] performed a simulation of II. MODEL GOMETRY heat transfer in the glass panel of an oven using an empirical This research has been treated as a two dimensional steady- correlation for natural convection heat transfer on the surface state problem. There was hot air on the inside surface of the of the glass panel. Bilgen [2] carried out a numerical and multi-layer door and room air on the outside surface of the experiment study of conjugate heat transfer by conduction and multi-layer door. natural convection on a heated vertical wall. Simulation of natural convection and radiation in the gas cavity and heat Furnace conduction in the walls of an enclosure was performed by Kuznetsov and Sheremet [3]. Nouanegue and Bilgen [4] Tsur carried out a numerical study on vertical channel solar chimneys with a practical geometry, which consist of two Room air Hot air International Science Index, Mechanical and Mechatronics Engineering Vol:5, No:1, 2011 waset.org/Publication/1781 dimensional channel boarded by two solid walls. Xaman et al. T∞,2 [5] study laminar and turbulent natural convection combined T∞,1 with surface thermal radiation in a square cavity with a glass y wall. Zrikem et al. [6] numerically investigated coupled heat Tw Multi-layer door transfers by conduction, natural convection and radiation in x hollow structures heated from below or above. Oztop and co- Fig. 1 Model geometry The top and bottom sides of the multi-layer door were U. Prasopchingchana is with the Mechanical Engineering Department, assumed to be insulated. The multi-layer door of the furnace is Burapha University, Chonburi 20131 THAILAND (phone: 66-3810-2222 Ext. 50 cm. wide, 50 cm. high and 20.3 cm. thick. The multi-layer 3385 Ext. 104; fax: 66-3874-5806; e-mail: [email protected]). International Scholarly and Scientific Research & Innovation 5(1) 2011 229 scholar.waset.org/1307-6892/1781 World Academy of Science, Engineering and Technology International Journal of Mechanical and Mechatronics Engineering Vol:5, No:1, 2011 door of the furnace consists of 3 layers: fire clay brick, glass Radiation heat transfer from the furnace wall surfaces to the fiber insulation and steel sheet. The thermal conductivity inner door surface and from the outer door surface to the values [13] and the thickness values of each layer are shown surrounding area was included in the problem. The nonlinear in Table 1. equations of radiation heat transfer were linearized. Thus, the radiation heat transfer equations are TABLE I ⎛ ()k ⎞ qr,w = hr,1A⎜Tw − T ⎟ (6) PROPERTIES OF THE MULTI-LAYER DOOR ⎝ s,1 ⎠ Layer Material Thickness Thermal and (m) conductivity ⎛ ()k ⎞ (W/m⋅K) qr,sur = hr,2 A⎜T − Tsur ⎟ (7) ⎝ s,2 ⎠ 1 Fire clay brick 0.1 1.09 2 Glass fiber insulation 0.1 0.038 where hr,1 is the radiation heat transfer coefficient between the 3 Steel sheet 0.003 28 furnace wall surfaces and the furnace door surface, and hr,2 is the radiation heat transfer coefficient between the furnace door surface and the surrounding area. The radiation heat The emissivity values [13] at the surfaces of the multi-layer transfer coefficients are door inside and outside of the furnace are 0.75 and 0.066 ⎛ ()k −1 ⎞ 2 ()k −1 2 respectively. The temperatures at the furnace wall surfaces hr,1 = ε1σ ⎜Tw + T ⎟(Tw + []T ) (8) ⎝ s,1 ⎠ s,1 equaled the hot air temperatures in the cavity. Also, the and temperatures of the surrounding area equaled the room air 2 temperatures. ⎛ ()k −1 ⎞⎛ 2 ()k −1 ⎞ hr,2 = ε 2σ ⎜Tsur + T ⎟⎜Tsur + []Ts,2 ⎟ (9) ⎝ s,2 ⎠⎝ ⎠ III. MATHEMATICAL MODELS where ε is emissivity and σ is the Stefan-Boltzmann constant. Newton’s law of cooling was employed to determine the All properties of hot air and room air around the inside and the outside of the furnace door surfaces were obtained from heat transfer rates at the surfaces of the furnace door. To NIST (National Institute of Standards and Technology) determine the average Nusselt number, the work used the Standard Reference Database 23, Version 7.0 by polynomial Churchill and Chu empirical correlation [12]: curve fitting: 2 ⎧ 1/ 6 ⎫ 2 n−1 ⎪ 0.387Ra y ⎪ x = a1 + a2T f + a3T f + K + anT f (10) Nu y = ⎨0.825 + ⎬ (1) 9 /16 8/ 27 where x are the air properties, a are the polynomial curve ⎩⎪ []1+ ()0.492 / Pr ⎭⎪ n fitting coefficients and T is the film temperature of air. The where Ra is the Rayleigh number and Pr is the Prandtl f y errors of the properties derived from polynomial curve fitting number at the film temperatures. The film temperatures are are less than 0.1 percent. For all properties, n is 10. T f = (()Ts ave + T∞ )/ 2 (2) The volumetric thermal expansion coefficient is where ()Ts is the average temperature at the door surfaces 1 ⎛ ∂ ρ ⎞ ave β = - ⎜ ⎟ and T is the ambient temperature. The average convection ⎜ ⎟ ∞ ρ ⎝ ∂ T ⎠P heat transfer coefficients were obtained from 1 2 n−2 β = − (a2 + 2a3T f + 3a4T f + K + (n −1)anT ) (11) k f ρ f P h = Nu y (3) y y Within the thickness of the furnace door, the heat diffusion equation was used to calculate the temperature where kf is thermal conductivity of air. The local convection heat transfer coefficients could be obtained from distribution: ∂ 2T ∂ 2T 1 ⎛ y y+δy ⎞ + = 0 (12) h = ⎜ h δy + h δy⎟ 2 2 y+δy ⎜ ∫ y ∫ δy ⎟ ∂x ∂y y + δy ⎝ 0 y ⎠ The finite difference method was used to discretize Eq. (12) ⎛ 1 y ⎞ ⎛ 1 y+δy ⎞ and it yielded ()y + δy h = y⎜ ∫ h δy⎟ + δy⎜ ∫ h δy⎟ y+δy ⎜ y ⎟ ⎜ δy ⎟ ⎡ ⎤ ⎝ y 0 ⎠ ⎝ δy y ⎠ T ()k = Δy 2 ⎛T (k−1 )+ T ()k ⎞+ Δx 2 ⎛T (k−1 )+ T ()k ⎞ / International Science Index, Mechanical and Mechatronics Engineering Vol:5, No:1, 2011 waset.org/Publication/1781 ⎢ ⎜ ⎟ ⎜ ⎟⎥ i, j ⎣ ⎝ i+1, j i−1, j ⎠ ⎝ i, j+1 i, j−1 ⎠⎦ ()y + δy hy+δy = yhy + δyhδy 2()Δy 2 + Δx 2 (13) ()y + δy hy+δy − yhy The finite difference equation for the nodes close to the hδy = (4) δy interfaces as shown in Fig. 2 is Newton’s law of cooling was used to determine the convection heat transfer rates at the furnace door surfaces. The convection heat transfer rates at the furnace door surfaces are qδy = hδywδy(T∞ − Ts ) (5) where w is the furnace door width. International Scholarly and Scientific Research & Innovation 5(1) 2011 230 scholar.waset.org/1307-6892/1781 World Academy of Science, Engineering and Technology International Journal of Mechanical and Mechatronics Engineering Vol:5, No:1, 2011 ⎛ T (k) T (k−1) T (k) T (k−1) ⎞ ()k ⎡ kΔx ()k kΔx ()k−1 kΔy ()k−1 (k) ⎜ i−1, j i+1, j i, j−1 i, j+1 ⎟ T = ⎢ T + T + T Ti, j = ⎜ + + + ⎟ / 1, j ⎣2Δy 1,J −1 2Δy 1,J +1 Δx 2,J ⎜ R1 R2 R3 R3 ⎟ ⎝ ⎠ ⎤ + ()h1 j ΔyT∞,1 + ()hr,1 ΔyTw / ⎛ 1 1 1 1 ⎞ j ⎦⎥ ⎜ + + + ⎟ (14) ⎝ R1 R2 R3 R3 ⎠ ⎡kΔx kΔy ⎤ ⎢ + + ()h1 j Δy + ()hr,1 Δy⎥ (17) where R1, R2 and R3 are the thermal resistances.