Study of heat transfer enhancement in a nanofluid-cooled miniature heat sink

This paper reports numerical solution for thermally developing temperature profile and analytical solution for fully developed velocity profile in a miniature plate fin heat sink with SiO₂–water nanofluid as coolant. The flow regime is laminar and Reynolds number varies between 0 and 800. The heat sink is modeled using porous medium approach. Modified Darcy equation for fluid flow and the two-equation model for heat transfer between the solid and fluid phases are employed to predict the local heat transfer coefficient in heat sink. Results show that the nanofluid-cooled heat sink outperforms the water-cooled one, having a considerable higher heat transfer coefficient. The effects of channel aspect ratio and porosity on heat transfer coefficient of the heat sink are studied in detail. Based on the results of our analysis, it is found that an increase in the aspect ratio or the porosity of the plate fin heat sink enhances the heat transfer coefficient.     

Analysis of laminar dissipative flow and heat transfer in a porous saturated circular tube with constant wall heat flux

A theoretical study of the thermal development of forced convection was performed using a circular tube filled with a saturated porous medium, with constant wall heat flux, and with the effect of viscous dissipation. The solution was obtained using the method of separation of variables. The Sturm–Liouville system was solved for the eigenvalues. Ordinary differential equations for the eigenfunctions were solved numerically by the fourth-order Runge–Kutta method. Results show that, in the presence of the viscous dissipation, both the level and distribution of temperature are altered remarkably, even for small values of the Brinkman number, Br, which is the ratio of heat generation caused by viscous dissipation to the value of heat flux at the wall. The value of the local Nusselt number, Nu, is demonstrably independent of Br, unlike the situation in which the wall temperature remains constant.     

The heat transfer and fluid flow of a partially heated microchannel heat sink

Numerical modeling of the conjugate heat transfer in microchannel heat sink is presented. As the most of the cooling applications deals with the partial heated sections, the influence of the heating position on the thermal and hydrodynamic behavior is analyzed. The laminar fluid flow regime and the water as a working fluid are considered. It is observed that partial heating together with variable viscosity has a strong influence on thermal and hydrodynamic characteristics of the micro-heat sink.     

Heat transfer of a parallel flow two-layered microchannel heat sink

A numerical study is conducted to predict the thermal performance of a parallel flow two-layered microchannel heat sink on heat transfer and compared to the case of counterflow for various channel aspect ratios. Findings reveal that the parallel flow configuration leads to a better heat transfer performance except for high Reynolds number and high channel aspect ratio. Further study on the horizontal rib thickness shows that lower thermal resistance can be achieved in a parallel flow two-layered microchannel heat sink with smaller thickness of middle rib.     

Numerical investigation of fluid flow and heat transfer under high heat flux using rectangular micro-channels

A 3D-conjugate numerical investigation was conducted to predict heat transfer characteristics in a rectangular cross-sectional micro-channel employing simultaneously developing single-phase flows. The numerical code was validated by comparison with previous experimental and numerical results for the same micro-channel dimensions and classical correlations based on conventional sized channels. High heat fluxes up to 130 W/cm² were applied to investigate micro-channel thermal characteristics. The entire computational domain was discretized using a 120 × 160 × 100 grid for the micro-channel with an aspect ratio of (α = 4.56) and examined for Reynolds numbers in the laminar range (Re 500–2000) using FLUENT. De-ionized water served as the cooling fluid while the micro-channel substrate used was made of copper. Validation results were found to be in good agreement with previous experimental and numerical data [1] with an average deviation of less than 4.2%. As the applied heat flux increased, an increase in heat transfer coefficient values was observed. Also, the Reynolds number required for transition from single-phase fluid to two-phase was found to increase. A correlation is proposed for the results of average Nusselt numbers for the heat transfer characteristics in micro-channels with simultaneously developing, single-phase flows.     

Bounds on natural convective heat transfer in a porous layer with fixed heat flux

Using the background field variational method, bounds on natural convective heat transfer in a porous layer heated from below with fixed heat flux are derived from the primitive equations. The enhancement of heat transfer beyond the minimal conduction value (the Nusselt number Nu) is bounded in terms of the non-dimensional forcing scale set by the ‘effective’ Rayleigh number (Rˆ) according to Nu ≤ 0.354Rˆ^1/2 and in terms of the conventional Rayleigh number (Ra) defined by the temperature drop across the layer according to Nu ≤ 0.125Ra. It is presented that fixing the heat flux at the boundaries does not change the linear dependence between Nusselt number and Rayleigh number at high Rayleigh number region.     

Numerical investigation on the heat transfer and flow in the mini-fin heat sink for CPU

Fluid flow and heat transfer in the mini-rectangular fin heat sink for CPU of PC using de-ionized water as working fluid are numerically investigated. Based on the real PC operating conditions, the three-dimensional governing equations for fluid flow and heat transfer characteristics are solved using finite volume scheme. The standard k–ε turbulent model is employed to describe the flow structure and behavior. The predicted results obtained from the model are verified by the measured data. There is a reasonable agreement between the predicted results and experiments. The results of this study are expected to lead to guidelines that will allow the design of the cooling system with improved cooling performance of the electronic equipments increasing reliable operation of these devices.     

Forced convection heat transfer of Williamson fluid flow in porous media over horizontal plate with constant heat flux

Forced convection of Williamson fluid flow in porous media under constant surface heat flux conditions is investigated numerically. A model of Darcy–Forchheimer–Brinkman is used and the corresponding governing equations are expressed in dimensionless forms and solved numerically using bvp4c with MATLAB package. Boundary layer velocity, shear stress, and temperature profiles, in addition to the local Nusselt number parameter over a horizontal plate, are found. The effects of the Forchheimer parameter, Nusselt number, Darcy parameter, porous inertia, and Williamson parameter on the velocity profiles, temperature profiles, coefficient of friction, and coefficient of heat transfer are investigated. The results showed that as the Darcy parameter increases, boundary layer velocity and shear stress increase, while the temperature and Nusselt number decrease. In addition, as Williamson's parameter increases, velocity within the boundary layer, shear stress, and Nusselt number decrease while the temperature profile increases. Also, with larger values of the Forchheimer parameter, the velocity of the boundary layer, shear stress, temperature, and Nusselt number increase. Furthermore, the Nusselt number and the coefficient of friction are obtained on the surface of the horizontal plate.     

Conjugate heat transfer in fractal-shaped microchannel network heat sink for integrated microelectronic cooling application

The hydrodynamic and thermal characteristics of fractal-shaped microchannel network heat sinks are investigated numerically by solving three-dimensional N–S equations and energy equation, taking into consideration the conjugate heat transfer in microchannel walls. It is found that due to the structural limitation of right-angled fractal-shaped microchannel network, hotspots may appear on the bottom wall of the heat sink where the microchannels are sparsely distributed. With slight modifications in the fractal-shaped structure of microchannels network, great improvements on hydrodynamic and thermal performance of heat sink can be achieved. A comparison of the performance of modified fractal-shaped microchannel network heat sink with parallel microchannels heat sink is also conducted numerically based on the same heat sink dimensions. It is found that the modified fractal-shaped microchannel network is much better in terms of thermal resistance and temperature uniformity under the conditions of the same pressure drop or pumping power. Therefore, the modified fractal-shaped microchannel network heat sink appears promising to be used for microelectronic cooling in the future.     

"水滴"型微针肋流动与传热特性数值研究

基于圆肋尾部流动传热强化的思想,对传统圆肋进行优化,进而设计出了新型的"水滴"型微针肋.采用Fluent模拟软件对不同角度"水滴"型针肋侧壁及整体流动传热特性进行了模拟分析.结果表明:a>90°时,流动没有显著改变,边界层涡依旧发生脱落;而当a<60°后,针肋尾部逆压梯度较小,边界层涡脱离现象消失;适当的"水滴"型结构...     

Local heat transfer distribution and effect of instabilities during flow boiling in a silicon microchannel heat sink

Flow boiling of the perfluorinated dielectric fluid FC-77 in a silicon microchannel heat sink is investigated. The heat sink contains 60 parallel microchannels each of 100 μm width and 389 μm depth. Twenty-five evenly distributed temperature sensors in the substrate yield local heat transfer coefficients. The pressure drop across the channels is also measured. Experiments are conducted at five flow rates through the heat sink in the range of 20–80 ml/min with the inlet subcooling held at 26 K in all the tests. At each flow rate, the uniform heat input to the substrate is increased in steps so that the fluid experiences flow regimes from single-phase liquid flow to the occurrence of critical heat flux (CHF). In the upstream region of the channels, the flow develops from single-phase liquid flow at low heat fluxes to pulsating two-phase flow at high heat fluxes during flow instability that commences at a threshold heat flux in the range of 30.5–62.3 W/cm² depending on the flow rate. In the downstream region, progressive flow patterns from bubbly flow, slug flow, elongated bubbles or annular flow, alternating wispy-annular and churn flow, and wall dryout at highest heat fluxes are observed. As a result, the heat transfer coefficients in the downstream region experience substantial variations over the entire heat flux range, based on which five distinct boiling regimes are identified. In contrast, the heat transfer coefficient midway along the channels remains relatively constant over the heat flux range tested. Due to changes in flow patterns during flow instability, the heat transfer is enhanced both in the downstream region (prior to extended wall dryout) and in the upstream region. A previous study by the authors found no effect of instabilities during flow boiling in a heat sink with larger microchannels (each 300 μm wide and 389 μm deep); it appears therefore that the effect of instabilities on heat transfer is amplified in smaller-sized channels. While CHF increases with increasing flow rate, the pressure drop across the channels has only a minimal dependence on flow rate once boiling is initiated in the microchannels, and varies almost linearly with increasing heat flux.     

Experimental investigation of heat transfer coefficient of R134a during condensation in vertical downward flow at high mass flux in a smooth tube

The two-phase heat transfer coefficients of pure HFC-134a condensing inside a smooth tube-in-tube heat exchanger are experimentally investigated. The test section is a 0.5 m long double tube with refrigerant flowing in the inner tube and cooling water flowing in the annulus. The inner tube is constructed from smooth copper tubing of 9.52 mm outer diameter and 8.1 mm inner diameter. The test runs are performed at average saturation condensing temperatures between 40–50 °C. The mass fluxes are between 260 and 515 kg m^− 2s^− 1 and the heat fluxes are between 11.3 and 55.3 kW m^− 2. The quality of the refrigerant in the test section is calculated using the temperature and pressure obtained from the experiment. The average heat transfer coefficient of the refrigerant is determined by applying an energy balance based on the energy transferred from the test section. The effects of heat flux, mass flux and condensation temperature on the heat transfer coefficients are also discussed. Eleven well-known correlations for annular flow are compared to each other using a large amount of data obtained from various experimental conditions. A new correlation for the condensation heat transfer coefficient is proposed for practical applications.     

Conjugate numerical analysis of flow and heat transfer with phase change in a miniature flat plate CPL evaporator

Two-dimensional numerical model for the global evaporator of miniature flat plate capillary pumped loop (CPL) is developed to describe heat and mass transfer with phase change in the porous wick, liquid flow and heat transfer in the compensation cavity and heat transfer in the vapor grooves and metallic wall. The governing equations for different zones are solved as a conjugate problem. The side wall effect heat transfer limit is introduced to estimate the heat transport capability of evaporator. The influences of liquid subcooling, wick material, metallic wall material and non-uniform heat flux on the evaporator performance are discussed in detail.     

Melting heat transfer analysis of electrically conducting nanofluid flow over an exponentially shrinking/stretching porous sheet with radiative heat flux under a magnetic field

Modern magnetic nanomaterial processing operations are progressing rapidly and require increasingly sophisticated mathematical models for their optimization. Stimulated by such developments, in this paper, a theoretical and computational study of a steady magnetohydrodynamic nanofluid over an exponentially stretching/shrinking permeable sheet with melting (phase change) and radiative heat transfer is presented. Besides, wall transpiration, that is, suction and blowing (injection), is included. This study deploys Buongiorno's nanofluid model, which simulates the effects of the Brownian motion and thermophoresis. The transport equations and boundary conditions are normalized via similarity transformations and appropriate variables, and the similarity solutions are shown to depend on the transpiration parameter. The emerging dimensionless nonlinear coupled ordinary differential boundary value problem is solved numerically with the Newton-Fehlberg iteration technique. Validation with special cases from the literature is included. The increase in the magnetic field, that is, the Hartmann number, is observed to elevate nanoparticle concentration and temperature, whereas it dampens the velocity. Higher values of the melting parameter consistently decelerate the boundary layer flow and suppress temperature and nanoparticle concentration. A higher radiative parameter strongly increases temperature (and thermal boundary layer thickness) and weakly accelerates the flow. The increase in the Brownian motion reduces nanoparticle concentrations, whereas a greater thermophoretic body force strongly enhances them. The Nusselt number and Sherwood number are observed to be decreased with an increasing Hartmann number, whereas they are elevated with a stronger wall suction and melting parameter.     

Analysis of heat transfer for compressible flow in two-dimensional porous media

Gas from a reservoir at constant pressure and temperature is forced through a two-dimensional porous region. The surface through which the gas exits is at a specified uniform temperature and pressure. The local gas and solid matrix temperatures are assumed equal. General solutions for the local temperature and pressure in the porous medium are found as a function of a potential. This potential can be determined by solving Laplace's equation in the porous region for a simple set of boundary conditions, and the temperature and pressure will then be known functions of position. Because of the nature of the boundary conditions it is particularly convenient to solve Laplace's equation by conformai mapping. By using this technique some illustrative heat and mass flow results were calculated for a porous wall with a step in thickness, a wall supplied with gas through periodic slots, and an eccentric annular region.     

Micropolar flow in a porous channel with high mass transfer

In this research the injective micropolar flow in a porous channel is investigated. The flow is driven by suction or injection on the channel walls, and the micropolar model is used to describe the working fluid. This problem is mapped into the system of nonlinear coupled differential equations by using Berman's similarity transformation. These are solved for large mass transfer via Optimal Homotopy Asymptotic Method (OHAM). Also the numerical method is used for the validity of this analytical method and excellent agreement is observed between the solutions obtained from OHAM and numerical results. Trusting this validity, effects of some other parameters are discussed.     

Hydromagnetic flow and heat transfer past a continuously moving porous boundary

The effect of magnetic field on the flow and heat transfer past a continuously moving porous plate in a stationary fluid has been analysed. The governing boundary layer equations have been reduced to a set of nonlinear ordinary differential equations using similarity transformations. The resulting boundary value problem has been solved numerically. The effects of magnetic and suction (or injection) parameters on the velocity and temperature profiles as well as on the skin friction and heat transfer coefficients have been studied. It has been observed that the effect of magnetic field is to increase the wall skin friction while the reverse occurs in the case of Nusselt number.     

Buoyancy-induced flow,heat, and mass transfer in a porous annulus

This paper reports on double-diffusive natural convection heat transfer in a porous annulus between concentric horizontal circular and square cylinders. A pressure-based segregated finite volume method is used to solve the problem numerically. The diffusion fluxes are discretized using the MIND fully implicit scheme. Furthermore, a modified pressure correction equation is derived that implicitly accounts for the nonorthogonal diffusion terms, which are usually neglected in the standard SIMPLE algorithm. Results indicate that convection effects increase with an increase in Rayleigh number, Darcy number, porosity, and enclosure aspect ratio. Further, at low Darcy values, porosity has no effect on the flow, temperature, and concentration fields.