Fluid Mechanics
Volume 2, Issue 2, November 2016, Pages: 28-32

Growth of Vapour Bubble Flow inside a Symmetric Vertical Cylindrical Tube

S. A. Mohammadein1, A. K. Abu-Nab2

1Department of Mathematics, Faculty of Science, Tanta University, Tanta, Egypt

2Department of Mathematics, Faculty of Science, Menoufia University, Shebin El-Koom, Egypt

Email address:

(A. K. Abu-Nab)

To cite this article:

S. A. Mohammadein, A. K. Abu-Nab. Growth of Vapour Bubble Flow inside a Symmetric Vertical Cylindrical Tube. Fluid Mechanics. Vol. 2, No. 2, 2016, pp. 28-32. doi: 10.11648/j.fm.20160202.12

Received: October 27, 2016; Accepted: November 9, 2016; Published: December 12, 2016


Abstract: The paper introduces the incompressible Newtonian fluid with heat transfer in a vertical cylindrical tube under the assumptions of long wavelength and low Reynolds number. The system of mass, momentum, and energy equations are solved analytically. The velocity and temperature field are obtained for two-phase densities. The growth of vapour bubble and its velocity between two-phase densities are obtained for first time under the effect of Grashof number and constant heating source. The obtained results are compared with experiment and Mohammadein at all model with good agreement.

Keywords: Heat Transfer, Two-Phase Flow, Bubbly Flow, Newtonian Fluid


1. Introduction

The study of temperature field and behavior vapour bubble dynamics between two-phase densities [1-3, 5-7, and 9] are most important physical phenomena because of its necessity in many technical processes. The mixed lubrication, chemical metallurgic, oil and gas processes are applications of heat exchange. It plays a major role in the industry of refrigerators, boilers and nuclear reactors which are used for generation of electrical current. The vapour bubble is considered as a finite sink growing inside a mixture (vapor and superheated liquid).There are three stages for bubble growth, inertial, thermal, and diffusion. In the inertial stage, the bubble nucleus depends strongly on the interfacial mechanical interactions such as acceleration, pressure force, and surface tension forces. The inertial stage takes a few milliseconds and thermal phenomena are negligible, therefore, this stage is called isothermal. In the thermal stage, the radius of the nucleus increases and growth becomes mostly dependent on the supply of heat that is consumed to vaporize the liquid on the bubble’s surface. The growth of vapour bubble is affected by heat transfer and pressure changes between two densities. The bubble radius grows within a superheated liquid has been studied by many workers [1-10]. The effect of heat transfer of a Newtonian fluid through asymmetric vertical cylindrical tube is studied by Raoet. al. [10]. The closed form solutions of velocity field and temperature are obtained. The influence of various physical parameters on the flow is observed. The temperature and the heat transfer are discussed through graphs. In Sec. 2, the mathematical model is presented and solving the velocity and temperature distribution, in addition to, the radius of vapour bubble and velocity of bubble are obtained. The discussion and results are presented in Sec. 3. In Sec. 4, the concluded remarks are indicated to the importance of study velocity and temperature distribution around the growing vapour bubble.

1.1. Analysis

Consider the flow of a viscous incompressible Newtonian fluid through a vertical tube. The flow is generated by sinusoidal wave trains propagating with constant speed along the wall of the outer rube. The axisymmetric cylindrical polar coordinate system is chosen such that the coordinate is along the center line of the tube and coordinate along the radial coordinate. The wall of the tube is maintained at a temperature and at the center we have used axisymmetric condition on temperature as in Fig. 1; which depicts the physical model of the problem. Where the radius of the tube is,  is the amplitude of the wave,  is the wavelength and  is the time. The flow is unsteady in the fixed frame . However, in a coordinate system moving with the propagation velocity .

Fig. 1. Physical Model.

The equations governing the flow of a Newtonian fluid are the linear momentum, the conservation of mass and the heat equation are given as follows:

(1)

(2)

(3)

(4)

Introducing the dimensionless variables as follows

(5)

Substituting by the equation (5) into the equations (1-4), we obtain

(6)

(7)

(8)

(9)

When the wavelength is large , the Reynolds number is quite small  and the equations (7-9) becomes

(10)

(11)

(12)

The dimensionless volume flow rate in the fixed frame of reference is given by

(13)

the corresponding dimensionless boundary conditions are

(14)

(15)

(16)

(17)

Solving Eq. (12) using the Eqs. (16) and (17), we get

(18)

substituting by Eq. (18) into the Eq. (11) and solving Eq. (11) with the boundary conditions of Eqs. (14) and (15), we get

(19)

The volume flow rate is given by

(20)

From the Eq. (20), we have

(21)

On the basis of continuity Eq. (1), we find that, the velocity of cylindrical coordinates of the vapour bubble, can be written as

(22)

From the Eq. (22), we can obtain the velocity of vapour bubble radius in a vertical cylindrical tube as the form

(23)

(24)

where,

,

,

(25)

1.2. Discussion and Results

The problem of Newtonian flow through a vertical cylindrical tube is described in Fig. 1, the physical problems is described by the system of equations (1-4), and solved analytically to calculate the temperature distribution in Eq. (18), velocity distribution is derived in Eq. (21), and the velocity and radius of vapour bubble are presented in Eqs. (23), (24). The bubble radius is derived in terms of r, z, and t. Heat equation surrounding the vapour bubbly flow is solved. The equation (18) represents temperature distribution in the mixture surrounding the vapour bubble. The equation (23) describes the vapour bubble flow in a vertical cylindrical pipe. The temperature distribution inside a superheated water under atmospheric pressure (=247kPa and Ts=400 K). The physical values are calculated by Haar [3] as given by Table 1. Moreover, by using Mathematica program (version 6.0), we get the following graphs that demonstrate the effect of the physical parameters on temperature, velocity distribution and the vapour bubble. Velocity distribution in terms of parameter r for two different values of Non-dimensional heat source parameter , is obtained in Fig. 2 (a), it is clearly the velocity distribution increases with the increasing of parameter . The velocity distribution taken the same behavior in the Fig. 2(b) under the effect of Grashof number.This behavior is in agreement with Refs. [10]. In the Fig. 3, we can show the velocity distribution is plotted as a function or and . We note that the velocity distribution  increases with the increasing of the volume rate flow .On the other hand, we can show the temperature distribution in terms of parameter r for two different values of amplitude ratio is obtained in Fig. 4, it is observed that the temperature distribution is proportional with amplitude ratio e. The behavior of vapour bubble radius for different values of two phase densities is shown in Fig. 5. It is observed that the bubble radius is increasing and shifted for the lower values with increasing of. This behavior is agreement with Ref. [7]. In the Fig. 6, the velocity of vapour bubble is plotted as a function of time  with two different values of physical parameter , which decreasing with the increasing of time, and shifted for lower values when the increasing of physical parameter .The velocity of vapour bubble is proportionally inversely with the radius of vapour bubble and shifted for the upper values with increasing the Grashof number , heat source parameter  as shown in Figs. 7, 8, and this is agreement with Ref. [4].The growth of vapour bubble compared with Mohammadein et. al. model is shows in Fig. 9. It is observed that the present model performs at lower values than Mohammadein model with fraction density  and initial velocity .

Table 1. Parameters’ values used in the present problem.

 

Fig. 2a. The velocity distribution is plotted as a function of  with thedifferent values of heat source parameter  ( =5.5, 5).

Fig. 2b. The velocity distribution is plotted as a function of with the different values of Grashof number G (G=3, 3.5).

Fig. 3. The velocity distribution is plotted in 3D as a function of  with thedifferent values of volume flow rate .

Fig. 4. Temperature distribution is plotted as a function of  with thedifferent Values of heat source parameter  ( =5, 5.5).

Fig. 5. The radius of vapour bubble is plotted as a function of time t for two different values of parameter .

Fig. 6. The Velocity of bubble radius is plotted as a function of time twith the different values of

Fig. 7. The Velocity of bubble radius is plotted as a function of bubble radius R with the different Values of heat source parameter ( =5, 5.5).

Fig. 8. The Velocity of bubble radius is plotted as a function of bubble radius R with the different Values of Grashof number G (G=3.5, 3).

Fig. 9. The comparison of present model with Mohammadein et al. [13].

2. Conclusion

The flow of a viscous incompressible Newtonian flow in a vertical tube with observing of vapour bubbles is formulated. The flow is generated by sinusoidal wave trains propagating with constant speed  along the wall of the outer tube. Temperature distribution in the mixture is proportional with different values of amplitude ratio . The velocity distribution in the mixture is proportional with different values of amplitude ratio , Grashof number , heat source parameter . The behavior of vapour bubble radius in the mixture increasing with the increasing values of physical parameter . The temperature, velocity distributions, and behavior of bubble radius in vertical cylindrical tube are obtained as special case when the amplitude ratio tends to zero.

3. Nomenclature

Radius of the tube

Amplitude of the wave

Wavelength

Specific heat of liquid at constant pressure Kg-1J/K

,

Density of vapour and liquid Kg m-3

Constant heat addition/absorption

Thermal conductivity

Radial coordinate
                  z Coordinate z
                  R Instantaneous bubble radius (m)

Instantaneous radial velocity of bubble boundary

Time of bubble growth (s)

Temperature of liquid (K0)

Initial temperature of liquid (K0)

Liquid velocity

Temperature distribution

Amplitude Factor

Constant defined by Eq. 25

Grashof number

Prandtl number

Non-dimensional heat source parameter

Dimensionless volume flow rat

Reynolds number

Wave number
           Subscript Liquid

Vapour

Pressure

References

  1. A. N. S. Srinivas, R. Hemadri Reddy, S. Srinivas, S. Sreenadh. Peristaltic transport of a casson fluid in a channel with permeable walls. International Journal of Pure and Applied Mathematics, Vol. 90, 2014, pp. 11-24.
  2. K. Das. Effect of slip and heat transfer on MHD peristaltic flow in an inclined a symmetric channel. Iranian Journal of Mathematical Sciences and Informatics. Vol. 7, 2012, pp. 35-52.
  3. Alexandr A. A. and Zudin B.Y. Inertial-thermal governed vapour bubble growth in highly superheated liquid. Heat Mass Transfer, Vol. 41, 2005, pp. 855-863.
  4. Batchelor G. K. The theory of Homogeneous turbulence. Cambridge University Press, London, 1959.
  5. Haar L., Callagher J. S., Kell G. S. NBS/NRV, Steam tables, 1984.
  6. Misra, J.C. and Pandey, S.K. "Peristaltic flow of a multi layered power-law fluid through a cylindrical tube", Int. J. Engng. Sci. Vol. 39, 2001, pp.387-402.
  7. S. A. Mohammadein and R.A. Gad El-Rab, The growth of vapour bubbles in superheating water between two finite boundaries, Can. J. Phys, Vol. 79, 2001, pp. 1021-1029.
  8. S.A. Mohammadein and Gorla S. R. The growth of vapor bubble and relaxation between two-phase bubble flow. Int. J. Heat and Mass transfer, Vol. 39, 2003, pp. 97-100.
  9. Mohammadein S.A. and Gouda Sh. A. Temperature distribution in a mixture surrounding a growing vapour bubble. Int. J. Heat Mass transfer, Vol. 42, 2006, pp. 359-363.
  10. Muthu P., B.V. Rathish Kumar, and Peeyush Chandra. Peristaltic motion of micropolar fluid in circular cylindrical tubes. Effect of wall properties Applied Mathematical Modeling, Vol. 32, 2008, pp. 2019–2033.
  11. Prosperetti A. and Plesset M. S. Vapour bubble growth in a superheated liquid. J. Fluid Mech. Vol. 85, 1978, pp. 349-360.
  12. Vasudev C., Rajeswara Rao U., Prabhakara Rao G. and M.V.Subba Reddy. Peristaltic flow of a Newtonian fluid through a porous medium in a vertical tube under the effect of a magnetic field. Int. J. Curr. Scie. Rese. Vol. 3, 2011, pp. 105-110.

Article Tools
  Abstract
  PDF(538K)
Follow on us
ADDRESS
Science Publishing Group
548 FASHION AVENUE
NEW YORK, NY 10018
U.S.A.
Tel: (001)347-688-8931