Abstract

A numerical study is conducted to clarify fluid-flow characteristics, pressure drop coefficient (), and the average skin friction coefficient () for wing-shaped-tubes bundle with longitudinal fins at downstream side of the tube using the commercial CFD FLUENT software package. The air-side ranges from 1800 to 9700. The tubes are employed with various fin heights () and fin thicknesses () such as 2 mm ≤ ≤ 12 mm and 1.5 mm ≤ ≤ 3.5 mm at the considered range. Results indicate that increases with for all . decreases with for all values for 1800 ≤ ≤ 4200 and then increases for 4200 ≤ ≤ 10000. Lowest values of , , and pumping power (PP) occurred at = 6 mm. Values of for = 6 mm are lower than those of NOF and = 2 mm by about 73 % and 32 %, respectively, at = 4000. decreases with . has negligible effect on . increases for 1.5 mm ≤ ≤ 2.5 mm while decreases with 2.5 mm < ≤ 3.5 mm for the considered range except for = 1850.

1. Introduction

Finned tube heat exchangers are extensively employed in chemical engineering, refrigeration, and HVAC applications. The dominant thermal resistance is usually on the air side in practical applications, and therefore the use of finned surfaces on the air side is very common to effectively improve the overall thermal performance of heat exchangers. Depending on the application, various fin patterns such as plate, louver, convex louver, and wavy fin are used.

There is, of course, an associated increase in due to the increase of friction and the drag contribution from the finite thickness of the interrupted fins. The typical tube geometries used in most of heat exchangers are circular and elliptical ones. A detailed survey on tube heat exchangers with different types of fin configurations is presented as follows.

Chen et al. [1] applied the inverse scheme of the finite difference method in conjunction with the least-squares scheme and experimental temperature data to estimate the air flow and heat transfer characteristics on vertical plate finned tubes heat exchanger in forced convection for various air speed () and fin spacing (). Results showed that the fin efficiency () values decreased with increasing and values. A combined numerical and experimental study of plate finned tubes heat exchangers was examined by Jang et al. [2]. They presented detailed numerical results of and , but experimental results were not enough. A systematic study on heat and friction characteristics of plate finned tubes heat exchangers was experimentally investigated by Wang et al. [3]. Results are limited to the cases with relatively large fin pitch heat exchangers.

Khudheyer and Mahmoud [4] carried out a numerical study on finned tubes heat exchanger to make CFD simulations and validated the results against experimental data. Ten different inlet ’s ranging from 330 to 7000 were simulated for turbulence models of laminar, -epsilon, and SST -omega. Reasonable agreement was found between the simulated and experimental data.

Three-dimensional CFD simulations were carried out to investigate heat transfer and fluid flow characteristics of three-row plain finned tubes heat exchangers using FLUENT software by Ghori and Kirar [5]. Pressure drop characteristics of the heat exchanger were investigated for Re ranging from 330 to 7000. Fluid flow was simulated and results compared for laminar flow model, turbulent flow models -, and SST -omega, with steady-state solvers to determine pressure drop , and temperature fields. Model validation was carried out by comparing simulated values of friction factor () with experimental results by Wang et al. [3]. Reasonable agreement between simulated and experimental data was found, thus proving that the FLUENT software is sufficient for simulating the flow fields in finned tubes heat exchangers.

Yan and Sheen [6] have carried out an experiment to investigate the heat transfer and pressure drop characteristics of finned tubes heat exchangers with plate, wavy, and louvered finned tubes. They found that, at the same , louvered finned heat exchanger shows larger values of compared with those for plate finned tubes one. Wolf et al. [7] studied the air flow and heat transfer performance of a wavy finned tubes heat exchanger by numerical and experimental methods. The air side pressure drop characteristics were successfully modeled using the CFD software FLUENT. The numerical results were validated with the experimental ones and the deviation was within 8%. Tang et al. [8] carried out an experimental and numerical study on the air side performance of finned tubes heat exchangers with various fin configurations, such as crimped spiral, plain, and slit finned tube heat exchanger. It was found that the heat exchanger with the crimped spiral fin has better performance than those for the other configurations. Also, it is found that heat exchanger with slit fins offers the best performance at higher .

Leu et al. [9] conducted a numerical simulation for louvered fin and tube heat exchangers having circular and oval tube configurations. The effects of the geometrical parameters such as louver angle (), louver pitches (), and louver length () were discussed. T’Joen et al. [10] worked on the interaction between the flow behavior (flow deflection and transition to unsteady flow) and the thermohydraulic performance of heat exchanger with an inclined louvered fin design. In this experiment, the impacts of fin pitch (), (), and were discussed in detail. Li and Wang [11] conducted an experimental study on the air side heat transfer and pressure drop characteristics of brazed aluminum heat exchangers, with multiregion louver fins and flat tubes. They found that tends to decrease with increasing and increase with the number of louvers.

Heat transfer and pressure drop characteristics for longitudinally finned tubes banks had been carried out by Sparrow and Kang [12]. The investigated geometrical parameters included the placement of the fins at the tube frontal side, at the rear, and at both. Results indicated that finned tube actually reduced the pressure drop relative to the corresponding unfinned tube. For the same pressure drop, finning yields significantly greater enhancements. At fixed mass flow, greater enhancements were attained with increased-diameter unfinned tubes, but at a high penalty in pressure drop and pumping power.

An experimental investigation had been conducted by Ibrahim and Moawed [13] to clarify heat transfer characteristics and entropy generation for individual elliptical tubes with LF. The investigated geometrical parameters included the placement of the fins at the front of the tube, at the rear of the tube, and at both. The air ranged from about 4.7 × 103 to 3.96 × 104. The analysis of the entropy generation is based on the principle of minimizing the rate of total entropy generation that includes the generation of entropy due to heat transfer and friction losses. Results indicated that the fin position on the elliptical tube has an effect on the results of and irreversibility ratio (). The correlations of average and with were presented.

Torabi and Zhang [14] studied, analytically, the temperature distribution within convective-radiative straight fins with simultaneous variation of , , surface emissivity (), and internal heat generation with temperature. They assumed nonzero convective and radiative sink temperatures for their analysis. Rectangular, triangular, convex, and exponential LF were considered. Calculations are carried out using the well-known, approximate analytical differential transformation method (DTM). From illustrative data, in this paper, one could find that fin with the rectangular profile has the highest comparing to the other three profiles. So, this paper showed the validity and great potential of the DTM for nonlinear problems in science and engineering applications.

There are few reports of implementation of longitudinal finned tubes heat exchanger (LFTH) for wing-shaped tubes. Notably, there are no reports on LFTH enhancement for multirow wing-shaped tubes heat exchangers. Therefore, the aim of the present study is to numerically investigate, by the aid of commercial code FLUENT version 6.3.26, the air flow, skin friction characteristics, and pressure contours through a wing-shaped tubes bundle in staggered arrangement with the placement of various longitudinal fins with different and at the downstream side of the tube. and are 2 mm ≤ ≤ 12 mm and 1.5 mm ≤ ≤ 3.5 mm for the considered range.

2. Description of the Computational Model

The numerical model is used to study the hydrodynamic and thermal performance of a staggered wing-shaped tubes bundle with LF mounted at the rear side of the tubes. The computational model considered, in the present work, is a simple design of air to water heat exchanger. The heat exchanger domain consists of three connected rectangular duct sections, entrance, tube bundle, and exit, as shown in Figure 1. The entrance section dimensions are 305 mm, 305 mm, and 320 mm as height, width, and length, respectively. The exit section is similar to the entrance one. The tube bundle section is fixed between the two above-mentioned sections.

Schematic plane of the numerical domain for the crossflow over the tube bundle, shown in Figure 2, consists of 22 wing-shaped tubes, which are distributed through three successive rows in addition to four half dummy ones. The cross-sectional dimensions of wing-shaped tube, drawn from 1 mm thick, 22.5 mm outer diameter circular copper tube with 305 mm length, are shown in Figure 3. LF of copper sheet having a width of 305 mm and varied and such as 2 mm ≤ ≤ 12 mm and 1.5 mm ≤ ≤ 3.5 mm are welded on the downstream side of the tube surface. Longitudinal and transverse tube pitches of 37 mm were kept constant. Hot air is forced to flow over the external cooled surface of the tubes with different . Flowing air velocity () is varied from 1.33 to 7 m/s, yielding ranging from 1.8 × 103 to 9.7 × 103. Tubes surface temperature is kept constant.

3. Numerical Investigation

3.1. Problem Description and Boundary Conditions

Since the tube length is much greater than its equivalent diameter, the flow across the tube bundle is considered two-dimensional flow one. The geometry of the numerical model including the boundary conditions is shown in Figure 2. The numerical solution is carried out by solving the governing equations of mass, momentum, and energy under the following assumptions; the flow is incompressible, steady, and turbulent. The fluid properties are constant; also, the effects of buoyancy forces and radiation are neglected. Based on these assumptions the governing equations can be set aswhere is a tensor indicating 1 and 2, is the viscous stress tensor, and is the fluid effective thermal conductivity.

Commercial CFD software FLUENT 6.3.26 is used to solve the governing equations. RNG - turbulence model is used to simulate the flow fluctuations because it is reported to be more responsive to the effects of flow strain rate, which significantly improves the accuracy for rapidly strained flows, such as the case in reactive swirling flows. Also the RNG - model shows better performance than the standard - model in the prediction of heat transfer as examined by Hu and Sun [15]. The RNG - model has shown an excellent agreement between numerical and experimental results for an isothermal flow over backward facing step (Yakhot et al.) [16]. The near wall region is treated using the enhanced wall method, which is capable of predicting the flow separation and recirculation behind bluff bodies, FLUENT [17].

3.2. Mesh Generation and Discretization

The geometry and mesh of the computational model described in the previous section are generated separately using GAMBIT 2.4.6. Quad/Map, Tri/Pave, and boundary layer meshing scheme is used in the present study with refining the mesh near walls and beside sharp edges.

Figure 4 shows the configuration of the computational mesh. Over the computational grid, the conservation equations as well as the turbulence equations are discretized using the second- order upwind scheme. The velocity pressure coupling is solved using SIMPLE algorism. The solution is considered converged when the scaled residual of the energy equation reaches 10−7 and the scaled residuals of other equations reach 10−4.

To test the dependence of the numerical results on the grid density, calculations are undertaken with different mesh densities in the and directions as recommended by FLUENT [17]. The first grid points adjacent to the walls were kept at values between 1 and 5. The grid sensitivity analysis is carried out, mainly, to obtain grid independent temperature distributions. The computational results of for the tube bundle are varied to give about 2.1% deviation when the number of grids ranges from 250,000 to 420,000. On the other hand, when the grid size ranges from 300,000 to 420,000, only 0.12% variation in is obtained as shown in Figure 5. It is clear that the effect of grid size on the computed results diminishes for grids of 300,000 nodes. Therefore, the computing time can be, drastically, cut down to save computing time without a penalty of poor accuracy.

4. Data Reduction

4.1. Fluid Flow Considerations

The pressure drop coefficient, , as defined below in (2) represents the ratio of the total pressure drop of the moving air through the tube bundle to its dynamic pressure:where is the pressure drop through the tube bundle and is the air film density and is the inlet air velocity.

and skin friction coefficient are denoted aswhere is the wall shear stress (N/m2).

5. Results

5.1. Validation of the Computational Model

The numerical results of the present study are validated with the corresponding experimental ones obtained by Sayed Ahmed et al. [18]. Figure 6 shows the validation of the experimental and numerical values of versus for wing-shaped tubes without fins (NOF). The average estimated error between the numerical and experimental results, at considered , does not exceed 9%.

5.2. Effect of and at Constant ( = 1.5 mm) on Fluid Flow Characteristics of Heat Exchanger Performance

From the earliest study, it has been found that the wing-shaped tubes heat exchangers have higher performance than those with circular and elliptical ones (Sayed Ahmed et al.) [18, 19]. Further performance improvement requires the elimination of heat transfer deficiency in certain regions over the wing-shaped tubes. This suggests the placement of extended surface and/or protrusions in the deficient regions, where the separation of flow makes the heat transfer, mainly, by convection.

Fins is conducting heat outside of this region to the main flow where forced convection is dominant. Thus, fins act as a heat transfer cascade where the heat transfer conduction through the fins is followed by convection to the main flow downstream.

5.2.1. Air Flow Characteristics

Changing for the wing-shaped tubes could affect the location at which separations and wakes were formed as well as the turbulence generation rate in the flow passing over the tubes. This in turn would have a direct effect on the resultant and on the heat exchange during the air cooling process. So, results will be portrayed by presenting the flow streamlines, velocity, and turbulent kinetic energy contours.

Flow-path-lines and velocity contours across the wing-shaped tubes of the studied bundle arranged for different fin heights are illustrated in Figures 710. Up to 12 mm fin height as shown in Figure 7, at = 1.8 × 103, for the arrangement of different fin heights there are three flow separation zones, two on the lateral sides and one at the rear surfaces of the tubes. This region is considered as a dead zone, where the separation of the flow makes the heat transfer mainly by natural convection. Flow separations occurred, as a result of travelling of the boundary layer far enough against an adverse pressure gradient that makes the velocity of the boundary layer fall, almost, to zero. On the other hand, flow separation can occur as a result of the increase of the pressure drag caused by the pressure differential of the flowing air between the frontal and rear surfaces of the tube.

On the other hand, for high , as shown in Figure 8, separations and secondary flows intensity occurred, mainly, on the rear surfaces of the tubes. It is seen from Figure 8 that, with increasing fin height, eddies and vortices at lateral tube surfaces disappeared. This is due to the fact that, for high , the flow becomes more energetic enabling the boundary layer to travel further along with the tube surface before separation occurs resulting, only, in a narrower wake at the rear surface of the tube.

As seen from Figures 9 and 10, the air flow is strongly accelerated in the passages between tubes. With increasing of fin heights and/or the flow is strongly guided between the tubes and the deficient flow separation zones are delayed and/or disappeared. Also, the levels of turbulence and formation of vortices through the passages of the bundle are increased.

Figures 11 and 12 present the turbulent kinetic energy contours through the wing-shaped tubes of the bundle for both low and high , which is defined as the mean kinetic energy per unit mass associated with eddies in turbulent flow. Physically, the turbulent kinetic energy produced by friction and buoyancy forces or through external forcing at low-frequency Eddie scales (integral scale) is characterized via root-mean-square method (RMS). Turbulence kinetic energy is transferred down the turbulence energy cascade and dissipated by viscous forces at the Kolmogorov scale (smallest scales in the spectrum that form the viscous sublayer range). It is found that the turbulent kinetic energy is larger in magnitude in the case of higher . This, in turn, enhances the fluid flow characteristics. It is, also, seen from the figures that there is a direct proportionality between and the values of turbulent kinetic energy () but the numerical results show that the values of at the exit section of the bundle are always greater than those at the inlet one for considered range, so the bundle works as a turbulence generator which leads to an increase in the heat transfer between tubes and flowing air. With the increasing of the values of increase at the fin tip leading to the enhancement of heat transfer over the following rows of the bundle. So, the fin acts as a turbulent generator.

5.2.2. Air Flow Pressure Drop through the Bundle

Figures 13 and 14 demonstrate the static pressure contours through the bundle. It can be seen that the pressure has the highest values at the stagnation point on the frontal portion of each tube, at which the flow velocity tends to be zero. As the flow passes over the tube surface, the pressure decreases until reaching its lowest value at the tube lateral surface.

Figure 15 shows the variation of versus at different . It is seen from the figure that increases with the increase of for different . Figure 16 shows the variation of versus various at different values of . decreases with for 2 mm ≤ ≤ 6 mm, while there is no significant change of for 6 mm < ≤ 12 mm except for slight increasing of at 6 mm < ≤ 12 mm for = 9700. Lowest values of occurred for arrangements of = 6 mm.

Figure 17 shows the variation of versus at different . It is seen from the figure that decreases with for NOF. has the highest values for all studied at = 2000 and decreases with the increase of in the range of 1800 ≤ ≤ 4200 and then increases for 4200 ≤ ≤ 10000 for all tested fins. Lowest values of occurred at = 6 mm. The values of for = 6 mm are lower than those of NOF and = 2 mm by about 73% and 32%, respectively, at = 4000.

Figure 18 shows the variation of versus various at different values of . decreases with for 2 mm ≤ ≤ 6 mm, while increases with 6 mm < ≤ 12 mm for the considered range.

These results can be attributed to the fact that the overall drag consists of two combined parts; one presents the pressure drag while the other presents the friction drag. The friction drag is more dominant than the pressure drag at the lower , thus leading to higher pressure drop while the opposite is true at higher . In the case of high , the effect of viscosity is less dominant and the total drag is rather affected by the inertia forces. As the flow tends to shift towards more turbulent region, the separation region, also, travels farther downstream, reducing the size of the wake and the magnitude of the pressure drag.

5.2.3. Skin Friction Coefficient ()

The distribution of average skin friction coefficient around wing-shaped tubes bundle is studied. Figures 19 and 20 show the effect of on average at different for wing-shaped tubes bundle. decreases as increases. It can be observed that there is no significant difference in the values of for different arrangements at considered .

This phenomenon may be explained as follows. For the baseline case NOF, the air flow resistance from the inlet to outlet of the heat exchanger arises from the local resistance of the wing-shaped tubes, which plays dominant part of the total pressure drop. When the LFTHs are installed, the air flow resistance will arise from the tubes and the local resistance of the LFTH. LFTH brings about additional form drag and delays the separation of boundary layer from the wing-shaped tubes and decreases the wake region behind the tube so that the form drag from the tubes decreases. Therefore, the LFTH with = 6 mm offers lower form drag than that with = 12 mm and that of NOF. Thus, the case of = 6 mm increases the heat transfer and decreases compared with the base case of NOF.

5.2.4. Selection of Best Arrangement for Lowest Pumping Power (PP)

One of the goals of this study is to predict optimal that would minimize the frictional losses of the flow. A typical plot is provided to compromise between and the pumping power (PP). A relationship (4) was obtained on multiplying the pressure drop and the volumetric flow rate aswhere is pressure drop (Pa). is air volume flow rate (m3/s).

Figure 21 shows the variation of PP versus at different . It is seen from the figure that PP increases with the increase of for arrangements with different . Figure 22 shows the effect of on PP at different . At certain , PP decreases with from 2 mm ≤ ≤ 6 mm, while the variation of PP is negligible for 6 mm < ≤ 12 mm. The figure shows, also, that the lowest values of the pumping power are achieved for arrangements of = 6 mm for the considered range.

This is due to the fact that the wing-shaped tubes bundle arrangement promoted turbulent mixing and lengthened the air flow path through the bundle. The size and the strength of the turbulence level as well as the reversed flow region are affected by and variations.

5.3. Effect of at Constant ( = 6 mm) on and

Figures 23 and 24 show the effect of on and for different at constant = 6 mm. It is clear from Figure 23 that has negligible effect on and the highest values of are obtained for = 1850 while the lowest ones belong to = 9700. Figure 24 shows a variation in values for the considered range, where increases for 1.5 mm ≤ ≤ 2.5 mm while decreases with 2.5 mm < ≤ 3.5 mm for the considered range except for = 1850 where there is no significant change of for 1.5 mm ≤ ≤ 3.5 mm.

6. Correlations of the Results

Correlations for and based on the numerical results obtained in terms of Re and / are predicted by (5) and (6), using power regression, as follows:where is the equivalent circular diameter.

Table 1 shows all constants with maximum percentage difference (max. diff. %) between the calculated and obtained results for (5) and (6). The excellent values confirm the statistical goodness of the fit. The obtained correlations are applicable for 1.8 × 103 ≤ 9.7 × 103.

7. Conclusions

Improvement of wing-shaped tubes bundle performance requires the elimination of fluid flow deficiency in certain regions over the tubes. This suggests the placement of extended surface in the deficient regions, where the separation of flow makes the heat transfer, mainly, by natural convection. A numerical study has been conducted to clarify fluid flow characteristics, pressure drop distributions, and over a wing-shaped tubes bundle in staggered arrangement with the placement of LF on downstream side of the tube. ranged from 1.8 × 103 to 9.7 × 103. The tubes bundles are equipped with various fins having 2 mm ≤ ≤ 12 mm and 1.5 mm ≤ ≤ 3.5 mm at considered range. The flow pattern around the staggered wing-shaped tubes bundle is predicted using the commercial CFD FLUENT 6.3.26 software package. Correlations of pressure drop coefficient and skin friction coefficient ()  in terms of , design parameters, for the studied cases are presented. Comparisons between the experimental and numerical results of the present study and those, previously, obtained for similar available studies showed good agreements. Results indicate that the attaching of LF to the tubes surfaces downstream acts as a heat transfer cascade where heat transfers firstly by conduction through the fins and then followed by convection to the main flow downstream. The following conclusions are preferred:(i) increases with of for LF with different . decreases with for 2 mm ≤ ≤ 6 mm, while there is no significant change of for 6 mm < ≤ 12 mm except for slight increasing of at 6 mm < ≤ 12 mm for = 9700. Lowest values of occurred for arrangements of = 6 mm.(ii) has the highest values for all studied at = 2000, and decreases with the increase of in the range of 1800 ≤ ≤ 4200 and then increases for 4200 ≤ ≤ 10000 for all tested fins. Lowest values of occurred at = 6 mm.(iii) decreases with for 2 mm ≤ ≤ 6 mm, while increases for 6 mm < ≤ 12 mm with considered . The values of for = 6 mm are lower than those of NOF and = 2 mm by about 73% and 32%, respectively, at = 4000.(iv) decreases as increases for all studied fins. The use of LF with = 6 mm results in decreasing average compared with these of the case of NOF.(v)PP increases with the increase of for arrangements with different . PP decreases with from 2 mm ≤ ≤ 6 mm, while the variation of PP is negligible for 6 mm < ≤ 12 mm. The lowest values of PP are achieved for arrangements of = 6 mm for the considered range.(vi) has negligible effect on and the highest values of are obtained for = 1850 while the lowest ones belong to = 9700.(vii) increases for 1.5 mm ≤ ≤ 2.5 mm while decreases with 2.5 mm < ≤ 3.5 mm for the considered range except for = 1850, where there is no significant change of for 1.5 mm ≤ ≤ 3.5 mm.(viii)Correlations for and based on the numerical results are obtained in terms of and .The Resultant Point of View. The above-mentioned discussion showed that the wing-shaped tubes bundle heat exchanger with = 6 mm has the lowest values of , , , and PP and hence the best performance comparing with the other bundles.

Symbols

Alphabet: Upper Case
PP:Pumping power, W
:Total outer surface area of the tubes, m2
:Equivalent circular diameter, m
:Fin pitch, m (fin spacing + fin thickness)
:Fin spacing, m
:Pressure drop coefficient, )/()
:Coefficient of determination ranged from 0 to 1
Re:Reynolds number, ()/μ
:Temperature, K
:Velocity, m/s
:Longitudinal tube pitch, m
:Transverse tube pitch, m
:Slit width, m
:Louver length, m
:Louver pitch, m
:Surface emissivity parameter, dimensionless.
Alphabet: Lower Case
:Specific heat at constant pressure, J/kgK
:Heat transfer coefficient, W/m2K
:Thermal conductivity, W/mK
:Mass flow rate, kg/s
:Tube thickness, m
:Axial coordinate
:The normal distance to the tube surface
:Dimensionless normal distance to the tube surface
:Fin height, m.
Greek Symbols
Absolute viscosity, Pas
:Louver angle, °
:Density, kg/m3
:Pressure drop across the bundle, Pa
:Irreversibility ratio
:Fin efficiency
:Fin thickness, m.
Subscripts
Air
:Circular
:Exit
:Film or fin
:Inlet or row number
:Outer or outlet
:Height
:Longitudinal.
Abbreviations
RNG:Reynolds normalized group
LFTH:Longitudinal finned tubes heat exchanger
NOF:Without fin
DTM:Differential transformation method.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.