Integrated Thermal Analysis of Natural Convection Air Cooled Electronic Enclosure

1999 ◽  
Vol 121 (2) ◽  
pp. 108-115 ◽  
Author(s):  
L. Tang ◽  
Y. K. Joshi

In the present paper, a methodology is described for the integrated thermal analysis of a laminar natural convection air cooled nonventilated electronic system. This approach is illustrated by modeling an enclosure with electronic components of different sizes mounted on a printed wiring board. First, a global model for the entire enclosure was developed using a finite volume computational fluid dynamics/heat transfer (CFD/CHT) approach on a coarse grid. Thermal information from the global model, in the form of board and component surface temperatures, local heat transfer coefficients and reference temperatures, and heat fluxes, was extracted. These quantities were interpolated on a finer grid using bilinear interpolation and further employed in board and component level thermal analyses as various boundary condition combinations. Thus, thermal analyses at all levels were connected. The component investigated is a leadless ceramic chip carrier (LCCC). The integrated analysis approach was validated by comparing the results for a LCCC package with those obtained from detailed system level thermal analysis for the same package. Two preferred boundary condition combinations are suggested for component level thermal analysis.

Author(s):  
Christoph Starke ◽  
Erik Janke ◽  
Toma´sˇ Hofer ◽  
Davide Lengani

Recent development in commercial CFD codes offers possibilities to include the solid body in order to perform conjugate heat transfer computations for complex geometries. The current paper aims to analyse the differences between a conjugate heat transfer computation and conventional uncoupled approaches where a heat transfer coefficient is first derived from a flow solution and then taken as boundary condition for a thermal conduction analysis of the solid part. Whereas the thermal analyses are done with a Rolls-Royce in-house finite element code, the CFD as well as the conjugate heat transfer computation are done using the new version 8 of the commercial code Fine Turbo from Numeca International. The analysed geometry is a turbine cascade that was tested by VKI in Brussels within the European FP6 project AITEB 2. First, the paper presents the aerodynamic results. The pure flow solutions are validated against pressure measurements of the cascade test. Then, the heat transfer from flow computations with wall temperature boundary conditions is compared to the measured heat transfer. Once validated, the heat transfer coefficients are used as boundary condition for three uncoupled thermal analyses of the blade to predict its surface temperatures in a steady state. The results are then compared to a conjugate heat transfer method. Therefore, a mesh of the solid blade was added to the validated flow computation. The paper will present and compare the results of conventional uncoupled thermal analyses with different strategies for the wall boundary condition to results of a conjugate heat transfer computation. As it turns out, the global results are similar but especially the over-tip region with its complex geometry and flow structure and where effective cooling is crucial shows remarkable differences because the conjugate heat transfer solution predicts lower blade tip temperatures. This will be explained by the missing coupling between the fluid and the solid domain.


1975 ◽  
Vol 97 (1) ◽  
pp. 47-53 ◽  
Author(s):  
R. E. Forbes ◽  
J. W. Cooper

Natural convection in horizontal layers of water cooled from above to near freezing was studied analytically. The water was confined laterally and underneath by rigid insulators, and the upper horizontal surface was subjected to: (1) a constant 0C temperature, rigid conducting boundary, and (2) a free, water to air convection boundary condition, in which the convective heat transfer coefficient was held constant at values of 5.68 W/m2 · K and 284 W/m2 · K (1.0 and 50.0 Btu/hr ft2F) and the temperature of the ambient air was maintained at 0C. The ratios of the width to the depth of the rectangular water layers under consideration were W/D = 1, 3, and 6. Initially the water is assumed to be at a uniform temperature of either 4C or 8C, and then the upper surface boundary condition was suddenly applied. It was observed in all cases for which the initial water temperature was 4C, that the water remained stagnant and became thermally stratified. Heat transfer application of either of the surface boundary conditions to water initially at 8C produced large convective eddies extending from the bottom to the top of the layer of water. As the liquid layer cooled further, two distinct horizontal regions appeared, the 4C isothermal line separating the two. This produces a region of hydrodynamic instability in the fluid since the maximum density fluid (4C) is physically located above the less dense fluid in the lower portion of the cavity. The large eddies which appeared initially were confined to the hydrodynamically unstable region bounded by the 4C isotherm and the bottom of the cavity. The action of viscous shearing forces upon the stable water above the 4C isotherm produced a second “layer” of eddies. An alternating direction implicit finite difference method was used to solve the coupled system of partial differential equations. The paper presents transient isotherms and streamlines and a discussion of the effect of maximum density on the flow patterns.


Author(s):  
Fakhreddine S. Oueslati ◽  
Rachid Bennacer ◽  
Habib Sammouda ◽  
Ali Belghith

The natural convection is studied in a cavity witch the lower half is filled with a porous media that is saturated with a first fluid (liquid), and the upper is filled with a second fluid (gas). The horizontal borders are heated and cooled by uniform heat fluxes and vertical ones are adiabatic. The formulation of the problem is based on the Darcy-Brinkman model. The density variation is taken into account by the Boussinesq approximation. The system of the coupled equations is resolved by the classic finite volume method. The numerical results show that the variation of the conductivity of the porous media influences strongly the flow structure and the heat transfer as well as in upper that in the lower zones. The effect of conductivity is conditioned by the porosity which plays a very significant roll on the heat transfer. The structures of this flow show that this kind of problem with specific boundary conditions generates a complex flow structure of several contra-rotating two to two cells, in the upper half of the cavity.


1991 ◽  
Vol 28 (3) ◽  
pp. 399-409 ◽  
Author(s):  
L. B. Smith ◽  
J. P. Graham ◽  
J. F. Nixon ◽  
A. S. Washuta

This paper presents a description of thermal analyses of forced-air ventilation and thermosyphon cooling systems, which were carried out in connection with the design of the concrete raft foundations that support hangars and other major structures to be constructed by the Department of National Defence adjacent to the existing airport near Inuvik, N.W.T. The cooling systems are required to prevent heat from the buildings from thawing the ice-rich permafrost present below the site. The analyses identified those parameters that have the most significant effect on the efficiency of each system. Based partially on the results of the analyses, it was decided to utilize air ventilation for cooling. The system is expected to perform satisfactorily under natural convection; however, the design includes a provision to install air blowers, if this should prove necessary in the future. A number of areas in which further research appears useful have been identified. Key words: permafrost, thermal analysis, raft foundation, hangar, ventilated slab, natural convection, thermosyphons.


Author(s):  
Jakob Hærvig ◽  
Anna Lyhne Jensen ◽  
Henrik Sørensen

Abstract Vertical smooth surfaces are commonly used for transferring heat by natural convection. Many studies have tried altering smooth surfaces in various ways to increase heat transfer. Many of these studies fail to increase global heat transfer. The problem commonly reported is dead zones appearing just upstream and downstream obstructions that effectively decrease wall temperature gradients normal to the surface. In this study, we simulate how changes geometry of forward facing triangular roughness elements affect local and global heat transfer for isothermal plates. We change the aspect ratio of the triangular elements from L/h = 5 to L/h = 40 at Grashof numbers of GrL = 8.0 · 104 and GrL = 6.4 · 105. In all cases the flow remains laminar. Even when accounting for the increase in surface area, we keep observing a decrease in global heat transfer compared to the smooth vertical plate. However, the results show by carefully selecting the aspect ratio and pitch distance of the triangular elements based on the Grashof number, the dead zone behind the horizontal part can be eliminated thereby significantly increasing local heat transfer. This observation could help to improve cooling of electronics with high localised heat fluxes.


Author(s):  
D. H. Greisen ◽  
V. P. Manno

Compact Thermal Models (CTMs) utilize a few connected thermal nodes to represent the thermal characteristics of electronic packages. These models are preferable to highly discretized models in preliminary design and system level analysis because of their computational efficiency. Surface heat flux non-uniformities often make it necessary to subdivide the package surfaces into multiple CTM nodes. This division is often quantified as the surface area ratio. This work assesses CTM performance sensitivity to area ratio changes and variation in heat transfer coefficient boundary conditions. CTMs for benchmark TQFP and BGA packages are developed using an admittance matrix approach. While optimum area ratios are identified, a direct correlation between these optimal values and the heat flux distributions computed from fully-discretized models was not obtained. CTM performance was found to be sensitive to changes in the heat transfer coefficient used to generate the CTM parameter values. A critical generating heat transfer coefficient was determined such that the resulting CTM, when optimized for a single boundary condition, was relatively accurate over the whole set of boundary conditions considered. This single boundary condition also provided an upper bound for error. This finding could be significant in future CTM development procedures.


Author(s):  
David M. Pratt ◽  
David J. Moorhouse

Aerospace vehicle design has progressed in an evolutionary manner, with certain discrete changes such as turbine engines replacing propellers for higher speeds. The evolution has worked very well for commercial aircraft because the major components can be optimized independently. This is not true for many military configurations which require a more integrated approach. In addition, the introduction of aspects for which there is no pre-existing database requires special attention. Examples of subsystem that have no pre-existing data base include directed energy weapons (DEW) such as high power microwaves (HPM) and high energy lasers (HEL). These devices are inefficient, therefore a large portion of the energy required to operate the device is converted to waste heat and must be transferred to a suitable heat sink. For HPM, the average heat load during one ‘shot’ is on the same order as traditional subsystems and thus designing a thermal management system is possible. The challenge is transferring the heat from the HPM device to a heat sink. The power density of each shot could be hundreds of megawatts. This heat must be transferred from the HPM beam dump to a sink. The heat transfer must occur at a rate that will support shots in the 10–100Hz range. For HEL systems, in addition to the high intensity, there are substantial system level thermal loads required to provide an ‘infinite magazine.’ Present models are inadequate to analyze these problems, current systems are unable to sustain the energy dissipation required and the high intensity heat fluxes applied over a very short duration phenomenon is not well understood. These are examples of potential future vehicle integration challenges. This paper addresses these and other subsystems integration challenges using a common currency for vehicle optimization. Exergy, entropy generation minimization, and energy optimization are examples of methodologies that can enable the creation of energy optimized systems. These approaches allow the manipulation of fundamental equations governing thermodynamics, heat transfer, and fluid mechanics to produce minimized irreversibilities at the vehicle, subsystem and device levels using a common currency. Applying these techniques to design for aircraft system-level energy efficiency would identify not only which subsystems are inefficient but also those that are close to their maximum theoretical efficiency while addressing diverse system interaction and optimal subsystem integration. Such analyses would obviously guide researchers and designers to the areas having the highest payoff and enable departures from the evolutionary process and create a breakthrough design.


Energies ◽  
2020 ◽  
Vol 13 (16) ◽  
pp. 4112
Author(s):  
Tian Zhao ◽  
Di Liu ◽  
Ke-Lun He ◽  
Xi Chen ◽  
Qun Chen

Optimization of heat transfer systems (HTSs) benefits energy efficiency. However, current optimization studies mainly focus on the improvement of system design, component design, and local process intensification separately, which may miss the optimal results and lack reliability. This work proposes a synergetic optimization method integrating levels of the local process, component to system, which could guarantee the reliability of results. The system-level optimization employs the heat current method and hydraulic analysis, the component level optimization adopts heuristic optimization algorithm, and the process level optimization applies the field synergy principle. The introduction of numerical simulation and iteration provides the self-consistency and credibility of results. Optimization results of a multi-loop heat transfer system present that the proposed method can save 16.3% pumping power consumption comparing to results only considering system and process level optimization. Moreover, the optimal parameters of component originate from the trade-off relation between two competing mechanisms of performance enhancement, i.e., the mass flow rate increase and shape variation. Finally, the proposed method is not limited to heat transfer systems but also applicable to other thermal systems.


Author(s):  
Koichi Hata ◽  
Katsuya Fukuda ◽  
Tohru Mizuuchi

Natural convection heat transfer from vertical 5×5 rod bundles in liquid sodium was numerically analyzed for two types of the bundle geometry (equilateral square and triangle arrays, ESA and ETA). The unsteady laminar three dimensional basic equations for natural convection heat transfer caused by a step heat flux were numerically solved until the solution reaches a steady-state. The PHOENICS code was used for the calculation considering the temperature dependence of thermophysical properties concerned. The 5×5 test rods for diameter (D = 7.6 mm), heated length (L = 200 mm) and L/d (= 26.32) were used in this work. The surface heat fluxes for each cylinder were equally given for a modified Rayleigh number, (Rf,L)ij and (Rf,L)5×5,S/D, ranging from 3.08 × 104 to 4.19 × 107 (q = 1 × 104∼7 × 106 W/m2) in liquid temperature (TL = 673.15 K). The values of S/D, which are ratios of the diameter of flow channel for bundle geometry to the rod diameter, for vertical 5×5 rod bundles were ranged from 1.8 to 6 on each bundle geometry. The spatial distribution of local and average Nusselt numbers, (Nuav)ij and (Nuav,B)5×5,S/D, on vertical rods of a bundle was clarified. The average value of Nusselt number, (Nuav)ij and (Nuav,B)5×5,S/D, for two types of the bundle geometry with various values of S/D were calculated to examine the effect of the bundle geometry, S/D, (Rf,L)ij and (Rf,L)5×5,S/D on heat transfer. The bundle geometry for the higher (Nuav,B)5×5,S/D value under the condition of S/D = constant was examined. The correlations for (Nuav,B)5×5,S/D for two types of bundle geometry above mentioned including the effects of (Rf,L)5×5,S/D and S/D were developed. The correlations can describe the theoretical values of (Nuav,B)5×5,S/D for two types of the bundle geometry for S/D ranging from 1.8 to 6 within −11.77 to 13.34 % difference.


Sign in / Sign up

Export Citation Format

Share Document