Transient Thermal Management of Microelectronics Using Free Liquid Jet Impingement

Author(s):  
Antonio J. Bula ◽  
Muhammad M. Rahman

The results of numerical simulation of a transient heat transfer process when a free jet of high Prandtl number fluid impinges perpendicularly on a solid substrate of finite thickness containing discrete electronics on the opposite surface are presented. The numerical model was developed considering both solid and fluid regions and solved as a conjugate problem. Equations for the conservation of mass, momentum, and energy were solved in the liquid region taking into account the transport processes at the inlet and exit boundaries as well as at the solid-liquid and liquid-gas interfaces. In the solid region, only heat conduction equation was solved. The shape and location of the free surface (liquid-gas interface) was determined iteratively as a part of the solution process by satisfying the kinematic condition as well as the balance of normal and shear forces at this interface. The number of elements in the fluid and solid regions were determined from a systematic grid-independence study. A non-uniform grid distribution was used to adequately capture large variations near the solid-fluid interface. Computed results included the local and average heat transfer coefficients at the solid-fluid interface. Computations were carried out to investigate the influence of different operating parameters such as jet velocity and plate material. It was found that the average heat transfer coefficient is maximum at early stages of the transient process and decreases gradually with time to the final steady state condition.

2000 ◽  
Author(s):  
Antonio J. Bula ◽  
Muhammad M. Rahman ◽  
John E. Leland

Abstract Transient conjugate heat transfer process during axial free jet impingement on a solid disk of finite thickness was considered. As the fluid reached steady state, power was turned on and a uniform heat flux was imposed on the disk at its opposite surface. The numerical model considered both solid and fluid regions. Equations for conservation of mass, momentum, and energy were solved in the liquid region taking into account the transport processes at the inlet and exit boundaries, as well as at the solid-liquid and liquid-gas interfaces. Inside the solid, only the heat conduction equation was solved. The shape and location of the free surface (liquid-gas interface) was determined iteratively as a part of the solution process by satisfying the kinematic condition as well as the balance of normal and shear forces at this interface. A non-uniform grid distribution, captured from a systematic grid-independence study, was used to adequately accommodate large variations near the solid-fluid interface. Computed results include the simulation of six different substrate materials namely, aluminum, constantan, copper, diamond, silicon, and silver, and three different impinging liquids, FC - 77, Mil - 7808, and water. The solids and fluids selected covered a wide range of possibilities of conjugate heat transfer phenomena. The analysis performed showed that the thermal storage capacity, defined as density times specific heat, is an important factor defining which material will attain steady state faster during conjugate heat transfer process, like the thermal diffusivity does it for pure conduction heat transfer.


2010 ◽  
Vol 132 (5) ◽  
Author(s):  
Jorge C. Lallave ◽  
Muhammad M. Rahman

Abstract This paper considers the transient conjugate heat transfer characterization of a partially confined liquid jet impinging on a rotating and uniformly heated solid disk of finite thickness and radius. A constant heat flux was imposed at the bottom surface of the solid disk at t=0, and heat transfer was monitored for the entire duration of the transient until the steady state condition was reached. Calculations were done for a number of disk materials using water as the coolant, covering a range of Reynolds numbers (225–900), Ekman numbers (7.08×10−5−∞), nozzle-to-target spacing (β=0.25–1.0), confinement ratios (rp/rd=0.2–0.75), disk thicknesses to nozzle diameter ratios (b/dn=0.25–1.67), and solid to fluid thermal conductivity ratios (36.91–697.56). It was found that a higher Reynolds number decreases the time to achieve the steady state condition and increases the local and average Nusselt number. The duration of the transient increases with the increment of the Ekman number and disk thickness, and the reduction in the thermal diffusivity of the disk material.


Author(s):  
Muhammad M. Rahman ◽  
Cesar F. Hernandez

Transient heating of a hemispherical solid plate of finite thickness during impingement of a free liquid jet is studied. A constant heat flux is imposed at the inner surface of the hemispherical plate at t = 0 and heat transfer is monitored for the entire duration of the transient until a steady state condition is reached. Calculations are done for Re = 500–1500 and b/dn = 0.083–1.5 using water (H2O) as the coolant and various solid materials such as silicon, Constantan, and copper. It was found that the time for the plate to achieve the steady-state condition decreases and Nusselt number increases with Reynolds number. A plate material with higher thickness provides higher average Nusselt number and longer transient period.


Energies ◽  
2021 ◽  
Vol 14 (14) ◽  
pp. 4327
Author(s):  
Min-Seob Shin ◽  
Santhosh Senguttuvan ◽  
Sung-Min Kim

The present study experimentally and numerically investigates the effect of channel height on the flow and heat transfer characteristics of a channel impingement cooling configuration for various jet Reynolds numbers in the range of 2000–8600. A single array consisting of eleven jets with 0.8 mm diameter injects water into the channel with 2 mm width at four different channel heights (3, 4, 5, and 6 mm). The average heat transfer coefficients at the target surface are measured by maintaining a temperature difference between the jet exit and the target surface in the range of 15–17 °C for each channel height. The experimental results show the average heat transfer coefficient at the target surface increases with the jet Reynolds number and decreases with the channel height. An average Nusselt number correlation is developed based on 85 experimentally measured data points with a mean absolute error of less than 4.31%. The numerical simulation accurately predicts the overall heat transfer rate within 10% error. The numerical results are analyzed to investigate the flow structure and its effect on the local heat transfer characteristics. The present study advances the primary understanding of the flow and heat transfer characteristics of the channel impingement cooling configuration with liquid jets.


Author(s):  
Fatih Selimefendigil ◽  
Hakan F. Oztop

Abstract Nanoliquid impingement heat transfer with phase change material (PCM) installed radial system is considered. Study is performed by using finite element method for various values of Reynolds numbers (100 ≤ Re ≤ 300), height of PCM (0.25H ≤ hpcm = 0.7H ≤ 0.75H) and plate spacing (0.15H ≤ hpcm = 0.7H ≤ 0.40H). Different configurations with using water, nanoliquid and nanoliquid+PCM are compared in terms of heat transfer improvement. Thermal performance is improved by using PCM while best performance is achieved with nanoliquid and PCM installed configuration. At Re=100 and Re=300, heat transfer improvements of 26% and 25.5% are achieved with nanoliquid+PCM system as compared to water without PCM. Height of the PCM layer also influences the heat transfer dynamic behavior while there is 12.6% variation in the spatial average heat transfer of the target surface with the lowest and highest PCM height while discharging time increases by about 76.5%. As the spacing between the plates decreases, average heat transfer rises and there is 38% variation.


1996 ◽  
Vol 118 (1) ◽  
pp. 21-26 ◽  
Author(s):  
David Copeland

Experimental measurements of multiple nozzle submerged jet array impingement single-phase and boiling heat transfer were made using FC-72 and 1 cm square copper pin fin arrays, having equal width and spacing of 0.1 and 0.2 mm, with aspect ratios from 1 to 5. Arrays of 25 and 100 nozzles were used, with diameters of 0.25 to 1.0 mm providing nozzle area from 5 to 20 mm2 (5 to 20% of the heat source base area). Flow rates of 2.5 to 10 cm3/s (0.15 to 0.6 l/min) were studied, with nozzle velocities from 0.125 to 2 m/s. Single nozzles and smooth surfaces were also evaluated for comparison. Single-phase heat transfer coefficients (based on planform area) from 2.4 to 49.3 kW/m2 K were measured, while critical heat flux varied from 45 to 395 W/cm2. Correlations of the single-phase heat transfer coefficient and critical heat flux as functions of pin fin dimensions, number of nozzles, nozzle area and liquid flow rate are provided.


Author(s):  
D. E. Metzger ◽  
R. S. Bunker ◽  
G. Bosch

All experimental technique has been developed for measurement of local convection heat transfer characteristics on rotating surfaces, utilizing thin liquid crystal surface coatings in a thermal transient test procedure. The encapsulated liquid crystal coatings used are sprayed directly on the test surface and their response is observed and processed during the transient with automated computer vision and data acquisition systems. Heat transfer coefficients are calculated from the thermal transient response of the test surface, as determined from the color indication from the thin coating. A significant advantage of the method, especially for convection in disk/shroud cavities that may contain recirculating fluid regions, is that appropriate thermal boundary conditions are naturally imposed on all of the boundary surfaces. The method is also relatively fast and inexpensive, and allows the geometry of the disk and stator surfaces to be changed easily, without the expenses of mounting discrete heat flux and temperature sensors and equipment to transfer information to a stationary frame of reference. Application of the experimental technique is demonstrated with detailed radially local surface Nusselt number distributions acquired for cases involving jet impingement onto a plane smooth disk, rotating in close proximity to a parallel plane stator disk. A single circular jet, with nozzle exit flush mounted in the stator, is oriented normal to the disk surface at various radii and flowrates. Local Nusselt numbers are presented nondimensionally as functions of both disk and flow Reynolds numbers. The results indicate that the local radial heat transfer distribution can be controlled by varying the impingement radius, but maximum radially averaged heat transfer is obtained with impingement at the disk center.


2000 ◽  
Author(s):  
Qiao Lin ◽  
Shuyun Wu ◽  
Yin Yuen ◽  
Yu-Chong Tai ◽  
Chin-Ming Ho

Abstract This paper presents an experimental investigation on MEMS impinging jets as applied to micro heat exchangers. We have fabricated MEMS single and array jet nozzles using DRIE technology, as well as a MEMS quartz chip providing a simulated hot surface for jet impingement. The quartz chip, with an integrated polysilicon thin-film heater and distributed temperature sensors, offers high spatial resolution in temperature measurement due to the low thermal conductivity of quartz. From measured temperature distributions, heat transfer coefficients are computed for single and array micro impinging jets using finite element analysis. The results from this study for the first time provide extensive data on spatial distributions of micro impinging-jet heat transfer coefficients, and demonstrate the viability of MEMS heat exchangers that use micro impinging jets.


1979 ◽  
Vol 101 (1) ◽  
pp. 120-125 ◽  
Author(s):  
K. R. Randall ◽  
J. W. Mitchell ◽  
M. M. El-Wakil

Heat transfer by natural convection in rectangular enclosures has been experimentally studied using interferometric techniques. The effects of Grashof number, tilt angle, and aspect ratio on both the local and average heat transfer coefficients have been determined. The Grashof number range tested was 4 × 103 to 3.1 × 105, and the aspect ratio (ratio of enclosure length to plate spacing) varied between 9 and 36. The angles of tilt of the enclosure with respect to the horizontal were 45, 60, 75 and 90 deg. Correlations are developed for both local and average Nusselt number over the range of test variables. The effect of tilt angle is found to reduce the average heat transfer by about 18 percent from the value of 45 deg to that at 90 deg. No significant effect of aspect ratio over the range tested was found. A method for characterizing the flow regimes that is based on heat transfer mechanisms is proposed.


1997 ◽  
Vol 119 (2) ◽  
pp. 381-389 ◽  
Author(s):  
M. E. Taslim ◽  
C. M. Wadsworth

Turbine blade cooling, a common practice in modern aircraft engines, is accomplished, among other methods, by passing the cooling air through an often serpentine passage in the core of the blade. Furthermore, to enhance the heat transfer coefficient, these passages are roughened with rib-shaped turbulence promoters (turbulators). Considerable data are available on the heat transfer coefficient on the passage surface between the ribs. However, the heat transfer coefficients on the surface of the ribs themselves have not been investigated to the same extent. In small aircraft engines with small cooling passages and relatively large ribs, the rib surfaces comprise a large portion of the passage heat transfer area. Therefore, an accurate account of the heat transfer coefficient on the rib surfaces is critical in the overall design of the blade cooling system. The objective of this experimental investigation was to conduct a series of 13 tests to measure the rib surface-averaged heat transfer coefficient, hrib, in a square duct roughened with staggered 90 deg ribs. To investigate the effects that blockage ratio, e/Dh and pitch-to-height ratio, S/e, have on hrib and passage friction factor, three rib geometries corresponding to blockage ratios of 0.133, 0.167, and 0.25 were tested for pitch-to-height ratios of 5, 7, 8.5, and 10. Comparisons were made between the rib average heat transfer coefficient and that on the wall surface between two ribs, hfloor, reported previously. Heat transfer coefficients of the upstream-most rib and that of a typical rib located in the middle of the rib-roughened region of the passage wall were also compared. It is concluded that: 1 The rib average heat transfer coefficient is much higher than that for the area between the ribs; 2 similar to the heat transfer coefficient on the surface between the ribs, the average rib heat transfer coefficient increases with the blockage ratio; 3 a pitch-to-height ratios of 8.5 consistently produced the highest rib average heat transfer coefficients amongst all tested; 4 under otherwise identical conditions, ribs in upstream-most position produced lower heat transfer coefficients than the midchannel positions, 5 the upstream-most rib average heat transfer coefficients decreased with the blockage ratio; and 6 thermal performance decreased with increased blockage ratio. While a pitch-to-height ratio of 8.5 and 10 had the highest thermal performance for the smallest rib geometry, thermal performance of high blockage ribs did not change significantly with the pitch-to-height ratio.


Sign in / Sign up

Export Citation Format

Share Document