# Energy and Exergy Analyses of a PWR-Type Nuclear Power Plant Coupled with an ME-TVC-MED Desalination System

Zakaria Triki, Rabah Menasri, Mohamed Najib Bouaziz, Hichem Tahraoui, Mohammed Kebir, Abdeltif Amrane, Jie Zhang, Lotfi Mouni

► To cite this version:

Zakaria Triki, Rabah Menasri, Mohamed Najib Bouaziz, Hichem Tahraoui, Mohammed Kebir, et al. Energy and Exergy Analyses of a PWR-Type Nuclear Power Plant Coupled with an ME-TVC-MED Desalination System. Sustainability, 2023, 15 (10), pp.8358. 10.3390/su15108358. [https://hal.science/hal-04164311](https://hal.science/hal-04164311)

HAL Id: hal-04164311

https://hal.science/hal-04164311

Submitted on 18 Jul 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Energy and Exergy Analyses of a PWR-Type Nuclear Power Plant Coupled with an ME-TVC-MED Desalination System

Zakaria Triki <fnref n="1" /> ID, Rabah Menasri <fnref n="1" />, Mohamed Najib Bouaziz <fnref n="1" />, Hichem Tahraoui <fnref n="1" />,<fnref n="2" />,* ID, Mohammed Kebir <fnref n="3" /> ID, Abdeltif Amrane <fnref n="4" />,* ID, Jie Zhang <fnref n="5" /> ID and Lotfi Mouni <fnref n="6" /> ID

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Triki, Z.; Menasri, R.; Bouaziz, M.N.; Tahraoui, H.; Kebir, M.; Amrane, A.; Zhang, J.; Mouni, L. Energy and Exergy Analyses of a PWR-Type Nuclear Power Plant Coupled with an ME-TVC-MED Desalination System. *Sustainability* **2023**, *15*, 8358. https://doi.org/10.3390/su15108358

Academic Editors: Muhammad Rizwan Haider, Jinglong Han and Muhammad Bilal Asif

Received: 30 January 2023

Revised: 18 May 2023

Accepted: 18 May 2023

Published: 21 May 2023

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Electricity-water cogeneration power plants are an important tool for advancing sustainable water treatment technologies because they provide a cost-effective and environmentally friendly solution for meeting the energy and water needs of communities. By integrating power and water production, these technologies can reduce carbon emissions and help mitigate the impact of climate change. This work deals with the energy and exergy analysis of a cogeneration plant for electrical power generation and water desalination using real operational data. The power side is a pressurized water reactor (PWR) nuclear power plant (NPP), while the desalination side is a multi-effect distillation (MED) system with a thermo-vapor compressor (TVC) plant coupled with a conventional multi-effect plant (ME-TVC-MED). A mathematical model was implemented in MATLAB software and validated through a comparison with previously published research. The exergy analysis was carried out based on the second law of thermodynamics to evaluate the irreversibility of the plant and the subsystems. In this study, the components of the sub-systems were analyzed separately to identify and quantify the component that has a high loss of energy and exergy. According to the energy and exergy analyses, the highest source of irreversibility occurs in the reactor core with 50% of the total exergy destruction. However, turbines, steam generators, and condensers also contribute to energy loss. Further, the thermodynamic efficiency of the cogeneration plant was obtained as 35.38%, which is more effective than other systems. In the ME-TVC-MED desalination unit, the main sources of energy losses are located in the evaporators and the thermo-compressor (about 50% and 36%, respectively). Moreover, the exergetic efficiency of the ME-TVC-MED unit was found to be low at 6.43%, indicating a high degree of technical inefficiency in the desalination process. Therefore, many opportunities exist to improve the performance of the cogeneration system.

**Keywords:** PWR nuclear power; ME-TVC-MED desalination; cogeneration plant; energy analysis; exergy evaluation; performance

## 1. Introduction

Many cogeneration plants for the simultaneous production of electrical and thermal energy for desalination plants still mainly use conventional fuels, such as natural gas and coal. However, the increasing environmental pollution caused by conventional fuels and the ineffectiveness and exhaustion of available fuel reserves has prompted many researchers to explore other sustainable and environmentally friendly sources of energy,

<footnote n="1">¹ Laboratory of Biomaterials and Transport Phenomena, University of Medea, Medea 26000, Algeria</footnote>

<footnote n="2">² Laboratoire de Génie des Procédés Chimiques, Department of Process Engineering, University of Ferhat Abbas, Setif 19000, Algeria</footnote>

<footnote n="3">³ Research Unit on Analysis and Technological Development in Environment (URADTE-CRAPC), BP 384, Bou-Ismail Tipaza 42004, Algeria</footnote>

<footnote n="4">⁴ Ecole Nationale Supérieure de Chimie de Rennes, Centre National de la Recherche Scientifique (CNRS), ISCR-UMR 6226, Université de Rennes, F-35000 Rennes, France</footnote>

<footnote n="5">⁵ School of Engineering, Merz Court, Newcastle University, Newcastle upon Tyne NE1 7RU, UK</footnote>

<footnote n="6">⁶ Laboratory of Management and Valorization of Natural Resources and Quality Assurance, SNVST Faculty, Akli Mohand Oulhadj University, Bouira 10000, Algeria</footnote>

<footnote>* Correspondence: tahraoui.hichem@univ-medea.dz (H.T.); abdeltif.amrane@univ-rennes1.fr (A.A.)</footnote>
such as renewable and nuclear energy. They are the most suitable and reliable candidates for industrial processes, such as desalination. Thus, nuclear energy, which avoids many of the problems associated with fossil fuels, is an efficient and advantageous option for the main desalination plants [1,2].

Among other clean and renewable resources, nuclear desalination is a sustainable source of energy and potable water that could assist in addressing the present issues with access to electricity and potable water. As such, it should be considered in public policy to protect the environment. However, the multiple subsystems of nuclear power plants (NPP) result in some irreversibilities and corresponding energy losses. Therefore, reducing the irreversibilities becomes a crucial duty in order to increase the amount of energy that is available in the plant. Undoubtedly, defining the irreversibility quantities in various components requires significant effort [2,3].

Exergy analysis is an essential diagnostic tool for energy conversion systems since it pinpoints where and how much energy is lost throughout a process. Additionally, the origins of the energy losses will be revealed. The minimum effort of separation is a crucial consideration during the design stage of thermal desalination processes in order to reduce entropy generation, which lowers freshwater costs, pollution, and greenhouse gas emissions. In addition, by reducing the entropy, the system can operate more efficiently, leading to a reduction in the amount of brine discharge into the sea [4,5].

Despite the significance of the electricity-water cogeneration principle, very few studies on the energetic analysis of such systems have been published. Additionally, despite having a high-performance ratio, several desalination technologies have not yet been researched. Kambiz Ansari et al. [6,7] performed an exergetic and exergo-economic analysis supplemented with multi-objective optimization of the coupling of a desalination plant by multi-effect thermal vapor compression distillation (MED-TVC) with a nuclear power plant with a pressurized nuclear reactor (PWR). A genetic algorithm (GA) was used to minimize the cost of system operation (cost of producing electricity and freshwater) and maximize the energy efficiency of the system. The results showed that the costs of thermodynamic inefficiencies (the cost of exergetic destruction) for the optimized system components have been reduced compared to the corresponding cost of exergetic destruction in the base scenario plant. In addition, it has been confirmed that their optimization improved the thermodynamic and economic characteristics of the system.

Khoshgoftar Manesh et al. [8,9] carried out a multi-objective optimization of a PWR nuclear power plant coupled to a multi-flash distillation (MSF) desalination system through GA and mixed-integer non-linear mathematical programming methods. In addition, a computer program has been developed in the MATLAB environment for a thermoeconomic analysis. The results demonstrated that the evolutionary algorithm (NSGA-II) could be systemically and elegantly applied to the PWR-MSF dual-purpose plant.

In addition, Khalid et al. [10] described a comparative evaluation of two nuclear desalination systems. The reverse osmosis (RO) system used for desalination is coupled to a CANDU6 nuclear reactor and a fast reactor cooled with SFR sodium. Exergetic analysis was used to evaluate the cogeneration systems' performance. The results showed that the exergetic efficiencies of the CANDU-6 and SFR systems are respectively 32.8% and 36.8%. In addition, the exergetic efficiencies of the RO process for the CANDU 6 and SFR systems are 49.2% and 36.3%, respectively. In addition, Khalid et al. [11] presented a parametric study of the coupling of a modular helium gas turbine reactor (GT-MHR) with the RO process performed to reveal the effect of certain parameters on energy efficiency. The analysis showed that the exergy efficiency of cogeneration processes increases by 10.3% if the waste heat from the reactor is used.

Priego et al. [12] presented a study on nuclear desalination using the SMART reactor coupled with two thermal desalination processes, multi-effect (MED) and multi-flash (MSF) distillation. The study was carried out using exergetic and thermo-economic analyses. In either case, the analyses were carried out for gain output ratios (GORs) of 5, 10, and 15. The results indicated that economic competitiveness and the feasibility of nuclear desalination
using the main stream line of the reactor. In addition, the greatest quantity of water for every extraction position in the MED and MSF processes as well is produced by the gain ratio of 15.

The number of reported studies concerned with nuclear desalination in the literature is increasing. A comprehensive overview of the most recent studies on the various aspects of nuclear desalination is provided by Al-Othman et al. [13]. The review highlighted the main advantages of this coupled technology and the principal challenges facing its development and extensive industrial use. However, despite the importance of the electricity-water cogeneration principle, until recently, few studies have been published on using exergetic analysis for such systems. In addition, some desalination technologies that exhibit high-performance ratios have not been explored yet. Hence, it is valuable to delve into this field to gain a more comprehensive understanding of current facilities and future undertakings.

In this paper, energy and exergy analysis was conducted on a multi-effect thermal vapor compression plant (ME-TVC-MED) combined with a PWR nuclear power plant based on real operational data for the simultaneous production of electricity and freshwater. A mathematical model of the PWR-ME-TVC-MED cogeneration plant was developed using the most recent and accurate correlations and developed values of the properties of seawater, including chemical exergy, which provide more novel approaches for the design and operation of the plant, leading to improved efficiency, performance, and sustainability. The model takes into account the behavior of the nuclear reactor, the thermo-vapor compression system, and the multi-effect distillation system, as well as their interactions with each other and with the environment. This approach allows for maximum efficiency and energy savings while reducing greenhouse gas emissions and water consumption compared to traditional power plants.

## 2. Cogeneration Plant Description

In this study, a cogeneration system based on a combined NPP and desalination plant is considered. The overall system is shown in Figure 1 with a schematic process flow shown in Figure 2. The energy side is the Angra (II) nuclear power plant in Brazil, while the water desalination side is based at the Al-Jubail desalination plant (ME-TVC-MED), which is located in the Kingdom of Saudi Arabia (KSA). The proposed NPP is equipped with a pressurized light-water reactor (PWR). A 3855 MW thermal reactor is used to generate the nominal power of 1300 MWe, operating under local climate conditions. The technical data of the PWR power plant are summarized in Appendix A. The Al-Jubail desalination plant's ME-TVC-MED consists of two Multi-Effect Distillation (MED) units equipped with thermo-vapor compressors (TVCs) integrated with a conventional MED unit to form a single entity.

Referring to Figure 1, the NPP plant consists of three main cycles, including the primary, secondary and tertiary cycles. In the primary cycle, the four similar cooling cycles (a reactor and four steam generators, four pumps, and four pressurizers) use water as the heat transfer fluid. This circuit is always kept under pressure to prevent vaporization of the heat transfer fluid at the reactor's operating temperature. The pressurizer is a component of the loop that maintains the pressure of the coolant loop. In the primary cycle, thermal energy is dissipated by the fission reaction of nuclear fuel at the heart of the reactor.

The cycle also includes several additional components, such as low-pressure (LP) heaters, high-pressure (HP) heaters, a feedwater tank (FWT), a water or moisture separator (MS), and a reheater (RH). The LP and HP heaters preheat water by recovering heat to increase its temperature before it enters the steam generator. The FWT stores and supplies water to the steam generator, while the MS removes any moisture from the steam to ensure its purity. The RH superheats steam before it enters the low-pressure turbine, resulting in improved overall system performance. Figure 1 provides an overview of the position of these devices in the nuclear power plant, as well as the thermodynamic cycle for the system.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Schematic flow diagram of the hybrid PWR-ME-TVC-MED system (adapted from [14]).</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Schematic flow diagram of the ME-TVC-MED desalination system (adapted from [15]).</figcaption></figure>

The heat generated by the fission reaction is transferred to the coolant (water). The coolant is recirculated in the steam generators and evaporates the working fluid (steam) from the second cycle. The secondary circuit is based on a Rankine cycle similar to that of most steam power plants. The steam that comes out of the steam generator is divided into two streams; one part of the outlet steam is directed to the HPT turbine, and the rest is used to heat the flows extracted from the humidity separator. The superheated steam from the pre-heater is directed to the BPT turbine. Finally, in the third circuit, the seawater cools the vapor in the condenser after it leaves the low-pressure turbine. This cooling process is
necessary to guarantee all the thermodynamic processes of the installation, including the cooling of the reactor, the production of steam, and the production of electricity.

The outlet flow to the first extraction line from the BPT turbine designated by flow number 44 is divided into two flows, numbered as 46 and 14. Flow 46 is directed to an intermediate insulation loop heat exchanger (designated as EX in Figure 1). The latent heat from flow 46 is transferred to the circulating water stream, increasing its temperature by about 5 °C. Some of the heated water then enters the expansion chamber, forming the vapor designated as Min. This vapor (Min) is utilized as the driving vapor in the first stage of the MED-TVC process. The heat exchanger (EX) is implemented to prevent direct contact between the steam from the PWR power plant and the water from the desalination plant, avoiding the risk of radioactive contamination in the freshwater. Flow 46 used to generate the driving steam for the ME-MED-TVC system is condensed in the EX and then returned to the second cycle of the PWR as flow 45.

Figure 2 illustrates a process flow diagram of the ME-TVC-MED, with the relevant operational data presented in Appendix B. The condensates of the first effects in the two MED-TVC units are divided into two flows (S and D<sub>r</sub>); the first returns to its source, while the second pass is to the product of desalinated water. The vapor generated (D) in the first effect of each MED-TVC unit is directed to the second effect as side effect and acts as source of heat in the next effect at lower pressure and temperature. The hot brine leaving the first effect in each MED-TVC unit passes into the second effect through the brine (B<sub>1</sub>). This iterative process is repeated through all the effects up to the final effect (N) of each MED-TVC unit.

The second flow (D<sub>v</sub>) is combined with the flow from the other MED-TVC unit and directed to the first effect of the MED unit to start the desalination process at a low temperature. This process is the same in the MED-TVC unit and repeated until the last effect (n) in the MED unit. The steam produced by the final effect is ultimately sent to the condenser, where it transfers its latent heat to the cooling water M<sub>c</sub>, raising its temperature from T<sub>sw</sub> to T<sub>f</sub>. The remainder of this cooling water is returned to the sea, while some of it is used as feed for the various effects (R<sub>rej</sub>).

### 3. Modeling and System Analysis

#### 3.1. Mathematical Laws and Assumptions

The mathematical model representing the thermodynamic and exergy analysis of PWR-MED-TVC-MED cogeneration includes three main equations: mass balance, energy balance, and exergy balance according to the first and second laws of thermodynamics [7,16,17]. These equations are as follows:

* The law of mass balance for each component of the power plant:

$$ \sum_i (\dot{m}_{in} - \dot{m}_{out}) = 0 \qquad (1) $$

* The first law of thermodynamics applies to each component of the power plant:

$$ \sum_i \dot{m}_i h_i - \sum_e \dot{m}_e h_e - \dot{W} + \dot{Q} = 0 \qquad (2) $$

* Exergy flow equation for each part of the power plant:

$$ \sum_i \dot{Q} \left( 1 - \frac{T_0}{T} \right) - \dot{W} + \sum_i \dot{m}_i e_i - \sum_e \dot{m}_e e_e = \dot{S}_g \qquad (3) $$

To simplify the analysis, some assumptions are made as follows:
* The hybrid system model is operating at a steady state.

* The potential and kinetic exergies are negligible.

* The pressure losses in all pipelines and heat exchangers are negligible.

* Heat losses in turbines and pumps are considered.

* The temperature difference between the effects is constant.

* The feed-flow rate in the effects is the same.

* Calculations have been made for the boiling point elevation (BPE), specific heat capacity, and other variables.

* The distillate is salt-free.

* Thermodynamic losses have been studied.

## 3.2. Mass and Energy Balance Model

For the evaluation of mass and energy conservation, they should be determined with some unknown parameters, such as the heat transfer coefficient and physical properties of water. These unknown variables are defined in the following section. Mass and salinity balance equations for all the effects, the condenser, and the distillate tank are given by Equations (4)–(18), as shown in Table 1.

<table><caption>**Table 1.** Mass and energy balance equations.</caption><thead><tr><th>Equation</th><th>Description</th><th>Eq.</th></tr></thead><tbody><tr><td>B<sub>1</sub> = F<sub>1</sub> - D<sub>1</sub></td><td>Mass balance of effect 1</td><td>(4)</td></tr><tr><td>B<sub>i</sub> = F<sub>i</sub> + B<sub>i-1</sub> + D<sub>i</sub></td><td>Mass balance of effects 2 to N</td><td>(5)</td></tr><tr><td>B<sub>i</sub> = F<sub>i</sub> + 2B<sub>i-1</sub> + D<sub>i</sub></td><td>Mass balance of effect N + 1</td><td>(6)</td></tr><tr><td>B<sub>i</sub> = F<sub>i</sub> + B<sub>i-1</sub> + D<sub>i</sub></td><td>Mass balance of effects N + 2 to n</td><td>(7)</td></tr><tr><td>X<sub>1</sub> = F<sub>1</sub>/B<sub>1</sub> × X<sub>f</sub></td><td>Salinity balance of effect 1</td><td>(8)</td></tr><tr><td>X<sub>i</sub> = F<sub>i</sub>/B<sub>i</sub> × X<sub>f</sub> + B<sub>i-1</sub>/B<sub>i</sub> × X<sub>i-1</sub>,</td><td>Salinity balance of effects 2 to N</td><td>(9)</td></tr><tr><td>X<sub>i</sub> = F<sub>i</sub>/B<sub>i</sub> × X<sub>f</sub> + 2B<sub>i-1</sub>/B<sub>i</sub> × X<sub>i-1</sub></td><td>Salinity balance of effect N + 1</td><td>(10)</td></tr><tr><td>X<sub>i</sub> = F<sub>i</sub>/B<sub>i</sub> × X<sub>f</sub> + B<sub>i-1</sub>/B<sub>i</sub> × X<sub>i-1</sub></td><td>Salinity balance of effects N + 2 to n</td><td>(11)</td></tr><tr><td>D<sub>1</sub> = 1/λ<sub>1</sub> [(S + D<sub>r</sub>)λ<sub>s</sub> - F<sub>1</sub> × C<sub>p</sub>(T<sub>1</sub> - T<sub>f</sub>)]</td><td>The vapor generated in the first effect</td><td>(12)</td></tr><tr><td>D<sub>i</sub> = 1/λ<sub>i</sub> [(D<sub>i-1</sub> + d'<sub>i-1</sub>)λ<sub>i-1</sub> - F<sub>i</sub>C<sub>p</sub>(T<sub>i</sub> - T<sub>f</sub>) - B<sub>i-1</sub>C<sub>p</sub>(T<sub>i-1</sub> - T<sub>i</sub>)]</td><td>The amount of steam released by the second up to the Nth effects</td><td>(13)</td></tr><tr><td>D<sub>N+1</sub> = 1/λ<sub>N+1</sub> [(2D<sub>v</sub>)λ<sub>N</sub> - F<sub>N+1</sub>C<sub>p</sub>(T<sub>N+1</sub> - T<sub>f</sub>) - 2B<sub>N</sub>C<sub>p</sub>(T<sub>N</sub> - T<sub>N+1</sub>)]</td><td>The vapor formed in the effect N + 1</td><td>(14)</td></tr><tr><td>D<sub>i</sub> = 1/λ<sub>i</sub> [(D<sub>i-1</sub> + d'<sub>i-1</sub>)λ<sub>i-1</sub> - F<sub>i</sub>C<sub>p</sub>(T<sub>i</sub> - T<sub>f</sub>) - B<sub>i-1</sub>C<sub>p</sub>(T<sub>i-1</sub> - T<sub>i</sub>)</td><td>The vapor formed in the effects (N + 2) to n</td><td>(15)</td></tr><tr><td>d'<sub>i</sub> = D<sub>i-1</sub>C<sub>p</sub>(T<sub>v,i-1</sub> - T<sub>i</sub>)</td><td>The amount of vapor removed from the flash boxes</td><td>(16)</td></tr><tr><td>M<sub>cw</sub> = D<sub>n</sub>λ<sub>n</sub>/C<sub>p</sub>(T<sub>f</sub> - T<sub>fw</sub>)</td><td>Cooling seawater flow rate</td><td>(17)</td></tr><tr><td>D<sub>ev</sub> = S/Ra</td><td>The amount of entrained steam</td><td>(18)</td></tr></tbody></table>

The temperature profile equations to determine the temperature difference across the effects, saturated vapor temperature, vapor condensation temperature, brine temperature, boiling point elevation, and non-equilibrium allowance are given by Equations (19)–(26), as shown in Table 2. The heat transfer area for each effect, the pre-heaters, the condenser, the total heat transfer area, the overall heat transfer coefficients, and the logarithmic mean temperature differences for the effects, pre-heaters, condenser, and the latent heat of the steam and distillate vapor can be obtained by using Equations (27)–(35), as shown in Table 3. The pressure equations to determine the compressed and entrained vapors can be obtained by using Equations (36) and (37), as shown in Table 4. In addition, the expansion ratio, compression, ratio, entrainment ratio, and heat capacity of water are presented in this table.
<table><caption>**Table 2.** Temperature profiles.</caption><thead><tr><th>Equation</th><th>Description</th><th>Eq.</th></tr></thead><tbody><tr><td>ΔT = <sup>T1-Tn</sup>/<sub>n-1</sub></td><td>Temperature difference across the effects</td><td>(19)</td></tr><tr><td>T<sub>s</sub> = T1 + ΔT</td><td>Temperature of compressed steam</td><td>(20)</td></tr><tr><td>T<sub>vn</sub> = T<sub>n</sub> - BPE</td><td>Vapor temperature in the last effect</td><td>(21)</td></tr><tr><td>T<sub>ci</sub> = T<sub>i</sub> - BPE - ΔT<sub>p</sub> + ΔT<sub>t</sub> - ΔT<sub>c</sub></td><td>Vapor condensation temperature of effects</td><td>(22)</td></tr><tr><td>T<sub>i+1</sub> = T<sub>i</sub> - ΔT<sub>t</sub></td><td>Saturated vapor temperature of effects</td><td>(23)</td></tr><tr><td>T'<sub>i</sub> = T<sub>v,i-1</sub> - NEA<sub>i</sub></td><td>Temperature of the vapor formed by flashing</td><td>(24)</td></tr><tr><td>BPE = X<sub>b</sub>[B + CX<sub>b</sub>] × 10<sup>-3</sup></td><td>Boiling point elevation (BPE)</td><td>(25)</td></tr><tr><td>B = [6.71 + (6.34 × 10<sup>-2</sup> × T<sub>n</sub>) + (9.74 × 10<sup>-5</sup>T<sub>n</sub><sup>2</sup>)] 10<sup>-3</sup></td><td></td><td></td></tr><tr><td>C = [22.238 + ((9.59 × 10<sup>-3</sup> × T<sub>n</sub>) + (9.42 × 10<sup>-5</sup>T<sub>n</sub><sup>2</sup>)] 10<sup>-8</sup></td><td></td><td></td></tr><tr><td>NEA<sub>i</sub> = 33 <sup>(T<sub>i-1</sub>-T<sub>i</sub>)<sup>0.55</sup></sup>/<sub>T<sub>v,i</sub></sub></td><td>Non-equilibrium tolerance (NEA<sub>i</sub>)</td><td>(26)</td></tr></tbody></table>

<table><caption>**Table 3.** Heat transfer area, heat transfer coefficient, logarithmic mean temperature difference, and latent heat steam equations.</caption><thead><tr><th>Equation</th><th>Description</th><th>Eq.</th></tr></thead><tbody><tr><td>A<sub>1</sub> = <sup>(S+D<sub>r</sub>)λ<sub>s</sub></sup>/<sub>U<sub>i</sub>(T<sub>s</sub>-T<sub>i</sub>)</sub></td><td>The heat transfer area of the 1st effect</td><td>(27)</td></tr><tr><td>A<sub>i</sub> = <sup>D<sub>i</sub>λ<sub>i</sub></sup>/<sub>U<sub>i</sub>(T<sub>v,i</sub>-T<sub>i</sub>)</sub></td><td>The heat transfer area of the 1st effect</td><td>(28)</td></tr><tr><td>A<sub>e</sub> = 2 ∑<sub>i=1</sub><sup>N</sup> A<sub>i</sub> + ∑<sub>i=N+1</sub><sup>n</sup> A<sub>i</sub></td><td>The total heat transfer area</td><td>(29)</td></tr><tr><td>A<sub>c</sub> = <sup>D<sub>n</sub>λ<sub>n</sub></sup>/<sub>U<sub>c</sub>(LMTD)<sub>c</sub></sub></td><td>The heat transfer area of the condenser</td><td>(30)</td></tr><tr><td>U<sub>c</sub> = 1.7194 + 3.2063 × 10<sup>-2</sup>T<sub>v,n</sub> - 1.5971 × 10<sup>-5</sup>T<sub>v,n</sub><sup>2</sup> + 1.9918 × 10<sup>-7</sup>T<sub>v,n</sub><sup>3</sup></td><td>The overall heat transfer coefficient of the condenser</td><td>(31)</td></tr><tr><td>U<sub>i</sub> = <sup>1939.4+1.40562T<sub>i</sub>-0.020752T<sub>i</sub><sup>2</sup>+0.0023186T<sub>i</sub><sup>3</sup></sup>/<sub>(LMTD)<sub>c</sub> = <sup>(T<sub>f</sub>-T<sub>sw</sub>)</sup>/<sub>ln[<sup>T<sub>v,n</sub>-T<sub>sw</sub></sup>/<sub>T<sub>c,n</sub>-T<sub>f</sub></sub>]</sub></td><td>The overall heat transfer coefficient</td><td>(32)</td></tr><tr><td>λ<sub>s</sub> = 2501.897149 - 2.407064037 × T<sub>s</sub> + 1.192217 × 10<sup>-3</sup> × T<sub>s</sub><sup>2</sup> - 1.5863 × 10<sup>-5</sup> × T<sub>s</sub><sup>3</sup></td><td>Latent heat of steam</td><td>(34)</td></tr><tr><td>λ<sub>i</sub> = 2501.897149 - 2.407064037 × T<sub>i</sub> + 1.192217 × 10<sup>-3</sup> × T<sub>i</sub><sup>2</sup> - 1.5863 × 10<sup>-5</sup> × T<sub>i</sub><sup>3</sup></td><td>Latent heat of distillate vapor</td><td>(35)</td></tr></tbody></table>

<table><caption>**Table 4.** Pressure and heat capacity of water equations [18].</caption><thead><tr><th>Equation</th><th>Description</th><th>Eq.</th></tr></thead><tbody><tr><td>P<sub>s</sub> = 1000 × exp(<sup>-3892.7</sup>/<sub>T<sub>s</sub>+273.15-42.6776</sub> + 9.5)</td><td>Pressure of the compressed vapor</td><td>(36)</td></tr><tr><td>P<sub>ev</sub> = 1000 × exp(<sup>-3892.7</sup>/<sub>T<sub>vn</sub>+273.15-42.6776</sub> + 9.5)</td><td>Pressure of the entrained vapor</td><td>(37)</td></tr><tr><td>ER = <sup>P<sub>n</sub></sup>/<sub>P<sub>ev</sub></sub></td><td>Expansion ratio (ER)</td><td>(38)</td></tr><tr><td>CR = <sup>P<sub>s</sub></sup>/<sub>P<sub>ev</sub></sub></td><td>Compression ratio (CR)</td><td>(39)</td></tr><tr><td>Ra = 0.235 <sup>P<sub>s</sub><sup>1.19</sup></sup>/<sub>P<sub>ev</sub><sup>0.04</sup></sub> ER<sup>0.015</sup></td><td>Entrainment ratio (Ra)</td><td>(40)</td></tr><tr><td>C<sub>p</sub> = [a + b × T<sub>1</sub> + c × T<sub>1</sub><sup>2</sup> + d × T<sub>1</sub><sup>3</sup>] × 10<sup>-3</sup></td><td></td><td></td></tr><tr><td>a = 4206.8 - (6.6197 × S) + (1.2288 × 10<sup>-2</sup> × S<sup>2</sup>)</td><td></td><td></td></tr><tr><td>b = -1.1262 + (5.4178 × 10<sup>-2</sup> × S) - (2.2719 × 10<sup>-4</sup> × S<sup>2</sup>)</td><td>Heat capacity of water</td><td>(41)</td></tr><tr><td>c = (1.2026 × 10<sup>-2</sup>) - (5.3566 × 10<sup>-4</sup> × S) + (1.8906 × 10<sup>-6</sup> × S<sup>2</sup>)</td><td></td><td></td></tr><tr><td>d = (6.8777 × 10<sup>-7</sup>) + (1.517 × 10<sup>-6</sup> × S) - (4.4268 × 10<sup>-9</sup> × S<sup>2</sup>)</td><td></td><td></td></tr></tbody></table>

### 3.3. Exergy Analysis

Exergy is the maximum theoretical useful work that can be obtained from a given form of energy in a system using the environmental parameters as the reference state. In order to arrive at Equation (3), a control volume undergoing a steady-flow process may
have multiple inlets and outlets and may exchange heat with the environment, which acts as a reservoir at a constant pressure $P_0$ and a constant temperature $T_0$ [19]. The second law of thermodynamics for this process can be expressed as follows:

$$ \dot{S}_g = \sum_e \dot{m}_e e_e - \sum_i \dot{m}_i e_i + \frac{\dot{Q}_J}{T_0} \qquad (42) $$

where $\dot{S}_g$ represents the entropy generation rate for the process due to the irreversibilities. Two terms, namely, $\sum_e \dot{m}_e e_e$ and $\sum_i \dot{m}_i e_i$, are considered as the entropy transfer. $\dot{Q}_J = -\dot{Q}$

denotes the heat transfer rate for the instantaneous temperature $T_0$. The last term, $\frac{\dot{Q}_J}{T_0}$, stands for the entropy transfer rate.

By neglecting the heat transfer as well as the kinetic and potential energies of the stream, the following expression is obtained:

$$ \dot{W}_{cv} = \dot{W}_u = \sum_i \dot{m}_i e_i - \sum_e \dot{m}_e e_e - T_0 \dot{S}_g \qquad (43) $$

This equation calculates the actual work accomplished during the process. This also defines the useful work, as steady-flow components have fixed boundaries and do not provide work performed by or against the environment. The reversible work can be determined by making the entropy generation term $\dot{S}_g = 0$. Hence, Equation (43) can be expressed as follows:

$$ \dot{W}_{rev} = \dot{W}_{u,max} = \dot{m}(e_i - e_e) \qquad (44) $$

The general exergy balance of a system can be written as follows:

$$ \dot{E}x = \dot{E}_{ph} + \dot{E}_{ch} + \dot{E}_{ke} + \dot{E}_{pe} \qquad (45) $$

where $\dot{E}_{ph}$, $\dot{E}_{ch}$, $\dot{E}_{ke}$, and $\dot{E}_{pe}$ are the physical, chemical, kinetic, and potential exergies, respectively. As for the selected reference of this study, the kinetic and potential exergy can be neglected from the exergy balance.

Given that the chemical composition of the flows in the PWR plant remains constant, the chemical exergy is disregarded in the energy balance equations applied to the PWR system. However, it is important to keep in mind that in the ME-TVC-MED desalination system, the salinity of the seawater changes during the process, making it necessary to consider the significance of the chemical exergy.

The physical and chemical exergies can be respectively provided as follows:

$$ \dot{E}_{ph} = \dot{m}e_{ph} = \dot{m}[(h - h_0) - T_0(S - S_0)] \qquad (46) $$

$$ \dot{E}_{ch,w} = \dot{m}e_{ch} = \dot{m} \sum_k w_k (\mu_k^s - \mu_k^o) \qquad (47) $$

where the subscript s refers to the initial state and o to the corresponding environmental or reference state, respectively. The terms $\mu$ and $w$ refer to the chemical potential and mass fraction respectively. The terms $e_{ph}$ and $e_{ch}$ refer to the specific physical exergy and specific chemical exergy respectively.

The reversible (maximum) work can be expressed as follows:

$$ \dot{W}_{u,max} = \sum_i \dot{E}x_i - \sum_e \dot{E}x_e \qquad (48) $$
The paper uses the term exergy flow rate instead of physical exergy and chemical exergy for simplicity. The rate of irreversibility (İ) is equal to exergy loss:

$$ \dot{I} = \dot{W}_{u,max} - \dot{W}_{u} = T_{0}\dot{S}_{g} \quad (49) $$

The energy and exergy balance for each component of the NPP is based on reference [1], and the mathematical model is presented in the supplementary material.

The ME-TVC-MED desalination unit has three streams: pure water, saline water, and steam. The thermodynamic properties of pure water and steam can be obtained either from thermodynamic tables or through the use of a specific equation [16]. In this study, the properties of seawater were calculated using empirical relationships that depend on temperature and concentration, under constant atmospheric pressure, as suggested by Sharqawy et al. [20]. These correlations can be adjusted for different pressures with equations developed by the same authors. These correlations have been corroborated by Sharqawy et al. [21].

### 3.4. System Performance

The system performance of the PWR-ME-TVC-MED model can be evaluated in terms of the following [15]:

* Total distillate ($D_t$)

The total distillate output of all effects is equal to the following:

$$ D_t = 2 \sum_{i=1}^{N} D_i + \sum_{i=N+1}^{n} D_i \quad (50) $$

* Gain Output Ratio (GOR)

Gain output ratio is one of the parameters used to evaluate the performance of thermal desalination processes. It is defined as the ratio of the total distilled water produced ($D_t$) and the motive steam supplied ($S$).

$$ \text{GOR} = \frac{D_t}{2S} \quad (51) $$

* Specific heat consumption (Q)

This stands as one of the most paramount characteristics of the systems of thermal desalination. It is defined as the thermal energy to be consumed by the system to produce 1 kg of distilled water, where $\lambda_m$ is the motive steam latent heat in kJ/kg.

$$ Q = \frac{2S\lambda_m}{D_t} \quad (52) $$

* Specific heat transfer area ($A_d$)

The specific total heat transfer area is equal to the sum of the heat transfer areas of the effects ($A_e$) and the condenser ($A_c$) per total distillate product (m²/kg/s).

$$ A_d = \frac{A_e + A_c}{D_t} \quad (53) $$

* Exergetic efficiency ME-TVC-MED

The exergetic efficiency represents the ratio of the minimum work of separation to fuel exergy supplied to the system and can be expressed as follows:

$$ \varepsilon_{MED} = \frac{\dot{W}_{min}}{\dot{E}_f} \quad (54) $$
The minimum work of separation ($\dot{W}_{min}$) refers to the product exergy in the desalination process, while the heat input represents the fuel exergy ($\dot{E}_f$).

* Exergy efficiency of the cogeneration plant [7]:

$$ \varepsilon_{PWR+MED} = \frac{(\dot{W}_{u,max})_{fis} - \sum_{PWR+MED} \dot{I}}{(\dot{W}_{u,max})_{fis}} \quad (55) $$

## 4. Results and Discussion

The ME-TVC-MED desalination unit performance results are shown in Table 5. The results have been validated based on a comparison with three units available in the literature [15,22], as illustrated in Table 6. In this study, the thermal balance diagram of the combined PWR-ME-TVC-MED system is shown in Figures 3 and 4. The diagram includes a total of 70 states, with the first three states located on the primary circuit of the PWR unit. The states from 4 to 60 are located along the flow stream on the secondary side piping unit of the PWR unit, while states 47 to 70 are found along the flow stream of the desalination unit at the inlet and outlet of its components. Streams 41 and 42 are located along the flow stream of the closed-loop cooling system piping unit at the inlet and outlet of the condenser. The properties of seawater are calculated using empirical correlations based on temperature and concentration at a constant atmospheric pressure.

<table><caption>Table 5. Performance data of the desalination unit.</caption><thead><tr><th>Description</th><th>Value</th><th>Unit</th></tr></thead><tbody><tr><td>Temperature difference, ΔT</td><td>3.57</td><td>°C</td></tr><tr><td>Entrainment ratio, Ra</td><td>0.806</td><td>-</td></tr><tr><td>Expansion ratio, ER</td><td>6.667</td><td>-</td></tr><tr><td>Compression ratio, CR</td><td>1.953</td><td>-</td></tr><tr><td>Distillate production, D<sub>t</sub></td><td>334.66</td><td>kg/s</td></tr><tr><td>GOR</td><td>9.56</td><td>-</td></tr><tr><td>Cooling water flow rate, M<sub>cw</sub></td><td>1266.8</td><td>kg/s</td></tr><tr><td>Specific heat consumption, q</td><td>197.83</td><td>kJ/kg</td></tr><tr><td>Specific heat transfer area, A<sub>d</sub></td><td>302.23</td><td>m<sup>2</sup>/kg/s</td></tr></tbody></table>

<table><caption>Table 6. Mathematical model validation.</caption><thead><tr><th>Desalination Plants</th><th>Model</th><th>Ref. [15]</th><th>Model</th><th>Ref. [22]</th><th>Model</th><th>Ref. [23]</th></tr></thead><tbody><tr><td><strong>Operating and design conditions</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Number of effects n</td><td>8</td><td>8</td><td>6</td><td>6</td><td>4</td><td>4</td></tr><tr><td>Motive pressure P<sub>m</sub>, kPa</td><td>2700</td><td>2700</td><td>2800</td><td>2800</td><td>21,000</td><td>21,000</td></tr><tr><td>Top brine temperature T(1), °C</td><td>63</td><td>63</td><td>63</td><td>63</td><td>63</td><td>63</td></tr><tr><td>Minimum brine temperature T(n), °C</td><td>42</td><td>42</td><td>44</td><td>44</td><td>48</td><td>48</td></tr><tr><td>Temperature drop per effect, °C</td><td>3</td><td>3</td><td>3.8</td><td>3.8</td><td>5</td><td>5</td></tr><tr><td>Feed seawater temperature T<sub>f</sub>, °C</td><td>37</td><td>37</td><td>40</td><td>40</td><td>43</td><td>43</td></tr><tr><td>Motive steam flow rate D<sub>m</sub>, kg/s</td><td>17.5 × 2</td><td>17.5 × 2</td><td>11 × 2</td><td>11 × 2</td><td>8.5 × 2</td><td>8.5 × 2</td></tr><tr><td><strong>TVC design</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Entrainment ratio Ra</td><td>0.67</td><td>0.98</td><td>0.65</td><td>NA</td><td>0.63</td><td>NA</td></tr><tr><td>Expansion ratio ER</td><td>17.6</td><td>18.7</td><td>17.1</td><td>NA</td><td>11.5</td><td>NA</td></tr><tr><td>Compression ratio CR</td><td>0.18</td><td>1.75</td><td>0.17</td><td>NA</td><td>1.59</td><td>NA</td></tr><tr><td><strong>System performance</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Distillate production D<sub>t</sub>, kg/s</td><td>378.94</td><td>349</td><td>188.03</td><td>184.38</td><td>114.8</td><td>127</td></tr><tr><td>Gain output ratio GOR</td><td>10.82</td><td>10</td><td>8.54</td><td>8.6</td><td>6.75</td><td>7.5</td></tr><tr><td>Specific heat consumption Q, kJ/kg</td><td>174.71</td><td>223</td><td>288.3</td><td>287.5</td><td>280.12</td><td>NA</td></tr><tr><td>Specific heat transfer area A<sub>d</sub>, m<sup>2</sup>/kg/s</td><td>360.09</td><td>452.2</td><td>221.32</td><td>310</td><td>217.90</td><td>NA</td></tr></tbody></table>

Abbreviation: NA, not available.
For all the state numbers shown in Figures 1 and 2, the reference state for the nuclear power plant is chosen as a sub-cooled liquid at an ambient pressure of *p* = 1.01 bar and a temperature of 25 °C for the water/steam, and for the environmental and dead state desalination unit, these were *T* = 298 K, *p* = 1.01 bar and *X*<sub>f</sub> = 39,000 ppm, which represents the intake conditions for the seawater for this analysis. Therefore, *h* = 104.8374 kJ/kg and *s* = 0.3661 kJ/kg.K. The results of the energy and exergy flow rates for the combined PWR and ME-TVC-MED plant are presented in Appendix C. Table 7 shows that the reversible (maximum) work that can be extracted from the system is 3854.7 MWth. The total thermodynamic efficiency of the PWR-ME-TVC-MED cogeneration plant was found to be 35.38%. According to the literature, the thermodynamic efficiency value of a PWR-ME-TVC-MED cogeneration plant of 1000 MWe was estimated as 33.83% [7]. In addition, the thermal efficiencies of CANDU6-RO and SFR-RO hybrid systems are 32.8% and 36.8%, respectively [10]. In addition, the thermodynamic efficiency of the GT-MHR-RO cogeneration plant was found to be 41% [11]. Therefore, the cogeneration plant presented in this study has an average efficiency value among those.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** The simplified energy balance of the PWR-ME-TVC-MED (in the units of MWth).</figcaption></figure>

In addition, the total thermodynamic efficiency of the PWR-ME-TVC-MED cogenera-
tion was found to be 35.38%. Comparatively, a previous study of a 1000 MWe PWR-ME-
TVC-MED cogeneration plant found a thermodynamic efficiency of 33.83% [7]. Moreover,
the thermal efficiencies of the CANDU6-RO and SFR-RO hybrid systems are 32.8% and
36.8% [10], respectively. The thermodynamic efficiency of the GT-MHR-RO cogeneration
plant was found to be 41% [11]. Thus, the cogeneration plant presented in the present study
has an average efficiency value among those.

The total amount of energy that is lost as a result of inefficiencies in the cogeneration
system is estimated to be 64%. The overall thermodynamic efficiency of the PWR-ME-TVC-
MED system is 35.38%, as reported in Table 7. The difference between the total exergy rate produced and the total rate of exergy destruction is 0.046% for the entire system.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** The simplified exergy balance of the PWR-ME-TVC-MED (in the units of MWth).</figcaption></figure>

To address the component irreversibility in the PWR-ME-TVC-MED cogeneration system, the following measures can be taken:

*   **Improve the efficiency of the PWR:** The PWR operates by generating high-pressure steam, which is used to drive a turbine and generate electricity. The efficiency of the PWR can be improved by using advanced materials for the reactor core, optimizing the fuel assembly design, and improving the heat transfer in the steam generators.

*   **Reduce heat losses:** Heat losses occur in various components of the cogeneration plant, including the PWR, the ME unit, the TVC, and the MED unit. These losses can be minimized by improving the insulation of the components, reducing the surface area of the heat exchangers, and using advanced heat transfer fluids.

*   **Optimize the operation of the system:** The operation of the cogeneration plant can be optimized by adjusting the operating parameters of the various components to ensure that they operate at their maximum efficiency. This includes optimizing the flow rates of the fluids, adjusting the temperature and pressure of the steam and water, and minimizing the pressure drops across the components.

In addition, there are several ways to reduce heat losses in the PWR-ME-TVC-MED cogeneration plant, including:

*   **Insulation:** Insulating the plant's pipes, tanks, and other equipment can reduce heat loss and improve energy efficiency.
* *Heat recovery:* Implementing heat recovery systems can capture waste heat generated by the plant and use it for other processes, such as heating water or generating electricity. This can help reduce heat losses and improve the overall efficiency of the plant.

* *Improved maintenance:* Regular maintenance and cleaning of the plant’s equipment can help ensure that it operates efficiently and reduces heat loss due to leaks, blockages, or other issues.

* *Upgraded equipment:* Upgrading older or inefficient equipment with newer, more energy-efficient models can reduce heat losses and improve the overall efficiency of the plant.

<table><caption>**Table 7.** Results of the second-law analysis for the PWR-ME-TVC-MED cogeneration plant.</caption><thead><tr><th></th><th>Actual Work (MWth)</th><th>Reversible Work (MWth)</th><th>Rate of Irreversibility (MWth)</th><th>Irreversibility Percentage (%)</th><th>Exergy Efficiency (%)</th></tr></thead><tbody><tr><td><strong>Steam generation section</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Reactor</td><td>1916.3</td><td>3854.7</td><td>1938.6</td><td>50.28</td><td>49.98</td></tr><tr><td>Steam generators (SG)</td><td>0</td><td>147.08</td><td>147.08</td><td>3.81</td><td>92.98</td></tr><tr><td>Reactor coolant pumps</td><td>-16.36</td><td>-10.51</td><td>5.85</td><td>0.15</td><td>64.18</td></tr><tr><td>Total</td><td></td><td></td><td>2091.53</td><td>54.25</td><td></td></tr><tr><td><strong>Power production section</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>High-pressure turbine (HPT)</td><td>474.14</td><td>545.61</td><td>71.47</td><td>1.85</td><td>86.91</td></tr><tr><td>Low-pressure turbine (LPT)</td><td>889.96</td><td>1014.8</td><td>124.84</td><td>3.23</td><td>87.69</td></tr><tr><td>Re-heater (RH)</td><td>0</td><td>20.21</td><td>20.21</td><td>0.52</td><td>77.81</td></tr><tr><td>Moisture separation (MS)</td><td>0</td><td>2.75</td><td>2.75</td><td>0.06</td><td>91.44</td></tr><tr><td>Total</td><td></td><td></td><td>219.27</td><td>5.64</td><td></td></tr><tr><td><strong>Condensation section</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Condenser (C)</td><td>0</td><td>92.42</td><td>92.42</td><td>2.39</td><td>30.71</td></tr><tr><td>Total</td><td></td><td></td><td>92.42</td><td>2.39</td><td></td></tr><tr><td><strong>Preheating section</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>LP feed water heater-1</td><td>0</td><td>0.691</td><td>0.691</td><td>0.017</td><td>27.02</td></tr><tr><td>LP feed water heater-2</td><td>0</td><td>1.63</td><td>1.63</td><td>0.042</td><td>48.48</td></tr><tr><td>LP feed water heater-3</td><td>0</td><td>5.02</td><td>5.02</td><td>0.13</td><td>72.54</td></tr><tr><td>LP feed water heater-4</td><td>0</td><td>13.08</td><td>13.08</td><td>0.339</td><td>70.44</td></tr><tr><td>LP feed water heater-5</td><td>0</td><td>5.78</td><td>5.78</td><td>0.149</td><td>86.66</td></tr><tr><td>HP feed water heater-1</td><td>0</td><td>9.89</td><td>9.89</td><td>0.256</td><td>78.46</td></tr><tr><td>HP feed water heater-2</td><td>0</td><td>14.22</td><td>14.22</td><td>0.368</td><td>66.27</td></tr><tr><td>HP feed water heater-3</td><td>0</td><td>2.71</td><td>2.71</td><td>0.07</td><td>65.54</td></tr><tr><td>Feed water tank (FWT)</td><td>-9.56</td><td>-5.29</td><td>8.27</td><td>0.214</td><td>34.51</td></tr><tr><td>Feed water pumps</td><td>-2.23</td><td>-1.71</td><td>0.52</td><td>0.013</td><td>50</td></tr><tr><td>Main condensate pumps</td><td></td><td></td><td>64.96</td><td>1.685</td><td></td></tr><tr><td>Intermediate insulation loop</td><td>0</td><td>1.88</td><td>1.88</td><td>0.048</td><td>96.1</td></tr><tr><td>Heat exchanger (EX)</td><td></td><td></td><td>1.88</td><td>0.048</td><td></td></tr><tr><td><strong>Desalination unit</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Thermo-compressor</td><td>0</td><td>6.92</td><td>6.92</td><td>0.179</td><td></td></tr><tr><td>Evaporators</td><td>0</td><td>9.99</td><td>9.99</td><td>0.259</td><td>17.56</td></tr><tr><td>Condenser</td><td>0</td><td>1.42</td><td>1.42</td><td>0.036</td><td>21.66</td></tr><tr><td>Seawater Pump</td><td>-0.23</td><td>-0.16</td><td>0.170</td><td>0.004</td><td>55.12</td></tr><tr><td>Distillate pump</td><td>-0.16</td><td>-0.11</td><td>0.133</td><td>0.003</td><td>30</td></tr><tr><td>Brine pump</td><td>-0.10</td><td>-0.086</td><td>0.014</td><td>0.0003</td><td>62.3</td></tr><tr><td>Reject pump</td><td>-0.085</td><td>-0.053</td><td>0.062</td><td>0.001</td><td>70</td></tr><tr><td>Total</td><td></td><td></td><td>18.99</td><td>0.49</td><td>37.5</td></tr><tr><td>Net work output</td><td></td><td></td><td>1364.1</td><td>35.38</td><td></td></tr><tr><td>Imbalance</td><td></td><td></td><td>3.38</td><td>0.046</td><td></td></tr><tr><td>Total</td><td></td><td></td><td>3855</td><td>100</td><td></td></tr></tbody></table>

The results estimated by the current study are presented briefly in Table 8. Isentropic efficiencies for low- and high-pressure turbines, pumps for reactor coolant, feed water, main
condensate, seawater, distillate, brine, and reject are estimated as 87.68%, 86.9%, 64.24%, 76.68%, 55.33%, 69.56%, 68.75%, 86%, and 62.35%, respectively. The simplified energy and exergy balance diagrams of the PWR-ME-TVC-MED cogeneration plant are illustrated in Figures 3 and 4, respectively. From Figure 3, it is obvious that the core energy makes a cycle of the steam generator, preheating unit, power generation unit, and energy flow through condensation. In addition, there is a desalination unit. In the case of the exergy balance diagram in Figure 4, power generation, steam generation, condensation, preheating, and desalination units should be the key components.

<table><caption>**Table 8.** Results of the first- and second-law analysis of the cogeneration plant.</caption><thead><tr><th>Estimated Powers from the Current Study</th><th>Value</th><th>Unit</th></tr></thead><tbody><tr><td>Reactor thermal power</td><td>3854.7</td><td>MW</td></tr><tr><td>Turbine work</td><td>1364.1</td><td>MW</td></tr><tr><td>Reactor coolant pump work</td><td>-16.36</td><td>MW</td></tr><tr><td>Pump work on the secondary side</td><td>-13.79</td><td>MW</td></tr><tr><td>Heat transfer by the condenser</td><td>2436.54</td><td>MW</td></tr><tr><td>Steam generator thermal power</td><td>3871.1</td><td>MW</td></tr><tr><td>Net work of the secondary side</td><td>1350.31</td><td>MW</td></tr><tr><td>Heating steam of ME-TVC-MED</td><td>22.07</td><td>MW</td></tr><tr><td>Minimum work of separation</td><td>1.42</td><td>MW</td></tr><tr><td>Total rate of the irreversibility of ME-TVC-MED</td><td>20.65</td><td>MW</td></tr><tr><td>Total rate of irreversibility</td><td>2492.17</td><td>MW</td></tr><tr><td>Net heat input to the feed water</td><td>1363</td><td>MW</td></tr><tr><td>Total rate of irreversibility and net work of the work side</td><td>3840.91</td><td>MW</td></tr><tr><td><strong>Rates obtained in this study</strong></td><td></td><td></td></tr><tr><td>Total irreversibility rate</td><td>64.51</td><td>%</td></tr><tr><td>Thermodynamic efficiency of PWR-ME-TVC-MED</td><td>35.38</td><td>%</td></tr><tr><td>Total error rate</td><td>0.046</td><td>%</td></tr><tr><td>Exergetic efficiency of ME-TVC-MED</td><td>6.43</td><td>%</td></tr></tbody></table>

The efficiency of the cogeneration plant and the rates of exergy destruction of the components are illustrated in Figures 5 and 6, respectively, for the components of the desalination unit, which are more detailed than the thermodynamic calculations. It was found that the reactor core exhibited an irreversible rate of 50.28% in the installation. In addition, there are evaporators at an irreversible rate of 51% in the desalination unit.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Thermodynamic efficiency and irreversibility rate for the cogeneration plant.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Rates of exergy destruction in the cogeneration plant components.</figcaption></figure>

Figure 5 displays the exergy destruction and efficiency for all components as a proportion of total exergy destruction. The main exergy destruction occurs in the reactor core due to the fission process. The transfer between the fuel rods and coolant, heat loss, and pressure drop also contribute. The high-pressure turbine has a relatively high exergy destruction rate compared to other components, due to the high fuel exergy. The steam generator, which has many components and experiences heat transfer loss, temperature differences, and pressure differences, is considered the third largest source of irreversibilities and inefficiencies in the cycle. However, the exergy destruction of other components is minor and there is a difference in the level of exergy destruction between PWR and ME-TVC-MED, due to differences in heat capacity, the inlet temperature, and the type of working fluid. Approximately 64% of the total exergy generated from the fission power in a cogeneration plant is lost. Around 54% of this loss occurs in the primary loop (core and steam generators) of the PWR due to heat losses and pressure drops in the primary circuit, while the secondary circuit components (turbines, condensers, heaters) have a loss of approximately 10%. Moreover, other components, such as desalination and EX, have a loss of 0.5%, while 66.2% of the exergy is lost in the cogeneration (PWR-MED-TVC) [7]. In addition, 62% of the exergy is lost in the VVER NPP [1], and 53.7% of the exergy of the fuel is diminished in the PWR NPP [17].

The exergetic efficiency of the desalination unit was found to be 6.43%. According to the literature, the exergetic efficiency value of the desalination unit was found to be 7.58% [24]. The results indicate that in a cogeneration plant, about 64% of the total exergy produced from fission power is lost, with about 54% of this exergy loss occurring in the primary loop (core and steam generators) of the PWR plant due to heat losses and pressure drops in the primary circuit. The secondary circuit components (turbines, condenser, heaters) have a loss of approximately 10%. Furthermore, the distillation process is extremely energy-intensive, occurring at low pressure, which increases the amount of latent heat of vaporization and the percent of exergy destruction of the component desalination unit, as shown in Figure 6.

The results show that the highest source of irreversibilities within a desalination unit is found in the evaporators and that the thermo-compressor has a relatively high rate of exergy destruction compared to the remaining components in the system. The exergy destruction of other components is rather small. Increasing energy efficiency and reducing waste could be achieved by designing optimal heat recovery systems through various methods, such as heat exchangers or preheaters to extract the waste energy from the production line [25,26]. Significant efficiency improvements could also be obtained
by using nanofluids, as they have been found to exhibit enhanced thermal conductivity compared to traditional fluids [27–30].

**5. Conclusions**

In this paper, the energy and exergy efficiency of a PWR-ME-TVC-MED cogeneration plant for the simultaneous production of electrical energy and freshwater has been analyzed. The second law was used in defining exergy to reveal contributions of the components. The most recent validated data in the literature were used to determine the thermodynamic characteristics of seawater. A mathematical model was implemented in MATLAB software and validated based on a comparison with previously published research. The following are the key findings of this investigation:

It has been shown that the nuclear reactor core is the main source of exergy loss in the cogeneration plant. Indeed, it is the worse component in the whole cogeneration system due to the irreversibility of the energy transformation. Therefore, a substantial improvement can be attained in the performance of the cogeneration plant by considering the core of the reactor components. While the thermodynamic efficiency of the PWR-ME-TVC-MED is found to be 35.38%, the irreversibilities of the core of the reactor and turbines have been calculated as 50.28% and 5.8% respectively. Low irreversibilities also exist in the components of steam generators, condensers, moisture separators, and desalination units. Additionally, in the ME-TVC-MED desalination unit, the main exergy losses occur in effects and in the thermos-compressor, which represents 85% of the total value of exergy destruction.

One way to address component irreversibility is through the use of heat exchangers, which are devices designed to transfer heat between two fluids at different temperatures. In the PWR-ME-TVC-MED cogeneration plant, heat exchangers can be used to capture the waste heat from the power generation process and use it to generate steam for the thermal desalination process. This can significantly increase the overall efficiency of the plant, as the waste heat that would otherwise be lost is put to productive use. In addition, the cogeneration plant can incorporate other technologies and strategies to improve efficiency and reduce component irreversibility. For example, the plant can use advanced materials and designs to reduce heat loss and improve the performance of key components. The plant can also use advanced control systems to optimize the operation of the various components and minimize waste.

Taking all of these conclusions into consideration, optimizing the operating conditions of the existing designs leads to the introduction of various possibilities to improve the design and performance of the cogeneration plant, based on the analysis conducted. In general, any improvement in thermodynamic performance has an economic and environmental impact. This will be published soon as the second part of the current study, using exergo-economic and exergo-environmental methodologies.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/su15108358/s1.

**Author Contributions:** Conceptualization, Z.T., R.M., M.N.B., H.T., M.K., A.A., J.Z. and L.M.; Methodology, Z.T., R.M., M.N.B., H.T., M.K., A.A., J.Z. and L.M.; Software, Z.T., R.M., M.N.B., H.T., M.K., J.Z. and L.M.; Validation, Z.T., R.M., M.N.B., H.T., M.K., A.A., J.Z. and L.M.; Formal analysis, Z.T., R.M., M.N.B., H.T., M.K., A.A., J.Z. and L.M.; Investigation, Z.T., R.M., M.N.B., H.T., M.K., A.A., J.Z. and L.M.; Resources, Z.T., R.M., M.N.B., M.K., A.A., J.Z. and L.M.; Data curation, Z.T., R.M., M.K., A.A., J.Z. and L.M.; Writing—original draft, Z.T.; Writing—review & editing, M.N.B., H.T., A.A., J.Z. and L.M.; Visualization, Z.T., H.T., M.K., A.A., J.Z. and L.M.; Supervision, M.N.B., A.A. and J.Z.; Project administration, Z.T., A.A., J.Z. and L.M. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported in entire part by the Biomaterials and Transport Phenomena Laboratory agreement N° 303 03-12-2003, at the University of Medea. The authors acknowledge and gratefully thank the financial support provided by DG-RSDT of Algeria.
**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** All subjects gave their informed consent for inclusion before they participated in the study.

**Data Availability Statement:** The authors confirm that the data supporting the findings of this study are available within the article and its supplementary materials.

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

## Nomenclature

A: heat transfer area, m²

B: brine water flow rate, kg.s⁻¹

BPE: boiling point elevation, °C

Cₚ: specific heat, kJ.kg⁻¹.K⁻¹

CR: compression ratio

D: distillate, kg.s⁻¹

e: exergy, kJ.kg⁻¹

Ė: total energy flow rate, kW

Ėx: total exergy flow rate, kW

ER: expansion ratio

F: feed seawater flow rate, kg.s⁻¹

GOR: gain output ratio

h: specific enthalpy, kJ.kg⁻¹

İ: rate of irreversibility, kW

i: number of effects

LMTD: logarithmic mean temperature
difference, °C

m: mass flow rate, kg.s⁻¹

Mₑ: total seawater inlet flow rate, kg.s⁻¹

P: pressure, Pa, bar

Q: specific heat consumption, kJ.kg⁻¹

Q̇: net rate of heat transferred, kW

Ra: entrainment ratio

S: specific entropy, kJ.kg⁻¹.K⁻¹

Ṡg: entropy generation rate

T: temperature, °C

TBT: top brine temperature, °C

U: heat transfer coefficient, kW/m² K

V: velocity, m/s

Ẇ: net rate of work, kW

w: work per unit mass, kJ.kg⁻¹

X: salinity balance

Z: height, m

### Greek letters

ε: exergetic efficiency

λ: latent heat, kJ.kg⁻¹

η: isentropic efficiency

μ: chemical potential

Δ: difference

### Subscripts

0: reference state

a: property that corresponds to the fission
fragments' kinetic energy

c: coolant

cv: control volume

e: exit, outflow

f: feed water

fis: fission

g: generation

i: inlet, inflow, number of effects

ke: kinetic energy

m: motive steam

max: maximum

p: pump

pe: potential

ph: physical

R: reactor

RCP: coolant pump of reactor

Rej: rejection

Rev: reversible

s: steam

t: turbine

### Abbreviations

CANDU: Canadian deuterium uranium

cond: condenser

EX: heat exchanger

FWT: feed water tank

GA: genetic algorithm

GT-MHR: gas turbine-modular helium reactor

H.P: high-pressure

H.P.T: high-pressure turbine

L.P: low-pressure

L.P.T: low-pressure turbine

M.S: moisture separator

n, N: number of effects

NPP: nuclear power plant

P: pump

ppm: parts per million

PWR: pressurized water reactor

R.H: re-heater

RO: reverse osmosis

S.G: steam generators

SFR: sodium-cooled fast reactor

TVC: thermo-vapor compressor
# Appendix A

<table><caption>**Table A1.** Technical data of the nuclear power plant.</caption><thead><tr><td>Item</td><td>Value</td></tr></thead><tbody><tr><td>Coolant inlet temperature to the reactor, °C</td><td>291.7</td></tr><tr><td>Coolant outlet temperature from the reactor, °C</td><td>324.9</td></tr><tr><td>Coolant</td><td>H₂O</td></tr><tr><td>Coolant flow rate, kg/s</td><td>19,643.8</td></tr><tr><td>Average reactor pressure, bar</td><td>157</td></tr><tr><td>Number of coolant cycles</td><td>4</td></tr><tr><td>Inlet water temperature in the steam generator, °C</td><td>211.13</td></tr><tr><td>Outlet water temperature in the steam generator, °C</td><td>282.3</td></tr><tr><td>Inlet pressure at steam generator, bar</td><td>66.07</td></tr><tr><td>Outlet pressure at steam generator, bar</td><td>66.41</td></tr><tr><td>Steam flow rate in the generator, kg/s</td><td>2069.35</td></tr><tr><td>Number of coolant steam generators</td><td>4</td></tr><tr><td>Cooling temperature at condenser inlet, °C</td><td>25</td></tr><tr><td>Cooling temperature at condenser outlet, °C</td><td>35</td></tr><tr><td>Cooling water flow rate at the condenser, kg/s</td><td>61,728.34</td></tr></tbody></table>

# Appendix B

<table><caption>**Table A2.** Technical data of the desalination unit.</caption><thead><tr><td>Description</td><td>Value</td><td>Unit</td></tr></thead><tbody><tr><td>Cooling water temperature, T<sub>cw</sub></td><td>25</td><td>°C</td></tr><tr><td>Feed sea water temperature, T<sub>f</sub></td><td>35</td><td>°C</td></tr><tr><td>Number of effects, n</td><td>8</td><td>-</td></tr><tr><td>Motive steam flow rate, S</td><td>17.5</td><td>kg/s</td></tr><tr><td>Maximum brine temperature, T<sub>1</sub></td><td>65</td><td>°C</td></tr><tr><td>Minimum brine temperature, T<sub>n</sub></td><td>40</td><td>°C</td></tr><tr><td>Motive pressure, P<sub>m</sub></td><td>101.31</td><td>kPa</td></tr><tr><td>Feed water salinity, X<sub>f</sub></td><td>39,000</td><td>Ppm</td></tr></tbody></table>

# Appendix C

<table><caption>**Table A3.** Total energy and exergy flow rates for the cogeneration plant.</caption><thead><tr><td>State No.</td><td>Total Energy Flow Rate (MW)</td><td>Total Exergy Flow Rate (MW)</td><td>State No.</td><td>Total Energy Flow Rate (MW)</td><td>Total Exergy Flow Rate (MW)</td></tr></thead><tbody><tr><td>1</td><td>29,118.7</td><td>8868.3</td><td>36</td><td>44.1</td><td>5</td></tr><tr><td>2</td><td>25,247.6</td><td>6941.7</td><td>37</td><td>47</td><td>3.6</td></tr><tr><td>3</td><td>25,264</td><td>6952.2</td><td>38</td><td>16.1</td><td>4</td></tr><tr><td>4</td><td>5739</td><td>2148.5</td><td>39</td><td>22.6</td><td>0.3</td></tr><tr><td>5</td><td>5383.5</td><td>2014.4</td><td>40</td><td>33</td><td>0.3</td></tr><tr><td>6</td><td>4489.1</td><td>1331</td><td>41</td><td>30,509.6</td><td>0</td></tr><tr><td>7</td><td>182</td><td>32.4</td><td>42</td><td>32,944.2</td><td>37.5</td></tr><tr><td>8</td><td>4045.2</td><td>1217.1</td><td>43</td><td>44</td><td>4.9</td></tr><tr><td>9</td><td>350.3</td><td>129.6</td><td>44</td><td>186</td><td>48.7</td></tr><tr><td>10</td><td>151.4</td><td>39</td><td>45</td><td>23.9</td><td>3.7</td></tr><tr><td>11</td><td>4239</td><td>1287.6</td><td>46</td><td>95.5</td><td>25</td></tr><tr><td>12</td><td>420.3</td><td>137.9</td><td>47</td><td>125.6</td><td>0</td></tr><tr><td>13</td><td>261.7</td><td>78.9</td><td>48</td><td>125.7</td><td>0.06</td></tr><tr><td>14</td><td>90.4</td><td>23.7</td><td>49</td><td>176.1</td><td>0.922</td></tr><tr><td>15</td><td>221.3</td><td>49.7</td><td>50</td><td>75.6</td><td>0.386</td></tr><tr><td>16</td><td>265.3</td><td>45.9</td><td>51</td><td>100.5</td><td>0.526</td></tr><tr><td>17</td><td>70.1</td><td>5.4</td><td>52</td><td>58.6</td><td>1.374</td></tr></tbody></table>
<table><caption>Table A3. *Cont.*</caption><thead><tr><th>State No.</th><th>Total Energy Flow Rate (MW)</th><th>Total Exergy Flow Rate (MW)</th><th>State No.</th><th>Total Energy Flow Rate (MW)</th><th>Total Exergy Flow Rate (MW)</th></tr></thead><tbody><tr><td>18</td><td>2606.4</td><td>123.1</td><td>53</td><td>58.6</td><td>1.446</td></tr><tr><td>19</td><td>200.2</td><td>01</td><td>54</td><td>36.2</td><td>0.832</td></tr><tr><td>20</td><td>202.5</td><td>2.1</td><td>55</td><td>75.6</td><td>0.41</td></tr><tr><td>21</td><td>201.9</td><td>2.7</td><td>56</td><td>75.5</td><td>0.37</td></tr><tr><td>22</td><td>208.4</td><td>1.7</td><td>57</td><td>56.3</td><td>4.723</td></tr><tr><td>23</td><td>233.4</td><td>3.3</td><td>58</td><td>47.6</td><td>10.977</td></tr><tr><td>24</td><td>288</td><td>7</td><td>59</td><td>102.8</td><td>12.239</td></tr><tr><td>25</td><td>505.2</td><td>36.2</td><td>60</td><td>56.3</td><td>4.723</td></tr><tr><td>26</td><td>682.5</td><td>75.2</td><td>61</td><td>47.6</td><td>10.977</td></tr><tr><td>27</td><td>726.6</td><td>80.2</td><td>62</td><td>102.8</td><td>12.239</td></tr><tr><td>28</td><td>1228</td><td>157.9</td><td>63</td><td>51.1</td><td>2.367</td></tr><tr><td>29</td><td>1239.5</td><td>161.2</td><td>64</td><td>4.9</td><td>0.252</td></tr><tr><td>30</td><td>1522.4</td><td>248.6</td><td>65</td><td>4.9</td><td>0.252</td></tr><tr><td>31</td><td>1839.7</td><td>358.7</td><td>66</td><td>3.2</td><td>0.08</td></tr><tr><td>32</td><td>1867.9</td><td>368.9</td><td>67</td><td>55.5</td><td>1.405</td></tr><tr><td>33</td><td>123.2</td><td>26</td><td>68</td><td>55.4</td><td>1.398</td></tr><tr><td>34</td><td>226.2</td><td>39.5</td><td>69</td><td>35</td><td>0.959</td></tr><tr><td>35</td><td>387</td><td>53.6</td><td>70</td><td>52</td><td>1.321</td></tr></tbody></table>

## References

1. Terzi, R.; Turkenmez, I.; Kurt, E. Energy and Exergy Analyses of a VVER Type Nuclear Power Plant. *Hydrog. Energy* **2016**, *41*, 12465–12476. [CrossRef]

2. Sahin, S.; Sahin, H.M.; Al-Kusayer, T.A.; Sefidvash, F. An Innovative Nuclear Reactor for Electricity and Desalination. *Energy Res.* **2011**, *35*, 96–102. [CrossRef]

3. Alonso, G.; Del Valle, E.; Ramirez, J.R. *Desalination in Nuclear Power Plants*; Woodhead Publishing: Sawston, UK, 2020. [CrossRef]

4. Terzi, R. Application of Exergy Analysis to Energy Systems. In *Application of Exergy*; IntechOpen: London, UK, 2018. [CrossRef]

5. Eshun, R.B. Energy and Exergy Based Performance Analysis of Westinghouse AP1000 Nuclear Power Plant. *Adv. Appl. Sci.* **2019**, *4*, 1–10. [CrossRef]

6. Khoshgoftar Manesh, M.H.; Amidpour, M. Multi-Objective Thermoeconomic Optimization of Coupling MSF Desalination with PWR Nuclear Power Plant Through Evolutionary Algorithms. *Desalination* **2009**, *249*, 1332–1344. [CrossRef]

7. Khoshgoftar Manesh, M.H.; Amidpour, M.; Hamedi, M.H. Optimization of the Coupling of Pressurized Water Nuclear Reactors and Multistage Flash Desalination Plant by Evolutionary Algorithms and Thermoeconomic Method. *Energy Res.* **2009**, *33*, 77–99. [CrossRef]

8. Ansari, K.; Sayyaadi, H.; Amidpour, M. A Comprehensive Approach in Optimization of a Dual Nuclear Power and Desalination System. *Desalination* **2011**, *269*, 25–34. [CrossRef]

9. Ansari, K.; Sayyaadi, H.; Amidpour, M. Thermoeconomic Optimization of a Hybrid Pressurized Water Reactor (PWR) Power Plant Coupled to a Multi Effect Distillation Desalination System with Thermo-Vapor Compressor (MED-TVC). *Energy* **2010**, *35*, 1981–1996. [CrossRef]

10. Khalid, F.; Dincer, I.; Rosen, M.A. Comparative Assessment of CANDU 6 and Sodium-Cooled Fast Reactors for Nuclear Desalination. *Desalination* **2016**, *379*, 182–192. [CrossRef]

11. Khalid, F.; Dincer, I.; Rosen, M.A. Analysis and Assessment of a Gas Turbine-Modular Helium Reactor for Nuclear Desalination. *Nucl. Eng. Radiat. Sci.* **2016**, *2*, 031014. [CrossRef]

12. Priego, E.; Alonso, G.; Valle, E.D.; Ramirez, R. Alternatives of steam extraction for desalination purposes using SMART reactor. *Desalination* **2017**, *413*, 199–216. [CrossRef]

13. Al-Othman, A.; Darwish, N.N.; Qasim, M.; Tawalbeh, M.; Darwish, N.A.; Hilal, N. Nuclear desalination: A state-of-the-art review. *Desalination* **2019**, *457*, 39–61. [CrossRef]

14. Marques, J.G.O.; Costa, A.L.; Pereira, C.; Fortini, Â. Energy and Exergy Analyses of Angra 2 Nuclear Power Plant. *Radiat. Sci.* **2019**, *7*, 1–18. [CrossRef]

15. Binamer, A. Second law and sensitivity analysis of large ME-TVC desalination units. *Desalination Water Treat.* **2016**, *53*, 1234–1245. [CrossRef]

16. Sayyaadi, H.; Sabzaligol, T. Exergoeconomic optimization of a 1000 MW light water reactor power generation system. *Energy Res.* **2009**, *33*, 378–395. [CrossRef]

17. Sayyaadi, H.; Sabzaligol, T. Various approaches in optimization of a typical pressurized water reactor power plant. *Appl. Energy* **2009**, *86*, 1301–1310. [CrossRef]
18. El-Dessouky, H.T.; Ettouney, H.M. *Fundamental of Salt-Water Desalination*; Elsevier Science: Amsterdam, The Netherlands, 2002; ISBN 0444508104.

19. Durmayaz, A.; Yavuz, H. Exergy analysis of a pressurized-water reactor nuclear power plant. *Appl. Energy* **2001**, *69*, 39–57.
[CrossRef]

20. Sharqawy, M.H.; Lienhard, J.H.; Zubair, S.M. Thermophysical properties of seawater: A review of existing correlations and data. *Desalination Water Treat.* **2010**, *16*, 354–380. [CrossRef]

21. Sharqawy, M.H.; Lienhard, J.H.; Zubair, S.M. On exergy calculations of seawater with applications in desalination systems. *Therm. Sci.* **2011**, *50*, 187–196. [CrossRef]

22. Darwish, M.A.; Al-juwayhel, F.; Abdulraheim, H.K. Multi-effect boiling systems from an energy viewpoint. *Desalination* **2006**, *194*, 22–39. [CrossRef]

23. AL-Habshi, S. Simulation and Economic Study of the MED-TVC Units at Umm Al-Nar Desalination Plant. Master's Thesis, UAE. 2002. Available online: https://scholarworks.uaeu.ac.ae/all_theses (accessed on 17 May 2023).

24. Almutairi, A.; Pilidis, P.; Al-mutawa, N.; Al-weshahi, M. Energetic and exergetic analysis of cogeneration power combined cycle and ME-TVC-MED water desalination plant: Part-1 operation and performance. *Appl. Therm. Eng.* **2016**, *103*, 77–91. [CrossRef]

25. Eshoul, N.; Almutairi, A.; Lamidi, R.; Alhajeri, H.; Alenezi, A. Energetic, Exergetic, and Economic Analysis of MED-TVC Water Desalination Plant with and without Preheating. *Water* **2018**, *10*, 305. [CrossRef]

26. Menasri, R.; Triki, Z.; Bouaziz, M.N.; Hamrouni, B. Energy and exergy analyses of a novel multi-effect distillation system with thermal vapour compression for seawater desalination. *Desalin. Water Treat.* **2022**, *246*, 54–67. [CrossRef]

27. AbdElhafez, S.E.; El-Shazly, A.H.; El-Maghraby, A. Improving the rate of solar heat recovery from nanofluids by using micro heat exchanger. *Desalin. Water Treat.* **2016**, *57*, 48–49. [CrossRef]

28. Ali, A.Y.M.; El-Shazly, A.H.; Elkady, M.F.; AbdElhafez, S.E. Evaluation of Surfactants on Thermo-Physical Properties of Magnesia-Oil Nanofluid. *Mater. Sci. Forum* **2018**, *928*, 106–112. [CrossRef]

29. AbdElhafez, S.E.; Abo-Zahhad, E.M.; El-Shazly, M.; El-Shazly, A.H.; El-Kady, M.F. Experimental Investigation of Heat Transfer for Graphene/Water Nanofluid in Micro Heat Exchanger. *AIP Conf. Proc.* **2017**, *1814*, 020014. [CrossRef]

30. Singh, T.; Atieh, M.A.; Al-Ansari, T.; Mohammad, A.W.; McKay, G. The Role of Nanofluids and Renewable Energy in the Development of Sustainable Desalination Systems: A Review. *Water* **2020**, *12*, 2002. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.