Study of the Inﬂuence of the Lack of Contact in Plate and Fin and Tube Heat Exchanger on Heat Transfer Efﬁciency under Periodic Flow Conditions

: Plate ﬁn-tube heat exchangers are widely used in air conditioning and refrigeration systems and other industry ﬁelds. Various errors made in the manufacturing process can result in the formation of an air gap between the tube and ﬁn. Several numerical simulations were carried out for a symmetric section of plate ﬁn-tube heat exchanger to study the inﬂuence of air gap on heat transfer under periodic ﬂow conditions. Different locations and sizes of an air gap spanning 1/2 circumference of the tube were considered for the range of airﬂow velocities. Velocity and temperature ﬁelds for cases with air gap were compared with ideal thermal contact cases. Blocking of heat ﬂow by the gap leads to the reduction of heat transfer rate. Fin discontinuity in the front of the tube causes the smallest reduction of the heat transfer rate in comparison to the ideal tube-ﬁn contact, especially for thin slits. The rear gap position is the worst in the smallest gap range. Therefore, reversing the ﬂow direction can lead to up to a 15% heat transfer increase, if mainly the rear gaps are present. The introduction of a thin slit in the front of the tube leads to convective heat transfer enhancement, which should be further investigated.


Introduction
Plate-fin and tube heat exchangers are most widely used as evaporators and condensers in refrigeration systems (commercial refrigeration) and in other applications where heat transfer between liquid or evaporating/condensing fluid and gaseous working fluid is required. Proven technology, reliability, and relatively low cost of manufacture make platefin and tube heat exchangers very popular, not only as components of refrigeration units but also as gas coolers/heaters, recuperative heaters, dryers or heat pipe, thermosyphon radiators [1][2][3]. Nowadays, when detailed numerical simulations of fluid flow with heat transfer are within reach, many of the researchers carry out Computational Fluid Dynamics (CFD) studies to optimize the geometry of plate-fin and tube heat exchangers [4,5]. The majority of these works assume an ideal tube fin contact that in practice has to be created during the manufacturing process. The contact is usually achieved by a tube expansion (plastic deformation of the tube by moving the expansion die or by applying pressure at its inner surface). These mechanical methods of expansion cause additional thermal resistance at the tube-fin interface, due to the gaps between the contact surfaces of the fin collars and tubes [4]. Formation of the gap can be also a result of the nonideal joint formation process or assembly requirements. For the egg-crate-type heat exchangers, which are used as evaporators in domestic refrigerators, slits between the tubes are introduced to allow one to pull the whole tube coil through the plate-fin packet [6]. In this case, the lack of contact is always present at approximately 1/2 perimeter of the tube, and it is a feature essential for the assembly of the heat exchanger.
The gaps can be also created because of the wear of the heat exchanger, e.g., as the result of corrosion by exposure to difficult weather conditions or by strong chemical interactions. In Figure 1, one can see a picture (courtesy of Blygold company [7]) of a corroded heat exchanger surface with gaps on the circumference of the tube. In this case, fin discontinuities are just one of the factors affecting heat transfer; there are also white corrosion products around the tubes that impede heat transfer (AlOH thermal conductivity is smaller than for Al).  [7]).
The quality of the joint can be improved by brazing the tube to the fin (the soldering process and its impact on heat transfer are discussed in detail in [8][9][10]). Because of the additional amount of energy needed for the soldering process and some technical difficulties (ensuring high purity of soldered surfaces), mechanical methods of expansion are most frequently chosen. Non-ideal tube fin connection can cause significant thermal contact resistance which deteriorates heat transfer. Therefore, some of the researchers have undertaken studies taking into account the realistic fin-tube thermal contact and its impact on the thermal efficiency of heat exchangers. Aside from the experimental research aimed at estimating the thermal contact resistance [11,12], its value is also obtained through semiempirical methods, partly based on numerical analysis. Taler and Oclon [13] and Taler and Cebula [14] changed iteratively the value of contact resistance in their CFD model in such a way that the thermal and flow quantities obtained from the numerical solution agreed with the measured ones. Less attention has been paid to cases where there is a lack of contact with the plate-fin along a part of the tube perimeter (no contact between the fin collar and tube). Recently, the topic was addressed by Singh et al. [15] who proposed a CFD model of a finned heat exchanger with a variable contact surface area between the tube and the fin. They considered several cases of gaps of different sizes. They found that the presence of gaps influenced heat transfer and flow (pressure drop) in a significant way. In contrast to [15], where several triangular gaps are distributed evenly along the joint circumference, our work presents a CFD analysis of plate-fin and tube heat exchanger, where there is a lack of contact between the tube and the fin at the half of the tube's circumference. The present study aims to numerically investigate the effect of the lack of contact on the energy efficiency of the heat exchanger for variable parameters, such as different air velocities and various positions of the gap relative to the flow direction. The CFD model is developed with the assumption of periodic flow boundary conditions. The novelty of the present work lies mostly in the consideration of the different locations and sizes of the gaps between the tube and fin relative to the flow direction. Because of the nonuniform heat transfer coefficients on the fin and tube surfaces, the expected degree of heat transfer deterioration is a function of the position of the cavity. Only in [15] was a similar CFD analysis done, but evenly distributed triangular voids excluded the possibility of dependency on the flow direction; additionally, the authors of [15] considered a slightly different geometry of the plate-fin tube heat exchanger than the typical one (with slitted fins). The present work assumes more realistic positions of gaps resulting from wear (corrosion) or erroneous tube expansion in the classic type of plate-fin tube heat exchanger. It contributes to the current area of research by expanding knowledge about the magnitude of heat transfer deterioration (fin efficiency) caused by the lack of contact between fin and tube depending on flow conditions. It can help to understand and predict the thermal effects of wear and manufacture errors on the functioning of compact heat exchangers.

Numerical Model of a Heat Exchanger
The heat exchanger under consideration works as a condenser in a refrigeration system. It is a plate-fin and tube heat exchanger, in which tubes are arranged in-line. It was recognized that the location of the gap relative to the flow direction is an important parameter because of the variation in the local heat transfer coefficients on the surface of the fin and the locally variable velocity field [16]. The gap location along tube-fin joint circumference should therefore have a significant impact on heat exchange and flow patterns. Figure 2 shows the analyzed configurations. Case I (without a gap-thickness of discontinuity δ = 0 mm) is treated as a reference. The slit position is nondimensionalized as θ/180 • , where θ is an angular position. In the second case, the gap spans the front side of the circumference (θ = 0 • ). In Case III, it is placed on the upper half perimeter (θ = 90 • ), in the last IV variant-in the back of the tube (θ = 180 • ). The angles indicate the centers of the gaps, taking into account the slit circumference span the corresponding angle ranges are front gap (90 and back gap (−90 • ≤ θ ≤ 90 • ). Four gap thicknesses were taken into account: 0.25, 0.5, 1.0, and 2 mm, corresponding to the following dimensionless values (relative to tube diameter): 0.03125, 0.0625, 0.125, and 0.25. Numerical simulations were carried out for airflow through the heat exchanger equal to 1, 2, 3, 4, and 5 m/s average velocities at the flow area directly upstream of the entrance to the device under consideration.

Geometry of the Computational Domain
The computational domain is a repeatable slice of the heat exchanger geometry, which is shown in Figure 3. The length and width of the domain are equal to the distance between axes of the tubes L f . The thickness of the fin in the computational domain is half of the actual fin thickness δ f , as the temperature field is symmetric at the opposite half of the fin (see the symmetry plane in Figure 3). Tube length in the domain equals 1 2 of the distance between fins. Taking the half of inter-fin spacing is a consequence of the flow symmetry assumption. The actual domain boundary is marked in Figure 3 by a dashed line and spans up to the symmetry plane. The specific dimensions of the domain are given in Table 1.

Boundary Conditions and Governing Equations
The shear stress transport (SST) turbulence model was chosen to solve Reynolds (time) averaged Navier-Stokes equations (RANS). SST combines the best features of k-ε and k-ω models. In the near-wall region, approximately up to half of the boundary layer k-ω formulation is used, which does not require any damping functions, while for the outer zone of the boundary layer k-ε is the governing turbulence model-it eliminates the problem of k-ω sensitivity to values of ω in the freestream (outside the boundary layer). The other advantage of SST is the accurate prediction of the boundary layer separation by application of the turbulence production limiter in stagnation regions [17]. The SST model is also characterized by very good accuracy of the numerical solutions for wall-bounded, complicated geometry flows, obtained with relatively low computational power. Recently, Lindqvist et al. have performed an extensive analysis of RANS turbulence models in the application to numerical simulation of airflow and heat transfer in plate-fin and tube heat exchangers [18]. The k-ω SST was chosen as optimal and used for further validation against experimental data. It was utilized in many works where CFD analysis of the heat transfer for finned and enhanced surfaces was carried out [4,13,[19][20][21].
The following simplifying assumptions were made in the present work: • Steady-state fluid flow and heat transfer; • Fluid flow and heat transfer are periodic (fully-developed), meaning that pattern of flow/thermal solution has a periodically repeating nature (this condition is generally fulfilled for tube rows greater than the fourth row); • Thermophysical properties of air are temperature-dependent (ideal gas); • Natural convection was not considered as the highest Richardson number calculated for simulation conditions is Ri = (g·β·∆T·L)/u 2 = 0.025 (for Ri < 0.1 mechanism of natural convection can be typically considered negligible).
Numerical calculations were carried out using the commercial ANSYS CFX code capable of solving the Navier-Stokes equations [22]: where τ is a stress tensor in the fluid, and ρu i u j are Reynolds stresses (turbulent), which depend on average velocity gradients (Boussinesq hypothesis): where k-turbulence kinetic energy and µ t -turbulent viscosity. The SST turbulence model allows for obtaining turbulent viscosity by using two transport equations of kinetic turbulence energy k and specific dissipation rate ω. Due to the extensive definition of the model, it will not be presented here; its formulation can be found in [17]. The turbulent viscosity is calculated using the equation: The energy equation for fluid is given by: where e-total fluid energy, T-temperature, λ-thermal conductivity, c p -specific heat at constant pressure, and Pr t -turbulent Prandtl number for energy (assumed value: 0.9). In the solid body (fin and tube domain), the heat is exchanged by means of steady conduction: The computational domain is divided into three parts: air (fluid), fin, and tube (solids). It was assumed that the thermal conductivity of the solid (aluminum) is constant in the tested temperature range, which is λ = 202.4 W/m/K. A periodic flow was solved with the following boundary conditions and parameters: On the other surfaces of the computational domain, a symmetry boundary condition was assigned, which can be expressed mathematically as: where n is a normal vector to the surface on which the boundary condition was set. Additionally, on the solid-fluid contact surfaces, a coupled boundary condition was present, ensuring the energy balance is satisfied between domains. The periodic condition of the velocity in the x-direction can be written as: where r-position vector and L-the periodic length vector (equations in the rest of the Cartesian directions are analogous). For viscous flows, the pressure field is considered periodic if the following condition is met: The pressure gradient is divided into two parts: a periodic gradient and a linear gradient: where p(r) is so-called periodic pressure and β L |L| -linearly variable pressure component.
Periodic pressure is the remainder after subtracting the linear pressure. On the other hand, the linearly variable component results in a force acting on the fluid in the momentum equations. During the solution process, the value of the coefficient β is iteratively changed until the desired convergence between the set and the calculated mass flow is obtained. The temperature field is analogously divided into two components for periodic simulation: where: Q . m·c p ·L (12) and .
Q-net heat transfer rate in a periodic domain. At the inlet and outlet, the periodic boundary condition has the form: The above-mentioned conditions for the periodicity of the airflow are practically implemented in the Ansys-CFX CFD software. The boundaries named in Figure 4 as periodic inlet and periodic outlet are coupled as two flow interfaces; the target mass flow can be specified, which is obtained by iteratively changing the pressure gradient β. Negative heat flux has to be fixed at the periodic outlet boundary to balance the heat transfer rate gained by air flowing through the domain, as the outflowing fluid is "returned" to the inlet (Figure 4). If a heat sink is not introduced, the air temperature rises infinitely during the solution process of the numerical simulation, causing the solver to crash after reaching some maximum temperature level. In the case of numerical simulations carried out within this research, the domain was properly thermally balanced to keep the inlet temperature at an approximately constant level of 22 • C. The symmetry boundary condition is applied at domain boundaries that were not marked in Figure 4 (for both solid and fluid surfaces).

Computational Grid
Numerical simulations of the fluid flow with coupled heat exchange process were carried out using ANSYS CFX CFD software. The computational grid for the case I is shown in Figure 5. It was created for the repetitive slice of heat exchanger with the use of ANSYS-Meshing application. The grid independence study was performed using a methodology similar to used in works [23,24]. The five grids were tested, namely, coarse, average, fine, very fine, and the finest corresponding to mesh sizes: 0.2 × 10 6 , 1.8 × 10 6 , 2.14 × 10 6 , 2.48 × 10 6, and 2.80 × 10 6 nodes. For all of the five meshes, the prism layer elements number at the wall was gradually increased, starting from 10 for the coarse grid by 25, 30, 35, up to 40 for the finest one. The Y + parameter and percent change of mean heat transfer coefficient h % were presented in Figure 6. The h % is the relative difference between the obtained average heat transfer coefficient for the ith mesh refinement level and i − 1 refinement (the one-step coarser mesh). The change does not exceed 0.5%; therefore, the average mesh (1.8×10 6 nodes) was chosen as optimal because it also satisfies the maximum Y + < 1 condition (maximum Y + = 0.44, average Y + = 0.24).

Validation of the Numerical Model
The computational model was validated against data published in the work of Saboya and Sparrow [16] as their paper is widely acknowledged as the standard for validation of the studies of heat transfer from finned surfaces. They used the naphthalene sublimation technique [25] to obtain local mass transfer coefficients (Sherwood numbers Sh) on the fin surface of the one-row plate-fin tube heat exchanger. The Sherwood number distribution can be converted to heat transfer coefficients by utilization of the heat-mass transfer analogy: where Sc-Schmidt number, d h -hydraulic diameter according to Kays and London [26], and n-heat-mass transfer analogy exponent. Rosman et al. [27] have chosen the value of the n = 0.4 based on his experiment and measurements available in the literature. However, the latest careful measurements have reaffirmed that the most suitable value of the exponent for both laminar and turbulent boundary layers is approximately n = 0.33 [23]. Prandtl number and thermal conductivity are assumed as for pure air (Pr = 0.7, λ = 0.026 W/m/K)the concentration of naphthalene is too low to significantly influence the thermophysical properties of the fluid. Saboya and Sparrow [16] assumed Sc = 2.5, yet it is an outdated value, and for room temperature the Sc = 2.28 is now recommended [25]. Local Sherwood numbers [16] were converted to heat transfer coefficients according to (14) and prepared to validate the CFD simulation. On the left side of Figure 7, the distribution of heat transfer coefficients on the fin surface is shown obtained from [16] for Re = 214. The distribution of local h is only known along the horizontal black lines where the mass transfer coefficients are measured-a 2D contour map is a result of an inverse distance interpolation made in Tecplot 360 EX software [28]. On the right side, there are results from numerical simulation for analogous flow and thermal conditions as in the experimental work. This setup allows a direct comparison between two contour maps. As can be seen, the ranges of heat transfer coefficient values are nearly the same. Near the front edge of the fin, the convective heat transfer seems to be more intensive in the case of the numerical results, but the lower h are a consequence of the interpolation of the Saboya and Sparrow's [16] data from the location significantly upstream of the edge of the fin. The accordance of the simulation's local h with the experiment is shown in Figure 8. Four lines were chosen, where the measurements were taken exactly-at the front (line 1), back (line 4), and middle segments (lines 2, 3) of the fin. Most of the experimental data is in very good agreement with the CFD solution (maximally ±20% and better). The greatest discrepancies are noted in the near tube region, but they can be caused by the fact that in the experiment [16] the tube surface was adiabatic. The uncertainty of the experimental h was also plotted in the form of error bars. It was based on the estimation of the uncertainty of Sherwood number measurements done in [29] (uncertainty within 7.3%). The validation presented above was made to test the ability of the model to predict local coupled heat and fluid flow phenomena. Because the experimental flow arrangement (one-row tube-fin) is not adequate to the simulation settings (periodical flow-more than four-row inline tube-fin geometry), additional validation was performed. Mean heat transfer obtained from the simulation was compared against VDI correlation [30], which is true for three rows inline heat exchangers: A comparison between simulation heat transfer coefficients and VDI correlation is presented in Figure 9. The two data sets' values are very close. The relative difference varies from 4.9% for the lowest Re to 1.8% for the biggest airflow.

Results and Discussion
The results of the numerical simulations were put together in a form of a series of graphs, allowing for the evaluation of the thermal efficiency of the cases with the presence of the gap to perfect fin-tube contact example. In Figure 10, the ratios of heat transfer rates (discontinuity to no-slit heat transfer rate) are presented in the function of the dimensionless gap thickness (gap thickness to diameter). For each Reynolds number, a curve approximating the numerical data is presented. For the gap in the back and front, the minimal ratio of heat transfer rate is approx. 0.75 for the considered airflow span, whereas the side gap attains a slightly lower minimum of 0.72. For the front and side gap localization, there is a bigger spacing between the curves (stronger influence of the airflow velocity) than for the rear position. This effect is the strongest for smaller gaps. For δ/D > 0.1, there is no significant difference between curves spacing for all examined cases. The shapes of the curves are similar for all the arrangements, although for the fin discontinuity in the rear of the tube the heat transfer ratio drop from unity (δ/D = 0) is much steeper than for the other locations. Heat transfer rate ratios are significantly lower for back orientation, especially for small gap thickness and small Re. Usually the higher the Reynolds number the lower value of . Q/ . Q 0 , although there is an exception for the lowest considered Reynolds number curve for nearly all gaps, excluding the biggest fin discontinuity. Q 0 values are below these for Re = 2013. One can see a nearly monotonic decrease of the heat transfer rate ratio with increasing gap thickness, being the sharpest in the lower δ/D range and reaching practically a plateau for δ/D > 0.1. The data points are connected by a piecewise linear function, instead of polynomial fit, because higher-order curve fitting produced slight inflection in the plateau range, which was a misleading artifact of the method.
A clearer picture of how the flow velocity influences heat transfer emerges from Figure 11. A decreasing trend of . Q/ . Q 0 is seen with increasing Reynolds number, excluding the smallest Re departing from this pattern. In the range of Re > 1342, there is a nearly linear negative slope. However, in the span of Re = 671-1342, the decrease of heat transfer rate ratio is smaller or even there is a reverse trend (increase), indicating there can be a change of the flow regime. The curve for the smallest gap is placed remarkably higher than the other ones for the side and front location. Other δ/D cases, for all slit positions, exhibit similar trends, with the δ/D = 0.125 curve approaching or even crossing δ/D = 0.25 with decreasing Re. Q 0 is addressed in Figure 12. In the case of the smallest gap, the front and side locations exchange nearly equal heat transfer rates, with a significant decrease that can be seen for slit in the back (θ/180 • = 1). For the rest of the considered gap thicknesses, there are very small differences between the tube circumference and gap locations. The weak dependence in the range of δ/D = 0.0625-0.25 is the most evident for the highest Re = 3355. The back and front slit cases indicate moderately higher, nearly equal thermal throughputs relatively similar to the side example. In Figure 13, the ratio of mean heat transfer coefficient with the gap presence to heat transfer ratio without any fin discontinuity is presented in the δ/D function. The front and side location results in enhancement of the mean convective heat transfer for the smallest discontinuity. This effect is not present in the rear location.   The dependence on the gap position is shown in Figure 15. At δ/D = 0.03125, there is the highest difference between the front, side, and back locations. The remaining slit thicknesses exhibit a very weak dependence on θ/180 • , typically slightly reducing, with a maximum for the front and minimum for the back case. In Figure 16, the contour map of a local heat transfer coefficient, which is defined as: where T in is periodic inlet air temperature, is presented for different gap sizes at the lowest Reynolds number (Re = 671) for front fin discontinuity location. Figure 16 was arranged to investigate the convective heat transfer enhancement relatively to no gap case. Intensification is present for the thinnest slit (δ/D = 0.03125)-for bigger gaps, this effect diminishes ( Figure 10). The contours are similar for all cases shown; therefore, the analysis starts with a description of common map characteristics. There is a high h x zone at the inlet edge of the fin-this is the region where the thermal boundary layer is the thinnest. The area of this region is reduced at the center because of the presence of the fluid circulation flow spanning from the leading edge of the fin to the surface of the tube ( Figure 17). This circulation is caused by the trace flow downstream of the upstream tubebeing a result of boundary layer detachment. In Figure 17, one can see the circulation regions visualized by streamlines. The circulation area is wider in the rear region of the fin (behind the tube) than upstream of the tube (air flows from the bottom to the top in Figure 17). Because the circulating fluid moves with reduced velocity, local heat transfer coefficients are significantly lower. On the sides of the tube, air flows with increased velocity (Figure 17)-this corresponds to the high h x region in Figure 16 (contour level: 4, h x = 18-24 W/m 2 /K).  At the fin surface zone near the front part of the tube circumference, there are two bands of intensified heat transfer (placed symmetrically on the sides). This local enhancement is caused by a horseshoe vortex, which is visible in Figure 17, where the velocity vector field is visualized at plane A. The swirl of the fluid can be observed at the base of the tube. For the one-row or staggered arrangement plate-fin tube heat exchangers, the horseshoe vortex spans the whole half of the tube's circumference. In the present analysis, periodic flow simulations are equivalent to the inline arrangement. Therefore, the wake behind the preceding tube divides the horseshoe vortex in two (Figure 16-δ/D = 0). There is also a thin region of enhanced heat transfer coefficient at the trailing edge, widening at the middle of the fin. This is the effect of the implementation of periodic heat transfer conditions in the simulations. Because the negative heat flux is set on the outlet, the small backflow of colder fluid causes this "artificial" heat transfer intensification. This effect is practically of the same magnitude for all cases-it cancels in the comparative analysis. Even the h x contour maps share the same characteristics; for δ/D = 0, the strength of the heat transfer enhancement caused by the horseshoe vortex is the greatest. With the increase of the gap thickness, the horseshoe intensification becomes weaker. The thinnest slit covers partially the area of vortex action, whereas in the δ/D = 0.125 case the vortex is present completely over the gap, which results in only moderately increased local h on the edges of the fin discontinuity. The introduction of the gap improves the heat transfer on the leading edge of the fin and changes somewhat the flow and heat transfer characteristics in the region downstream of the tube.
For δ/D = 0.03125 and δ/D = 0.0625, this alteration is not significant, yet for δ/D = 0.125 there is a dramatic change-low heat transfer coefficient region (contour 3 in Figure 16) doubles its size and merges with the smaller one at the upstream part of the tube. The combined effect of the enhancement on the leading edge with still strong horseshoe vortex intensification can explain the observed augmentation of mean HTC for the thinnest slit relatively to no-gap example ( Figure 13). The mean of heat transfer enhancement with thin slits cut in the front of the tube (for low airflow velocities) should be further studied. It could be a promising method of heat transfer intensification in compact heat exchangers.
The enhancement of h x , by the presence of the slit, relative to ideal fin-tube contact example is not enough to cancel the negative effects of the fin discontinuity at its base (tube-fin junction)-in Figure 11 one can observe at least 6% lower heat transfer rate in the respect to δ/D = 0. In Figure 18, the most obvious consequence of the introduction of the gap can be seen-the blockage of heat flow in the fin area adjacent to the discontinuity. This effect is visualized by dimensionless temperature excess Θ contour map, where: The fin is the most isothermal for the no-gap example, and the contour map is the most symmetrical. In Figure 18, the heat flux blockage is evident for the cases with the gap-Θ is significantly lower for nearly half of the fin on the side of the slit. For the gap in the front, the reduction of temperature is the smallest, while for the side location it is the highest (lowest fin efficiency). The mean h on the sides of the fin is the highest, and that is the reason for the severe temperature decrease. The Θ is the biggest on the average for the front case because the heat flux is directed backward (low h x region) and the enhancing effect of the horseshoe vortex in the front of the tube is partially blocked. Analyzing the back gap position, the Θ is significantly lower on average. Heat flux is directed towards the front where the heat transfer intensification regions cause effective cooling of the fin. The result is substantially lower fin temperature in the back. In a conclusion, the gap placement in the front blocks the heat flow in the direction of the enhanced convective heat transfer area, which is not beneficial, although it results in higher fin efficiency than for the gap in the back. This can explain why there is a sudden heat transfer rate ratio drop for the rear slit position. The effect of low fin efficiency seems to combine with poorer convective heat transfer relatively similar to other discontinuity placement cases ( Figure 10). The side gap position is the worst of the considered positions in terms of the fin temperature contour. The colder fin area is slightly bigger than the gap in the back. As the heat flux is directed to the side where high-velocity flow increases the intensity of heat transfer, the fin is efficiently cooled there. As the mean heat transfer coefficient is nearly as high as for the gap in the front, . Q/ . Q 0 is slightly lower than the front gap for thin slits. The enhanced convective heat transfer compensates for the low fin efficiency. Figure 19 shows the influence of the flow Reynolds number on the h x . The h x contour maps are presented at four Reynolds numbers for a 1 mm thick gap in the back. The heat transfer intensification caused by the presence of the horseshoe vortex is less significant for the lowest Re, and it gets increasingly distinctive as the airflow velocity increases. The thermal boundary layer becomes thinner at the leading edge of the fin with the raising Re, which is visible as the wider high h x area. As the Re increases, the h x contours change their values proportionally, but the overall topology (shape of the contours) is analogous for all airflow velocities. In Figure 14, one can see that for Re = 671 h/h 0 is the lowest-this is mainly caused by the low activity of the horseshoe vortex ( Figure 19).

Conclusions
To sum up the key findings from the numerical simulation, the fin discontinuity along the circumference reduces the heat transfer rate relatively to the perfect fin-tube contact case. The placement of the gap reduces the skin friction locally to zero, which can promote beneficial flow structures, which causes enhancement concerning the no-gap case, especially for small slit sizes. However, the convective heat transfer enhancement is overpowered by the heat flux blockage effect which leads to . Q/ . Q 0 < 1 for all the cases. Fin discontinuity in front of the tube causes the smallest reduction of the heat transfer rate in comparison to the ideal tube-fin contact, especially for thin slits. Even the heat transfer rate ratio minimum is obtained for the gap at the side position (approx. 0.72), it performs slightly worse for the small gap range than the front location. The rear gap is characterized by the mean . Q/ . Q 0 = 0.8 for all the gap sizes and therefore is considered the worst scenario, without the advantage of the smallest gap range. For gaps wider than δ/D > 0.1, the heat transfer rate reduction becomes nearly constant for a specified Reynolds number. Generally, the higher the airflow, the . Q/ . Q 0 becomes lower. For multirow plate-fin tube heat exchangers with small fin discontinuities, reversing the flow direction can result in an even 15% heat transfer increase if most of the gaps are in the rear of the tube (low airflow velocity). Realistically, obtaining this maximum is not very probable, although the present study shows that reversing the flow can have a positive effect on the plate-fin tube heat exchanger performance, and its longevity depends on the nature of the origin of the discontinuities (are they, and how fast they enlarge). The convective heat transfer enhancement for thin slits in the front should be further investigated to test the possibility of the introduction of fin cut-outs as heat transfer intensification structures in compact heat exchangers.

Conflicts of Interest:
The authors declare no conflict of interest.