Open access peer-reviewed chapter

Bubble Dynamics in a Narrow Rectangular Channel

Written By

Xu Jianjun, Xie Tianzhou, Chen Bingde and Bao Wei

Submitted: November 29th, 2017 Reviewed: January 29th, 2018 Published: March 20th, 2018

DOI: 10.5772/intechopen.74608

Chapter metrics overview

842 Chapter Downloads

View Full Metrics

Abstract

It is very important to study the bubble dynamics in order to understand the physical process of boiling heat transfer in a narrow channel. Experimental and theoretical studies on bubble dynamics in a narrow rectangular channel are proposed. The cross section of the narrow rectangular channel is 2 mm × 8 mm. A high speed digital camera is applied to capture bubble behaviors from the narrow side and wide side of the narrow rectangular channel. Bubble growth rate, bubble departure diameter, bubble interface parameter and others are obtained according to the observation. A force balance analysis on a growing bubble is proposed to predict the bubble departure diameter and sliding bubble velocity, and the predicted results agree with the experimental data. Thus, the mechanism of bubble departure, slide and lift-off behavior in a narrow rectangular channel can be explained by the analysis of forces.

Keywords

  • narrow channel
  • bubble dynamics
  • force balance

1. Introduction

Thermal hydraulics and safety analysis are important in nuclear reactor design and operation. There have been active investigations on nuclear thermal hydraulics and safety analysis in different core fuel elements. The rod bundles with grid spacers and plate-type fuel elements are frequently encountered at the aspect of developing new type nuclear reactor. The coolant channels of the plate-type fuel element are composed of some typical narrow rectangular channels, and the flow and heat transfer in a narrow rectangular channel is also widely employed in the compact evaporators and heat exchangers. So it is important to perform the thermal hydraulic study of boiling heat transfer in the narrow rectangular channels.

During the last decades, there were a lot of investigations focusing on two-phase flow and boiling heat transfer in narrow rectangular channels in published literatures [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14]. The results show that two-phase flow and boiling heat transfer in a narrow channel is a very complex physical process. Some researchers think that boiling heat transfer in a narrow channel are different from that in an ordinary sized channel and the boiling heat transfer mechanism in a narrow rectangular channel is not fully understood.

In order to understand the physical process of boiling heat transfer in a narrow channel, and it is important to study the bubble growth, departure and lift-off behavior by using the visualized technology. Over the last half century, a lot of studies of bubble dynamics have been conducted, and this kind of study is significant to understand the boiling heat transfer mechanism. It is the fundamental for understanding and predicting the boiling heat transfer correctly. However, most of these studies have been focused on observation of bubble behaviors in an ordinary sized channel. In our previous work, bubble characteristics in narrow rectangular channel were visually studied, and the results showed that the bubbles departing from the nucleation sites always slid along the heating wall [15]. Some researchers presented the importance of the sliding bubble to enhance heat transfer in nucleate boiling [16, 17, 18, 19, 20], the results showed that turbulence caused by sliding bubbles could be an important mechanism in enhancing heat transfer. The effect of sliding bubbles on boiling heat transfer in a narrow channel will be more important than that in an ordinary sized channel because of the smaller cross section size.

It is known that the analysis of forces acting on a single bubble is a promising approach to sufficiently understand the mechanism of bubble motion. Levy [21] first proposed a model for the prediction of bubble departure diameter based on the balance of forces, including the buoyancy, surface tension force and drag force. Klausner et al. [22] built the models of force balance to predict the bubble departure diameter and lift-off diameter, and there have been active investigations on the analysis of forces balance [14, 23, 24, 25], and these researches implies that bubble interface parameters were required as the input parameters in order to solve this model, including the bubble inclination angle, upstream contact angle, downstream contact angle and bubble contact diameter. It is also important to develop the wall heat flux partitioning model and the interfacial area transport equation for using computational fluid dynamics (CFD) code.

These studies in the above literatures have not provided an adequately detailed understanding of boiling heat transfer in a narrow rectangular channel. The present study showed bubble growth, departure, sliding and lift-off behaviors in a narrow rectangular channel, and the mechanism of bubble departure and lift-off from the heating wall was discussed based on the analysis of forces acting on the bubble.

Advertisement

2. Experimental facilities and measuring system

Figure 1 shows a schematic of the test section of the narrow rectangular channel. The cross section of the narrow rectangular channel was 2 × 8 mm2, which was fabricated in an optical quartz glass. Both the wide side and the gap side of the optical quartz glass were polished to improve the transparency. Bubble behavior was visually observed from the wide side and gap side of the rectangular channel. The observation windows were made of optical quartz glass. The width of the narrow rectangular channel is 8 mm and height is 2 mm. The same experimental apparatus was described to investigate the bubble dynamics in the previous study of the authors [13].

Figure 1.

Schematic diagram of test section.

In the present work, the experiment pressure was atmospheric pressure. The range of the liquid subcooling of the inlet of the test section was 10–40°C. The range of mass flux was 129–706 kg/m2 s. Temperatures at the wall, inlet, and outlet of the test section were registered by 1 mm diameter T-type thermocouples, which were calibrated and a uncertainty of 0.5°C for all thermocouples. The uncertainty in the heat supplied to the test section was less than ±0.9%. The heat loss of the test section was estimated to within 7% of the total of the heating power in the present study. Flow rate was measured by a venturimeter and the uncertainty was less than ±0.8% of the full scale range. The pressure of test section was measured by a pressure transmitter, and the uncertainty was less than ±0.5 kPa of the full scale range.

A high speed digital camera (Micro-NIKKOR lens 200 mm 1:1) was used to capture bubble dynamics. The successive typical bubble images were obtained by using a high speed digital camera. The bubble dynamics can be obtained from the wide side and gap side of the narrow rectangular channel, respectively. A picture of the vernier caliper was captured with the high speed digital camera, and the scale factor was defined by the image of vernier caliper. The real dimension of bubble was obtained with the scale factor. The bubble dynamics was captured with a speed of 3800 fps (frames per second). The 576 × 576 pixel picture corresponds to a 7.2 mm × 7.2 mm recorded field in real dimension. Figure 2 defined the bubble equivalent diameter as following:

π6d3=π6dbxdbz2d=dbxdbz23ifdbx>dbzE1
π6d3=π6dbxdby2d=dbxdby23ifdbx>dbyE2

Figure 2.

Definition of bubble equivalent diameter.

The ratio Rof the bubble contact diameter to the bubble departure diameter is expressed as:

R=dw/d0E3
Advertisement

3. Experimental observation of bubble dynamics

3.1. Bubble growth and departure

A series of images of a typical bubble growing at the nucleation site is obtained by observation from the narrow side of the rectangular channel, as shown in Figure 3. The bubble shape is almost spherical when the bubble is growing at the nucleation sites. The contact diameter between the bubble base and heating wall is observed when the bubble is growing at the nucleation site. The interface parameters of typical bubbles just departing from the nucleation site, such as the bubble inclination angle, upstream contact angle, downstream contact angle, and bubble contact diameter, are measured by observation [14], as shown in Figure 4. The bubble upstream contact angle varies from 44° to 45°, the bubble downstream contact angle varies from 37° to 40°, the inclination angle varies from 5° to 6°, the range of Rin the current experiment is about 0.41–0.5, and the average of Ris about 0.45.

Figure 3.

Visualization of bubble growth and departure from narrow side.

Figure 4.

Bubble contact angle and inclination angle [14].

Figure 5 shows the effect of heat flux on bubble growth at the same nucleation site [14]. The bubble growth rate increases with increasing heat flux. The bubble departure time decreases with increasing heat flux, and so the bubble departure diameter decreases with increasing heat flux. According to Zuber’s model [26], the bubble growth rates can be approximated by a power law curve fit, and the range of the power law varies from about 0.34 to 0.684.

Figure 5.

Bubble growth diameter with time [14].

3.2. Sliding bubble behavior in vertical channel

The images of typical sliding bubbles motion at the same observation window are obtained from the narrow side of the rectangular channel, as shown in Figure 6. After a bubble departs from the nucleation site, it always slides along the heating wall. The typical bubble sliding phenomenon in vertical upflow boiling is observed at low heat flux in the isolated bubble region [13]. During the bubble sliding along the heating wall, the bubble contact diameter between the bubble base and heating wall is observed. The sliding bubble diameter increases with increasing heat flux. The number of sliding bubbles increases gradually with increasing heat flux. The sliding bubbles start to coalescence, and sliding bubbles become the larger after coalescence, kenning et al. [19] considered it as the main model for influence of the thermal boundary layer reformation near the wall due to sliding bubble motion. The phenomenon of bubble lift-off from the heating wall is not observed, which shows the effect of the sliding bubble interaction among the sliding bubbles is also little in the isolated bubble region with low heat flux. It is known that the turbulent effect is strong with increasing number of bubbles for the high heat flux, which maybe results in bubble lift-off from the heating wall.

Figure 6.

Typical sliding bubble with increasing heat flux.

The shape of the small sliding bubble appears spherical, but it is elongated in the direction normal to the heating wall with increasing bubble diameter. When the bubble slides along the heating wall, the upstream contact angle and downstream contact angle of the sliding bubble is almost equal, and so the inclination angle of the sliding bubble approaches to zero. In the present work, the confined bubble growing and sliding motion in the narrow rectangular channel is not observed.

Figure 7 shows the sliding bubble velocity under different bubble diameter with time, as is seen from Figure 7, the sliding bubble velocity increase obviously at the initial moment, but the sliding bubble velocity approach asymptote, ultimately. The influence of the sliding bubble diameter on the velocity is larger. The sliding bubble ultimate velocity increase with increasing bubble diameter, which shows the importance of buoyancy in the flow direction. The velocity of sliding bubble (d0 = 0.389 mm) approaches to the bulk liquid velocity.

Figure 7.

Sliding bubble velocity with time.

3.3. Bubble lift-off behavior

In inclined 45° upward facing upflow boiling, the phenomenon on bubble lift-off from the heating wall is observed from the narrow side of the rectangular channel at the same observation window, as is shown in Figure 8. After a bubble slides some distance, it tends to lift-off from the heating wall, and the buoyancy at the direction normal to the surface is responsible for the bubble lift-off in inclined 45° upward facing heating wall. The bubble velocity will increase when the bubble lifts off from the heating wall. The condensation phenomenon on the interface of the lift-off bubble occurs when the bulk flow is subcooled, which results in the change of bubble shape.

Figure 8.

Single bubble motion in inclined 45° upward facing upflow boiling.

In vertical downflow boiling, the images of bubbles are obtained from the same observation window with the different heat flux, as is shown in Figure 9. After a bubble departs from the nucleation site, it tends to lift-off from the surface. When the bubble lifts off from the heating wall, the shadow of the bubble is observed on the heating wall. The bubble size is small due to the condensation on the interface of the lift-off bubble. The bubble velocity in vertical downflow boiling is small due to the opposite direction of the buoyancy and drag forces acting on the bubble. The local fluid velocity is always higher than that of the bubble velocity, and so the shear lift force is the driving force to promote bubble lift-off from the heating wall.

Figure 9.

Bubble motion in vertical downflow boiling.

Advertisement

4. Analysis of force acting on the bubble

4.1. Balance of forces acting on a single bubble

The forces acting on a single bubble are schematically shown in Figure 10. Different forces acting on a bubble in the directions parallel and normal to a heating wall are analyzed, and the force in the xand ydirections are evaluated as following:

Fx=Fbx+Fsx+Famx+Fqs=mbdvxdtE4
Fy=Fbx+Fsy+Famy+Fsl+Fh+Fcp=mbdvydtE5

where Fbis the buoyancy force, Fsis the surface tension force, Famis the added-mass force, Fqsis the quasi-steady drag force, mb is the mass of the sliding bubble, vxis the velocity of displacement of the center of mass of the bubble in the xdirection, tis the time, Fslis the shear lift force, Fhis the force due to the hydrodynamic pressure, Fcpis the contact pressure force, θis the incline angle of the bubble, αis the upstream contact angle, βis the downstream contact angle, vyis the velocity of displacement of the center of mass of the bubble in the ydirection, dwis the contact diameter, u(y)is the local liquid velocity profile near the wall, and Φis the incline angle of the heating wall.

Figure 10.

Force acting on a bubble.

4.1.1. Buoyancy

The buoyancy acting on a bubble due to the surrounding liquid is expressed as

Fb=ρlρvVbgE6

where ρlis the density of liquid, ρvis the density of the vapor bubble, gis the gravitational acceleration, and Vbis the bubble volume.

4.1.2. Surface tension force

The surface tension forces at xand ydirections are given as follows [22]:

Fsx=dwσlvπαβcosβcosαE7
Fsy=1.25dwσlvπαβπ2αβ2sinα+sinβE8

where σlvis surface tension.

4.1.3. Added-mass force

In order to estimate the added-mass force, Thorncroft and Klausner [23] developed the following expression by solving the inviscid flow problem for a growing sphere in a uniform, unsteady flow.

Fam=1243πρla3dUdtdvxdt+2πρla2UvxȧE9

where the first term is the conventional added-mass force. The second term is the added-mass force due to bubble expansion. In the steady-state flow, dU/dtapproaches zero. Uis bulk liquid velocity, vx is bubble velocity and ais the bubble growth radius.

When the bubble is growing at the nucleation site, the bubble’s center of mass moves away from the heating wall and the added-mass force leads to the unsteady force in the direction normal to the wall. In keeping with convention, the added-mass force is referred to as the bubble growth force. Zeng et al. [27] suggested the growth force can be estimated from:

Fam=ρlπa232CSȧ2+aa¨E10

where CS=20/3.

4.1.4. Quasi-steady drag force

Delnoij et al. [28] suggested that the quasi-steady drag force can be estimated from:

Fqs=12CDρlπa2uvxuvxE11
CD=240Reb0.124Reb1+0.15Reb0.6870.1<Reb10000.44Reb>1000E12
Reb=2atuvxγE13

where γis the liquid kinematic viscosity, uis the velocity of local liquid, and Rebis the bubble Reynolds number.

The local liquid velocity of the mass center of the bubble is estimated by using universal single-phase flow profile. The dimensionless velocity of different region is expressed by the following correlations:

u+=y+y+55lny+3.055<y+302.5lny++5.5y+>30E14
u+=uu=uτw/ρlE15
y+=yuγ=yτw/ρlγE16
u=τw/ρlE17

The wall shear stress can be calculated by the following expression

τw=12CfρlU2E18

where Cfis the friction coefficient, which is calculated by the following expression

Cf=λ/4E19

where λis the friction factor. For a smooth surface, the friction factor is expressed by the following expression

λ=64/ReRe<20000.3163/Re0.252000Re100000E20
Re=UDhγE21

where Dh isthe hydraulic equivalent diameter.

4.1.5. Shear lift force

Saffman [29] developed an expression for the lift force on a particle in a low Reynolds number shear flow. Recently, Mei and Klausner [30] modified the Saffman model for a bubble as follows:

Fsl=12Clρlπa2uvxuvxE22

where Clis the shear lift force coefficient and is defined as

Cl=3.877Gs1/2Reb2+0.014Gs21/4E23

where

Gs=dudyatuvx=1uvx1k+u+E24

4.1.6. Contact pressure force

The contact pressure force due to the pressure difference inside and outside the bubble at the reference point over the contact area is given as follows [22]:

Fcp=πdw242σlvrE25

4.1.7. Hydrodynamic pressure force

Klausner et al. [22] suggested that the hydrodynamic force can be estimated from

Fh=1294ρluvxuvxπdw24E26

The criteria of the bubble departing and sliding from the nucleation site was proposed by Klausner et al. [22]. The bubble remains attached to the nucleation site when the conditions Fx<0and Fy<0. When the condition Fx=0is violated prior to Fy=0, the bubbles will slide along the heating wall, and the criterion for bubble departure from the nucleation site is defined according to the condition Fx=0, at which point the bubble will departure from the nucleation site. When the condition Fy=0is violated prior to Fx=0, the bubbles will lift-off from the heating wall without sliding, and the criterion for bubble lift-off from the heating wall is defined according to the condition Fy=0, at which point, the bubble will lift-off from the heating wall.

4.2. Analysis of mechanism of bubble departure

In order to predict the bubble departure diameter by analysis of force acting on the bubble, the interface parameters of a typical bubble, such as bubble contact diameter, dw, bubble inclination angle, θ, upstream and downstream contact angle, α and β, must be measured. Based on the experimental observation of bubble behavior by Xu et al. [14], dw, θ, α and βcan be set as 0.45db, 45°,40°, and 10°, respectively. The bubble departure diameters under different thermal hydraulic conditions are predicted by solving the balance model of forces. The predictions of the bubble departure diameter are compared by 48 experimental data in both vertical and inclined conditions. Figure 11 shows the predicted result agrees with experimental data.

Figure 11.

Comparison between predicted and measured results [14].

The forces acting on the bubble in the xdirection and ydirection are calculated by using the balance model of forces in order to analyze the mechanism of bubble departing from the nucleation sites. Table 1 shows the value of different forces acting on the bubble just departing from the nucleation sites.

Table 1.

The force acting on the typical bubble (unit: N).

The prediction of results indicates that a bubble will slide along the heating wall before lift-off because the condition Fx=0is violated prior to the condition Fy=0, which is accorded with the visual experimental observation. As is seen from Table 1, the bubble growth force is much less than that of other forces due to the lower growth rate of the bubble, and it is negligible in the current case. The buoyancy force, quasi-steady drag force, and surface tension force are the main forces controlling the bubble departure.

4.3. Analysis of mechanism of sliding bubble motion

According to the experimental results, the inclination angle of the bubble approached 0 during the bubble sliding along the heating wall, and so the xcomponents of the growth and surface tension forces are approximately 0. The balance of forces acting on a single bubble at the xdirection is expressed as

ρlρvVbgsinΦ+12CDρlπa2uvxuvx1243πρla3dvxdt+2πρla2Uvxȧ=mbdvxdtE27

In order to solve the above differential equation, the initial condition is defined as v(tdepart= 0, i.e., vis 0 when the vapor bubble just departs from the nucleation site and tdepart is the departure time of the vapor bubble. In our previous work [31], the sliding bubble velocities for pool and flow boiling in an ordinary size channel were obtained and the predicted results agreed with the experimental data of Maity [32]. However, the above balance of force is (27) was directly used to calculate the sliding bubble velocity in a narrow rectangular channel, it is found that the predicted result is not reasonable, and the force acting on the bubble is need to re-analyzed. The added-mass force is associated with the bubble growth rate and rate of change of bubble velocity. Some researchers showed that flow and heat transfer in a narrow rectangular channel were different with that in an ordinary size channel. In this experiment, the results show that the bubble grows slowly in a narrow rectangular channel and the bubble departure time is longer. This indicates that the size of a narrow channel inhibits bubble growth and departure. While the added-mass force were developed in an ordinary size channel by Thorncroft et al. [23]. Based on 72 experimental data on the sliding bubble velocity in a narrow rectangular channel, the value of the empirical constant Cused for the evaluation of added-mass force is proposed, and the added-mass forces is expressed as:

Fam=C1243πρla3dUdtdvxdt+2πρla2UvxȧE28

When the empirical constant Cis 13.4, the prediction results agree with the experimental results.

Figure 12 shows the sliding bubble velocity against time for vertical flow boiling under the conditions of a bulk velocity of 0.143 m/s. As is seen from Figure 12, the predicted trends of the sliding bubble velocity agree with the experimental data. The sliding bubble velocity increases with increasing time, but the trend of an increase in the sliding bubble velocity decreases gradually as time increases. At about 14 ms, the sliding bubble velocity is higher than the local liquid velocity of the center of mass of the bubble. At about 50 ms, the sliding bubble velocity is almost equal to the bulk velocity. Figure 13 shows the forces acting on the sliding bubble against time. The buoyancy gradually increases with increasing time, but the trend of an increase in the buoyancy is small, this is because the bubble growth rate is low. The buoyancy is the driving force to promote the bubble to slide along the surface in vertical flow boiling. At initial moment, the quasi-steady drag force is the driving force to promote the bubble to slide along the surface. At about 14 ms, the sliding bubble velocity exceeds the local liquid velocity and therefore, the quasi-steady drag force will become the resistance, which prevents the bubble from sliding along the surface. The quasi-steady drag force gradually increases because of an increasing difference between the sliding bubble velocity and the liquid velocity. During this process of the bubble sliding motion, the added-mass force is always a resistance to prevent the bubble from sliding along the surface. As the time increases, since the sliding bubble velocity exceeds the local liquid velocity, and therefore the shear lift force in the ydirection will push the bubble against the wall, maybe this is cause of the bubble sliding along the surface in a vertical channel.

Figure 12.

Sliding bubble velocity for vertical flow.

Figure 13.

Forces acting on the sliding bubble for vertical flow.

Figure 14 shows the sliding bubble velocity against time for in inclined 45° upward facing upflow boiling under the conditions of a bulk velocity of 0.143 m/s. As is seen from Figure 14, the predicted trends of the sliding bubble velocity agree with the experimental data. The trends of the sliding bubble velocity against time are similar to the above in Figure 14. The sliding bubble velocity increases quickly at the initial moment, but the trend of an increase in the sliding bubble velocity decreases gradually as the time increases. The buoyancy in the xdirection is the driving force to promote the bubble to slide along the surface, as shown in Figure 15. At about 5 ms, the sliding bubble velocity is higher than the local liquid velocity of the center of mass of the bubble, and therefore the quasi-steady drag force will become the resistance, which prevents the bubble from sliding along the surface. At initial moment, the added-mass force is a resistance to prevent the bubble from sliding along the surface. At about 9 ms, the added-mass force will become a driving force to promote the bubble to slide along the surface. As the time increases, since the sliding bubble velocity exceeds the local liquid velocity of the center of mass of the bubble, and therefore the shear lift force in the ydirection will push the bubble against the wall. The buoyancy in the ydirection promotes the bubble lift-off from the surface.

Figure 14.

Sliding bubble velocity in inclined 45° upward facing upflow.

Figure 15.

Force acting on the sliding bubble in inclined 45° upward facing upflow.

Advertisement

5. Conclusions

Experimental and theoretical studies on bubble dynamics in a narrow rectangular channel are proposed in this chapter. A high speed digital camera is applied to capture bubble behaviors from the narrow side and wide side of the narrow rectangular channel. Bubble growth rate, bubble departure diameter, bubble interface parameter and others were obtained according to the observation. The bubbles always slide along the heating wall after departing from the nucleation sites. During the process of the sliding bubble motion along the heating wall, the upstream and downstream contact angle is almost equal. The phenomenon on bubble lift-off from the heating wall is not observed in vertical flow boiling with low heat flux in the isolated bubble region. The bubble tends to lift-off from the heating wall after sliding some distance in inclined upward facing upflow boiling and in vertical downflow boiling.

An analysis of force balance on a growing bubble is performed to analyze the mechanism of bubble departure, slide and lift-off behavior in the narrow channel. The bubble growth force is much less than the other forces acting on a bubble, and the buoyancy force, surface tension force, and quasi-steady drag force are the main forces controlling the bubble departure. The buoyancy, quasi-steady drag force, and added-mass force in the flow direction are the main forces controlling the sliding bubble motion along the heating wall. As the time increases, since the sliding bubble velocity exceeds the local liquid velocity of the center of mass of the bubble, and therefore the shear lift force in the ydirection will push the bubble against the wall.

Advertisement

Acknowledgments

This work is supported by the National Natural Science Foundation of China (No.11475161 and 51106142).

Advertisement

Nomenclature

Fxcomposition of forces in the x direction
Fycomposition of forces in the y direction
Fqsquasi-steady drag force
Fbbuoyancy
Fslshear lift force
Fssurface tension force
Famadded-mass force
Fcpcontact pressure force
Fhhydrodynamic pressure force
Vbbubble volume
qheat flux
Gmass flux
ulocal liquid velocity of the center of mass of the bubble
vbubble velocity
ttime
Ttemperature
dwbubble contact diameter
τwwall shear stress
ggravitational acceleration
rbubble radius
a(t)bubble growth rate
Rebbubble Reynolds number
Ubulk velocity,
Cffriction coefficients
λfriction factor
Dhhydraulic equivalent diameter
Rratio of the bubble contact diameter to bubble departure diameter
Cempirical constant
aupstream contact angle
βdownstream contact angle
θinclination angle
ρldensity of liquid,
ρvdensity of bubble
γliquid kinematic viscosity
Φincline angle of the heating wall
xdirection parallel to the heating wall
ydirection normal to the heating wall
bxdirection parallel to bulk flow
bydirection normal to the heating wall
bzdirection perpendicular to flow
ininlet

References

  1. 1. Peng XF, Wang BX. Forced-convection and flow boiling heat transfer for liquid flow through micro-channels. International Journal of Heat and Mass Transfer. 1993;36(14):3421-3427. DOI: 10.1016/0017-9310(93)90160-8
  2. 2. Mishimal K, Mibiki T, Nishihara H. Some characterstics of gas-liguid flow in narrow rectangular duct. International Journal of Multiphase Flow. 1993;19(1):l15-124. DOI: 10.1016/0301-9322(93)90027-R
  3. 3. Tran TN, Wambsganss MW, France DM. Small circular- and rectangular-channel boiling with two refrigerants. International Journal of Multiphase Flow. 1996;22(3):485-498. DOI: 10.1016/0301-9322(96)00002-X
  4. 4. Jiang L, Wong M, Zohar Y. Forced convection boiling in a micro-channel heat sink. Journal of Microelectromechanical Systems. 2001;10:80-87. DOI: 10.1109/84.911095
  5. 5. Kandikar SG. Fundamental issues related to flow boiling in minichannels and microchannels. Experimental Thermal and Fluid Science. 2002;26(2–4):389-407. DOI: 10.1016/S0894-1777(02)00150-4
  6. 6. Piasecka M, Hozejowska S, Poniewski ME. Experimental evaluation of flow boiling incipience of subcooled fluid in a narrow channel. International Journal of Heat and Fluid Flow. 2004;25(2):159-172. DOI: 10.1016/j.ijheatfluidflow.2003.11.017
  7. 7. Lie YM, Lin TF. Saturated flow boiling heat transfer and associated bubble characteristics of R-134a in a narrow annular duct. International Journal of Heat and Mass Transfer. 2005;48(25–26):5602-5615. DOI: 10.1016/j.ijheatmasstransfer.2005.05.013
  8. 8. Lie YM, Lin TF. Subcooled flow boiling heat transfer and associated bubble characteristics of R-134a in a narrow annular duct. International Journal of Heat and Mass Transfer. 2006;49(13–14):2077-2089. DOI: 10.1016/j.ijheatmasstransfer.2005.11.032
  9. 9. Cheng P, Wang GD, Quan XJ. Recent work on boiling and condensation in microchannels. Journal of Heat Transfer. 2009;131(4):1-15. DOI: 10.1115/1.3072906
  10. 10. Yu W, France DM, Wambsganss MW, Hull JR. Two-phase pressure drop, boiling heat transfer, and critical heat flux to water in a small-diameter horizontal tube. International Journal of Multiphase Flow. 2002;28(6):927-941. DOI: 10.1016/S0301-9322(02)00019-8
  11. 11. Chen DQ, Pan LM, Yuan DW, Wang XJ. Dual model of bubble growth in vertical rectangular narrow channel. International Communications in Heat and Mass Transfer. 2010;37(8):1004-1007. DOI: 10.1016/j.icheatmasstransfer.2010.06.023
  12. 12. Chen DQ, Pan LM, Ren S. Prediction of bubble detachment diameter in flow boiling based on force analysis. Nuclear Engineering and Design. 2012;243:263-271. DOI: 10.1016/j.nucengdes.2011.11.022
  13. 13. Xu JJ, He JS, Chen BD, Huang YP, Yan X, Yuan DW. Experimental visualization of sliding bubble dynamics in a vertical narrow rectangular channel. Nuclear Engineering and Design. 2013;261(8):156-164. DOI: 10.1016/j.nucengdes.2013.02.055
  14. 14. Xu JJ, Chen BD, Xie TZ. Experimental and theoretical analysis of bubble departure behavior in narrow rectangular channel. Progress in Nuclear Energy. 2014;77(4):1-10. DOI: 10.1016/j.pnucene.2014.06.002
  15. 15. Xu JJ, He JS, Chen BD, Wang XJ. Visualization of behavior of subcooled boiling bubbles in narrow rectangular slits. Journal of Power Engineering. 2007;27(3):389-392. DOI: 10.3321/j.issn:1000-6761.2007.03.020
  16. 16. Cornwell K. The influence of bubbly flow on boiling from a tube in a bundle. International Journal of Heat and Mass Transfer. 1990;33(23):2579-2584. DOI: 10.1016/0017-9310(90)90193-X
  17. 17. Thorncroft GE, Klausner JF, Mei R. An experimental investigation of bubble growth and detachment in vertical upflow and downflow boiling. International Journal of Heat and Mass Transfer. 1998;41:3857-3871. DOI: 10.1016/S0017-9310(98)00092-1
  18. 18. Thorncroft GE, Klausner JF. The influence of vapor bubble sliding on forced convection boiling heat transfer. Journal of Heat Transfer. 1999;121(1):73-79. DOI: 10.1115/1.2825969
  19. 19. Kenning DBR, Bustnes OE, Yan Y. Heat transfer to a sliding vapour bubble. Multiphase Science Technology. 2000;14(1):75-94. DOI: 10.1615/MultScienTechn.v14.i1.20
  20. 20. Li X, Hollingsworth K, Witte LC. Vapor bubble rise under a heated inclined plate. Experimental Thermal and Fluid Science. 2007;32(2):529-544. DOI: 10.1016/j.expthermflusci.2007.06.003
  21. 21. Levy S. Forced convection subcooled boiling–prediction of vapor volumetric fraction. International Journal of Heat and Mass Transfer. 1967;10(7):951-965. DOI: 10.1016/0017-9310(67)90071-3
  22. 22. Klausner JF, Mei R, Bernhard DM, Zeng LZ. Vapor bubble departure in forced convection boiling. International Journal of Heat and Mass Transfer. 1993;36(6):651-662. DOI: 10.1016/0017-9310(93)80041-R
  23. 23. Thorncroft GE, Klausner JF, Mei R. Bubble force and detachment models. Multiphase Science and Technology. 2001;13(3&4):35-76. DOI: 10.1615/MultScienTechn.v13.i3-4.20
  24. 24. Situ R, Hibiki T, Ishii M, Mori M. Bubble lift-off size in forced convective subcooled boiling flow. International Journal of Heat and Mass Transfer. 2005;48(24–25):5536-5548. DOI: 10.1016/j.ijheatmasstransfer.2005.06.031
  25. 25. Cho Yun J, Yum SB, Lee JH, Park GC. Development of bubble departure and lift-off diameter models in low heat flux and low flow velocity conditions. International Journal of Heat and Mass Transfer. 2011;54(15–16):3234-3244. DOI: 10.1016/j.ijheatmasstransfer.2011.04.007
  26. 26. Zuber N. The dynamics of vapor bubbles in nonuniform temperature fields. International Journal of Heat and Mass Transfer. 1961;2(1):83-98. DOI: 10.1016/0017-9310(61)90016-3
  27. 27. Zeng LZ, Klausner JF, Mei R. A unified model for the prediction of bubble detachment diameters in boiling systems—II. Flow boiling. International Journal of Heat and Mass Transfer. 1993;36(9):2271-2279. DOI: 10.1016/S0017-9310(05)80112-7
  28. 28. Delnoij E, Kuipers JAM, van Swaaij WPM. Dynamic simulation of gas–liquid two-phase flow: Effect of column aspect ratio on the flow structure. Chemical Engineering Science. 1997;52(21–22):3759-3772. DOI: 10.1016/S0009-2509(97)00222-4
  29. 29. Saffman PG. The lift on a small sphere in a slow shear flow. Journal of Fluid Mechanics. 1965;22(2):385-400. DOI: 10.1017/S0022112065000824
  30. 30. Mei R, Klausner JF. Shear lift force on spherical bubbles. International Journal of Heat and Fluid Flow. 1994;15(1):62-65. DOI: 10.1016/0142-727X(94)90031-0
  31. 31. Xu JJ, Chen BD, Wang XJ. Prediction of sliding bubble velocity and mechanism of sliding bubble motion along the surface. Journal of Enhanced Heat Transfer. 2010;17(2):111-124. DOI: 10.1615/JEnhHeatTransf.v17.i2.10
  32. 32. Maity S. Effect of velocity and gravity on bubble dynamics [thesis]. University of California. 2000

Written By

Xu Jianjun, Xie Tianzhou, Chen Bingde and Bao Wei

Submitted: November 29th, 2017 Reviewed: January 29th, 2018 Published: March 20th, 2018