Fluid flow and heat transfer in a two-dimensional wavy channel

A numerical study is presented for the laminar fully developed flow and heat transfer in a two-dimensional wavy channel. The effects of the geometry, Reynolds and Prandtl number on the flow field and heat transfer are investigated. The channel is characterized by a wavy wall, heated at uniform heat flux, and an opposite wall, being plane and adiabatic. The extent of the wall waviness and the distance between the channel walls are found to significantly affect the streamlines contours as well as the heat transfer coefficients. Comparisons with the straight channel, in the same flow rate and heat transfer conditions, have been performed. Pressure drop of the wavy channel is found to be always larger than the value characteristic of a straight channel, while heat transfer performance decreases or increases depending on the values of the parameters (geometry, Reynolds and Prandtl numbers).     

Visualization and determination of local heat transfer coefficients in shell-and-tube heat exchangers for in-line tube arrangement by mass transfer measurements

The local heat transfer coefficients on the shell-side of shell-and-tube heat exchangers for in-line tube arrangement are visualized and determined from mass transfer measurements. The mass transfer experiments are carried out using a technique based on absorption, chemical and coupled colour reaction. Local mass transfer coefficients are measured for fully developed flow conditions on each tube surface. These coefficients were transformed to heat transfer coefficients by employing the analogy between heat and mass transfer. The averaged heat transfer coefficients and the pressure drop are compared with the predictions from the literature. Received on 2 May 1997     

Visualisation of flow boiling heat transfer in a microtube

This paper presents the results of the flow boiling patterns and heat transfer coefficients of FC-72 in a small tube. The internal diameter of the tube is 0.48 mm, with a heated length of 73 mm. The mass flow rate varies from 50 to 3,000 kg/m²-s. The microtube is made of Pyrex in order to obtain the visualisation of the flow pattern along the heated channel. Different types of flow pattern have been observed: bubbly flow, deformed bubbly flow, bubbly/slug flow, slug flow, slug/annular flow, and annular flow. The experiments show the presence of flow instabilities in a large portion of the tests at low mass flow rates and low subcooling. Flow patterns in presence of flow instabilities are mainly characterized by bubbly/slug flow and slug/annular flow. Heat transfer rates have been studied in all flow pattern conditions. The two groups of results, with flow instabilities and without flow instabilities, show similar heat transfer behaviour. The heat transfer characteristics of the pipes have been studied in comparison with mass flux and vapour quality.     

Investigation of temperature statistic characteristics in turbulent water flow with unsteady heat transfer

The results of an experimental study of a temperature field and its statistical characteristics in turbulent water flow upon a sudden change of heat generation in the channel wall are reported. Measurements were performed in 5 mm × 40 mm, 10 mm × 40 mm, and 20 mm × 40 mm channels in the regions of thermal stabilization and stabilized heat transfer at Reynolds numbers of (0.8–6.8) × 10⁴. The measurement results are generalized using a dimensionless time scale. The results of the calculation of heat transfer coefficients at unsteady heat transfer are presented.     

Macro-to-microchannel transition in two-phase flow: Part 2 - Flow boiling heat transfer and critical heat flux

This part of the paper presents the current experimental flow boiling heat transfer and CHF data acquired for R134a, R236fa and R245fa in single, horizontal channels of 1.03, 2.20 and 3.04 mm diameters over a range of experimental conditions. The aim of this study is to investigate the effects of channel confinement, heat flux, flow pattern, saturation temperature, subcooling and working fluid properties on the two-phase heat transfer and CHF. Experimentally, it was observed that the flow boiling heat transfer coefficients are a significant function of the type of two-phase flow pattern. Furthermore, the monotonically increasing heat transfer coefficients at higher vapor qualities, corresponding to annular flow, signifies convective boiling as the dominant heat transfer mechanism in these small scale channels. The decreasing heat transfer trend at low vapor qualities in the slug flow (coalescing bubble dominated regime) was indicative of thin film evaporation with intermittent dry patch formation and rewetting at these conditions. The coalescing bubble flow heat transfer data were well predicted by the three-zone model when setting the dryout thickness to the measured surface roughness, indicating for the first time a roughness effect on the flow boiling heat transfer coefficient in this regime. The CHF data acquired during the experimental campaign indicated the influence of saturation temperature, mass velocity, channel confinement and fluid properties on CHF but no influence of inlet subcooling for the conditions tested. When globally comparing the CHF values for R134a in the 0.51-3.04 mm diameter channels, a peak in CHF peak was observed lying in between the 0.79 (Co ≈ 0.99) and 1.03 (Co ≈ 0.78) mm channels. A new CHF correlation has been proposed involving the confinement number, Co that is able to predict CHF for R134a, R236fa and R245fa in single-circular channels, rectangular multichannels and split flow rectangular multichannels. In summary, the present flow boiling and CHF trends point to a macro-to-microscale transition as indicated by the results presented in Ong and Thome (2011) [1].     

The developing heat transfer and fluid flow in micro-channel heat sink with viscous heating effect

The numerical modeling of the conjugate heat transfer and fluid flow through the micro-heat sink was presented in the paper, considering the viscous dissipation effect. Three different fluids with temperature dependent fluid viscosity are considered: water, dielectric fluid HFE-7600 and isopropanol. The square shape of the cross-section is considered with D _h  = 50 μm with a channel length L = 50 mm. As most of the reported researches dealt with fully developed fluid flow and constant fluid properties in this paper the thermal and hydro-dynamic developing laminar fluid flow is analyzed. Two different heat transfer conditions are considered: heating and cooling at various Br. The influence of the viscous heating on local Nu and Po is analyzed. It was shown that for a given geometry the local Po and Nu numbers are strongly affected by the viscous heating. Moreover the Po number attains the fully developed value as the external heating is equal with the internal viscous heating.     

Two-dimensional wall temperature measurements and heat transfer enhancement for top-heated horizontal channels with flow boiling

Two-dimensional (circumferential and axial) wall temperature distributions were measured for top-heated coolant channels with internal geometries that include smooth walls, spiral fins and both twisted tape and spiral fins. Freon-71 was the working fluid. The flow regimes studied were single-phase, subcooled flow boiling, and stratified flow boiling. The inside diameter of all test sections was near 10.0 mm. Circumferentially averaged heat transfer coefficients at several axial locations were obtained for selected coolant channels for a volumetric flow rate of 4.738 x 10⁻⁵m³/s, 0.19 MPa (absolute) exit pressure, and 22.2°C inlet subcooling. Overall (averaged over the entire channel) heat transfer coefficients were compared for the various channel geometries. This comparison showed that the channel with large-pitch spiral fins had higher heat transfer coefficients at all power levels. However, the results appear to indicate that if the twist ratio (ratio of the twisted tape period to the inside diameter) is decreased, the configuration employing both fins and a twisted tape will have had greater enhancements.     

Convective heat transfer of aqueous alumina nanosuspensions in a horizontal mini-channel

This work reports an experimental study of convective heat transfer of aqueous alumina nanofluids in a horizontal mini-channel under laminar flow condition 40 < Re < 1,000. The variation of local heat transfer coefficients, in both entrance and developed flow regimes, was obtained as a function of axial distance. The heat transfer coefficient of nanofluids was found to be dependent on not only nanoparticle concentration but also mass flow rate. Different to the behavior in conventional-sized channels, the major heat transfer coefficient enhancement is shown in the fully developed regime in the minichannel where up to 40% increase is observed. Discussions of the results suggest that apart from the need of a careful assessment of different thermo-physical properties of nanofluids, i.e., viscosity, specific heat and thermal conductivity, the heterogeneous nature of nanoparticle flow should be considered especially under high flow rate conditions.     

An updated three-zone heat transfer model for slug flow boiling in microchannels

This work proposes a novel physics-based model for the fluid mechanics and heat transfer associated with slug flow boiling in horizontal circular microchannels to update the widely used three-zone model of Thome et al. (2004). The heat transfer model has a convective boiling nature and predicts the time-dependent variation of the local heat transfer coefficient during the cyclic passage of a liquid slug, an evaporating elongated bubble and a vapor plug. The capillary flow theory, extended to incorporate evaporation effects, is applied to estimate the bubble velocity along the channel. A liquid film thickness prediction method also considering bubble proximity effects, which may limit the radial extension of the film, is included. The minimum liquid film thickness at dryout is set to the channel wall roughness. Theoretical heat transfer models accounting for the thermal inertia of the liquid film and for the recirculating flow within the liquid slug are utilized. The heat transfer model is compared to experimental data taken from three independent studies. The 833 slug flow boiling data points cover the fluids R134a, R245fa and R236fa, and channel diameters below 1 mm. The proposed evaporation model predicts more than 80% of the database to within ±30%. It demonstrates a stronger contribution to heat transfer by the liquid slugs and correspondingly less by the thin film evaporation process compared to the original three-zone model. This model represents a new step towards a complete physics-based modelling of the bubble dynamics and heat transfer within microchannels under evaporating flow conditions.     

Boiling heat transfer,pressure drop,and flow pattern in a horizontal square minichannel

This study experimentally investigated the flow boiling heat transfer, pressure drop, and flow pattern in a horizontal square minichannel with a hydraulic diameter of 2.0 mm, and the effects of mass flux, vapor quality, heat flux, and refrigerant properties on the flow boiling characteristics were clarified. The heat transfer coefficient and pressure drop of R32 and R1234yf were measured in a mass flux range of 50–400 kgm⁻²s⁻¹ at a saturation temperature of 15 °C. The flow pattern of the square minichannel outlet was observed and was classified as plug, wavy, churn, and annular flows. The heat transfer coefficients in the square minichannel were larger than those in the circular minichannel with a similar hydraulic diameter at low mass flux conditions. The heat transfer coefficients of R32 indicated higher values compared with those of R1234yf at same mass flux and qualities. An empirical heat transfer model taking into account the forced convection, nucleate boiling, and thin liquid film evaporation was developed for horizontal square and circular minichannels. The frictional pressure drop of R32 was 1.5–2 times higher than that of R1234yf at same mass flux and vapor quality condition, and the effect of channel shape on the frictional pressure drop was small unlike the boiling heat transfer.     

Numerical investigation of heat transfer and fluid flow characteristics inside a wavy channel

The periodically fully developed laminar heat transfer and fluid flow characteristics inside a two-dimensional wavy channel in a compact heat exchanger have been numerically investigated. Calculations were performed for Prandtl number 0.7, and Reynolds number ranging from 100 to 1,100 on non-orthogonal non-staggered grid systems, based on SIMPLER algorithm in the curvilinear body-fitted coordinates. Effects of wavy heights, lengths, wavy pitches and channel widths on fluid flow and heat transfer were studied. The results show that overall Nusselt numbers and friction factors increase with the increase of Reynolds numbers. According to the local Nusselt number distribution along channel wall, the heat transfer may be greatly enhanced due to the wavy characteristics. In the geometries parameters considered, friction factors and overall Nusselt number always increase with the increase of wavy heights or channel widths, and with the decrease of wavy lengths or wavy pitches. Especially the overall Nusselt number significantly increase with the increase of wavy heights or channel widths, where the flow may become into transition regime with a penalty of strongly increasing in pressure drop. An erratum to this article can be found at     

Measurements and high-speed visualizations of flow boiling of a dielectric fluid in a silicon microchannel heat sink

Experiments were conducted to investigate flow boiling heat transfer to a dielectric fluid in a silicon chip-integrated microchannel heat sink. Twenty-four microchannels, each 389 μm × 389 μm in cross-section, were fabricated into the 12.7 mm × 12.7 mm silicon substrate. High-speed visualizations (at 12,500 frames per second) were performed simultaneously with heat transfer and pressure drop measurements to investigate the physics of flow boiling in parallel microchannel arrays. At low heat fluxes, bubbly flow is dominant, with the bubbles coalescing to form vapor slugs as the heat flux is increased. At high heat fluxes, the flow regimes in the downstream portion of the microchannels are characteristic of alternating wispy-annular flow and churn flow, while flow reversal is observed in the upstream region near the microchannel inlet. Local heat transfer measurements, obtained at three flow rates ranging from 35 to 60 ml/min, show that at lower heat fluxes, the heat transfer coefficient increases with increasing heat flux. The heat transfer coefficient in fully developed boiling is seen to be independent of flow rate in this range. At higher heat fluxes (exceeding 542, 673, 730 kW/m², respectively, for flow rates of 35, 47 and 60 ml/min), this trend is reversed, and the heat transfer coefficient decreases with further increases in heat flux due to partial dryout in some of the microchannels. Heat fluxes at which fully developed boiling is achieved depend on the flow rate. The pressure drop in fully developed boiling increases with increasing heat flux and is independent of flow rate for the test conditions considered in this work.     

Two phase flow and heat transfer characteristics of a separate-type heat pipe

Two phase flow and heat transfer characteristics of a separate-type heat pipe have been studied experimentally and theoretically. The experimental apparatus have the same geometry for the evaporator and the condenser which consist of 5-tube-banks, with working temperature ranges of 80–125°C. The experimental working fluid is dual-distilled water with corrosion-resistant agents. Heat transfer coefficients for boiling and condensation along with heat flux and working temperature are measured at different filling ratio. According to the results of the experiments, the optimized filling ratio ranges from 16 to 36%. Fitted correlations of average heat transfer coefficients of the evaporator and Nusselt numbers of the condenser at the proposed filling ratio are obtained. Two phase flow characteristics of the evaporator and the condenser as well as their influence on heat transfer are described on the basis of simplified analysis. Reasons for the pulse-boiling process remain to be studied.     

Laminar flow heat transfer for simultaneously developing velocity and temperature profiles in ducts of rectangular cross section

Summary A numerical method is used to solve the heat transfer equations for laminar flow in ducts of rectangular cross section with simultaneously developing temperature and velocity profiles, both for constant wall temperature and for constant heat input per unit length of the duct. Like the solutions for a fully developed velocity profile, the Nusselt number for each aspect ratio is found to increase from a limiting value at large distances from the entry plane to a maximum at the entry plane. The results also show a strong effect of the Prandtl number on the heat transfer coefficients with uniform and fully developed velocity profiles representing the upper and lower limits respectively. Comparisons are made with analytical solutions for circular ducts and parallel plates and with experimental data.     

Numerical and experimental study of the heat transfer and fluid flow by thermocapillary convection around gas bubbles

 At liquid–gas or liquid–liquid interfaces thermocapillary or Marangoni convection develops in the presence of a temperature or concentration gradient along the interface. This convection was not paid much attention up to now, because under terrestrial conditions it is superimposed by the strong buoyancy convection. In a microgravity environment, however, it is the remaining mode of natural convection. During boiling in microgravity it was observed at subcooled conditions. Therefore the question arises about its contribution to the heat transfer. Thus the thermocapillary convection was intensively studied at single gas bubbles in various liquids both experimentally and numerically. Inside a temperature gradient chamber, the overall heat transfer around single bubbles of different volume was measured with calorimetry and the liquid flow with PIV and LDV. In parallel to the experiment, a 2-dimensional mathematical model was worked out and the coupled heat transfer and fluid flow was simulated with a CV-FEM method both under earth gravity level and under microgravity. The results are described in terms of the dimensionless Nusselt-, Peclet-, Marangoni-, Bond- and Prandtl-number. Received on 23 August 1999     

Adverse and favorable mixed convection heat transfer in a two-side heated square channel

Experimental heat transfer measurements and analysis for mixed convection in a vertical square channel are presented. Water flow directions are selected such that buoyancy assists or opposes the bulk flow pressure gradient. Unlike most previous experiments with symmetrically heated circular tubes, the present configuration uses an asymmetric heating condition (two sides heated and two sides insulated) and shows significant increase in the Nusselt number for both assisted and opposed flow conditions. Observed heat transfer coefficient distributions are different from the symmetrically heated channels; and this difference in heat transfer coefficient is attributed to the formation of buoyancy driven large-scale flow structures. In general, opposed flow shows higher heat transfer coefficients, and the Nusselt number ratio is observed to increase as Gr/Re or Gr/Re² ratios increase for both assisted and opposed flow conditions. A correlation based on the buoyancy parameter predicts the heat transfer pattern well in both assisted and opposed mixed convection. The range of Reynolds numbers discussed (Re=400–10,000) is of importance for direct numerical simulations and the details provided here can serve as the benchmark data required for complicated buoyancy affected turbulence simulations.     

Seed bubbles stabilize flow and heat transfer in parallel microchannels

Seed bubbles are generated on microheaters located at the microchannel upstream and driven by a pulse voltage signal, to improve flow and heat transfer performance in microchannels. The present study investigates how seed bubbles stabilize flow and heat transfer in micro-boiling systems. For the forced convection flow, when heat flux at the wall surface is continuously increased, flow instability is self-sustained in microchannels with large oscillation amplitudes and long periods. Introduction of seed bubbles in time sequence improves flow and heat transfer performance significantly. Low frequency (∼10 Hz) seed bubbles not only decrease oscillation amplitudes of pressure drops, fluid inlet and outlet temperatures and heating surface temperatures, but also shorten oscillation cycle periods. High frequency (∼100 Hz or high) seed bubbles completely suppress the flow instability and the heat transfer system displays stable parameters of pressure drops, fluid inlet and outlet temperatures and heating surface temperatures. Flow visualizations show that a quasi-stable boundary interface from spheric bubble to elongated bubble is maintained in a very narrow distance range at any time. The seed bubble technique almost does not increase the pressure drop across microsystems, which is thoroughly different from those reported in the literature. The higher the seed bubble frequency, the more decreased heating surface temperatures are. A saturation seed bubble frequency of 1000–2000 Hz can be reached, at which heat transfer enhancement attains the maximum degree, inferring a complete thermal equilibrium of vapor and liquid phases in microchannels. Benefits of the seed bubble technique are the stabilization of flow and heat transfer, decreasing heating surface temperatures and improving temperature uniformity of the heating surface.     

二维水平通道内流动沸腾换热的格子Boltzmann模拟

利用格子Boltzmann方法模拟二维水平通道内水的流动沸腾过程,获得不同壁面过热度下流型特点和不同因素对换热过程的影响规律。结果表明,随着壁面过热度升高,流道内流型依次经历从泡状流、弹状流到反环流的转变,平均热流密度和平均换热系数先增大后减小。入口流速降低会使流道内出现受限气泡流,核态沸腾受到抑制。提高入口流速能够有效促进气泡脱离,壁面平均换热系数随入口流速增大而增大,但增长速率有所减小。减小通道宽度有利于汽化现象发生,核态沸腾得到强化,壁面平均换热系数有所提高。     

Investigation of flow and heat transfer in corrugated-undulated plate heat exchangers

 An experimental and numerical investigation of heat transfer and fluid flow was conducted for corrugated-undulated plate heat exchanger configurations under transitional and weakly turbulent conditions. For a given geometry of the corrugated plates the geometrical characteristics of the undulated plates, the angle formed by the latter with the main flow direction, and the Reynolds number were made to vary. Distributions of the local heat transfer coefficient were obtained by using liquid-crystal thermography, and surface-averaged values were computed; friction coefficients were measured by wall pressure tappings. Overall heat transfer and pressure drop correlations were derived. Three-dimensional numerical simulations were conducted by a finite-volume method using a low-Reynolds number k–ɛ model under the assumption of fully developed flow. Computed flow fields provided otherwise inaccessible information on the flow patterns and the mechanisms of heat transfer enhancement. Received on 5 February 1999     

A self-consistent,physics-based boiling heat transfer modeling framework for use in computational fluid dynamics

Computational Fluid Dynamics (CFD) offers the opportunity to investigate physically and geometrically complex systems with high fidelity. Its applicability to multiphase flow, and particularly boiling heat transfer, is currently limited by the lack of appropriate closure models to describe all relevant phenomena. In this paper, we present an original subcooled flow boiling modeling framework for CFD, which aims at consistently and accurately characterizing the key physics that affect heat transfer at the boiling surface. The new framework introduces a fully mechanistic representation of heterogeneous boiling that improves numerical robustness and reduces sensitivity to closure coefficients. The proposed formulation is inspired by new experimental insight, and significantly extends the existing boiling models by capturing the effects of (i) the microlayer on surface evaporation, (ii) the boiling surface, and (iii) bubbles sliding along the boiling surface. A new statistical treatment of the location and mutual interactions of bubbles on the surface allows for mechanistic prediction of the dry surface area, an important quantity that affects the boiling heat transfer coefficient. This approach lends itself naturally to extension to very high heat fluxes, potentially up to the critical heat flux. An assessment and sensitivity study of the model is presented for a range of mass fluxes (500–1250 kg/m²/s), heat fluxes (100–1600 kW/m²), inlet subcoolings (5, 10, 15 K), and pressures (1, 1.5, 2 bars), demonstrating improved robustness and predictive accuracy at all tested conditions in comparison to traditional heat partitioning approaches, including high heat fluxes, where classic models often fail to converge. Lastly, the framework proposed here should not be viewed as another heat partitioning model, but rather as a general platform that allows incorporation of advanced models for each physical phenomenon considered, leveraging the growing insight generated by modern experimental diagnostics for boiling heat transfer.