Study on the natural air cooling design of electronic equipment casings: Effects of the height and size of outlet vent on the flow resistances

Similar documents
Natural Convection from Horizontal Rectangular Fin Arrays within Perforated Chassis

Conjugate heat transfer from an electronic module package cooled by air in a rectangular duct

NUMERICAL SIMULATION OF CONJUGATE HEAT TRANSFER FROM MULTIPLE ELECTRONIC MODULE PACKAGES COOLED BY AIR

Experimental Analysis of Rectangular Fin Arrays with Continuous Fin and Interrupted Fins Using Natural and Forced Convection

Experimental Study of Convective Heat Transfer and Thermal Performance in the Heat-Sink Channel with Various Geometrical Configurations Fins

International Journal of Engineering Research and General Science Volume 3, Issue 6, November-December, 2015 ISSN

PHYSICAL MECHANISM OF NATURAL CONVECTION

Free and Forced Convection Heat Transfer Characteristics in an Opened Box with Parallel Heated Plates

CHME 302 CHEMICAL ENGINEERING LABOATORY-I EXPERIMENT 302-V FREE AND FORCED CONVECTION

Experiment 1. Measurement of Thermal Conductivity of a Metal (Brass) Bar

THERMO-FLOW CHARACTERISTICS OF A PIN-FIN RADIAL HEAT SINKS ACCORDING TO THEIR FIN HEIGHT PROFILE

EXPERIMENTAL AND NUMERICAL STUDIES OF A SPIRAL PLATE HEAT EXCHANGER

Experimental Analysis of Natural Convective Heat Transfer Performance From Rectangular Fin Array with Perforations

FLUID FLOW AND HEAT TRANSFER INVESTIGATION OF PERFORATED HEAT SINK UNDER MIXED CONVECTION 1 Mr. Shardul R Kulkarni, 2 Prof.S.Y.

Theoretical and Experimental Studies on Transient Heat Transfer for Forced Convection Flow of Helium Gas over a Horizontal Cylinder

EFFECT OF BAFFLES GEOMETRY ON HEAT TRANSFER ENHANCEMENT INSIDE CORRUGATED DUCT

Heat transfer coefficient of near boiling single phase flow with propane in horizontal circular micro channel

NATURAL CONVECTION HEAT TRANSFER CHARACTERISTICS OF KUR FUEL ASSEMBLY DURING LOSS OF COOLANT ACCIDENT

Heat Transfer Characteristics of Square Micro Pin Fins under Natural Convection

Experimental Study on Heat Transfer from a Perforated Fin Array with Cross Perforations

UNIT II CONVECTION HEAT TRANSFER

The Effect of Solid and Perforated Pin Fin on the Heat Transfer Performance of Finned Tube Heat Exchanger

Heat Transfer Enhancement Through Perforated Fin

Table of Contents. Foreword... xiii. Preface... xv

An experimental investigation on condensation of R134a refrigerant in microchannel heat exchanger

Keywords: Spiral plate heat exchanger, Heat transfer, Nusselt number

Ali S. Sharbuddin, International Journal of Advance Research, Ideas and Innovations in Technology.

IEEE TRANSACTIONS ON COMPONENTS, PACKAGING, AND MANUFACTURING TECHNOLOGY PART A, VOL. 20, NO. 4, DECEMBER

Analysis of the Cooling Design in Electrical Transformer

Maximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection

Heat Transfer Studies on Structured Metal Sheets

Computational Fluid Dynamics Based Analysis of Angled Rib Roughened Solar Air Heater Duct

CHAPTER-4 EVALUATION OF NANOFLUIDS PROPERTIES

Effect of roughness shape on heat transfer and flow friction characteristics of solar air heater with roughened absorber plate

Heat Transfer Enhancement of Solar Air Heater By Using Artificial Roughness double inclined ribs

Enhanced MicroChannel Heat Transfer in Macro- Geometry using Conventional Fabrication Approach

International Journal of Modern Trends in Engineering and Research e-issn No.: , Date: April, 2016

Microelectronics Heat Transfer Laboratory

Transitional Flow and Heat Transfer Characteristics in a Rectangular Duct with Stagger-arrayed Short Pin Fins

Available online at ScienceDirect. Procedia Engineering 90 (2014 )

Enhancement of Natural Convection Heat Transfer within Closed Enclosure Using Parallel Fins F. A. Gdhaidh, K. Hussain, H. S. Qi

COMPUTATIONAL ANALYSIS OF LAMINAR FORCED CONVECTION IN RECTANGULAR ENCLOSURES OF DIFFERENT ASPECT RATIOS

Experimental Analysis of Heat Transfer Enhancement over Dimpled Surface on One Side of Plate

Simple Analytical Calculations for

Review of Natural Convection from Rectangular Interrupted Fins.

Thermo-Fluid Performance of a Vapor- Chamber Finned Heat Sink

EFFECT OF THE INLET OPENING ON MIXED CONVECTION INSIDE A 3-D VENTILATED CAVITY

Natural and Mixed Convection Heat Transfer Cooling of Discrete Heat Sources Placed Near the Bottom on a PCB

Chapter 3 NATURAL CONVECTION

STUDY OF A PASSIVE SOLAR WINTER HEATING SYSTEM BASED ON TROMBE WALL

NUMERICAL HEAT TRANSFER ENHANCEMENT IN SQUARE DUCT WITH INTERNAL RIB

Constructal design: geometric optimization of micro-channel heat sinks

IJSRD - International Journal for Scientific Research & Development Vol. 3, Issue 06, 2015 ISSN (online):

Numerical Study of Conjugate Natural Convection Heat Transfer Using One Phase Liquid Cooling

CFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel

Investigation of Jet Impingement on Flat Plate Using Triangular and Trapezoid Vortex Generators

Flow and heat transfer characteristics of tornado-like vortex flow

HEAT TRANSFER CAPABILITY OF A THERMOSYPHON HEAT TRANSPORT DEVICE WITH EXPERIMENTAL AND CFD STUDIES

Enhancement of Heat Transfer Effectiveness of Plate-pin fin heat sinks With Central hole and Staggered positioning of Pin fins

Lecture 28. Key words: Heat transfer, conduction, convection, radiation, furnace, heat transfer coefficient

Introduction to Heat and Mass Transfer. Week 5

Heat Transfer Enhancement with Different Square Jagged Twisted Tapes and CuO Nano fluid

Numerical investigation to study effect of radiation on thermal performance of radiator for onan cooling configuration of transformer

HEAT TRANSFER THERMAL MANAGEMENT OF ELECTRONICS YOUNES SHABANY. C\ CRC Press W / Taylor Si Francis Group Boca Raton London New York

Atrium assisted natural ventilation of multi storey buildings

Thermo-Hydraulic performance of Internal finned tube Automobile Radiator

Experimental Investigations on the Local Distribution of wall static pressure coefficient Due To an Impinging Slot Air Jet on a Confined Rough Surface

The Electrodynamics of a Pair of PV Modules with Connected Building Resistance

Natural convection heat transfer around a horizontal circular cylinder near an isothermal vertical wall

of the heat is dissipated as a result of the flow of electrical current in various conductors. In order to predict temperatures

Finite-Element Evaluation of Thermal Response Tests Performed on U-Tube Borehole Heat Exchangers

Assessment of the Baffle Effects on the Mixed Convection in Open Cavity Nabil Jamil Yasin 1, Kadhum Audaa Jehhef 2, Asraa Mahdi Shaker 3 1,3

Experimental Analysis for Natural Convection Heat Transfer through Vertical Cylinder

World Journal of Engineering Research and Technology WJERT

Numerical Study of the Three-dimensional Flow in a Flat Plate Solar Collector with Baffles

Thermocouples calibration and analysis of the influence of the length of the sensor coating

Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions

A study of Heat Transfer Enhancement on a Tilted Rectangular Stainless Steel Plate

A Numerical Analysis of Indoor Thermal Environment and Human Thermophysiological Responses under Natural Ventilation S. Iizuka 1,*, T. Sakoi 2, T. Sai

Fluid Flow and Heat Transfer of Combined Forced-Natural Convection around Vertical Plate Placed in Vertical Downward Flow of Water

Forced Convection Heat Transfer Augmentation from Ribbed Surface

International Journal of Advanced Engineering Technology E-ISSN

Convection Heat Transfer. Introduction

Pin Fin Lab Report Example. Names. ME331 Lab

CFD study for cross flow heat exchanger with integral finned tube

Optimization of Plate Fin Heat Sinks Using Entropy Generation Minimization

Available online at ScienceDirect. International Conference On DESIGN AND MANUFACTURING, IConDM 2013

Natural Convection Heat Loss from A Partly Open Cubic Enclosure Timothy N Anderson 1,a * and Stuart E Norris 2,b

Numerical investigation of entropy generation in microchannels heat sink with different shapes

"Experimental investigation of heat transfer enhancement by u-shaped Turbulator"

Piping Systems and Flow Analysis (Chapter 3)

Heat Transfer Convection

Thermo Hydraulic Performance of Solar Air Heater by Using Double Inclined Discrete Rib Roughened Absorber Plate

CFD AND CONJUGATE HEAT TRANSFER ANALYSIS OF HEAT SINKS WITH DIFFERENT FIN GEOMETRIES SUBJECTED TO FORCED CONVECTION USED IN ELECTRONICS COOLING

Numerical Investigation of Convective Heat Transfer in Pin Fin Type Heat Sink used for Led Application by using CFD

W-Discrete Rib for Enhancing the Thermal Performance of Solar Air Heater

FLOW BOILING HEAT-TRANSFER IN PLATE MICRO- CHANNEL HEAT SINK

Experimental Study of Natural Convective Heat Transfer from Vertical Rectangular Fin Array at Different Angle of Inclination

Mathematical Calculations Of Heat Transfer For The CNC Deposition Platform Based On Chemical Thermal Method

NUMERICAL STUDY OF MIXED CONVECTION AND THERMAL RADIATION IN A SQUARE CAVITY WITH AN INSIDE INCLINED HEATER

Transcription:

Journal of Physics: Conference Series Study on the natural air cooling design of electronic equipment casings: Effects of the height and size of outlet vent on the flow resistances To cite this article: M Ishizuka et al 2012 J. Phys.: Conf. Ser. 395 012122 View the article online for updates and enhancements. Related content - A Test Bench for Validation of CFD Analysis for Air Flow in Thin Electronic Casing Model Tomoyuki Hatakeyama, Masaru Ishizuka, Risako Kibushi et al. - First average and local heat transfer measurements on a forced air-flow at low Re-numbers through a rectangular channel with ribbed surfaces D Fustinoni, P Gramazio, L P M Colombo et al. - Thermal management of electronics using phase change material based pin fin heat sinks R Baby and C Balaji This content was downloaded from IP address 46.3.193.154 on 25/02/2018 at 05:57

Study on the natural air cooling design of electronic equipment casings: Effects of the height and size of outlet vent on the flow resistances M Ishizuka 1, T Hatakeyama, R Kibushi and M. Inoue Department of Mechanical Systems Engineering, Toyama Prefectural University, 5180 Kurokawa, Imizu, Toyama 939-0398, Japan E-mail: ishizuka@pu-toyama.ac.jp Abstract. This paper describes the effects of the outlet vent size and the distance between the outlet vent location and the power heater position on the flow resistance in natural-air-cooled electronic equipment casings. An experiment was carried out using a simple model casing simulated for the practical natural-air-cooled casing which is composed of 4 side walls, a top plate and bottom plate which has an inlet opening. A power heater to served as a power dissipation unit was placed at its open bottom. An outlet opening was set on one of the side walls. The opening area, the height of the outlet and the heater location were varied. The experimental results were analyzed using the flow resistance coefficient K which was related to the distance between the outlet vent and the power heater position and the heat removal from the outlet vent, and K values were plotted against a pair of Reynolds numbers Re and the outlet vent porosity β which is defined as the ratio of outlet vent open area to the top surface area of the casing. 1. Introduction The power dissipation density of electronic equipment has increased in recent years and, as a result, it is necessary to consider the cooling design of electronic equipment in order to develop suitable cooling techniques. The literature [1] includes a number of studies on air cooling, which is an important subject because almost all electronic equipment are cooled by air convection. Of the various cooling systems of electronic equipment, natural air cooling is often utilized for those applications, such as telecommunication, for which high reliability is essential. The main advantage of natural convection is that no fan or blower is required, because air movement is simply generated by density differences in the presence of a gravity field. The optimum thermal design of electronic devices, cooled by natural convection, depends on an accurate choice of geometrical configuration and on heat sources distribution to be able to promote the thermo-circulation flow rate that minimizes the temperature rise (mean and maximum values) inside the casings. Guglielmini et al. [2] have reported on the natural air cooling of electronic cards in ventilated enclosures. Ishizuka et al.[3] presented a simplified set of equations which represents the cooling capability on the basis of data on natural air cooling of electronic equipment casing. However, there are insufficient data for detailed discussion regarding the design of practical electronic equipment. For example, the simplified set of equations was introduced using a ventilation model like a chimney Published under licence by Ltd 1

which was formed by a heater at the duct bottom and an outlet vent on the top. However, in practical electronic equipment, the outlet vent is located at the upper part of the side walls, not on the top, and a casing is not composed of a circular duct. Therefore, in the present work, an attempt is made to study the effects of the distance between the outlet vent location and the power heater position on the cooling capability of natural-air-cooled electronic equipment casings. 2. Set of equations Ishizuka et al.[3] proposed the following set of equations for engineering applications in the thermal design of electronic equipment o 1.78 1.25 0.5 1.5 eq. m 300 o ( / ) o ( 1) 2.5(1 1.3 2 β ) / β Tm ( 2) ( 3) = + + 1/ 2 Q = S T + A h K T K = T = S S S S eq. top side bottom 4 where, Q denotes the total heat generated by the components, S eq. is the equivalent total source of the casing, T m is the average temperature rise in the casing, A o is the outlet vent area of the casing, h is the distance from the heater position to the outlet, K is the flow resistance coefficient arising from the air path blockage at the outlet, T o is the air temperature rise at the outlet vent and β is the porosity coefficient of the outlet vent. The flow resistance coefficient K was obtained approximately as a function of porosity coefficient of the outlet vent in the casing. The following relations for wire nets have been obtained in a low Reynolds number range by Ishizuka[4]. K = Re β β 2 0.95 40( (1 ) / ) 5 where Re is defined on the basis of the wire diameter used in the wire nets. However, since the effect of the Reynolds number on K is less prominent than that of the porosity coefficient β, K can be approximated as a function of the porosity coefficient only. Therefore, Eq.(2) is considered to be a reasonable expression for practical applications. The h value in Eq.(1) is called chimney height and Eq.(1) was introduced assuming a ventilation model like a chimney which is composed of a heater at the duct bottom and outlet vent on the top as shown in Fig.1. However, in practical electronic equipment, the outlet is located at the upper part of the side walls, not on the top, and casings are not circular ducts. Following experiments were carried out by taking those factors into account. Figure 1 Ventilation model 3. Experiments 3.1Experimental apparatus 2

The experimental apparatus is shown in Fig.2 and consists of a casing having the following external dimensions: length x width x height = 220 x 230 x 310 mm 3. The casing surfaces were made of plastic whose thickness and thermal conductivity were 10 mm and 0.01 W/(mK), respectively. A wire heater which is shown in Fig.3, served as heat dissipation unit and is placed at a position 25 mm above the bottom of the casing. The diameter of the heater was 1.6 mm. On one of the side walls of the box, one rectangular shaped opening is left and the size and location of which were varied. The opening dimensions have six patterns and are showed in Fig.4. The experiments were performed for three different distances between the wall bottom and the center of the opening of the outlet vent (H v = 275, 200, and 150 mm) and for three different vertical locations of the heater from the bottom of the wall( H h = 25, 125, and 225 mm). However, the experiments were carried out in such a way that the heater was always placed below the outlet vent. The cooling air enters through the casing bottom and is exhausted through the upper outlet opening. There is an inlet vent whose area is 150 mm x 130 mm in the central region of the casing bottom. The air temperature distribution inside the casing, the room temperature and the wall temperature have been measured by means of 14 calibrated K-type thermocouples (±0.1K and 0.1 mm in diameter). Inside the casing, they are arranged at the different horizontal locations ( 30, 115 and 115 mm) from the inside of the left side wall for each three different height values from the wall bottom (75, 175, and 275 mm) as shown in Fig.2. Figure 2 Experimental Casing Figure 3 Heater shape On the walls, 4 thermocouples were placed at the center of the top surface and at the 3 locations of the central line of each of the side walls, and one thermocouple was placed outside the casing for the measurement of room temperature. The mean inner air temperature rise T m has been calculated from the locally measured values over the whole volume of the casing. 3.2 Estimation of heat removed from casing surfaces At first, the amount of heat removed from the casing surfaces was estimated from the experiment. For the purpose of this estimation, the casing in Fig.2 was considered to be closed unit casing with no vents except the open bottom. That is, it consisted of a top and four side walls. Using the natural convective heat transfer equations for individual surfaces presented by Ishizuka et al[4], the amount of heat removed from the casing surfaces, Qs, was expressed as follows. 1.25 Q = D T s 1 m 6 3

Figure 4 Outlet vent Openings The Eq.(6) was introduced in the same way as that of the first term on the right hand side of Eq.(1).In Eq.(5), constant D 1 includes a radiative heat transfer factor as well. The amount of heat generation was within the Q s = 6 W- 40 W range, and the room temperature T a was 298K. Results obtained from this experiment are presented in Fig.5. It is found from this figure that the results are well-expressed by the following Eq.(7),in which the arithmetic mean temperature rise T m was obtained from measuring the temperature of the 9 internal locations in the casing. 1.25 Q = 0.445 T 7 s m where, the constant D 1 was determined as 0.445 to meet the measured values. Temperature distribution in the casing was relatively uniform and the maximum ratio of the temperature difference value to the temperature in the casing was 5 %. Hereafter, the amount of heat removed from the outlet vent Q v was calculated on the basis of the following equation. Q = Q Q v s 8 3.3 Influence of outlet vent size on the temperature rise in the casing The experiment was carried out using a reference casing in Fig.2 by varying the outlet vent size. The reference casing was defined as that outlet vent position was at H v =275 mm and the heater position was at H h =25 mm. Porosity coefficient βo, presented in Fig.4, is defined as the ratio of the open area of each individual vent to the area of the reference vent (150 mm x 50 mm). The results are shown in Fig.6. 4

Figure 5 Relationship between Q s and T m Figure 6 Influence of vent porosity on the mean temperature rise in the casing In Fig.6, mean temperature rise T m in the casing is logarithmically plotted against the porosity coefficient of the openings, β o in log scale by varying the input power, Q value. As the β o decreases, T m increases linearly in the logarithm plot. T m also increases, as the input power increases. 3.4 Influence of outlet vent position on the mean temperature rise in the casing The relationship between internal temperature rise T m and outlet vent position H was investigated by varying opening size under the condition that the heater position was fixed at the bottom. Fig.7 shows the experimental results. The outlet vent height was varied by preparing 3 kinds of side walls with 3 different outlet vent height values. T m decreases as H value increases. However, the slope which was a little larger in the case of small H value, becomes smaller in the range of H>200 mm. The relation between T m and H for two opening size values exhibited a similar pattern. 3.5 Influence of the distance between the outlet vent position and the heater position on the mean temperature rise in the casing The relationship between internal temperature rise T m and the distance between the outlet vent position and the heater position h was investigated by varying opening size and input power value under the condition that the outlet was fixed at H v =275 mm. Fig.8 shows the experimental results. As h value increases, T m decreases. T m was defined as the average of the temperature values measured at the locations above the heater position in the casing. That is, the temperature values at the temperature locations under the heater position were not considered. The slopes of three point groups are a little different from each other. 5

Figure 7 Influence of outlet vent position on the mean temperature rise in the casing 4. Correlation of the data using non-dimensional parameters 4.1 Flow Resistance Coefficient K The flow resistance coefficient K was used to correlate the data obtained in the experiments, since the K value was related to Q v and the distance h as follows. A ventilation model, as shown in Fig.1, is considered. Assuming a uniform temperature distribution and a one-dimensional steady state flow, the two expressions below can be written, one for the overall energy balance and the other for the balance between flow resistance and buoyancy force. Q = ρc Au T 9 v p gh= K u 2 ( ρa ρ) ρ / 2 10 where, Q is dissipated power, ρ is air density, c p is specific heat of the air at constant pressure, u is airflow velocity, A is the cross-sectional area of the duct, T is temperature rise, g is acceleration due to gravity, h is the distance between the outlet and the heater and K is the flow resistance coefficient for the system. The subscript, a, is used to represent the atmosphere condition. Since the pressure change in the system is small, the expression is obtained assuming a perfect gas. ( ρ ρ) / ρ = T T / T 11 a a a and the equation for K is obtained as follows, 3 2 K = 2g h T / ( T ( ρc A / Q ) ) 12 a p v Eq.(12) indicates the relationship among K, h instead of H in Fig.7 and Q v. Figure 8 Influence of the distance between the outlet vent position and the heater position on the mean temperature rise in the casing 4.2 Porosity coefficient β In Fig.3, although the porosity coefficient of opening β o was defined as the ratio of the area of the opening to the reference opening, its value was generally unknown, since the data on reference opening was generally unavailable.. Thus, the porosity coefficient β is defined as follows in this study. β = open area in the outlet vent / inner casing top surface area. (13) 6

4.3 Reynolds number Re The Reynolds number was also used in order to correlate the data. The velocity u was obtained from Eq.(9) and the hydrodynamic equivalent diameter of the openings was used as a reference length as follows with the height A and the width B of the opening. L = 4 A B 2 A + B 14 Thus, the Reynolds number Re is defined as follows. Re = u L / ν 15 4.4 Relationship among K, Re and β In Figure 9, the variation of K with X is plotted. The parameter X was defined on the basis of the method of correlation of the K values for wire nets and perforated plates reported the author et al[3,-4]. X = Re β 2 / (1 β ) 16 Almost all the K values lie on the line obtained by the empirical correlation. The K values obtained in the case of the reference condition that the outlet vent was at the upper position and the heater was kept at the bottom, lie on the line plotted by the empirical correlation. However, the K values in the non-reference case indicate slightly smaller values. This reason behind this discrepancy is expected to be due to measurement uncertainty, other than the change of basic phenomena. And the error due to measurement uncertainty composed of mean temperature rise values. The line was obtained from the best fit to the measure data as follows. 1.5 K = 0.4 X 17 where, the coefficient having a value of 0.4 is inherent in this apparatus and not a general value. This indicates that the distance between the outlet vent and the heat dissipation unit position h can be considered as a chimney height in practical equipment as well and the definition of β is reasonable. A more detailed discussion requires 3-dimensional thermo-fluid analysis and more precise measurement. However, it is considered that a useful relationship among the flow resistance coefficient, Reynolds number and porosity coefficient can be clarified from a practical point of view. Figure 9 Relationship among K, Re and β 7

5. Conclusion The effects of the size of outlet vent opening and the distance between the outlet vent and the power heater location on the flow resistance in natural air cooled electronic equipment casings were studied experimentally. From the results, the relationship among flow resistance coefficient K, Reynolds number Re and the outlet vent porosity coefficient β defined on the basis of the top surface are obtained as follows. 1.5 K = B X X = Re β 2 / (1 β ) B is the inherent coefficient of the casing in question. 6. References [1]Bergles AE,1990, Heat Transfer in Electronic and Microelectronic Equipment, Hemisphere. New York. [2] Guglielmini D., Milano G., and Misale M., 1988, Second UK National Conference on Heat Transfer,.1 199. [3] Ishizuka M., 1995, ASME.HTD-Vol.303, National Heat Transfer Conference-1 65. [4] Ishizuka M., Miyazaki Y. and Sasaki T., 1987, ASME J.of Heat Transfer,.109 540 8