E. BARTULI, M. RAUDENSKY: NUMERICAL INVESTIGATION OF HEAT TRANSFER ON THE OUTER SURFACE ... 459–463 NUMERICAL INVESTIGATION OF HEAT TRANSFER ON THE OUTER SURFACE OF POLYMERIC HOLLOW FIBERS NUMERI^NA RAZISKAVA PRENOSA TOPLOTE NA ZUNANJI POVR[INI VOTLIH POLIMERNIH VLAKEN Erik Bartuli, Miroslav Raudensky Brno University of Technology, Faculty of Mechanical Engineering, Heat Transfer and Fluid Flow Laboratory, Technicka 2896/2, 616 69, Brno, Czech Republic Erik.Bartuli1@vut.cz Prejem rokopisa – received: 2016-07-25; sprejem za objavo – accepted for publication: 2018-02-02 doi:10.17222/mit.2016.221 Polymeric hollow-fiber heat exchangers (PHFHEs) are heat exchangers that use polymeric microchannels as the heat-transfer surface. The main purpose of this work was to determine the optimal arrangement of polymeric hollow fibers in a PHFHE. The object of study was a bundle of parallel polymeric hollow fibers with an outer diameter of 0.8 mm. The heated working fluid flew inside the hollow fibers. The heat was transferred by a cooling medium from the outside of the hollow fibers. Heat-transfer coefficients on the outer surface of the hollow fibers were obtained using computer modelling. For the case of water cooling, a sharp heat-transfer-coefficient increase was observed as the angle between the hollow-fiber position and the direction of the cooling-medium flow increased to 30°; however, its further increase to 45° had no significant effect. With the water cooling, the total heat-transfer coefficient of the PHFHE increased by 35 % and with the air cooling, it increased by 15 %. The angle of 30° can be considered as the optimal angle for manufacturing such PHFHEs. Keywords: heat transfer, heat exchanger, polymeric hollow fibers, numerical modelling Toplotni izmenjevalniki z votlimi polimernimi vlakni (PHFHE; angl.: Polymeric hollow fiber heat exchangers) so izmenje- valniki toplote, ki uporabljajo polimerne mikrokanale kot povr{ine za prenos toplote. Glavni namen pri~ujo~ega dela je bil dolo~iti optimalno razporeditev votlih polimernih vlaken v PHFHE. Predmet {tudije je bil sve`enj vzporednih votlih polimernih vlaken z zunanjim premerom 0,8 mm. Ogrevano delovno sredstvo (voda ali zrak) se je pretakalo znotraj votlih vlaken. Toplota se je prena{ala z ohlajanjem s hladilnim medijem, ki se je nahajal na zunanji strani vlaken. Koeficiente prenosa toplote na zunanji povr{ini votlih vlaken so avtorji tega ~lanka dobili z ra~unalni{kim modeliranjem. V primeru vodnega hlajenja so zaznali mo~no nara{~anje koeficienta prenosa toplote, ko je kot med polo`ajem votlih vlaken in smerjo pretoka hladilnega medija narastel na 30°, nadaljnje nara{~anje kota do 45° pa ni imelo pomembnega vpliva. Z vodnim hlajenjem je celotni koefi- cient prenosa toplote PHFHE narastel za 35 %, pri zra~nem hlajenju pa za 15 %. Avtorji so kot optimalen za proizvodnjo tak{nih PHFHE ocenili kot 30°. Klju~ne besede: prenos toplote, toplotni izmenjevalnik, votla polimerna vlakna, numeri~no modeliranje 1 INTRODUCTION Nowadays the use of polymeric heat exchangers is becoming more and more popular in different fields, such as heat-recovery system, evaporative cooling sys- tem, desiccant-cooling system, electronic-device cooling and water-desalination system. 1 D. Rousse et al. 2 pre- sented the theoretical and experimental analysis of using a plastic heat exchanger as a dehumidifier in a green- house for the agriculture industry. Due to the corrosive resistance of the polymer, heat exchangers can be used in the chemical industry where aggressive media are used with metals. L. Jia et al. 3 pre- sent an experimental study on the heat-transfer perfor- mance of a wet flue-gas heat-recovery system, using a plastic, longitudinal, spiral-plate heat exchanger. Another important advantage of a polymer heat ex- changer is that its production uses twice less energy than the production of metal heat exchangers. W. Liu et al. 4 tested two types of heat exchangers (tube-in-shell and immersed tube) made either from high-temperature nylon or cross-linked polyethylene for solar water-heat- ing systems. Authors reported that the polymer heat exchangers can provide thermal outputs equivalent to the conventional copper heat exchangers at a lower cost. The cost of a nylon tube-in-shell heat exchanger is about 80 % of that of a copper tube-in-shell heat exchanger. Other significant advantages of plastic heat ex- changers are low weight, smooth surface, simplicity of shaping and producing, resistance to fouling. 5 Most commercially available polymer heat exchangers can be used in low-temperature applications. 6 One type of polymer heat exchanger is the polymeric hollow-fiber heat exchanger (PHFHE). 7 Using thin-wall polymeric hollow fibers as the heat-exchanger tubes, it was first proposed by D. Zarkadas et al. 8 as a new type of heat exchanger for lower temperature/pressure applica- tions. This heat exchanger utilizes polymeric micro- channels as the heat-transfer surface. The outer diameter of these microchannels is smaller than 1 mm. The heat exchanger is made of hundreds of such fibers, resulting in a very large heat-transfer area compared to the size of Materiali in tehnologije / Materials and technology 52 (2018) 4, 459–463 459 UDK 67.017:677.494.7:536.24 ISSN 1580-2949 Original scientific article/Izvirni znanstveni ~lanek MTAEC9, 52(4)459(2018) the entire heat exchanger (Figure 1). This design allows the PHFHE to compete with metal heat exchangers in terms of heat-transfer intensity. I. Krasny et al. 9 con- ducted comparative tests of PHFHE and a conventional automobile aluminum finned-tube radiator and observed that the heat-transfer rates (up to 10.2 kW) and the overall heat-transfer coefficients (up to 335 W/(m 2 K)) of the PHFHE were competitive with the metal heat exchanger. The theory of the heat transfer inside a polymeric hollow fiber has been well-described. 10 However, the description of the heat-transfer processes on the outer surface of the fibers is more complex and currently still has not been sufficiently explored. This work was devoted to the numerical modelling of the heat transfer on the external surface of polymeric hollow fibers. The main purpose of this work was to determine the optimal configuration of polymeric hollow fibers in a PHFHE allowing it to achieve the maximum heat-transfer coefficients on the external surface of poly- meric hollow fibers and of the total device as a whole. The hollow-fiber position with respect to the direc- tion of the cooling-medium flow, determined by the angle , was examined as a variable parameter. Figure 2 shows a real PHFHE with the angle of about 20°. 2 EXPERIMENTAL PART The polymeric hollow-fiber heat exchanger studied in this paper was designed and manufactured in our Heat Transfer and Fluid Flow Laboratory (Figure 3a). It is based on the tube-in-tube type heat exchanger. Inside the outer tube is a polymeric hollow-fiber bundle made of polypropylene. Individual fibers were weaved with a thread that made it possible to put them in order (Figure 3b). Hot water is used as the working fluid inside the fibers. Outside of these fibers, heat is removed by a cold medium, which can be liquid or gaseous and should be compatible with polypropylene. In this work, pairs of working fluids of water-water and water-air were considered. The tested heat exchanger contained about 1800 hollow fibers. Due to the large number of fibers, the simulation of the whole internal volume of the heat exchanger took a lot of computing power. Therefore, a rectangular parallelepiped domain with a square cross- section of 10 mm × 10 mm and a length of 100 mm was selected as the object of investigation. The fibers in the model, which were considered to be hollow cylinders, passed through the cooling medium’s domain. The outer diameter of each fiber was 0.8 mm and the interaxial distances between adjacent fibers were 2.4 mm. As the fiber wall temperature at this length changes slightly in reality, in the model, it was set as a constant temperature of 60 °C. Air and water were chosen as the cooling media with the velocities equal to 0.3 m/s and 2 m/s, respectively. The pressure inside the domain was 1 bar. Plastic fibers have a low surface roughness (less than 1 μm); therefore, it was not considered in this model. E. BARTULI, M. RAUDENSKY: NUMERICAL INVESTIGATION OF HEAT TRANSFER ON THE OUTER SURFACE ... 460 Materiali in tehnologije / Materials and technology 52 (2018) 4, 459–463 Figure 1: Twisted polymeric hollow-fiber bundle Figure 3: Polymeric hollow-fiber heat exchanger Figure 2: Angle in a real PHFHE Angle in the model was set according to Figure 4. The largest value of the angle that can be set in a real heat exchanger is 45°. Therefore, in this simulation, this parameter varied from 0° to 45°. The heat-exchange process was simulated using the finite-volume method in ANSYS CFX. 11,12 Because the aim of the work was to compare the heat-transfer coeffi- cients for the hollow-fiber surfaces, the domain size and boundary conditions for each case were identical. The mesh was set using the ANSYS meshing prog- ram. For = 0°, the mesh had quadro elements. The number of these elements was equal to 3.7 mL. One main parameter for the mesh quality is the Y + coeffi- cient. It represents the non-dimensional wall distance for the wall-bounded flow and can be described with Equation (1): Y uy + = ⋅ (1) where u is the friction velocity at the wall, y is the dist- ance to the wall and is the local kinematic viscosity of the fluid. For a heat-transfer calculation, Y + near the wall must usually be less or equal to 1. When = 0°, Y + max is equal to 0.94. When > 0°, a mesh with tetra elements was used in the volume and the near-wall layer was defined by the quadro elements. For different values of , the amount of elements varied from 5 to 7 million and the Y + max never exceeded 0.98. General views of each mesh are shown in Figure 4. At angle = 0°, R e < 1000 and the water-cooling flow regime was laminar. A laminar model for modelling this flow was used for this case. The heat from the fiber outer surface was transferred to the cooling medium only by the thermal conductivity of the water. For > 0°, the SST-turbulence model was selected as a combination of the k- and k- models. The k- model was used for simulating the near-wall layer turbulence. The turbulence in the flow volume is described by the k- model. 13 The LES and DES tabulating models facilitate a simulation of the flow turbulence in more details, but they require significant computing-power resources and for the case under discussion, they are re- dundant. 14 For all the angle- values, the air-cooling flow regime was turbulent and the SST-turbulence model was used. 3 RESULTS 3.1 Total thermal resistance of the heat exchanger Equation (2) was used for the calculation relating to the heat-transfer coefficient and the thermal resistance was inversely proportional: R h = 1 (2) where h is the heat-transfer coefficient. Total thermal resistance of the heat exchanger can be presented as a sum of the resistances: R total = R in + R wall + R out (3) where R in is the thermal resistance of the heat transfer from the hot fluid inside the fiber to its wall, R wall is the thermal resistance of the fiber wall, R out is the thermal resistance of the heat transfer from the fiber wall to the cooling medium outside the fiber. Each fiber is a capillary tube and the heat exchange inside the tubes with a small diameter is very intensive (Figure 5), so R in is a relatively small part of R total and attempts for its further reduction do not have a signifi- cant impact on R total . R wall depends on the wall thickness. For this PHFHE, hollow fibers with a wall thickness of 0.1 mm were used. This value was determined on the basis of the strength properties of the fibers. A further reduction in the wall E. BARTULI, M. RAUDENSKY: NUMERICAL INVESTIGATION OF HEAT TRANSFER ON THE OUTER SURFACE ... Materiali in tehnologije / Materials and technology 52 (2018) 4, 459–463 461 Figure 4: Angle between hollow fibers and the cooling-medium flow Figure 5: Computational mesh for the cases of: a) = 0° and b) = 45° thickness would have adversely affected the reliability of the heat exchanger. So, for this type of PHFHE, value R wall could be considered as a constant. That is why the main attention in this paper is focused on the intensity of the heat transfer on the outer surface of the fibres (R out ). The contribution of R out to the total thermal resistance of the heat exchanger at =0°is about 35–40 %. 3.2 Water-cooling medium The computer-simulation results indicating the dependence of the heat-transfer coefficient on the outer- fiber surface on the angle between the cooling-flow direction and the fibers is shown on Figure 6. The graph shows that even an increase in angle by up to 5° results in a twofold increase in the heat-transfer coefficient. This fact is explained with the change of the laminar cooling-flow regime in the heat exchanger with the angle of = 0 to the turbulent cooling-flow regime for the exchangers with > 0. Hollow fibers in the water flow and the water start to intermix, leading to an intensifi- cation of the heat exchange on the outer surface of the hollow fibers. It can be seen from the graph that heat-transfer coefficients increase linearly when angle changes from 0° to 30°. Further increases in , from 30° to 45°, lead to a very slight increase in the heat-transfer coefficients. A decrease in the growth rate of the heat-transfer coeffi- cients can be explained with the regime of advanced turbulence. The maximum heat-transfer coefficient was 10,563 W/(m 2 K). Figure 7 illustrates the dependence of the total heat- transfer coefficient of the whole heat exchanger on angle . Even an increase in by 5° increases the heat-transfer coefficient by nearly 20 %. At the maximum angle of = 45°, the heat-transfer coefficient is increased by 35 % and is 1267 W/(m 2 K). 3.3 Air-cooling medium For the air-cooling medium, the main contribution to the total thermal resistance of the heat exchanger is the resistance on the outer wall of the fibers (R out ) due to poor air conductivity. Figure 9 shows that the heat-trans- fer coefficient on the outer surface of the fibers is almost equal to the total heat-transfer coefficient of the heat exchanger. At = 0°, the air flow is already turbulent, so E. BARTULI, M. RAUDENSKY: NUMERICAL INVESTIGATION OF HEAT TRANSFER ON THE OUTER SURFACE ... 462 Materiali in tehnologije / Materials and technology 52 (2018) 4, 459–463 Figure 8: Total heat-transfer coefficient of the hollow-fiber heat ex- changer vs angle Figure 6: Dependence of heat-transfer coefficients inside the tubes on the tube diameter Figure 9: Heat-transfer-coefficient dependence on angle at the air- cooling mediu Figure 7: Heat-transfer coefficient on the outer surface of hollow fibers vs angle an increase in angle does not greatly increase the heat-transfer coefficient. When angle changes from 0° to 45°, the total heat-transfer coefficient changes from 41.1 W/(m 2 K) to 47.3 W/(m 2 K). 4 CONCLUSIONS A computer simulation was used for determining the optimal arrangement of the polymeric hollow fibers in a PHFHE. Air and water were examined as cooling me- diums. It was found that with the water cooling, even an in- crease by 5° in the angle between the fiber axis and the cooling-medium flow direction leads to a twofold in- crease in the heat-transfer coefficient on the outer surface of the hollow fibers. A sharp heat-transfer-coefficient increase was ob- served as angle increased to 30°. Changing angle to 45° had no significant effect. The angle of 30° can be considered as optimal for manufacturing such PHFHEs. The maximum value of the heat-transfer coefficient on the outer surface of the fibers with water as the cool- ing medium was 10,563 W/(m 2 K). The maximum heat- transfer coefficient of the full PHFHE was 1267 W/(m 2 K). These values are comparable to the ones achieved in metal heat exchangers. For the air cooling, a linear increase in the heat-trans- fer coefficient was observed when angle was increased to 45°. The total heat-transfer coefficients of the PHFHE increased by 15 % when angle changed from 0° to 45°. However, the process of heat exchange with air is non-intensive, so the heat-transfer coefficients increased from 41.1 W/(m 2 K) to 47.3 W/(m 2 K). Acknowledgement The research leading to these results received funding from the MEYS under the National Sustainability Programme I (Project LO1202) and project EPSILON No.TH01020139, Heat exchangers with polymeric hollow fibres for energetic systems of buildings, granted by the Technology Agency of the Czech Republic. 5 REFERENCES 1 X. Chen, Y. Su, D. Reay, S. Riffat, Recent research developments in polymer heat exchangers, Renewable and Sustainable Energy Reviews, 60 (2016), 1367–1386, doi:10.1016/j.rser.2016.03.024 2 D. Rousse, D. Martin, R. Theriault, F. Leveillee, R. Boily, Heat reco- very in greenhouses: a practical solution, Applied Thermal Engineering, 80 (2000) 8, 687–706, doi:10.1016/S1359-4311(99) 00048-4 3 L. Jia, X. Peng, J. Sun, T. Chen, An experimental study on vapor condensation of wet flue gas in a plastic heat exchanger, Heat Trans- fer – Asian Research, 30 (2001) 7, 571–580, doi:10.1002/htj.10000 4 W. Liu, J. Davidson, S. Mantell, Thermal analysis of polymer heat exchangers for solar water heating: a case study, Journal of Solar Energy Engineering – ASME, 122 (2000) 2, 84–91, doi:10.1115/ 1.1288027 5 I. Astrouski, M. Ruadensky, M. Dohnal, Fouling of polymeric hollow fiber heat exchanger by wastewater, Chemical Engineering Transac- tions, 45 (2015), 949–954, doi:10.3303/CET1545159 6 T. Malik, C. Bullard, Suitability of polymer heat exchangers for air conditioning applications, College of Engineering, University of Illinois at Urbana-Champaign, 2005 7 I. Astrouski, M. Ruadensky, The study of polymeric hollow fiber heat exchangers, 18th International Conference Engineering Mecha- nics, Svratka, Czech Republic, 2012, 47–57 8 D. Zarkadas, B. Li, K. Sirkar, Polymeric hollow fiber heat ex- changers (PHFHEs): a new type of compact heat exchanger for lower temperature applications, Proceedings of the ASME Summer Heat Transfer Conference, (2005), 429–38 9 I. Krasny, I. Astrouski, M. Raudensky, Polymeric hollow fiber heat exchanger as an automotive radiator, Applied Thermal Engineering, 108 (2016), 798–803,doi:10.1016/j.applthermaleng.2016.07.181 10 H. Herwig, Flow and heat transfer in micro systems: is everything different or just smaller, ZAMM Journal of applied mathematics and mechanics, 82 (2002) 9, 1–10, doi:10.1002/1521-4001(200209) 82:93.0.CO;2-V 11 J. Ferziger, M. Peric, Computational methods for fluid dynamics, 2002, 423 12 R. Lewis, P. Nithiarasu, K. Seetharmu, Fundamentals of the finite element method for heat and fluid flow, 2004, 341 13 H. Versteeg, An introduction to computational fluid dynamics: the finite volume method, Longman, 1995, 257 14 T. Ganesan, M. Awang, Large eddy simulation (LES) for steady-state turbulent flow prediction, Engineering applications of computational fluid dynamics, (2015), 17–32, doi:10.1007/978-3-319-02836-1_2 E. BARTULI, M. RAUDENSKY: NUMERICAL INVESTIGATION OF HEAT TRANSFER ON THE OUTER SURFACE ... Materiali in tehnologije / Materials and technology 52 (2018) 4, 459–463 463