Impact of dryness degree x on parameters of two-phase flow of water in slick vertical tube.
Aliu, Mehush Rifat ; Terziqi, Avni Kahriman ; Mulliqi, Ismet Sejdi 等
1. INTRODUCTION
The motion of two phase flow during the heat transfer depends on a
number of different characteristics. All characteristics are mainly
related to the mutual hydrodynamic action of phases among themselves and
to the frozen wall as well as to the changes that contribute to the
hydrodynamic flow of phase transition (Sacadura, 1993). For theoretical
analyses made in this study is used the overall linkage between
variables and qualitative interdependences, respectively between unknown
specific operands and those determined by different authors (Gorodecki,
1987), (Michael & Howard, 2004), and (Kutepov & Sterman, 1983).
For this particular case the correlation of the tube wall temperature
and heat transfer coefficient of saturated vapours is determined for
different levels of dryness degree x. From the analysis of this
correlation it will be observed that for greater values of dryness
degree x, the heat transfer coefficient of saturated vapours and its
temperature will be lower. The obtained results can be used for
comparison with experimental ones.
2. MATHEMATICAL MODEL
During the motion of flow containing of vapour-water, the speed of
vapour and water phase is different. In elevating tubes the speed of
mixture of vapour phase is greater than the speed of liquid phase,
whereas in release tubes it is smaller. Hence, to continue with a
dimensional analysis, we avoid the flow speed at entrance and instead we
rely on superficial mass flux, G, along the tube:
G = [??]/[A.sub.pipe] (kg/[m.sup.2]s (1)
This mass flow per surface unit is constant throughout the tube if
the flow is stable. As a result, Reynold's number can be determined
"only for fluid"
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)
Which would be Reynold's number if the entire mass flow was
fluid.
Friction factor for slick tubes is provided in the equation
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)
Whereas Nuselt's number for tubes with slick walls according
to Gnielinski (John, L.IV & John, L.V, 2005) can be calculated with
equation:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)
"D 1.07+127^/ f /8(Pr2'3-i)
For [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]
Accordingly, the physical arguments suggest that the functional
dimensional equation for heat transfer coefficient, should take the
following shape for the flow of saturated vapour in vertical tubes:
[h.sub.fb] = [f.sub.n]([h.sub.lo,G,x,[h.sub.fg],[q.sub.w],[[rho].sub.f][[rho].sub.g],D) (5)
It should be noted that other characteristics of fluids, such as
viscosity and conductivity, are indirectly represented through hlo
according to equation
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)
Functional equation (5) has eight dimensional variables (and one
dimensional variable, x) per unit (m, kg, s, J, K).
These ways are obtained three times more dimensional groups under
x, respectively:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (7)
In fact, the problem is simpler than this since the arguments
related to the pressure gradient show that the quality and rapport of
density can be joined in a single group, known as convection number:
[C.sub.0] = [(1 - x/x).sup.0.8]
[([[rho].sub.g]/[[rho].sub.f]).sup.0.5] (8)
Other group of dimensions in equation (7) is called boiling number:
[B.sub.0] [equivalent to] [q.sub.w]/G [h.sub.fg] (9)
Therefore
[h.sub.fo]/[h.sub.lo] [equivalent to] [f.sub.n] ([B.sub.0],
[C.sub.0]) (10)
According to Kandlikar (John, L.IV & John, L.V, 2005) two
correlations [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] can be
calculated with sufficient accuracy, as follows:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11a)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11b)
Where: "nbd"--means "nuclear boiling
domination" and "cbd" means "convection boiling
domination", [f.sub.0]-orientation factor, F--fluid dependent
parameter for which is recommended F=1 for non corroded tubes for
regimens with dryness degree 0 < x < 0.8. From above correlations
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] the greater value is
selected.
The wall temperature is determined according to the following
equation:
[T.sub.w] = [T.sub.b] + [q.sub.w]/[h.sub.b] (12)
3. CALCULATIONS AND RESULTS
Model for analyzing the dependence of the tube wall temperature and
enthalpy of saturated vapour for different levels of dryness degree x
requires the following entry data: m=0.5 kg/s; Tb=227 [degrees]C;
D=0.0445 m; [q.sub.w]=176000 W/[m.sup.2]; x=5-35 %;
[[mu].sub.f]=0.0001177 kg/ms; [P.sub.r]=0.853; k=0.6439 W/mK; F=1;
[[rho].sub.g]=13.2 kg/[m.sup.3]; [[rho].sub.f]=831.3 kg/[m.sup.3];
[h.sub.fg]=1828000 J/kg;[f.sub.0]=1.
[FIGURE 1 OMITTED]
[FIGURE 2 OMITTED]
[FIGURE 3 OMITTED]
[FIGURE 4 OMITTED]
4. CONCLUSION
As provided in fig.1, convection number [C.sub.0] is in exponential
correlation with the dryness degree x. For lower values of dryness
degree x convection number [C.sub.0] has smaller values whereas for
greater values of dryness degree x convection number [C.sub.0] has
greater values.
Whereas for greater values of x values [MATHEMATICAL EXPRESSION NOT
REPRODUCIBLE IN ASCII] and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN
ASCII] are lower. Whereas for greater values of x the difference
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is lower.
The internal temperature of tube wall is high for greater values of
the dryness degree x and vice versa.
5. REFERENCES
John, L.IV & John, L.V (2005). A heat transfer textbook, third
edition, Cambridge, Massachusetts, U.S.A
Gorodecki, A., (1987). Steam cooling of metallurgical aggregates,
Metallurgy, ISBN5-229-00711-7, Moskva
Michael, J. & Howard, N. (2004). Fundamental of Engineering
Thermodynamics, John Wiley, ISBN 0-471-27471-2, Hoboken
Kutepov, M. & Sterman, S. (1983). Hydrodynamics and heat
transfer in the creation of steam, BBK 31. 31 K95 YDK 621.1.7, Moskva
Sacadura, J. (1993). Knowledge for thermal transfer, TEC&DOC,
ISBN 2-85206-618-1, Paris