Effect of Nanofluid Thermophysical Properties on the Performance Prediction of Single-Phase Natural Circulation Loops

: Specifying nanofluids’ thermophysical properties correctly is crucial for correct interpretation of a system’s thermo-hydraulic performance and faster market-uptake of nanofluids. Although, experimental and theoretical studies have been conducted on nanofluids’ thermophysical properties; their order-of-magnitude change is still a matter of debate. This numerical study aims to reveal the sensitivity of single phase natural circulation loops (SPNCL), which are the passive systems widely used in solar thermal and nuclear applications, to thermophysical property inputs by evaluating the effects of measured and predicted nanofluid thermophysical properties on the SPNCL characteristics and performance for the first time. Performance and characteristics of an SPNCL working with water-based-Al 2 O 3 nanofluid (1–3 vol.%) for heating applications is evaluated for different pipe diameters (3–6 mm). The thermal conductivity effect on SPNCL characteristics is found to be limited. However, viscosity affects the SPNCL characteristics significantly for the investigated cases. In this study, Gr m ranges are 1.93 × 10 7 –9.45 × 10 8 for measured thermophysical properties and 1.93 × 10 7 –9.45 × 10 8 for predicted thermophysical properties. Thermo-hydraulic performance is evaluated by dimensionless heat transfer coefficients which is predicted within an error band of ±7% for both the predicted and measured thermophysical properties of the data. A Nu correlation is introduced for the investigated SPNCL model, which is useful for implementing the SPNCL into a thermal system. nuclear cooling systems in which nanofluids’ potential is also comprehensively studied. In this study, the influence of thermophysical property definition on the performance of nanofluid based SPNCmLs has been studied numerically. Results of measured thermophysical properties (MTP) and predicted thermophysical properties (PTP) have been compared by using DIW and DIW based Al 2 O 3 nanofluid (1, 2 and 3 vol.%). Numerical study was conducted for different pipe diameters to broaden the operating range. Several major conclusions are presented as follows:


Introduction
In a general sense, an increase in thermal conductivity by addition of nanoparticles into a base fluid (distilled water (DIW), ethylene glycol (EG), etc.), producing a nanofluid, results in higher heat transfer performance. Enhanced thermal conductivity is the main reason for the improved heat transfer performance, but the change in other thermophysical properties (especially viscosity) may reduce the overall efficiency of the system. Therefore, the correct interpretation of the convective heat transfer and overall efficiency is crucial for nanofluid applications to obtain a realistic conclusion [1,2].
The measurement of the thermophysical properties of nanofluids involves some difficulties due to the complex physical interaction between fluid and nano-particles [3]. The thermal conductivity and viscosity of nanofluids are the most studied properties since they highly influence the performance of the system [4][5][6][7][8][9][10]. Flow type and drag type viscometers are generally used for viscosity measurements and transient techniques, such as transient hot wire (THW) and 3-ω method, and are widely used to measure the thermal conductivity of nanofluids [3]. Measurements by commercial devices that are tuned for ordinary fluids results in less accurate and/or misleading property values for nanofluids [11,12]. For example, Antoniadis et al. [11] published the necessary conditions for accurate transient hot wire measurements of nanofluid thermal conductivity. Non-standardized production and measurement procedures of nanofluids also introduced confusing empirical thermophysical property models, which are vital for the implementation of nanofluids into real-life applications for faster and cheaper design of engineering systems [13]. The accuracy of such models is improved by introducing parameters such as concentration, size and shape of particles, micro and nanoscale physical interactions, etc. [14,15]. Maxwell and Hamilton-Crosser, and Einstein and Brinkman models are widely used for thermal conductivity and viscosity prediction, respectively. The difference between Measured Thermophysical Properties (MTP) and Predicted Thermophysical Properties (PTP) by existing models has been studied widely. Table 1 presents a comparison of the data obtained from predictions and measurements for selected studies. It is seen that different models could give different results compared to each other and measurements for a nanofluid. Table 1. Comparison of nanofluids' thermophysical properties obtained from predictions and measurements for selected studies.
Natural circulation loops (NCLs) are well-known passive components in heat transfer systems. NCLs have been used for many applications such as cooling of nuclear reactors, solar water heaters and electronic chip cooling [48][49][50]. Their reliability, stability and performance are widely studied [51][52][53][54][55]. In nuclear reactors, decay heat and core heat removal are provided by using NCL [50,51,56,57]. NCLs are also getting attention for solar thermal applications. Utilization of NCL in solar water heating system avoids freezing and corrosion problems and nighttime heat losses [58], and when a solar still system is designed as NCL, water vapor is transferred from evaporator to condenser without direct contact [59,60]. Moreover, Michael and Iniyan [61] compared the performance of CuO/DIW nanofluid in a flat plate solar water heater under natural and forced circulations. Although the highest efficiency is observed for forced circulation, they found that the highest efficiency improvement occurs for natural circulation.
In single phase natural circulation loops (SPNCLs), temperature difference causes a density gradient (without phase change) and resulting buoyancy forces are the driving mechanism of SPNCLs. They are getting attention due to their advantages of simplicity, easy control, enhanced safety and reliability. In order to obtain a high-performance SPNCL together with robust characteristics (i.e., highly efficient heat transfer and stable flow) within a desired operational range, the effect of different geometries and working fluids has previously been investigated [48]. The geometry of the loop (pipe diameter, loop shape, inclination angle, heater and cooler orientation, scale and aspect ratio of the loop) is considered as the primary factor that affects the loop performance. Recently, nanofluids have also been used in SPNCLs as a working fluid to increase their performance [62][63][64][65][66], improve the stability [62,66] of the system and reduce the required time to reach a steady state [64].
Investigations on mini SPNCLs (SPNCmL) have been increasing during the last decade. Although there is not any certain criterion for defining an SPNCL as mini, SPNCLs having loop dimensions of tens of centimeters and pipe diameters of a few millimeters are assumed as mini loops. The principle driving mechanisms are the same for the large-scale SPNCLs and SPNCmLs, however, miniaturizing affects system behavior strongly [67]. With the knowledge of enhanced stability by miniaturization of SPNCLs, nanofluids were also used to improve thermo-hydraulic performance and stability of SPNCmLs. Misale et al. [68] studied DIW-based Al2O3 nanofluid in an SPNCmL experimentally for the first time. Although using nanofluid improves the thermal performance of SPNCmL slightly for inclined loop orientations, a significant enhancement of thermo-hydraulic performance of the nanofluid-based SPNCmL was not observed compared to DIW. However, results for the mass flow rate inside the loop were in good agreement with the Vijayan's correlation [69], which is one of the important loop characteristics. Turgut and Doganay [70] introduced the effectiveness factor (ε) for the evaluation of the thermal performance of SPNCmLs. ε increases with an increase in nanofluid concentration and decreases with heater power. Karadeniz et al. [71] developed a steady 3D numerical model for predicting the performance and characteristics of the SPNCmLs and validated the model by comparing it with the results of Turgut and Doganay [70]. The deviation from the experimental results was up to 10% for all cases. They concluded that the proposed numerical model is a powerful tool for comparative parametric studies about the performance and characteristics of SPNCmLs.
Doganay and Turgut [72] introduced and evaluated the effectiveness (ε) of Al2O3 nanofluid based SPNCmL experimentally at different combinations of inclination angles and heat sink temperatures. ε was calculated higher for lower heat sink temperatures and higher inclination angles. Karadeniz et al. [73] analyzed the experimental results of Doganay and Turgut [72] with their previous model [71]. Koca et al. [74] used a DIW based Ag nanofluid in an SPNCmL and compared the results with a DIW based Al2O3 nanofluid. SPNCmL using a lower concentration of Ag nanofluid (0.25 wt.%) had the same thermal stability with DIW as a working fluid. However, the Al2O3 nanofluid had an almost constant ε ratio at different heater powers while for higher concentrations of Ag-nanofluid, ε increased significantly with increasing heater power.
Experimental and numerical attempts to represent the performance and characteristics were important for understanding and improving a single SPNCmLs behavior. However, their in-situ characteristics must be well predicted to match them with other components for designing a system properly. SPNCLs can be used for heating and cooling purposes. Therefore, Nu correlations for the hot and cold ends are needed for correctly placing them into an application. Ho et al. [75] used a minichannel heat sink and heat source to investigate a nanofluid based SPNCmL and proposed Nu correlations for both heating and cooling ends. They showed that an increase in concentration increases the average Nu number for both ends. In addition, at constant hot wall temperature, cold wall temperature rise causes a reduction of Nu in the heating section unlike the increments in the cooling section. Although an increase in Ra number increases the Ress number; at constant Ra, Ress decreases with the addition of nanoparticles. Moreover, they also developed a correlation for Ress as a function of nanofluid concentration and Grm to generalize the characteristics of the system. On the other hand, Mohan et al. [76] used mineral oil based-hexagonal boron nitride nanofluid in SPNCmL. Using nanofluid increased Theater, flow rate and maximum Nu number. They also compared the Ress predicted by Vijayan's correlation [69] and experiments and they reported the reason behind the deviation in Ress as the incapability of existing thermophysical property models for oils and nano-oils.
As mentioned above, thermophysical properties of nanofluids are defined either by measurements or theoretical models in the literature. Correct definition of thermophysical properties of nanofluids is crucial for the realistic evaluation of a heat transfer systems' performance. Duangthongsuk and Wongwises investigated both the effect of thermophysical property models [77] and the determination by either MTP or PTP [78] on the calculation of convective heat transfer of a double-tube counter-flow heat exchanger under turbulent flow conditions, experimentally. It is found that different models gave different results, and various models had no significant effect on the heat transfer coefficient and Nu number. Although, MTP had higher values compared to PTP of thermal conductivity and viscosity, the discrepancy between PTP and MTP for the calculation of Nu was found as a few percentages. Recently, Yıldız et al. [79] numerically compared the MTP and PTP in a natural convection cavity under laminar flow conditions and found that PTP utilization underestimated the heat transfer enhancement.
To the best of the authors' knowledge, there are a few studies in the literature on the effect of the determination of thermophysical properties (MTP or PTP) on the overall system performance but there is no study so far for SPNCLs. An SPNCmL is used in this study, which is one of the candidates for the acceleration of nanofluid's market uptake because of its passive nature. As SPNCLs are used to manage a heat transfer need robustly, their in-situ characteristics must be predicted before assembly. This is important for predicting the overall system performance both at design-range and off-design conditions. Within this scope, the main question of this study is how the definition of nanofluids' viscosity and thermal conductivity affect this natural convection-driven heat transfer system's performance prediction. Performance and characteristics of an SPNCmL working with DIW based-Al2O3 nanofluid (1, 2 and 3 vol.%) are evaluated for different pipe diameters (3-6 mm). Furthermore, the effects of nanofluids' thermophysical property prediction method on SPNCmL characteristics and performance has been investigated for the first time.

Numerical Model
Computer simulations of this study have been done by using a pre-validated numerical methodology [73]. Nanofluid is assumed as a single-phase fluid in the steady 3D numerical model, and validation of the model with experimental results of Doganay and Turgut [72] is also re-reported briefly in Section 3.1. The schematic of geometry, as well as, the parameter table including the dimensions and boundary conditions are given in Figure 1. In the heater part, constant heat flux (corresponding to 10, 30 and 50 W of heat rates for the heater surface area) was applied as the heater power input, whereas in the cooler part, a constant temperature was defined as 20 °C, which is the cooling fluid temperature. A commercial software (ANSYS CFX) was used for building the geometry and solving the numerical model (a steady laminar full buoyancy model with temperature-dependent thermophysical properties). The general form of the governing equations for steady-state 3-D incompressible flow with negligible viscous dissipation can be expressed as [80]: Mass conservation equation: Momentum equation: Energy equation: The thermal diffusion in the streamwise/axial direction is also neglected in the calculations. In the governing equations, h is enthalpy (m 2 /s 2 ) and ⃗ is velocity vector (m/s). S denotes the source of momentum (kg m 2 /s 2 ) and energy (kg m/s 3 ). and k are the density (kg/m 3 ) and thermal conductivity (W/m K) of the fluid, respectively. τ is stress tensor and given as: where is viscosity (kg/m s) and is the Kronecker delta. The full buoyant method is applied in the numerical method as: Here, is the reference density of the working fluid, which is calculated at the average fluid temperature in the system and given in Appendix A for all working fluids and investigated cases.
Since the stability of the SPNCL was improved by miniaturization, the steady numerical model can be used to reduce the computational power (see Section 3.1 for the method validation results). Mesh dependency of the solution was investigated for three different meshes (745,456, 565,681 and 348,087 elements, respectively), all including boundary layer elements to correctly resolve thermal and hydrodynamic boundary layers, in the previous study [73] by comparing temperature difference between hot leg and cold leg (ΔTheater) and maximum temperature (Tmax), which are the main parameters for comparing experimental and numerical results. The difference is negligible, therefore the coarse mesh number (346,876) was used in the analysis to reduce computational power. The details of the mesh at different parts of the model are presented in Figure 1. Heater and cooler sections of the system have finer boundary layer elements and finer meshes with 0.6 mm of grid size compared to the rest of the system (0.9 mm and 0.8 mm of grid size for elbows and other pipes, respectively) to model precisely. A high-resolution scheme is used to solve advection terms. Convergence criteria were set to 10 -7 for all of the governing equations in all of the analyses.

Thermophysical Properties
One of the goals of this study is to understand how the determination of nanofluids' thermophysical properties affects the performance of a representative natural convection heat transfer system. The emphasis is given to the effect of measured and predicted nanofluid thermophysical properties. However, both methods agreed well on some properties (specific heat, density and thermal expansivity). Therefore, only predicted data is used as the input to the numerical simulations for these properties. There is still a discussion on the order-of-magnitude change of the viscosity, and thermal conductivity data, which would lead researchers to misleading results for the prediction of overall system performance. Therefore, both measured and predicted data is used as the input to the numerical simulations for these properties to investigate their effect on overall system performance.
Temperature dependence of thermophysical properties for modeling SPNCmLs is also crucial [73]. Tabular values were used to define the properties for DIW [81] as a function of temperature. However, the existing models to predict nanofluid data do not include temperature effects. Therefore, the temperature dependency of the thermophysical properties was derived by assuming the DIW's temperature dependency as the dominant parameter and using the corresponding equation repetitively for calculating the values at different temperatures. Thermophysical properties of Al2O3 nanoparticles are presented in Table 2 [82]. Table 2. Thermophysical properties of Al2O3 nanoparticles [82].

Thermal Expansion (K −1 )
Thermal Conductivity (W m −1 K −1 ) Al2O3 3970 765 0. 85 25 The effective specific heat (C ) values were calculated by Equation (6) [83]. Specific heat does not change with respect to the operational temperature range as seen in Figure 2. Therefore, it is taken as constant for each nanofluid and DIW for the simulations. For Al2O3-DIW nanofluids at different volumetric concentrations (φ), models in the literature were used to determine the effective density (e) [84] and effective thermal expansion (e) [85] as expressed by Equations (7) and (8), respectively ( Figure 2). Specific heat does not change within the operational temperature range as seen in Figure  2. Therefore, it is taken as constant for each nanofluid and DIW for the simulations.
( , , , , , ) = The discussion on the viscosity and thermal conductivity of nanofluids is the main focus of this study. Well-known models such as the Maxwell Equation (9) [17] and Einstein Equation (10) [18] correlations have been used for predicting thermal conductivity and viscosity, respectively. These models have been preferred for many nanofluid based system simulations previously.
( , , , ) = The measured values reported in [16] were used to prepare the thermal conductivity and viscosity data of MTP. A 3-ω setup and a Brookfield rheometer were used for the measurements of thermal conductivity and viscosity, respectively. Measured and predicted thermal conductivity and viscosity values for working fluids are given in Figure 3. Predicted thermal conductivity data by Maxwell model is slightly higher compared to 3-ω system measurements but they are in good agreement by means of tendency. However, the results of the Einstein viscosity model of nanofluids gives almost concentration independent results compared to measured data. Nevertheless, an increase in temperature closes the gap between measured and predicted viscosity data. Second-order polynomial functions were fitted to all temperature-dependent thermophysical data in order to define them in the numerical model except for the Einstein correlation for which a fourth-order polynomial was used to increase the accuracy. Thermophysical properties for all working fluids are presented in Figures 2 and 3 and corresponding polynomials are given in Appendix B.

Performance Indicators
The density gradient due to the temperature difference between heater and cooler sections is the driving mechanism of SPNCLs. Therefore, Tmax and ΔTheater can be considered as the main characteristics of SPNCmL that are used to compare different operational cases. However, it is possible to observe more information and accurate results on the performance of SPNCLs by calculating mass flow rate and non-dimensional indicators such as ε, Nu, Ress, Grm and Ra. These parameters are defined below. Researchers tend to characterize their systems by using some of these performance indicators but there is not an agreement on which represents the performance more accurately.
Steady-state mass flow rate is a crucial parameter for evaluation of the SPNCLs' performance. Due to the difficulties in the measurement of mass flow rate compared to temperature in real life applications, a theoretical approach that is proposed by Vijayan [69] is mostly used (Equation 11). This correlation is valid for the natural circulation loops having a uniform diameter (ranging between 4-40 mm), under steady-state conditions with an incompressible Boussinesq fluid, and where the local pressure losses can be neglected. The corresponding mass flow rate in the loop is calculated by [69]: Here, ΔT is the average temperature difference and N is effective loss coefficient and can be calculated as [69]: where ξ is local friction loss and it is obtained from [82].
During the calculation of the mass flow rate, Ress can be written as [81]: For the evaluation of a natural convection intensity of the working fluid, the modified Grashoff number (Grm) must be determined. Grm is the ratio of gravitational forces to shear stresses. For an SPNCL with a stable flow, Grm can be determined by [68]: Although Grm is used to understand flow characteristics in natural convection, the working fluid that is forced to move while passing along the heating section carries momentum through the loop. Therefore, the system also exhibits forced convection characteristics. For laminar and turbulent flow, Ress and Grm can be correlated by this theoretical approach on mass flow rate [69].
where B and C depend on the loop geometrical parameters and working fluid's thermophysical properties. Vijayan defined B = 0.1768 and C = 0.5 for fully developed laminar flow [69].
is defined originally for the heat exchangers, which is a ratio of the actual heat transfer to maximum possible heat transfer. It is also used to understand SPNCL heat transfer performance. is defined as: Moreover, Num is also used for the evaluation of the heat transfer performance of the SPNCLs and is calculated by [86]: The overall heat transfer coefficient (U) in Equation (17) was evaluated on the basis of the logarithmic mean temperature difference and it refers to the internal area of the cooler [68]. Overall heat transfer coefficient at the cooler was calculated by using the equations below: Num can be calculated in the same manner. Ra number is also considered as a non-dimensional performance indicator due to being useful for investigation of the buoyancy-driven heat transfer and it is introduced as [87]: where Pr is Prandtl number and defined as [87]:

Results and Discussion
In this section, firstly validation of the numerical model using experimental results of Doganay and Turgut [72] is reported briefly. Moreover, the thermo-hydraulic performance of SPNCmL for PTP and MTP cases, and for different pipe diameters are discussed in terms of performance parameters to reveal the sensitivity of the system to thermophysical property inputs. As SPNCmLs are passive components that are used to manage a heat transfer process robustly, their in-situ characteristics must be predicted before assembly. This is important for predicting the overall system performance both at design-range and off-design conditions. Therefore, the SPNCmL performance is generalized by introducing Nu correlations for both thermophysical property inputs and the applicability and accuracy of the Nu correlations are discussed afterward.
A legend system that is valid for all the results given in the following sections is introduced as the high number of parameters makes it difficult to conduct the results. Table 3 compiles colors for nanofluids, and pipe diameters and heater powers are presented in Table 4. Moreover, filled symbols are used for 50 W, patterned symbols for 30 W and empty symbols for 10 W. Symbol "✕" denotes numerical results of 4.75 mm pipe diameter for the validation study.

Validation
Results of the 3D numerical model for PTP and MTP cases are compared with the experimental results of a well-proven setup [70,72,74] (pipe diameter is 4.75 mm) by means of Tmax, ΔTheater and ε. Considering MTP results, Figure 4 shows that results obtained from the numerical model have a good agreement with the experimental data. ε, Tmax and ΔTheater, varied between ±10%, ±20% and ±20% error bands, respectively. As shown in Figure 4, ε has the same tendency with the experimental results and increases with the increase in concentration and decreases with the increasing heater power for both PTP and MTP. Numerical ε results are getting closer to experimental results at lower powers. Heater power doesn't affect the accuracy of the Tmax results, as it is always over-predicted by the simulation. However, ΔTheater is better predicted as the heater power increases opposite to ε results.  Increase in concentration results in a significant increase of ΔTheater for both MTP and PTP, MTP having a higher slope. As concentration increases, Tmax increases significantly for MTP but for PTP a minor change is observed. An increase in viscosity means a higher resistance to flow, which results in lower mass flow rates. The slower motion of the fluid causes higher heat dissipation to the inner layers of the fluid which results in higher average temperatures in cross-sections. Moreover, for all heater powers, DIW results are in good agreement with the experimental results of Tmax and ΔTheater. But for ε, nanofluids are well-predicted compared to DIW at lower heater powers. To conclude, compared to MTP, PTP results fit well with the experimental study in terms of loop temperatures, which are mainly used to characterize the system. However, for ε, which is a performance indicator, as heater power increases PTP results diverge more from the experimental data compared to that of MTP. Therefore, it is not easy to choose one or the other by means of a validation study.

MTP and PTP Comparison for Different Pipe Diameters
A set of pipe diameters (3, 4, 5 and 6 mm) was selected to extend the operating range with the aim of the generalization of the results for the SPNCmL performance and improve the discussion on MTP and PTP cases. Loop temperatures (Tmax and ΔTheater) for different pipe diameters are reported in Figure 5 together with MTP and PTP comparison with ±15 deviation bands. For MTP, Tmax and ΔTheater increase with the increase in heater power and concentration of nanofluid as shown in Figure  5. Power-dependent groups (i.e., empty, striped and filled) of data are observed for both Tmax and ΔTheater. Moreover, the addition of nanoparticles into DIW increases the Tmax and ΔTheater, which means an increase in SPNCmL implementation potential in thermal systems due to higher temperatures. The same behavior has been reported in the literature [70,71,73,76]. However, Bejjam and Kumar [88] said that nanofluids decrease the temperature difference between hot and cold legs. In another study [89], it was reported that for up to 5% vol. concentration temperature difference at the heater increases with the increase in heater power but at 6% vol. it decreases. By using PTP, Figure 5 shows that the concentration effect on Tmax and ΔTheater is small compared to MTP and concentration-dependent difference decreases with the increase in pipe diameter. This effect can be seen slightly in Figure 4 but is more apparent in Figure 5 due to the increased operational range by changing the pipe diameter. The main difference between MTP and PTP is at viscosity values (Figure 3), namely; viscosity is almost unaffected by the particle content when Einstein correlation was used. Therefore, it is thought that this difference is directly related to the viscosity data for the PTP. As pipe diameter increases, the total mass in the SPNCmL increases and results in a higher heat capacity of the working fluid. For constant heat flux boundary conditions, higher heat capacity decreases for both Tmax and ΔTheater as shown in Figure 5. A smaller ΔTheater closes the viscosity range during operation ( Figure 3) and thus, the difference between the Tmax data of PTP and MTP decreases.
ε totally depends on the temperatures of the SPNCmL (Equation 16). Therefore, similar trends with temperatures are observed for ε, too. As presented in Figure 6, ε decreases with the increase in heater power for different loop diameters for both cases. For higher heater powers, ΔTheater and Tmax have higher values, and as a consequence ε increases. The previous studies [70][71][72][73] in the literature have reported that ε increases with an increase in nanofluid concentration. In this study, it is also found that increments in concentration for all diameters enhance effectiveness as a result of higher temperatures. Contrary to the effect of concentration, a decrease in pipe diameter results in an enhancement in ε. Sahu and Sarkar [90] reported that ε decreases with the increased pipe diameter within the range of 25-30 mm, which is much higher than the pipe diameter range considered in this study. Due to the insignificant concentration dependency of viscosity to temperatures for the PTP case (Figure 3), change in ε with concentration is smaller than MTP. However, as shown in Figure 6, decreasing the pipe diameter to below 3 mm would not enhance ε at lower powers for both MTP and PTP based systems because the limit for ε is unity. To conclude, SPNCmLs with smaller pipe diameters at higher heater power have an improved potential to be used in thermal applications. As seen from loop temperature and effectiveness results, the definition of the thermophysical properties with existing correlations or experiments affects the performance evaluation of the SPNCmL. ε is used to predict overall performance of SPNCmL and it is related to temperature differences. In this case, using MTP is thought to be more accurate than PTP.

Generalization of Results
Characterization parameters and effectiveness are not sufficient to understand the performance of the SPNCmLs and to conclude on the PTP and MTP selection. Therefore, non-dimensional parameters are used to generalize the discussion on SPNCmL performance for both thermophysical property inputs. In this study, widely used dimensionless parameters have been used to evaluate pipe diameter, working fluid, and MTP and PTP cases' effects on thermo-hydraulic performance. Although the flow inside SPNCLs is driven by density difference arising from the temperature change, a representative Reynold's number for the steady flow (Ress) can be calculated as given previously in Equation (13). Mass flow rate values obtained from the numerical study were used to calculate Ress in this study. A correlation developed by Vijayan [69] that is valid for fully developed laminar flows for SPNCLs (not mini) and for a wide pipe diameter range (4-40 mm) is used for comparison. Effective loss coefficient (NG) and modified Grashoff number (Grm) were calculated by Equations (12) and (14) by handling the data obtained from the simulations for this purpose. The results were compared with Vijayan's correlation (Equation 15) to check the validity of the correlation for different pipe diameters of SPNCmL ( Figure 7). As seen in Figure 7, the discrepancy between the simulation results and Vijayan's correlation is up to 50% for the MTP case and up to 45% for the PTP case. As seen in Figure 7a,b the data is grouped by the pipe diameters (i.e., circles, triangles, diamonds and squares) firstly and then heater power groups are distinguished (i.e., empty, striped and filled). Ress increases (because of the mass flow rate increase) with the increase in heater power for both cases and decreases as the concentration increases for MTP cases. However, concentration-dependent distribution of PTP data points is less compared to MTP but PTP data in Figure 7b shows that the change of Ress with concentration is insignificant. The main reason behind this behavior could be the minor change of viscosity with concentration for the PTP case. Bejjam and Kumar [89] (studied numerically and used the mixture rule for both thermal conductivity and viscosity data) have found that Re increases with the heater power but an increase in concentration also increases the Ress contrary to MTP and PTP results obtained in this study. Nayak et al. [62,66] (studied experimentally and used the Hamilton-Crosser model for thermal conductivity and Brinkman for viscosity data), Thomas et al. [91] and Thomas and Sobhan [63] (studied experimentally and used the THW method such as KD2 Pro data for thermal conductivity and unknown for viscosity) have reported that Ress increased with the concentration. Moreover, Ho et al. [75] (studied numerically and used KD2 Pro data for thermal conductivity and Brookfield rheometer data for viscosity) reported that concentration of the nanofluid has a higher impact on Ress.
Besides concentration, Ress increases with pipe diameter for all heater powers in this study. Since the system has a higher total mass for higher pipe diameters, mass flow rate and thus, Ress increases. For higher pipe diameters, the reduction of Ress compared to DIW has increased for MTP. As the heater power increases, Grm increases. In addition to heater power, Grm also increases with the increase in loop diameter as Ress but it also increases with the concentration in contrary to Ress for both MTP and PTP. For the MTP case, Grm range is 1.93 × 10 7 -9.45 × 10 8 and Ress range is 30-355; and for the PTP case, Grm varies in the range of 1.93 × 10 7 -7.97 × 10 8 and Ress range is 34-358. Although, concentration effects are smaller on Grm as well as Ress, the minor dependency of viscosity on concentration for the PTP case (Figure 3) results in a smaller change of Grm compared to MTP. Figure 7 shows a significant difference between the Vijayan's correlation (Equation 10) and the numerical results of this study. The deviation is more for the MTP case. Vijayan's correlation (Equation 10) is valid for fully developed laminar flows for SPNCLs (not mini) and for a wide pipe diameter range (4-40 mm). It is clear that pipe dimensions and heater power values are at a smaller range for this study and they affect the flow characteristics together with the concentration, increasing the deviation from the Vijayan's correlation (Equation 10).
The addition of nanoparticles or an increase in concentration decreases the mass flow rate and Re . Therefore, results shifted below the correlation as shown in Figure 7a,b. In order to correlate this shift for DIW and different nanofluid concentrations, a concentration-dependent Ress model for nanofluid based SPNCmLs has been developed by using the numerical results. The correlation which is valid for the working fluids within a concentration range of φ =0-3% and 10 4 ≤ ≤ 10 7 is given in Equation (24) and the plot is given in Figure 8. This correlation is more accurate (data scattering is reduced to 20% bandwidth, which was previously up to 50%) for both the MTP and PTP cases compared to the Vijayan's correlation. As NCLs are passive components that are used to manage a heat transfer need robustly, their in-situ characteristics must be predicted before assembly. This is important for predicting the overall system performance both at design-range and off-design conditions. A wide operational range description of the characteristics of NCLs can be predicted by the dimensionless heat transfer coefficient, the modified Nusselt number (Num). SPNCmLs (also all NCLs) can be used for heating or cooling purposes. Therefore, a Num can be defined either for the heater or for the cooler section. Most of the previous studies on SPNCLs applied constant temperature boundary condition in the heater section and developed Num correlations as a function of Ra [92] and Ress [93]. Different from the related literature, in this study, a nanofluid based SPNCmL is aimed to be used for thermal applications with the heating purpose. Figure 9 presents the relationship of Num for the cooler section with Ra for constant temperature boundary condition. It is obvious that higher heater powers cause a higher Ra and Num as expected for both PTP and MTP. Num decreases as the nanoparticle concentration increases for constant Ra for both cases. Reduction in Num with concentration could be due to neglecting the Brownian motion and thermophoresis effects in the numerical simulations [94]. Moreover, as the pipe diameter increases Num decreases and Ra increases for all working fluids (for both MTP and PTP cases) and heater powers. Changing working fluid and heater powers result in changes of thermophysical properties and thus Ress, Grm and Pr. Figure 9 shows that evaluating Num with respect to Ra results in groups with regard to pipe diameter and heater power. Therefore, numerical results are used to develop a correlation for Num at the cooler section considering a shape factor ( ) for generalization (Equation 25), consistent with the previous literature [92].
As shown in Figure 10, , = , is applicable for SPNCmLs working with different fluids at various pipe diameters. Results show that the correlation given in Equation (25), has a discrepancy of ±7% for both MTP and PTP cases. Considering the deviation of the numerical study from the experimental results (Figure 4), the overall error of the proposed correlations is not more than 20%, which is acceptable for the design purposes of a system including an SPNCmL.

Conclusions
Correct prediction of the thermophysical properties has great importance for investigating the nanofluid's effect on the performance of the heat transfer systems. Due to restrictions arising from the increased viscosity of nanofluids, passive heat transfer systems have great potential in nanofluid's market uptake. SPNCmLs, as passive systems, are widely used in the solar thermal applications and nuclear cooling systems in which nanofluids' potential is also comprehensively studied. In this study, the influence of thermophysical property definition on the performance of nanofluid based SPNCmLs has been studied numerically. Results of measured thermophysical properties (MTP) and predicted thermophysical properties (PTP) have been compared by using DIW and DIW based Al2O3 nanofluid (1, 2 and 3 vol.%). Numerical study was conducted for different pipe diameters to broaden the operating range. Several major conclusions are presented as follows: Although nanofluids have smaller Num at the cooler section compared to DIW, ε of the overall system is higher for nanofluids. The reason behind the reduction of Num could be the Brownian motion and thermophoresis effects of particles;  The concentration-dependent mass flow rate model was introduced for nanofluid based SPNCmLs, which is valid for working fluids with concentrations between = 0-3% and 10 4 ≤ ≤ 10 7 , and more accurate than Vijayan's correlation;  Num correlation for the cooler section is defined as = ( , ) to generalize the heat transfer performance of SPNCmL;  The difference between predicted thermal conductivity and viscosity of nanofluids by wellknown models and measured data does not affect the heat transfer performance significantly.
To conclude, the effect of the thermal conductivity on SPNCmL characteristics is found to be limited for the investigated cases. However, viscosity affects the SPNCmL characteristics significantly as the viscosity values differ greatly for both the measured and Einstein correlation's data. Thermohydraulic performance is also evaluated by calculating the dimensionless heat transfer coefficient. A Nu correlation is introduced for the investigated SPNCmL model, which is useful for implementing the SPNCmL into a solar thermal system. The methodology proposed in this study can be used to derive Nu correlations for different geometries.

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

3D
Three Appendix A.

Reference Density Values of Working Fluids
Buoyancy calculations in the solution require density difference ( − ) and using results with the exclusion of hydrostatic gradient of the pressure in the momentum equation. is calculated at the average fluid temperature for all investigated cases.
values of the simulations are presented in Table A1.

Polynomials for Thermophysical Properties of Working Fluids
Thermophysical properties of working fluids were applied to the numerical model as temperaturedependent polynomial functions of the data presented in Figures 2 and 3.
Temperature-dependent thermophysical properties of DIW used in the numerical model are presented as follows: = 0.6065 a + a + a , Here, a1 = −1.48445, a2 = 4.12292 and a3 = −1.63866 with the unit of (W/(m K)) for all constants and the reference temperature (Tref) is 298.15 K.