A CFD Study on Heat Transfer Performance of SiO2-TiO2 Nanofluids under Turbulent Flow

A CFD model was performed with commercial software through the adoption of the finite volume method and a SIMPLE algorithm. SiO2-P25 particles were added to water/ethylene glycol as a base fluid. The result is considered a new hybrid nanofluid (HN) for investigating heat transfer (HT). The volume concentrations were 0.5, 1.0, and 1.5%. The Reynolds number was in the range of 5000–17,000. The heat flux (HF) was 7955 W/m2, and the wall temperature was 340.15 K. The numerical experiments were performed strictly following the rules that one should follow in HT experiments. This is important because many studies related to nanofluid HT overlook these details. The empirical correlations that contain the friction factor perform better with higher Reynolds numbers than the correlations based only on Reynolds and Prandtl numbers. When temperature differences are moderate, researchers may consider using constant properties to lower computational costs, as they may give results that are similar to temperature-dependent ones. Compared with previous research, our simulation results are in agreement with the experiments in real time.


Introduction
Common liquids such as water, ethylene glycol, thermal oil, etc., are used in HT devices for different industrial applications. In the last few years, many studies have improved the HT to maintain more suitable HT in different thermal systems [1,2]. The hybrid nanofluid is a novel fluid with a hybrid nanofluid for high thermal conductivity (k) and stability. Much research has been conducted on HT through a tube with different nanofluids [3,4].
Hussein et al. investigated the impact of the cross area on HT and friction [5]. They found that, compared with circular and elliptical pipes, flat pipes have a notable increment of HT and pressure drop (PD). Additionally, they studied the augmentation of forced convection HT in an automotive cooling system [6].

Theoretical Background
The HT research was conducted in at least two directions. The first was to perform experiments and then find and construct empirical and semi-empirical formulae (correlations) for solving engineering problems. The second was to analyse the HT phenomenon, investigating its physics, discovering the effects of working fluid properties, the wall-fluid interaction, and the validity of the known equations applied to describe the process. The results that yield the solutions for HT problems are based on a combination of experimental data and theoretical analysis. Earlier, physical experiments were the only way to determine HT subprocesses. Recent sophisticated CFD software offers reasonable alternatives, the so-called numerical experiments, for investigating HT problems related by determining the velocity and temperature fields in different arrangements. Independent of the type of experiment, the dimensionless numbers or groups resulting from the similarity theories (Buckingham pi theorem or differential similarity) prove crucial in HT calculations. Boussinesq assumptions are generally applied in solving the governing equations of the HT (mass, energy, momentum conservation). The matter properties variations are indirectly set up, usually by the proportion of the Prandtl numbers at two characteristic temperatures (average fluid and wall). Neither the differential similarity nor the pi theorem applies if the variation in properties is included in the fundamental equations.
The most investigated and basic arrangement is the horizontal tube with an isothermal wall or steady wall HF. One can also use these configurations in engineering formulae for HT calculation. It should be mentioned that the isothermal wall temperature, in most cases, can be maintained only in strict experimental conditions, and in practice, this assumption rarely exists. Nevertheless, the difference between the two basic approaches is nearly negligible for the Nu number in a turbulent flow. Promoting the numerical experiments (making CFD calculation), the boundary conditions can easily be established, and the results of the calculations are according to the desired conditions. That is one possible explanation for the observed variance, compared with the real experimental results.

Forced Flow in Circular Tubes
Forced flow occurs in multiple configurations, and one of the important engineering arrangements is the forced flow in the duct. The duct's cross section may appear in various shapes, and in the engineering calculation, it is maintained through an equivalent circular diameter. Therefore, investigations generally deal with only the circular tubes.
Although the tests were carried out on a circular duct, it is advisable to recall that, in practice, the non-circular ducts cases can be handled well using the hydraulic (D eqH ) and thermal (D eqT ) equivalent diameter even in complex permanent cross-sectional structures. The definitions are as follows: The hydraulic equivalent diameter should be used in the Reynolds number and the thermal in the Nusselt number.
According to the Re number, the flows in a tube are characterised as laminar flow (Re < 2300), transitional flow (2300 < Re < 4000), and turbulent flow (4000 < Re). However, transitional flow is often considered part of turbulent flow in engineering practice. Moreover, depending on the distance from the entrance, there are different flow formations such as developing and (fully) developed flow regimes, as shown in Figure 1. Therefore, the local HT varies alongside the tube length, and some correlation formulae count this phenomenon. As both the velocity and the temperature vary from those of the entrance, there are hydraulically developing, thermally developing, and simultaneously developing flows. Nanomaterials 2022, 12, x FOR PEER REVIEW  3 of 17 cases, can be maintained only in strict experimental conditions, and in practice, this assumption rarely exists. Nevertheless, the difference between the two basic approaches is nearly negligible for the Nu number in a turbulent flow. Promoting the numerical experiments (making CFD calculation), the boundary conditions can easily be established, and the results of the calculations are according to the desired conditions. That is one possible explanation for the observed variance, compared with the real experimental results.

Forced Flow in Circular Tubes
Forced flow occurs in multiple configurations, and one of the important engineering arrangements is the forced flow in the duct. The duct's cross section may appear in various shapes, and in the engineering calculation, it is maintained through an equivalent circular diameter. Therefore, investigations generally deal with only the circular tubes.
Although the tests were carried out on a circular duct, it is advisable to recall that, in practice, the non-circular ducts cases can be handled well using the hydraulic (DeqH) and thermal (DeqT) equivalent diameter even in complex permanent cross-sectional structures. The definitions are as follows: The hydraulic equivalent diameter should be used in the Reynolds number and the thermal in the Nusselt number.
According to the Re number, the flows in a tube are characterised as laminar flow (Re < 2300), transitional flow (2300 < Re < 4000), and turbulent flow (4000 < Re). However, transitional flow is often considered part of turbulent flow in engineering practice. Moreover, depending on the distance from the entrance, there are different flow formations such as developing and (fully) developed flow regimes, as shown in Figure 1. Therefore, the local HT varies alongside the tube length, and some correlation formulae count this phenomenon. As both the velocity and the temperature vary from those of the entrance, there are hydraulically developing, thermally developing, and simultaneously developing flows.

Turbulent Flow
The inner perturbation of fluid causes a turbulent flow. Different perturbations always exist in the fluid. When velocity is low, the viscous forces stabilise the flow pattern, but this stabilisation is not effective when the velocity increases and the local vortex formations result in a turbulent flow pattern. The laminar-turbulent transition never means a sharp borderline. There is always a particular transition zone. This phenomenon explains why the HT coefficient (Nu number) calculations perform less when the Re number is not much larger than the laminar Re number limit. There is always a particular length

Turbulent Flow
The inner perturbation of fluid causes a turbulent flow. Different perturbations always exist in the fluid. When velocity is low, the viscous forces stabilise the flow pattern, but this stabilisation is not effective when the velocity increases and the local vortex formations result in a turbulent flow pattern. The laminar-turbulent transition never means a sharp borderline. There is always a particular transition zone. This phenomenon explains why the HT coefficient (Nu number) calculations perform less when the Re number is not much larger than the laminar Re number limit. There is always a particular length after the fully developed flow pattern is generated. Since this study subject is the turbulent forced flow, the Re number is greater than 5000 in all investigated cases.
In Handbook of Heat Transfer [23], there are 12 formulae for calculating the Nu number in the smooth tube for turbulent flow. The wide variation of the starting range for the Re numbers shows that generating a formula that accounts for the transition and the fully developed regime is not a simple task. The formulae can be characterised in two groups-either (a) the ξ friction factor is not included in the formulae or (b) it is, as shown in Tables 1 and 2. Analysing all the formulae is beyond the scope of this study; therefore, we selected only three formulae to validate our numerical experiments for which the application range should fit 5000 ≤ Re and 3 ≤ Pr. The Dittus-Boelter (1) is generally used, while Sieder-Tate (6) accounts for the property variation, and Gnielinski (13) is based on the friction factor and also counts the tube length diameter ratio. The arrangements of the experiments should follow the limitations stated in the application ranges of the correlations considering flow development (hydraulic, thermal, and simultaneous). The other listed correlations might be useful for readers in further investigation or application of Nu number calculations.  [23].

Correlations (a) Application Range
Dittus-Boelter (1)       According to the second law of thermodynamics, any methods or efforts made to intensify the HT result in an entropy generation, which manifests in the increased amount of work that must be performed to maintain flow configuration. The calculation of that increased work-practically the pumping work-is based on the determination of the friction factor of the flow arrangement, smoothness of the tube, and viscosity of the working fluid. When the working fluid property is subject to change, increasing the HT, this phenomenon-namely, the increased friction, always appears. The increased viscosity is the cost for the increased HT as the second law of thermodynamics states. The friction factor calculation might also be carried out according to several formulae. As the Prandtl-Kármán-Nikuradse friction factor definition is implicit, several value calculations are recommended and listed in Table 3 to overcome that difficulty. Table 3. Correlations for friction factor, [23].

Correlations
Application Range

Fluid Properties
SiO 2 -P25 nanoparticles and water/ethylene glycol mixture as base fluid were considered hybrid nanofluids (HNs), and the authors investigated these properties. The HNs in point were considered incompressible, Newtonian, and single-phase fluids. Table 4 illustrates the density (ρ), k, dynamic viscosity (µ), and specific heat (c p ) of HNs obtained from the authors' previous study [24]. The measurements were conducted at the temperature range from 20 • C to 60 • C. Some empirical formulae in Tables 1 and 2 use the Pr numbers ratio as temperature correction. Therefore, the prediction was applied for the HN's Pr numbers extension to 100 • C, as Table 5 reports. The prediction method can be summarised as follows: The ratio of the nano fluid's Pr number against the water's Pr number in the measurement range was calculated and extended by linear regression. Figure 2 visualises the measured and extended values. As shown, the extended Pr numbers are the natural continuity of the measured data as the third-order polynomial curve R 2 values confirm. If needed, this prediction method works for other properties as well.  The empirical formulae should always follow the specific order related to the formulae. Generally, the characteristic length is specified, and the characteristic temperature is defined so that one should use it to determine fluid properties. When the formula has an extended application range, the Re and Pr numbers appear, and specific corrections may also count the specialties of the arrangement. Such correction could be the length-to-diameter ratio or the Pr numbers ratio. Figure 3 illustrates the effect that the properties variation may mean. The empirical formulae should always follow the specific order related to the formulae. Generally, the characteristic length is specified, and the characteristic temperature is defined so that one should use it to determine fluid properties. When the formula has an extended application range, the Re and Pr numbers appear, and specific corrections may also count Nanomaterials 2022, 12, 299 7 of 16 the specialties of the arrangement. Such correction could be the length-to-diameter ratio or the Pr numbers ratio. Figure 3 illustrates the effect that the properties variation may mean. The empirical formulae should always follow the specific order related to the formulae. Generally, the characteristic length is specified, and the characteristic temperature is defined so that one should use it to determine fluid properties. When the formula has an extended application range, the Re and Pr numbers appear, and specific corrections may also count the specialties of the arrangement. Such correction could be the length-to-diameter ratio or the Pr numbers ratio. Figure 3 illustrates the effect that the properties variation may mean. Since the mean fluid temperature plays a central role in the HT coefficient calculation, Equation (3) shows its calculation.
where Tmean,f practically is the so-called completely stirred temperature. Since the mean fluid temperature plays a central role in the HT coefficient calculation, Equation (3) shows its calculation.
where T mean,f practically is the so-called completely stirred temperature.

Arrangements of the Numerical Experiment
The inlet temperature of the flow was 30 • C. The water was used for validating the CFD analysis. The inlet velocity was uniform. Static pressure conditions at the entrance and stagnation pressure conditions at the outlet were applied. The fluid's velocity and temperature were varied so that the Reynolds number varied from 5000 to 17,000, and the flow was turbulent.
The CFD study applied (a) constant HF of 7955 W/m 2 and (b) a constant tube wall temperature of 67 • C. Figure 4 represents the arrangement, flow in a pipe with L = 1.5 m length, d = 16 mm hydraulic diameter, and wt = 2 mm wall thickness. We assumed that the flow in the tube is symmetrical and steady. The following equations were used for calculating the mass, energy, and momentum conservation with variable properties [26]: Nanomaterials 2022, 12, x FOR PEER REVIEW 8 of 17

Arrangements of the Numerical Experiment
The inlet temperature of the flow was 30 °C. The water was used for validating the CFD analysis. The inlet velocity was uniform. Static pressure conditions at the entrance and stagnation pressure conditions at the outlet were applied. The fluid's velocity and temperature were varied so that the Reynolds number varied from 5000 to 17,000, and the flow was turbulent.
The CFD study applied (a) constant HF of 7955 W/m 2 and (b) a constant tube wall temperature of 67 °C. Figure 4 represents the arrangement, flow in a pipe with L = 1.5 m length, d = 16 mm hydraulic diameter, and wt = 2 mm wall thickness. We assumed that the flow in the tube is symmetrical and steady. The following equations were used for calculating the mass, energy, and momentum conservation with variable properties [26]:   ANSYS Workbench 16 (ANSYS Inc., Canonsburg, PA, USA) was used for this simulation research by solver strategy. The control volume approach provided the governing equations for the one-phase flow of the HNs. The geometry of the studied object was simulated as a circular duct (pipe), then the meshing was created based on the geometry, and boundary conditions were set up before the solving phase started. After the convergence reached the given criteria in the post-process phase, the HT and pressure drop were obtained at the computational domain. The applied mesh configuration is illustrated in Figure 5.
= 0 (4) ANSYS Workbench 16 (ANSYS Inc., Canonsburg, PA, USA) was used for this simulation research by solver strategy. The control volume approach provided the governing equations for the one-phase flow of the HNs. The geometry of the studied object was simulated as a circular duct (pipe), then the meshing was created based on the geometry, and boundary conditions were set up before the solving phase started. After the convergence reached the given criteria in the post-process phase, the HT and pressure drop were obtained at the computational domain. The applied mesh configuration is illustrated in Figure 5.

Grid Independence Test
The grid selection was tested using the result of the grid independence investigation. The investigated Re range was from 5000 to 17,000. The Nu numbers by the ANSYS calculation were determined for the Re = 5000 and Re = 17,000, with the mesh sizes of 7.0 × 10 5 , 8.0 × 10 5 , 9.0 × 10 5 , 1.0 × 10 6 , 1.1 × 10 6 , and 1.2 × 10 6 . Figure 6 shows the calculated ratio (Nu divided by Nu of the highest mesh) variation against the mesh density. We assumed that if mesh independence is proved at the end of the investigated range, the same mesh used is validated inside the Re range. The grid selection was tested using the result of the grid independence investigation. The investigated Re range was from 5000 to 17,000. The Nu numbers by the ANSYS calculation were determined for the Re = 5000 and Re = 17,000, with the mesh sizes of 7.0 × 10 5 , 8.0 × 10 5 , 9.0 × 10 5 , 1.0 × 10 6 , 1.1 × 10 6 , and 1.2 × 10 6 . Figure 6 shows the calculated ratio (Nu divided by Nu of the highest mesh) variation against the mesh density. We assumed that if mesh independence is proved at the end of the investigated range, the same mesh used is validated inside the Re range.

Fully Developed Flow
In the formulae that account for the fully developed flow HT, the diameter-to-length ratio did not appear, except in (7) and (13), but (7) was not applicable due to its validity for low Pr numbers. The numerical experiment had to be performed for a long tube. The local HT variation was localised for a shorter length in turbulent flow configuration than the laminar flow; for calculating the entry length (Le), we used Equation (6) from [27]. However, in planning the experiment, the tube should be long enough to neglect that part. Figure 7 represents the local Nu number variation decrease from [23]. According to the investigated NFs, the Pr = 5 curve was assumed to determine that part of the length where the local Nu number varied. Practically, this means that our investigated zone started at =~23 · = 0.37 m from the tube entrance, as shown in Table 6, which lists the calculation results based on [27]. It is also worth mentioning that in [28], the proposed equation for hydrodynamic entry length is = • 1.359 • . and in [29], =~10 • . The smallest approximation for entry length is in [27][28][29] offers the largest one. We accepted

Fully Developed Flow
In the formulae that account for the fully developed flow HT, the diameter-to-length ratio did not appear, except in (7) and (13), but (7) was not applicable due to its validity for low Pr numbers. The numerical experiment had to be performed for a long tube. The local HT variation was localised for a shorter length in turbulent flow configuration than the laminar flow; for calculating the entry length (L e ), we used Equation (6) from [27]. However, in planning the experiment, the tube should be long enough to neglect that part. Figure 7 represents the local Nu number variation decrease from [23]. According to the investigated NFs, the Pr = 5 curve was assumed to determine that part of the length where the local Nu number varied. Practically, this means that our investigated zone started at L = ∼ 23·D = 0.37 m from the tube entrance, as shown in Table 6, which lists the calculation results based on [27]. It is also worth mentioning that in [28], the proposed equation for hydrodynamic entry length is Le = D·1.359·Re 0.25 and in [29], Le = ∼ 10·D. The smallest approximation for entry length is in [27][28][29] offers the largest one. We accepted this latter value.

Velocity and Temperature Profiles
The velocity and temperature profiles of the HNs at the fully developed flow regime are represented in Figure 8 as an example for Re = 10,000. Part (a) shows wall subjects of constant HF, while part (b) shows wall subjects of constant temperature. According to sources in the literature, for example [23], in turbulent flow, there are no countable differences between part (a) and part (b) regarding the velocity profiles, but the temperature values are different; meanwhile, the shapes of the temperature profiles are also very similar.

Results
The following sections introduce our results based on the ANSYS calculation. The results are valid according to the limitations of the previous sections' reports. Here, we mention again that the developing flow was not included in the results, compared with the correlations, because most are valid only for infinite length tubes. Figure 9 represents the numerical simulation results for a constant wall HF of Re = 5000…17,000. In agreement with sources in the literature, the 1.5% concentration HN has the highest Nu numbers.

Velocity and Temperature Profiles
The velocity and temperature profiles of the HNs at the fully developed flow regime are represented in Figure 8 as an example for Re = 10,000. Part (a) shows wall subjects of constant HF, while part (b) shows wall subjects of constant temperature. According to sources in the literature, for example [23], in turbulent flow, there are no countable differences between part (a) and part (b) regarding the velocity profiles, but the temperature values are different; meanwhile, the shapes of the temperature profiles are also very similar.

Velocity and Temperature Profiles
The velocity and temperature profiles of the HNs at the fully developed flow regim are represented in Figure 8 as an example for Re = 10,000. Part (a) shows wall subjects constant HF, while part (b) shows wall subjects of constant temperature. According sources in the literature, for example [23], in turbulent flow, there are no countable diff ences between part (a) and part (b) regarding the velocity profiles, but the temperatu values are different; meanwhile, the shapes of the temperature profiles are also very si ilar.

Results
The following sections introduce our results based on the ANSYS calculation. T results are valid according to the limitations of the previous sections' reports. Here, w mention again that the developing flow was not included in the results, compared w

Results
The following sections introduce our results based on the ANSYS calculation. The results are valid according to the limitations of the previous sections' reports. Here, we mention again that the developing flow was not included in the results, compared with the correlations, because most are valid only for infinite length tubes.
3.5.1. Constant HF with Non-Variable Properties, Developed Flow Figure 9 represents the numerical simulation results for a constant wall HF of Re = 5000 . . . 17,000. In agreement with sources in the literature, the 1.5% concentration HN has the highest Nu numbers.      do not contain the friction factor, while the vonK and Gn include the friction factor. Since the tendencies are the same for all Re numbers in the investigated range, only the Re = 5000 and Re = 17,000 results are shown. The major observation can be summarised as follows:

•
The Nu numbers of the constant wall temperature are 7-10% higher than the Nu numbers of the constant HF, and the difference increases as the Re number increases. The difference for constant HF and constant wall temperature is considered small in the literature, for example, in [23], but due to the greater temperature difference alongside the tube length, the HT is more intensive in the case of constant wall temperature. Moreover, the observed differences in the context of the accuracy of HT measurements are satisfactory; • The DB correlation underestimates, while the ST correlation gives practically the same Nu numbers for constant HF; • The ST correlation values are less than the numerical simulation result till Re = 14,000 and are above the simulation result only when Re = 17,000. It should also be mentioned that the ST correlation by [23] is valid only to Re = 10,000, but according to our results, it can be used at least to Re = 17,000; • The correlation that contains the friction factor performs better at higher Re numbers. It should be mentioned that the application range for vonK starts from Re = 10,000, and for Gn, from Re = 2300, but vonK performs well below Re = 10,000 as well.
Nanomaterials 2022, 12, x FOR PEER REVIEW 12 of 17 alongside the tube length, the HT is more intensive in the case of constant wall temperature. Moreover, the observed differences in the context of the accuracy of HT measurements are satisfactory; • The DB correlation underestimates, while the ST correlation gives practically the same Nu numbers for constant HF; • The ST correlation values are less than the numerical simulation result till Re = 14,000 and are above the simulation result only when Re = 17,000. It should also be mentioned that the ST correlation by [23] is valid only to Re = 10,000, but according to our results, it can be used at least to Re = 17,000; • The correlation that contains the friction factor performs better at higher Re numbers. It should be mentioned that the application range for vonK starts from Re = 10,000, and for Gn, from Re = 2300, but vonK performs well below Re = 10,000 as well.

Constant HF with Variable Properties
The properties of the fluids vary by temperature, but when the difference between the mean fluid temperature and the tube wall temperature is not high, the differences that the property variations may cause are nearly negligible. Table 7 illustrates the calculation results with CFD simulations using constant and variable properties. By considering the variable properties, the Nu number is generally higher by 1% in the investigated cases.  Table 8 illustrates the calculation results with constant and variable property CFD simulations. By considering the variable properties, the Nu number is generally higher by

Constant HF with Variable Properties
The properties of the fluids vary by temperature, but when the difference between the mean fluid temperature and the tube wall temperature is not high, the differences that the property variations may cause are nearly negligible. Table 7 illustrates the calculation results with CFD simulations using constant and variable properties. By considering the variable properties, the Nu number is generally higher by 1% in the investigated cases.  Table 8 illustrates the calculation results with constant and variable property CFD simulations. By considering the variable properties, the Nu number is generally higher by 2% in the investigated cases. The higher differences are due to the higher differences between the wall and mean fluid temperatures.

HT in Flow Development
Comparing the constant HF and constant wall temperature, one might find only limited differences according to the practical approach, as reported in [29]. Figure 12 reports the relative local Nusselt number variations (Nu/Nu∞) in constant HF for Re = 5000 and Re = 17,000 for 1.5 vol% concentrate NF. These results are in agreement with those shown in Figures 8 and 9, based on [29]. The decrease in the Nu number is instant, and the development of the flow ends around Le = 10·D when the Reynolds number is above 10 4 , but for smaller Reynolds numbers, the flow development occurs at a somewhat greater distance. The authors advise the use of short tubes in the practical engineering application, considering flow development. A more intensive HT would occur than expected. Generally, this is not considered a problem and can be ignored. 2% in the investigated cases. The higher differences are due to the higher differences between the wall and mean fluid temperatures.

HT in Flow Development
Comparing the constant HF and constant wall temperature, one might find only limited differences according to the practical approach, as reported in [29]. Figure 12 reports the relative local Nusselt number variations (Nu/Nu ∞ ) in constant HF for Re = 5000 and Re = 17,000 for 1.5 vol% concentrate NF. These results are in agreement with those shown in Figures 8 and 9, based on [29]. The decrease in the Nu number is instant, and the development of the flow ends around Le = 10·D when the Reynolds number is above 10 4 , but for smaller Reynolds numbers, the flow development occurs at a somewhat greater distance. The authors advise the use of short tubes in the practical engineering application, considering flow development. A more intensive HT would occur than expected. Generally, this is not considered a problem and can be ignored.

Friction Factors
As mentioned in 2.2.1, the friction factor has a significant role in determining the pumping power of the working fluids, and its value also appears in correlations for calculating the HT coefficient. Before easy access to the math software, several formulae were worked out to find the implicit Prandtl-Kármán-Nikuradse correlation. Table 9 shows all friction factor results we made using the formulae in Table 3, and the CFD calculation pressure drops in different cases. Figure 13 shows the percentage differences from the Prandtl-Kármán-Nikuradse friction factor calculation results for Re = 5000 and Re = 17,000. While the seven other correlation results are symmetrical, the CFD simulation results are "one sided". When the Re value is 5000, the friction factors are less, and when the Re value is 17,000, the friction factors are higher than the references calculated by the

Friction Factors
As mentioned in Section 2.2.1, the friction factor has a significant role in determining the pumping power of the working fluids, and its value also appears in correlations for calculating the HT coefficient. Before easy access to the math software, several formulae were worked out to find the implicit Prandtl-Kármán-Nikuradse correlation. Table 9 shows all friction factor results we made using the formulae in Table 3, and the CFD calculation pressure drops in different cases. Figure 13 shows the percentage differences from the Prandtl-Kármán-Nikuradse friction factor calculation results for Re = 5000 and Re = 17,000. While the seven other correlation results are symmetrical, the CFD simulation results are "one sided". When the Re value is 5000, the friction factors are less, and when the Re value is 17,000, the friction factors are higher than the references calculated by the Prandtl-Kármán-Nikuradse correlation. Our observation led us to not analyse these results in this study and leave it for future studies of the scientific community. The same results for No. 9-16 are due to the non-variable fluid properties. Prandtl-Kármán-Nikuradse correlation. Our observation led us to not analyse these results in this study and leave it for future studies of the scientific community. The same results for No. 9-16 are due to the non-variable fluid properties.   The friction factor is more general than the pressure drop; Figure 14 shows the pressure drop of HNs obtained from simulation for non-variable properties. The results present how the pressure drop rises with higher Re numbers. If one substitutes the ξ(Re) function into the pressure drop calculation formula, the result is as follows: where the ξ(Re) is considered as C·Re m , according to our results, and the trendlines exponent 1.767, one obtains m = 1.767 − 2 = −0.233, which is close to the Blasius correlation exponent (−0.25), as shown in Table 3. The friction factor is more general than the pressure drop; Figure 14 shows the pressure drop of HNs obtained from simulation for non-variable properties. The results present how the pressure drop rises with higher Re numbers. If one substitutes the ( ) function into the pressure drop calculation formula, the result is as follows: where the ( ) is considered as • , according to our results, and the trendlines exponent 1.767, one obtains m = 1.767 − 2 = −0.233, which is close to the Blasius correlation exponent (−0.25), as shown in Table 3.
The pressure drop relatively increases by 56% and 137% at 0.5 and 1.5 vol%, respectively, because of the rise in viscosity.

Conclusions
In this research, the HT performance of SiO2-P25 nanofluids passing through a circular tube was investigated by CFD with constant HF and constant wall temperature. During the analysis, a comprehensive treatment of turbulent HT was offered by mentioning different known but often ignored details of the phenomenon. The following conclusions were drawn:  The pressure drop relatively increases by 56% and 137% at 0.5 and 1.5 vol%, respectively, because of the rise in viscosity.

Conclusions
In this research, the HT performance of SiO 2 -P25 nanofluids passing through a circular tube was investigated by CFD with constant HF and constant wall temperature. During the analysis, a comprehensive treatment of turbulent HT was offered by mentioning different known but often ignored details of the phenomenon. The following conclusions were drawn:

•
Most of the simplifications in the context of the turbulent flow are acceptable, but one has to refer to the borders of the simplifications; • When the real conditions do not meet the simplified circumstances, one should refer to the complex phenomenon; • In the case of moderate temperature differences, the constant properties give similar results to the temperature-dependent case with lower computational cost; • The Nusselt number and pressure drop increase with increasing concentrations of hybrid nanofluids and flow rate, in agreement with many experimental and theoretical studies. • Compared with previous research, the simulated results are acceptable.