Performance Evaluation of a Solar-Powered Regenerative Organic Rankine Cycle in Different Climate Conditions

A model to evaluate the performance of a solar powered regenerative Organic Rankine Cycle (R-ORC) using five dry organic fluids: RC318, R227ea, R236ea, R236fa, and R218, is presented in this paper. The system is evaluated in two locations in the U.S.: Jackson, MS and Tucson, AZ. The weather data for each location is used to determine the heat available from the solar collector that could be used by the R-ORC to generate power. Results from the R-ORC performance are compared with a basic ORC using first and second law criteria as well as primary energy consumption (PEC) and carbon dioxide emission (CDE) savings for both locations. An economic analysis to determine the maximum capital cost for a desired payback period is presented in this paper. A parametric analysis is also performed to study the effect of the turbine efficiency as well as the open feed organic fluid heater intermediate pressure on the system performance. Results indicate that the R-ORC is able to generate more power than the basic ORC for some of the selected working fluids. For the R-ORC, R236ea is the working fluid that show the best performance among the evaluated fluids under the modeled conditions. On the other hand, the basic ORC with R236ea as the working fluid outperformed three of the fluids in the R-ORC. Also, the R-ORC evaluated in Tucson, AZ is able to generate more power, to provide more PEC and CDE savings, and had a higher available capital cost than the R-ORC evaluated in in Jackson, MS.


Introduction
Organic Rankine Cycles (ORCs) have become an area of increased interest and study because of their ability to use heat transfer from low to medium temperature sources to generate power.Various sources of heat can be used to power ORCs such as waste heat, geothermal energy, and solar energy [1][2][3][4][5].ORCs typically generate small amounts of electricity, so they are ideal for small scale applications.The working fluid greatly affects the performance of ORCs, therefore fluid selection and performance has been widely studied by several researchers [6][7][8][9][10][11]. Rayegan and Tao [12] evaluated 34 fluids for several solar powered ORC configurations based on thermal efficiency, exergetic efficiency, vapor expansion ratio and network output.Mago et al. [13] modeled a regenerative ORC using four dry fluids and performed a first and second law analysis.They found that of the evaluated fluids, R113 had the highest thermal efficiency and that the regenerative cycle increased the thermal efficiency while decreasing the irreversibility.
Regenerative ORCs (R-ORC) have been studied for various applications as well [14][15][16][17][18]. Roy and Misra [19] evaluated a R-ORC using R134a and R123 as the working fluid for waste heat recovery.They found that for a turbine inlet pressure of 2.5 MPa, R123 performed better than R134a.solar collector in line with an ORC and validated their proposed model with experimental results.They also performed a parametric study to optimize the system configuration and showed that the optimized configuration can produce net mechanical power with an efficiency up to 11%.In another study by Wang et al. [34] a flat plate solar R-ORC was analyzed using R123, R245fa, R134a, and isobutane as working fluids.In addition, they used a thermal storage tank to prevent the fluctuations in heat supplied to the ORC from the solar collector and used a regenerator to transfer heat from the working fluid as it left the turbine to the working fluid before it enters the evaporator.Pei et al. [35] modeled a solar powered R-ORC using a two stage compound parabolic concentrator system with thermal storage using phase change materials.They used a two stage compound parabolic concentrator and evaporator system to heat the ORC working fluid first to a saturated liquid in the first stage and then vaporize the working fluid in the second stage.They found that the regenerator increased the efficiency of the ORC and slightly increased the overall system efficiency.The authors found the regenerator temperature had a significant effect on the solar collector efficiency as well as the ORC, i.e., as the regenerator temperature increased, the ORC efficiency increased, but the solar collector efficiency decreased.
In a previous study, Spayde and Mago [36] studied a basic solar powered ORC using five dry organic fluids in Jackson (Mississippi) and Tucson (Arizona).In their study, hourly solar irradiation values were obtained from ASHRAE by latitude for the 21st day of each month.They performed a parametric study to investigate the effect of temperature, pressure, solar collector area, and turbine efficiency on the overall system performance.The effect of producing electricity by the ORC versus purchasing electricity from the grid on PEC and CDE was also investigated as well as the available capital cost.They reported that R236ea performed the best of the evaluated fluids under the modeled conditions and generated the most PEC and CDE savings and available capital cost.For all the evaluated fluids, the system showed the best performance in the middle of the day.They also reported that increasing the solar collector area increased the net energy produced and total exergy destroyed whereas increasing solar collector pressure and condensing temperature decreased both the net energy generated and total exergy destroyed.
The present study expands the previous study by the authors by investigating the performance of a solar powered R-ORC using several dry working fluids, RC318, R227ea, R236ea, R236fa, and R218, in two different locations in the U.S. The results are compared with a solar powered basic ORC to establish the benefits of the proposed R-ORC.A thermo-economic analysis is performed to determine the fluid that provides the best performance under the modeled conditions for the evaluated fluids.In this study a flat plate solar collector is modeled as the evaporator of the R-ORC system.Because the flat plate collector is in line with the ORC, the irradiation collected from the flat plate collector, which is therefore available to the ORC, will vary with the local weather.The solar irradiation data is calculated hourly from local typical meteorological year 3 (TMY3) weather data available from National Renewable Energy Laboratory.A parametric analysis is also performed to study the effect of the turbine efficiency as well as the open feed organic fluid heater intermediate pressure on the system performance.In addition to the parametric analysis, the proposed system is evaluated in two locations because weather has a direct effect on the performance of the proposed R-ORC.Jackson, MS and Tucson, AZ were chosen because they have approximately the same latitude but are in different climate zones.The performance of the proposed R-ORC is analyzed based on the PEC and CDE savings, and the available capital cost for both locations.

System Model
This section describes the model used to evaluate the performance of the solar powered R-ORC.Figure 1 shows a schematic of the evaluated R-ORC and the associated T-s diagram.As can be seen in Figure 1, the system is a direct vapor generation (DVG), since the organic working fluid is directly heated to the vapor state at the evaporator [37].a. Pump 1 (Process 1-2): The pump power can be determined as: where is the ideal power of Pump 1, is the isentropic efficiency of Pump 1, is the mass flow rate of the working fluid at State 5, is the fraction of the extracted fluid from the first stage of the turbine, and ℎ , ℎ , and ℎ are the enthalpy values of the working fluid at the pump inlet, the ideal enthalpy value at the pump exit, and the actual value at the pump exit, respectively.The exergy destruction rate for Pump 1 is: where is the pump exergy rate and and are the exergy rates for State 2 and State 1, respectively.Pump 1 (Process 1-2): The pump power can be determined as: .
W p1s is the ideal power of Pump 1, η p1 is the isentropic efficiency of Pump 1, .m 5 is the mass flow rate of the working fluid at State 5, X is the fraction of the extracted fluid from the first stage of the turbine, and h 1 , h 2s , and h 2 are the enthalpy values of the working fluid at the pump inlet, the ideal enthalpy value at the pump exit, and the actual value at the pump exit, respectively.
The exergy destruction rate for Pump 1 is: where .The exergy rate of Pump 1 is given by: .
The exergy rate change from State 1 to State 2 can be estimated as: .
where T o , s 1 , and s 2 are the temperature at the dead state (298 K), and entropy values of the working fluid for States 1 and 2, respectively.
b. Pump 2 (Process 3-4): The power required for Pump 2 is given by .
W p2s is the ideal power of Pump 2, η p2 is the isentropic efficiency of Pump 2, and h 3 and h 4 are the enthalpy values of the working fluid at State 3 and 4, which are the inlet and outlet of Pump 2, respectively.
The exergy destruction rate for Pump 2 can be determined as: where .
E p2 is the exergy rate of Pump 2 and  The exergy rate of Pump 2 is given by: .
The change in exergy rate between States 4 and 3 is: .
where s 4 and s 3 are the entropy values of the working fluid at States 4 and 3, respectively.
c. Solar Collector (Process 4-5): This is an isobaric process where heat is supplied to the organic working fluid before the turbine inlet after the fluid exits the second pump.The flat plate solar collector replaces the evaporator in a typical R-ORC.The solar collector heat transfer rate into the working fluid follows: .
where h 5 is the enthalpy of the working fluid at the exit of the solar collector.
The heat transfer rate from the solar collector can also be expressed as a function of irradiation: .
where η solar is the solar collector efficiency, I is irradiation, and A is the collector area.
The solar collector efficiency is determined using the relationship below: where y int is the y-intercept and m is slope.These two terms are provided by the manufacturer or a third party certification.In the proposed model, m = 4.910 W/m 2 • • C and y int = 0.706 [38].
The equation for solar collector efficiency is the Hottel-Whiller-Bliss equation [39] where y int and m correspond to: where F R is the collector heat removal factor, τ is the transmissivity of the glass cover plates, α is the absorptivity of the absorber plate, and U L is the losses due to conduction and radiation.
The irradiation values can be determined using the following equation: where I t is the total irradiation, I DN is the direct normal irradiation, θ is the incidence angle, I dH is the diffuse horizontal irradiation, Σ is the surface tilt angle, I tH is the total horizontal irradiation, and ρ is the ground reflectance.Direct normal irradiation, diffuse horizontal irradiation, and total horizontal irradiation can be found in TMY3 data from the National Renewable Energy Laboratory [40].The value for ground reflectance used in the paper is 0.2 which was taken from literature [41].The incidence angle and surface tilt angle are dependent on the solar collector configuration (placement).In this study the solar collectors were modeled to be 2 axis tracking solar collectors.This gives the maximum solar irradiation.A two axis tracking system maintains the incidence angle at zero.The two axis tracking system leads to the following surface tilt equation: where β is the solar altitude which is given by [41]: where L is latitude, δ is declination, and H is the hour angle.Declination can be found by using the following equation [42]: where n is the day of the year.
The exergy destruction rate of the solar collector is: where .

E .
Q s is the exergy rate due to the heat input to the solar collector and .
E 5 is the exergy rate of State 5.
The change in exergy rate across the solar collector is found by: .
where s 5 is the entropy of the working fluid at the collector exit.
The exergy rate of the solar collector can be estimated as: .
where T is the solar radiation temperature which is assumed to be 6000 K [43].
d. Turbine (Process 5-6, 7): The power of the two stage turbine is determined by: .
W ts is the power of the ideal turbine, η t is the turbine isentropic efficiency, h 6 and h 6s are the enthalpies of the working fluid for the exit of the first stage of the turbine for the real and ideal cases respectively, and h 7 and h 7s are the enthalpies of the working fluid of the outlet of the second stage of the turbine for the real and ideal cases, respectively.
The exergy destruction rate of the turbine is expressed as: where .
E t is the exergy rate of the turbine and .The change in exergy rates from the inlet to the outlets of the turbine is: .
where s 6 and s 7 are the entropy values at State 6 and 7, respectively.
The turbine exergy rate is: . e.
Open feed organic fluid heater (Process 6, 2-3): The extraction fraction is defined as: where .The open feed organic fluid heater exergy destruction rate is determined by: The exergy balance is: .
The exergy destruction rate of the condenser is: where .
E .

Q c
is the exergy rate of the condenser.
The exergy balance across the condenser is given by: .
The exergy rate of the condenser is given by the following equation: .
where T L is the low temperature heat sink which is assumed to be 303 K.

g. R-ORC Net Power:
The equation for the net power of the R-ORC is: .
h. R-ORC Efficiencies: The R-ORC thermal efficiency and R-ORC exergetic efficiency can be expressed as: where .

E
.

W net
and .

E .
Q in are the exergy of the products and the exergy input to the ORC. .

E .
W net can be estimated as: .
i. R-ORC Component Exergetic Efficiencies: The exergetic efficiency of each component of the R-ORC can be found by defining exergetic efficiency as the used exergy divided by the available exergy for each component.Equations ( 36) through (41) define each component's exergetic efficiency: Q in (38) Energies 2017, 10, 94 9 of 20 The contribution of each component to the total exergy destruction rate of the R-ORC can be estimated as: where Π i represents the exergy destruction for each component and Π total is the total exergy destruction rate of the R-ORC system: Primary Energy Consumption (PEC) Savings: Generating electricity on site using solar energy instead of purchasing electricity from the grid results in PEC savings and can be estimated by subtracting the PEC from the R-ORC (PEC R−ORC ) from the PEC from the conventional system (PEC Conv ) as follows: where PEC savings is the PEC savings, ECF PEC is the site-to-source conversion factor for electricity (grid purchase), which varies depending on location [23,44], and SCF PEC is the site-to-source conversion factor for electricity (solar), which was chosen to have a value of 1 [45,46].

k.
Carbon Dioxide Emission (CDE) Savings: Since the R-ORC utilizes solar energy, there are no CDE associated with the electricity generation.This results in CDE savings when compared to purchasing electricity from the grid.The CDE savings can be found by subtracting the CDE from the R-ORC (CDE R−ORC = 0) from the CDE from the conventional system (CDE Conv ) as follows: where CDE reduction is the CDE savings and ECF CDE is the carbon dioxide emissions factor for electricity which is location dependent [47].

l.
Cost Savings and Available Capital Cost: Available yearly capital cost (ACC) is a means to determine the maximum capital cost for a given payback period based on the cost of purchased electricity from the grid.It is based on the purchased electricity savings that result from using the solar powered R-ORC to generate electricity.The available capital cost can indicate the feasibility of implementing a solar powered ORC for a specific location.The following equation determines the cost savings and maximum capital cost available depending on the desired payback period: where Cost e is the yearly forecasted cost of electricity and PBP is the desired payback period.The cost of electricity is forecasted by plotting the average retail price of residential electricity from 2004 to 2015 [48] and performing a linear regression to estimate the future average retail price.

Results and Discussion
The model described in Section 2 was employed to simulate an hourly solar powered R-ORC in Jackson, MS and Tucson, AZ.The solar collector was modeled as a two axis tracking collector, so that the maximum irradiation values would be used.The hourly irradiation and ambient temperature data that was used to determine the solar collector efficiency and the heat transferred from the solar collector was found from TMY3 data for both locations selected for this study.Jackson and Tucson were chosen because they are at approximately the same latitude but are in different climate zones.Figure 2 depicts the yearly irradiation data used for both Jackson, MS and Tucson, AZ calculated from TMY3 weather data.Figure 2 illustrates that Tucson receives more irradiation than Jackson throughout the year.Tucson also had in general higher ambient temperatures which affected the solar collector efficiency.
Energies 2017, 10, 94 9 of 19 from 2004 to 2015 [48] and performing a linear regression to estimate the future average retail price.

Results and Discussion
The model described in Section 2 was employed to simulate an hourly solar powered R-ORC in Jackson, MS and Tucson, AZ.The solar collector was modeled as a two axis tracking collector, so that the maximum irradiation values would be used.The hourly irradiation and ambient temperature data that was used to determine the solar collector efficiency and the heat transferred from the solar collector was found from TMY3 data for both locations selected for this study.Jackson and Tucson were chosen because they are at approximately the same latitude but are in different climate zones.Figure 2 depicts the yearly irradiation data used for both Jackson, MS and Tucson, AZ calculated from TMY3 weather data.Figure 2 illustrates that Tucson receives more irradiation than Jackson throughout the year.Tucson also had in general higher ambient temperatures which affected the solar collector efficiency.Five dry organic working fluids were used in this paper: R236ea, R236fa, RC318, R227ea, and R218.As mentioned before, in this study, the state of the working fluid was assumed to be fixed as a saturated vapor as the fluid leaves the solar collector; therefore, the maximum available energy from the solar collector was harvested by changing the mass flow rate of the working fluid.Thus the changing irradiation values and ambient temperatures affected the amount of power generated by the R-ORC.In addition, for each of the selected fluids, the fluid enters both pumps were assumed to be saturated liquid.For all fluids, the condensing temperature was assumed to be 30 °C, the solar collector pressure was 2 MPa, and the pressure of the open feed organic fluid heater was set as the average pressure between the solar collector and the condenser.The pump and turbine isentropic efficiencies were assumed to be 80%.Since the mass flow rate of the organic working fluid changes depending on the amount of irradiation from the solar collector, the efficiencies of the pump and turbine will vary during the process, however, they were assumed constant to simplify the proposed model.To account for this, the effect of turbine efficiency on the net energy produced and the total exergy destroyed is investigated in a parametric study presented later in this paper.The pressure and temperature ranges for each fluid are given in Table 1 [49].Five dry organic working fluids were used in this paper: R236ea, R236fa, RC318, R227ea, and R218.As mentioned before, in this study, the state of the working fluid was assumed to be fixed as a saturated vapor as the fluid leaves the solar collector; therefore, the maximum available energy from the solar collector was harvested by changing the mass flow rate of the working fluid.Thus the changing irradiation values and ambient temperatures affected the amount of power generated by the R-ORC.In addition, for each of the selected fluids, the fluid enters both pumps were assumed to be saturated liquid.For all fluids, the condensing temperature was assumed to be 30 • C, the solar collector pressure was 2 MPa, and the pressure of the open feed organic fluid heater was set as the average pressure between the solar collector and the condenser.The pump and turbine isentropic efficiencies were assumed to be 80%.Since the mass flow rate of the organic working fluid changes depending on the amount of irradiation from the solar collector, the efficiencies of the pump and turbine will vary during the process, however, they were assumed constant to simplify the proposed model.To account for this, the effect of turbine efficiency on the net energy produced and the total exergy destroyed is investigated in a parametric study presented later in this paper.The pressure and temperature ranges for each fluid are given in Table 1 [49].

System Performance
The R-ORC was compared with a basic ORC using the same irradiation data and locations for each of the five evaluated fluids.Table 2 presents the thermal and exergetic efficiencies for each fluid for both the R-ORC and basic ORC.For all of the evaluated fluids, the R-ORC had higher thermal and exergetic efficiencies than the basic ORC when both systems operated using the same working fluid.In Table 3, the net energy generated and total exergy destroyed are presented for the five fluids in both the R-ORC and basic ORC.The R-ORC model was able to generate more net energy and had less total exergy destruction than the modeled basic ORC for each of the evaluated fluids.Among all the evaluated fluids, R236ea performed the best under the modeled conditions for both the R-ORC and basic ORC with respect to the thermal and exergetic efficiencies, the net energy produced, and the total exergy destroyed.To study the performance of the proposed solar R-ORC, the system was evaluated using four collectors with an area of 3.696 m 2 each.Under the modeled conditions, the R-ORC using R236ea as the working fluid produced 2852 kWh/year in Jackson, MS.The R-ORC using R236fa as the working fluid had slightly higher efficiencies, higher net energy produced, and lower total exergy destroyed than the basic ORC with R236ea as the working fluid.Whereas the remaining fluids performed worse than R236ea in the basic ORC regardless of the chosen cycle under the modeled conditions.RC318 had the largest percent increase, 13.27%, of net energy produced per year, and R236ea had the largest percent decrease, 1.84%, of total exergy destruction per year when comparing the R-ORC to the basic ORC.Table 4 shows the percent contribution to the total exergy destruction rate for each of the components in the R-ORC for each of the evaluated fluids.The solar collector had the highest percentage contribution to the exergy destruction rate for all of the evaluated fluids.R236ea had the lowest percentage contribution for the solar collector among the evaluated fluids with a value of 91.76%.The turbine had the next highest percent contribution to the total exergy destruction rate with a range of 4.22% to 1.95% for the evaluated fluids.R236ea had the highest percent contribution for the turbine at 4.22%.The open feed organic fluid heater was the third largest contributor to the total exergy destruction rate with a percentage range of 3.12% to 0.67% for the evaluated fluids.Again R236ea had the highest percent contribution for the open feed organic fluid heater among the evaluated fluids.
Since R236ea performed the best of the evaluated fluids under the modeled conditions, it is used as the organic working fluid for the remainder results presented in this paper.Figures 3 and 4 compare the basic ORC versus the R-ORC over the course of a year in Jackson, MS.The R-ORC produced more net energy and slightly less total exergy destruction per month than the basic ORC.The R236ea R-ORC produced 2852 kWh/year of net energy and destroyed 16,301 kWh/year of exergy while the R236ea basic ORC generated 2546 kWh/year of net energy and destroyed 16,331 kWh/year of exergy resulting in a 12.01% increase of net energy production and a 1.84% reduction of total exergy destruction.
Furthermore, the R-ORC was evaluated in Tucson, AZ and the results compared with those presented in Figures 3 and 4 for Jackson, MS.In Figure 5, the net energy produced and total exergy destroyed are compared for Jackson and Tucson.Since there was more solar irradiation in Tucson, the net energy produced and total exergy destruction was higher in Tucson when compared to Jackson.The increase in both net energy production and total exergy destruction was more pronounced in the summer months.The net energy generated per year for Tucson, AZ was 4162 kWh/year which was 46% increase from the net energy produced in Jackson, MS.Table 4 shows the percent contribution to the total exergy destruction rate for each of the components in the R-ORC for each of the evaluated fluids.The solar collector had the highest percentage contribution to the exergy destruction rate for all of the evaluated fluids.R236ea had the lowest percentage contribution for the solar collector among the evaluated fluids with a value of 91.76%.The turbine had the next highest percent contribution to the total exergy destruction rate with a range of 4.22% to 1.95% for the evaluated fluids.R236ea had the highest percent contribution for the turbine at 4.22%.The open feed organic fluid heater was the third largest contributor to the total exergy destruction rate with a percentage range of 3.12% to 0.67% for the evaluated fluids.Again R236ea had the highest percent contribution for the open feed organic fluid heater among the evaluated fluids.
Since R236ea performed the best of the evaluated fluids under the modeled conditions, it is used as the organic working fluid for the remainder results presented in this paper.Figures 3 and 4 compare the basic ORC versus the R-ORC over the course of a year in Jackson, MS.The R-ORC produced more net energy and slightly less total exergy destruction per month than the basic ORC.The R236ea R-ORC produced 2852 kWh/year of net energy and destroyed 16,301 kWh/year of exergy while the R236ea basic ORC generated 2546 kWh/year of net energy and destroyed 16,331 kWh/year of exergy resulting in a 12.01% increase of net energy production and a 1.84% reduction of total exergy destruction.
Furthermore, the R-ORC was evaluated in Tucson, AZ and the results compared with those presented in Figures 3 and 4 for Jackson, MS.In Figure 5, the net energy produced and total exergy destroyed are compared for Jackson and Tucson.Since there was more solar irradiation in Tucson, the net energy produced and total exergy destruction was higher in Tucson when compared to Jackson.The increase in both net energy production and total exergy destruction was more pronounced in the summer months.The net energy generated per year for Tucson, AZ was 4162 kWh/year which was 46% increase from the net energy produced in Jackson, MS.

Primary Energy Consumption and Carbon Dioxide Emissions
The effect of replacing electricity purchased from the grid with electricity generated using solar power on the PEC and CDE was also investigated.Producing electricity from solar power has the potential for PEC and CDE savings as compared to purchasing electricity form the utility grid.Table 5 presents the site-to-source conversion factor for electricity (grid purchase) for Jackson, MS and Tucson, AZ, the site-to-source conversion factor for electricity (solar), and the carbon dioxide emissions factor for electricity for both locations.

Primary Energy Consumption and Carbon Dioxide Emissions
The effect of replacing electricity purchased from the grid with electricity generated using solar power on the PEC and CDE was also investigated.Producing electricity from solar power has the potential for PEC and CDE savings as compared to purchasing electricity form the utility grid.Table 5 presents the site-to-source conversion factor for electricity (grid purchase) for Jackson, MS and Tucson, AZ, the site-to-source conversion factor for electricity (solar), and the carbon dioxide emissions factor for electricity for both locations.

Primary Energy Consumption and Carbon Dioxide Emissions
The effect of replacing electricity purchased from the grid with electricity generated using solar power on the PEC and CDE was also investigated.Producing electricity from solar power has the potential for PEC and CDE savings as compared to purchasing electricity form the utility grid.Table 5 presents the site-to-source conversion factor for electricity (grid purchase) for Jackson, MS and Tucson, AZ, the site-to-source conversion factor for electricity (solar), and the carbon dioxide emissions factor for electricity for both locations.
Table 5. Site-to-source conversion factor for electricity and carbon dioxide emissions factor for electricity.
ECF PEC (Jackson, MS) [44] 3.14 kWh/kWh ECF PEC (Tucson, AZ) [44] 3.06 kWh/kWh SCF PEC [45] 1 kWh/kWh ECF CDE (Jackson, MS) [47] 0.467 kg/kWh ECF CDE (Tucson, AZ) [47] 0.534 kg/kWh Figure 6 presents the PEC savings for both the R-ORC and basic ORC in Jackson, MS and Tucson, AZ.The PEC savings were higher in Tucson, AZ regardless of which cycle is used than Jackson, MS even though the ECF PEC was slightly higher in Jackson.The higher savings in Tucson can be contributed from the higher amount of net energy produced by the ORC due to the increased solar irradiation available and the decreased amount of electricity purchased from the grid.The R-ORC produced more PEC savings than the basic ORC for each location because there was an increase in electricity produced in the R-ORC compared to the basic ORC.The yearly PEC savings for the R-ORC in Jackson and Tucson were 6104 kWh/year and 8574 kWh/year, respectively.
Implementing a solar powered R-ORC has also the potential of reducing the CDE since the R-ORC emits practically zero CDE. Figure 7 presents the CDE savings for the basic ORC and R-ORC in Jackson and Tucson.As with the PEC savings, the possible CDE savings were higher in Tucson, AZ versus Jackson, MS.The R-ORC also produced more potential CDE savings than the basic ORC.The possible yearly CDE savings for the R-ORC in Jackson and Tucson were 1332 kg/year and 2224 kg/year, respectively.The reason that more CDE savings are obtained in Tucson, AZ is that the carbon dioxide emission factor for electricity are higher in this location and emissions which indicates that higher carbon content is present in the flue mix used to generate electricity in that area.
Table 5. Site-to-source conversion factor for electricity and carbon dioxide emissions factor for electricity.
ECFPEC (Jackson, MS) [44] 3.14 kWh/kWh ECFPEC (Tucson, AZ) [44] 3.06 kWh/kWh SCFPEC [45] 1 kWh/kWh ECFCDE (Jackson, MS) [47] 0.467 kg/kWh ECFCDE (Tucson, AZ) [47] 0.534 kg/kWh Figure 6 presents the PEC savings for both the R-ORC and basic ORC in Jackson, MS and Tucson, AZ.The PEC savings were higher in Tucson, AZ regardless of which cycle is used than Jackson, MS even though the ECFPEC was slightly higher in Jackson.The higher savings in Tucson can be contributed from the higher amount of net energy produced by the ORC due to the increased solar irradiation available and the decreased amount of electricity purchased from the grid.The R-ORC produced more PEC savings than the basic ORC for each location because there was an increase in electricity produced in the R-ORC compared to the basic ORC.The yearly PEC savings for the R-ORC in Jackson and Tucson were 6104 kWh/year and 8574 kWh/year, respectively.
Implementing a solar powered R-ORC has also the potential of reducing the CDE since the R-ORC emits practically zero CDE. Figure 7 presents the CDE savings for the basic ORC and R-ORC in Jackson and Tucson.As with the PEC savings, the possible CDE savings were higher in Tucson, AZ versus Jackson, MS.The R-ORC also produced more potential CDE savings than the basic ORC.The possible yearly CDE savings for the R-ORC in Jackson and Tucson were 1332 kg/year and 2224 kg/year, respectively.The reason that more CDE savings are obtained in Tucson, AZ is that the carbon dioxide emission factor for electricity are higher in this location and emissions which indicates that higher carbon content is present in the flue mix used to generate electricity in that area.

Available Capital Cost
Available Capital Cost (ACC) determines the maximum capital cost for a given payback period based on the cost saving potential from grid electricity price.The ACC for the proposed system is shown in Figure 8 for a payback period ranging from 1 to 10 years.The ACC values were determined based on the current and forecasted national average purchased electricity cost until 2025.The predicted national average for 2016 is 0.13 $/kWh and for 2025 is 0.16 $/kWh.Applying these predicted electricity costs resulted in ACCs of $4198 for Jackson, MS and $6127 for Tucson, AZ given a 10 year payback period.Tucson had a higher ACC because the modeled R-ORC produced more electricity in Tucson than in Jackson.This indicates that for the R-ORC to have a payback period of less than 10 years in Tucson the capital cost must be less than $6127 based on an increasing national average electricity purchase price.

Available Capital Cost
Available Capital Cost (ACC) determines the maximum capital cost for a given payback period based on the cost saving potential from grid electricity price.The ACC for the proposed system is shown in Figure 8 for a payback period ranging from 1 to 10 years.The ACC values were determined based on the current and forecasted national average purchased electricity cost until 2025.The predicted national average for 2016 is 0.13 $/kWh and for 2025 is 0.16 $/kWh.Applying these predicted electricity costs resulted in ACCs of $4198 for Jackson, MS and $6127 for Tucson, AZ given a 10 year payback period.Tucson had a higher ACC because the modeled R-ORC produced more electricity in Tucson than in Jackson.This indicates that for the R-ORC to have a payback period of less than 10 years in Tucson the capital cost must be less than $6127 based on an increasing national average electricity purchase price.

Available Capital Cost
Available Capital Cost (ACC) determines the maximum capital cost for a given payback period based on the cost saving potential from grid electricity price.The ACC for the proposed system is shown in Figure 8 for a payback period ranging from 1 to 10 years.The ACC values were determined based on the current and forecasted national average purchased electricity cost until 2025.The predicted national average for 2016 is 0.13 $/kWh and for 2025 is 0.16 $/kWh.Applying these predicted electricity costs resulted in ACCs of $4198 for Jackson, MS and $6127 for Tucson, AZ given a 10 year payback period.Tucson had a higher ACC because the modeled R-ORC produced more electricity in Tucson than in Jackson.This indicates that for the R-ORC to have a payback period of less than 10 years in Tucson the capital cost must be less than $6127 based on an increasing national average electricity purchase price.

Parametric Analysis
The effect of the intermediate pressure (extraction pressure) and the turbine efficiency on the net energy produced, the total exergy destroyed, and the system mass flow rate was investigated.Figure 9 shows the effect of the intermediate pressure on the net energy produced, total exergy destroyed, and the average mass flow rate for the months of January and July when the intermediate pressure varies from 0.3 MPa to 1.9 MPa.The net energy produced, total exergy destroyed, and average mass flow rate were higher in the month of July compared to those in January.The average mass flow rates for both January and July increased as the intermediate pressure increased.The maximum net energy produced and the minimum exergy destroyed occurred at the same intermediate pressure value of 0.75 MPa.Changing the turbine efficiency only affected the net energy produced and the total exergy destroyed.The effect of the turbine efficiency on the net energy produced and the total exergy destroyed for the months of January and July is displayed in Figure 10.The efficiency of the turbine was varied from 50% to 80%.As the turbine efficiency increased, the net energy increased and the total exergy destroyed decreased slightly.The percent increase of net energy produced did not change for the months of January and July.Likewise, the percent decrease of the total exergy destroyed did not change for January and July.As the turbine efficiency increased from 70% to 80%, the percent increase of net energy produced was 15%, and the percent decrease of exergy destruction was 2.2%.

Parametric Analysis
The effect of the intermediate pressure (extraction pressure) and the turbine efficiency on the net energy produced, the total exergy destroyed, and the system mass flow rate was investigated.Figure 9 shows the effect of the intermediate pressure on the net energy produced, total exergy destroyed, and the average mass flow rate for the months of January and July when the intermediate pressure varies from 0.3 MPa to 1.9 MPa.The net energy produced, total exergy destroyed, and average mass flow rate were higher in the month of July compared to those in January.The average mass flow rates for both January and July increased as the intermediate pressure increased.The maximum net energy produced and the minimum exergy destroyed occurred at the same intermediate pressure value of 0.75 MPa.Changing the turbine efficiency only affected the net energy produced and the total exergy destroyed.The effect of the turbine efficiency on the net energy produced and the total exergy destroyed for the months of January and July is displayed in Figure 10.The efficiency of the turbine was varied from 50% to 80%.As the turbine efficiency increased, the net energy increased and the total exergy destroyed decreased slightly.The percent increase of net energy produced did not change for the months of January and July.Likewise, the percent decrease of the total exergy destroyed did not change for January and July.As the turbine efficiency increased from 70% to 80%, the percent increase of net energy produced was 15%, and the percent decrease of exergy destruction was 2.2%.
The organic fluid at State 3 is pumped to State 4 and enters the flat plate solar collector where heat is transferred to the fluid.The organic working fluid at State 5 enters a two stage turbine, as a saturated vapor, where power is generated.In the first stage of the turbine, the fluid is expanded to State 6, where a fraction of the organic fluid is extracted into an open feed organic fluid heater operating at the extraction pressure (Intermediate pressure).The remainder of the working fluid expands through the second stage of the turbine to State 7 and then passes through a condenser (State 1).This portion of the fluid is pumped to the intermediate pressure and introduced into the open feed organic fluid heater at State 2. For this model, this fluid is then mixed and heated with the fraction of the fluid that exits the first stage of the turbine so that the working fluid leaving the open feed organic fluid heater is a saturated liquid at the intermedia pressure.After the two fluid streams combine, the organic working fluid enters the second pump, at State 3, to increase the pressure before entering the solar collector at State 4 and repeating the cycle.Assumptions for this model include: a constant irradiation rate for each hour, a steady state system, no pressure losses, the intermediate pressure is the average of the high and low pressures for each working fluid, and constant isentropic efficiencies for the pumps and the turbine.In the model presented in this paper dry organic fluids are used since the saturated vapor curve has a negative slope which has been proven to provide improved results verses fluids that have a positive slope for the saturated vapor curve (wet fluids)[13].Energies 2017, 10, 94 4 of 19 operating at the extraction pressure (Intermediate pressure).The remainder of the working fluid expands through the second stage of the turbine to State 7 and then passes through a condenser (State 1).This portion of the fluid is pumped to the intermediate pressure and introduced into the open feed organic fluid heater at State 2. For this model, this fluid is then mixed and heated with the fraction of the fluid that exits the first stage of the turbine so that the working fluid leaving the open feed organic fluid heater is a saturated liquid at the intermedia pressure.After the two fluid streams combine, the organic working fluid enters the second pump, at State 3, to increase the pressure before entering the solar collector at State 4 and repeating the cycle.Assumptions for this model include: a constant irradiation rate for each hour, a steady state system, no pressure losses, the intermediate pressure is the average of the high and low pressures for each working fluid, and constant isentropic efficiencies for the pumps and the turbine.In the model presented in this paper dry organic fluids are used since the saturated vapor curve has a negative slope which has been proven to provide improved results verses fluids that have a positive slope for the saturated vapor curve (wet fluids) [13].

Figure 1 .
Figure 1.(a) Schematic of modeled solar powered R-ORC; (b) T-s diagram for the solar powered R-ORC.

Figure 1 .
Figure 1.(a) Schematic of modeled solar powered R-ORC; (b) T-s diagram for the solar powered R-ORC.

E 3
are the exergy rates for State 4 and State 3, respectively.

E 7
are the exergy rates for State 6 and State 7, respectively.

m 6
is the mass flow rate of the working fluid at the exit of the first stage of the turbine.

Figure 2 .
Figure 2. Yearly irradiation data for Jackson, MS and Tucson, AZ calculated from TMY3 weather data [40].

Figure 2 .
Figure 2. Yearly irradiation data for Jackson, MS and Tucson, AZ calculated from TMY3 weather data [40].

Figure 3 .
Figure 3. Net energy produced per month in Jackson, MS for the basic ORC and R-ORC configurations using R236ea as the working fluid.

Figure 3 .
Figure 3. Net energy produced per month in Jackson, MS for the basic ORC and R-ORC configurations using R236ea as the working fluid.

Figure 4 .
Figure 4. Total exergy destroyed per month in Jackson, MS for the basic ORC and R-ORC using R236ea as the working fluid.

Figure 5 .
Figure 5. Net energy produced and total exergy destroyed in Jackson, MS and Tucson, AZ by the R-ORC using R236ea as the working fluid.

Figure 4 .
Figure 4. Total exergy destroyed per month in Jackson, MS for the basic ORC and R-ORC using R236ea as the working fluid.

Figure 4 .
Figure 4. Total exergy destroyed per month in Jackson, MS for the basic ORC and R-ORC using R236ea as the working fluid.

Figure 5 .
Figure 5. Net energy produced and total exergy destroyed in Jackson, MS and Tucson, AZ by the R-ORC using R236ea as the working fluid.

Figure 5 .
Figure 5. Net energy produced and total exergy destroyed in Jackson, MS and Tucson, AZ by the R-ORC using R236ea as the working fluid.

Figure 6 .
Figure 6.PEC savings from the basic ORC and R-ORC in Jackson, MS and Tucson, AZ using R236ea as the working fluid.

Figure 6 .
Figure 6.PEC savings from the basic ORC and R-ORC in Jackson, MS and Tucson, AZ using R236ea as the working fluid.

Figure 7 .
Figure 7. CDE savings from the basic ORC and R-ORC in Jackson, MS and Tucson, AZ using R236ea as the working fluid.

Figure 8 .
Figure 8. Available capital cost for Jackson, MS and Tucson, AZ for a payback period of up to ten years using a forecasted national average electricity price using R236ea as the working fluid.

Figure 7 .
Figure 7. CDE savings from the basic ORC and R-ORC in Jackson, MS and Tucson, AZ using R236ea as the working fluid.

Energies 2017, 10 , 94 14 of 19 Figure 7 .
Figure 7. CDE savings from the basic ORC and R-ORC in Jackson, MS and Tucson, AZ using R236ea as the working fluid.

Figure 8 .
Figure 8. Available capital cost for Jackson, MS and Tucson, AZ for a payback period of up to ten years using a forecasted national average electricity price using R236ea as the working fluid.

Figure 8 .
Figure 8. Available capital cost for Jackson, MS and Tucson, AZ for a payback period of up to ten years using a forecasted national average electricity price using R236ea as the working fluid.

Figure 9 .
Figure 9.The effect of intermediate pressure on the net energy produced, total exergy destroyed, and average mass flow rate for January and July in Jackson, MS.

Figure 9 .
Figure 9.The effect of intermediate pressure on the net energy produced, total exergy destroyed, and average mass flow rate for January and July in Jackson, MS.
p1 is the pump exergy rate and .

Table 1 .
[49]ical pressure[49]and pressure and temperature ranges for each of the evaluated fluids.

Table 2 .
Thermal and exergetic efficiencies for each of the evaluated fluids in both the basic ORC and R-ORC.

Table 3 .
The net energy produced and total exergy destroyed for each of the evaluated fluids for the basic ORC and R-ORC.

Table 4 .
Percentage of contribution for each device to the total exergy destruction rate for the R-ORC.

Table 4 .
Percentage of contribution for each device to the total exergy destruction rate for the R-ORC.