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
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
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. . 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.
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 .
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:
Introducing the dimensionless variables as follows
Substituting by the equation (5) into the equations (1-4), we obtain
When the wavelength is large , the Reynolds number is quite small and the equations (7-9) becomes
The dimensionless volume flow rate in the fixed frame of reference is given by
the corresponding dimensionless boundary conditions are
Solving Eq. (12) using the Eqs. (16) and (17), we get
substituting by Eq. (18) into the Eq. (11) and solving Eq. (11) with the boundary conditions of Eqs. (14) and (15), we get
The volume flow rate is given by
From the Eq. (20), we have
On the basis of continuity Eq. (1), we find that, the velocity of cylindrical coordinates of the vapour bubble, can be written as
From the Eq. (22), we can obtain the velocity of vapour bubble radius in a vertical cylindrical tube as the form
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  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. . 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. . 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. .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 .
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).
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.
|Radius of the tube|
|Amplitude of the wave|
|Specific heat of liquid at constant pressure Kg-1J/K|
|Density of vapour and liquid Kg m-3|
|Constant heat addition/absorption|
|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)|
|Constant defined by Eq. 25|
|Non-dimensional heat source parameter|
|Dimensionless volume flow rat|