Transient Data Processing of Flow Boiling Local Heat Transfer in a Multi-Microchannel Evaporator Under a Heat Flux Disturbance

Multi-microchannel evaporators are often used to cool down electronic devices subjected to continuous heat load variations. However, so far, rare studies have addressed the transient flow boiling local heat transfer data occurring in such applications. The present paper introduces and compares two different data reduction methods for transient flow boiling data in a multi-microchannel evaporator. A transient test of heat disturbance from 20 to 30W cm(-2) was conducted in a multi-microchannel evaporator using R236fa as the test fluid. The test section was 1 x 1 cm(2) in size and had 67 channels, each having a cross-sectional area of 100 x 100 mu m(2). The micro-evaporator backside temperature was obtained with a fine-resolution infrared (IR) camera. The first data reduction method ( referred to three-dimensional (3D)-TDMA) consists in solving a transient 3D inverse heat conduction problem by using a tridiagonal matrix algorithm ( TDMA), a Newton-Raphson iteration, and a local energy balance method. The second method ( referred to two-dimensional (2D)controlled) considers only 2D conduction in the substrate of the micro-evaporator and solves at each time step the well-posed 2D conduction problem using a semi-implicit solver. It is shown that the first method is more accurate, while the second one reduces significantly the computational time but led to an approximated solution. This is mainly due to the 2D assumption used in the second method without considering heat conduction in the widthwise direction of the micro-evaporator.

Numerical Modeling of Annular Laminar Film Condensation in Circular and Non-Circular Micro-Channels under Normal and Micro-Gravity

Stefano Nebuloni

A theoretical and numerical model to predict film condensation heat transfer in mini, micro and ultra micro-channels of different internal shapes is presented in this thesis. The model is based on a finite volume formulation of the Navier-Stokes and energy equations and it includes the contributions of the unsteady terms, surface tension, axial shear stresses, gravitational forces and wall thermal conduction. Notably, interphase mass transfer and near-to-wall effects (disjoining pressure) are also included. This model has been validated versus various benchmark cases and versus published experimental results from three different laboratories, predicting micro-channel heat transfer data with an average error of 20 % or better. The conjugate heat transfer problem arising from the coupling between the thin film fluid dynamics, the heat transfer in the condensing fluid and the heat conduction in the channel wall has been studied and analyzed. The work has focused on the effects of three external wall boundary conditions: a uniform wall temperature, a non uniform wall heat flux and single-phase convective cooling. The thermal axial and peripheral conduction occurring in the wall of the channel can affect the behavior of the condensate film, not only because it redistributes the heat, but also because the annular laminar film condensation process is dependent on the local saturation to wall temperature difference. When moving from mini to micro and ultra-micro channels, the results shows that the axial conduction effects can become very important in the prediction of the wall temperature profile and they can not be ignored. Under these conditions, the overall performances of the heat exchanger become dependent not only on the fluid properties and the operative conditions but also on the geometry and wall material. Results obtained for steady state conditions are presented for circular, elliptical and flattened shape cross sections for R-134a and ammonia, for hydraulic diameters between 10 µm and 3 mm. Microscale condensation finds applications in heat pipes and compact heat exchangers for electronic equipment or spacecraft thermal control, in automotive condensers, in residential air conditioning and in refrigeration applications: the influence of the steady or unsteady gravitational field and of the inertia forces on the flow field and consequently on the heat transfer performances is investigated allowing the model to be applied as a design and optimization tool for enhanced heat exchangers.

EPFL2010

Cooling High Heat Flux Micro-Electronic Systems using Refrigerants in High Aspect Ratio Multi-Microchannel Evaporators

Etienne Costa-Patry

Improving the energy efficiency of cooling systems can contribute to reduce the emission of greenhouse gases. Currently, most microelectronic applications are air-cooled. Switching to two-phase cooling systems would decrease power consumption and allow for the reuse of the extracted heat. For this type of application, multi-microchannel evaporators are thought to be well adapted. However, such devices have not been tested for a wide range of operating conditions, such that their thermal response to the high non-uniform power map typically generated by microelectronics has not been studied. This research project aims at clarifying these gray areas by investigating the behavior of the two-phase flow of different refrigerants in silicon and copper multi-microchannel evaporators under uniform, non-uniform and transient heat fluxes operating conditions. The test elements use as a heat source a pseudo-chip able to mimic the behavior of a CPU. It is formed by 35 independent sub-heaters, each having its own temperature sensor, such that 35 temperature and 35 heat flux measurements can be made simultaneously. Careful measurements of each pressure drop component (inlet, microchannels and outlet) found in the micro-evaporators showed the importance of the inlet and outlet restriction pressure losses. The overall pressure drop levels found in the copper test section were low enough to possibly be driven by a thermosyphon system. The heat transfer coefficients measured for uniform heat flux conditions were very high and typically followed a V-shape curve. The first branch was associated to the slug flow regime and the second to the annular flow regime. By tracking the minimum level of heat transfer, a transition criteria between the regimes was established, which included the effect of heat flux on the transition. Then for each branch, a different prediction method was used to form the first flow pattern-based prediction method for two-phase heat transfer in microchannels. A non-uniform heat flux creates important temperature gradients in the evaporator, such that the data reduction procedure needs to be adapted to include heat spreading within the evaporator. To do so, a robust multi-dimensional thermal conduction scheme was developed. Once these effects were taken into consideration, the local heat transfer coefficients provided by two-phase flow were found to be the same for uniform and non-uniform heat fluxes, allowing the flow pattern-based method to be extended to non-uniform heat flux conditions. Lastly, with proper control of the mass flow, transient heat flux situations were well handled by the micro-evaporators.

EPFL2011

Numerical Modeling of the Conjugate Heat Transfer Problem for Annular Laminar Film Condensation in Microchannels

Stefano Nebuloni, John Richard Thome

This paper presents numerical simulations of annular laminar film condensation heat transfer in microchannels of different internal shapes. The model, which is based on a finite volume formulation of the Navier-Stokes and energy equations for the liquid phase only, importantly accounts for the effects of axial and peripheral wall conduction and nonuniform heat flux not included in other models so far in the literature. The contributions of the surface tension, axial shear stresses, and gravitational forces are included. This model has so far been validated versus various benchmark cases and versus experimental data available in literature, predicting microchannel heat transfer data with an average error of 20% or better. It is well known that the thinning of the condensate film induced by surface tension due to gravity forces and shape of the surface, also known as the "Gregorig" effect, has a strong consequence on the local heat transfer coefficient in condensation. Thus, the present model accounts for these effects on the heat transfer and pressure drop for a wide variety of geometrical shapes, sizes, wall materials, and working fluid properties. In this paper, the conjugate heat transfer problem arising from the coupling between the thin film fluid dynamics, the heat transfer in the condensing fluid, and the heat conduction in the channel wall has been studied. In particular, the work has focused on three external channel wall boundary conditions: a uniform wall temperature, a nonuniform wall heat flux, and single-phase convective cooling are presented. As the scale of the problem is reduced, i.e., when moving from mini- to microchannels, the results show that the axial conduction effects can become very important in the prediction of the wall temperature profile and the magnitude of the heat transfer coefficient and its distribution along the channel. [DOI: 10.1115/1.4005712]

2012