Influences of Interfacial Shear in Turbulent Film Boiling on a Horizontal Tube with External Flowing Liquid ()
1. Introduction
The research of film boiling on a horizontal tube was conducted by the pioneering investigator, Bromley [1] . After Bromley’s research, many related researches had been reported. In 1966, Nishikawa et al. [2] analyzed two-phase boundary-layer treatment of free-convection film boiling. The theoretical study had been made of film boiling from an isothermal vertical plate and a horizontal cylinder without considering radiative effects. Jordan [3] investigated the laminar film boiling and transition boiling, and the separated region had also been discussed. Sakurai et al. [4] presented the pool film boiling on a horizontal cylinder with theoretical solutions. The analytical heat transfer model is based on laminar boundary theory including radiation effects. For the forced convection film boiling, Huang et al. [5] investigated the flow film boiling across a horizontal cylinder with uniform heat flux. The numerical results agreed with those experimental data where the wall temperature did not vary a lot around the heater at high heat fluxes, within ±18%.
Laminar film boiling on horizontal tubes has been widely discussed in published literature, and there is also some development on the researches of turbulent film boiling. For example, Sarma et al. [6] presented turbulent film boiling under uniform heat flux condition on a horizontal cylinder. In the research, the assumption of equal shear condition both at the wall and the vapor-liquid interface is reasonable. About the turbulent film boiling on a horizontal isothermal circular cylinder, Sarma et al. [7] presented some theoretical results. The analysis compared the theoretical results with previous experimental results, and found that their results were in a good agreement with the experimental data. In addition, Hu [8] investigated the free convection in turbulent film boiling under the quiescent liquid on a sphere with variable wall temperature.
Even though there were many researches about laminar film boiling and turbulent film boiling on horizontal tubes, there was little publication about the turbulent film boiling with high velocity liquid. Predicting interfacial shear in a turbulent film boiling system under high velocity liquid is not easy. However, the present paper successfully predicts the vapor-liquid interfacial shear by using Colburn analogy. The present study applies the interfacial shear into the forced balance equation, and then combines the forced balance equation with the energy equation and thermal energy balance equation. At last, both the film thickness and Nusselt number are obtained. Besides, the present analysis also includes eddy diffusivity, radiation effects and temperature ratio. Finally, a comparison between the results of the present study and those reported in previous theoretical and experimental studies is provided. It is found that a good agreement exists between the two sets of results.
2. Formulations
Consider a horizontal tube immersed in flowing liquid with the high velocity and saturated temperature Ts. The wall temperature is assumed high enough to occur turbulent film boiling on the surface of the tube, and then a continuous film of vapor runs upward over the tube. The physical model and the coordinate system adopted in the present study are shown in Figure 1.
Boiling under the forced convection, the viscosity component and the buoyancy effect are assumed more significant than the inertia force. Then the force balance equation for the vapor film can be expressed as:
(1)
It is assumed the thickness of vapor film is much thinner than the radius of the tube. And it’s further assumed the turbulent conduction term across the vapor layer is more significant than the convective term, and hence the convective term can be neglected. The energy equation can be expressed as:
(2)
Figure 1. Physical model and coordinate system.
The boundary conditions of energy equation under isothermal condition are as follows:
at
at for (3)
For a pure substance, the thermal energy balance equation of the vapor film can be expressed as:
(4)
Under the high velocity liquid, the equation, which describes heat transfer in the flow around the tube, can be expressed as following equation:
(5)
where C is a constant in flow configuration. The constant value C is listed in Table 1.
According to Colburn analogy, the friction factor can be expressed as the following equation:
(6)
The mean friction coefficient in the stream wise direction may then be written as:
(7)
Furthermore, the local friction can be obtained as:
(8)
The turbulent boundary layer exerts a friction force on the liquid-vapor boundary. The shear stress is estimated by considering the external flowing liquid across the surface of the tube when there is no vapor film on the surface. The local friction coefficient is defined as:
(9)
According to potential flow theory, when the uniform liquid flow of velocity passing a tube, the liquid velocity at the edge of the boundary is as follows:
(10)
Combining Equation (8-10), the local shear stress can be expressed as:
(11)
Incorporating the interfacial vapor shear stress given by Equation (11) into the elemental forced balance equation enables Equation (1) to be rewritten in the following form:
(12)
Table 1. Values of C and n in Equation (5) (*data are used in the present paper).
The forced balance equation Equation (12) yields the following dimensionless equation:
(13)
It’s further assuming the pressure across the boundary layer is constant and the density variation across the boundary layer is given by the following equation:
(14)
The energy equation Equation (2) yields the following dimensionless energy equation:
(15)
The dimensionless boundary conditions of Equation (15) are:
(16)
where the absolute viscosity equation in dimensionless energy equation Equation (15) is expressed as the nitrogen at the saturation temperature corresponding to a system pressure under 1 atm. i. e.
(17)
Besides, the thermal energy balance equation Equation (4) can be rewritten in dimensionless form as follows:
(18)
where the absolute conductivity equation in thermal energy balance equation Equation (18) is expressed as the nitrogen the saturation temperature corresponding to a system pressure under 1 atm. :
(19)
Furthermore, the dimensionless thermal energy balance equation Equation (18) requires the velocity profile in the vapor film. And can be obtained by following equation [8] :
(20)
The boundary condition is:
(21)
The eddy diffusivity distribution presented by Kato et al. [9] is expressed as:
(22)
The heat transfer of turbulent film boiling can be given by the following equation:
(23)
Obviously, the local Nusselt number can be expressed as:
(24)
The mean Nusselt number for the entire surface of the tube can be written as:
(25)
3. Numerical Method
The dimensionless governing Equations (13), (15)-(22) and (24)-(25) subject to the relevant boundary conditions given can be used to estimate, and Nu for the vapor film by means of the following procedures by using C++:
1) Suitable dimensionless parameters, such as, S, NR, Fr and are specified.
2) The boundary conditions of velocity and temperature are as follows:
(26)
3) At the bottom of the tube, , i = 0, the dimensionless film thickness δ+ is also zero
. At the next node, i.e. i = i + 1, the value of is given by, where.
4) Guess an initial value of; substitute Equations (17), (22) into Equations (15) and (16) and then get the
value of.
5) Substitute, Equations (19), (20) and (22) into Equation (18), and get the value of. After calculating, the values of can be gotten, and then substitute the values of into Equation (13).
6) The criterion for the accuracy of is assessed by Equation (13), and the equation can be rewritten in the
following form:. The equation can be written more clearly as the following form:
Furthermore, the error of the numerical calculation is less than, and the precision is
good enough. It can be expressed as the following unequal equation:
If the calculation is a convergence, process the film thickness of next angular position. If the calculation is not a convergence, guess a new thickness and repeat processes (4)-(6).
7) The process above is repeated at the next node position, i.e., and then subsequently at all nodes within the range.
8) The local Nusselt number and mean Nusselt number are then calculated.
9) The flow chart for calculating the vapor thickness is expressed in Figure 2.
4. Results and Discussion
Figure 3 presents dimensionless local vapor velocity. The dimensionless velocity of the tube wall is zero under the no-slip condition. Then, the velocity increases along the y-direction, and because of the shear stress of vapor-liquid interfacial, it becomes larger and larger. Besides, the velocity will increase along with the increase of the angle. And the velocity will reach the maximum value at the top of the tube.
Figure 4 presents the dimensionless temperature distribution for entire angular positions along with the vapor film thickness direction y+. Because assuming the stagnation point is at the bottom of the tube, the figure shows a linear temperature distribution at the bottom of the tube. Besides, the present paper considers the interfacial
Figure 2. The flow chart for calculating the vapor thickness and Nusselt number.
Figure 3. Local velocity distribution in vapor film.
Figure 4. Local temperature distribution in vapor film.
shear with high velocity liquid under the effects of turbulence; thus, when the angular position increases, the non-linear temperature distribution of the profile is becoming more and more significant.
Figure 5 presents the variation of the dimensionless vapor thickness of the tube. Specifically, the film thickness increases continuously from a minimum value at the bottom of the tube as the value of θ increases. It can be seen that the film thickness reaches its maximum value at the top of the tube. Furthermore, the increase of the Froude number, Fr will bring out the increase of the interfacial shear and lead to the increase of the evaporative rate. And the increase of the evaporative rate also leads to an increase of the vapor film thickness.
Figure 6 presents the relationship between the mean Nusselt number and the Froude numbers for five values of Grashof number. The figure shows the results of the forced convection film boiling. A higher Froude number will bring an increase in the mean Nusselt. Besides, the Grashof number Gr is also one of the dominant factors, and therefore increasing the Grashof number will bring out an increase in the mean Nusselt.
Figure 7 provides a comparison between the present results for the mean Nusselt number and the results generated in a previous study [10] . The range of Ra in Pomerantz’s research was. Ra in the present research is defined as modified Rayleigh number,. And the range of Ra of the present study is. The figure shows that the mean Nusselt number calculated in the present study is broadly similar to the results generated by Sarma et al. [8] for the range of Ra of . When Fr = 0, = 0, there is a good agreement between Sarma et al.’s results which is under the condition of quiescent liquid and the present data.
Besides, the figure also compares the present results with the experimental results [9] obtained in an earlier study for laminar film boiling. It can be seen that there is a similar trend between the two sets of results at both
Figure 5. Dimensionless film thickness on tube surface.
Figure 6. Effects of the Fr on mean Nusselt number.
low and mid Ra number. Besides, the increase of the Ra values leads to an increase of the mean Nusselt number. The increase in the Froude number will also bring out an increase in the mean Nusselt number.
5. Conclusions
The following conclusions can be drawn from the results of the present theoretical study:
Figure 7. Effects of the Ra on mean Nusselt number.
1) With the help of Colburn analogy, the present research successfully predicts the shear stress of the vapor- liquid interface in a film boiling system under high velocity of liquid.
2) The increase in the Grash of number will lead to an increase in the mean Nusselt number. Besides, turbulent film boiling under the high velocity of external flowing liquid, the increase in the Froude number will bring out an increase in the mean Nusselt number.
3) It shows a good agreement between the present result which is under the condition of Fr = 0 and the previous studies with free convection.
Acknowledgements
The authors gratefully acknowledge the support provided to this projects by the Ministry of Science and Technology of Taiwan under Contract Number MOST 104-2221-E-019-052.
Nomenclature
a acceleration due to graviton force (m/sec2)
specific heat capacity (J/kg×K)
D diameter of tube, 2R
shear Reynolds,
wall shear parameter,
Fr Froude number,
modified Grashof number,
g acceleration due to gravity (m/sec2)
h heat transfer coefficient, w/(m2k)
hfg latent heat (J/kg)
k thermal conductivity (W/m×k)
dimensionless thermal conductivity,
Local Nusselt number,
Num mean Nusselt number
NR radiation parameter,
Pr Prandtl number,
R Radius of tube (m)
Ra modified Rayleigh number,
Reynolds number of vapor,
S heat capacity parameter,
St Stanton number,
T temperature (K)
temperature ratio,
T+ dimensionless temperature,
u vapor velocity in x-direction (m/s)
shear velocity,
dimensionless velocity,
v velocity normal to the direction of flow (m/s)
x peripheral coordinate (m)
y coordinate measured distance normal to tube surface (m)
y+ dimensionless distance,
Greek Symbols
vapor film thickness (m)
dimensionless film thickness,
absolute viscosity (kg/m×s)
dimensionless absolute viscosity,
kinematic viscosity (m2/s)
density (kg/m3)
shear stress (N/m2)
angle measured from bottom of tube
interfacial shear parameter,
eddy diffusivity for momentum
emissivity
Stefan-Boltzmann constant,
Subscripts
l liquid
s vapor at saturation temperature
v vapor
w tube wall
x x-direction
vapor-liquid interface