RETRACTED ARTICLE: Numerical assessment of the influence of helical baffle on the hydrothermal aspects of nanofluid turbulent forced convection inside a heat exchanger

This study is devoted to the numerical assessment of the influence of helical baffle on the hydrothermal aspects and irreversibility behavior of the turbulent forced convection flow of water-CuO nanofluid (NF) inside a hairpin heat exchanger with 100 mm length, 10 mm inner tube internal diameter, and 15 mm outer diameter internal diameter. The variations of the first-law and second-law performance metrics are investigated in terms of Reynolds number (Re = 5000–10,000), volume concentration of NF (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi =0-4\mathrm{\%}$$\end{document}φ=0-4%) and baffle pitch (B = 25–100 mm). The results show that the NF Nusselt number grows with the rise of both the Re and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document}φ whereas it declines with the rise of B. In addition, the outcomes depicted that the rise of both Re and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document}φ results in the rise of pressure drop, while it declines with the increase of B. Moreover, it was found that the best thermal performance of NF is equal to 1.067, which belongs to the case B = 33.3 mm, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document}φ=2%, and Re = 10,000. Furthermore, it was shown that irreversibilities due to fluid friction and heat transfer augment with the rise of Re while the rise of B results in the decrease of frictional irreversibilities. Finally, the outcomes revealed that with the rise of B, the heat transfer irreversibilities first intensify and then diminish.


Abbreviations
Heat flux (W/m 2 ) S gen Local entropy generation rate (W/m 3  It is well known that the turbulent flow has a higher heat exchange rate and pumping power than the laminar flow, the former being desirable and the latter undesirable 1-3 . So, this idea came into the researchers' minds to put equipment in the path of laminar flow and create local turbulence [4][5][6][7][8][9] . The idea was very successful and was widely used in the industry today [10][11][12] . This equipment is called turbulator, and so far, various types of turbulators were introduced and their performance was studied experimentally and numerically. Also, the use of nanostructures was developed to increase the efficiency of engineering systems [13][14][15][16] . Although the use of turbulators has led to an improvement in the performance of heat transfer systems, this has not prevented researchers from looking for ways to further improve the performance of these systems. One of these amazing techniques, which are originated from the low thermal conductivity of heat transfer fluids, is the use of nanofluids (NFs) instead of common heat transfer fluids [17][18][19] . Choi 20 first made these modern fluids and called them NFs. After the introduction of NFs and their amazing thermal properties, much research was done on their performance in diverse applications such as thermal management of electronic components 21,22 and photovoltaic panels 23,24 , nuclear applications 25,26 , and performance improvement of solar thermal collectors [27][28][29] , heat exchangers 30-32 and car radiators [33][34][35] . The literature inspection shows that the performance of thermal systems with NF coolants equipped with turbulators was investigated by various researchers. Bellos et al. 36 analyzed the efficacy of oil-CuO NF in a parabolic trough collector equipped with turbulators. They found that using the combination of NF and turbulator causes a 1.54% thermal efficiency improvement. Nakhchi and Esfahani 37 inspected the efficacy of aqueous Cu NF inside a heated tube equipped with the perforated conical rings in a turbulent regime. It was reported that using compound NF and turbulator results in considerable heat transfer intensification. Akyurek et al. 38 experimentally evaluated the forced convection flow of water-Al 2 O 3 NF inside a horizontal tube equipped with a wire coil turbulator. They utilized two turbulators with different pitches and found that the performance metrics of a tube filled with NF without any turbulator are superior to that of the cases with turbulator. Xiong et al. 39 simulated the efficacy of aqueous CuO NF flowing inside a tube having a compound turbulator. The outcomes portrayed that using turbulator elevates the rate of heat exchange between tube wall and NF. Xiong et al. 40 numerically investigated the forced convection of NF through a pipe equipped with a complex-shaped turbulator. They inspected the consequence of elevating width ratio, flow rate, and pitch ratio on the performance features. It was found that the Nusselt number intensifies with boosting pitch ratio of turbulator. In a numerical investigation, Ahmed et al. 41 explored the forced convection of water-Al 2 O 3 and water-CuO NFs inside a triangular duct with a delta-winglet pair of turbulator under a turbulent flow regime. They reported the significant effect of using both NF and turbulator on the performance aspects. The design of a heat transfer system can be done both according to the first law and the second law of thermodynamics. If the first law applies, the system must have the highest overall hydrothermal performance, and if the second law applies, the system performance must have the least irreversibility. The efficacy of NF flow in turbulator-equipped thermal units has rarely been investigated from a second-law perspective [42][43][44][45] . Sheikholeslami et al. 42 inspected the irreversibility features of turbulent flow of aqueous CuO NF inside a pipe equipped with complex turbulators. Li et al. 43 simulated the NF irreversibility in a tube with helically twisted tapes. The thermal irreversibility was found to be declined with elevating the height ratio of turbulator, while the opposite was true for the frictional irreversibility. Farshad and Sheikholeslami 44 analyzed the irreversibility aspects for aqueous Al 2 O 3 NF flow in a solar collector having a twisted tape. The outcomes revealed that the irreversibility diminishes with the increase of diameter ratio. Al-Rashed et al. 45 examined the influence of nanomaterial type on the irreversibility production of an NF in a heat exchanger. It was reported that the maximum irreversibility belongs to the platelet shape nanomaterials.
This numerical work aims to evaluate the features of turbulent flow of water-CuO NF through a hairpin heat exchanger equipped with helical baffles in the annulus side from both the first and second-law perspectives. The www.nature.com/scientificreports/ impacts of B , Re , and ϕ on the NF efficacy are assessed. This investigation is the first work on the consequences of using a helical baffle on the irreversibility production inside the annulus of a hairpin heat exchanger filled with NF. Figure 1 gives a demonstrative sketch of the geometry under investigation. It is a hairpin heat exchanger with 100 mm length, 10 mm inner tube internal diameter, and 15 mm outer diameter internal diameter. Additionally, the wall thickness for both the inner and outer tubes is 1 mm. Moreover, the considered B values are 25 mm, 33 mm, 50 mm, and 100 mm. The purpose of using this device is to cool the water-CuO NF passing through the annulus with the help of water passing through the inner tube. In both water and NF streams, time is of the essence and both streams are in a turbulent regime. Both streams enter the device at a uniform velocity and temperature and are discharged into the atmosphere and as a result, the relative pressure at the outlets is zero. Also, the exterior wall of the heat exchanger is considered insulated and the no-slip condition is utilized on the walls.

Governing equations and problem parameters
To simulate the steady, incompressible, and turbulent flow of nanofluid and water in the annulus and inner tube of the hairpin heat exchanger studied in the present study, the conservation equations of mass, momentum, and energy must be solved. These equations are as follows 46-48 : • Continuity equation where µ t and ρ nf u ′ j u ′ i are defined as follows; The k-ε turbulence scheme is employed to establish the turbulence. To this end, two additional equations for turbulent kinetic energy ( k ) and the turbulent dissipation rate ( ε ) are added which are as follows 42 : The irreversible generation in forced convection flow of NF flow is due to two sources, namely fluid friction and heat exchange. Therefore, the total irreversibility is computed as 47 : By integrating the whole computational domain, the total irreversibility can be obtained as 47 :  www.nature.com/scientificreports/ The characteristic of an NF includes its thermophysical properties that establish the relationship between the base fluid and the nanoparticles. These basic relationships for the water-CuO NF are defined as follows [50][51][52][53] : The heat transfer rates for hot and cold fluids are calculated as follows: • Average Nusselt number 10 : • The heat transfer coefficient 54-56 : where where q ′′ is the total heat flux that the fluid receives over the entire computational domain and is calculated as follows: The pressure drop (�P) is defined as follows: Thermal performance can be computed as:

Numerical scheme, validation, and grid independency
In the current numerical investigation, the simulations were conducted using the ANSYS Fluent 18.1 software. The governing equations were discretized using the second-order upwind technique. Besides, the SIMPLE scheme was utilized to perform the pressure-velocity coupling. The convergence metric was set to 10-6 [57][58][59][60][61][62][63] . Four different types of grids were used to make the outputs of the problem independent of the mesh. The tetrahedral mesh was used. The walls had a boundary layer mesh with a factor of 5%. The Nusselt number was used to perform the mesh study. Finally, it was found that the most appropriate mesh is the one with a 650,538 element number. To ensure the accuracy and validity of the results of the present work, the experimental findings of Wen and Ding 64 and the numerical outcomes of Goktepe et al. 65 for the local convection coefficient of NF flowing through a uniformly heated tube were employed. The outcomes are reported in Fig. 2. As can be observed, there is a good (10) Be = S gen,t S tot www.nature.com/scientificreports/ consistency between the outcomes, and the largest discrepancy between our results and the data from the other works (i.e. Refs. 64,65 ) is 10%. To further validate the simulations, the average Nusselt number of NF flow inside a channel with rectangular rib (in-line) reported by Vanaki and Mohammed 66 was compared with our results. As can be seen in Fig. 3, at low velocities, a good consistency is observed, and as the flow velocity elevates, the discrepancy of the results is elevated and reaches 7%.

Results and discussion
The main focus of this numerical investigation is to analyze the turbulent flow of water-CuO NF inside a hairpin heat exchanger with the helical baffle on the annulus side. The effects of ϕ (0-4%), Re (5000-10,000), and B (25-100 mm) on the performance metrics are investigated. Figure 4 depicts the contour plots of NF velocity and temperature for different values of B at Re = 5000 and ϕ = 4% . It is seen that as B augments, the NF velocity diminishes. Figure 5 gives the consequences of elevating Re on the velocity and temperature contours of NF for B=25 mm. As is seen, intensifying the Re entails an elevation in the NF velocity and a decrement in the NF temperature. Figure 6 gives the Nusselt number in terms of B for various ϕ . The Nusselt number improves by elevating both the Re and ϕ , while it declines by augmenting the B . For instance, at Re = 5000 and ϕ = 0% , elevating the B from 25 to 100 mm results in a 19.05% decrease in the Nusselt number of NF, while this amount for ϕ = 4% is 15.62%. In addition, at B=25 mm and ϕ = 4% , intensification of Re from 5000 to 10,000 causes a 61.83% decline in the Nusselt number of NF. Elevating the Re entails a rise in the thickness of the velocity and thermal boundary layer, which elevates the temperature gradient and, consequently, elevates the convective heat transfer coefficient and Nusselt number. Moreover, intensification of ϕ causes augmentation in the k nf and, thereby, an augmentation in the convective heat transfer and Nusselt number. Furthermore, intensifying the baffle pitch entails a decrease in the NF velocity and, therefore, a decrease in the convective heat transfer coefficient and Nusselt number.
One of the important issues when choosing a working fluid is that the fluid performance should be examined from both the heat transfer and pressure drop (pumping power) aspects. Figure 7 gives the variations of the pressure drop of NF versus ϕ and B for different Re . It is observed that the pressure drop intensifies by boosting both the ϕ and Re , while it reduces by boosting the B . For example, at Re = 5000 and ϕ = 0% , boosting the B from 25 to 100 mm results in a 73.59% reduction in the pressure drop of NF, while this amount for ϕ = 4% is 73.55%.   The results presented so far have shown that intensifying the ϕ entails an elevated Nusselt number and pressure drop, which the first is a favorable effect and the latter is undesirable. For the final decision on the usefulness of using NF in the heat exchanger under investigation, the performance index should be examined. Figure 8 demonstrates the variations of the performance index versus ϕ and B for different Re values. The outcomes reveal that the hydrodynamic performance of the water-CuO NF in the considered heat exchanger is superior to that of the pure water in just the following four cases: • B = 100 mm, ϕ = 4% and Re = 5000 where η = 1.029. In the remainder of this section, the flow of water-CuO NF with ϕ = 2% in the considered heat exchanger is examined from the perspective of irreversibility production. Figure 9 displays the influences of Re and B on www.nature.com/scientificreports/ hence a decrease in the average NF temperature which results in elevated frictional irreversibility. Elevating the B results in a decrease in the flow mixing, which results in a decrease in the velocity gradient and a decrease in the NF temperature, which in turn decreases and elevates the frictional irreversibility. It can be concluded that the decreasing effect of velocity gradient on the frictional irreversibility is overcome by the increasing effect of NF temperature and, therefore, the frictional irreversibility decreases with increasing B. Figure 10 gives the changes of thermal irreversibility of NF with ϕ = 2% versus B in terms of Re . It is seen that the thermal irreversibility declines with the elevation of Re . For example, at B = 100 mm, the intensification of www.nature.com/scientificreports/ Re from 5000 to 10,000 causes a 52.57% decrease in thermal irreversibility. As mentioned before, the NF temperature decreases with increasing Re at a constant ϕ and B , which results in elevated thermal irreversibility. Also, the augmentation of Re results in improved mixing of the flow and consequently, a decrease in the temperature gradient which ultimately results in a decrease in the thermal irreversibility. It can be said that the impact of temperature gradient on the thermal irreversibility outweighs the impact of NF temperature and therefore, the thermal irreversibility declines with the rise of Re . Moreover, Fig. 10 shows that with the rise of B at a constant Re and ϕ , the thermal irreversibility first elevates and then reduces. The highest thermal irreversibility occurs   www.nature.com/scientificreports/ at B=50 mm. Boosting of B at a constant ϕ and Re reduces the flow mixing which results in the decrease of both the NF temperature and temperature gradient which respectively elevates and diminishes the rate of thermal irreversibility. Figure 10 reveals that for the B lower than 50 mm, the increasing impact of temperature outweighs the decreasing of the temperature gradient, and ultimately, the thermal irreversibility declines while the opposite is true for the B higher than 50 mm. The results presented in Figs. 9 and 10 show that increasing the baffle pitch results in a decrease in both the thermal and frictional irreversibilities, and therefore, it can be easily deduced that increasing the B entails a declined total irreversibility. However, the influence of Re on the total irreversibility cannot be predicted because the thermal irreversibility elevates with decreasing Re and then decreases.

Conclusion
In this study, the first law and the second law of thermodynamics are employed to investigate the turbulent flow of aqueous CuO NF through a hairpin heat exchanger equipped with the helical baffle on the annulus side. The influence of Re (5000-10,000), ϕ (0-4%), and B (25-100 mm) on the performance metrics are assessed. The following results can be deduced from this simulation: • With the rise of B at a constant Re and ϕ , the thermal irreversibility first elevates and then reduces.
• Intensifying the Re entails an elevated NF velocity and a declined NF temperature.
• Pressure drop intensifies by boosting both the ϕ and Re , while it reduces by boosting the B.
• Frictional irreversibility declines and rises with increasing B and Re , respectively. • Re effect on the frictional irreversibility at lower B is greater.
In future studies, the effect of using perforated and porous turbulators on the results presented in the present study will be investigated.  www.nature.com/scientificreports/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4.