Computational heat transfer and two-phase flow topology in miniature tubes

Detailed computational multi-fluid dynamics simulations have been performed to study the effect of two-phase flow regime on heat transfer in small diameter pipes. Overall the heat removal rate in two-phase flow is higher than in single phase. Subtle differences in thermal removal rates are revealed when the flow-regime transitions from bubbly to slug and slug-train configurations. It is found that the wall thermal layer is affected by two separate mechanisms: an early-stage compression due to gas-jet fragmentation into slugs or bubbles, and a background inclusion-induced flow superimposed on the equivalent single-phase fully developed flow far downstream. The first mechanism resembles a confinement or blockage effect, and is shown to directly influence radial temperature gradients. The downstream mechanism is a cell-based developed flow (rather than fully developed), and is shown here to increase the wall shear in the vicinity of the cell, leading to higher heat transfer rates. The mean Nusselt number distribution shows that the bubbly, slug and slug-train patterns transport as much as three to four times more heat from the tube wall to the bulk flow than pure water flow. A mechanistic heat transfer model is proposed, based on frequency and length scale of inclusions.     

Finite-element analysis of a rotary oven for bread

The paper deals with a finite-element analysis of the turbulent forced-convection air flow inside a rotary oven for bread. The code ADINA-F is used for the numerical finite-element simulation. To describe the flow and the heat transfer inside the oven, a two-dimensional model has been adopted. The numerical heat-transfer convection coefficient is in good agreement with the experimental data thus providing a validation of the approximate model formulated in the numerical analysis. A structural change intended to enhance the convective heat transfer has been suggested and numerically tested. The results show a mean Nusselt number enhancement of about 30%. The finite-element analysis is, therefore, a helpful tool in estimating the performances of this kind of oven.     

Analytical solution of mixed electromagnetic/pressure driven gaseous flows in microchannels

The influences of wall-slip/jump conditions on the fluid flow and heat transfer for hydrodynamically and thermally fully developed electrically conducting gaseous flow subject to an electromagnetic field inside a parallel plate microchannel with constant heat flux at walls are studied under the assumptions of a low-magnetic Reynolds number. The governing equations are non-dimensionalized and then analytical solutions are derived for the friction and the heat transfer coefficients. The fluid flow and the heat transfer characteristics obtained in the analytical solutions are discussed in detail for different parameters such as the Knudsen, Hartmann, and Brinkman numbers. The velocity profiles verify that even with a constant Knudsen number, applying a stronger electromagnetic field gives rise to an increase in the slip velocity. The results also reveal that on increasing the Hartmann number, the heat transfer rate as well as the friction factor is enhanced, whereas it tends to suppress the movement of the fluid. Further, it is found that the Nusselt and the Poiseuille numbers are less sensitive to the electromagnetic field effects with increase in rarefaction.     

Numerical study of seed bubble-triggered evaporation heat transfer in a single microtube

High boiling incipience temperature and flow instabilities in silicon-based microchannels with smooth surface are challenging issues. This work numerically investigates the seed bubble-triggered evaporation heat transfer in a microtube, with a length of 5.0 mm and diameter of 106 μm. Acetone was the working fluid. Seed bubbles were assumed to be generated periodically at the microtube upstream. The fixed grid allocation technique was proposed to successfully perform the parallel computation via a set of computer core solvers. It is found that the seed bubble-guided heat transfer consists of a start-up stage and a steady operation stage. The start-up time equals to the residence time of the first seed bubble growing and traveling in the microtube. The seed bubble frequency is a key parameter to influence the performance. Low-frequency seed bubbles cause alternative flow patterns of liquid flow and elongated bubble flow, corresponding to the apparent spatial-time oscillations of wall and bulk fluid superheats. High-frequency seed bubbles result in quasi-stable elongated bubble flow, corresponding to quasi-uniform and stable wall and fluid superheats. There is a saturation seed bubble frequency beyond which no further performance improvement can be made. There are residual fluid superheats specifying the required minimum superheats to sustain the evaporation heat transfer between the two phases. Elongated bubbles with thin liquid films are responsible for the heat transfer enhancement. Contrary to wall temperatures, the transient local Nusselt numbers are slightly changed due to the fact that heat transfer is more closely related to the dynamic elongated bubble flow evolution within millisecond timescale in the microchannel. The heat transfer coefficients can be 2.0 to 3.5 times of that for the superheated liquid flow before seed bubble injections.     

Airside fin efficiencies for finned-tube heat exchangers with forced convection

The heat transfer and mass transfer fin efficiencies were analyzed numerically to show that popular models for heat transfer fin efficiency for circular fins are not always reasonable.The numerical results show that the effective heat transfer area of a circular fin increases several times faster than that of a straight fin for the same tube radius.Then,a simple but accurate heat transfer fin efficiency model was developed and verified by numerical results for a wide range of fin designs.This model predicts...     

Interlayer cooling potential in vertically integrated packages

The heat-removal capability of area-interconnect-compatible interlayer cooling in vertically integrated, high-performance chip stacks was characterized with de-ionized water as coolant. Correlation-based predictions and computational fluid dynamic modeling of cross-flow heat-removal structures show that the coolant temperature increase due to sensible heat absorption limits the cooling performance at hydraulic diameters ≤200 μm. An experimental investigation with uniform and double-side heat flux at Reynolds numbers ≤1,000 and heat transfer areas of 1 cm² was carried out to identify the most efficient interlayer heat-removal structure. The following structures were tested: parallel plate, microchannel, pin fin, and their combinations with pins using in-line and staggered configurations with round and drop-like shapes at pitches ranging from 50 to 200 μm and fluid structure heights of 100–200 μm. A hydrodynamic flow regime transition responsible for a local junction temperature minimum was observed for pin fin in-line structures. The experimental data was extrapolated to predict maximal heat flux in chip stacks having a 4-cm² heat transfer area. The performance of interlayer cooling strongly depends on this parameter, and drops from >200 W/cm² at 1 cm² and >50 μm interconnect pitch to <100 W/cm² at 4 cm². From experimental data, friction factor and Nusselt number correlations were derived for pin fin in-line and staggered structures.     

Gaseous slip flow forced convection through ordered microcylinders

This is a theoretical study dealing with longitudinal gaseous slip flow forced convection between a periodic bunch of microcylinders arranged in regular array. The selected geometry has applications in microscale pin fin heat sinks used for cooling of microchips. The flow is considered to be hydrodynamically and thermally fully developed. The two axially constant heat flux boundary conditions of H1 and H2 are considered in the analysis. The velocity and temperature discontinuities at the boundary are incorporated into the solutions using the first order slip boundary conditions. The method considered is mainly analytical in which the governing equations and three of the boundary conditions are exactly satisfied. The remaining symmetry condition on the right-hand boundary of the typical element is applied to the solution through the point matching technique. The results show that both the Poiseuille number and the Nusselt number are decreasing functions of the degree of rarefaction characterized by the Knudsen number. While an increase in the blockage ratio leads to a higher Poiseuille number, the functionality of the Nusselt number on this parameter is not monotonic. At small and moderate values of the blockage ratio, the Nusselt number is higher for a higher blockage ratio, whereas the opposite may be right for higher values of this parameter. It is also observed that the angular variations of the parameters are reduced at smaller blockage ratios. Accordingly, the H1 and H2 Nusselt numbers are the same for small and moderate blockage ratios.     

Forced convection in electro-osmotic/Poiseuille micro-channel flows of viscoelastic fluids: fully developed flow with imposed wall heat flux

The analytical solution for heat transfer in a dynamic and thermally fully developed channel flow of the simplified Phan-Thien–Tanner fluid induced by combined electro-osmosis and pressure gradient was obtained assuming that material properties are independent of temperature. The flow forcing was quantified by an appropriate dimensionless parameter and its effect and that of all other relevant dimensionless numbers is presented and discussed. Specifically, the forced convection occurs under conditions of constant wall heat flux and the solution includes the effects of Weissenberg number, electric double layer (EDL) thickness, forcing ratio parameter, viscous dissipation as well as of Joule heating due to the electric currents and was obtained under the simplifying Debye–Hückel approximation. Generally speaking, the Joule effect is stronger than the viscous dissipation except in very narrow channels, but these fall outside the validity of the Debye–Hückel conditions. For pure electro-osmosis, viscous dissipation is restricted to the near-wall region and virtually nonexistent elsewhere, so it is irrelevant for thin electric double layers and Joule heating is more relevant. As the EDL thickens and/or the pressure gradient contribution increases, the role of viscous dissipation grows and shear-thinning effects also appear more clearly on the Nusselt number. Generally speaking, an increase in internal heating results in lower Nusselt numbers and this effect is stronger than the effect of shear-thinning, which is responsible for a slight increase in the Nusselt number.     

降膜蒸发过程的数值模拟和传热传质分析

利用CFD软件对逆流降膜蒸发过程进行了实验模拟研究,研究了速度边界层、热边界层和浓度边界层的变化对降膜蒸发传热传质特性的影响规律;通过建立一维逆流降膜蒸发的数学方程编程求解出了对流传热传质的Nu数和Sh数,利用Fluent软件模拟出的实验结果采用回归分析得出了气液流量比Raw、流道的长宽比αL、空气进口无量纲温度θai以及空气进口Re数与Nu数、Sh数之间的无量纲关系式,可为降膜蒸发换热器的设计提供参考。     

Laminar fluid flow and heat transfer in an internally corrugated tube by means of the method of fundamental solutions and radial basis functions

The paper shows application of the method of fundamental solutions in combination with the radial basis functions for analysis of fluid flow and heat transfer in an internally corrugated tube. Cross-section of such a tube is mathematically described by a cosine function and it can potentially represent a natural duct with internal corrugations, e.g. inside arteries. The boundary value problem is described by two partial differential equations (one for fluid flow problem and one for heat transfer problem) and appropriate boundary conditions. During solving this boundary value problem the average fluid velocity and average fluid temperature are calculated numerically. In the paper the Nusselt number and the product of friction factor and Reynolds number are presented for some selected geometrical parameters (the number and amplitude of corrugations). It is shown that for a given number of corrugations a minimal value of the product of friction factor and Reynolds number can be found. As it was expected the Nusselt number increases with increasing amplitude and number of corrugations.     

Numerical investigation of heat transfer enhancement by carbon nano fibers deposited on a flat plate

Numerical simulations of flow and heat transfer have been performed for flow over a plate surface covered with carbon nano fibers (CNFs). The CNFs influence on fluid flow and heat transfer has been investigated. Firstly, a stochastic model for CNFs deposition has been explained. Secondly, the lattice Boltzmann model for simulating fluid flow and heat transfer is described. Finally, results from calculations for the heat transfer are presented. The results show substantial heat transfer enhancement for a densely covered surface with CNFs of varying length. From a detailed analysis of the results it is concluded that the enhancement is caused by boundary layer regeneration.     

不同错开位置锯齿翅片热力特性的三维仿真

对通道壁面沿流动方向周期性地错开而在横向上通道壁面有几组不同错开位置的强化换热翅片的层流流动和传热情况进行了三维仿真分析。这类翅片经常被应用在汽车和电子设备的冷却系统中。三维有限元方法用来模拟Re数在100-1500范围时翅片的速度场和温度场，考虑了通道壁面的不同错开位置(SSP)对翅片传热性能的影响。仿真结果与以前的试验数据也进行了比较研究。     

MPTL板翅式换热器性能的数值分析

研究机械泵驱动的两相冷却回路(MPTL)换热器传热效能优化问题,由于换热器的应用环境极为复杂,目前的实验数据极度匮乏,而且结构特殊,现有的标准叉流板式换热器的设计方法并不适用。为了提出可行的设计优化方法,针对某典型MPTL板翅式换热器,建立一种简化的三维物理数学模型,采用有限体积法及流固耦合方法深入研究该换热器在典型结构及设计工况下的流场特性,随后研究其在不同翅高及进口雷诺数时的换热性能变化规律。数值仿真结果表明,换热器换热效率较高,压力损失较小,且随着翅高的降低、进口雷诺数的增加,换热性能逐渐提高。结果证明,可为换热器的进一步设计与优化提供有益的参考。     

Numerical simulation of heat transfer enhancement by elastic turbulence in a curvy channel

Although many investigations on elastic turbulence have been conducted in recent years, two major research topics still call for in-depth mechanistic investigations. Specifically, one is heat transfer performance affected by elastic turbulence; the other is so-called high Weissenberg number problem (HWNP) in numerical simulation of viscoelastic fluid flow. Taking these two topics into account simultaneously, the coupled problem becomes heat transfer characteristic of viscoelastic fluid in elastic turbulence at high Weissenberg number (Wi) and very low Reynolds number (Re). In this work, we implement numerical simulations by embedding log-conformation reformulation algorithm into the open-source software OpenFOAM. The heat transfer process of viscoelastic fluid flow in a three-dimensional (3D) curvy channel is simulated over a wide range of Wi. For the first time, significant heat transfer enhancement induced by elastic turbulence in a curvy channel at high Wi was identified numerically. When Wi is above the critical value of O(1), the heat transfer performance is found to be dramatically improved by elastic turbulence and then approaches a saturation. From the transient analysis of flow motions in the axial and cross sections, it can be seen that the flow twists and wiggles in the curvy channel and the field synergy effect of viscoelastic fluid flow becomes more intensive than that of Newtonian fluid flow. These effects give rise to the extremely irregular flow motions in the cross section and consequently lead to heat transfer enhancement.     

Direct computation of perturbed impinging hot jets

The unsteady flow and temperature fields of an impinging hot jet at a Reynolds number of 1000 and a nozzle-to-plate distance of 6 jet diameters have been obtained by direct numerical solution of the compressible time-dependent three-dimensional Navier-Stokes equations using highly accurate numerical methods. Effects of an external perturbation on the flow and heat transfer characteristics of the transitional impinging jet have been examined. Oscillatory behaviour induced by the external perturbation has been observed for the impinging jet. The external perturbation leads to the large-scale vortical structures in the primary jet stream, which subsequently lead to the strong oscillatory behaviour of the impinging jet. The vortical structures lead to flow transitional behaviour that enhances mixing of the hot jet with the ambient fluid. It has been observed that the wall boundary layer of the impinging jet is thin. Both the instantaneous and time-averaged wall shear and normal stresses and Nusselt number are examined. Although the external perturbation strongly affects the flow structures in the primary jet stream, it does not have significant effects on the wall stresses and heat transfer characteristics of the impinging jet due to the re-laminarization effect of the wall.     

Effect of particle migration on heat transfer in suspensions of nanoparticles flowing through minichannels

Recent work has shown that suspensions of highly thermally conducting nanoparticles with a size considerably smaller than 100 nm have great potential as a high-energy carrier for small channel systems. However, it is also known that particles in a suspension under certain conditions may migrate. This indicates that the efficiency of heat transfer in the small channels may not be as superior as expected, which bears significance to the system design and operation. This work aims at addressing this issue by examining the effect of particle migration on heat transfer under a fully developed laminar flow regime in small channels. This involves the development of both flow and heat transfer models, and a numerical solution to the models. The flow model takes into account the effects of the shear-induced and viscosity-gradient-induced particle migration, as well as self-diffusion due to Brownian motion, which is coupled with an energy equation. The results suggest a significant non-uniformity in particle concentration and, hence, thermal conductivity over the tube cross-section due to particle migration, particularly for large particles at high concentrations. Compared with the constant thermal conductivity assumption, the non-uniform distribution due to particle migration leads to a higher Nusselt number, which depends on the Peclet number and the mean particle concentration. Further improvement of the model is needed to take into account other factors such as entrance effects, as well as the dynamics of particles and particle–wall interactions.     

Enhancement of Heat Transfer in a Tube with Swirl Generators

Influence of strip self-rotating plastic spiral elements inserted in a tube on heat transfer enhancement is studied experimentally. The strip consists of 100 self-rotating spiral elements made of plastic polycarbonate inserted in the inner tube of a concentric tube heat exchanger with a view to generating swirl flow that helps to increase the heat transfer rate of the tube. Cold water flows in the annulus whereas the hot water flows in the inner tube. The obtained experimental data are compared with the data obtained from plain tube .Experimental results confirmed that the use of self-rotating plastic spiral elements leads to higher overall heat transfer coefficient than the plain tube. This technology is not only useful for heat transfer enhancement but also can be used for self cleaning deposition fouling in heat exchanger when the flow velocity is higher than 0.2 m/sec.     

汽车空调平流式冷凝器性能仿真分析

为分析汽车空调平流式冷凝器的换热、流动性能,假设制冷剂沿管长方向做一维流动,空气侧流动视为零维流动,忽略制冷剂加速压降,对制冷剂两相区采用均相模型.使用AMESim建立平流式冷凝器仿真模型,并通过与试验对比验证模型的准确性.改变冷凝器结构参数,分析对冷凝器的性能影响,发现合理的制冷剂回路流程布置可以改善平流式冷凝器性能;增加流程数可以增加换热量,但是压降也会增大;制冷剂侧总横截面积相等时,微通道数目增加,换热量增加;空气速度较小时,减小翅片间距可以增加换热量.     

Efficient heat transfer enhancement by elastic turbulence with polymer solution in a curved microchannel

Heat and mass transfer in microscale flows are limited due to extremely low Reynolds number (Re). In a curved microchannel, however, complex flow behaviors, such as elastic instability and elastic turbulence, can be induced via viscoelastic fluid at vanishingly low-Re conditions, which is of great potential to enhance the heat transfer performance. The influence of elastic instabilities and turbulence on heat dissipation of exothermic components is experimentally investigated in this study. The heat transfer performance of both viscoelastic (polymer solutions) and Newtonian (sucrose solutions) fluid flows in a curved microchannel with a square cross section is experimentally characterized. Titanium–platinum (Ti–Pt) thin films embedded at the bottom wall of the polydimethylsiloxane (PDMS) microchannel serve as both microheater and temperature sensor. For viscoelastic fluids, the spectrum of outlet temperature fluctuation in broad frequency (f) region fits the power law of f ^?1.1. Heat transfer enhancement due to the elastic turbulence in a curved microchannel is thereby identified by the drastic growth of the Nusselt number (Nu, the ratio of convective to conductive heat transfer normal to the boundary) with the increase in the Weissenberg number (Wi, the ratio of elastic stress to viscous stress). The mechanism of heat transfer enhanced by the convection effect of elastic turbulence is also elucidated.