Multi-Objective Optimization of Organic Rankine Cycle (ORC) for Tractor Waste Heat Recovery Based on Particle Swarm Optimization

: Diesel engines are widely used in agricultural tractors. During ﬁeld operations, the tractors operate at low speed and high load for a long time, the fuel efﬁciency is only about 15% to 35%, and the exhaust waste heat accounts for 38% to 45% of the energy released from the fuel. The use of tractor exhaust waste heat can effectively reduce fuel consumption and pollutant emissions, of which the organic Rankine cycle (ORC)-based waste heat recovery conversion efﬁciency is the highest. First, the diesel engine map is achieved through the test rig, a plate-ﬁn evaporator is trial-produced based on the tractor size, and the thermodynamic and economic performance model of the ORC are established. Then, taking the thermal efﬁciency of ORC and the speciﬁc investment cost (SIC) as the objective function, the particle swarm optimization (PSO) algorithm and the technique for order of preference by similarity to ideal solution (TOPSIS) decision method were used to obtain the optimal operating parameter set under all working conditions. Finally, the results showed that the ORC thermal efﬁciency could reach a maximum of 12.76% and the corresponding SIC value was 8539.66 $/kW; the ORC net output power could be up to 8.31 kW compared with the system without ORC; and the maximum brake speciﬁc fuel consumption (BSFC) could be reduced by 8.3%. The improvement in the thermodynamic performance will lead to a sacriﬁce in economic performance, and at high speeds, the economic beneﬁts and thermal efﬁciency reach a balance and show a better thermal economic performance. Recovering exhaust heat energy through ORC can reduce tractor fuel consumption and pollution emissions, which is one of the effective technical means to achieve “carbon neutrality” in agricultural production. At the same time, through the PSO algorithm, the optimal combination of ORC operating parameters is obtained, which ensures that the exhaust heat energy can be effectively recovered during the tractor ﬁeld operation, and provides a basis for the adjustment of real-time work strategies for future research.


Introduction
Since the 21st century, the demand for tractors has increased with the continuous development of agriculture in China. As a significant part of agricultural production and having the characteristics of high production efficiency and low unit cost, the number of tractors has increased from 15.80 million in 2004 to 23.17 million in 2020, and the total power of the country's agricultural machinery reached 1.056 billion kilowatts [1]. China's tractor market is a typical stock market, and the increase in tractor ownership and total agricultural machinery power indicates that China's agricultural production has shifted from relying mainly on human and animal power to mechanical power, entering a new period of mechanization. China is continuing to promote the full mechanization of agricultural production, which shows that agricultural modernization is in a critical period of changing the development mode, and the use and ownership of tractors will continue to increase in the coming period.
As agricultural machinery with the longest average annual working time (more than 500 h) [2], tractors are increasingly used in agricultural production, and the problem of their energy consumption and pollutant emissions is gradually attracting attention [3,4]. The diesel engine is widely applied in tractors and has accounted for more than 95% of the agricultural machinery power source; tractors operate at a constant speed most of the time with long working hours and heavy load, so a large amount of exhaust waste heat energy is lost directly to the air, failed to be effectively used, and has been recognized as a notable source of pollutant emissions, contributing to poor air quality and negative impacts on human health [5].
What calls for special attention is how to reuse the exhaust waste heat energy to improve the fuel efficiency and reduce emissions. In recent years, with the increasing environmental pollution and energy dependence, the traditional diesel engine energy saving, and emissions reduction have been of wide concern, but in the field of agricultural machinery, tractor diesel engine waste heat recovery (WHR) and energy reuse have not yet attracted sufficient attention. Diesel engine fuel efficiency during tractor operation in the field is only 15% to 35%, and exhaust energy accounts for 38% to 45% of the energy released from the fuel, so recycling the reuse of exhaust waste heat can help improve fuel utilization and reduce emissions. Research has shown that the WHR based on the organic Rankine cycle (ORC) has the highest exhaust heat energy conversion rate [6]. Punov et al. [7] presented the possibilities of exhaust heat recovery for a 110 kW tractor engine and the simulations results revealed that the engine output power and efficiency increased within the range of 3.9-7.5%. At full load, the calculated exhaust enthalpy was within the range of 32.1 kW to 103.4 kW, as the corresponding temperature varied from 693 K to 785 K. At full engine load, 8.36-14.8% of the fuel energy can be converted into mechanical work. ORC is considered as a promising heat-to-power technology to utilize waste heat and renewable energy. However, since the thermodynamic, economic, and environmental performance is usually in conflict, a single objective design can no longer meet the ORC system requirements, and the urgent requirements for multi-objective optimization have attracted increasing attention [8]. As the most widely used means of transportation today, automobiles have caused a lack of energy due to the rapid increase in their numbers. How to improve energy utilization to achieve energy saving and emissions reduction has attracted a series of scholarly research on automobiles. ORC technology has received extensive attention due to its superior performance. Although the current research on ORC for diesel engine WHR is now concentrated in the field of on-road vehicle engineering, the WHR of tractors, which also use diesel engines as a power source, can utilize similar research methods and technical means.
In the cases of off-highway vehicles, it is possible to observe stable high speed and torque profiles during operations for agricultural tractors. Lion et al. [9] assessed the possibility of applying an ORC system in a commercial agricultural tractor, recovering waste heat from a 300 kW power heavy-duty diesel engine. A maximum fuel consumption reduction of 10.6% was obtained using heat from the tailpipe and EGR. Different WHR technologies including ORC to be used in agricultural applications were analyzed [10]. More than 80% of the time, the engine operated at high load and almost constant speed when the tractors are run on the fields; at full engine load, the exhaust heat is lost to the environment in the form of power from 42.5 kW to 139.9 kW with the temperature varying from 700 K to 800 K. The exergetic analysis showed that a maximum of 15% of fuel low heating value can be reused by employing ORC, and at the most typical operation points, this value is from 10% to 12.8% [11].
The existing research is mainly on tractor steady-state operating conditions or for a single, stable heat source condition; there has been less analysis of tractors in off-design performance or variable operating conditions of the ORC operating state as well as the impact on system stability, where system optimization has not been achieved. Therefore, the focus and difficulty of the research have changed from the energy conversion efficiency of the ORC at a single, stable design operating point to how to dynamically adjust the ORC operating parameters and the synergistic matching with the exhaust heat energy to ensure the operational stability and effectiveness of the system under variable operating conditions [12,13]. A good design of an ORC system for tractor applications should be lightweight, a small size, low cost, evaporator pressure drop, high net power output, better off-design performance, and high controllability [14]. The most appropriate way to solve these problems is to perform multi-objective optimization. At present, the solutions in relevant research are roughly divided into two types. One is to predict through simulation, and the other is to optimize through intelligent algorithms combined with neural networks. In addition, due to the actual operation process, the collection of data such as flow, temperature, speed, and other data at each operating point consumes a lot of manpower, financial resources, and energy. Therefore, in the research, it is often necessary to process the data and conduct correlations based on the existing data. In recent years, due to the advantages of self-adaptation, self-learning, nonlinear mapping, and fault tolerance, neural networks with optimization algorithms such as particle swarm optimization (PSO) have proven to be an effective tool for design and performance optimization in the field of thermal engineering.
System optimization based on the PSO algorithm was performed to indicate the development potential of the ORC system, and the PSO was employed to optimize the parameters of pulsating flow when a diesel engine operates under the motorway road condition [15,16]. The multi-objective genetic algorithm (MOGA) was used to optimize three combined heat and power based on ORC systems, respectively, to achieve higher exergy efficiency and profit ratio of investment [17]. Additionally, the PSO was used to optimize the net output power as the objective function, and results revealed that the overall thermal and exergy efficiencies were obtained at 18.23% and 62.37%, respectively, for the base case [18]. A reliable and time-efficient optimization tool was developed in the MATLAB environment that was able to optimize ORC, the genetic algorithm (GA) toolbox was used, and the fluids' thermophysical properties were acquired from the CoolProp [19] and REFPROP [20] databases [21]. A machine learning prediction model was developed and applied to predict the ORC power output, GA was used to optimize the initial weights and thresholds of the different structural parameter ORC model to further improve the model generalization ability [22]. A simulation model of an ORC fin-and-tube evaporator via the computational fluid dynamic method was established, where the evaporator with a nonuniform structure obtained by GA optimization could further enhance the heat transfer and improve the flow state of the evaporator [23]. Based on machine learning, the GA and PSO were combined into a GA-PSO hybrid algorithm to predict and optimize the pump isentropic efficiency under full operating conditions in an ORC system [24]. The performances of GA, PSO, and repulsive particle swarm optimization (RPSO) based on the optimization of thermo-economic indicators were compared as well as in the optimization of the thermodynamic process of an ORC energy recovery system [25]. A multi-objective optimization using the PSO algorithm aimed at the proposed heat exchanger was conducted and the optimal design scheme had the potential to achieve a 71.46% decrease in the total annual cost [26]. A novel thermo-economic performance indicator, namely, maximum net power output with the constraint (MPC), was proposed, and the organic and steam Rankine cycle systems driven by the flue gas were optimized to maximize the net power output with the constraint of electricity production cost by using the non-dominated sorting genetic algorithm-II [27]. Performance prediction and multi-objective optimization using a back propagation neural network (BPNN) were investigated, an ORC experiment platform was used to obtain 950 sets of basic experimental data, the results demonstrated that the prediction error of the BP-ORC model was very low, and the system performance could be improved by adjusting several operating parameters experimentally according to the model prediction [28]. A comprehensive method to achieve a reasonable application of machine learning into ORC research for the prediction and optimization of the ORC's parameters and performance was presented. BPNN and support vector regression (SVR) prediction models for ORC were built by predicting error analysis with part of the database [29]. A preliminary design method was developed to optimize the eight critical parameters of the ORC based on the PSO algorithm; the isotropic efficiency was the objective function and the results indicated that the turbine efficiencies increased with the reduction in the pressure ratio and turbine inlet temperature and the increase in the power output. The exergy efficiency and isotropic efficiency dropped slightly with the increase in the turbine inlet temperature. Therefore, the turbine with R245fa exhibited good off-design performance [30]. PSO with parallel computation was adopted to optimize the performance of the nonlinear ORC system, where the result showed that the optimal system could output 6.87% more power than the initial system and 20.08% more power [31]. A fuzzy nonlinear dynamic model of the ORC evaporator using renowned finite volume (FV) was developed to capture the transient effects accurately, and the results showed that the fuzzy-based model was able to capture the transient effects at a data fitness of over 90% [32]. A control-oriented neurofuzzy model of brazed-plate evaporators for use in ORC was built to assess the dynamic performance of the ORC system. The proposed adaptive neuro-fuzzy inference system (ANFIS) model consisted of two separate neuro-fuzzy sub-models to predict the evaporator output temperature and evaporating pressure. The effect of training the models using gradient-descent least-square estimate (GD-LSE) and PSO techniques was investigated [33]. A PSO algorithm was used to optimize the operating parameters of the regenerative ORC system under various engine operating conditions [34]. A multi-objective optimization approach on the PSO for parametric optimization was applied to find a set of Pareto optimal solutions [35]. A multi-condition performance prediction method with the PSO algorithm for ORC centrifugal pumps was proposed and compared with the experimental results, where the performance under multiple operating conditions was well-predicted by introducing performance constraints [36]. To reduce the computation time, a new fuzzybased evaporator model was developed and the performance compared with the finite volume method. The results showed that the fuzzy technique can be applied to predict the output of the supercritical evaporator in the WHR system and can significantly reduce the required computation time [37]. Considering the great potential of deep reinforcement learning (DRL) for solving complex control problems, a DRL-based energy management system (EMS) was proposed for an HEV-ORC. The simulation results demonstrated that the DRL-based EMS could save 2% more fuel energy than the rule-based EMS [38]. A nonlinear model predictive controller (NMPC) was designed to provide optimal control input for maximum turbine power generation in an ORC-WHR system. The power optimization-based NMPC utilized an extended Kalman filter (EKF) as a state estimator. The working fluid evaporation pressure was controlled by an external PID control loop. The experimental study showed that the augmented control scheme outperformed the baseline multi-loop PID control in terms of both the tracking error and settling time during transient engine operation [39,40]. Most of the existing research in the literature on machine learning has been to collect single indicators such as the ORC thermal efficiency with different ORC operating parameters including the working fluid mass flow rate, pressure, etc. at a single engine operating condition, then build a neural network to predict and find the optimal value in a given range of ORC operating parameters, which means the amount of data needed to be collected and trained for each of the engine operating conditions. However, this research method often only considers a single indicator and not all of the engine operating conditions are available in all cases.
The current research literature on exhaust waste heat recovery based on ORC has mainly focused on the technical economy and thermodynamic performance of the systems under engine on-design operating conditions by means of experimental data and simulations. However, in the operation of the tractor, there are differences in the engine speed and load under the switching of different working modes such as field tillage, headland steering, and idle parking, which lead to system instability due to low degree superheating of the working fluid at the evaporator outlet. Therefore, the study of the tractor exhaust heat transfer law with the ORC combined and how to determine the operating range and predictive the performance of the ORC parameters under off-design conditions was the research question of this paper. In order to ensure the efficient and stable operation of the system, the parameters of the working fluid and coolant water such as the pressure, degree of superheat and subcooling, and mass flow rate were selected as the optimization parameters, and the thermal efficiency of ORC and the economic performance indicators were used as the objective functions. Compared to the several implementations presented with machine learning, the PSO algorithm combined the technique for order of preference by similarity to ideal solution (TOPSIS) decision-making were adopted to solve the Pareto frontier and optimize the parameter combinations, respectively. The optimal parameter combinations and distributions were determined under all operating conditions and the validity of the model was verified by comparing the experimental data with the PSO optimization model and the prediction results of the BPNN. The rest of this paper is organized as follows. In Section 2, the universal characteristic data of the diesel engine are obtained through the test rig, the detailed thermodynamic and economic performance model of each component of the ORC are established, and a multi-objective optimization model is built on this basis. In Section 3, single and multi-operating conditions optimization and thermo-economic optimization are analyzed. The conclusions are presented in Section 4.

System Description
The aim of this research was to study the exhaust heat potential of the engine under all operating conditions, and how the thermodynamic parameters of the ORC working fluid are adapted to the engine operating conditions and exhaust parameters to ensure that the exhaust heat energy can be effectively utilized. Therefore, it is necessary to measure the brake power, brake specific fuel consumption (BSFC), exhaust temperature, pressure, and mass flow rate of the engine under all operating conditions. In this paper, an eddy current dynamometer was used to measure the 4-cylinder, 4-stroke, 1.9 L turbocharged diesel engine brake power, torque, and speed; a NI 4-slot, Ethernet chassis cDAQ-9184 combined with four modules was employed to build a test data acquisition system. The exhaust data of the diesel engine at 210 operating points were collected. Figure 1 shows the brake power and exhaust temperature of the diesel engine under all operating conditions, the maximum torque and rated power were 361.4 N-m at 2250 rpm and 118.7 kW at 3500 rpm, respectively, the exhaust temperature increased with the engine speed and load, and the highest exhaust temperature was 950 K with the engine speed of 4000 rpm at full load. According to the variation law of the exhaust energy of a vehicle diesel engine, a set of the ORC system was designed, the diesel engine was used as the top cycle, the working cycle of the ORC system was used as the bottom cycle, and the combined system of the engine-ORC was constructed. The conceptual schematic of the engine-ORC system is shown in Figure 2, where R245fa was selected as the organic working fluid for circulation. As a key component of the ORC, the analysis of the thermal performance of the evaporator  According to the variation law of the exhaust energy of a vehicle diesel engine, a set of the ORC system was designed, the diesel engine was used as the top cycle, the working cycle of the ORC system was used as the bottom cycle, and the combined system of the engine-ORC was constructed. The conceptual schematic of the engine-ORC system is shown in Figure 2, where R245fa was selected as the organic working fluid for circulation. As a key component of the ORC, the analysis of the thermal performance of the evaporator under the space constrained condition of the tractor can provide a theoretical basis for the optimal design of the evaporator parameters as well as effectively improve the utilization of the exhaust heat of the tractor. Therefore, in this paper, after the actual measurements of the front space of the hood and the dimensions of the cab roof of three different horsepower tractors such as Dongfanghong 804/854/954, a plate fin evaporator was made on a trial basis. The experimental evaporator was constructed from aluminum heat-sink stock, which consists of three separate heat exchangers; the evaporator structural sketch and dimensional parameters are shown in Figure 3 and Table 1. It consisted of three separate heat exchangers, the thicker middle section as the hot-side for the engine exhaust and the two outer sections as the cold-side designed for the ORC working fluid flow. The thicker middle section was constructed using two identical sections of heat-sink material, the fins of one side, while facing the fins of other, were inserted into the gaps between the fins of the opposing side, resulting in a fin spacing one-half that of the heat sink material. The fin height is therefore equal to the fin height of the original heat sink stock.         Process 1-4: The working fluid starts a cycle as the liquid leaving the pump at state point 1, then is preheated to saturation and evaporated through the evaporator. The heat absorbed from the exhaust gases is calculated in Equation (1), the specific enthalpy i, and the temperature of each state point can be calculated from the energy balance equation.
where r m  and e m  are the mass flow rate of R245fa and the exhaust gases, respectively.
The evaporator, as a counter-flow finned heat exchanger, was modeled using log mean temperature difference (LMTD), considering the three-zone method, as shown in Figure 5. The thermodynamic evaporator model was developed for this paper to establish the thermodynamic model of the evaporator, and the following assumptions were made to simplify the numerical calculation procedure.
(1) Ignore the heat dissipated by the evaporator to the environment; (2) Ignore evaporator fouling factor and material thermal resistance; (3) Ignore the changes in heat conduction, potential energy and kinetic energy of the evaporator fluid along the axis of the flow channel and pipeline. Process 1-4: The working fluid starts a cycle as the liquid leaving the pump at state point 1, then is preheated to saturation and evaporated through the evaporator. The heat absorbed from the exhaust gases is calculated in Equation (1), the specific enthalpy i, and the temperature of each state point can be calculated from the energy balance equation.
where . m r and . m e are the mass flow rate of R245fa and the exhaust gases, respectively. The evaporator, as a counter-flow finned heat exchanger, was modeled using log mean temperature difference (LMTD), considering the three-zone method, as shown in Figure 5. The thermodynamic evaporator model was developed for this paper to establish the thermodynamic model of the evaporator, and the following assumptions were made to simplify the numerical calculation procedure.
(1) Ignore the heat dissipated by the evaporator to the environment; (2) Ignore evaporator fouling factor and material thermal resistance; (3) Ignore the changes in heat conduction, potential energy and kinetic energy of the evaporator fluid along the axis of the flow channel and pipeline. the thermodynamic model of the evaporator, and the following assumptions were made to simplify the numerical calculation procedure.
(1) Ignore the heat dissipated by the evaporator to the environment; (2) Ignore evaporator fouling factor and material thermal resistance; (3) Ignore the changes in heat conduction, potential energy and kinetic energy of the evaporator fluid along the axis of the flow channel and pipeline. The thermodynamic properties of the R245fa refrigerant and diesel exhaust gases such as the specific heat and the convection heat transfer coefficient were sourced from an online property calculator developed and maintained by the National Institute of Standards and Technology (NIST) [20]. Thermodynamic properties of the diesel engine exhaust gases were derived from known air and flue mass flow rates and assuming the complete combustion of the fuel having assumed the molecular structure of dodecane (C12H26).
The cross-sectional area of the evaporator pipe was circular, while the cross-sectional area of the flow channel was rectangular. Considering the convective heat transfer in the The thermodynamic properties of the R245fa refrigerant and diesel exhaust gases such as the specific heat and the convection heat transfer coefficient were sourced from an online property calculator developed and maintained by the National Institute of Standards and Technology (NIST) [20]. Thermodynamic properties of the diesel engine exhaust gases were derived from known air and flue mass flow rates and assuming the complete combustion of the fuel having assumed the molecular structure of dodecane (C 12 H 26 ).
The cross-sectional area of the evaporator pipe was circular, while the cross-sectional area of the flow channel was rectangular. Considering the convective heat transfer in the tube with different cross-sectional shapes, the Nusselt number Nu of R245fa and exhaust gas was calculated by Equation (2) [41]. The respective heat transfer coefficients h r and h e of R245fa and exhaust gas, the total heat transfer coefficient U of the evaporator, and the heat transfer area A can be calculated from the structural parameters and the LMTD method, as shown in Equations (3)-(5), where the two-phase heat transfer coefficient h ev is determined by Equation (6).
where a 0 is the scaling factor, 8.24 for the exhaust gases, and 5.39 for the working fluid, respectively. Re is the Reynolds number based on hydraulic diameter, Pr is the Prandtl number, k is the thermal conductivity of the fluid, D is the hydraulic diameter of the flow channel calculated by the parameters in Table 1, F and S are the forced convective heat transfer enhancement factor and suppression factor, respectively, h pool is the pooling boiling transfer coefficient [42]. Process 4-5: Saturated vapor from the evaporator at state 4, having an elevated temperature and pressure, expands through the expander to produce work and is then discharged to the condenser at state 5 with relatively low pressure. Neglecting heat transfer with the surroundings, the mass and energy rate balances for a control volume around the expander reduce at the steady state to give: where η exp,iso is the isentropic efficiency of the expander. Process 5-8: In the condenser, there is heat transfer from the working fluid to the cooling water. The working fluid condenses and the temperature of the cooling water increases; heat emitted from the work fluid to the cooling water is expressed in Equation (8). The modeling approach for the thermodynamic model of the condenser is similar to that of the evaporator, while the overall structure is half the size of the evaporator.
Process 8-1: The working fluid condensate leaving the condenser at 8 is pumped from the condenser into the higher-pressure evaporator. Taking a control volume around the pump and assuming no heat transfer with the surroundings, the mass and energy rate balances give where η pump,iso is the isentropic efficiency of the pump. The thermal efficiency gauges the extent to which the energy input to the working fluid passing through the evaporator is converted to the net output power. Using the quantities and expressions introduced above, the thermal efficiency η ORC of the ORC system is

Economic Model
The thermodynamic model of the ORC described in the previous section was expanded by considering the costs of the main elements in the cycle. These costs were estimated from the work of Lecompte et al. [43] and Quoilin et al. [44], and are presented in Equation (11). The total cost (TC) of the ORC system is mainly dependent on the component costs, C cp , and labor cost, C lab , which is half of the C cp . The cost calculation of common components should include the ORC heat exchanger, expander, pump, electric motor, and working fluid, and refer to the NETL (National Energy Technology Laboratory) and the results provided by Pantaleo [45] TC = C cp + C lab (11) The cost correlation is based on the heat exchanger heat transfer area A: The heat exchanger cost, C hx , is determined from: where where F p is the pressure factor and F m is the material factor, F m is set to 2.
As screw compressors can be modified to operate in reverse mode as expanders, the correlation from Astolfi [46] was used, where the size parameter is the volumetric flow rate at the outlet of the expansion, and the expander cost C exp is expressed in Equation (15). The pump cost expression C pump is shown in Equation (16) W p 300 0.25 (16) Energies 2022, 15, 6720 10 of 24 The cost of the electric motor that drives the pump is C motor : The working fluid cost C wf is determined by the mass: Therefore, the estimate TC is equal to the result of the component costs and the labor costs, and the specific investment cost (SIC) is used as the thermo-economic objective function.

Multi-Objective Optimization Model
Based on the above analysis, the size of the diesel exhaust heat is mainly determined by the operating conditions, and the distribution characteristics of the exhaust energy will change with the engine operating conditions. Since the ORC system involves multiple operating parameters and the coupling between each parameter, it is difficult to analyze the performance of the ORC under variable engine operating conditions and the relationship between the operating parameters simply through the sensitivity analysis. Therefore, in this paper, a multi-objective optimization model was established to optimize the main operating parameters of the ORC system, and the solution was solved by the PSO algorithm. Figure 6 shows the flow chart of the PSO algorithm established in this paper.
The particle swarm algorithm simulates a bird in a flock by designing a massless particle with only two properties: velocity and position, where velocity represents the speed of movement and position represents the direction of movement. Each particle searches for the optimal solution independently in the search space, which is recorded as the current individual extreme value, and shares the individual extreme value with other particles in the entire particle swarm. All particles in the particle swarm find the current individual extreme value according to their own. The current global optimal solution shared with the entire particle swarm adjusts its speed and position.
In this paper, the ORC thermal efficiency η ORC and SIC were used as the objective function, and the evaporation pressure, condensation pressure, mass flow rate, degree of superheating and subcooling of the working fluid, and the mass flow rate of the coolant water were selected as the decision variables. The multi-objective optimization mode can be expressed as: change with the engine operating conditions. Since the ORC system involves multiple operating parameters and the coupling between each parameter, it is difficult to analyze the performance of the ORC under variable engine operating conditions and the relationship between the operating parameters simply through the sensitivity analysis. Therefore, in this paper, a multi-objective optimization model was established to optimize the main operating parameters of the ORC system, and the solution was solved by the PSO algorithm. Figure 6 shows the flow chart of the PSO algorithm established in this paper.  The particle swarm algorithm simulates a bird in a flock by designing a massless particle with only two properties: velocity and position, where velocity represents the speed of movement and position represents the direction of movement. Each particle The value ranges of the optimization variables were set as shown in Table 2.  Since the objectives in multi-objective optimization are often mutually constrained, optimizing one of the objectives must be performed at the expense of degrading the performance of the other objectives. The difference with single-objective optimization is that instead of a unique optimal solution, there is an optimal solution set, which is explained by a concept in economics-the Pareto solution set. The frontier formed by the Pareto solution is also called the Pareto frontier. For practical problems, it is often necessary to pick one or more solutions from the Pareto solution set as the optimal solution to the problem to be solved. The setting of the algorithm parameters of the optimization model has a great influence on the optimization results. When applying the particle swarm algorithm for solving multi-objective optimization problems, the optimized model needs to be calibrated. The algorithm parameters involved mainly include the population size, number of iterations, mutation probability, learning factor, etc. Taking the rated working condition of the diesel engine as an example, to verify the optimization model, the optimal parameter setting should make the Pareto front obtained by the multi-objective optimization evenly distributed in a small area. Figure 7 shows the Pareto frontier results of different population numbers. The population numbers were set as 10, 20, 30, and 40, respectively. It can be seen that the population number had a greater impact on the optimization results. The Pareto frontier gradually dispersed with an increasing population number, but the thermal efficiency converged to about 14%, and the SIC value gradually increased. Therefore, the population size was set to 30. The Pareto frontier with different iterations is shown in Figure 8, and the maximum iterations were set to 30, 300, 500, and 1000, respectively. When the maximum number of iterations was 30, it had not converged enough to form a Pareto frontier, because the number of iterations was too small and the position and velocity of the particles had not been updated to the optimal position. The increase in the number of iterations will lead to an increase in the running time of the program, and does not change the shape of the Pareto frontier, so the maximum number of iterations was set to 300. The Pareto frontier for different archive sizes are shown in Figure 9. The archive sizes were set to 20, 200, 300, and 500, respectively. The archive size was used to save the best particle situation, similar to the maximum number of iterations, and the shape of the Pareto frontier was basically unchanged after exceeding 200, and the archive size was set to 200. Therefore, the specific parameter settings are shown in Table 3. the maximum iterations were set to 30, 300, 500, and 1000, respectively. When the maximum number of iterations was 30, it had not converged enough to form a Pareto frontier, because the number of iterations was too small and the position and velocity of the particles had not been updated to the optimal position. The increase in the number of iterations will lead to an increase in the running time of the program, and does not change the shape of the Pareto frontier, so the maximum number of iterations was set to 300. The Pareto frontier for different archive sizes are shown in Figure 9. The archive sizes were set to 20, 200, 300, and 500, respectively. The archive size was used to save the best particle situation, similar to the maximum number of iterations, and the shape of the Pareto frontier was basically unchanged after exceeding 200, and the archive size was set to 200. Therefore, the specific parameter settings are shown in Table 3.

Multi-Objective Optimization Decision Method
Unlike single-objective optimization problems, the solution set of multi-objective optimization problems often does not have an optimal solution, and the solution set has a Pareto frontier. The multi-objective optimization objectives often interact with each other or even contradict each other, and the decision-making process of selecting an optimal solution from the Pareto frontier according to the actual situation is very complicated. The results studied by Basílio [47] showed, in the multi-objective decision approach, that the analytic hierarchy process(AHP) and TOPSIS approaches had the most publications in the four decades (1977-2022) and were also the most commonly employed methods in solving multi-objective related issues. Computer science stands out among others because 47% of the methods address issues related to these areas, with the TOPSIS method being used the most. Engineering follows, with 35% of the methods, with the AHP method being the second most commonly used method. Therefore, TOPSIS was used to make a decision on the multi-objective optimization.
The basic principle of TOPSIS is based on the concept that the chosen alternative should have the shortest geometric distance from the positive ideal solution and the longest geometric distance from the negative ideal solution. If the evaluation object is closest to the positive ideal solution and farthest away from the negative ideal solution at the same time, it is the best, otherwise it is not optimal. Due to the different dimensions, it is often necessary to perform dimensionless processing on the optimization parameters. In this paper, the maximum and minimum method was used to perform the dimensionless processing of the decision-making objective, and the TOPSIS method was used to perform the decision processing on the Pareto solution set.
The distance between the target alternative and the positive ideal solution and the distance between the target alternative and the negative ideal solution were calculated as shown in Equations (23)-(24), respectively.

Multi-Objective Optimization Decision Method
Unlike single-objective optimization problems, the solution set of multi-objective optimization problems often does not have an optimal solution, and the solution set has a Pareto frontier. The multi-objective optimization objectives often interact with each other or even contradict each other, and the decision-making process of selecting an optimal solution from the Pareto frontier according to the actual situation is very complicated. The results studied by Basílio [47] showed, in the multi-objective decision approach, that the analytic hierarchy process(AHP) and TOPSIS approaches had the most publications in the four decades (1977-2022) and were also the most commonly employed methods in solving multi-objective related issues. Computer science stands out among others because 47% of the methods address issues related to these areas, with the TOPSIS method being used the most. Engineering follows, with 35% of the methods, with the AHP method being the second most commonly used method. Therefore, TOPSIS was used to make a decision on the multi-objective optimization.
The basic principle of TOPSIS is based on the concept that the chosen alternative should have the shortest geometric distance from the positive ideal solution and the longest geometric distance from the negative ideal solution. If the evaluation object is closest to the positive ideal solution and farthest away from the negative ideal solution at the same time, it is the best, otherwise it is not optimal. Due to the different dimensions, it is often necessary to perform dimensionless processing on the optimization parameters. In this paper, the maximum and minimum method was used to perform the dimensionless processing of the decision-making objective, and the TOPSIS method was used to perform the decision processing on the Pareto solution set. The distance between the target alternative and the positive ideal solution and the distance between the target alternative and the negative ideal solution were calculated as shown in Equations (23) and (24), respectively.
Equation (25) calculates the closeness of the optimization result to the positive ideal solution. The higher the value of C i , the closer it is to the positive ideal solution.
As shown in Figure 10, it describes how to find the optimal solution in a Pareto frontier. Since the goal in this paper was to have a greater efficiency and a smaller SIC value, the point in the lower right corner was the positive ideal solution and the upper left corner was the negative ideal solution. The closeness of each Pareto solution to the ideal solution was calculated through the equations, and the closest one was the optimal solution determined by the decision.
d + d (25) As shown in Figure 10, it describes how to find the optimal solution in a Pareto frontier. Since the goal in this paper was to have a greater efficiency and a smaller SIC value, the point in the lower right corner was the positive ideal solution and the upper left corner was the negative ideal solution. The closeness of each Pareto solution to the ideal solution was calculated through the equations, and the closest one was the optimal solution determined by the decision.  Figure 11 shows the effect of the ORC operating parameters on the thermal efficiency and SIC with the engine speed of 2250 rpm at full load. From Figure 11a, it can be seen that the thermal efficiency of ORC and SIC increased with the increase on the evaporating pressure. This is mainly because the enthalpy of the inlet of the expander increases with the increase in the evaporating pressure, which leads to an increase in the work capacity of the expander, so both the heat transfer area and the net output power increase, making the thermal efficiency and SIC value increase. In Figure 11b, the net ORC output power decreased as the condensing pressure increased, resulting in a continuous decrease in the thermal efficiency while the SIC value increased. In Figure 11c, the degree of superheating had little effect on the thermal efficiency of the cycle, but has a greater impact on the SIC value. The reason for the analysis is that the degree of superheating is a physical quantity used to ensure that the working fluid is in a superheated state at the inlet of the expander. In Figure 11d, as the degree of subcooling increased, the thermal efficiency decreased and the SIC value increased. This is because the temperature of the working fluid at the inlet of the pump decreased, and when the evaporating and condensing pressures remained unchanged, the power consumed by the pump during the isentropic process increased,  Figure 11 shows the effect of the ORC operating parameters on the thermal efficiency and SIC with the engine speed of 2250 rpm at full load. From Figure 11a, it can be seen that the thermal efficiency of ORC and SIC increased with the increase on the evaporating pressure. This is mainly because the enthalpy of the inlet of the expander increases with the increase in the evaporating pressure, which leads to an increase in the work capacity of the expander, so both the heat transfer area and the net output power increase, making the thermal efficiency and SIC value increase. In Figure 11b, the net ORC output power decreased as the condensing pressure increased, resulting in a continuous decrease in the thermal efficiency while the SIC value increased. In Figure 11c, the degree of superheating had little effect on the thermal efficiency of the cycle, but has a greater impact on the SIC value. The reason for the analysis is that the degree of superheating is a physical quantity used to ensure that the working fluid is in a superheated state at the inlet of the expander. In Figure 11d, as the degree of subcooling increased, the thermal efficiency decreased and the SIC value increased. This is because the temperature of the working fluid at the inlet of the pump decreased, and when the evaporating and condensing pressures remained unchanged, the power consumed by the pump during the isentropic process increased, resulting in a decrease in the net output power and thermal efficiency, and the SIC value increased.

Multi-Operating Conditions Optimization Results Analysis
In this paper, PSO was used to optimize the thermal efficiency and economic index SIC of the ORC system as the objective function under the operating conditions of the diesel engine from 850 rpm to 4500 rpm with 25%, 50%, 75% and 100% load, respectively. The optimization calculation results are presented in Figure 12. Figure 12a shows the variation in the optimal values of the evaporating and condensation pressure of the working fluid with the diesel engine operating conditions. The optimized evaporating pressure decreased with the increase in the torque and engine speed, the operating conditions had little effect on the condensing pressure and its value basically kept the minimum value of about 0.1 MPa within the given range, while the evaporating pressure ranged from 1.4 MPa to 3 MPa. Higher evaporating pressure will increase the power consumed by the pump, so in order to ensure the optimal SIC value target, when the output net power does not change much, the cost can be reduced by decreasing the pump power consumption. Figure 12b shows the variation in the degree of superheating and subcooling with engine speed and torque for the optimized ORC system. The optimized superheat range was between 1 and 10 K and the subcooling range was between 0 and 4.5 K under variable engine

Multi-Operating Conditions Optimization Results Analysis
In this paper, PSO was used to optimize the thermal efficiency and economic index SIC of the ORC system as the objective function under the operating conditions of the diesel engine from 850 rpm to 4500 rpm with 25%, 50%, 75% and 100% load, respectively. The optimization calculation results are presented in Figure 12. Figure 12a shows the variation in the optimal values of the evaporating and condensation pressure of the working fluid with the diesel engine operating conditions. The optimized evaporating pressure decreased with the increase in the torque and engine speed, the operating conditions had little effect on the condensing pressure and its value basically kept the minimum value of about 0.1 MPa within the given range, while the evaporating pressure ranged from 1.4 MPa to 3 MPa. Higher evaporating pressure will increase the power consumed by the pump, so in order to ensure the optimal SIC value target, when the output net power does not change much, the cost can be reduced by decreasing the pump power consumption. Figure 12b shows the variation in the degree of superheating and subcooling with engine speed and torque for the optimized ORC system. The optimized superheat range was between 1 and 10 K and the subcooling range was between 0 and 4.5 K under variable engine operating conditions. The graph shows that in most cases, the degree of superheating was the maximum value in a given range and the subcooling was the minimum value. This is because a higher superheat makes the working fluid at the evaporator outlet have a larger enthalpy, which makes the expander work harder, while a smaller subcooling makes the pump consume less power. Figure 12c shows the mass flow rate of the working fluid increased with the increase in the torque and engine speed. The exhaust mass flow rate was small under low and medium engine speed and torque, and obtained between 0 and 0.05 kg/s; when the engine speed was at 4500 rpm, the mass flow rate of the working fluid could reach a maximum of 0.23 kg/s. In Figure 12d, the variation law of the net output power of the ORC system was similar to the brake power of the diesel engine, and both increased with the increase in the engine torque and speed. This is mainly because the exhaust temperature and mass flow rate increase with the increase in the diesel engine speed and torque, so the exhaust energy that can be utilized by the ORC system also increases and leads to the net output power increases accordingly. When the engine speed was 4500 rpm with a full load at 216.7 N·m, the net output power of the ORC system reached the maximum power of 8.31 kW. operating conditions. The graph shows that in most cases, the degree of superheating was the maximum value in a given range and the subcooling was the minimum value. This is because a higher superheat makes the working fluid at the evaporator outlet have a larger enthalpy, which makes the expander work harder, while a smaller subcooling makes the pump consume less power. Figure 12c shows the mass flow rate of the working fluid increased with the increase in the torque and engine speed. The exhaust mass flow rate was small under low and medium engine speed and torque, and obtained between 0 and 0.05 kg/s; when the engine speed was at 4500 rpm, the mass flow rate of the working fluid could reach a maximum of 0.23 kg/s. In Figure 12d, the variation law of the net output power of the ORC system was similar to the brake power of the diesel engine, and both increased with the increase in the engine torque and speed. This is mainly because the exhaust temperature and mass flow rate increase with the increase in the diesel engine speed and torque, so the exhaust energy that can be utilized by the ORC system also increases and leads to the net output power increases accordingly. When the engine speed was 4500 rpm with a full load at 216.7 N·m, the net output power of the ORC system reached the maximum power of 8.31 kW.  In Figure 13a, the thermal efficiency of the ORC system increased with the increase in the engine speed and torque, the optimized range was from 10.72% to 12.76% and could reach the maximum efficiency at full load and an engine speed of 3500 rpm. In Figure 13b, the SIC value decreased with the increase in speed and torque. This is because the net output power increased with an increase in the speed and torque, and from Equation (20), the SIC value decreased. In Figure 13a, the thermal efficiency of the ORC system increased with the increase in the engine speed and torque, the optimized range was from 10.72% to 12.76% and could reach the maximum efficiency at full load and an engine speed of 3500 rpm. In Figure 13b, the SIC value decreased with the increase in speed and torque. This is because the net output power increased with an increase in the speed and torque, and from Equation (20), the SIC value decreased.

Thermo-Economic Optimization Results Analysis
In the process of TOPSIS multi-objective decision-making, when the optimal situation is selected under a single operating condition, it is actually a balanced state, and one of the objectives is often not the optimal situation. In Figure 10, the thermal efficiency of the optimal solution did not reach the maximum value, which is the ideal solution in the figure, and the SIC was not taken to the minimum value, which is the worst solution, that can be seen as a loss of part of the thermal efficiency to improve the economic efficiency in this case. Taking the data shown in Figure 10, considering the maximum thermal efficiency under this operating condition and the maximum SIC as the initial situation, the optimal solution obtained by the decision can be regarded as reducing the thermal efficiency by 7.93%, but increasing the thermal efficiency by 16.72% in terms of economic benefits. Figure 14 shows the ORC thermal efficiency loss rate and economic efficiency improvement rate of the optimization decision results under all of the operating conditions. Figure 14a represents the loss rate of thermal efficiency relative to the optimal solution under this operating condition, and Figure 14b indicates the economic benefit relative to the worst solution at this operating condition. The improvement rate in the economic benefits gradually decreased with the increase in the engine speed and torque. In the medium and low speed and load ranges, reducing the thermal efficiency can appropriately improve the economic benefits. When the engine was 1250 rpm and the load rate was 50%, the largest thermal efficiency loss rate was 0.12, and the corresponding economic improvement rate also reached the maximum value of 0.83; while at a 4000 rpm full load, the economic benefit improvement rate reached the minimum value of 0.03, which indicates that the economic and thermal benefits of the ORC system reached a relative balance when the system was under heavy load conditions, and it was difficult to improve the economic benefit by changing the thermal efficiency. Therefore, in the practical application of the ORC system, the economic benefit can be improved by reducing a part of the thermal

Thermo-Economic Optimization Results Analysis
In the process of TOPSIS multi-objective decision-making, when the optimal situation is selected under a single operating condition, it is actually a balanced state, and one of the objectives is often not the optimal situation. In Figure 10, the thermal efficiency of the optimal solution did not reach the maximum value, which is the ideal solution in the figure, and the SIC was not taken to the minimum value, which is the worst solution, that can be seen as a loss of part of the thermal efficiency to improve the economic efficiency in this case. Taking the data shown in Figure 10, considering the maximum thermal efficiency under this operating condition and the maximum SIC as the initial situation, the optimal solution obtained by the decision can be regarded as reducing the thermal efficiency by 7.93%, but increasing the thermal efficiency by 16.72% in terms of economic benefits. Figure 14 shows the ORC thermal efficiency loss rate and economic efficiency improvement rate of the optimization decision results under all of the operating conditions. Figure 14a represents the loss rate of thermal efficiency relative to the optimal solution under this operating condition, and Figure 14b indicates the economic benefit relative to the worst solution at this operating condition. The improvement rate in the economic benefits gradually decreased with the increase in the engine speed and torque. In the medium and low speed and load ranges, reducing the thermal efficiency can appropriately improve the economic benefits. When the engine was 1250 rpm and the load rate was 50%, the largest thermal efficiency loss rate was 0.12, and the corresponding economic improvement rate also reached the maximum value of 0.83; while at a 4000 rpm full load, the economic benefit improvement rate reached the minimum value of 0.03, which indicates that the economic and thermal benefits of the ORC system reached a relative balance when the system was under heavy load conditions, and it was difficult to improve the economic benefit by changing the thermal efficiency. Therefore, in the practical application of the ORC system, the economic benefit can be improved by reducing a part of the thermal efficiency at low and medium engine speeds and torques, while the thermal efficiency target is dominant at high speeds and high torques.  Figure 15 represents the variation in the BSFC optimization. The BSFC had a similar variation pattern before and after the addition of the ORC system, but the value after the addition of the ORC system decreased significantly compared to that before the optimization, with a maximum decrease of 8.3%. This is because the BSFC of the engine-ORC combined system is not only affected by parameters such as fuel consumption and brake power, but also by the net output power of the ORC system. When the engine speed was 1750 rpm and the torque was 354.90 N-m, the BSFC value of the combined system was the lowest, which was 238.9 g/(kW·h) −1 .

Bench Test and Model Validation
The mass flow rate of the working fluid, the net output power, and thermal efficiency of the ORC system under the full load condition of the engine speed from 850 rpm to 4500 rpm were measured through the test rig. Taking the inlet temperature, mass flow rate, and pressure of the exhaust and working fluid as input parameters of the PSO program, Figure 14. The optimization results of the two objectives: (a) thermal efficiency loss rate; (b) economic benefit improvement rate. Figure 15 represents the variation in the BSFC optimization. The BSFC had a similar variation pattern before and after the addition of the ORC system, but the value after the addition of the ORC system decreased significantly compared to that before the optimization, with a maximum decrease of 8.3%. This is because the BSFC of the engine-ORC combined system is not only affected by parameters such as fuel consumption and brake power, but also by the net output power of the ORC system. When the engine speed was 1750 rpm and the torque was 354.90 N-m, the BSFC value of the combined system was the lowest, which was 238.9 g/(kW·h) −1 .
Energies 2022, 15, x FOR PEER REVIEW 20 of 25 efficiency at low and medium engine speeds and torques, while the thermal efficiency target is dominant at high speeds and high torques.
(a) (b) Figure 14. The optimization results of the two objectives: (a) thermal efficiency loss rate; (b) economic benefit improvement rate. Figure 15 represents the variation in the BSFC optimization. The BSFC had a similar variation pattern before and after the addition of the ORC system, but the value after the addition of the ORC system decreased significantly compared to that before the optimization, with a maximum decrease of 8.3%. This is because the BSFC of the engine-ORC combined system is not only affected by parameters such as fuel consumption and brake power, but also by the net output power of the ORC system. When the engine speed was 1750 rpm and the torque was 354.90 N-m, the BSFC value of the combined system was the lowest, which was 238.9 g/(kW·h) −1 .

Bench Test and Model Validation
The mass flow rate of the working fluid, the net output power, and thermal efficiency of the ORC system under the full load condition of the engine speed from 850 rpm to 4500 rpm were measured through the test rig. Taking the inlet temperature, mass flow rate, and pressure of the exhaust and working fluid as input parameters of the PSO program,

Bench Test and Model Validation
The mass flow rate of the working fluid, the net output power, and thermal efficiency of the ORC system under the full load condition of the engine speed from 850 rpm to 4500 rpm were measured through the test rig. Taking the inlet temperature, mass flow rate, and pressure of the exhaust and working fluid as input parameters of the PSO program, the degree of superheating was set to 10, the degree of subcooling was set to 0, the mass flow rate of the coolant water was 0.15 kg/s, and the net output power and thermal efficiency of the ORC system were obtained. The comparison of the test rig, PSO numerical calculation, and BPNN fitting prediction result data are shown in Figure 16. It can be seen from the figure that the PSO numerical calculation program and the BPNN fitting results were basically consistent with the test rig results. Except at 3500 rpm, the relative errors of the PSO and BPNN results and the test rig were not more than 5%. the degree of superheating was set to 10, the degree of subcooling was set to 0, the mass flow rate of the coolant water was 0.15 kg/s, and the net output power and thermal efficiency of the ORC system were obtained. The comparison of the test rig, PSO numerical calculation, and BPNN fitting prediction result data are shown in Figure 16. It can be seen from the figure that the PSO numerical calculation program and the BPNN fitting results were basically consistent with the test rig results. Except at 3500 rpm, the relative errors of the PSO and BPNN results and the test rig were not more than 5%.
(a) (b) Figure 16. The comparison of the performance prediction models for the 10 test conditions with full load at different speeds: (a) net output power; (b) thermal efficiency.

Conclusions
1. The engine-ORC combined system can recover the exhaust energy of a tractor diesel engine. When the engine speed was 3500 rpm and full load torque was 323.8 N-m, the ORC thermal efficiency reached a maximum of 12.76% and the corresponding SIC value was 8539.66 $/kW. When the engine speed was 4500 rpm and the full load torque was 216.7 N m, the net output power of the ORC system was up to 8.31 kW. When the engine speed was 1750 rpm, the full load torque was 354.90 N·m, it had the minimum BSFC value of 238.9 g/(kW·h) −1 with the addition ORC system, and compared to the engine without the ORC system, the BSFC improvement was 8.3%. 2. The analysis of a single-case operating condition shows that with the increase in the evaporating pressure, the ORC thermal efficiency and SIC value will continue to increase, and when the condensation pressure is increased, the ORC thermal efficiency will decrease and the SIC value increase; increasing the ORC working fluid superheating had little effect on the thermal efficiency. The SIC value first increased and then decreased while the continuous increase of the degree of subcooling led to a decrease in the thermal efficiency and an increase in the SIC value. 3. The multi-objective optimization results showed that the optimal value of the evaporation pressure varied from 1.4 MPa to 3 MPa, the condensing pressure was less affected by the change in the working conditions, and was basically maintained at about 0.1 MPa. The optimized degree of superheating and subcooling were in the ranges of 1-10 K and 0-4.5 K, respectively. The mass flow of the working fluid was mainly affected by the operating conditions of the diesel engine. The optimized ratio of the mass flow rate of the working fluid to the exhaust was between 1.2 and 1.5. Under the conditions of low and medium speed and torque, the exhaust mass flow was small and the temperature was relatively low, and the working fluid mass flow

1.
The engine-ORC combined system can recover the exhaust energy of a tractor diesel engine. When the engine speed was 3500 rpm and full load torque was 323.8 N-m, the ORC thermal efficiency reached a maximum of 12.76% and the corresponding SIC value was 8539.66 $/kW. When the engine speed was 4500 rpm and the full load torque was 216.7 N m, the net output power of the ORC system was up to 8.31 kW. When the engine speed was 1750 rpm, the full load torque was 354.90 N·m, it had the minimum BSFC value of 238.9 g/(kW·h) −1 with the addition ORC system, and compared to the engine without the ORC system, the BSFC improvement was 8.3%.

2.
The analysis of a single-case operating condition shows that with the increase in the evaporating pressure, the ORC thermal efficiency and SIC value will continue to increase, and when the condensation pressure is increased, the ORC thermal efficiency will decrease and the SIC value increase; increasing the ORC working fluid superheating had little effect on the thermal efficiency. The SIC value first increased and then decreased while the continuous increase of the degree of subcooling led to a decrease in the thermal efficiency and an increase in the SIC value.

3.
The multi-objective optimization results showed that the optimal value of the evaporation pressure varied from 1.4 MPa to 3 MPa, the condensing pressure was less affected by the change in the working conditions, and was basically maintained at about 0.1 MPa. The optimized degree of superheating and subcooling were in the ranges of 1-10 K and 0-4.5 K, respectively. The mass flow of the working fluid was mainly affected by the operating conditions of the diesel engine. The optimized ratio of the mass flow rate of the working fluid to the exhaust was between 1.2 and 1.5. Under the conditions of low and medium speed and torque, the exhaust mass flow was small and the temperature was relatively low, and the working fluid mass flow varied between 0 and 0.05 kg/s, which increased with the increase in the engine speed and torque, and reached a maximum value of 0.23 kg/s at 4500 rpm under full load conditions. In most operating conditions, the ORC thermal efficiency remained around 12%. When the engine ran in the medium and high load region, the ORC system showed better thermal economy, and the maximum value was 12.76% under a full load at 3500 rpm.

4.
For the ORC, the improvement in the thermodynamic performance will inevitably lead to a sacrifice in economic performance, so economic benefits can be improved by reducing a part of the thermal efficiency. The analysis of the Pareto frontier under each operating condition showed that the loss rate of thermal efficiency was between 0 and 0.12, and the improvement rate in the economic benefits was 0 to 0.83. When the loss rate of thermal efficiency was the largest, the corresponding economic improvement rate could also reach the maximum value. At low and medium speeds, there was an obvious conflict between the thermal efficiency and economic benefits, and the economic benefits could be improved by appropriately reducing the thermal efficiency in the process of practical application, while at high speeds, the economic benefits and thermal efficiency reached a balance and showed a better thermal economic performance.

5.
Through the test rig, the performance parameters of the engine at 850-4500 rpm including 10 different speeds and loads, a total of 210 operating conditions were measured. The net output power and thermal efficiency of the ORC system with all of the operating conditions were obtained by adjusting the mass flow rate and pressure of the working fluid to verify the validity of the developed thermodynamic model. Comparing the test rig collected data with the PSO model optimization and the BPNN prediction results, it was found that, except at 3500 rpm, the maximum relative errors of the thermal efficiency and net power output were 5.3% and 5.585%, respectively, indicating that the PSO optimization model established in this paper has certain validity and usability. Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: Not applicable.

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