Development of Decision-Making Tool and Pareto Set Analysis for Bi-Objective Optimization of an ORC Power Plant

Power plants based on organic Rankine cycle (ORC) are known for their capacity in converting low and medium-temperature energy sources to electricity. To find the optimal operating conditions, a designer must evaluate the ORC from different perspectives including thermodynamic performance, technological limits, economic viability, and environmental impact. A popular approach to include different criteria simultaneously is to formulate a bi-objective optimization problem. This type of multi-objective optimization (MOO) allows for finding a set of optimal design points by defining two different objectives. Once the optimization is completed, the decision-making analysis shall be carried out to identify the final design solution. This study aims to develop a decision-making tool for facilitating the choice of the optimal design point. The proposed procedure is coded in MATLAB based on the commonly used Technique for Order Preference by Similarity to Ideal Solution (TOPSIS). By providing the capability to graphically identify the decisions taken, the tool developed in the study is called Tracking and Recognizing Alternative Design Solutions (TRADeS). Analysis of our data shows that certain regions of Pareto set points should be excluded from the design space. It was noted that in these regions a high rate at which one of the objectives moves away from its ideal value coincides with a low rate at which the second criterion approaches its ideal solution. Hence, it was recommended that the criteria weights corresponding to excluded regions of the Pareto set should be discarded when selecting the final design point. By comparing the results obtained using the proposed model to those of existing decision-making techniques, it was concluded that while the known approaches are appropriate for an easy and fast selection of the final design point, the presented procedure allows for a more comprehensive and well-rounded design. It was shown that our design tool can be successfully applied in the decision-making analysis for problems that aim at optimizing the ORC using two design criteria. Finally, the proposed software benefits from a generic structure and is easy to implement which will facilitate its use in other industrial applications.


Introduction
The depletion of fossil fuels combined with increasing electricity demand is the main reason for the efforts aiming at the utilization of alternative energy sources. As a technology capable of converting renewable and waste heat to electricity, the organic Rankine cycle (ORC) power plant has been widely investigated in the last decade. The growing interest in developing ORC systems is partly because it lends itself well to applications including geothermal energy [1], solar thermal energy [2], biomass energy [3], and low-grade waste heat recovery [4]. Furthermore, the system can be cycle (BFC) considering the net power output and heat exchanger area as two conflicting criteria. By applying NSGA-II, they found the Pareto optimal solutions. To examine various design points, the authors used TOPSIS with a weight factor for the net power output varying in a range from 0.0 to 1.0. The values between 0.1 and 0.6 have been reported as optimal. By limiting the range of the weights, the subjective decision-making process was facilitated. Detchusananard et al. [25] optimized a steam biomass gasification system coupled with the solid oxide fuel cell by solving a bi-objective optimization problem with two economic criteria. To comprehensively evaluate the obtained results, the authors applied ten selection methods, including LINMAP and TOPSIS approaches. Different decision-making results have been compared by varying the weights of the objectives in a range between 0.1 and 0.9. It was reported that for the same weight distributions, LINMAP and TOPSIS provided similar outcomes.
Previous studies on the decision-making process show that it is as essential as the optimization procedure and it should be conducted thoroughly to design the system comprehensively and reliably. Hence, this study aims to develop a decision-making tool to facilitate the choice of the final design point in ORC power plants. The proposed software is based on the well-established TOPSIS procedure. The decisions that can be made during the program operation include covering the desired decision-makers' preferences, selecting the dominant criterion, or specifying a superiority of the selected objective using a rating scale. New recommendations concerning the decision-making process are presented. In particular, exclusion of weight distributions that correspond to adverse Pareto set points is put forward. These design points are located in regions where a steep rate at which the first objective moves away from its ideal value is combined with a gradual rate at which the second criterion approaches its ideal solution. The consideration of design points that are getting closer to the ideal value of one of the criteria is questionable since a slight improvement of this criterion is accompanied by a significant deterioration of the second objective. The literature review clearly reveals that the discussed effect was not addressed in previous studies. Overall, the main contributions of this study can be listed as (1) Developing an intuitive tool that can be applied in the decision-making process for a bi-objective optimization of the ORC or other industrial systems, (2) Making the developed software code available for public use (3) Making recommendations concerning the criteria weights and the selection of final system design.

System Description
The system under consideration is modeled as a basic subcritical ORC power plant. As depicted in Figure 1a, the simple ORC consists of a condenser, pump, vapor generator, and turbine coupled with an electric generator. In each of the system devices, the working fluid of organic origin undergoes thermodynamic processes presented in Figure 1b. During the expansion process (1)(2)(3)(4)(5) in a turbine, the energy of the working medium is converted into mechanical work of the shaft which is used to drive the electric generator. Then, the working fluid is directed to the condenser in which it is cooled down (5-5") and condensed (5"-6), by transferring the heat with cold water. In a pump, the medium is pressurized (6-7) to achieve a high-pressure state before entering the vapor generator. Utilizing the heat of geothermal water in the vapor generator, the working fluid is preheated (7)(8), evaporated (8)(9) and superheated (9-1) to re-enter the turbine and start a new cycle.

Modeling of ORC Components
The expansion device of the examined system is modeled as a radial-inflow turbine (RIT). RIT converts low enthalpy of the working fluid to shaft work in a compact single-stage structure. Hence, the RIT prevailed as an excellent choice for small-scale ORC systems [27]. As seen in Figure 1b, the expansion of the working fluid is divided into sub-processes that take place within the RIT including volute, nozzle, interspace, and rotor. The modeling of the RIT is based on a mean-line (onedimensional) model [28] which assumes that the working fluid parameters change only with respect to mean streamline through the individual stage components.
In the preliminary design of the RIT, most attention is focused on the analysis of the flow through the rotor [29]. The cross-section of the rotor with mean streamline drawn between its inlet (4) and outlet (5) stations is depicted in Figure 2b. Furthermore, the velocity triangles at the rotor inlet (parameters with subscript 4) and outlet (parameters with subscript 5) are also presented. The common practice in designing the RIT is to assume the radial blades at the rotor inlet [30] which is equivalent to set the blade angle β4,blade to 0° (the positive sign of an angle is measured according to the positive direction of peripheral velocity u). For this reason, the relative flow angle and the

Modeling of ORC Components
The expansion device of the examined system is modeled as a radial-inflow turbine (RIT). RIT converts low enthalpy of the working fluid to shaft work in a compact single-stage structure. Hence, the RIT prevailed as an excellent choice for small-scale ORC systems [27]. As seen in Figure 1b, the expansion of the working fluid is divided into sub-processes that take place within the RIT including volute, nozzle, interspace, and rotor. The modeling of the RIT is based on a mean-line (one-dimensional) model [28] which assumes that the working fluid parameters change only with respect to mean streamline through the individual stage components.
In the preliminary design of the RIT, most attention is focused on the analysis of the flow through the rotor [29]. The cross-section of the rotor with mean streamline drawn between its inlet (4) and outlet (5) stations is depicted in Figure 2b. Furthermore, the velocity triangles at the rotor inlet (parameters with subscript 4) and outlet (parameters with subscript 5) are also presented. The common practice in designing the RIT is to assume the radial blades at the rotor inlet [30] which is equivalent to set the blade angle β 4,blade to 0 • (the positive sign of an angle is measured according to the positive direction Energies 2020, 13, 5280 5 of 27 of peripheral velocity u). For this reason, the relative flow angle and the incidence angle are equal (β 4 = i) since the latter is defined as i = β 4 − β 4,blade . To ensure no swirl flow at the rotor outlet, the absolute flow angle α 5 is assumed to be equal to 0 • [31]. Figure 2a,c, this type of heat exchanger geometry is applied to the vapor generator and condenser. Regarding the modeling assumptions, the vapor generator is divided with respect to the thermodynamic processes of the working fluid, including preheating, evaporation, and superheating. Evaporation is a two-phase process that is characterized by the rapid changes in the thermophysical properties of the fluid. For this reason, the evaporation is discretized into n sub-processes and the value of n = 10,000 is assumed as it provides high accuracy of the calculation procedure. Similarly, the condenser is divided into sections in which de-superheating and condensation processes take place. As the condensation is also a two-phase process, it is modeled in the same manner as the evaporation of the working fluid.
As a common practice in designing the ORC, the pump of a system is modeled as a device working with a constant efficiency [33].

Assumptions and Calculation Steps
The calculation model is built in a form of the originally developed code written in MATLAB environment [34]. The detailed description of the modeling is out of the scope of this study. Nevertheless, the individual calculation steps are briefly explained using a simplified flowchart  Because of the high heat transfer performance, flexible thermal sizing, and compact dimensions, plate heat exchangers proved to be particularly suitable for ORC power plants [32]. As depicted in Figure 2a,c, this type of heat exchanger geometry is applied to the vapor generator and condenser. Regarding the modeling assumptions, the vapor generator is divided with respect to the thermodynamic processes of the working fluid, including preheating, evaporation, and superheating. Evaporation is a two-phase process that is characterized by the rapid changes in the thermophysical properties of the fluid. For this reason, the evaporation is discretized into n sub-processes and the value of n = 10,000 is assumed as it provides high accuracy of the calculation procedure. Similarly, the condenser is divided into sections in which de-superheating and condensation processes take place. As the condensation is also a two-phase process, it is modeled in the same manner as the evaporation of the working fluid.
As a common practice in designing the ORC, the pump of a system is modeled as a device working with a constant efficiency [33].

Assumptions and Calculation Steps
The calculation model is built in a form of the originally developed code written in MATLAB environment [34]. The detailed description of the modeling is out of the scope of this study. Nevertheless, the individual calculation steps are briefly explained using a simplified flowchart presented in Figure 3.

Validation of the Calculation Model
The verification of the algorithm which is applied in this work was presented in [37]. By comparing the results generated by the developed procedure with those presented in the literature, it was shown that the presented calculation model leads to reliable findings.

Bi-Objective Optimization
It is possible to include different performance aspects simultaneously. Hence, a multi-objective approach has attracted a lot of attention in optimizing the operation of the ORC. As mentioned before, a bi-objective optimization has become especially popular in that field. The mathematical formulation The procedure starts by assigning values to the input variables including the volumetric flow rate of the geothermal water . V gw , the inlet temperature of the geothermal water T gw1 , the inlet temperature of the cold water T cw1, and thermophysical properties of the working fluid applied in the ORC. At this stage of the modeling, the decision variables are specified as well and they include the evaporation Energies 2020, 13, 5280 7 of 27 temperature T eva , the condensation temperature T con , the specific speed n s, and the temperature difference at the saturation liquid point ∆T eva (see Figure 1b). The first, the second, and the last parameter strongly affect the thermodynamic performance of the ORC. The third variable is closely related to the turbine efficiency and it is defined as follows: The detailed description of the decision variables and their limits can be found in [26]. Furthermore, the assumed operating parameters of the ORC are the same as those in [26] as listed in Table 1. To calculate basic thermodynamic parameters, including the mass flow rates or the temperature distributions of the fluids (geothermal water, working fluid, cold water), each component (or part of the component) of the ORC is considered as a control volume. Assuming that each process in the components (or subcomponents) is a steady-state, the mass, and energy balance equations can be written as follows: . m in = . m out (2) .
For determination of the working fluid properties (such as the specific enthalpy h or specific entropy s), REFPROP 9.0 [35] database is employed.
Using the parameters determined in the previous step, the radial-inflow turbine is designed according to the procedure partly discussed in [26]. The key output variable at this stage is a turbine total-to-static efficiency which is calculated iteratively using the following formula: The value of the η T given in Table 1 serves as an initial guess and it is updated until the convergence of the iterative procedure is ensured. The sum of the enthalpy losses Σ∆h loss involves exit velocity loss, incidence loss, passage loss, tip clearance loss, trailing edge loss, windage loss, and nozzle loss. The individual losses are calculated based on the empirical models available in the literature. While detailing such models is beyond the scope of this study, loss correlations similar to those in [36] were applied with detailed explanation for the specific enthalpies h 01 and h 5s given in [37].
In the next step, the thermodynamic and heat transfer analyses are conducted to determine the ORC performance indicators including the exergy efficiency η ex and total heat transfer area A tot . The η ex is calculated as [38]: The A tot is determined as the sum of the heat transfer areas of the vapor generator and condenser. To calculate the heat transfer areas corresponding to the individual sections of the vapor generator (preheating, evaporation, superheating) and condenser (de-superheating, condensation), the following formula is applied [39]: The detailed explanation concerning the indicators given in Equations (5) and (6) is presented in [40].
In the next stage, the optimization problem is formulated by defining the objective functions and constraints imposed on the decision variables (see Optimization problem section). Then, the optimal solutions can be found applying the Non-dominated Sorting Genetic Algorithm-II (NSGA-II) [41] which is a tool embedded in MATLAB. When the solutions are determined, they are displayed in a form of the Pareto front graph, a well-known graphical representation of the bi-objective optimization results. Additionally, the numerical values of the outcomes are given as the matrices of the objective functions F and decision variables DV, which completes the procedure.

Validation of the Calculation Model
The verification of the algorithm which is applied in this work was presented in [37]. By comparing the results generated by the developed procedure with those presented in the literature, it was shown that the presented calculation model leads to reliable findings.

Bi-Objective Optimization
It is possible to include different performance aspects simultaneously. Hence, a multi-objective approach has attracted a lot of attention in optimizing the operation of the ORC. As mentioned before, a bi-objective optimization has become especially popular in that field. The mathematical formulation of bi-objective optimization (BOO) can be written as follows: subject to The BOO can be described as minimizing the objective functions f 1 and f 2 (Equation (7)) over a vector → X, while satisfying the constraints expressed in Equations (8) and (9). The case in which the criterion is to be maximized can be treated by negating the objective function. The vector → X of the decision variables are expressed as: Energies 2020, 13, 5280 9 of 27 In the literature, there are many methods for dealing with multi-objective optimization (MOO) problems. As one of the most commonly applied optimization techniques [42], Non-dominated Sorting Genetic Algorithm-II (NSGA-II) is used in this study.

Optimization Problem
The optimization problem which was analyzed in [26] is recalled in this paper. Following the notation presented in Equations (7)-(9), the discussed BOO problem can be formulated as follows: subject to 60.0 3.00 The vector → X consists of the following decision variables: In Equation (11) the problem of maximizing the exergy efficiency η ex while minimizing the total heat transfer area A tot is formulated. The bounds for the decision variables are presented in Equations (12)- (15). The constraints providing a rational design of the radial-inflow turbine are expressed in Equations (16)- (22). The practical limitation imposed on the outlet temperature T gw2 of the geothermal water, which allows avoiding the silica oversaturation, is presented in Equation (23). A more detailed description of the optimization problem can be found in [26].

General Concept and Weighting Technique
The main feature of the proposed tool is the possibility of tracking the decisions taken by graphically indicating the extent to which the selected Pareto set points deviate from the ideal solution. For this reason, the developed software is called TRADeS (Tracking and Recognizing Alternative Design Solutions) which is developed in MATLAB (the program code is available under the link given in Supplementary Materials). The core of the developed procedure is the TOPSIS method which proved to be an effective technique for the decision-making process. The detailed description of the TOPSIS procedure was presented in [26].
It is worth noting that TOPSIS allows for specifying user preferences by setting weights to the criteria. The problem of assigning weights to the objectives is solved by applying a rating scale between 1.0 and 6.0. For this purpose, the weight superiority WS parameter is introduced. The WS is defined by assuming the linear relationship between the rating that is selected by the user and the weight that is assigned to the objective. Thereby, the weight w is determined using the following formula: where low and up are the lower (1.0) and upper (6.0) bounds of the adopted rating scale. The graphical illustration of the results from Equation (25) is presented in Figure 4. As seen, the WS = 1.0 is equivalent to assign w = 0.5, while by assigning WS = 6.0, the weight is equal to w = 1.0. In the case of the BOO, two criteria weights, w 1 and w 2 , must be specified. By specifying one of the weights, the second one is always determined using the following relationship: Energies 2020, 13, x FOR PEER REVIEW 10 of 28 solution. For this reason, the developed software is called TRADeS (Tracking and Recognizing Alternative Design Solutions) which is developed in MATLAB (the program code is available under the link given in Supplementary Materials). The core of the developed procedure is the TOPSIS method which proved to be an effective technique for the decision-making process. The detailed description of the TOPSIS procedure was presented in [26]. It is worth noting that TOPSIS allows for specifying user preferences by setting weights to the criteria. The problem of assigning weights to the objectives is solved by applying a rating scale between 1.0 and 6.0. For this purpose, the weight superiority WS parameter is introduced. The WS is defined by assuming the linear relationship between the rating that is selected by the user and the weight that is assigned to the objective. Thereby, the weight w is determined using the following formula: where low and up are the lower (1.0) and upper (6.0) bounds of the adopted rating scale. The graphical illustration of the results from Equation (25) is presented in Figure 4. As seen, the WS = 1.0 is equivalent to assign w = 0.5, while by assigning WS = 6.0, the weight is equal to w = 1.0. In the case of the BOO, two criteria weights, w1 and w2, must be specified. By specifying one of the weights, the second one is always determined using the following relationship: Therefore, assigning the WS to one of the objectives is equivalent to specifying to what extent the selected objective is superior to the other.

Description of the Procedure
The procedure (see Figure 5) starts from providing the matrices of the objective functions F and decision variables DV which are the outputs of the bi-objective optimization (see Figure 3). The user should also provide the information which objective function is the cost (to be minimized) and benefit (to be maximized) criterion by assigning the column numbers of the F to indexmin and indexmax parameters.
After running the program, two figures are displayed (see Figure 6a). The first one is a graph with the Pareto set points (Pareto front) which conventionally are called design points (DPs). To have a deeper insight into the individual DPs, the second figure shows a relative change between a certain Therefore, assigning the WS to one of the objectives is equivalent to specifying to what extent the selected objective is superior to the other.

Description of the Procedure
The procedure (see Figure 5) starts from providing the matrices of the objective functions F and decision variables DV which are the outputs of the bi-objective optimization (see Figure 3). The user should also provide the information which objective function is the cost (to be minimized) and benefit (to be maximized) criterion by assigning the column numbers of the F to index min and index max parameters.   After running the program, two figures are displayed (see Figure 6a). The first one is a graph with the Pareto set points (Pareto front) which conventionally are called design points (DPs). To have a deeper insight into the individual DPs, the second figure shows a relative change between a certain i-th DP and the ideal solution for the first (left y-axis) and the second (right y-axis) objective function. By referring to the examined optimization problem in which the objectives are the benefit (the exergy Energies 2020, 13, 5280 12 of 27 efficiency η ex ) and the cost (the total heat transfer area A tot ) criteria, the corresponding relative changes can be defined as follows: Energies 2020, 13, x FOR PEER REVIEW 13 of 28   By convention, the relative changes defined in Equations (27) and (28) will be named as deviations from the ideal solution for easier identifying their meaning and intended use. Having defined index min and index max , the program will automatically adjust the appropriate formula for the other configurations of the objective function types (cost or benefit).
Next, the user decides what range of decision-making preferences is to be covered (question 1). Technically, it is done by introducing the preference range PR parameter which is used in place of the WS in Equation (25). By assigning a certain value to the PR, the program applies Equation (25) twice to generate two-weight distributions, i.e., one that favors the first objective and the equivalent one that is beneficial for the second criterion. As an example, assigning PR = 5.0 is equivalent to specifying the criteria weights as: w 1 = 0.9, w 2 = 0.1 and w 1 = 0.1, w 2 = 0.9. By applying the TOPSIS procedure, the design points corresponding to the selected weight distributions are indicated. Additionally, the design point calculated by assigning w 1 = 0.5, w 2 = 0.5 is highlighted as being the characteristic Pareto solution reflecting equal importance of the criteria. The solutions which are distributed between the equal weights point and design points indicated using PR parameter are graphically indicated (see Figure 6b) as solutions favoring the first and the second objective. The discussed step of the procedure can be perceived as being equivalent to deciding what range of criteria weights are to be considered [24,25]. In other words, by assigning a certain value to the PR parameter, the set of design points corresponding to the selected criteria weights is indicated.
In the next step, the user decides which objective function is a priority criterion (question 2). By simply assigning the value (1 or 2) to the priority function PF parameter, one of the criteria (f 1 or f 2 ) is selected as primary. Then, the user must decide to what extent the final design point should differ from the point of the equal weight (question 3). This is done by direct assigning a value to the WS parameter. As depicted in Figure 5, the choice is limited to the range between 1.0 and PR. In other words, the user cannot choose the design point which is out of the range covered by assigning PR a few steps earlier. As an example, by assigning PF = 1 and WS = 2.0, the criteria weights are specified as: w 1 = 0.6, w 2 = 0.4. To highlight graphically (see Figure 6c) the proposed design point, the TOPSIS is applied for the second time. Thereby, the WS parameter is used (after choosing the priority criterion by applying the PF parameter) for indicating a single, and potentially the final design point.
The user can decide whether to consider other design points or mark the last choice as the final solution (see Figure 6d). To provide deeper insight into the decisions taken, the graph depicted in Figure 7 is displayed. It presents the design points corresponding to various weights w 1 which were specified indirectly by assigning values to the PF and WS parameters. Obviously, for a certain value of the w 1 , the weight distribution is fully defined since the second weight w 2 is calculated using Equation (26).   If the user opts for conducting the entire procedure from the beginning, the program will clear the marked points and the decision-making process will be repeated. If not, the procedure is completed. In the end, the program displays graphs with decision variable values corresponding to each design point. The decisions taken by the user are marked, indicating what values of the decision variables favor the examined objectives and what is the decision variable value corresponding to the final design point. As an example, the graph generated for one of the examined decision variables, i.e., the temperature difference ∆T eva (see Figure 1b), is presented in Figure 8. Being the fourth element of a vector written in Equation (24), it is marked with a number 4.

Results and Discussion
For a better understanding of the developed tool, an illustrative decision-making procedure applied for the results of bi-objective optimization is presented. As mentioned, the exergy efficiency ηex (objective function f1) and total heat transfer area Atot (objective function f2) were defined as the

Results and Discussion
For a better understanding of the developed tool, an illustrative decision-making procedure applied for the results of bi-objective optimization is presented. As mentioned, the exergy efficiency η ex (objective function f 1 ) and total heat transfer area A tot (objective function f 2 ) were defined as the criteria optimizing the operation of the ORC power plant. The set of optimal design points generated by applying NSGA-II is depicted in Figure 9 (obtained for the salinity of geothermal water equaling 0.00 kgkg −1 , for more details, see [26]). At this stage, the designer of the ORC may choose one of the Pareto set points and consider it as the final solution. The description of this process by analyzing the data generated by the TRADeS is given below.

Decision-Making Analysis
Defining the preference range of design points. It was decided to assign a value of 5.0 to the preference range PR parameter. In other words, the marked solutions (see Figure 9) are indicated for the criteria weights with their numerical values between 0.1 and 0.9 (see Table 2), similar to [25]. Based on this decision, the design points are divided into solutions favoring the first and the second criterion. The point corresponding to w 1 = 0.5, w 2 = 0.5 is indicated automatically and its role is to separate these solutions. The design points marked with the black dots are excluded from the analysis as solutions corresponding to weight distributions (WDs) which are out of the selected range. Examining weights' influence on the selected design points. The design points marked with the black stars are selected based on the different values assigned to the prior function PF and weight superiority WS parameters (see Table 2). By analyzing the graphs depicted in Figure 10 and the outcomes listed in Table 2, it is noted that both objective functions f 1 and f 2 are very sensitive to weights in a range between 0.4 and 0.6. As an example, by changing the weight w 1 from 0.4 to 0.6, the deviation (δf 1 ) id from the ideal solution of the exergy efficiency increased (towards (δf 1 ) id = 0%) by 36.3 pp, while by increasing the w 1 from 0.6 to 0.8, the (δf 1 ) id increased by only 10.2 pp. For a more detailed analysis of this relationship, the graph depicted in Figure 11 is examined. Specifically, the steepest slopes of the (δf 1 ) id and (δf 2 ) id are reported for the weights w 1 ranging between 0.4 and 0.6. Thereby, by increasing w 1 , within the aforementioned range, a significant improvement in the exergy efficiency η ex can be obtained. Similarly, by decreasing w 1 , from 0.6 to 0.4, the total heat transfer area A tot can be substantially reduced. Therefore, the decision-maker should keep in mind that the change of weights in the interval of 0.4 to 0.6 is associated with significant improvement of the first objective and substantial deterioration of the second one. For the remaining weight distributions (0.6 < w 1 ≤ 0.9 and 0.1 ≤ w 1 < 0.4), the objectives are less sensitive, or their numerical values do not change (see design point numbered as 16 in Figure 9).
From the user's perspective, the discussed sensitivity of the solutions to the criteria weights can be easily identified by tracking the locations of design points marked on graphs presented in Figures 9  and 10a. In particular, if the selected solutions are distributed close to each other, they are less sensitive to the assigned weights. For points located far from each other, an inverse relationship is observed, i.e., they are more sensitive to the selected weights.
Excluding adverse design point regions. Another decision-making aspect which can be examined by applying the proposed approach relates to the rate at which the individual solutions of the objectives f 1 and f 2 change for different regions of design points. To obtain a more complete picture of the problem, the results presented as points in Figure 10a are approximated by the polynomial curves (see Figure 10b) using the Curve Fitting toolbox embedded in MATLAB software. It was assessed that the regions I and III correspond to the substantial difference of rates at which the objectives approach their ideal solutions. The assessment of which Pareto set points belong to which specific region is subjective. However, the decision is not made unreasonably, but it is based on a careful analysis which is explained in the following paragraphs.   Energies 2020, 13, x FOR PEER REVIEW 21 of 28 Figure 11. Deviations from an ideal solution for different sets of weights. Figure 11. Deviations from an ideal solution for different sets of weights.
As presented in Table 3, for region I, the highest average rate of change of the (δf 1 ) id is combined with the lowest average rate of change of the (δf 2 ) id . By following the direction in which (δf 2 ) id Energies 2020, 13, 5280 19 of 27 approaches 0%, the above relationship means that a slight improvement of the objective f 2 is associated with a significant deterioration of the criterion f 1 . By examining the results covering region III, an inverse relationship is reported, i.e., when the (δf 1 ) id approaches 0%, a small beneficial increase in f 1 is combined with a substantial unfavorable increase in f 2 . For this reason, the design points located in regions I and III are excluded, limiting the choice of the decision-maker to solutions covering the region II. Thereby, the weights w 1 ≤ 0.3 and w 1 ≥ 0.7 (see Table 2), as values corresponding to design points which are placed in regions I and III, should not be considered for the selection of the final operating point. For further investigating the effects described above, the instantaneous rates at which the objectives approach their ideal values were determined. This was done by calculating the first derivatives of the deviations (δf 1 ) id and (δf 2 ) id (presented in Figure 10b) with the design point as a variable. As depicted in Figure 10c, regions I and III are indeed the areas in which an intensified contradicting relation between the objectives is observed. Specifically, for region I, the low rates at which the objective f 2 approaches the ideal value, are combined with the high rates at which f 1 moves away from its ideal solution. Similarly, for region III, the low rates at which the objective f 1 approaches the ideal value, are associated with the high rates at which the f 2 moves away from its ideal solution. On the other hand, the design points covering region II are characterized by the intermediate rates at which the objectives approach/move away from their ideal values. Thereby, the design points located in region II may be associated with the most sustainable relation between deterioration and improvement of the criteria, while approaching the ideal value of one of them.
Technically, regions I and III can be perceived as sets of Pareto solutions similar to Pareto weak optimal points [43]. Considering such points, it is possible to improve/aggravate one of the objectives without changing the value of the other. For these cases, the rate of change of the deviation from the ideal solution would be equal to 0.00% for one of the objectives. As depicted in Figure 10c, the rates of the (δf 1 ) id and (δf 2 ) id are close to 0.00%, being less than 0.55% for the first criterion (region III) and less than 0.30% for the second objective (region I).
Summarizing, it should be noted that a discussed relationship between the objectives, emerging in regions strongly favoring one of the criteria, occurs in the case of the Pareto front convex to the ideal solution [44]. For the non-convex Pareto curves, the opposite relation would be observed. Specifically, by approaching the ideal value of the first objective, a significant improvement of this objective would be accompanied by a slight deterioration of the second criterion. For such regions, selecting design points that are closer to the ideal value of the first objective would be preferable. In practical applications, the non-convex shapes of the Pareto curves are considered less frequently [45], thereby making the above-mentioned recommendations applicable for the ORCs and other industrial systems.
Selecting the final design point. After an initial overview of the design points, the designer can choose the potential final optimal solution. Since the size of the system is usually not as important as the power capacity for the geothermal ORC, it was decided to favor the first objective (PF = 1) with the WS parameter equaling 2.0 which is equivalent to the w 1 = 0.6, w 2 = 0.4 weight distribution.
Overview of the considered design points with respect to weights. By negating the selection of other design points, the user is directed to the next decision-making step. For deeper identifying the decisions taken, the results depicted in Figure 12 are examined. The considered design points are presented with respect to the selected weights w 1 . As in the case of the graph shown in Figure 11, it is seen that the objective functions are very sensitive to weights in a range between 0.4 and 0.6. Therefore, a detailed analysis of the impact of the weight on the objectives is provided, thus satisfying one of the requisites of the robust decision-making process [46].   Restarting/ending the decision-making process. At this stage, the designer can opt for reconsidering the weight distributions and thereby analyzing new solutions, including the final design point. For the studied case, it was decided not to consider other distributions of the potential design points, thus ending the decision-making process.
Analysis of the decisions taken with respect to decision variables. The graphs which complete the entire procedure are presented in Figure 13. The final decisions taken by the user are marked on scatter distribution graphs which are generated for each decision variable of the examined optimization problem.
Taking into consideration the bounds specified in Equations (12) and (14), it can be stated that the evaporation temperature T eva (see Figure 13a) and specific speed n s (see Figure 13c) are distributed in a narrow range. Specifically, the T eva values close to 90 • C are reported as optimal, while for the n s most outcomes are hovering around 0.55. For this reason, it is not possible to divide the obtained outcomes into values favoring the first (η ex ) or the second (A tot ) objective. In the case of the condensation temperature T con (see Figure 13b), it is clear that the values close to 40 • C (upper bound, see Equation (13)) favor the second criterion, while the values close to 25 • C (lower bound, see Equation (13)) are beneficial for the first objective. However, most of the solutions (favoring both f 1 and f 2 ) are obtained for 25 • C. As seen in Figure 13d, the temperature difference ∆T eva is monotonically decreasing with respect to design points and the values closer to 20 K (upper bound, see Equation (15)

Comparative Analysis with other Decision-Making Methods
The design point selected using the proposed procedure is marked along with solutions indicated by applying well-established decision-making approaches (see Figure 14). As mentioned in the Introduction, the LINMAP technique is based on the assumption that the Pareto set point which has the shortest distance from the ideal point is the best candidate for the final solution. Similarly, by applying the TOPSIS approach, the design point which has the shortest distance from the ideal point and the longest distance from the non-ideal point is considered as the definitive choice. As discussed previously, the scholars utilizing these techniques rarely took the opportunity to specify the decisionmaking "preferences" by assigning weights to the criteria. Having analyzed the LINMAP or TOPSIS calculation procedure, it is clear that omitting this step is equivalent to specifying an equal importance of the objectives by implicitly assigning w1 = 1.0 and w2 = 1.0. The design points indicated by applying this approach for the LINMAP and TOPSIS techniques are numbered as 53 and 60, respectively.
If one decides to specify the preferences with respect to the optimization criteria, one can opt for the subjective or objective methods for assigning weights. An application of the objective approach is illustrated in Fig 14. Specifically, the design point numbered as 101 is indicated by employing the TOPSIS technique combined with Shannon's entropy method. The latter was applied for determining the criteria weights. The detailed description of Shannon's entropy procedure can be found in [23]. With regard to the design point numbered as 83, it was determined by applying the developed Summarizing, the graphs depicted in Figure 13 are complementary to the decision-making procedure. Even though they do not affect the decision-making process itself, they provide information as to what values of the decision variables correspond to the final design point and what ranges of these parameters favor the individual objectives.

Comparative Analysis with other Decision-Making Methods
The design point selected using the proposed procedure is marked along with solutions indicated by applying well-established decision-making approaches (see Figure 14). As mentioned in the Introduction, the LINMAP technique is based on the assumption that the Pareto set point which has the shortest distance from the ideal point is the best candidate for the final solution. Similarly, by applying the TOPSIS approach, the design point which has the shortest distance from the ideal point and the longest distance from the non-ideal point is considered as the definitive choice. As discussed previously, the scholars utilizing these techniques rarely took the opportunity to specify the decision-making "preferences" by assigning weights to the criteria. Having analyzed the LINMAP or TOPSIS calculation procedure, it is clear that omitting this step is equivalent to specifying an equal importance of the objectives by implicitly assigning w 1 = 1.0 and w 2 = 1.0. The design points indicated by applying this approach for the LINMAP and TOPSIS techniques are numbered as 53 and 60, respectively.
The unquestionable advantage of the examined approaches (LINMAP, TOPSIS, Shannon's entropy + TOPSIS) is a simple structure and fast indication of the final design point for the examined system. These features make them particularly useful for inexperienced decision-makers. The tool developed in this study allows for conducting a more comprehensive decision-making process, providing the opportunity for specifying the decision-maker preferences, tracking the decisions taken and analyzing the results of the procedure.

Conclusions
This study aimed to propose a decision-making tool (called TRADeS) for a bi-objective optimization of the geothermal ORC power plant. By maximizing the exergy efficiency and minimizing the total heat transfer area, the Pareto set of optimal design points was determined. To conduct an in-depth decision-making study and select the final optimal solution, the developed tool was utilized. A thorough analysis which was performed by examining the data generated by created software allowed to draw the following conclusions: If one decides to specify the preferences with respect to the optimization criteria, one can opt for the subjective or objective methods for assigning weights. An application of the objective approach is illustrated in Figure 14. Specifically, the design point numbered as 101 is indicated by employing the TOPSIS technique combined with Shannon's entropy method. The latter was applied for determining the criteria weights. The detailed description of Shannon's entropy procedure can be found in [23]. With regard to the design point numbered as 83, it was determined by applying the developed decision-making tool. As already discussed, the proposed procedure is based on the application of the TOPSIS technique. However, the criteria weights are determined subjectively based on the detailed analysis described in the decision-making process. Thereby, the proposed approach is not solely based on the mathematically defined procedures, making the selection of the final operating point a more conscious and comprehensible process.
In the "Decision-making process" section, it was suggested that the final design point should be located in a range between 20 and 85 design point (see Figure 10b). As seen in Figure 14, the solutions obtained using the LINMAP and TOPSIS procedures are located inside the suggested region (in region II), whereas the design point indicated by applying the Shannon's entropy + TOPSIS approach is placed out of the acceptable range (in region III). The choice of design points belonging to the region III was not recommended since, in that range of solutions, an intensified adverse relationship between the objectives was reported.
The unquestionable advantage of the examined approaches (LINMAP, TOPSIS, Shannon's entropy + TOPSIS) is a simple structure and fast indication of the final design point for the examined system. These features make them particularly useful for inexperienced decision-makers. The tool developed in this study allows for conducting a more comprehensive decision-making process, providing the opportunity for specifying the decision-maker preferences, tracking the decisions taken and analyzing the results of the procedure.

Conclusions
This study aimed to propose a decision-making tool (called TRADeS) for a bi-objective optimization of the geothermal ORC power plant. By maximizing the exergy efficiency and minimizing the total heat transfer area, the Pareto set of optimal design points was determined. To conduct an in-depth decision-making study and select the final optimal solution, the developed tool was utilized. A thorough analysis which was performed by examining the data generated by created software allowed to draw the following conclusions: • while analyzing the rates at which the criteria approach their ideal solutions, it was found that there are regions of design points in which an intensified deterioration of one of the objectives is combined with a slight improvement of the second criterion. For this reason, from among the 105 design points, these numbered as lower than 20 and higher than 85 have been excluded from the set of potential final design points; • while considering the final design point, it was highlighted to discard the weights w 1 ≥ 0.7 and w 1 ≤ 0.3 as values corresponding to excluded Pareto set points; • most commonly applied forms of decision-making techniques, such as LINMAP or TOPSIS, are appropriate for inexperienced designers as they do not require specifying the decision-maker preferences. As a tool allowing for a careful examination of the criteria weights, the proposed program provides a more conscious and comprehensible choice of the final optimal configuration for the examined system.
The tool which was developed in this study should be perceived as a useful decision-making support for practical bi-objective optimization problems. By using the software and implementing the proposed recommendations, a more detailed analysis of the Pareto solutions is provided, giving a deeper insight into the decision-making process and selection of the final design point. Simultaneously, intuitive operation and generic structure allow for applying the tool in other industrial applications. Funding: The work has been partially financed by the Polish National Agency for Academic Exchange for implementation the project entitled "R&D in the field of increasing the efficiency of low-temperature geothermal power plants" under the agreement PPN/BEK/2018/1/00408/DEC/1.

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