Multiobjective Optimization of a Plate Heat Exchanger in a Waste Heat Recovery Organic Rankine Cycle System for Natural Gas Engines

A multiobjective optimization of an organic Rankine cycle (ORC) evaporator, operating with toluene as the working fluid, is presented in this paper for waste heat recovery (WHR) from the exhaust gases of a 2 MW Jenbacher JMS 612 GS-N.L. gas internal combustion engine. Indirect evaporation between the exhaust gas and the organic fluid in the parallel plate heat exchanger (ITC2) implied irreversible heat transfer and high investment costs, which were considered as objective functions to be minimized. Energy and exergy balances were applied to the system components, in addition to the phenomenological equations in the ITC2, to calculate global energy indicators, such as the thermal efficiency of the configuration, the heat recovery efficiency, the overall energy conversion efficiency, the absolute increase of engine thermal efficiency, and the reduction of the break-specific fuel consumption of the system, of the system integrated with the gas engine. The results allowed calculation of the plate spacing, plate height, plate width, and chevron angle that minimized the investment cost and entropy generation of the equipment, reaching 22.04 m2 in the heat transfer area, 693.87 kW in the energy transfer by heat recovery from the exhaust gas, and 41.6% in the overall thermal efficiency of the ORC as a bottoming cycle for the engine. This type of result contributes to the inclusion of this technology in the industrial sector as a consequence of the improvement in thermal efficiency and economic viability.


Introduction
Natural gas is one of the most profitable fuels to replace conventional fuels, such as diesel and gasoline, around the world [1]. One of the reasons is a large number of reserves of fossil fuels around the world that are becoming considered as a suitable alternative in the industrial energy sector [2]. All devices such as heat exchangers implicate energy losses during the process, which is equivalent to increasing the total operational costs and reducing the energy performance of the system. Therefore, improvement in the thermal performance of a gas engine was proposed by incorporating an exhaust gas heat recovery system using an organic Rankine cycle (ORC) [3]. However, this proposed configuration requires a thermal oil circuit, in which the residual exhaust gas is used to evaporate the organic fluid with the help of a thermal oil that regulates the amount of heat transferred and decreases the ORC evaporation temperature.
Wang et al. [16] performed multiobjective optimization of the ORC condenser using a genetic algorithm, where the pressure drop and the heat transfer area were minimized under constant heat transfer conditions, resulting in a series of optimal solutions presented in a Pareto frontier. However, the author did not base their optimization on a thermo-economic analysis, which indicates that the effect of equipment costs on the results obtained was not considered. Consequently, minimization of equipment costs and the entropy generation number, considering geometrical parameters of the parallel plate heat exchanger (ITC2) without its effect on the heat transfer capacity of the equipment and reducing system performance, had not been developed considering some energy indicators, such as global energy conversion efficiency and energy efficiency, which decrease due to the irreversibility of the process and limits its application and commercial feasibility [17].
From these results and their limitations, the main contribution of this research is the multiobjective optimization of a plate heat exchanger, which was used to evaporate the organic fluid in an ORC as a WHR system from a natural gas engine using the NSGA-II genetic algorithm, and a detailed thermodynamic model of the heat exchanger. The most sensitive variables of the system are determined, and the optimal configuration is selected to obtain the minimum acquisition cost and entropy generation number. In addition, the influence of the design parameters of the thermodynamic cycle, such as the evaporator and condenser pinch temperatures, the turbine and pump efficiency, and different working fluids were studied.

Description of The Waste Heat Recovery System
The configuration shown in Figure 1a [18] operates as the T-s diagram in Figure 1b is presented. An air-fuel mixture (1) is delivered to the internal combustion engine (ICE), which is compressed in the compressor stage of the turbocharger (6) to achieve the conditions required for combustion in the cylinders. The gases at the outlet of the exhaust manifold (9) expand in the expansion stage of the turbocharger to reach the stream (10) and are disposed the environment (11) after transferring the energy in the heat exchanger tube and shell (ITC1) with the thermal oil (3 AT), which flows through the thermal oil circuit pumped by Pump 1 (B1).
The thermal oil (1AT) transfers the energy to evaporate the ORC working fluid in the (ITC2) to obtain the superheated steam in (1ORC) to generate thermal power in the turbine (T1) and electric energy in the generator (G). In the evaporator (ITC2), three zones are presented: the preheating zone (zone 1), the evaporation zone (zone 2), and the superheating zone (zone 3). The working fluid expands until the low pressure (2ORC) is inputted in the condenser heat exchanger (ITC3), where the fluid is fully condensed to the saturated liquid phase (3ORC). Condensation is conducted in two successive stages called condensation (1A to 1gA) and cooling (1gA to 2A). Subsequently, the working fluid pump (B2) increases the fluid pressure to the evaporating pressure of the ITC2, which finishes the WHR system from the natural gas engine by ORC using toluene as the working fluid [19][20][21][22].

Energy Analysis
The components of the WHR system were studied under the assumption of an open system operating under steady state conditions. Consequently, the mass and energy balance must be calculated, as shown in Equation (1) and Equation (2), respectively. ∑ ṁi n − ∑ ṁo ut = 0 (1)

Energy Analysis
The components of the WHR system were studied under the assumption of an open system operating under steady state conditions. Consequently, the mass and energy balance must be calculated, as shown in Equation (1) and Equation (2), respectively.
Some performance indicators were calculated to assess the WHR systems, such as the thermal efficiency of the configuration (η I, c ), heat recovery efficiency (ε hr ), and overall energy conversion efficiency (η I, global ) [23]. The thermal efficiency of the configuration (Equation (3)) is the ratio of the net power output of the ORC ( · W net ) recovered, with respect to the waste heat. In this case, the net power output of the ORC is the turbine power output (T1) less the energy consumption of the thermal oil pump (B1) and organic fluid pump (B2).
In addition, the waste heat recovery efficiency (Equation (4)) is the heat recovered from the exhaust gas engine line, with respect to the maximum heat available to be removed, and the overall energy conversion efficiency (Equation (5)) is the relation between the net power output and the available waste heat.
Also, to consider global indicators involving engine operation conditions, the absolute increase of engine thermal efficiency (Equation (6)) was calculated as a measurement of the net power output to the energy supply by the natural gas in the engine.
Because additional power is delivered in the heat recovery system over the engine power, there is a break-specific fuel consumption (BSFC), which is determined by Equation (7), and generating more power with the same fuel consumption decreases the specific fuel consumption of the engine, which is calculated as follows in Equation (8) [24].

Exergy Analysis
The specific exergy for the process states is calculated by neglecting the variations of kinetic and potential energy, resulting as shown in Equation (9).
The chemical exergy of the exhaust gases as a product of combustion (stream 10), is defined by Equation (10). One has a mixture of gas products of the combustion of natural gas. The chemical exergy of the gas mixture is given by Equation (10).
Also, the exergetic efficiency based on the second law of thermodynamics (η II, ORC ) is calculated with Equation (11), additionally expressed as a function of the destroyed exergy ( · Ex Destr ) using Equation (12).

Modeling of The Plate Heat Exchanger
The parallel plate heat exchangers are devices designed to provide a large surface area of heat transfer per unit volume [25] in addition to achieving high heat transfer rates between two fluids in a small volume [26,27]. The design of this heat exchanger takes into account the phase change in the organic fluid at low temperatures [28] to select the optimal configurations that guarantee the functional operation of the equipment within the thermal oil circuit.
For the heat exchanger studied, Therminol 75 was used as thermal oil and hot fluid for the ITC2 design. This oil recovers a fraction of the heat available in the exhaust gas line of the MCI and transfers the energy to heat the organic fluid, which presents a phase change during the process. For this reason, the thermal design considers different stages of the organic fluid, represented by zones, to facilitate the description of the process and the design of the device, as shown in the T-s diagram in Figure 2 of the heat transfer process between the two working fluids in the ITC2.
The increase of plate spacing has a direct effect on the channel's cross-section area, channel velocity, equivalent diameter, thermal oil, and the organic fluid Reynolds number in the ITC2. Generally, these parameters affect the total heat transfer area, total pressure drop, and exchanger performance [29]. In this section, the determination of geometric parameters is presented, as shown in Figure 3, which affects the number of plates and the heat transfer area of the equipment in each zone. The equations are analyzed under three zones: preheating (zone 1), evaporation (zone 2), and overheating (zone 3).
For the heat exchanger studied, Therminol 75 was used as thermal oil and hot fluid for the ITC2 design. This oil recovers a fraction of the heat available in the exhaust gas line of the MCI and transfers the energy to heat the organic fluid, which presents a phase change during the process. For this reason, the thermal design considers different stages of the organic fluid, represented by zones, to facilitate the description of the process and the design of the device, as shown in the T-s diagram in Figure 2 of the heat transfer process between the two working fluids in the ITC2. The increase of plate spacing has a direct effect on the channel's cross-section area, channel velocity, equivalent diameter, thermal oil, and the organic fluid Reynolds number in the ITC2. Generally, these parameters affect the total heat transfer area, total pressure drop, and exchanger performance [29]. In this section, the determination of geometric parameters is presented, as shown in Figure 3, which affects the number of plates and the heat transfer area of the equipment in each zone. The equations are analyzed under three zones: preheating (zone 1), evaporation (zone 2), and overheating (zone 3).    Figure 2 shows the temperature limits in the mixing phase, where the temperature of the thermal oil at the beginning of the phase is calculated using Equation (13) with the pinch (known as the minimum temperature difference), and the temperature at which the organic fluid leaves ITC2 is calculated with Equation (14) [30].
The heat transfer flow is calculated for each particular zone. Equation (15) allows calculation of heat flow in zone 1, and the heat of each zone is a function of the different enthalpies presented in the process.
The enthalpy of the organic fluid at the ITC2 inlet is calculated as Equation (16).
where h 4ORC and h ORCf are in kJ/kg. The logarithmic mean temperature difference is an essential factor in the heat exchanger design, which is calculated for zone 1 as Equation (17).
where LMTD Z1 is given in Kelvin. Other relevant factors in designing the heat exchanger are the plate area and the heat transfer area in m 2 , which are calculated using Equations (18) and (19),  Figure 2 shows the temperature limits in the mixing phase, where the temperature of the thermal oil at the beginning of the phase is calculated using Equation (13) with the pinch (known as the minimum temperature difference), and the temperature at which the organic fluid leaves ITC2 is calculated with Equation (14) [30].
The heat transfer flow is calculated for each particular zone. Equation (15) allows calculation of heat flow in zone 1, and the heat of each zone is a function of the different enthalpies presented in the process.
The enthalpy of the organic fluid at the ITC2 inlet is calculated as Equation (16). where h 4ORC and h ORCf are in kJ/kg. The logarithmic mean temperature difference is an essential factor in the heat exchanger design, which is calculated for zone 1 as Equation (17).
where LMTD Z1 is given in Kelvin. Other relevant factors in designing the heat exchanger are the plate area and the heat transfer area in m 2 , which are calculated using Equations (18) and (19), respectively.
where plate width (W) and plate height (H) are in mm, and U Z1_A represents the overall heat transfer coefficient for zone 1 in W m 2 K . Through the application of Equations (20)- (31) it is possible to calculate the most relevant geometric characteristics of the equipment, such as the plate thickness, number of plates per zone, number of passes, the velocities, pressure drops in the system, and the exergy destroyed in each zone during the process The velocities determine the fluid flow regimen and their characteristics. Therefore, the design considers the calculation for each zone of the dimensionless numbers of Reynolds (Equation (23)), Prandtl (Equation (24)), and Nusselt (Equation (25)) [31]. This allows us to obtain valuable information about the fluid properties, the flow characteristics, and the fluid heat transfer capacity.
where the hydraulic diameter (Dh) is in m, the fluid velocity (v Z1_AT ) is in m/s, the kinematic viscosity (υ Z1_AT ) in m 2 /s, density (ρ Z1_AT ) is in kg/m 3 , the thermal conductivity (k Z1_AT ) is in W/m·K, and the specific heat at constant pressure Cp Z1 AT is in J/kg·K. To calculate the drop pressure at the inlet, and because of the flow within zone 1 for the thermal oil side, Equations (27) to (28) are used, which requires the velocity calculation at the ITC2 input of the thermal oil using Equation (26), where vin Z1_AT is in m/s. where the thermal oil mass flow ( · m AT ) is in kg/s, hydraulic diameter (Dh) is in m, and the drop pressure (∆P) is obtained in mbar.
The total pressure drop in the heat exchanger is the total pressure loss on the thermal oil side (Equation (29)) and the total pressure drop on the organic fluid side (Equation (30)), involving the pressure drops due to flow (∆PF) and the pressure drops at the fluid inlet (∆Pin) in each of the zones.
Consumption of the potential useful work in zone 1 of ITC2 is calculated by Equation (31), which represents the exergy destroyed by the system in this zone expressed in watts. Similarly, the destroyed exergy is calculated for the other ITC2 zones, facilitating the computation of one of the objective functions considered in the optimization of this equipment.

Methodology of Design
The objective of the design of the parallel plate heat exchanger is to provide the lowest amount of entropy generated during the process [32] with the smallest heat transfer area and, therefore, the best cycle efficiencies [25]. For this reason, it is necessary to determine a geometric configuration that accomplishes the thermodynamic requirements with the lowest possible acquisition cost [33]. It is for this reason that the methodology implemented for this design is multiobjective optimization applied in various disciplines to minimize or maximize two or more functions simultaneously [34], which can be expressed as shown in Equation (32), and is subject to restrictions g i (X) ≤ 0; i = 1, . . . , m and h i (X) = 0; and j = 1, . . . , n within a range of criteria for X k, min ≤ X k ≤ X k, max . For these systems, the calculation of the entropy generation number (EGN) is proposed using Equation (33), and the heat exchanger acquisition cost is expressed as in Equation (34) [35].
The results obtained by this technique display a series of geometric combinations that provide the functionality of the proposed objectives, which form a set of points called a Pareto front. In determining the most optimal points, however, a multicriteria decision technique is required to facilitate decision-making and, therefore, the optimization of the parallel plate heat exchanger.

Multi-Criteria Decision
The multicriteria decision method called the technique of order preference for similarity to ideal solution (TOPSIS) [36][37][38], is applied to select a point for the ideal geometrical configuration, which is a suitable method for this type of applications and is used to classify the alternatives to the Pareto solutions obtained. This technique mathematically identifies the point with the nearest distance to the positive ideal solution with the longest distance to the negative ideal solution, conditions which are calculated by Equation (35) where d ix and d iy are the distances from the points to the ideal positive and negative solution, respectively; t ij is the reference value of alternative i for objective j; and t xj and t yj are the ideal and nonideal values, respectively. The relative proximity to the ideal solution (S iy ) is calculated using Equation (36) S iy = d iy where the best solution is the one whose S iy is the closest to 1.

Working Fluid Selection
This research mainly focused on improving the energetic and exergetic performance of the WHR system, and selection of the working fluid was oriented to the increase net power [39]. For this case, 14 preselected working fluids were evaluated in simulation to study the performance of the WHR system for the Jenbacher JMS 612 GS-N.L engine, using the simple ORC configuration under typical engine operating conditions. The results of the performance indicators studied are shown in descending order (Table 1). The results show the MD4M and D6 fluids had the lowest net power with values of 30.84 kW and 32.04 kW. These results were near to those obtained in the ASPEN HYSYS ® simulation with the CATERPILLAR C32 ATAAC diesel engine [40], where a simple ORC was used to make use of the residual heat available in the engine exhaust gases.
Based on these results, acetone, cyclohexane, and toluene were the organic fluids that presented the best values in the analysis performed. Therefore, they were selected to analyze the influence they had on the performance of the simple ORC configuration.

Parametric Study
In this section a parametric study was developed to visualize the performance of the system with the selected organic fluids and the performance of the parameters in the presence of vibration; this allowed us to determine the fluids that provided better results and had the best properties for the working conditions. Figures 4-7 show the influence of variation on the performance parameters in the ORC cycle using toluene, cyclohexane, and acetone. The parameters selected to visualize the behavior of the ORC cycle were the net power, the absolute increase of the thermal efficiency, and the overall exergetic efficiency. Toluene was the fluid that guaranteed the best performance of the system in all parametric cases, as it maintained the highest values compared to other organic fluids.
In this section a parametric study was developed to visualize the performance of the system with the selected organic fluids and the performance of the parameters in the presence of vibration; this allowed us to determine the fluids that provided better results and had the best properties for the working conditions. Figures 4-7 show the influence of variation on the performance parameters in the ORC cycle using toluene, cyclohexane, and acetone. The parameters selected to visualize the behavior of the ORC cycle were the net power, the absolute increase of the thermal efficiency, and the overall exergetic efficiency. Toluene was the fluid that guaranteed the best performance of the system in all parametric cases, as it maintained the highest values compared to other organic fluids. Figure 4 shows the performance of the simple ORC cycle with the three fluids selected at different turbine efficiency values (T1). It can be noted that the maximum values of power, the absolute increase in thermal efficiency, and overall exergetic efficiency were obtained when the turbine efficiency was 85%. That is, as the efficiency of this equipment increased, increases in net power and all system efficiencies were observed. Figure 4a shows how the net power values changed concerning the variation of efficiencies for each fluid used. In the case of toluene, an increase from 75% to 85% in turbine efficiency allowed an 11.9% increase in the net power of the system, from 101.02 kW to 113.12 kW, specifically.  The values allow us to visualize the performance of the system corresponding to the variations in efficiency of the pump, as shown in Figure 5. It can be observed in Figure 4b,c that the values of ∆η th and η II.overall , respectively, were similarly affected, indicating that the variation of efficiencies in the turbine and the pump had a similar influence on the thermal parameters of the system. On the contrary, net power was influenced by the variations, as shown in Figure 5a, which reached a maximum value of 106.57 kW with an 85% efficiency in the pump (B2). Effects of the variations of the condenser pinch temperature are shown as in Figure 6. It can be observed that toluene was the fluid that guaranteed the best performance of the ORC cycle, supplying a better net power and higher efficiencies of the cycle. In addition, it is indicated that variations in condenser pitch values were not significant in the performance of the ORC. Figure 6a evinces that, with a condenser pinch temperature of 60 °C and toluene as the working fluid, the net power reached  The results obtained in the parametric study allow us to affirm that the organic fluid with the most suitable properties for the established working conditions is toluene. Therefore, for thermodynamic calculations, the values and properties corresponding to toluene as an organic fluid were taken for this research.

Results of The Optimization
The multiobjective optimization of the parallel plate heat exchanger was performed considering two objective functions with four decision variables, which are detailed as in Table 2. The selected variables correspond to the geometric parameters of the ITC2, all to obtain the optimal geometric  The results obtained in the parametric study allow us to affirm that the organic fluid with the most suitable properties for the established working conditions is toluene. Therefore, for thermodynamic calculations, the values and properties corresponding to toluene as an organic fluid were taken for this research.

Results of The Optimization
The multiobjective optimization of the parallel plate heat exchanger was performed considering two objective functions with four decision variables, which are detailed as in Table 2. The selected variables correspond to the geometric parameters of the ITC2, all to obtain the optimal geometric  Figure 4 shows the performance of the simple ORC cycle with the three fluids selected at different turbine efficiency values (T1). It can be noted that the maximum values of power, the absolute increase in thermal efficiency, and overall exergetic efficiency were obtained when the turbine efficiency was 85%. That is, as the efficiency of this equipment increased, increases in net power and all system efficiencies were observed. Figure 4a shows how the net power values changed concerning the variation of efficiencies for each fluid used. In the case of toluene, an increase from 75% to 85% in turbine efficiency allowed an 11.9% increase in the net power of the system, from 101.02 kW to 113.12 kW, specifically.
The values allow us to visualize the performance of the system corresponding to the variations in efficiency of the pump, as shown in Figure 5. It can be observed in Figure 4b,c that the values of ∆η th and η II.overall , respectively, were similarly affected, indicating that the variation of efficiencies in the turbine and the pump had a similar influence on the thermal parameters of the system. On the contrary, net power was influenced by the variations, as shown in Figure 5a, which reached a maximum value of 106.57 kW with an 85% efficiency in the pump (B2).
Effects of the variations of the condenser pinch temperature are shown as in Figure 6. It can be observed that toluene was the fluid that guaranteed the best performance of the ORC cycle, supplying a better net power and higher efficiencies of the cycle. In addition, it is indicated that variations in condenser pitch values were not significant in the performance of the ORC. Figure 6a evinces that, with a condenser pinch temperature of 60 • C and toluene as the working fluid, the net power reached a maximum value of 106.5 kW, which is 18.86% higher than the value obtained using cyclohexane and 18.09% higher than that obtained with acetone. The influence of the working fluid and the variations of the evaporator pinch temperature on the performance of the ORC cycle can be observed as in Figure 7. The best results were obtained when AP increased to 35 • C, and toluene was used as the working fluid. It is asserted that the variations in temperatures pinch either in the condenser or the evaporator did not represent a determinant influence on the performance of the cycle because both variables maintained an inversely proportional relation. On the other hand, the working fluid did influence the performance of the ORC cycle, as is the case with net power.
The results obtained in the parametric study allow us to affirm that the organic fluid with the most suitable properties for the established working conditions is toluene. Therefore, for thermodynamic calculations, the values and properties corresponding to toluene as an organic fluid were taken for this research.

Results of The Optimization
The multiobjective optimization of the parallel plate heat exchanger was performed considering two objective functions with four decision variables, which are detailed as in Table 2. The selected variables correspond to the geometric parameters of the ITC2, all to obtain the optimal geometric configuration that guarantees the operational conditions and the thermodynamic requirements of the secondary thermal circuit. For this case, 130 commercial configurations were available to obtain the reference values of the objective functions. Positive exergy destruction in the three zones of the equipment was considered as a restriction · e Z1 + · e Z2 + · e Z3 > 0. In addition, the angle of inclination of the plates ranged from 10 • to 80 • . On the other hand, the limits of C2, C3, and C4 were given by the geometric dispositions of the heat exchangers that can be purchased in the market with a determined cost. Figure 8 shows the variations of the EGN and the cost of acquisition with the decision criteria, where it can be observed that the angle of inclination, the height of the plate, width of the plate, and length between plates affected the selected objective functions (cost and EGN). Figure 8 shows the influence of design parameters on both the cost of equipment acquisition and the amount of entropy generated within ITC2. An increasing angle of inclination in the grooves of the plates resulted in a lower cost of the equipment, while, at the same time, the entropy generated in the system increased, as in Figure 8a. This is caused by an increase in corrugation on the surface of the plate. Likewise, a similar behavior was shown with an increase in the height of the plate in the ITC2, as in Figure 8b, for variations in the cost of the equipment and variations in the amount of entropy generated. Increasing the width of the ITC2 plate created oversizing, which generated high manufacturing costs and, therefore, acquisition of equipment, causing a significant increase in its values. Similar behavior occurred in the generation of entropy, almost with the same inclination as in Figure 8.  The optimal points have a linear tendency; this is due to the order of the related objective functions in the optimization process, which are of the first order because of the correlation that exists On the other hand, an increasing length between the plates decreased the number of plates required for the heat transfer area. Therefore, an increase in this parameter represents a decrease in the pressure drop per flow, which minimizes the amount of entropy generated during the process (Figure 8d). In our case, it represents a critical design parameter because the acquisition cost increases due to an oversizing of the equipment, which represents a problem for the optimization of this objective.
The multiobjective optimization solution does not have a global optimal point because no single response exists that simultaneously optimizes each objective variable. Therefore, a procedure of selecting the final optimal values by applying the NSGA II [41,42] method is included, where 32 points are obtained from the objective functions after performing 230 possible iterations. Figure 9 shows the Pareto frontier for the cost of acquisition with the exchanger EGN, plus five possible solution points.
The optimal points have a linear tendency; this is due to the order of the related objective functions in the optimization process, which are of the first order because of the correlation that exists between the geometric parameters.
The five points selected at the Pareto frontier meet the optimization requirements, as shown in Figure 9. Applying TOPSIS gives the ideal positive and negative solutions to find the optimal points (Table 3).  The optimal points have a linear tendency; this is due to the order of the related objective functions in the optimization process, which are of the first order because of the correlation that exists between the geometric parameters.   distribution with a tendency close to the lower limit of the domain, as shown in Figure 10a, which suggests that the variable plays a significant role in the compensation of the objective variables. The plate width value tended to approximate the lower values around 200 mm, as shown in Figure 10c, indicating that decreasing this parameter improves the optimization result. The plate height and the length between plates had scattered distributions with tendencies to a higher value, as shown in Figure 10b,d. The observations obtained suggest that these design variables are determinant for identifying the critical points of the objective functions proposed for equipment optimization. The multicriteria decision technique indicated that point D is nearest to the ideal solution. From this point, the optimal values of the decision criteria are obtained and are shown in Table 4.

Conclusions
This work presents a model to optimize a parallel plate heat exchanger (ITC2), which corresponds to the evaporator of a secondary thermal circuit, through the NSGA-II approach. The acquisition cost of the equipment and the entropy generated were selected as objective functions, whereas the geometric parameters of the exchanger were considered as decision variables. As a single The multicriteria decision technique indicated that point D is nearest to the ideal solution. From this point, the optimal values of the decision criteria are obtained and are shown in Table 4. The optimal values of the geometric parameters of ITC2 form a transfer area of 22.04 m 2 , which ensures a heat transfer flow from the thermal oil to the toluene of 693.87 kW. An overall thermal engine efficiency of 41.6% was achieved for the ORC cycle as a result of the minimum amount of possible exergy destroyed by both temperature and pressure drops.

Conclusions
This work presents a model to optimize a parallel plate heat exchanger (ITC2), which corresponds to the evaporator of a secondary thermal circuit, through the NSGA-II approach. The acquisition cost of the equipment and the entropy generated were selected as objective functions, whereas the geometric parameters of the exchanger were considered as decision variables. As a single parameter value cannot satisfy both functions to be optimized, a series of optimal points are presented in the form of the Pareto frontier, which represents the equilibrium curve between both functions. The working fluid has a determinant role in the performance of the simple ORC cycle, so it is necessary to select the organic fluid that provides more significant benefits to the system. A study on performance parameters of the system is carried out before varying the efficiency of both the turbine and the pump as well as the temperature pinches of both the evaporator and the condenser of the system. With this, we attempted to identify which organic fluid had the highest net power values, the highest absolute increase in the thermal efficiency, and highest overall exergetic efficiency. Effects of the decision criteria on the objective functions are also studied by means of a sensitivity analysis, which showed that the length between plates is the most promising criteria, since its increase causes an elevation in the costs of the equipment up to a maximum value of USD 18,000 and a decrease in the entropy generation number.
The results found that, when applying the methodology proposed for this evaporator design through multiobjective optimization and selecting the best configuration of the five possible solutions through the TOPSIS method, point D was the best solution according to the established criteria. It was possible to minimize the entropy generation number (NGE = 0.058) and the acquisition cost of the equipment (USD 10,385.55), with an inclination angle (20.44 • ), plate height (2070.32 mm), plate width (205.16 mm), length between plates (800.49 mm), and heat transfer area of 22.04 m 2 . This guarantees that heat transferred from the thermal oil (Therminol 75) to the toluene is 693.87 kW, overall thermal motor efficiency of the ORC cycle is 41.6%, and pressure drop is 980.32 mbar, which is within the permissible backpressure range of the engine.
The results of optimization may vary with the change of some configurations in the genetic algorithm, for instance: the fraction of the population, the population number, and the number of allowed iterations. In this case, optimization was performed with a total of 200 particles in the iterative space fractionated, achieving a value of 0.16 for a total of 32 optimal points, as shown in the Pareto frontier. In conclusion, the optimal geometry will be different for each case because there are infinite combinations in the iteration space, and adjustments were made in the configuration of the optimization model. This proposed methodology can be applied to the thermodynamic and economic optimization of plate heat exchangers in any type of heat recovery system with indirect evaporation of the organic fluid. This methodology is always more relevant for cases where there are limitations in heat source backpressure, such as industrial engines with medium and high exhaust gases temperatures, and is applicable in cases where ORC technology has not been widely applied commercially.

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

Abbreviations
The following abbreviations are used in this manuscript: Preheating zone in ITC2 Z2 Evaporation zone in ITC2 Z3 Overheating zone in ITC2