Measurement of liquid film thickness in micro square channel

In micro channels, slug flow becomes one of the main flow regimes due to strong surface tension. In micro channel slug flow, elongated bubble flows with the thin liquid film confined between the bubble and the channel wall. Liquid film thickness is an important parameter in many applications, e.g., micro heat exchanger, micro reactor, coating process etc. In the present study, liquid film thickness in micro square channels is measured locally and instantaneously with the confocal method. Square channels with hydraulic diameter of D_h = 0.3, 0.5 and 1.0 mm are used. In order to investigate the effect of inertial force on the liquid film thickness, three working fluids, ethanol, water and FC-40 are used. At small capillary numbers, liquid film at the channel center becomes very thin and the bubble interface is not axisymmetric. However, as capillary number increases, bubble interface becomes axisymmetric. Transition from non-axisymmetric to axisymmetric flow pattern starts from lower capillary number as Reynolds number increases. An empirical correlation for predicting axisymmetric bubble radius based on capillary number and Weber number is proposed from the present experimental data.     

Measurement of the liquid film thickness in micro tube slug flow

Slug flow is one of the representative flow regimes of two-phase flow in micro tubes. It is well known that the thin liquid film formed between the tube wall and the vapor bubble plays an important role in micro tube heat transfer. In the present study, experiments are carried out to clarify the effects of parameters that affect the formation of the thin liquid film in micro tube two-phase flow. Laser focus displacement meter is used to measure the thickness of the thin liquid film. Air, ethanol, water and FC-40 are used as working fluids. Circular tubes with five different diameters, D = 0.3, 0.5, 0.7, 1.0 and 1.3 mm, are used. It is confirmed that the liquid film thickness is determined only by capillary number and the effect of inertia force is negligible at small capillary numbers. However, the effect of inertia force cannot be neglected as capillary number increases. At relatively high capillary numbers, liquid film thickness takes a minimum value against Reynolds number. The effects of bubble length, liquid slug length and gravity on the liquid film thickness are also investigated. Experimental correlation for the initial liquid film thickness based on capillary number, Reynolds number and Weber number is proposed.     

Experimental investigation of vapor bubble growth during flow boiling in a microchannel

Experiments were conducted to analyze flow boiling characteristics of water in a single brass microchannel of 25 mm length, 201 μm width, and 266 μm depth. Different heat flux conditions were tested for each of two different mass flow rates over three different values of inlet fluid temperature. Temporal and spatial surface temperature profiles were analyzed to show the relative effect of axial heat conduction on temperature rise along the channel length and the effect of flow regime transition on local surface temperature oscillation. Vapor bubble growth rate increased with increasing wall superheat. The slower a bubble grew, the further it was carried downstream by the moving liquid. Bubble growth was suppressed for increased mass flux while the vapor bubble was less than the channel diameter. The pressure spike of an elongating vapor bubble was shown to suppress the growth of a neighboring bubble by more than 50% of its volume. An upstream progression of the Onset of Bubble Elongation (OBE) was observed that began at the channel exit and progressed upstream. The effects of conjugate heat transfer were observed when different flow regime transitions produced different rates of progression for the elongation sequence. Instability was observed at lower heat fluxes for this single channel experiment than for similar studies with multiple channels.     

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

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

生物芯片微通道周期性电渗流特性

以双电层的Poisson-Boltzmann方程和黏性不可压缩流体运动的Navier-Stokes方程为 基础,提出二维均匀微通道周期电渗流的解析解. 分析结果表明,周期电渗流速度大 小不但与双电层特性和外电场有关, 而且与流动雷诺数(Re = \omega h^2/\nu )密切相关. 随雷诺数增加,双电层滑移速度下降. 当离开固壁距离增加时,双电层以外区域流动速度快 速衰减,速度滞后相位角明显增加. 研究发现在微通道有波浪状速度剖面. 给出在低雷 诺数时的周期电渗流渐近解,它的速度振幅与定常电渗流速度相同,并具有柱栓式速度分布 形态. 还得到在微通道宽对双电层厚的比值(\kappa h)很小时,Debye-H\"{u}ckel近似 的周期电渗流解, 并与解析解进行分析比较 微通道,双电层,周期电渗流,雷诺数     

Laminar flow and mass transfer in channels with a porous bottom wall and with fins attached to the top wall

 The steady incompressible, viscous, two- dimensional flow of a solution in a channel was considered. The bottom wall was porous and the fins were attached to the top wall. Employing control volume approach, a computer program based on SIMPLE algorithm was developed. Computations were carried out to investigate the effects of the inlet Reynolds number, the fin length, the suction Reynolds number and the slip coefficient on the flow structure and the concentration distribution. It was observed that the thickness of concentration boundary layer increases in the flow direction. The concentration on the porous wall and the concentration boundary layer thickness decrease with increasing fin length, the slip coefficient and the inlet Reynolds number. These results show that fins attached to the upper wall of the channel can be utilized to reduce the concentration polarization and hence improve the effectiveness of the separation process. Received on 24 February 1999     

Binary diffusion and heat transfer in laminar mixed convection channel flows with uniform wall heat flux: extremely thin film thickness

A detailed numerical study has been performed to investigate the combined heat and mass transfer in laminar mixed convection channel flows with uniform wall heat flux. In an initial effort the liquid film on the channel wall is assumed to be extremely thin in thickness. Major dimensionless groups governing the present problem areGr _T,Gr _Mx,Pr,Sc, φ andRe. Results are specifically presented for an air-water system under various conditions. The effects of wall heating flux, the Reynolds number and the relative humidity of the moist air in the ambient on the momentum, heat and mass transfer in the flow are investigated in great detail.     

Numerical study of condensing a small concentration of vapour inside a vertical tube

The purpose of this study is to analyse the combined heat and mass transfer of liquid film condensation from a small steam–air mixtures flowing downward along a vertical tube. Both liquid and gas stream are approached by two coupled laminar boundary layer. An implicit finite difference method is employed to solve the coupled governing equations for liquid film and gas flow together with the interfacial matching conditions. The effects of a wide range of changes of three independent variables (inlet pressure, inlet Reynolds number and wall temperature) on the concentration at exit tube, local Nusselt and Sherwood numbers, film thickness, accumulated condensate rate and temperature are carefully examined. The numerical results indicate that in the case of condensing a small concentration of vapours from a mixture, the resistance to heat and mass transfer by non-condensable gas becomes very intense. The comparisons of average Nusselt number and local condensate heat transfer coefficient with the literature results are in good agreement.     

Experimental studies on heat and mass transfer performance of a coiled tube absorber for R134a-DMAC based absorption cooling system

Absorber is an important component in vapor absorption refrigeration system and its performance has greater influence in overall efficiency of absorption machines. Falling film heat and mass transfer in an absorber is greatly influenced by fluid properties, geometry of heat exchanger and its operating parameters. This paper presents on the results of experimental studies on the heat and mass transfer characteristics of a coiled tube falling film absorber, using 1,1,1,2-Tetrafluroethane(R-134a) and N-N Dimethyl Acetamide (DMAC) as working fluids. The effects of film Reynolds number, inlet solution temperature and cooling water temperature on absorber heat load, over all heat transfer coefficient and mass of refrigerant absorbed are presented and discussed. Normalized solution and coolant temperature profiles and refrigerant mass absorbed along the height of absorber are also observed from the experimental results. The optimum over all heat transfer coefficient for R-134a–DMAC solution found to be 726 W/m²K for a film Reynolds number of 350. The R-134a vapour absorption rate is maximum in the normalized coil height of 0.6 to 1.     

Direct contact condensation of vapor on turbulent, falling liquid film

The problem of condensation of pure vapor on turbulent falling liquid film of the same species is analytically solved. The gradual change in enthalpy of the coolant liquid film in the flow direction is considered to take place in three successive stages. The study brings out the influence of inlet Reynolds number, Prandtl number and degree of subcooling of the coolant on condensation heat transfer coefficients. The heat transfer coefficients predicted from the theoretical analysis are in reasonable agreement with the experimental data available in literature.     

Axial static pressure measurements of water flow in a rectangular microchannel

An experimental investigation of water flow through an aluminum rectangular microchannel with a hydraulic diameter of 169 μm was conducted over a Reynolds number (based upon mean velocity and hydraulic diameter) range from 230 to 4,740. Pressure measurements were simultaneously acquired at eight different axial locations within the channel along with pressure measurements in the inlet and outlet ports. The 27 μm pressure taps were more densely packed near the channel entrance in order to study the developing flow region. The average Poiseuille number for laminar flows was 86.4, which is in excellent agreement with the theoretical value of 86.9. The average critical Reynolds number was found to be 2,370. The limited turbulent friction factor data were in good agreement with the Haaland equation. The inlet to the channel was not well rounded and pressure distributions near the channel entrance show a region of pressure recovery. Entrance length and some minor loss coefficient data were not in agreement with theory, but the cause of these deviations were primarily a function of the inlet geometry and pressure recovery in the microchannel rather than a microscale effect.     

Numerical solutions of pulsating flow and heat transfer characteristics in a channel with a backward-facing step

The incompressible laminar flow of air and heat transfer in a channel with a backward-facing step is studied for steady cases and for pulsatile inlet conditions. For steady flows the influence of the inlet velocity profile, the height of the step and the Reynolds number on the reattachment length is investigated. A parabolic entrance profile was used for pulsatile flow. It was found with amplitude of oscillation of one by Re=100 that the primary vortex breakdown through one pulsatile cycle. The wall shear rate in the separation zone varied markedly with pulsatile flows and the wall heat transfer remained relatively constant. The time-average pulsatile heat transfer at the walls was greater as with steady flow with the same mean Reynolds number.     

A two-phase model for laminar film condensation from steam-air mixtures in vertical parallel-plate channels

Detailed results are presented for laminar film condensation from steam-air mixtures flowing downward in vertical flat-plate channels. The mixture flow is laminar and saturation conditions prevail at the inlet. A fully coupled implicit numerical approach is used that achieves excellent convergence behavior, even for high inlet gas mass fractions. The detailed results include velocity, temperature, and gas mass fraction profiles, as well as axial variations of film thickness, pressure gradient and Nusselt number. The effects of a wide range of changes in the four independent variables (the inlet-to-wall temperature difference and the inlet values of gas concentration, Reynolds number, and pressure) on the film thickness, axial pressure gradient, and the local and average Nusselt numbers are carefully examined. It was found that increases in inlet concentration of noncondensable gas caused significant decreases in the film thickness, local Nusselt number, and axial pressure gradient. An analytical solution for the film thickness and velocity field at the end of condensation path was developed and shown to be the asymptotic value of the numerical results for large distances along the channel.     

Evaporation of pure liquids with increased viscosity in a falling film evaporator

The present study investigated fluid dynamics and heat transfer of viscous pure liquids in a falling film evaporator. This is of special benefit as it avoids mass transfer effects on the evaporation behaviour. Experiments at a single-tube glass falling film evaporator were conducted. It allowed a full-length optical film observation with a high-speed camera. Additionally the evaporator was equipped with a slotted weir distribution device. Test fluids provided viscosities ranging from μ = 0.3 to 41 mPa s. The Reynolds number was between 0.7 and 1,930. Surface evaporation and the transition to nucleate boiling were studied to gain information about the film stability at maximum wall superheat. A reliable database for laminar and laminar-wavy viscous single component films was created. The experimental results show a significant enhancement in the wave development due to the film distribution. A wavy flow with different wave velocities was superposed to the film in each liquid load configuration without causing a film breakdown or dry spots on the evaporator tube. It was found that nucleate boiling can be allowed without causing film instabilities over a significant range of wall superheat.     

Two-phase annular flow and evaporative heat transfer in a microchannel

A physical and mathematical model has been developed to predict the two-phase flow and heat transfer in a microchannel with evaporative heat transfer. Sample solutions to the model were obtained for both constant wall temperature and constant wall heat flux conditions. Results are provided for evaporation rate, liquid film thickness, liquid and vapor phase pressure and temperature distributions. In addition to the sample calculations that were used to illustrate the transport characteristics, computations based on the current model were performed to generate results for comparisons with the experimental results of Qu and Mudawar (2004) where two different mass flow rates of the working fluid were used in the experiment. The comparisons of total pressure drops with the experimental data of Qu and Mudawar (2004) cover the wall heat flux range of 142.71-240 W/cm² with a total channel mass flux of 400.1 kg/m² s and also the wall heat flu range of 99.54-204.39 W/cm² with total channel mass flux of 401.9 kg/m² s. The calculated results from the current model match closely with those of Qu and Mudawar (2004).     

Conjugate forced convection and heat conduction in circular microchannels

The effects of axial heat conduction in the solid walls of microchannels of circular cross-sections are analyzed here. A systematic approach is adopted, with the aim of pointing out the influence of geometrical parameters and of solid wall thermal conductivity on microchannel heat transfer. The reliability of a commonly adopted criterium, based on the so-called axial conduction number, to assess the relevance of axial heat conduction is also discussed. Numerical simulations concern the simultaneously developing laminar flow of a constant property fluid in microchannels of different length, wall thickness and wall material, heated with a uniform heat flux at the outer surface, for different values of the Reynolds number. Moreover, since often in experimental tests the two end sections of the microchannel wall are not perfectly insulated, the effects of heat losses through these sections are also considered. A hybrid finite element procedure, which implies the step-by-step solution of the parabolized momentum equations in the fluid domain, followed by the solution of the energy equation in the entire domain, corresponding to both the solid and the fluid parts, is used for the numerical simulations.     

Heat transfer in droplet-laden turbulent channel flow with phase transition in the presence of a thin film of water

We present results of a numerical study of turbulent droplet-laden channel flow with phase transition. Previous studies of the same system did not take into account the presence of gravity. Here, we do so introducing a thin film of water at the bottom wall and permitting droplets to fall into and merge with it. We treat the carrier phase with the Eulerian approach. Each droplet is considered separately in the Lagrangian formulation, adopting the point–particle approximation. We maintain the film thickness constant by draining water from the bottom wall to compensate for (a) the droplets that fall onto the film and (b) evaporation/condensation. We also maintain on average the total mass of water in the channel by inserting new droplets at the top wall to compensate for the water that has been drained from the bottom wall. We analyze the behavior of the statistically averaged gas and droplet quantities focusing on the heat exchange between the two phases. We increase (a) the initial droplet diameter keeping the same initial droplet volume fraction and (b) the initial number of droplets in the channel keeping their diameter the same. In both parameter studies we find that droplets grow less than in the reference case. In case (a) this is explained by the larger velocity with which they travel to the bottom wall and in case (b) by the lower rate of condensation of vapor due to the presence of neighboring droplets.     

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.     

Comparative study on performance of a zigzag printed circuit heat exchanger with various channel shapes and configurations

Comparative study has been performed with various channel cross-sectional shapes and channel configurations of a zigzag printed circuit heat exchanger (PCHE), which has been considered as a heat exchanging device for the gas turbine based generation systems. Three-dimensional Reynolds-averaged Navier–Stokes equations and heat transfer equations are solved to analyze conjugate heat transfer in the zigzag channels. The shear stress transport model with a low Reynolds number wall treatment is used as a turbulence closure. The global Nusselt number, Colburn j-factor, effectiveness, and friction factor are used to estimate the thermal–hydraulic performance of the PCHE. Four different shapes of channel cross section (semicircular, rectangular, trapezoidal, and circular) and four different channel configurations are tested to determine their effects on thermal–hydraulic performance. The rectangular channel shows the best thermal performance but the worst hydraulic performance, while the circular channel shows the worst thermal performance. The Colburn j-factor and friction factor are found to be inversely proportional to the Reynolds number in cold channels, while the effectiveness and global Nusselt number are proportional to the Reynolds number.     

Convective boiling in a parallel microchannel heat sink with a diverging cross section and artificial nucleation sites

To develop a highly stable microchannel heat sink for boiling heat transfer, three types of diverging microchannels (Type 1, Type 2 and Type 3) were designed to experimentally investigate the effect of different distributions of artificial nucleation sites (ANS) on the enhancement of flow boiling heat transfer, in 10 parallel diverging microchannels with a mean hydraulic diameter of 120 μm. Water was used as the working fluid with mass flux, based on the mean cross section area, ranging from 99 to 297 kg/m² s. The Type-1 system did not contain any ANS; the Type-2 system contained ANS distributed uniformly along the downstream half of the channel; and the Type-3 system contained ANS distributed uniformly along the entire channel. The ANS are laser-etched pits on the bottom wall of the channel and have a mouth diameter of approximately 20-22 μm, as indicted by the heterogeneous nucleation theory. The results of the present study reveal that the presence of ANS for flow boiling in parallel diverging microchannels significantly reduces the wall superheat and enhances the boiling heat transfer performance. The Type-3 system shows the best boiling heat transfer performance.