Computers, Materials & Continua DOI:10.32604/cmc.2022.019407 

Article 
Simulation of NonIsothermal Turbulent Flows Through Circular Rings of Steel
1Department of Mathematics and Social Sciences, Sukkur IBA University, Sukkur, 65200, Sindh, Pakistan
2Department of Mathematics, Faculty of Science, Khon Kaen University, Khon Kaen, 40002, Thailand
3Department of Mathematics, College of Science AlZulfi, Majmaah University, AlMajmaah, 11952, Saudi Arabia
*Corresponding Author: Kamsing Nonlaopon. Email: nkamsi@kku.ac.th
Received: 12 April 2021; Accepted: 15 June 2021
Abstract: This article is intended to examine the fluid flow patterns and heat transfer in a rectangular channel embedded with three semicircular cylinders comprised of steel at the boundaries. Such an organization is used to generate the heat exchangers with tube and shell because of the production of more turbulence due to zigzag path which is in favor of rapid heat transformation. Because of little maintenance, the heat exchanger of such type is extensively used. Here, we generate simulation of flow and heat transfer using nonisothermal flow interface in the Comsol multiphysics 5.4 which executes the Reynolds averaged Navier stokes equation (RANS) model of the turbulent flow together with heat equation. Simulation is tested with Prandtl number (Pr = 0.7) with inlet velocity magnitude in the range from 1 to 2 m/sec which generates the Reynolds number in the range of 2.2 × 105 to 4.4 × 105 with turbulence kinetic energy and the dissipation rate in ranges (3.75 × 10−3 to 1.5 × 10−2) and (3.73 × 10−3−3 × 10−2) respectively. Two correlations available in the literature are used in order to check validity. The results are displayed through streamlines, surface plots, contour plots, isothermal lines, and graphs. It is concluded that by retaining such an arrangement a quick distribution of the temperature over the domain can be seen and also the velocity magnitude is increasing from 333.15% to a maximum of 514%. The temperature at the middle shows the consistency in value but declines immediately at the end. This process becomes faster with the decrease in inlet velocity magnitude.
Keywords: Incompressible turbulent flow; heat transfer; nonisothermal; finite element method; energy equation
The implication of the circular rings or cylinders in the channel of any shape to execute the efficient heat exchanger is extensively used in the field of engineering, science and industry where the improvement of the heat transfer or heat energy is important or beneficial [1–3]. Due to cooling or heating purpose, rectangular heat exchangers have significant roles in the field of spacecraft, power generating industries, refrigeration and cryogenic procedures, food supplier industries, and petroleum industries. To magnify the heat transfer in the rectangular heat exchanger fitted with the circular tubes several fluid dynamics factors are utilized for example pressure drop, flow speed of the particles, different shapes of the geometry, obstacles of regular boundaries and curved shape, different types of flow (laminar or turbulent), different materials of fluid, etc. It is also reasonable to say that convective process, which depends on the heat transfer coefficient possesses, has a direct relationship with the fluid velocity [4–7]. In the cases when flow is examined through turbulence, velocity of the fluid remains increased due to the addition of kinetic energy of the fluid with dissipation rates [8–10]. Due to the turbulence the amount of airflow as well as the blare of the flow are enhanced which stimulates the heat conduction and thermal collaboration. This indicates to control the heat transfer to make the efficient rectangular heat exchanger, the turbulence modeling in the channel plays a vital role [11–13]. Along with the turbulence modeling the making of the heat exchanger with the air as testing fluid with addition of the nonisothermal flow properties carries a lot of importance in numerous applications like making coolant apartments, maintaining the residential heat up, producing ventilated air conditioning rooms, solar collectors, production of the turbine blades running through gaseous, etc. [14–16]. Because of having variety of applications in civil, electrical, and mechanical engineering, different scientists and researchers have put their energy to study (numerically and experimentally) the heat transfer enhancement of the physical problems.
The turbulent flow simulation along with the heat transfer over the circular cylinder was studied by [17] with the large eddy simulation using immersed boundary procedure. To get the experimental data for the turbulence near the cylinder, particle image velocimetry and hotwire anemometer was used. A large data for the turbulent flow was gathered near the premises of the cylinder. Another experimental study [18] was performed using the laminar region for the heat transfer enhancement of circular cylinder. They intended to find a certain reference temperature that yields the kinematic viscosity that can be used to predict the forces due to convection through the bluff objects. Leading the previous experimental study [19] used the numerical as well as experimental approach by using the different input temperatures with different Reynolds numbers. The study suggested the equation that shows the relationship between the temperature coefficient on the surface of the cylinder and input temperature. On observing the phenomena of vortices [20] due to the alteration of the momentum equations as well as the heat transfer, hot anemometry was used to measure the heat fluxes in the recirculating region and average momentum. With the amplification of the same study [21] observed in detail the alteration of the velocity correlations, Prandtl number for the turbulent and fluctuation of the temperature in the region. It was also added that, due to the momentum created in the region, transport of the heat was greatly affected. From the ancient, many experiments and numerical approaches [22–26] were applied to insert the circular tubes as the tabulator creator or vortex generator to enhance the heat transfer in the rectangular channel in such a way a valuable heat exchanger can be made. Also, much literature [27–29] can be found that developed the correlations among the fluid dynamics parameter and the nondimensional parameter which might help design the heat exchanger of any shape. Many of the researchers who had focused on the shape of the heat exchanger for example [30–32] have suggested that increasing the heat transfer in the tube dimpled channel proved to be the most beneficial in terms of transferring the heat than the plain tubes. From many of the articles [33–35], it was also discussed that to fit elliptical shells in the tube made the heat exchanger more efficient about 20% than fitting the circular shells in the tube! The heat exchanger is widely used to make diesel engines in mechanical engineering because of its heat energy transformation efficiency. Engineers are laboring hard to make efficient heat exchangers that resist a lot the emissions and raise the fuel capacity of the engines. There are several articles found [36–40] that use the turbulent models with different materials to work out on the heat controlling in heat exchangers.
2.1 Implementing Physics in the Model Construction for the Channel
The simulation has been expanded to explore the rectangular heat exchanger enclosed by the three circular cylinders of steels for temperature configuration and analysis of fluid flow over these cylinders. The inner radius of all cylinders and thickness is 0.25 [m] and 0.05 [m] respectively. The cylinder at the lower boundary is at the middle of the channel whereas the other two cylinders are placed at 25% and 75% of the total length of the channel at the upper boundary. Air as experimenting fluid is practiced and Flow is considered to be nonisothermal flow where the thermal properties are essentially time and temperaturedependent. For this goal, a rectangular tube is conjectured to have a length of 4 [m] and a height of 1 [m] see Fig. 1. A trial inlet velocity with magnitude in the range from 1 to 2 m/sec with the interval of 0.1 is conceded to enter with the initial temperature Tin = 232 [K] whereas the three solid circular cylinders are initially kept at constant temperature of Tp = 278 [K]. The simulation for such type of fluid flow is based on the questioning to check just for the pattern of velocity field and the temperature distribution on the circular tubes and other variables of interest. The channel is open on the right whereas the pressure is kept zero with stifling the backflow condition. This implies that the fluid and heat can enter and leave from the exit of the channel. Moreover, the symmetrical walls are well insulated.
2.2 The Meshing Process and Convergence Analysis
The selected rectangular channel has been meshed with irregular, triangular and the quad meshes. In Fig. 2a the coarse mesh is viewed which is composed of 3126 elements with edge elements of 288. To get the better convergence Fig. 2b, fine mesh is used which is consisting of 12429 elements that are connected through 641 edges and 22 vertex elements. Here, the geometry is composed of 11433 triangles and 996 quads. Here 6950 mesh vertices are used with minimum element quality of 0.1479 and average element quality of 0.8938. Near the circular cylinders we have used more fine mesh in order to have higher accuracy. For simulation using Comsol Multiphysics 5.4, 41913 degrees of freedom are solved with the internal degrees of freedom 841. Tab. 1 represents the convergence plot for the last iteration for some parameter Uin. We achieved excellent results when compared with the available literature.
2.3 Reynolds Navier Stokes Equation and the Standard
Turbulent flows are difficult to analyze due to addition of energy terms along with Navier Stokes equations. We might state that the turbulence is another form of energy. Due to the basis of complexity and addition of the energies in the fluid, the simulation through the Navier Stokes equation is controlled by time average quantities. The Reynolds averaged Navier Stokes equation is essential to predict the fluid flow for the turbulence, the idea was given by Osborne Reynolds. The equations are based on the experience for the turbulence flow which divines the approximate solution of the Navier Stokes equation. The equations are acquired with the use of the basic Navier Stokes equations when to write in the Einstein notation in the scalar form:
The derivation of the RANS model is accompanied by the Reynolds decomposition. According to his idea, each variable used to predict the flow can be disintegrated into timeaveraged component
Now substituting (3) into (1) and (2) and averaging the quantities we get the following reduced equations.
Due to the assumption, the fluctuating terms in the RANS model would be omitted. The term in
Using (7) in (6) and manipulating the terms we get:
where
In these Eq. (5) variables are chosen to be free or to optimize the results to give better accuracy. For a particular problem we pick them as:
Finally, due to the steadystate computation for the fluid flow and the heat transfer, the time derivative for all the independent quantities would be taken as zero. i.e.,
To examine the heat transfer problem, the governing steadystate heat equation is accompanied by the
Here, k is the thermal conductivity.
Parameter Computation Formula
The heat transfer coefficient will be computed by the following correlation (a famous DittusBoelter equation) for the turbulent flows:
where,
d = Hydraulic Diameter of the channel
k = Thermal conductivity
j = Mass flux of the flow per unit area
The local Nusselt numbers and the local Reynolds numbers will be computed by the formula:
3 Validation of the Existing Work with Correlations
To approve the numerical results obtained with the least Square scheme of the finite element method, the two correlations Bejan et al. [41] and Lienhard [42] are used to associate with the present results. According to the provided literature, the local Nusselt number is related to local Reynolds number and with nondimensional Prandtl number. For the current simulation of the rectangular heat exchanger fitted with three cylinders, the outcomes are estimated and compared with the correlations given by Figs. 3a and 3b at the downstream of the cylinder connected at the lower boundary. The association has been confirmed for the inlet velocities of 1, 1.5 and 2 m/sec see Figs. 3a–3c. It can be inferred that the results gained are in better accuracy with experimental results or with the correlations. It has been asserted that using the Least Square technique of Galerkin's strategy consumes less time for the computation and gives more reliable results. We assume that by practicing FEM via COMSOL Multiphysics 5.4 the results are reliable to some length. On behalf of the comparison with the correlations, we are going to provide some more relationships for the fluid dynamics parameters and the parameters of the thermodynamics.
The current fluid flow problem along with heat distribution is checked by using the turbulent
4.1 Streamlines Surface Plots of the Velocity Field, Contours of Pressure and Isothermal Lines of Temperature
The velocity field streamlines with the arrow and surface plots are exhibited in Figs. 4a–4d. For each input velocity in the domain, the velocity at the outlet is improved by a maximum value of 333.5% see Fig. 4d. The rate of flow is enhanced and move in the recurrent direction of flow. Such an adjustment of the flow will be convenient where the heat energy at the outlet is maximized. The circular cylinders generate restraint in the way of flow and create the hydraulic jumps. Hydraulic jumps after one and another period generate more kinetic energy in the flow. Thus, a dynamic growth in the velocity of the flow in terms of percentage obtained of 333.5% when to elaborate on the basis of inlet velocity. It is the part of the turbulent flows that they are often depending on the velocity given at the inlet of the channel. It is manifest in Fig. 4a. with the given velocity of the 1 m/sec, the velocity at the outlet is extended by 4.26 m/sec. In Fig. 4b. with the given inlet velocity of 1.4 m/sec, the velocity at exit is increasing to 6.14 m/sec. That is the gain of about 514% of the inlet velocity. The minimum and the maximum values in the domain is given by the Tab. 2.
In Figs. 5a–5d the contour plots for the pressure distribution over the domain is presented for the particular velocities at the inlet of 1 m/sec, 1.4 m/sec, 1.8 m/sec and 2 m/sec. For each input velocity in the simulation the flow starts with maximum magnitude in the domain. As the fluid moves to the exit of the channel the pressure is declining down up to a negative pressure. Production of the negative pressure will support the controlling of the fluid inside the pipe. The negative pressure is mostly producing at the end of the channel on the upper corner beside the cylinder. The rapid variation of the pressure can be seen on the space bounded by the three circles due to the rapid increase of the flow rate around the corners of cylinders. Tab. 3 presents the numerical computed results for the minimum and maximum pressure on the whole domain. It is clear that with the increase in inlet velocity the minimum pressure around the domain is decreased whereas the maximum pressure is increased with the increase in inlet velocity magnitude. The portion where the minimum pressure exits may located at the end of the corners of the channel whereas the maximum pressure can be located portion near the inlet of the channel before the first circular cylinder.
The isothermal contours are exhibited in Figs. 6a–6d for some selected inlet velocities. They are determining that the distribution of the temperature is not uniform across all the domains. But with the rise in inlet velocity, the extent of the temperature intensity is increasing over the domain. Originally, the cylinders are kept at the temperature of 278 [k] whereas the air infiltrates the region bearing the temperature of 323 [k]. As soon as the fluid is passed to enter the region the distribution of heat flows quickly over the domain. The method of distribution becomes quick with increasing the velocity magnitude at the exit of the channel. The minimum and maximum temperatures at the domain of the channel are given in Tab. 4. With the rise of the input velocity, the minimum temperature is growing but there will be no impact on the maximum temperature. The minimum temperature can be perceived at the region surrounded by the circle or at the region after the circle attached at the end see Figs. 6a–6d.
4.2 Variation of the Temperutre on the Surface of the Cylinder
In the portion we are working to explain the changing of temperature over the surface circular cylinder attached at 75% of the length of the channel against the velocity magnitude, pressure, local Nusselt number and the Reynolds number see Figs. 7a–7d. The temperature of the surface is increasing with the increase in the flow rate of the fluid. The speed of the flow at the center of the cylinder attained maximum and then the temperature declined slowly with the reduction of velocity field magnitude see Fig. 7a. With the increase of pressure on the surface of the circular cylinder, the temperature confers a positive response. With the increase of inlet velocities, the connection between the temperature and pressure bestows more stability and response in the terms of enrichment. The pressure on the surface of the circular cylinder diversifies in a wide range when the velocity magnitude at the inlet is a little boosted.
The convective process is lagging over the surface of the cylinder for each inlet velocity. See Fig. 7c which concludes that the temperature over the cylinder is declining with the increase in local Nusselt number. With enhancing the velocity magnitude at the inlet, the local Nusselt number shows the direct correlation with the temperature. For the current problem of fluid flow and heat transfer, the local Nusselt number varies from 50 to 400 when the inlet velocities are varying in range from 1 to 2. The local Reynolds number behaves the same as the local Nusselt number at the surface of the cylinder. The influence on the temperature of the local Reynolds is identical and with the raising the range of Reynolds number the temperature on the surface decays in the same way as the local Nusselt number does see Figs. 7c and 7d. With the increase in the velocity at the opening of the channel, the range in the local Reynolds is increased with the new derivation of the initial temperature on the surface of the cylinder. It is also apparent that due to the declination of local Reynolds numbers the temperature at the appropriate inlet velocity also operates down with the common pattern for all Reynolds numbers.
4.3 Pattern of the Flow in the Middle Horizontal Line of the Channel
It is to assert that while the fluid is departing the region the whole model of the flow can be analyzed in the best fashion in the middle of the channel. Because the motion of the fluid is adjusted in that circumstantial direction. In the part we are going to perceive the pattern of flow for u and v components of velocity pressure and the temperature at the middle of the channel in terms of presented input velocities see Figs. 8a–8d. It is shown in Fig. 8a when the air is deducted to flow the ucomponent of the velocity is equal to inlet velocity because of the horizontal orientation of the channel. For a particular inlet velocity Uin = 1–2 m/sec, the ucomponent of velocity is raised while passing the narrow paths enclosed by the three cylinders. It endeavors some optimum values and then declines.
It is also noticeable that the greater the velocity greater the achievement of optimum values of the flow rate over the domain. While inspecting the velocity in ydirection it starts from zero for all given inlet velocities. The ycomponent of velocity moves in oscillating paths and attempts some minimum values before arriving at the exit of the channel. From Figs. 8a–8b it is inferred that the fluid is mostly stretched towards horizontal paths due to the contraction of the circles. It is apparent from Fig. 8c, the greater input velocity will start the flow with greater pressure but the pressure wanes to zero at the end of the channel. Also, the rate of drop of the pressure is directly proportional to the velocity magnitude at the inlet of the channel. It can be recommended that to start a flow with maximum acceleration a pressure of very less value must be forced at the end of the channel. The temperature distribution acts at the middle with the constant temperature of 323 [k] unless the fluid particles passed by the region molded by the three circles see Fig. 8d. The rapid declination of the temperature can be observed near the exit of the channel refer Fig. 8d. Fig. 8d is also intimating that the less the inlet velocity greater the rapid declination in the temperature down the exit of the channel.
The heat flow with the help of air flow was discussed with the use of nonisothermal interface in the finiteelement based software Comsol Multiphysics 5.4 which implements the Reynolds averaged Navier Stokes model for the
(1) The discussion of the streamlines pattern and the surface plots indicated that the velocity magnitude is increasing in the range of 333.5% to 514% and due to such arrangement of the cylinders creates the more hydraulic jumps to accelerate the fluid motion. The Rapid variation in the pressure can be seen in the region which is surrounded by the cylinder. For the current problem of the turbulent flow quick distribution of the temperature was seen in the domain.
(2) When to express the distribution of the temperature over the circular cylinder against velocity it is found that temperature is maximum at the middle of the cylinder. With the increase in the pressure the temperature of the cylinder is increasing for each input velocity. The range of the pressure is also increasing with the increase in input velocity field. Temperature on the surface of the circular cylinder declines due to increase in local Nusselt number as well as the Reynolds number for each input velocity. But for each input velocity the range for local Reynolds number as well as the local Nusselt number is increased.
(3) Since the channel is horizontal therefore the xcomponent of the velocity plays a major role in developing the flow at the start of the channel in the middle line and increase in magnitude during the whole flow. Whereas the ycomponent of the velocity is starting from zero position and remain fluctuate in the domain. The pressure at the middle line is increasing with inlet velocity magnitude and will be zero at the exit of the channel. The temperature in the middle line remains at the same value given at the inlet and declines quickly before reaching at the domain.
Acknowledgement: Authors are very much thankful to Sukkur IBA University, Sukkur, Pakistan for providing conducive environment for research.
Funding Statement: The authors received no specific funding for this study.
Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.
References
1. A. M. Jacobi and R. K. Shah, “Heat transfer surface enhancement through the use of longitudinal vortices: A review of recent progress,” Experimental, Thermal and Fluid Science, vol. 11, no. 3, pp. 295–309, 1995.
2. M. H. Mousa, N. Miljkovic and K. Nawaz, “Review of heat transfer enhancement techniques for single phase flows,” Renewable and Sustainable Energy Reviews, vol. 137, pp. 110566, 2021.
3. C. Maradiya, J. Vadher and R. Agarwal, “The heat transfer enhancement techniques and their thermal performance factor,” BeniSuef University Journal of Basic and Applied Sciences, vol. 7, no. 1, pp. 1–21, 2018.
4. Y. Sun, W. J. Jasper and E. A. DenHartog, “Effects of air velocity, air gap thickness and configuration on heat transfer of a wearable convective cooling system,” Journal of Textile Science & Engineering, vol. 5, no. 6, pp. 1000227, 2015.
5. A. S. Yang, Y. C. Shih, C. L. Lee and M. C. Lee, “Investigation of flow and heat transfer around internal channels of an air ventilation vest,” Textile Research Journal, vol. 84, no. 4, pp. 399–410, 2014.
6. S. Okamoto and Y. Sunabashiri, “Vortex shedding from a circular cylinder of finite length placed on a ground plane,” Journal of Fluids Engineering, vol. 114, no. 4, pp. 512–521, 1992.
7. W. M. Abed, R. D. Whalley, D. J. C. Dennis and R. J. Poole, “Experimental investigation of the impact of elastic turbulence on heat transfer in a serpentine channel,” Journal of NonNewtonian Fluid Mechanics, vol. 231, pp. 68–78, 2016.
8. D. Zhang, L. Cheng, H. An and S. Draper, “Flow around a surfacemounted finite circular cylinder completely submerged within the bottom boundary layer,” European Journal of MechanicsB/Fluids, vol. 86, pp. 169–197, 2021.
9. L. RedjemSaad, M. OuldRouiss and G. Lauriat, “Direct numerical simulation of turbulent heat transfer in pipe flows: Effect of prandtl number,” International Journal of Heat and Fluid Flow, vol. 28, no. 5, pp. 847–861, 2007.
10. M. Piller, “Direct numerical simulation of turbulent forced convection in a pipe,” International Journal for Numerical Methods in Fluids, vol. 49, no. 6, pp. 583–602, 2005.
11. P. Constantin and C. R. Doering, “Heat transfer in convective turbulence,” Nonlinearity, vol. 9, no. 4, pp. 1049, 1996.
12. J. C. Han, J. S. Park and C. K. Lei, “Heat transfer enhancement in channels with turbulence promoters,” Journal of Engineering in Gas Turbines and Power, vol. 107, no. 3, pp. 628–635, 1985.
13. G. W. Lowery and R. I. Vachon, “The effect of turbulence on heat transfer from heated cylinders,” International Journal of Heat and Mass Transfer, vol. 18, no. 11, pp. 1229–1242, 1975.
14. K. Boulouchos, “Energy Research in the ETH Domain: Science and Technology for Sustainable Energy,” ETH Zurich, 2005. https://doi.org/10.3929/ethza004953211.
15. T. Arghand, “Direct ground cooling systems for office buildings,” Ph.D. dissertation, Chalmers University of Technology, Gothenburg, Sweden, 2021.
16. M. Sheikholeslami, M. GorjiBandpy and D. D. Ganji, “Experimental study on turbulent flow and heat transfer in an air to water heat exchanger using perforated circularring,” Experimental Thermal and Fluid Science, vol. 70, pp. 185–195, 2016.
17. P. Parnaudeau, J. Carlier, D. Heitz and E. Lamballais, “Experimental and numerical studies of the flow over a circular cylinder at reynolds number 3900,” Physics of Fluids, vol. 20, no. 8, pp. 085101, 2008.
18. A. B. Wang, and Z. Trávnıček, “On the linear heat transfer correlation of a heated circular cylinder in laminar crossflow using a new representative temperature concept,” International Journal of Heat and Mass Transfer, vol. 44, no. 24, pp. 4635–4647, 2001.
19. L. Baranyi, S. Szabó, B. Bolló and R. Bordás, “Analysis of low reynolds number flow around a heated circular cylinder,” Journal of Mechanical Science and Technology, vol. 23, no. 7, pp. 1829–1834, 2009.
20. M. Matsumura and R. A. Antonia, “Momentum and heat transport in the turbulent intermediate wake of a circular cylinder,” Journal of Fluid Mechanics, vol. 250, pp. 651–668, 1993.
21. R. A. Antonia, Y. Zhou and M. Matsumura, “Spectral characteristics of momentum and heat transfer in the turbulent wake of a circular cylinder,” Experimental, Thermal and Fluid Science, vol. 6, no. 4, pp. 371–375, 1993.
22. A. I. ElSherbini and A. M. Jacobi, “The thermalhydraulic impact of deltawing vortex generators on the performance of a plainfinandtube heat exchanger,” HVAC&R Research, vol. 8, no. 4, pp. 357–370, 2002.
23. C. C. Wang, Y. J. Chang, C. S. Wei and B. C. Yang, “A comparative study of the airside performance of winglet vortex generator and wavy Finandtube heat exchangers,” ASHRAE Transactions, vol. 110, no. 1, pp. 53–57, 2004.
24. W. W. Cheng and C. Madhusudana, “Effect of electroplating on the thermal conductance of fin–tube interface,” Applied Thermal Engineering, vol. 26, no. 17–18, pp. 2119–2131, 2006.
25. D. Tang, D. Li, Y. Peng and Z. Du, “A new approach in evaluation of thermal contact conductance of tube–fin heat exchanger,” Applied Thermal Engineering, vol. 30, no. 14–15, pp. 1991–1996, 2010.
26. M. C. Gentry and A. M. Jacobi, “Heat Transfer Enhancement Using tip and Junction Vortices,” Air Conditioning and Refrigeration Center, College of Engineering, University of Illinois, Urbana, 1998.
27. M. E. Nakhchi, M. Hatami and M. Rahmati, “Experimental investigation of heat transfer enhancement of a heat exchanger tube equipped with doublecut twisted tapes,” Applied Thermal Engineering, vol. 180, pp. 115863, 2020.
28. H. Cai, S. Liang, C. Guo, T. Wang, Y. Zhu et al., “Numerical investigation on heat transfer of supercritical carbon dioxide in the microtube heat exchanger at low reynolds numbers,” International Journal of Heat and Mass Transfer, vol. 151, pp. 119448, 2020.
29. P. Ocłoń, S. Łopata, T. Stelmach, M. Li, J. F. Zhang et al., “Design optimization of a hightemperature finandtube heat exchanger manifold–a case study,” Energy, vol. 215, pp. 119059, 2021.
30. S. Xie, Z. Liang, L. Zhang and Y. Wang, “A numerical study on heat transfer enhancement and flow structure in enhanced tube with cross ellipsoidal dimples,” International Journal of Heat and Mass Transfer, vol. 125, pp. 434–444, 2018.
31. S. Xie, Z. Liang, L. Zhang, Y. Wang, H. Ding et al., “Numerical investigation on heat transfer performance and flow characteristics in enhanced tube with dimples and protrusions,” International Journal of Heat and Mass Transfer, vol. 122, pp. 602–613, 2018.
32. S. Xie, Z. Liang, J. Zhang, L. Zhang, Y. Wang et al., “Numerical investigation on flow and heat transfer in dimpled tube with teardrop dimples,” International Journal of Heat and Mass Transfer, vol. 131, pp. 713–723, 2019.
33. R. S. Matos, T. A. Laursen, J. V. C. Vargas and A. Bejan, “Threedimensional optimization of staggered finned circular and elliptic tubes in forced convection,” International Journal of Thermal Sciences, vol. 43, no. 5, pp. 477–487, 2004.
34. R. S. Matos, J. V. C. Vargas, T. A. Laursen and A. Bejan, “Optimally staggered finned circular and elliptic tubes in forced convection,” International Journal of Heat and Mass Transfer, vol. 47, no. 6–7, pp. 1347–1359, 2004.
35. R. S. Matos, J. V. C. Vargas, T. A. Laursen and F. E. M. Saboya, “Optimization study and heat transfer comparison of staggered circular and elliptic tubes in forced convection,” International Journal of Heat and Mass Transfer, vol. 44, no. 20, pp. 3953–3961, 2001.
36. B. Farajollahi, S. Gh. Etemad and M. Hojjat, “Heat transfer of nanofluids in a shell and tube heat exchanger,” International Journal of Heat and Mass Transfer, vol. 53, no. 1–3, pp. 12–17, 2010.
37. E. Ozden and I. Tari, “Shell side CFD analysis of a small shellandtube heat exchanger,” Energy Conversion and Management, vol. 51, no. 5, pp. 1004–1014, 2010.
38. A. A. A. Arani and R. Moradi, “Shell and tube heat exchanger optimization using new baffle and tube configuration,” Applied Thermal Engineering, vol. 157, pp. 113736, 2019.
39. S. G. Etemad, B. Farajollahi, M. Hajipour and J. Thibault, “Turbulent convective heat transfer of suspensions of ΓAl2O3 and CUO nanoparticles (nanofluids),” Journal of Enhanced Heat Transfer, vol. 19, no. 3, pp. 191–197, 2012.
40. D. Kim, Y. Kwon, Y. Cho, C. Li, S. Cheong et al., “Convective heat transfer characteristics of nanofluids under laminar and turbulent flow conditions,” Current Applied Physics, vol. 9, no. 2, pp. e119–e123, 2019.
41. A. Bejan and A. D. Kraus, Heat Transfer Handbook, vol. 1. Hoboken, New Jersey: John Wiley & Sons, Inc., 2003.
42. J. H. Lienhard, “Exterior shape factors from interior shape factors,” Journal of Heat Transfer, vol. 141, no. 6, pp. 0613011–06130110, 2019.
This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 