Thermo-Economic Analysis on Integrated CO 2 , Organic Rankine Cycles, and NaClO Plant Using Liqueﬁed Natural Gas

: The thermal energy conversion of natural gas (NG) using appropriate conﬁguration cycles represents one of the best nonrenewable energy resources because of its high heating value and low environmental effects. The natural gas can be converted to liqueﬁed natural gas (LNG), via the liquefaction process, which is used as a heat source and sink in various multigeneration cycles. In this paper, a new conﬁguration cycle is proposed using LNG as a heat source and heat sink. This new proposed cycle includes the CO 2 cycle, the organic Rankine cycle (ORC), a heater, a cooler, an NaClO plant, and reverse osmosis. This cycle generates electrical power, heating and cooling energy, potable water (PW), hydrogen, and salt all at the same time. For this purpose, one computer program is provided in an engineering equation solver for energy, exergy, and thermo-economic analyses. The results for each subsystem are validated by previous researches in this ﬁeld. This system produces 10.53 GWh electrical energy, 276.4 GWh cooling energy, 1783 GWh heating energy, 17,280 m 3 potable water, 739.56 tons of hydrogen, and 383.78 tons of salt in a year. The proposed system energy efﬁciency is 54.3%, while the exergy efﬁciency is equal to 13.1%. The economic evaluation showed that the payback period, the simple payback period, the net present value, and internal rate of return are equal to 7.9 years, 6.9 years, 908.9 million USD, and 0.138, respectively


Introduction
Energy demand is envisaged to increase by approximately 50% by 2070 [1].Whilst renewable energy sources are exponentially growing, the stable high-power generation of nonrenewable methods is still required for industrialization.Due to the shortage and the environmental impacts of nonrenewable energy resources, the development of new methods of energy utilization is necessary [2].For achieving this target, the use of natural gas (NG) is a prime candidate due to its high lower heating value (LHV), high H/C ratio (hydrogen/carbon), and low environmental pollution, making it the cleanest fossil fuel [3,4].For easier transportation and storage, NG liquefaction is an appropriate approach, especially for remote regions.
LNG is a liquid state of natural gas which has a maximum amount of methane (CH 4 ), with other components such as ethane (C 2 H 6 ) and propane (C 3 H 8 ).This mixture is cooled to below 113.1 K, changing its state from gas to liquid through a cryogenic process, thereby facilitating shipping and storage.The volume of LNG is 600 times smaller than that of natural gas.Upon reaching its destination, LNG is reverted to natural gas and distributed via pipelines to the final consumers, comprising industries, residential buildings, and commercial buildings [5].
LNG can be used for electrical power production as a heat source or sink.For heat source applications, it can be used as fuel upon conversion to NG.Moreover, it can be used as a heat sink to absorb the cycle's heat dissipation via a condenser [15].

Literature Review
The thermodynamic performance of the hybrid integration of a two-stage ORC with R116 and R227ea as the working fluid, driven by low-grade heat and LNG as the cryogenic energy storage, was studied.The results of this study stated that a two-stage ORC system was able to produce 1776.44 kW of power with 25.64% thermal efficiency and 31.02%exergy efficiency.Moreover, the cost of power production for this system was 6.3 USD/W, since the LNG temperature was about 283.15 K [24].
The CO 2 cryogenic flue gas capture of an LNG power generation system was investigated.In this system, the flue gas was compressed to remove CO 2 from the exhaust gas.The results of this study showed that this system could obtain a 90% CO 2 recovery rate [25].
The performance of a combined gas system coupled with cold energy released during the regasification process of LNG was investigated.The results showed that the electrical efficiency of this system could be increased by about 8% [26].
The performance of an advanced thermal power plant was compared with a conventional combined power plant.In this comparison, the usage of intercooling and inlet air cooling was investigated.Furthermore, the latent heat of steam from a steam turbine and the heat rejected from compressed air were used for power generation and to heat the LNG.The results showed that the overall output work and electrical efficiency of the combined cycle increased by 2.8% up to 76.8 MW [27,28].
An ammonia/water-based combined power cooling cycle with LNG as cold storage, driven by low-temperature waste heat, was employed for electrical generation and as a district cooling system in coastal hot climate regions.At 70 • C heat source temperature, the system's energy efficiency was 39% and the system's exergy efficiency was 36% [29].
A novel system comprising a combination of biomass-based gas and steam cycles, a cascade ORC, an absorption refrigeration system, a PEM fuel cell, and LNG as a thermal sink was applied for power, cooling, and hydrogen production.The results showed that the mass flow rate of the fuel was the most important factor affecting the energy efficiency and the total cost.As a result, an increase of 4 kg•s −1 to 10 kg•s −1 in the mass flow rate of the fuel caused a decrease of 8.5% in the total cost and an increase of 122.8% in the system's energy efficiency [30].
A mathematical model for the integration of a solid oxide fuel cell, gas and CO 2 cycles, and ORC using LNG cold energy storage was applied to achieve cascade energy exploitation and CO 2 capture.The results showed that an energy efficiency of 79.48% and an exergy efficiency of 62.3% were obtained [31].
A summary of the studies related to LNG-powered multigeneration systems is presented in Table 1.The system thermal efficiency was equal to 6.75% The net output power was 108.7 kW [26] The novelty of this study, with respect to the references discussed in Table 1, is the investigation of a system that includes the ORC, the CO 2 cycle, RO, and an NaClO plant using LNG as a heat source and sink.Furthermore, the proposed cycle has various products, i.e., electrical power, heating and cooling, potable water (PW), hydrogen, and sodium chloride (NaCl), which are produced simultaneously when the system operates.
In this paper, integration of the CO 2 cycle and ORC is proposed using LNG as a heat source and heat sink.The LNG absorbs the dissipated heat of condensers of both cycles (heat sink).Then, it expands in turbine III for electrical production.Moreover, it is used in the cooler for cooling load production.At this stage, the LNG is converted to NG, and it reacts with air in the boiler to warm the water and convert it to steam.The produced steam is used in the steam generator and evaporator of the CO 2 cycle and ORC to meet the heat needs of both cycles.Additionally, it produces heating load in the heater.Some of the electrical power in both cycles is consumed in the RO and NaClO plants to simultaneously produce potable water (PW), hydrogen, and NaCl.
The novelties of this paper are as follows: • A new multigeneration system is proposed using LNG as a heat source and sink; • Energy, exergy, and economic investigations of the new configuration system are carried out; • Various products are generated, i.e., electrical power, heating and cooling loads, PW, hydrogen, and NaCl.

Materials and Methods
Figure 1 shows a schematic diagram of the cycle.The working fluids in the CO 2 cycle and the organic Rankine cycle are carbon dioxide and R134a, respectively.The physical and thermodynamic properties of tetrafluoroethene (R134a) are similar to R12 with lower ozone depletion.This cycle is powered by LNG.Furthermore, the LNG is used as a heat sink.The LNG, after pressurizing in pump I (points 5, 6), absorbs the dissipated heat from the condensers of the CO 2 cycle and ORC (points 6, 7, 8).Then, it is expanded in turbine III and the generator to produce electricity (points 8, 9), while it is used as a coolant fluid in the cooler (points 9, 10), where it is also converted to NG (point 10).The NG is burned in the boiler to warm the water and convert it to steam (points 16,17).The steam provides the heat for the steam generator (CO 2 cycle) and the evaporator of the ORC (points 17,18,36,37).Some of the steam is used in the heater to produce heating energy (points 19,20).The exhaust hot water of the ORC evaporator and heater is mixed, and it is transferred and pressurized in pump IV (points 15,16,20,37).Some of the electricity produced in this system is consumed in the RO system to produce PW.In the RO system, seawater goes throughout pumps V and VI, and it is pressurized (points 21,22,23,24,25).The pressurized seawater is transferred to membranes I and II, and it is separated into brain water (points 27, 29) and potable water (points 26,28,30,31).The brain water rotates the recovery turbine to recover energy and produce electrical power.Then, it is transferred to the NaClO plant, where it produces NaCl and hydrogen.
In summary, in this cycle, four components produce electricity (turbines I, II, III, and the recovery turbine), while six components (pumps I to VI) and the NaClO plant consume electricity.Figure 2 shows the energy flow between the subsystems of the proposed system.The following assumptions are made in this study [14,21,[32][33][34]: (f) The turbine and pump polythrophic efficiencies are assumed to be 80%; (g) The heat exchanger effectiveness factor is assumed to be 80%; (h) The salt concentration in the electrolyzer is assumed constant; (i) The inlet CO 2 and LNG of the pump are in a liquid state; (j) The RO recovery ratio is 0.3.

Mass and Energy Balance
In general, the mass and energy balance equations are described as shown below [35]. .
W and .
Q are the power and heat transfer rate, while .m and h are the mass flow rate and enthalpy, respectively.Subscripts R, f, P, and 0 denote reactant, formation, product, and dead state, respectively.
The mass and energy balance relationships for the CO 2 cycle, ORC, LNG and water loops, and boiler are shown in Table 2 [23,36,37].  .
The mass balance equations in the RO system are shown below [14,38,39]. . .
where x is the salt concentration.Subscripts PW, SW, and BW denote potable water, seawater, and brain water, respectively.The sea and potable water relationship is shown below [14,38]. .
where RR denotes the recovery ratio.
The osmosis pressure for the sea, potable, and brain water streams is shown below [14,38].
where R defines the universal gas constant.The membrane net pressure equation is shown below [14,38].
The RO pump net pressure is calculated as shown below [14,38].
where A m is the membrane area.
The power consumption of the RO pump can be calculated as shown below [14,38]. .
where ρ is the density.The mass, energy, and concentration balance relationships for the RO system are shown in the Table 3.The following reaction can be considered in the NaClO plant: In the NaClO plant, the following relationships can be considered [14,38]: x NaClO = 1 6 x BW .
The NaClO plant concentration, mass, and energy balance relationships are shown in Table 4.The electrical power production of the system is calculated as follows: .
The system energy efficiency is expressed as

Exergy Analysis
The specific exergy is expressed below, neglecting the kinetic and potential exergies [40,41].
where x and e are the mass fractions and specific exergy, respectively, T, h, and s are the temperature, specific enthalpy, and specific entropy, respectively, and y denotes the mole fraction.Subscript ch denotes chemical.Abbreviations i and 0 represent the species and the dead state condition.Table 5 shows the exergy destruction rates and exergy efficiencies for the system components [23,42,43].
The system exergy efficiency is expressed below.

Thermo-Economic Analysis
The multigeneration annual income and outcome are calculated as follows [44,45]: where c is the specific cost, and A is the annual capacity of system production.The costs of products are shown in Table 6.The system investment cost can be calculated as shown below [44,45].
where C is the investment cost of each component, as shown in Table 7.The operation and maintenance cost is assumed to be 3% of the initial cost [51][52][53][54].NaClO NaClO (HD:6000) 45,000 [14] In Table 7, A, H, D, and V denote the surface area, head of pump or turbine, diameter, and volume, respectively.
The logarithmic method is applied to calculate the surface area of the heat exchanger, as shown below [51]. .
Q is the heat transfer rate, ∆T In is the logarithmic mean temperature difference, A and U denote the surface area and overall heat transfer coefficient, respectively, and F is the correction factor.The overall heat transfer coefficient values for components are shown in Table 8 [52].The effect of inflation can be seen in the following relationship [64]: where i and n are the inflation rate (3%) and the number of years, respectively [65].
The simple payback period (SPP) can be calculated as shown below [44,45].
The internal rate of return (IRR) can be calculated as follows [55,67,68]:

Simulation Method Description
For the simulation, a computer code provided in the engineering equation solver (EES) software was used.The computer code input information is shown in Table 9.

Model Validation
Since the proposed plant had a new configuration, its complete validation was not achievable.Hence, each subsystem was considered.To validate the CO 2 cycle and LNG loop, the results in [32] were considered.Four key parameters (CO 2 and LNG turbine power production and pump power production) were compared.Table 10 shows a comparison between the present work and the results in [32].The source of error was related to the calculation of thermodynamic properties, which features slight differences in each software.To validate the ORC results, Figure 2 in [56] was considered.R600a was selected as the working fluid.The data for the ORC presented in [56] were used.Figure 3 shows a comparison between the net output power of this work and the results in [56].The maximum error was 3.5%, which is compatible with engineering calculations.The source of error was similar to that of the CO 2 cycle.
For RO validation, the results in [14] were used.Table 11 shows the comparison.The error ranged from 0.7% to 7% due to some information about membrane performance not being exactly clarified in [14] For NaClO plant validation, the results in [57] were used.The electrical power consumption of the NaClO plant was found to be 4 kW, while this value was theoretically calculated as 3.78 kW by the computer code.The error of 5.5% was related to the type and model of the NaClO plant.

Energy and Exergy Analyses
Table 12 shows the thermodynamic properties at every point of the system.Table 13 shows the values of the system products.This system produces 10.53 GWh of electrical energy, 276.4 GWh of cooling energy, 1783 GWh of heating energy, 17,280 m 3 of potable water, 739.56 tons of hydrogen, and 383.78 tons of salt in a year.The ratios of cooling and heating energies to electrical energy are 26.25 and 169.33, respectively.Figure 4 shows the system energy and exergy efficiencies.According to the figure, the energy and exergy efficiencies are 54.3% and 13.1%, respectively.The system energy efficiency is much higher than the system exergy efficiency, as, in the exergy efficiency equation (Equation ( 20)), both the heat and the cooling production rates are multiplied by the ratios, which are lower than the unit.The variation in RO and NaClO plant power consumption and in recovery turbine and system power productions versus seawater inlet mass flow rate ( ṁ21 ) is shown in Figure 5. Upon increasing the ṁ21 from 1 to 5 kg/s, the power consumption of the RO and NaClO plant can be increased from 1.486 and 134.8 kW to 35.7 and 673.9 kW, respectively.According to the equations presented in Table 3, the power consumption of the RO pumps has a linear relationship with the mass flow rate of feed water.According to Equation ( 16), the power consumption of the NaClO plant also has a linear relationship with feed brain water.The net electrical power production of the total proposed system is reduced from 1456.7 kW to 876.7 kW following this increase in ṁ21 .
The system energy and exergy efficiencies versus seawater mass flow rate ( ṁ21 ) are shown in Figure 6.Increasing the seawater mass flow rate did not have a major effect on the system energy and exergy efficiencies.Upon increasing the seawater mass flow rate from 1 to 5 kg/s, the system energy and exergy efficiencies were reduced from 54.25% and 13.11% to 54.23% and 13.01%, respectively.It can be concluded that the system energy and exergy efficiencies remained constant.Since the extra electrical power consumed by the NaClO and RO systems is compensated for by the extra energy and exergy flows of PW, NaCl, and hydrogen, according to Equations ( 18) and (20).Figure 7 shows the changes in pump electrical power consumption and turbine electrical power production, as well as the net output power of the CO 2 cycle, with the variation in outlet pressure of the pump (P 2 ) in the CO 2 cycle.Upon increasing the outlet pressure of the pump, the electrical power consumption of the pump in the CO 2 cycle is increased.Moreover, the electrical power production of the turbine is decreased.This is because the superheat temperature at point 3 is fixed at 493.1 K (above the supercritical temperature), which is about 293.1 K less than the superheated steam temperature at point 17 according to pinch analysis.Thus, by increasing the outlet pressure of the pump in the CO 2 cycle (P 2 ) at a fixed temperature (493.1 K), the enthalpy at point 3 is reduced, and the electrical power production in the turbine is reduced.It is clear that, by increasing the electrical power consumption of the pump and the electrical power production of the turbine, the net electrical power of the CO 2 cycle is decreased.
Figure 8 shows the variation in CO 2 cycle energy and exergy efficiencies with outlet pump pressure in the CO 2 cycle.According to Figure 7, the CO 2 cycle net output electrical power is reduced upon increasing the outlet pressure of the pump, whereas the inlet heat rate of the CO 2 cycle ( ṁ (h 17 -h 18 )) remains constant.
The change in net output power of the ORC with pump outlet pressure at different inlet turbine temperatures is shown in Figure 9. Upon increasing the ORC pump outlet pressure, the cycle net output power is reduced, since it decreases the enthalpy value at point 13.Furthermore, the power consumption of the pump increases.Upon increasing the turbine inlet temperature in the ORC, the net output electrical power is increased by increasing the enthalpy at point 13.
The contribution to exergy destruction rate of each subsystem is presented in Figure 10.The highest contribution comes from the water loop, which includes a heater, burner/boiler, and pump IV.In the burner, a chemical reaction (combustion) warms the water in the boiler, thereby converting water to steam.Thus, this loop has the highest exergy destruction rate.In distant second place, the CO 2 cycle has a considerable exergy destruction rate due to the heat transfer taking place in the steam generator and condenser.The lowest exergy destruction rate occurs in the RO system, as the pressures and temperatures in this system are all close to the dead state.

Thermo-Economic Analyses
Table 14 shows the key economic parameters for the system.The system's NPV is 908.9 million USD, the system's PP and SPP are equal to 7.9 and 6.9 years, and the system's IRR is 0.138.Figure 11 shows the variation in PP and SPP versus feed seawater mass flow rate.Upon increasing the seawater mass flow rate, both PP and SPP are reduced, which is beneficial.Specifically, upon increasing the feed seawater mass flow rate from 1 kg/s to 5 kg/s, the PP is reduced from 8.7 to 6.2 years and the SPP is reduced from 7.56 to 5.59 years.
In general, increasing the seawater mass flow causes the following effects: (a) Increase in the initial cost of the NaClO and RO plants (negative effect); (b) Increase in the system product costs (NaCl, PW, and hydrogen) due to an increase in these products; (c) Decrease in the electrical power product cost due to an increase in RO and NaClO plant power consumption.
The magnitude of the second effect exceeds that of the other two effects; thus, both PP and SPP are reduced upon increasing the feed seawater mass flow rate.The variation in system NPV with the mass flow rate of seawater is depicted in Figure 12.The NPV of the system is increased upon increasing the mass flow rate of seawater.Similar to Figure 11, it is clear that, if more electrical power is consumed to produce PW, NaCl, and hydrogen, it is more beneficial.

Conclusions
A new integrated system was proposed in this paper which simultaneously uses LNG as a heat source and sink.Through this method, the ratio of the heat sink and source temperatures is high, thereby improving system performance.
The CO 2 cycle and ORC are combined to generate electricity in this proposed system.In both cycles, LNG is used to absorb the heat dissipated by condensers, as well as to produce cooling load.Subsequently, the LNG is converted to NG.The NG is used to convert water into steam in the burner and boiler, thus meeting the energy needs of the steam generator and evaporator in the CO 2 cycle and the ORC, in addition to producing heating load.A portion of the electrical power generated in the two mentioned cycles is used to produce hydrogen, salt, and PW in the RO and NaClO plants.The main conclusions are summarized below.
(a) In comparison with a system featuring LNG only as a heat sink [14], which uses solar energy through a flat plate collector as the heat source, the system energy and exergy efficiencies are further improved from 12.4% and 4.45% to 54% and 13.1%, respectively.(b) The NPV of this system is equal to 908.9 million USD.(c) The PP and SPP of this system are 7.9 and 6.9 years, respectively.(d) The IRR value of this system is equal to 0.138.(e) Greater desalination of seawater to produce PW, salt, and hydrogen is beneficial according to the economic evaluation.(f) Greater seawater desalination does not have a considerable effect on the system energy and exergy efficiencies.(g) The highest and lowest contributions to the exergy destruction rates were presented by the water loop and RO system.
In the future, the integration of various cycles such as the Kalina and Goswami cycles can be examined as alternatives.

Figure 1 .
Figure 1.Schematic diagram of the system.

Figure 2 .
Figure 2. Energy flow between the subsystems.
(a) The system works in steady-state conditions; (b) The ambient pressure and temperature are 288 K and 1 bar, respectively; (c) Pressure loss in the heat exchanger is assumed to be 2%; (d) The kinetic and potential energies are ignored; (e) The pressure loss in the cycles is ignored;

Figure 3 .
Figure 3.Comparison of the net output power in the present work with that in [14] versus the evaporator temperature of the ORC.

Figure 4 .
Figure 4.The system energy and exergy efficiencies.

Figure 5 .
Figure 5.The variation in RO, NaClO plant, recovery turbine, and system power with inlet seawater mass flow rate ( ṁ21 ).

Figure 6 .
Figure 6.The system energy and exergy efficiencies versus seawater mass flow rate ( ṁ21 ).

Figure 7 .
Figure 7.The changes in pump and turbine electrical power consumption/production and net output power of the CO 2 cycle with the variation in pump outlet pressure (P 2 ) in the CO 2 cycle.

Figure 8 .
Figure 8.The variation in CO 2 cycle energy and exergy efficiencies with outlet pump pressure in the CO 2 cycle.

Figure 9 .
Figure 9.The variation in net output power of the ORC with pump outlet pressure at different inlet turbine temperatures.

Figure 10 .
Figure 10.Contributions of each system to the exergy destruction rate.

Figure 11 .
Figure 11.The variation in PP and SPP versus feed seawater mass flow rate.

Figure 12 .
Figure 12.The variation in system NPV with the feed seawater mass flow rate.

Table 1 .
Studies related to multigeneration systems.

Table 2 .
Mass and energy balance relationships for the CO 2 cycle, ORC, LNG and water loops, and boiler. .

Table 3 .
Mass, energy, and concentration balance relationships for the RO system.32 (h 32 − h 33 ) x 32 = x 33 m

Table 4 .
Mass, energy, and concentration balance relationships for the NaClO plant.

Table 5 .
The exergy destruction rates and exergy efficiencies for system components.

Table 6 .
The specific costs of fuel and products.

Table 7 .
Values of investment and installation costs for different components.

Table 8 .
The overall heat transfer coefficient values for components.

Table 9 .
Input information of the simulation code.

Table 10 .
[32]comparison of key parameters in the present work and in[32].

Table 11 .
Validation results of the RO system.

Table 12 .
The thermodynamic properties at every point of the system.

Table 13 .
The values of system products.

Table 14 .
The key economic parameters for the system.