Optimization of the Adsorption/Desorption Contribution from Metal-Organic-Heat-Carrier Nanoparticles in Waste Heat Recovery Applications: R245fa/MIL101 in Organic Rankine Cycles

: The efﬁcient recovery of low temperature waste heat, representing from 25% up to 55% of the energy losses in industrial processes, still remains a challenge and even Organic Rankine Cycles (ORCs) experience a strong efﬁciency decay in such a low temperature operating range ( T < 150 ◦ C). In similar heat transfer processes, several nanoﬂuids have been proposed as a solution for increasing heat transfer efﬁciency, but they produced only moderate enhancements of the heat transfer efﬁciency in comparison with pure ﬂuids. This paper aims at numerically assessing the potential gain in efﬁciency deriving from the application of an unconventional type of nanoparticles, the metal-organic heat carriers (MOHCs), in the ORC ﬁeld. In comparison with standard nanoparticles, these MOHCs make it possible to extract additional heat from the endothermic enthalpy of desorption, with a theoretically high potential for boosting the heat transfer capacity of ORC systems. In this paper a numerical model was developed and customized for considering the adsorption/desorption processes of the pure ﬂuid R245fa (pentaﬂuoropropane) combined with a crystal structure for porous chromium terephthalate (MIL101). The R245fa/MIL101 nanoﬂuid behavior was experimentally characterized, deﬁning proper semi-emipirical correlations. Then, an optimization procedure was developed, combining the numerical model with a PSO algorithm, to optimize the thermodynamic conditions in the ORC so as to maximize the contribution of desorption/absorption processes. The results conﬁrm the increase in net power output (+2.9% for 100 ◦ C) and in expander efﬁciency (+2.4% for 100 ◦ C) at very low heat source temperature. The relevance of tuning the operating cycle and the nanoﬂuid properties is also demonstrated.


Introduction
One of the pillars of the energy transition of the industrial sector is certainly to make industrial processes more circular and more energy efficient. However, nowadays, the process efficiency in this sector is still too low [1,2] and a non-negligible share of the energy losses is represented by low temperature heat, wasted due to the lack of internal heat demand [3].
Even if there are promising technologies for recovering waste heat (such as organic Rankine cycles, ORC) [4,5], their effective and sustainable application is still a challenge due to the strong efficiency decay for low-grade waste heat (T < 150 • C). In particular, as regards ORC, the few available modules in the market have low efficiencies (between 5% and 9%) and high costs (5000 to 8000 €/kW [6][7][8]) with negative consequences in terms of sustainability and profitability.
One of the main reasons for these poor efficiency values is the strong performance decay experienced by pure organic fluids in the case of heat sources characterized by temperatures lower than 150 • C and especially around 100 • C [9].
Several studies have been carried out to identify a working fluid with characteristics suitable for low grade heat sources. Screening criteria based on the fluid thermodynamic and chemical properties [10][11][12][13], the Jakob number [14] and on tuning the heat source temperature and the working fluid critical temperature [15][16][17] have all been proposed. However, it is plain to see that the best working fluid is still far from being identified.
To overcome this issue, molecular engineers have studied how to increase the amount of heat extracted from the heat source and concluded that the contribution of the nanoparticles in terms of sensible heat is too limited at low temperature to have a positive impact on the recovery efficiency, so to increase the performance at these temperature values, they developed an innovative type of nanoparticles, the so-called metal-organic heat carriers (MOHCs), whose structure presents porosity, engineered on purpose to allow the molecule to absorb and desorb the pure fluid in which it is suspended depending on temperature and pressure values [38]. Since the desorption process is endothermic, it allows to extract an additional amount of heat from the heat source with a potential significant increase of the recovery efficiency of the nanofluid in comparison with the pure fluid.
The assessment of the potential increase in ORC performance deriving from the use of these MOHC-based nanofluids in low-grade heat recovery applications started with some preliminary analyses carried out on the pure fluid R245fa (pentafluoropropane) paired with nanoparticles of the MOHC family termed MIL101 [39], a crystal structure for porous chromium terephthalate whose characteristics are suitable for adsorption/desorption processes [40]. In these studies, semi-empirical correlations were determined in order to model the reversible adsorption/desorption process of the MIL101 in the R245fa [41,42]. These correlations were included in the numerical model of an ORC in order to estimate the impact of the adsorption/desorption process on the system performance. The results highlighted a promising potential of the MOHC nanoparticles in terms of area reduction in the heat exchanger but not significantly in terms of net power output. This limited performance improvement achieved by the R245fa/MIL101 nanofluid was due to the non-optimal mass fraction and the non-optimal thermodynamic conditions preventing from the maximization of the adsorption/desorption processes.
The primary aim of this paper is to move on from the results in [41,42], overcoming the limits of the non-optimal operating conditions and maximizing the ORC efficiency gain by optimally tuning the desorption process for the MIL101/R245fa nanofluid and the cycle parameters (pressure and temperature above all). This goal required the development of a new numerical model allowing to include the nanofluid behavior and in particular the desorption/adsorption processes within the optimization procedure of the ORC cycle parameters.
Moreover, to address future research in nanofluids applied to energy conversion processes, the paper also investigates possible paths for molecular improvements of the MIL101 characteristics. In particular, the paper analyses the benefits deriving from a theoretical increase in MIL101 uptake capacity.
The paper is organized as follows: Section 2 presents the numerical model of ORC developed for pure fluids, whereas Section 3 enters the details of the nanofluid characterization and behaviour, explaining the main equations adopted for developing the numerical model. In Section 4 the enhancements achieved by the use of the nanofluid in ORC are critically discussed and a potential path for the MOHC molecule performance enhancement is investigated.

Numerical Model of ORC for Pure Fluids
The ORC system was investigated in a simple and sub-critical configuration, without any regeneration process, since it was widely demonstrated both in literature and by the market that this is the most effective configuration for low grade heat recovery applications (Figures 1 and 2) [43,44].  The low grade heat source was assumed to be a mass flow rate . m s of water equal to 5 kg/s at the temperature T s,in of 100 • C, whereas for the cold source in the condenser, water at the temperature T c,in of 10 • C was fixed [9,45].
As regards the numerical model, the following assumptions were made: • the pump isentropic efficiency: was assumed equal to 0.8 [46].

•
The expander efficiency was set equal to 0.8 for the pure fluid R245fa [47]. In the nanofluid case, as it will be better explained later on, the expander efficiency was determined within the optimization procedure since it resulted to be positively affected by the adsorption process of the R245fa molecule into the MIL101 porous structure [38]. More details can be found in Section 3.4.

•
As regards the vapour quality at the end of the expansion process, the lower limit was fixed equal to 0.85 in order to avoid damages in the expander due to the excessive presence of liquid droplets [48]. • Pressure drops and heat losses were not considered [49].
To evaluate the cycle performance, the efficiency η syst of the ORC system was determined, considering both the thermal efficiency η th and the heat recovery efficiency χ: where P net is the net power output P net over the thermal power absorbed in the evaporator Q evap and Q s is the heat flux which could be transferred to the working fluid if the heat source was cooled to the ambient temperature (T0 = 20 • C). The model was developed in Matlab environment, whereas CoolProp was used for the thermodynamic data of the pure fluid R245fa [50]. More details about the equations can be found in [41,42].

Numerical ORC Model for the R245fa-MIL101 Nanofluid
Unlike pure organic fluids and conventional fluid mixtures, the properties of a nanofluid cannot be obtained by standard commercial software like CoolProp and hence to numerically investigate its behavior is necessary to determine proper correlations modelling the characteristics as a function of pressure and temperature values.

Experimental Determination of Adsorption Characteristics of the Nanofluid R245fa-MIL101
As explained in the introduction section, the potential of MOHC particle is related to the possibility of extracting/releasing additional heat from the heat source and to the cold source as a consequence of the adsorbing/desorbing process in their porous structure [38]. Since the uptake capacity of the porous structure depends not only on the operating conditions (pressure and temperature values) but also on the pure fluid characteristics and on its matching with the structure of the adopted MOHC [51], the adsorption/desorption process has to be experimentally investigated in order to numerically model it, so in previous studies, the MIL101 [52] was selected as a MOHC nanoparticle and its behaviour in combination with the pure fluid R245fa was experimentally analysed, as summarized in the following sub-sections [41,42] In literature, the adsorption behaviour of MOHC-based nanofluids can be traced to a few theoretical models, the most representative of which are summarized in Figure 3 [51]. As it can be seen in this figure, each behaviour is characterized by a different shape of the isotherm, describing the amount of gas adsorbed in the molecule as a function of its pressure at constant temperature. To determine the behaviour of the MIL101 in R245fa, gas adsorption experiments were carried out at four different temperatures (10, 25, 40, 60 • C) [41], resulting in the isotherms shown in Figure 4. The uptake is reported as a function of the ratio between the pressure p and the saturation MIL101 pressure p 0 at different adsorption temperatures. Similar experimental analyses were also carried to analyse the desorption behaviour and the resulting curves were almost superimposed with the adsorption ones, confirming the lack of hysteresis in the nanofluid behaviour. The comparison between the isotherm shapes in Figure 4 and the adsorption models of Figure 3 identified the R245fa-MIL101 as a type I at all temperatures, as shown for comparison in the Figure 4 on the right. This type is characterized by a monolayer adsorption process [51] which can be modelled by the Langmuir equation [53]: where a is the R245fa uptake at pressure p, C is a constant and a m is the maximum monolayer mass capacity per adsorbent mass unit, which is also constant. More details on the major assumptions of this model can be found in [53].
The value of the constants in the Langmuir equation (C and a m ) were determined by fitting the experimental data at the different temperature values, resulting in the curves reported in Figure 4.

The R245fa/MIL101 Enthalpy of Desorption/Adsorption
The added enthalpy gain ∆h a related to the desorption/adsorption process can be determined by means of the Clausius−Clapeyron equation: and it depends on the correlation between pressure p and temperature T, at a certain value of the R245fa molar uptake in the MIL101 molecule w (R = 62.02 J/kgK-specific 245fa constant).
According to the Langmuir's model, representing the R245fa/MIL101 nanofluid (Section 3.1.1), the uptake depends on pressure p and temperature T according to the following equation: where the correlation coefficients were determined starting from the experimental data [41,42]. The R245fa mass uptake m R245 f a by weight of m MIL101 dry mass can hence be easily determined with the following equation: where MM R245 f a is the R245fa molar mass. Then, for a number of values of R245fa molar uptake, the enthalpy of desorption ∆h a was determined by means of the Equation (5), resulting in an average value of about 8.7 kcal/mol.
In the case of adsorption, the average absolute value remains the same because of the lack of hysteresis. However, the adsorption process cools down the fluid and hence the enthalpy of adsorption has a negative value (−8.7 kcal/mol) [38,54].

Nanofluid Modelling Equations: The Heat Exchangers and the Expander
To model the performance of the R245fa/MIL101 nanofluid in an ORC system, it is necessary to adopt equations for estimating the total enthalpy variation in the different components of the system, adding to the conventional enthalpy variation of the pure organic fluid ∆h 2,3 , with the nanofluid contribution in terms of sensible heat and desorption/adsorption process.
As regards the nanofluid enthalpy increase ∆h 2,3 n f in the High Temperature Heat Exchanger (HTHE), the following literature correlation [38] was adopted: where: is the fraction between the MIL101 mass m MIL101 and the total fluid mass m f ; • ∆m R245 f a is the R245 mass uptake difference between the inlet and the outlet of the HTHE determined exploiting Equation (7). • T 2 and T 3 n f are the inlet and outlet nanofluid temperature in the HTHE; • c p,MIL101 = 1.1 J g K is the average MIL101 specific heat capacity in the HTHE defined according to literature values [38,55]; • ∆h a = 8.7 kcal/mol is the average enthalpy of the R245fa/MIL101 desorption process caused by the temperature increase in the HTHE.
Equation (8) was also applied in the low temperature heat exchanger (LTHE), in which the decrease in temperature favors an exothermic adsorption process.
As regards the expander, Equation (8) can also be applied for the definition of the enthalpy variation ∆h 3,4 n f : However, its evaluation is more complex due to the conflicting influence of the pressure and temperature variation on the desorbing/adsorbing behaviour of the MIL101 nanoparticles. During the expansion phase, the nanofluid temperature decreases favoring the adsorption process ( Figure 1). Simultaneously, the process is accompanied by a pressure decrease which favors the opposite desorption process (Figure 1). Depending on the pressure and temperature variations of the ORC during the expander process, one process prevails on the other, causing a positive/negative enthalpy variation due to the desorption/adsorption process and also affecting the turbine isoentropic efficiency. Section 3.4 will show how the process is handled in the optimization procedure.
Equations (8) and (9) clearly show the significant influence of the R245fa mass uptake ∆m R245 f a on the benefits provided by the nanoparticle suspension. If the mass uptake ∆m R245 f a is smaller than the mass m MIL101 of MIL101 suspended nanoparticles, the term ∆m R245 f a m MIL101 − 1 in Equations (8) and (9) is negative, reducing or even annulling the benefits related to the nanoparticle suspension (second and third terms in Equations (8) and (9). But the R245fa mass uptake ∆m R245 f a depends on the operating conditions of the cycle (pressure and temperature values), clearly confirming the need of the optimal tuning of the operating conditions for the nanofluid adoption, as demonstrated in the preliminary analyses [41,42].

The R245fa-MIL101 Numerical Model: Assumptions and Iterative Procedures
The numerical model developed for the pure fluid (Section 2) was properly modified, introducing the correlations developed for the ORC components (Equations (8) and (9)) in order to take into account the MIL101 nanoparticle suspension in the R245fa fluid.
The complexity of the model requires some simplifying assumptions: • The R245fa/MIL101 specific heat can be determined by adopting the following equation [56][57][58][59]: However, since the mass load fraction ϕ was fixed to be smaller or at most equal to 2% for tackling stability problems and since the MIL101 specific heat (c p,MIL101 = 1.1 kJ kg K ) is smaller than that of the R245fa, the specific nanofluid heat deriving from eq. 7 resulted to have a negligible difference with the R245fa specific heat. For this reason, the following assumption was made: resulting in an error smaller than 1%.
• Since the specific heat affects the isobar shape in the T-s diagram, the previous assumption (Equation (8)) also allows to consider, for the nanofluid operating points, the same isobars of the pure R245fa fluid during the heat exchanges in the HTHE and in the LTHE.

•
In the pumping process, the enthalpy contribution due to the exothermic adsorption process was considered negligible.
In addition to the previous assumptions, it was also necessary to tackle the interdependence between some parameters in the Equations (8) and (9).
In particular, for the HTHE, the outlet temperature T 3 n f can be determined once the enthalpy at the outlet h 3 n f is known: However, the enthalpy gain ∆h 2,3 n f depends in turn on the outlet temperature T 3 n f directly and indirectly via m f ads = f p evap , T 3 n f (Equation (7)). To overcome the problem related to this interdependence, an iterative procedure was developed and included in the numerical model.
The main steps of the procedure are here briefly summarized: (1) The R245fa mass uptake at the inlet of the HTHE is determined by applying Equations (6) and (7): [%] = w 2 ·MM R245 f a where (T 2 , p evap ) are the known value of inlet temperature and pressure.
(2) As regards the outlet, in the first iteration, the outlet nanofluid temperature (T 3 n f ) in the HTHE is assumed equal to the outlet temperature in case of pure fluid R245fa T 3 : (3) The R245fa mass uptake at the outlet of the HTHE is hence determined: h 3,n f = ∆h 2,3 n f + h 2 (6) Starting from the enthalpy and pressure values, it is hence possible to determine the temperature of the nanofluid at the HTHE outlet: If this value differs more than 0.1 from the value of the previous iteration i: the procedure starts again from step 3 by assuming: Otherwise, it ends. The same interdependence characterizes the outlet conditions in the expander (T 4 n f , h 4 n f ) and in the LTHE (T 5 n f , h 5 n f ). In both cases, once determined the inlet conditions from the previous iteration process, the procedure previously described was adapted and applied for determining the component outlet conditions.

Optimization Procedure for the Nanofluid
The results obtained by preliminary studies [41] highlighted the promising potential of the MIL101 nanoparticles particularly in terms of area reduction, but did not allow to fully appreciate their potential since the results were affected by two main factors: • the cycle parameters (pressures, temperatures, etc.) used for the nanofluid were the optimal for the pure fluid and were not tuned to optimize the nanoparticles contribution; • the influence of the mass fraction loading ϕ on the nanofluid contribution was not considered.
Therefore, to overcome these issues, the nanofluid model, described in Section 3.4, was included in the optimization procedure to tune the cycle parameters and the nanofluid behaviour. The optimization procedure was aimed at maximizing the power output, through the minimization of the loss in terms of system efficiency η syst : The optimization variables in the procedure (independent variables in Figure 5), were the evaporation pressure p evap , the pinch point temperature differences ( ∆T pp evap , ∆T pp cond ) in both the heat exchangers, the approach point in the HTHE ∆T ap evap and the MIL101 mass fraction ϕ. Whereas the approach point at the condenser ∆T ap cond was fixed equal to 25 • C and was included in the optimization procedure together with other input data, listed in Figure 5.
Knowing the input data and by assuming starting values of the optimization variables, the algorithm estimates the first value of the objective function (Equation (13)) following the steps summarized in Figure 5. Then, the optimization algorithm continues searching for the optimal values of the independent variables till the maximum system efficiency is reached. As regards the optimization algorithm, the ASD-PSO, already used in previous ORC optimization analyses [17], was adopted [60].
The limits in the search domain for the optimization variables were fixed as follows: • the evaporator approach point was assumed to vary from 10 •C to 25 •C and both the evaporator and condenser pinch point varies from 5 • C to 20 • C; • the evaporation pressure boundaries are driven by several constraints, fixed as follows: the lower bound was defined to ensure the pressure to be always greater than the condensing pressure. Since the approach point at the condenser ∆T ap cond and the cold source temperature T c,in were known, the limit was fixed equal to the pressure at a saturation temperature T 1,max = T c,in + ∆T ap cond the upper bound has to guarantee two different conditions: (i) the ORC cycle must be sub-critical and hence the evaporation pressure should be smaller than the critical pressure by at least 1 bar, as suggested in literature [61]; (ii) the vapour characteristic at the HTHE outlet must be saturated or superheated with an evaporation pressure smaller or at least equal to that at the saturation temperature of T 3,min = T s,in − ∆T ap evap . In order to guarantee both the conditions, the upper bound was fixed equal to the minimum between the two values.
This choice ensures that, even when p evap reaches its maximum value, state point 3 is in the vapour region. The working fluid in state 3 is saturated if p evap = p evap,max and superheated in the other cases.
• the mass fraction load (ϕ) of nanoparticles was supposed to vary between 0.05% and 2%: the minimum value of the mass fraction load was fixed to have a minimum impact of the MIL101 on the pure fluid performance, whereas the maximum value was fixed to avoid stability problems according to previous studies on the topic [62].

Results
To analyse the influence of the heat source temperature on the R245fa/MIL101 contribution, optimization analyses were carried out at five different heat source temperatures (from 100 • C to 150 • C). The cycle parameters were optimized to maximize the system efficiency and the results achieved by the nanofluid R245+MIL101 were then compared with those achieved by the pure fluid.
First, as expected, the different behaviours of the two fluids (pure fluid vs. nanofluid) result in optimized operating cycles with slightly different values of some parameters. This is clear, for example, in Figure 6, comparing the optimized cycles in the T-s diagram, and confirms the need of the developing a numerical model for optimizing the nanofluids performance. As regards the numerical results (Tables 1 and 2), for the same heat source temperatures, the introduction of nanoparticles in R245fa resulted in higher net power output especially when considering ultra-low heat source temperature (+2.9% for 100 • C, +2.3% for 110 • C, +1.7% for 120 • C). This positive contribution is achieved by maximizing the mass load fraction ϕ which is always equal to the maximum value (2%) in the allowed range (0.5-2%). However, this contribution becomes lower and lower for increasing heat source temperatures and, at a certain point, is of the same order of the model uncertainty (0.2% for 150 • C). To explain these results, it is necessary to enter into the details of the behaviour of the working fluid.
For the pure fluid, at a fixed heat source temperature, the optimization algorithm searches for the best compromise between two conflicting aspects: on one side, the need of limiting the value of the evaporation pressure in order to maximize the heat recovery efficiency and, on the other side, the need of maximizing the enthalpy drop (by increasing the evaporation pressure) to maximize the thermal efficiency.
In the case of the nanofluid, there are other two contributions to account for: • the endothermic desorption process of the MOHCs, which decreases for increasing evaporation pressure values; • the increase in expander efficiency for increasing evaporation pressure.
As a consequence of the thermodynamic behaviour of MOHC, the entropy of the system is affected by the occurrence of desorption/adsorption processes [40]. In particular, when the adsorption process takes place, the total entropy of the system is reduced. This could happen in the expander, where the temperature decrease favours the R245fa adsorption. If the pressure conditions, also affecting the R245fa uptake, are favourable, the adsorption process limits the entropy increase during the expansion, with a slight increase in expander isentropic efficiency and benefits for the overall efficiency [38]. This is clear in Table 3 reporting the resulting values of the expander efficiencies because of the adsorption process in the optimization analyses for the different heat source temperatures.
Among these two further contributions, the second one seems to be predominant, and this explains why, in the nanofluid case, the optimized cycle is characterized by evaporation pressure values from 1% to 4% higher than those of the pure fluid (Tables 1 and 2). Therefore, the enthalpy drops increase, allowing to achieve higher thermal efficiency η syst but lower heat recovery efficiency χ (1% to 2%). However, at increasing heat source temperatures, the MOHC desorption contribution further decreases because of the increase in the evaporation pressure but this decrease is not compensated by the MOHC adsorption contribution due to the unfavourable pressure and temperature values in the first part of the expansion process. This is clearly seen in Figure 7, which presents the uptake of R245fa in MIL101 at the inlet of the HTHE (point 2) and at the expander inlet and outlet (points 3 and 4) for heat sources of 100 • C, 130 • C and 150 • C. In all the three cases, at the HTHE the uptake always decreases (Figure 7): indeed, during the heat transfer at constant pressure, the temperature raises and thus, considering the experimental adsorption isotherms (Figure 4), desorption occurs allowing to extract extra heat at the heat source.
At the expander, pressure and temperature induce opposite effects on the uptake. While the pressure drop induces desorption, the temperature decrease induces adsorption, which is preferable in the expander since it makes the expansion process more isentropic and the enthalpy drop greater.
As it can be seen, for heat source temperatures of 100 • C and 130 • C (Figure 8a,b), adsorption occurs (m f ads,3 < m f ads,4 ), allowing enthalpy drops and thus power outputs greater than those with the pure fluid. For 140 • C and higher temperatures (Figure 8c), because of the unfavourable pressure and temperature variations, desorption occurs at the expander (m f ads,3 > m f ads,4 ). Because of this, the positive contribution of the MOHC nanoparticles becomes almost negligible, with a performance increase of the same order of the model uncertainty (0.2% for 150 • C).
This behaviour is also reflected in the expander efficiency values ( Table 3). The increase in the expander efficiency compared to the assumed starting value of pure fluid in expander design operating condition (80%) is not negligible at ultra-low heat source temperatures (+2.4% for 100 • C, +1.63% for 110 • C), but becomes less significant for 140 • C (+0.59%) and almost negligible for 150 • C (+0.34%). Thus, the results obtained shows that the nanofluid influence is more effective when the source temperature is particularly low (<130 • C).

Double Uptake Scenario
Research on nanofluids is still ongoing and new improved nanoparticles are and will be synthesized. Since the year 2000s, where the first experimental studies were conducted more than 2000 publications on experimental studies with nanofluids have been presented, of which more than 50% have appeared during the last 5 years [63].
In this section, a scenario considering a theoretical nanoparticle doubling the amount of working fluid adsorbed or desorbed compared to MIL101, is not impossible and is here considered in order to appreciate the influence of enhanced uptake capability in the improvement of the ORC technology performance.
To account for this scenario, the same numerical model presented in Section 3.4 has been adopted, but a double amount of desorbed/adsorbed mass was hypothesized, by modifying Equation (6) as follows: where c = 2.
With these enhanced uptake properties of the MOHC nanoparticles, the optimization was carried out again and the results are presented in Table 4. The positive influence of the nanofluid on the ORC efficiency is confirmed up to the heat source temperature of 120 • C, with a maximum increase in system performance for 100 • C (50.83 kW, +3.51% compared to optimal pure fluid). However, the enhanced adsorption/desorption nanofluid processes are limited at low temperatures (+0.11% compared to the standard MIL101 for 100 • C) and become even detrimental for heat source temperatures greater than 120 • C (system efficiency always smaller than that of the standard MIL101). To understand the reason of these results, the uptake of the enhanced MOHC particles is compared with that of the standard MIL101 during the evaporation and expansion processes ( Figure 8).
It is clear that the general trend of the adsorption/desorption process is similar for both the theoretical and the enhanced nanofluid. In both cases, in the HTHE the nanofluid experiences a desorption process, whereas in the expander an adsorption process only occurs up to a heat source temperature of 130 • C.
However, despite this similar trend, a non-negligible difference emerges: the boosted desorption process increases the temperature values at the end of the evaporation process and at the beginning of the expander process. Despite these increased values, for heat source temperatures lower than 120 • C, the temperature and pressure evolution during the expansion phase favours a strong adsorption process, resulting in enthalpy drops and isentropic efficiencies higher than those of the standard MIL101. As a result, in these cases the performance of the theoretical nanofluid is slightly better or at least comparable with those of the standard MIL101.
On the contrary, for heat source temperatures greater than 120 • C, the desorption process in the HTHE causes a too significant increase in the outlet temperature, negatively affecting the following adsorption process in the expander. Because of this, at the outlet of the expander the temperature T 4 is still high and the final uptake value is smaller than that of the standard MIL101. So, for these heat source temperatures, the enhanced adsorption/desorption properties of the nanofluid, combined with the unfavorable high evaporation pressure values (characterizing these operating conditions), negatively affect the ORC performance.
In order to mitigate this negative influence, the optimization algorithm minimizes the mass loading fraction (ϕ = 0.5%) and lowers the evaporation pressure p ev even below that of the pure fluid (8.78 bar vs. 10.48 bar, −13.37% for T s,i = 130 • C) so as to push on the increase in the heat recovery efficiency. Nevertheless, the power output turns out to be lower than that of the standard MIL101 (−1.15% for T s,i = 130 • C; −0.20% for T s,i = 140 • C; −0.20% for T s,i = 150 • C).
In conclusion, the doubled uptake scenario is not promising in terms of enhancement of the nanofluid performance in an ORC system and further enhancement of the nanofluid properties are not supposed to pursue an increase in the uptake capability, at least for heat source temperatures greater than 120 • C.

Conclusions
The possibility of recovering the large amount of low-grade heat, which is wasted nowadays, certainly represents a significant target for supporting the transition towards a climate neutral society. However, despite the huge potential, the most promising technologies, such as ORCs, are still not able to successfully exploit it due to the negative influence of the low grade of the heat source temperature on the system efficiency.
To deal with this low grade, nanofluids have been more than once suggested as a promising solution. However, in similar applications, standard nanoparticles have proven not to perform as well as expected.
To overcome this issue, molecular engineers have suggested the possible use of the metal organic heat carrier (MOHC) nanoparticles, capable of extracting not only the nanoparticle sensible heat but also a second contribution due to the endothermic desorption process from the MOHC porous structure.
In this paper, to assess the potential gain deriving from the MOHC application, a MOHC-based nanofluid was investigated, i.e., the R245fa/MIL101. Since the uptake capacity of the porous structure depends not only on the operating conditions (pressure and temperature values) but also on the pure fluid characteristics and on its matching with the structure of the adopted MOHC, the adsorption/desorption process of the R245fa in the MIL101 was experimentally investigated and proper correlations were defined for modelling the nanofluid behaviour in the ORC components. Then, the nanofluid model was included in the optimization procedure to tune the cycle parameters and the nanofluid behaviour with the aim of maximizing the net power output.
Optimization analyses were carried out at five different heat source temperatures (from 100 • C to 150 • C). The cycle parameters were optimized to maximize the system efficiency and the results achieved by the nanofluid R245+MIL101 were then compared with those achieved by the pure fluid.
The comparison brought to the following conclusions: • The introduction of the MIL101 nanoparticles in the pure fluid R245fa resulted in higher net power output especially when considering ultra-low heat source temperature (+2.9% for 100 • C, +2.3% for 110 • C, +1.7% for 120 • C).

•
The main reason of the higher net power output is the increase in expander efficiency due to the adsorption process limiting the entropy increase during the expansion (+2.04% for 100 • C, +1.63% for 110 • C, +1.25% for 120 • C, +1.41% for 130 • C).

•
The need of tuning the system thermodynamic operating conditions and the nanofluid behaviour is confirmed by the different evaporation pressure values characterizing the nanofluid case: from 1% to 4% higher than those of the pure fluid. • For heat source temperatures higher than 120 • C, the constraints of the thermodynamic cycle force to have higher evaporation pressure. This creates unfavorable pressure and temperature values in the expansion, resulting in a desorption process and an almost negligible contribution of the MOHC nanoparticles (+0.2% of net power output for 130 • C). To overcome the constraints.

•
As regards the MIL101 mass volume fraction ϕ, the maximum allowed value (2%) turned out to be beneficial for the system performance. However, mass volume fractions higher than 2% are at the time being not feasible due to technical limits. Above all, stability and deposit issues are still far from being fully solved in a challenging environment as an ORC system, in which the nanofluid passes through not only heat exchangers but also fluid machines for thousands of operating cycles.
To trace the path for further MIL101 enhancement by molecular engineering, further analyses considering a theoretical enhancement of the adsorption/desorption attitude of the MOHC were carried out. However, the improvements achieved were almost negligible at 100 • C, so discouraging to pursue this research direction.
In terms of future developments, even if MOHC-based nanofluids have a clear potential of increasing the ORC efficiency in waste heat recovery applications, their successful application requires an optimal tuning between the heat source temperature and the nanofluid behavior. To do this, it will be necessary to: • Experimentally analyze the adsorption/desorption processes for different combination of pure fluids and MOHC nanoparticles allowing to model the each MOHC-based nanofluid by means of semi-empirical correlations; • Experimentally investigate the stability and the performance of MOHC-based nanofluid in fluid machines, i.e., expanders and pumps, so as to improve the numerical models.