The Influence of Nusselt Correlation on Exergy Efficiency of a Plate Heat Exchanger Operating with TiO 2 :SiO 2 /EG:DI Hybrid Nanofluid

: This paper studies how the correlation with the Nusselt number affects the final result of the efficiency, ε , and exergy efficiency, η ex , of a chevron-type gasketed plate heat exchanger, which is installed in a typical small solar installation dedicated to single-family housing; the solar fluid is a TiO 2 :SiO 2 /EG:DI hybrid nanofluid with concentrations from 0% to 1.5% vol. The experimental model assumes constant flow of the solar fluid and varies on the domestic hot water side—from 3 lpm to 6 lpm. The inlet temperatures are 30 ◦ C and 60 ◦ C on the cold and hot sides of the heat exchanger, respectively. Of the six analysed correlations that showed similar trends, it is concluded that for the assumed flow conditions, geometry, and chevron angle of the plate heat exchanger, one model is the most accurate. The largest difference between the η ex values for a given concentration is 3.4%, so the exergy efficiency is not affected by the chosen Nusselt model by very much. However, the choice of correlation with the Nusselt number significantly affects the efficiency, ε ; the difference between the values obtained within a given concentration is more than 40% and depends on the Reynolds number and flow. Most research discusses the scenario with the nanofluid as a coolant. This paper considers the opposite situation in which the solar fluid is a hotter working medium that transfers heat to domestic hot water installation.


Introduction
Nowadays, much attention is paid to environmental aspects, sustainable development, and energy efficiency [1].Therefore, the role of heat exchangers, especially in HVAC (heating, ventilation, and air conditioning) installations, is becoming more important.These devices make it possible to improve the efficiency of systems by using and reusing, for example, waste heat.Additionally, heat exchange between two fluids with different thermodynamic parameters takes place in a non-contact manner.Intensive research is also carried out on working fluids whose thermodynamic potential and heat capacity allow for even more effective heat exchange or cold recovery [2].The new generation of working fluids cooperating with heat exchangers undoubtedly includes mono-and hybrid nanofluids of various concentrations.The most famous and tested heat exchangers include double-tube heat exchangers [3], plate heat exchangers (PHEs) [4], and shell and tube heat exchangers [5].By entering the two terms 'nanofluids' + 'heat exchangers' in the search engine on the ScienceDirect platform, a total of 12,859 publications of various types are obtained; the first dates back to 1998 [6].Among the types of publications available in Scopus, the most common are research articles, 78.6%; followed by review articles, 10%; short communications, 4.5%; and book chapters, 3.3%.The remainder constitutes approximately 3%, including conference abstracts, encyclopaedia, case reports, editorials, mini reviews, correspondence, news, and other types (see Figure 1).Furthermore, up to 63% of the publications concern PHEs.And, this type of exchanger is also the subject of this paper.
Inventions 2024, 9, x FOR PEER REVIEW 2 of 16 reviews, correspondence, news, and other types (see Figure 1).Furthermore, up to 63% of the publications concern PHEs.And, this type of exchanger is also the subject of this paper.The most common analysis in the literature concerns the efficiency of PHEs operating with nanofluids on the secondary side of the installation [2][3][4][5].This means that the heatreceiving factor is a suspension of various concentrations, produced on the basis of deionised water [7], and mixtures of water with glycols [8] or oils [9].The methods of producing and stabilising nanofluids are described in detail in [10], and the basic limitations and requirements include, among others, economic considerations or problems with increasing pumping power in [11][12][13].The use of nanofluids is extremely effective due to their proven superior thermal and rheological properties [14].The thermal conductivity coefficient can be from several (0.25% CaCO3(50:50%)/DI [15] or 0.3% Fe3O4-CNT/water [16]) to even several dozen percents higher (1% TiO2 + MWCNT/EG [17] and 1.5% Al2O3 + MWCNT/thermal oil [9]) than in the case of the base liquid.Nanofluids are also eagerly tested in closed-loop solar systems [18], in which solar radiation energy is converted into thermal [19] or electrical energy [20].The efficiency of exchanger installations can be effectively determined using the exergy method, as in, for example, [21][22][23].Iranmanesh et al. in [24] showed that graphene nanoplatelets (water-based) nanofluids with a concentration of 0.1% wt. in a solar installation with a PHE will increase the exergy and thermal efficiencies by approximately 40% and 20%, respectively, compared to distilled water; the flow of working fluid was (1.5 dm 3 /min).
The algorithm to determine the exergy efficiency of the plate heat exchanger requires, among others, assuming the appropriate correlation to the Nusselt number, Nu.It is directly related to other similarity numbers such as the Reynolds number, Re, or Prandtl number, Pr, with Nu = f(Re, Pr), and allows us to calculate the heat transfer coefficient, h (see Equation ( 1)) [25].The Nu number is the quotient of the thermal energy exchanged by convection with the conduction energy, which is generally written as Equation (2).Table 1 presents the universal, most known and proven functional relationships for the Nu number, which are used by researchers around the world to analyse laminar [26] and turbulent [27] flows, including nanofluids, but without taking into account their concentration and the type of nanoadditives.Some correlations take into account the influence of The most common analysis in the literature concerns the efficiency of PHEs operating with nanofluids on the secondary side of the installation [2][3][4][5].This means that the heat-receiving factor is a suspension of various concentrations, produced on the basis of deionised water [7], and mixtures of water with glycols [8] or oils [9].The methods of producing and stabilising nanofluids are described in detail in [10], and the basic limitations and requirements include, among others, economic considerations or problems with increasing pumping power in [11][12][13].The use of nanofluids is extremely effective due to their proven superior thermal and rheological properties [14].The thermal conductivity coefficient can be from several (0.25% CaCO 3 (50:50%)/DI [15] or 0.3% Fe 3 O 4 -CNT/water [16]) to even several dozen percents higher (1% TiO 2 + MWCNT/EG [17] and 1.5% Al 2 O 3 + MWCNT/thermal oil [9]) than in the case of the base liquid.Nanofluids are also eagerly tested in closed-loop solar systems [18], in which solar radiation energy is converted into thermal [19] or electrical energy [20].The efficiency of exchanger installations can be effectively determined using the exergy method, as in, for example, [21][22][23].Iranmanesh et al. in [24] showed that graphene nanoplatelets (water-based) nanofluids with a concentration of 0.1% wt. in a solar installation with a PHE will increase the exergy and thermal efficiencies by approximately 40% and 20%, respectively, compared to distilled water; the flow of working fluid was (1.5 dm 3 /min).
The algorithm to determine the exergy efficiency of the plate heat exchanger requires, among others, assuming the appropriate correlation to the Nusselt number, Nu.It is directly related to other similarity numbers such as the Reynolds number, Re, or Prandtl number, Pr, with Nu = f(Re, Pr), and allows us to calculate the heat transfer coefficient, h (see Equation ( 1)) [25].The Nu number is the quotient of the thermal energy exchanged by convection with the conduction energy, which is generally written as Equation (2).Table 1 presents the universal, most known and proven functional relationships for the Nu number, which are used by researchers around the world to analyse laminar [26] and turbulent [27] flows, including nanofluids, but without taking into account their concentration and the type of nanoadditives.Some correlations take into account the influence of liquid viscosity [27][28][29], while others do not [30,31].The Nusselt number should be calculated separately for each side of the plate heat exchanger, i.e., for the working medium with different flows and temperatures, in the hot and cold channels, respectively.In [27], Akturek et al. performed a new Nusselt correlation for the gasket type of PHE with a chevron angle of 30 • and a heat transfer/extended surface area of 0.14 m 2 .The counterflow water is maintained in the range of 0.57-6.6m 3 /h, on both the cold and hot sides, and the respective temperatures are T c = 9-25 • C and T h = 53-90 • C.An uncertainty analysis of the basic experimental parameters are as follows: for the Reynolds number, u(Re) = 4.16; for the Nusselt number, u(Nu) = 1.83; for the Prandtl number, u(Pr) = 2.19; for heat transfer, u(Q) = 0.40; for the overall heat transfer coefficient, u(U) = 0.41; for the mean logarithmic temperature, u(LMTD) = 0.01; and for the friction factor, u(f ) = 4.18.Proposed by Akturek et al. [27], the correlation for Nu is compared with other similar patterns, and the best fit is noted for Kumar et al. [29].
A very old correlation from 1936 given by Sieder and Tate in [28] originally for tubulartype water heat exchangers has also been used in nanofluids, e.g., by Alklaibi et al. [4], who investigated the cooling properties of hybrid nanofluids when flowing through the PHE.
For new and more accurate models, the Kumar et al. [29] correlation for the Nusselt number is recommended by [32,33].
The correlation given by Okada et al. [31] gives Nu values with even more than a 30% difference from other studies [27,34].
where k is the conductive heat transfer coefficient of the fluid, Nu is the Nusselt number describing the phenomenon of dynamic similarity, and D refers to the hydrodynamic channel diameter.
where E conv , E cond are thermal energy exchanged by convection and conduction, respectively.It should also be noted that for precisely defined experimental conditions of flow and temperature, correlations with the Nu number are also suggested, which are true for a specific nanofluid and concentration range.Models like this take into account very precisely defined thermodynamic parameters of the working medium, such as the conductivity of nanoparticles and the base liquid, as well as their density and specific heat (see Equation (3)).Yashawantha et al.,in [35], developed a new correlation (see Equation ( 4)), which was checked for nanofluid temperatures, with T nf,in in the range from −5 up to +20 • C. On the other side of the heat exchanger, the temperature of the working medium was T h,in = 60 • C. The experiment used a 0.2-2% wt.Al 2 O 3 nanofluid produced on the basis of a mixture of deionised water and ethylene glycol in a weight ratio of 65:35.Ultimately, a good agreement with the experimental data was obtained.The presence of nanoadditives resulted in an increase in heat transfer rate by a maximum of 5.07% (for a nanofluid concentration of 2%) compared to the base liquid; the overall heat transfer coefficient, U, increased by a maximum of 14.99%; the effectiveness, e, increased by 5.18%; and the power requirement for pumping, dp, increased by 13.74% (for a concentration and a temperature of nanofluid at 2% and T nf,in = −10 • C).The performance index was above unity, hence the conclusion that nanofluids are useful for practical and commercial use.The coefficient of determination was given as R 2 = 0.98.
where φ is the volume or weight concentration of the suspension; k is the thermal conductivity; c is the specific heat; and indexes nf and bf refer to the nanoparticles and the base liquid, respectively.
where, a, b, m, and n are the constants received from the regression analyses and depend on the nanofluid inlet temperature, φ is concentration, and the range of validated Reynolds number is Re nf,c = 100-200 (2-4 lpm).
In [36], local Nusselt number correlations for the heat transfer of nanofluids using genetic programming were proposed.It turns out that the variables that determine the Nusselt number are primarily the Re and Pr numbers, and the presence of nanoparticles affects the flow inertia forces and the thermal diffusivity.Many research and numerical works on nanofluids include Nusselt number correlations, which were validated for a specific flow range, i.e., the Re number.Most publications are concerned with the turbulent flow of nanofluids, e.g., in a circular tube [37][38][39].Relatively few works have concerns about the low Reynolds number range, e.g., during the flow of nanofluids in micro heat sinks [8], through the PHE [35], or just in the pipe [40].Lai et al. in [40] investigated an Al 2 O 3 -water nanofluid during laminar flow (Re < 270) in a millimetre-sized stainless steel tube.They showed an improvement in Nusselt number of about 8%.This research was conducted with a nanoparticle volume fraction of 0-1% and a size of 20 nm.Table 2 lists exemplary and more well-known correlations of the Nusselt number for nanofluids and the validation range.
Optimising the operation of exchanger systems involves selecting the output parameters so that their operation is highly efficient and energy-saving [41,42].This often requires checking and applying specific solutions provided in the literature, such as correlations with the Nusselt number.Indirectly, it enables the calculation of many parameters related to heat transfer and the effective design of a given process.Therefore, this work attempts to analyse the impact of the adopted correlation on the Nusselt number on the final result of the efficiency, ε, and exergy efficiency, η ex , of the PHE working in a typical solar installation.The atmospheric conditions of the experimental system also deserve attention.Most research studies that calculate the efficiency of the PHE assume that the nanofluid works on the secondary side of the plate heat exchanger as a coolant.This article assumes the opposite situation.The TiO 2 :SiO 2 hybrid nanofluid with a concentration in the range of 0.5-1.5% is a solar fluid that circulates in the hot water installation on the primary side (that is, outside) of the PHE.The nanofluid is produced on the basis of ethylene glycol.This scenario occurs for installations located in countries where external air temperatures are low enough to cause the solar fluid to freeze, e.g., in central European countries.In this work, it is assumed that the heat from the solar installation will be transferred to the internal domestic water installation with an initial temperature of T c,in = 30 • C; the initial temperature of the working medium on the primary side is T nf,in = 60 • C.This paper also validates the research methodology for other publications.

Thermal Performance Analysis
Due to that this paper being a numerical study, the operation of the heat exchanger was modelled in MATLAB.This experimental modelling concerns determining the efficiency of heat exchange with countercurrent flow of heating and heated liquids by a plate heat exchanger installed in a solar installation for preparing domestic hot water for a singlefamily residential building.It is assumed that the device is made of stainless steel with the dimensions given in Table 3 and Figure 2.
The technical parameters of the plate heat exchanger and its area are given in Table 3, and they are initial conditions.Heat exchangers of similar size and structure are commonly used in industry.The purpose of this work is not strictly to determine the size of the heat exchanger, but of course, it is possible.On the one hand, the total heat transfer surface area (A) can be estimated based on the known relationship, assuming the logarithmic mean temperature (∆T LM ) and the overall required heat transfer coefficient (U); A = Q/(U ∆T LM ).On the other hand, the plate surface area increases with corrugations on plates, which further leads to high heat transfer.The ratio of effective area to projected area of the plate is defined as the surface enlargement factor (Equations ( 5)-( 9)), and for commercial plates, it generally ranges from 1.15 to 1.25.
The enlargement factor of corrugation can also be determined as follows [43,44]: where γ is a channel profile aspect ratio.The solar fluid in this case is a TiO2 + SiO2 hybrid nanofluid with a concentration in the range of 0.5-1.5% vol.; the base liquid is a mixture of deionised water and bioethylene glycol in the DI:EG ratio 60:40.The thermodynamic properties of the working medium are presented in Table 4; the values come from our own research and are very close to that presented in the literature [45].Moreover, the initial conditions are also given in Figure 2. On the solar liquid side, under the operating conditions of the installation, the flow should be constant and the value recommended by [46] is 3 dm 3 /min; while on the hot water side, it is the result of the current hot water consumption and ranges from 3 to 6 dm 3 /min (which corresponds to the water outflow from the draw valve).The heat transfer was performed by modelling within MATLAB.
A P = L P w P (6) where L and w are the plate length and plate width, respectively; L p and w p are the projected length and projected width between the inlet and outlet ports; and d p is the port diameter.
The enlargement factor of corrugation can also be determined as follows [43,44]: where γ is a channel profile aspect ratio.The solar fluid in this case is a TiO 2 + SiO 2 hybrid nanofluid with a concentration in the range of 0.5-1.5% vol.; the base liquid is a mixture of deionised water and bioethylene glycol in the DI:EG ratio 60:40.The thermodynamic properties of the working medium are presented in Table 4; the values come from our own research and are very close to that presented in the literature [45].Moreover, the initial conditions are also given in Figure 2. On the solar liquid side, under the operating conditions of the installation, the flow should be constant and the value recommended by [46] is 3 dm 3 /min; while on the hot water side, it is the result of the current hot water consumption and ranges from 3 to 6 dm 3 /min (which corresponds to the water outflow from the draw valve).The heat transfer was performed by modelling within MATLAB.

Methodology
The exergy efficiency [47] of the PHE is defined as the ratio of useful exergy output to total exergy input and can be written as follows: where T a and T w,ex are the ambient and external surface temperatures of the PHE, respectively; . S gen is the total entropy generation of the fluid (see Equation ( 12)) that is a sum of thermal .
S gen,th and friction .

S gen, f r entropy; and
.
Q h is the heat transfer rate on the hot fluid side. .
Q av is the actual, average heat flux; T in , T out are the temperatures of the working fluid during input and output to and from the heat exchanger; L v is the height of the PHE; ρ is the density; and f is the friction factor.
Then, based on Equations ( 13) and ( 14), the Reynolds and Prandtl numbers were calculated separately for fluids in the cold and hot channels of the studied PHE.

Re
where G is the mass velocity; D h is the hydraulic diameter that includes a mean channel gap b and enlargement factor ϕ, with D h = 2b ϕ ; and µ is the dynamic viscosity of the fluid.Now, in the next step, the Nusselt similarity number can be calculated on the basis of a few correlations, listed in Table 1 [27][28][29][30][31]; additionally, calculations were performed for the model proposed by Xiu and Liu in [39].
Moreover, for the countercurrent flow, the PHE effectiveness can be calculated based on the following: where NTU is the number of transfer units, C is the heat capacity, C min = min(C h , C c ), and where U is the heat transfer coefficient and A is the area of heat transfer.The heat capacity rate ratio C min /C max in the solar system is assumed to be a constant value.On the hot water circuit, the rate is variable and depends on the flow; it takes values from 0.97 for V = 0.05 L/s to 0.486 for V = 0.1 L/s.

Results and Discussion
In the hot channel, due to the constant flow of the working medium at 3 dm 3 /min, the value of the Re number is slightly variable and depends primarily on the nanofluid concentration in the installation and its thermodynamic parameters and ranges from 430.25 to 328.67 for the concentration of nanofluid TiO 2 :SiO 2 from 0% to 1.5% vol., respectively.In the cold channel, on the side of the domestic hot water installation, the Re number varies from 719 for a flow of 3 dm 3 /min to 1438 for a flow of 6 dm 3 /min.
Figure 3 shows how the Reynolds number changes with the concentration of the hybrid nanofluid and is inversely proportional to it.Next, for ~329 < Re nf < ~430, the Nu can be expressed for models [27][28][29][30][31]39].Furthermore, the dynamic viscosity of the nanofluid changes with concentration and temperature so it is safe and accurate to give a dependency in the form of Equation (13).
Figure 4a-d show how the Nusselt number changes with the Re and nanofluid concentration.According to data published in the literature, nanofluids consisting of two types of nanoparticles are characterised by a higher Nu than mono-nanofluids.An important observation is that ethylene glycol-based nanofluids also increase the Nu number from 10 to 32% [48].Moreover, the Nu increases proportionally to the concentration of the suspension.
The range of Re numbers on the primary side of the installation is Re = 329-430 and refers to laminar flow.Figure 4 shows how the correlations available in the literature can affect the Nu number and, consequently, the result of the final analysis of the energy efficiency of a given system (further, Figures 5 and 6).The models defined by Akturk et al. [27], Sieder and Tate [28], Kumar et al. [29], Focke et al. [30], Okada et al. [31], and Xiu and Liu [39] were used for the analysis.
Figure 3 shows how the Reynolds number changes with the concentration of the hybrid nanofluid and is inversely proportional to it.Next, for ~329 < Renf < ~430, the Nu can be expressed for models [27][28][29][30][31]39].Furthermore, the dynamic viscosity of the nanofluid changes with concentration and temperature so it is safe and accurate to give a dependency in the form of Equation ( 13). Figure 4a-d show how the Nusselt number changes with the Re and nanofluid concentration.According to data published in the literature, nanofluids consisting of two types of nanoparticles are characterised by a higher Nu than mono-nanofluids.An important observation is that ethylene glycol-based nanofluids also increase the Nu number from 10 to 32% [48].Moreover, the Nu increases proportionally to the concentration of the suspension.The range of Re numbers on the primary side of the installation is Re = 329-430 and refers to laminar flow.Figure 4 shows how the correlations available in the literature can affect the Nu number and, consequently, the result of the final analysis of the energy efficiency of a given system (further, Figures 5 and 6).The models defined by Akturk et al.The exergy efficiency, ηex, of the PHE in each case increases with concentration and the nanofluid flow through the PHE.However, if identical flow is assumed on the primary and secondary sides of the PHE, in this case 3 lpm, some of the correlations, e.g., Sieder and Tate [28], indicate very low values of ηex (see Figure 5a).This means that the usefulness of the process is relatively low.When DHW is distributed and the flow through the installation is increased to 6 lpm, a visible increase in exergy efficiency is also observed, as shown in Figure 5b.The largest difference between the exergy efficiency, ηex, obtained within a given concentration is 3.4% for Re = 329; the difference decreases with the flow of nanofluid through the solar installation to 2.4% for Re = 430.Therefore, it is concluded that the adoption of the correlation model on the Nusselt number does not significantly affect the final result of the PHE exergy efficiency.The efficiency, ε, of the PHE is enhanced with the concentration of nanofluid in the solar installation and the flow on the primary side of the heat exchanger.The largest difference between the ε values obtained within a given concentration is 49% for Re = 329 but decreases with the flow of nanofluid in the solar installation to 42% for Re = 430.It is  The exergy efficiency, ηex, of the PHE in each case increases with concentration and the nanofluid flow through the PHE.However, if identical flow is assumed on the primary and secondary sides of the PHE, in this case 3 lpm, some of the correlations, e.g., Sieder and Tate [28], indicate very low values of ηex (see Figure 5a).This means that the usefulness of the process is relatively low.When DHW is distributed and the flow through the installation is increased to 6 lpm, a visible increase in exergy efficiency is also observed, as shown in Figure 5b.The largest difference between the exergy efficiency, ηex, obtained within a given concentration is 3.4% for Re = 329; the difference decreases with the flow of nanofluid through the solar installation to 2.4% for Re = 430.Therefore, it is concluded that the adoption of the correlation model on the Nusselt number does not significantly affect the final result of the PHE exergy efficiency.The efficiency, ε, of the PHE is enhanced with the concentration of nanofluid in the solar installation and the flow on the primary side of the heat exchanger.The largest difference between the ε values obtained within a given concentration is 49% for Re = 329 but decreases with the flow of nanofluid in the solar installation to 42% for Re = 430.It is V n f = 6 dm 3 /min [28][29][30][31]39].
As it turns out, the models of Akturk et al. [27] and Okada et al. [31] are characterised by the highest accuracy; the standard uncertainty of type A is below u A < 2% [27].The average difference between the obtained results is below 10% for Re < 1000.However, the model [27] is validated for Re > 450 and [31] for Re > 400.For lower values, Re < 430, it is suggested to use the model of Kumar et al. [29], which takes average values.In terms of flows resulting from this experiment, it differs from [28] and [31] on average by 80% and 48%, respectively, and its validation has already been carried out for Re > 10.It should be noted that the correlation proposed by Focke et al. [24] gives results up to three times higher for the laminar range, 150 < Re < 600; the discrepancies increase as the turbulence brews.The model of Xiu and Liu [39], for the 1.5% nanofluid, differs on average by 55% from [28], by 58% from [29], and by 87% from [30]; validation of the model [39] was performed for CuO nanofluids in the concentration range up to 2%.
The correlations by Okada [28] and Focke et al. [30] also do not take into account the influence of the liquid viscosity, which may be important in the case of nanofluids.Nanofluids are characterised by even several times higher viscosity than that of water (see Table 4).The chevron angle of the PHE [43] should also be taken into account.However, not all researchers include it as a validation parameter.
Finally, for low flows and the full range of Re = 329-430 covered by this work, it is assumed that the model of Kumar et al. [29] is proven and safe when modelling heat transfer in PHEs, as also suggested by researchers Unverdi and Islamoglu in [49].Furthermore, the geometry of the plate and the chevron angle required in this work (which is β = 30 • ) coincides with the model; in [29], the influence of a viscosity is also taken into account.
The following figures show how the adopted correlation to the Nusselt number affects the final results of the exergy efficiency (see Figure 5) and the effectiveness (see Figure 6) of the PHE.The results were compared for five models, as above.
The exergy efficiency, η ex , of the PHE in each case increases with concentration and the nanofluid flow through the PHE.However, if identical flow is assumed on the primary and secondary sides of the PHE, in this case 3 lpm, some of the correlations, e.g., Sieder and Tate [28], indicate very low values of η ex (see Figure 5a).This means that the usefulness of the process is relatively low.When DHW is distributed and the flow through the installation is increased to 6 lpm, a visible increase in exergy efficiency is also observed, as shown in Figure 5b.The largest difference between the exergy efficiency, η ex , obtained within a given concentration is 3.4% for Re = 329; the difference decreases with the flow of nanofluid through the solar installation to 2.4% for Re = 430.Therefore, it is concluded that the adoption of the correlation model on the Nusselt number does not significantly affect the final result of the PHE exergy efficiency.
The efficiency, ε, of the PHE is enhanced with the concentration of nanofluid in the solar installation and the flow on the primary side of the heat exchanger.The largest difference between the ε values obtained within a given concentration is 49% for Re = 329 but decreases with the flow of nanofluid in the solar installation to 42% for Re = 430.It is therefore concluded that the adoption of the correlation model based on the Nusselt number in this case significantly affects the efficiency of the PHE; the coefficient of determination resulting from the analysis of variance is R 2 = 1.
Moreover, based on the analysis of variance, the strength of the relationship between the exergy efficiency, η ex , and efficiency of the PHE, ε, and the flow expressed by the Re number for the exemplary correlation of Kumar et al. [29] are determined.
Based on the statistical parameters calculated in Table 5, and the coefficient of determination of R 2 > 0.99, it is concluded that the exergy efficiency, η ex , and the efficiency of the PHE, ε, strongly depend on the flow conditions in the PHE.However, residuals sum of squares assume much lower values for ε than in the case of η ex , which may indicate a better fit of the model.In each case, the sum of squares of the residuals is much lower than the sum of squares of the regression model inputs.As a supplement and validation of the adopted methodology, the friction factor is also computed.As the flow is laminar, the basic Hagen-Poiseuille relationship, f = 64/Re, was used.A graph of the relationship between the friction factor and the flow is shown below.The Hagen-Poiseuille correlation was compared with the one proposed by Alklaibi et al. [4]-see Equation (14).As is shown in Figure 7, there is a 8% difference in the obtained values.
Adj. R-Square 0.99555 0.99532 0.99544 0.98812 0.98826 As a supplement and validation of the adopted methodology, the frictio also computed.As the flow is laminar, the basic Hagen-Poiseuille relationship was used.A graph of the relationship between the friction factor and the flow below.The Hagen-Poiseuille correlation was compared with the one proposed et al. [4]-see Equation (14).As is shown in Figure 7

Conclusions
Plate heat exchangers are modern and effective heat transfer devices that ca heat recovery and energy efficiency.To optimise the operation of such systems factors need to be selected so that the operation is very efficient and energy-sav fore, this study analyses how the correlation with the Nusselt number affec

Conclusions
Plate heat exchangers are modern and effective heat transfer devices that can increase heat recovery and energy efficiency.To optimise the operation of such systems, the input factors need to be selected so that the operation is very efficient and energy-saving.Therefore, this study analyses how the correlation with the Nusselt number affects the final result of the efficiency, e, and the exergy efficiency, η ex , of a chevron-type gasketed plate heat exchanger which was installed in a typical small solar installation dedicated to single-family housing.Six Nusselt correlations were considered.The heat exchanger with six plates and with a solar fluid as the TiO 2 :SiO 2 /EG:DI hybrid nanofluid on one side and domestic water on the other side was investigated.The experimental model assumes constant flow of the solar fluid and varies on the secondary side of the PHE-between 3 lpm and 6 lpm-and the inlet temperatures are 30 • C and 60 • C on the cold and hot sides.Most research studies that calculate the efficiency of the plate PHE assume that the nanofluid works on the secondary side of the plate heat exchanger as a coolant.This paper assumes the opposite situation.In the hot channel, due to the constant flow of the working medium, the value of the Re number is slightly variable and depends primarily on the concentration of nanofluids in the installation and its thermodynamic parameters and ranges from 430.25 to 328.67 for concentrations of TiO 2 :SiO 2 nanofluid at 0-1.5% vol., respectively.Of the six analysed correlations that showed similar trends, it was concluded that for the assumed flow conditions, geometry, and chevron angle of the PHE, the model of Kumar et al. is the most accurate [29].Furthermore, the adoption of the correlation model in the Nusselt number does not significantly affect the final result of the PHE exergy efficiency, which, according to the correlation [29], is approximately 98% for a 1.5% nanofluid.The largest difference between the η ex values for a given concentration is 3.4% for Re = 329 and decreases with the flow of nanofluids in the solar installation to 2.4% for Re = 430.However, the choice of correlation on the Nusselt number significantly affects the result of the PHE efficiency result, ε; the difference between the values obtained within a given concentration is 49% for Re = 329 and decreases with the flow of nanofluid in the solar installation to 42% for Re = 430.
Funding: This research received no external funding.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The author declares no conflicts of interest.

Figure 1 .
Figure 1.The number of publications registered in the Scopus database for the entries 'nanofluids' + 'heat exchangers' from 1999 to 22 September 2023.

Figure 1 .
Figure 1.The number of publications registered in the Scopus database for the entries 'nanofluids' + 'heat exchangers' from 1999 to 22 September 2023.

Figure 2 .
Figure 2. (a) The hydraulic scheme of the solar system including the plate heat exchanger and operating conditions; (b) initial conditions of the solar system and basic dimensions of a chevron-type PHE with its cross-section dimensions normal to the direction of troughs.

Figure 2 .
Figure 2. (a) The hydraulic scheme of the solar system including the plate heat exchanger and operating conditions; (b) initial conditions of the solar system and basic dimensions of a chevron-type PHE with its cross-section dimensions normal to the direction of troughs.φ = E f f ective area Projected area = A E A P (5)

Figure
Figure4a-dshow how the Nusselt number changes with the Re and nanofluid concentration.According to data published in the literature, nanofluids consisting of two types of nanoparticles are characterised by a higher Nu than mono-nanofluids.An important observation is that ethylene glycol-based nanofluids also increase the Nu number from 10 to 32%[48].Moreover, the Nu increases proportionally to the concentration of the suspension.

Figure 5 .
Figure 5. Exergy efficiency of PHE for TiO 2 + SiO 2 /DI:EG hybrid nanofluid versus Reynolds number for varying Nusselt models and volumetric flow of (a) .V n f = 3 dm 3 /min and

Figure 6 .
Figure 6.Effectiveness of PHE for TiO 2 + SiO 2 /DI:EG hybrid nanofluid versus Reynolds number for varying Nusselt models and volumetric flow of (a) .V n f = 3 dm 3 /min and

Table 1 .
General Nusselt number correlations for the chevron type of PHE.

Table 3 .
Technical parameters of the studied plate heat exchanger.

Table 4 .
Thermodynamic parameters of working fluids.

Table 4 .
Thermodynamic parameters of working fluids.

Table 5 .
Results of regression analysis examining the influence of Re number on the efficiency, ε, and exergy efficiency, η ex , of the PHE for .V h,c = 3 dm 3 /min.