Simplified analysis of heat transfer through a finned tube bundle in air cooled condenser

Simplified analysis of heat transfer through a finned tube bundle in air cooled condenser

Yanán Camaraza-MedinaÁngel M. Rubio-Gonzales Oscar M. Cruz-Fonticiella Osvaldo F. García-Morales 

Center of Energy Studies and Environmental Technology, Universidad Central de Las Villas, Santa Clara 54440, Cuba

Technical Sciences Faculty. Universidad de Matanzas.Matanzas 44440, Cuba

Corresponding Author Email: 
ycamaraza1980@yahoo.com
Page: 
237-242
|
DOI: 
https://doi.org/10.18280/mmep.050316
Received: 
11 July 2018
|
Revised: 
12 July 2018
|
Accepted: 
27 July 2018
|
Available online: 
30 September 2018
| Citation

OPEN ACCESS

Abstract: 

In this paper. is present a new model for heat transfer calculations during airflow through the finned tubes bank in air cooled condenser systems (ACC). The model was correlated with a total of 736 sets of available experimental data provided by ten authors of recognized prestige in the research area. The air flow bathes a package of finned tubes with an inclination of 450 to 600 with respect to the horizontal line. The studies were performed for ST/SL intervals between 0.4 and 2, air flow velocity of 0.1 to 20 m/s, an ambient temperature of 15 to 43 °C, the incident wind speed over the installation between 0 and 45 km/h, the outer equivalent diameter of the bare tubes between 0.019 and 0.035 m, a height of the fins between 2.7 and 7.5 mm. thicknesses of the fins ranging between 2.3 and 3 mm and a number of fins per unit tube length between 315 and 394. The mean deviation found was 6.5% in 84.8% of the correlated experimental data.

Keywords: 

airflow, heat transfer coefficient, fins tube bank

1. Introduction

It is well known that in a power cycle the heat exchange is a decisive process in its efficiency, since approximately 90% of the heat extracted from the cycle is done through the condensation system. The waste heat from the steam turbine is released to the atmosphere from the cooling system, which, depending on the environmental conditions, makes this exchange from water circulation systems or direct cooling with the environment. [1-3].

In the air cooled condenser (ACC), the refrigerant used is the air. The general theory of heat transfer states that in the transverse cross of a solitary tube or tube package is strongly dependent on the thickness of the film that contours the posterior portion of the tube [4-6, 18]. This problem is known in the literature as a paradox of Schlichting, however, in the available and reviewed sources, there are scarcely a dozen works available that have treated the subsequent bathing of the tubes, according to what is proposed by [5-6].

The first known attempts to solve this problem date back to the year 1927. Several authors of recognized prestige, using combinations of dimensionless groups, empirical coefficients and additional assumptions, took into account approximately the boundary layer thickness, and the drag coefficient generated in the turbulence and the amount of heat transfer generated by the arrangement of the rear rows of tubes.

According to [7-9], the work of Griminson and Zhukauskas, are the reference works of the subject. But in both methods are required to divide the zones of applicability in a dozen regions and despite this complexity the obtained values compute an average error in the best case of 15%. All of this is due to the lack of prediction of the film thickness of the lamellar viscose sublayer [10-12] and to the non-inclusion of its effect on the expressions used.

This problem is further aggravated when the angle of incidence takes components in two directions, that is, it is not perpendicular to the tube package, since a rigorous analysis requires a two-dimensional treatment of the problem [13-15]. Of all the methods known and proposed is the method of Zukauskas who provides the best results in these cases. However, can compute errors of the order of 40%. [7-9].

Given that, there is an appreciable group of experimental quantities reported by different researchers, and the limitations of the current methods in the thermal evaluation of ACC systems by their external portion being known. The authors set themselves the goal of developing an expression that partially solves the constraints of the traditional methods and that allows in an amicable way the obtaining of the average coefficient of external heat transfer in this type of installations.

2. Methods and Validation

2.1 Transverse crossing of a finned tube. Local variation of Nusselt number

In general, transverse batting of solitary cylinders (or a package of these) results in flow separation, which is difficult to handle analytically. Therefore, flows of this type should be studied experimentally or numerically. In fact, flow through cylinders and cylinder packs has been studied experimentally by numerous researchers and several empirical correlations have been developed for the determination of heat transfer coefficient [6-8]. The complicated pattern of the flow through a single finned cylinder greatly influences the heat transfer.

Many investigations have been carried out in order to clarify this problem. In figure 1 are given the average values of the experimental results reported by [2-3]. This figure provides the variation of the local Nusselt number $\mathrm{V} u$ depending on the angle $\varphi$ , on the periphery of a finned cylinder bathed by a cross and perpendicular airflow.

A detailed analysis of figure 1 and the available experimental quantities allows us to establish that in the transversal bathing of a cylinder finned by a turbulent airflow, two basic regions are formed depending on the dimensionless number of. The first in the interval $4 \cdot 10^{3} \leq \operatorname{Re} \leq 10^{5}$  and the second in $10^{5}<\operatorname{Re} \leq 4 \cdot 10^{5}$ .

In the first interval the value of $N u$  starts relatively high at the stagnation point $\left(\varphi=0^{\circ}\right)$ , reaching at this point its maximum value $\left(N u_{M A Y}\right)$ . For $0^{0}<\varphi \approx 60^{\circ}$ , $N u$  presents a decreasing tendency very close to a linear function with negative slope $m \approx-1.15$ . In interval $60^{\circ}<\varphi \approx 80^{\circ}$  the downward trend undergoes a sharp change in slope taking now $-5<m \approx-6$ . In both cases this decrease is caused by the thickening of the laminar boundary layer.

In the point $\varphi \approx 80^{\circ}$ separation of the laminar flow occurs, and therefore it is logical that the number of $\mathrm{N} u$  reach its minimum value, causing this to also coincide in $\varphi \approx 80^{\circ}$  the presence of a positive inflection point. At this point it is true that $N u \approx\left(N u_{M \ell X}\right)^{0.84}$ , For $80^{\circ}<\varphi \leq 180^{\circ}$ , $N u$ present a clear tendency to grow as a result of intense mixing in the region of the separate stream (the wake). This growth follows approximately a positive slope line $m \approx 1.15$ until $\varphi \approx 180^{\circ}$ .where it is true that $N u \approx 0.9 \mathrm{Nu}_{\mathrm{MAX}}$ .

However the second region $10^{5}<\mathrm{Re} \leq 4 \cdot 10^{5}$  its behavior differs from the first. In the first interval, value of $N u$  starts relatively high at the point of stagnation $\left(\varphi=0^{\circ}\right)$ , reaching at this point a high value. $\left(N u_{0}\right)$  but not maximum as happened in the first region. For $0^{0}<\varphi \approx 90^{\circ} \mathrm{N} u$ number present a decreasing tendency very close to a linear function with negative slope $m \approx-1.8$ . The point $\varphi \approx 90^{\circ}$  is where the laminar flow separation occurs, and therefore it is logical that the Nusselt number reaches its minimum value and this point is at once an inflection point. At this point it is true that $N u \approx\left(N u_{0}\right)^{0.92}$ .

For $90^{\circ}<\varphi \leq 120^{\circ}$ , $N u$ number present a clear tendency to growth being asymptotically guided by a positive slope whose value varies in the intervals. 12$\approx m \approx 14$ . This increase in the local number is due to the transition from laminar to turbulent. In $\varphi=120^{\circ}$  is where a maximum $\left(N u_{M A Y}\right)$  and a second point of inflection are located, what is fulfilled $N u_{M A X}=\left(N u_{0}\right)^{1.05}$ .

For the $N u$ number present a pronounced decreasing behavior, not ruled by a single slope, until reaching $\varphi \approx 140^{\circ}$ , at which point the second minimum of the local $N u$ number is presented, being fulfilled that $N u=\left(N u_{0}\right)^{0.98}$ .

This point is precisely where the separation of the flow occurs in the turbulent flow. For $140^{\circ}<\varphi \leq 180^{\circ}$ , a new growth of the local number as a result of the intense mixing in the turbulent region of the wake, this growth presents an asymptotic behavior to a linear function with a positive slope that oscillates in the interval 4.1$\approx m \approx 4.8$ , is also fulfilled in that, $\varphi=180^{\circ}$  that $N u=(16 / 15) \cdot\left(N u_{0}\right)$ .

The discussions given in the previous paragraphs about the local heat transfer coefficients provide insight into the behavior of this at any point in the finned cylinder and for any operant flow regime, however, have little value in the transfer calculations of heat, since in these the average heat transfer coefficient over the entire surface is required. The main author provides a detailed discussion of this problem in the reference [10].

Figure 1. Local Nusselt number variation as a function of the angle in air cross flow over finned cylinder

2.2 Transverse crossing of a bundle of finned tubes

The tubes bank is usually arranged in a row or triangle in the direction of flow, as shown in Figure 2. The outside tube diameter is taken as the characteristic length. The arrangement of the tubes in the bench is characterized by the transverse passage ST. the longitudinal passage SL and the diagonal step SD between the centers of the tubes (see figure 2).

The diagonal step is determined from the following equation:

$S_{D}=\sqrt{S_{L}^{2}+\left(\frac{S_{T}}{2}\right)^{2}}$              (1)

As the fluid enters the bank, the flow area decreases from $A_{T}=S_{T} L$ to $A_{T}=\left(S_{T}-d\right) L$  between the tubes, and as a consequence, the flow velocity increases [17]. In the stepped arrangement the speed can increase further in the diagonal region if the rows of tubes are very close together. In tube banks, the characteristics of the flow are dominated by the maximum velocity $V_{M \alpha x}$ , which is held within the bank, rather than by the approximate velocity V. Therefore, the Reynolds number is defined on the basis of the maximum velocity as:

$\operatorname{Re}_{d}=\frac{\rho V_{\operatorname{Max}} d}{\mu}=\frac{V_{\operatorname{Max}} d}{v}$            (2)

The maximum velocity is determined based on the mass conservation requirement for the incompressible stationary flow. For alignment arrangement, the maximum velocity is in the minimum flow area between the tubes and the mass conservation can be expressed as (see Figure 2):

$\rho V A_{1}=\rho V_{M C \alpha} A_{T}=V_{M \alpha x}\left(S_{T}-d\right)$             (3)

Then from Equation (3) it follows that the maximum velocity will be as:

$V_{M c x}=\frac{S_{T}}{S_{T}-d} V$              (4)

In the stepwise arrangement the fluid approaching through the area A1 (See Figure 2) passes through the area AT and before for area 2AD, as it is wound around the tube in the next row. If it is true that 2AD > AT, the maximum speed still occurs in AT between the tubes and consequently, the VMax in Equation (4) can be used for staggered bundles tube.

Figure 2. Tubes arrangement in a bank

If it happened that $2 A_{D}<A_{T}\left[\text { that is if } 2\left(S_{D}-d\right)<\left(S_{T}-d\right)\right]$, the maximum velocity will be taken in the diagonal cross-sections and. in this case, that maximum velocity is given by:

$V_{M \alpha x}=\frac{S_{T}}{2\left(S_{D}-d\right)} V$          (5)

The conservation of mass can be expressed as (see Figure 2):

$\rho V A_{1}=\rho V_{M \alpha x}\left(2 A_{D}\right)$           (5.1)

Or what is the same.

$V S_{T}=2 V_{M \alpha x}\left(S_{D}-d\right)$           (5.2)

The peculiarities of heat transfer in the transverse flow in a tube bank are very similar to those found in a solitary tube. The variation of heat transfer around a tube bank is determined by the flow pattern, which depends largely on the arrangement of the tubes. Thus in line layout banks, there are two collision points, which results in two thermal transfer peaks. However in banks with staggered arrangement, the process of heat transfer is to some extent similar to that of a solitary tube. Figure 3 provides a comparison of the quotient arising between the local heat transfer coefficient and the average heat transfer coefficient for a solitary tube and for the third row of in-line and stepped arrangements.

Local heat transfer coefficients plotted in Figure 3 correspond to the average value of a measurable group of measurements reported in [2, 11] for $\mathrm{Re}=7 \cdot 10^{4}$ and arrangements $(1,5 x 2) D$ , that is, the transverse passage $S_{T}=1.5 \cdot D$ and the longitudinal passage $S_{L}=2 \cdot D$ .

Figure 3 shows that in the aligned arrays the maximum local heat transfer is located at the point $\varphi=50^{\circ}$ , while the stepped arrangement has two maxima, one at the point of fluid stagnation $\varphi=0^{0}$ and another in the separation of the flow in. It is expected that the presence of two maxima will cause the area under the curve to be larger in the staggered arrangements and therefore, taking into account also the low corrosion and deposition properties of the air, make preferable the stepped arrangement for ACC systems, because with this a higher average heat transfer coefficient is achieved. An inspection of ACC manufacturer’s catalogs confirms that the staggered arrangement is the one used in all marketed equipment.

In any type of arrangement the heat transfer coefficient increases from the first row of tubes to the third, and takes a constant value from the latter. In figure 4 and 5 are plotted for staggered and in-line arrangements, respectively, the average experimental values reported by the authors [9] for the first three rows of a tube package of an ACC. for $\mathrm{Re}=1.4 \cdot 10^{4}$ and one dimension $(2 \times 2) D$ .

The row variation in aligned arrangements does not produce significant changes in local heat transfer coefficient. The authors [8-11], who examined experimentally arrangements with different dimensions, confirm this criterion and for a fixed value of $\operatorname{Re}=3 \cdot 10^{4}$ .The, results obtained confirm that the local heat transfer coefficient follows approximately a common pattern for all the arrangements examined. In figure 6 are plotted the averages available values.

An analysis of 738 sets of available experimental data, reported by 10 authors linked to centers of high prestige [13-15], allowed developing a dimensionally non-homogeneous expression, which was generated by an integral analysis of residues by cross-methods of breaks at intervals, describing with a mean error of the order $\pm 6.5 \%$  in the 84.8% of the available experimental samples.

This equation is given by the following expression:

$\alpha=\frac{\left(T_{T B S}\right)^{0.0064}\left(V_{\max }\right)^{0,6}(e \cdot h)^{0.01}}{\left[0.15\left(S_{T}-d\right)^{0,4}\right] \cdot 0.17 \operatorname{Ln}(F)}\left(\mathrm{W} / \mathrm{m}^{2} \mathrm{K}\right)$           (6)

In Equation (6):

$T_{T B S}$ is the dry bulb temperature, in 0C

$V_{\max }$ is the velocity in the narrowest section of the tube package, in m/s

$e$ is the fins thickness , in mm

$h$ is the fins height, in mm

$F$ is the number of fins per linear meter in length of the finned tube.

The value of $V_{m}$ is determined by the use of the following relations:

if $\quad 2\left(S_{D}-d\right)>\left(S_{T}-d\right) \quad ; \quad V m=\frac{S_{T}}{S_{T}-d} V_{0}$

if $\quad 2\left(S_{D}-d\right) \leq\left(S_{T}-d\right) \quad ; \quad V m=\frac{S_{T}}{2\left(S_{D}-d\right)} V_{0}$          (7)

While the diagonal step is given by:

$S_{D}=\sqrt{S_{L}^{2}+\left(\frac{S_{T}}{2}\right)^{2}}$(m)              (8)

where:

$V_{0}$  It is the airflow velocity at entry to the finned tube pack, in m/s

$S_{T}$ is the transverse step, in m

$S_{L}$ is the longitudinal step, in m

$S_{D}$ is the diagonal step, in m

D is the external diameter of the tubes + fins, in m

d is the outer diameter of the bare tubes (without fins), in m

Figure 3. Local heat transfer coefficient variation for different arrangements of finned tubes

Figure 4. Variation of a local heat transfer coefficient in ACC with staggered arrangement

Expression (6) is valid for $15 \leq T_{T B S} \leq 43^{\circ} C$ , for a refrigerant speed between 0.1 and 20 m/s, as well as ratios $0.4<S_{T} / S_{L}<2$ ,external diameters of plain tubes $0.019<d<0.035 m$ ,a wind velocity incident on installation between 0 and 45 km/h, a height of the fins between 2.7 and 7.5 mm, thicknesses of the fins ranging between 2.3 and 3 mm and a number of fins per unit of pipe length between 315 and 394.

A total of 738 sets of available experimental data found a mean error of the order $\pm 6.5 \%$  to the 84.8% of the correlated experimental samples. Table 1 provides a summary of the main parameters of the experimental quantities available and the mean deviation values found in each case. Figure 7 gives the correlation of available experimental samples with Equation (6).

Figure 5. Variation of a local heat transfer coefficient in ACC with aligned arrangement

Figure 6. Local heat transfer coefficient variation for different dimensions in ACC with aligned arrangement

Figure 7. Correlation of available experimental samples with equation (6)

Table 1. Comparison between Equation (6) and experimental values available for staggered arrangements in ACC systems

Source

Number

of data

Fluid

(ST /SL)

V0

(m/s)

TTBS

(oC)

θ

VW

(km/h)

de (m)

h (mm)

E (mm)

Fins per

meter of

tube

Deviation

percent

(%)

Kumar et al.

(2017)

36

Air

0.4

0.1,16

18, 35

450

0,25

0.019

2.7

2.8

354

3.2, -0.8

44

Air

1

0.3, 12

20,36

600

0,42

0.0254

2.7

2.5

354

5.1, -1.4

38

Air

1.5

0.9, 8

20,32

600

0,36

0.0254

4.7

2.8

354

9.8, 3.2

21

Air

2

0.3, 18

18,30

600

0,24

0.0254

4.7

3.0

354

11.1, 5.4

VDI (2016)

112

Air

0.4

0.5, 19

15,42

600

2,28

0.0254

2.7

2.8

315

5.3, -4.7

14

Air

1

0.1, 11

15,43

450

4,40

0.019

2.7

2.5

315

8.9, 3.7

78

Air

1.5

0.9, 20

18,40

600

0,32

0.0254

2.7

2.8

315

10.2, 4.4

44

Air

2

0.4, 13

18,32

600

0,40

0.0254

2.7

2.8

315

7.7, 1.6

O’Donovan et al.

(2015)

14

Air

1

0.1, 8

17,35

600

0,35

0.0254

2.7

2.8

354

8.3, 0.7

18

Air

1.5

1.1, 7.6

17,30

600

0,35

0.0254

4.7

2.3

354

-2.4, -9.2

11

Air

2

2, 7

15,38

600

0,35

0.019

4.7

2.5

354

11.3, 6.1

Morthensen et al.

(2013)

40

Air

0.4

0.1, 8

18,35

600

0,42

0.0254

4.7

2.8

354

7.5, 4.3

31

Air

1

0.1, 10

18,35

600

0,42

0.035

5.5

2.8

354

8.2, 1.9

27

Air

1.5

0.1, 10

18,35

450

0,42

0.019

5.5

2.5

354

-0.5, -4.8

33

Air

2

0.1, 9

18,35

600

0,42

0.0254

4.7

2.3

354

10.8, 7.2

Heyns- Kröger

(2015)

7

Air

1.5

1.1, 9

26,35

600

0,45

0.0254

4.7

2.5

354

9.2, -5.4

4

Air

2

2.3, 5

28,32

600

0,40

0.019

2.7

2.3

354

11.8, 9.1

Owen-Kröger

 (2013)

8

Air

0.4

0.3, 7

16,34

600

0,28

0.0254

2.7

2.5

354

12.3, 6.1

13

Air

1

1.1, 8

15,31

600

0,28

0.0254

3.8

2.8

354

11.2, 7.4

19

Air

1.5

2.3, 6

15,40

600

0,30

0.0254

42.7

3

354

5.1, 4.4

21

Air

2

1.8, 5

24,36

450

0,35

0.035

7.5

2.8

354

4.5, -6.3

Ingadottir Bara

(2015)

13

Air

1

1.4, 7

21,27

600

6,24

0.0254

4.7

2.3

354

6.1, -4.8

19

Air

1.5

1.2, 9

22,29

600

7,32

0.035

7.5

2.8

315

5.3, 2.2

21

Air

2

1.8, 8.1

24,30

600

4,44

0.035

7.5

2.5

354

5.1, -4.2

 

N. H. Khandy

(2015)

18

Air

1

1.2, 6

18,30

450

0,35

0.0254

4.7

2.5

354

6.1, 1.9

23

Air

1.5

1.8, 5

18,30

450

0,35

0.0254

4.7

2.8

354

9.7, -3.2

11

Air

2

1.1, 7.1

18,30

450

0,35

0.035

7.5

3

354

-3.2, -8.1

Resume

738

 

0.4,2

0.1, 20

15,43

450,600

0,45

0.019,0.035

2.7,7.5

2.3,3

315,354

7.3, 1.1

3. Conclusions

In the study presented, a model has been developed for the determination of the average heat transfer coefficient in transversal cross of airflow to a package of finned tubes used in an ACC system. Could be verified that the expression obtained (6) is valid for $15 \leq T_{T B S} \leq 43^{\circ} C$ , an refrigerant speed between 0.1 and 20 m/s, as well as ratios $0.4<S_{T} / S_{L}<2$ , external diameters of plain tubes $0.019<d<0.035 m$ , a wind speed incident on the installation between 0 and 45 km/h, a fins height between 2.7 and 7.5 mm, fins thicknesses between 2.3 and 3 mm and a number of fins per tube unit length between 315 and 394.

738 sets of available experimental data found a mean error of the order $\pm 6.5 \%$  for 84.8% of the correlated experimental samples. The correlation of the available experimental data with expression (6) was shown in Table 1, and it can be verified that the proposed new model yields satisfactory results in the field of action for which it was developed, reducing in more than half error rates obtained with the use of the Zukauskas or Griminson equations in ACC systems [16-18].

A dimensionally non-homogeneous expression was obtained that responds to the following equation:

$\alpha =\frac{{{\left( {{T}_{TBS}} \right)}^{0.0064}}{{\left( {{V}_{\max }} \right)}^{0,6}}{{\left( e\cdot h \right)}^{0.01}}}{\left[ 0.15{{\left( {{S}_{T}}-d \right)}^{0,4}} \right]\cdot 0.17Ln\left( F \right)}$

Acknowledgment

This work was supported by Doctoral Research Program of Universidad Central “Marta Abreu” de las Villas, Cuba.

Nomenclature

F

number of fins per linear meter of tube

h

Fins height, mm

e

Fins thickness, mm

TTBS

Dry bulb temperature, oC

VMax

Maximum velocity in the minimum flow area, m.s-1

VW

Wind velocity over ACC installation, m.s-1

de

Inside equivalent diameter of tubes, m

ST

Transverse step, m

SL

Longitudinal passage, m

SD

Diagonal step, m

D

external diameter of the tubes with fins, m

V0

Velocity of refrigerant entry to the bundle, m.s-1

Greek symbols

 

$\alpha_{L}$

Local heat transfer coefficient. kg.m-2.s-3.K-1

$\alpha$

Mean heat transfer coefficient. kg.m-2.s-3.K-1

$\varphi$

Angle in air cross flow over finned cylinder, degrees

$\theta$

ACC tubes inclination respect to horizontal line

Subscripts

 

Eq.

Equation

  References

[1] Dorao CA, Fernandhino M. (2017). Dominant dimensionless groups controlling heat transfer coefficient during flow condensation inside pipes. International Journal of Heat and Mass Transfer 112: 465-479. https://doi.org/10.1016/j.ijheatmasstransfer.2017.04.104

[2] Boyko LD, Kruzhilin GN. (1967). Heat transfer and hydraulic resistance during condensation of Steam in a horizontal tube and in a bundle of tubes. International Journal of Heat and Mass Transfer 10(3): 361-373. https://doi.org/10.1016/0017-9310(67)90152-4

[3] Kim SM, Mudawar I. (2013). Universal approach to predicting heat transfer coefficient for condensing mini/micro-channel flow. International Journal of Heat and Mass Transfer 56(1–2): 238–250. http://dx.doi.org/10.1016/j.ijheatmasstransfer.2012.09.032 

[4] Medina YC, Khandy NH, Carlson KM, Fonticiella OMC, Morales OFG. (2018). Mathematical modeling of two-phase media heat transfer coefficient in air cooled condenser system. International Journal of Heat and Technology 36(1): 319-324. https://doi.org/10.18280/ijht.360142

[5] Medina YC, Khandy NH, Fonticiella OMC, Morales OFG. (2017). Abstract of heat transfer coefficient modelation in single-phase systems inside pipes. Mathematical Modelling of Engineering Problems 4(3): 126-131. https://doi.org/10.18280/mmep.040303

[6] Rifert VG, Sereda VV. (2015). Condensation inside Smooth horizontal tubes: Part 1. Survey of the methods of heat-exchange prediction. Thermal Science 19(5): 1769-1789. https://doi.org/10.2298/TSCI140522036R

[7] Thome JR. (2005). Condensation in plain horizontal tubes: Recent advances in modelling of heat transfer to pure fluids and mixtures. Journal of the Brazilian Society of Mechanical Sciences and Engineering 27(1): 23-30. http://dx.doi.org/10.1590/S1678-58782005000100002

[8] Cttani L, Bozzoli F, Raineri S. (2017). Experimental study of the transitional flow regime in coiled tubes by the estimation of local convective heat transfer coefficient. Intl Journal of Heat and Mass Transfer 112: 825-836. https://doi.org/10.1016/j.ijheatmasstransfer.2017.05.003 

[9] Bhagwat SM, Ghajar AJ. (2016). Experimental investigation of non-boiling gas-liquid two phase flow in upward inclined tubes. Exp Thermal and Fluid Science 79: 301–318. https://doi.org/10.1016/j.expthermflusci.2016.08.004

[10] Derby M, Joon H, Peles Y, Jensen MK. (2011). Condensation heat transfer in square, triangular, and semi-circular mini-channels. International Journal of Heat and Mass Transfer 55(1-3): 187–197. https://doi.org/10.1016/j.ijheatmasstransfer.2011.09.002

[11] Borishanskiy VM, Volkov DI, Ivashenko NI, Makarova GA, Illarionov T, Yu LA, Vorontsova LA, Alekseyev NI, Kretunov IOP. (1978). Heat transfer in steam condensing inside vertical pipes and coils. Heat Transfer Soviet Research 10(4): 44–58.

[12] Medina YC, Fonticiella OMC, Morales OFG. (2017). Design and modelation of piping systems by means of use friction factor in the transition turbulent zone. Mathematical Modelling of Engineering Problems 4(4): 162-167. https://doi.org/10.18280/mmep.040404

[13] Kezzar M. Tabet I. Chieul M. Nafir N. Khentout A. (2018). Analytical investigation of heat transfer of solar air collector by Adomian decomposition method. Mathematical Modelling of Engineering Problems 5(1): 40-45. https://doi.org /10.18280/mmep.050106

[14] Zhang W, Du X, Yang L, Yang Y. (2016). Research on performance of finned tube bundles of indirect air-cooled heat exchangers. Mathematical Modelling of Engineering Problems 3(1): 47-51. https://doi.org /10.18280/mmep.030108

[15] Camaraza Y. (2017). Introducción a la termotransferencia. Ed. Universitaria, La Habana, Cube. ISBN: 978-959-16-3286-9

[16] Camaraza-Medina Y, Khandy NH, Carlson KM, Cruz-Fonticiella OM, García-Morales OF, Reyes-Cabrera D. (2018). Evaluation of condensation heat transfer in air-cooled condenser by dominant flow criteria. Mathematical Modelling of Engineering Problems 5(2): 76-82. https://doi.org /10.18280/mmep.050204

[17] Rahim NA, Azudin NY, Shukor SRA. (2017). Computational fluid dynamic simulation of mixing in circular cross sectional microchannel. Chemical Engineering Transactions 56: 79-84. https://doi.org/10.3303/CET1756014

[18] Koh j ZZ, Veerasamy D. (2017). Overview of the use of hydrocarbon refrigerants in air conditioning systems. Chemical Engineering Transactions 56: 1849-1854. https://doi.org/10.3303/CET1756309