Investigating the Impact of Economic Uncertainty on Optimal Sizing of Grid-Independent Hybrid Renewable Energy Systems

: One of the many barriers to decarbonization and decentralization of the energy sector in developing countries is the economic uncertainty. As such, this study scrutinizes economics of three grid-independent hybrid renewable-based systems proposed to co-generate electricity and heat for a small-scale load. Accordingly, the under-study systems are simulated and optimized with the aid of HOMER Pro software. Here, a 20-year average value of discount and inﬂation rates is deemed a benchmark case. The techno-economic-environmental and reliability results suggest a standalone solar/wind/electrolyzer/hydrogen-based fuel cell integrated with a hydrogen-based boiler system is the best alternative. Moreover, to ascertain the impact of economic uncertainty on optimal unit sizing of the nominated model, the ﬂuctuations of the nominal discount rate and inﬂation, respectively, constitute within the range of 15–20% and 10–26%. The ﬁndings of economic uncertainty analysis imply that total net present cost (TNPC) ﬂuctuates around the benchmark value symmetrically between $478,704 and $814,905. Levelized energy cost varies from an amount 69% less than the benchmark value up to two-fold of that. Furthermore, photovoltaic (PV) optimal size starts from a value 23% less than the benchmark case and rises up to 55% more. The corresponding ﬁgures for wind turbine (WT) are, respectively, 21% and 29%. Eventually, several practical policies are introduced to cope with economic uncertainty.


Introduction
Energy acts as the foremost ingredient of thriving macro-and micro-economies in the world.In spite of this crucial role of energy in the communities, there have been controversial views, discussions, and debates surrounding the way that energy is yielded.This is due to the fact that conventional methods of electricity generation have long caused major damages to the environment, accordingly leading to the deteriorating issue of climate change [1,2].As such, viable alternative means and environmentally benign solutions of energy production, by which the detrimental impacts can be alleviated to a very low level, have been introduced and implemented over the past few decades [3,4].Among several methods available to generate clean, sustainable, and renewable energy, wind turbine (WT) and photovoltaic (PV) systems have been the most appealing ones and have been applied worldwide in suitable regions, mostly owing to the merits of cost-effectiveness, zero-or low-carbon emissions, and mature industry [5][6][7][8][9].Apart from the aforementioned benefits of utilizing renewables, all nations should consider that the global energy demand is projected to touch the value of 38,700 TWh by the year 2050 [10].This colossal amount of required energy implies and urges that the exploitation of renewables is of utmost prominence and priority.
The main shortcoming of renewables is the inherent issue of intermittency which can considerably decline the reliability of renewable-based systems.Fortunately, the concept of hybrid renewable energy (HRE) solutions has solved the problem to a great extent.Besides, HRE systems may impose less costs compared to single-source energy generation systems [11][12][13].Further, to raise the reliability and also lessen the power outage, they are usually coupled with energy storage units such as a diesel generator (DG), a battery, and a fuel cell (FC) [14].Few critical drawbacks, such as high operational and maintenance cost and harmful effects to the environment, are associated with DG power plants [15,16].Among two other above-mentioned backups, FC can constitute an environmentally sustainable solution if the required hydrogen (H 2 ) is provided by renewable means.
One step towards the blueprint for energy security, decarbonization, decentralization of energy sector, fulfilment of the Paris Agreement emissions targets, and sustainable energy production is to electrify remote island regions using autonomous renewable-based power generation plants.Hence, if this step is taken, then the growth of the nation and high living standards of those areas can be guaranteed [17].These places usually suffer from frequent power shortages and are very vulnerable to disruptions due to their far-off location and high dependency on fossil fuel power sites.Additionally, connecting them to the national grid will impose significant expenses linked to either capital or regular maintenance and review costs [18,19].Therefore, harvesting renewable resources of energy existing in the vicinity of these communities can noticeably aid corresponding countries in terms of reaching sustainable energy development.However, optimal design of HRE systems will remain a big question which requires a careful and thorough investigation.In this regard, the very first step, and the most important, is to conduct the techno-economic feasibility studies in order to avert any losses and to minimize the potential risks [20].
In this study, three configurations for a grid-independent HRE system comprising PV, WT, renewable H 2 -based FC, renewable H 2 -based boiler, natural gas (NG)-fueled boiler, and biomass generator (BMG) are scrutinized to ascertain optimal sizing and their technoeconomic performance for co-supplying a touristy village, Mesr in Isfahan province of Iran, with electricity and heat.To increase sustainability and efficiency of the proposed models, the excess electricity is sent to electrolyzer for generating H 2 .Then, during the times when solar and wind electricity cannot meet electric and thermal demands, this stored renewable H 2 is utilized to react with oxygen (O 2 ) for electricity generation.This transition occurs in a electrochemical process within FC-engaging electrodes and electrolytes, releasing water as a by-product [21].This effective usage of the surplus electricity implies the necessity of adding FC to the system, otherwise a great deal of electricity would be dumped without ever being exploited [22].Additionally, in each of other two proposed systems, it is supposed that boiler consumes renewable H 2 or NG for meeting thermal energy demands.Furthermore, to comprehensively examine the available resources in the area, a HRE system, which includes a BMG, is investigated.Finally, among these three systems, the most suitable one, with respect to technical, economic, environmental, and reliability aspects, is selected to be further analyzed under economic uncertainty.
Due to the volatile nature of economy in developing countries, unpredictable high fluctuations are likely to happen in the discount rate and inflation.For this reason, the study looks into the effect of economic uncertainty on optimization and economics of the chosen system in such countries.To this end, the fluctuations in the nominal discount rate and inflation are considered based on a 20-year average from mid-2000 to mid-2020, as represented by Ref. [23], which are in the range of 15-20% and 10-26%, respectively.Due to the volatile nature of economy in developing countries, un fluctuations are likely to happen in the discount rate and inflation.Fo study looks into the effect of economic uncertainty on optimization and chosen system in such countries.To this end, the fluctuations in the n rate and inflation are considered based on a 20-year average from midas represented by Ref. [23], which are in the range of 15-20% and 10-2 Figure 1 illustrates the waves of the nominal discount rate and inflatio the above-mentioned period.The main contribution of this study is that the proposed approac sizing of an HRE system to be applied in remote rural areas under econ which is inherent mostly in under-developed countries, can be deploy lar nations with islands and far-off villages.Additionally, a rigorous su has indicated that the the impact of economic uncertainty on optimal u systems has not been explored, which means this work could be a mar to the field.Thus, interested researchers can undertake the approach o ure 2) as a guideline to scrutinize the potential of electricity and heat the nominated regions, as well as to optimally size the units when eco The benefits of such a scheme include (I) making the remote places sel ergy generation, (II) decentralizing electricity coverage, and finally (III carbonization of energy sector set to meet the Paris Agreement targets.The main contribution of this study is that the proposed approach to optimal unit sizing of an HRE system to be applied in remote rural areas under economic uncertainty, which is inherent mostly in under-developed countries, can be deployed by other similar nations with islands and far-off villages.Additionally, a rigorous survey of literature has indicated that the the impact of economic uncertainty on optimal unit sizing of HRE systems has not been explored, which means this work could be a marked contribution to the field.Thus, interested researchers can undertake the approach of this study (Figure 2) as a guideline to scrutinize the potential of electricity and heat co-generation in the nominated regions, as well as to optimally size the units when economy is volatile.The benefits of such a scheme include (I) making the remote places self-sufficient in energy generation, (II) decentralizing electricity coverage, and finally (III) accelerating decarbonization of energy sector set to meet the Paris Agreement targets.

Figure 2.
The approach of the study.

Literature Review
In terms of incorporating economic uncertainty into the process of configuration of an off-grid HRE system, a thorough survey of literature has demonstrated that economic uncertainty has not yet been investigated for optimal sizing of units for electricity and heat co-generation.Nonetheless, some studies have examined the influence of uncertainty associated with other aspects of renewable-related projects; the most recent ones are as follows.The impact of economic policy uncertainties on energy-growth-emissions nexus was assessed for countries regarding the application of renewable and non-renewable energy resources [24].Furthermore, the difference in the impact of economic policy uncertainty on the investment of traditional and renewable energy enterprises was studied by [25].Uncertainties in wind velocity, solar radiation, water flow speed, load demand, and electricity price were investigated for planning a stand-alone renewable energy-based micro-grid [26].Bouchekara et al. [27] utilized a multi-objective approach to optimally configure an HRE system consisting of PV, WT, and DG for a micro-grid under load uncertainty.Blazquez et al. [28] scrutinized how policy instruments would work under different market conditions when considering the impact of price volatility and uncertainty of investments.Ebrahimi et al. [29] explored the uncertainty of renewable energy resources on optimal design of a multiple form energy in the presence of electric vehicle charging stations and the impact of this.

Literature Review
In terms of incorporating economic uncertainty into the process of configuration of an off-grid HRE system, a thorough survey of literature has demonstrated that economic uncertainty has not yet been investigated for optimal sizing of units for electricity and heat co-generation.Nonetheless, some studies have examined the influence of uncertainty associated with other aspects of renewable-related projects; the most recent ones are as follows.The impact of economic policy uncertainties on energy-growth-emissions nexus was assessed for countries regarding the application of renewable and non-renewable energy resources [24].Furthermore, the difference in the impact of economic policy uncertainty on the investment of traditional and renewable energy enterprises was studied by [25].Uncertainties in wind velocity, solar radiation, water flow speed, load demand, and electricity price were investigated for planning a stand-alone renewable energy-based micro-grid [26].Bouchekara et al. [27] utilized a multi-objective approach to optimally configure an HRE system consisting of PV, WT, and DG for a micro-grid under load uncertainty.Blazquez et al. [28] scrutinized how policy instruments would work under different market conditions when considering the impact of price volatility and uncertainty of investments.Ebrahimi et al. [29] explored the uncertainty of renewable energy resources on optimal design of a multiple form energy in the presence of electric vehicle charging stations and the impact of this.
Similarly, Li et al. [30] performed an economic analysis on micro-grids considering renewable energy uncertainty.In another research, a model was proposed for distributed generation expansion planning problem from the viewpoint of private investors considering uncertainties of renewable generations [31].Mirzamohammadi et al. [32] also investigated the role of renewable resources uncertainty for energy supply planning for greenhouse farming.Further, the role of uncertainty in wind and solar power was evaluated for power grid peak shaving [33].Optimal operation of a combined cooling, heating, and power micro-grid system was studied under the impact of wind power output uncertainty using an improved two-stage robust optimization model [34].Likewise, the potential risk and benefit of complementary operation of a large PV/WT/hydropower system were assessed under the consideration of forecast uncertainty in wind and solar power [35].Ezbakhe and Pérez-Foguet [36] introduced a modified multi-criteria decisionmaking technique to incorporate uncertainty of decision environments in the prioritization of alternatives for renewable energy planning.Fitiwi et al. [37] conducted an assessment on renewable energy utilization for the Irish power system under community acceptance policies considering the uncertainty of storage costs.Scott et al. [38] carried out a study on the role of long-term uncertainty in electricity market modelling for renewable energy support policy evaluation.There are other studies that consider uncertainty in different aspects of renewable-related projects, such as load demand uncertainty [39,40]; renewable energy resources uncertainty [41,42]; uncertainties of the inflow of the hydropower plant and the PV power output [43]; wind energy uncertainty [44]; uncertainty in the costs of renewable energy [45]; uncertainties of wind speed, solar radiation intensity, and power load [46]; and uncertainties of renewable energy production and load demand [47].
Since Hybrid Optimization Model for Electric Renewables (HOMER) software has proved applicability and reliability in simulation, analysis, and optimization of renewable power generation systems leading to precise results, it has been utilized in the past few years to explore HRE systems from the techno-economic-environmental perspective.Olatomiwa et al. [48] studied the viability of different configurations of PV, WT, battery, and DG units to be used in a villages in Nigeria using HOMER software.Considering total net present cost (TNPC), levelized cost of electricity (LCOE), and renewable fraction, the combination of PV/DG/battery was identified as the best architecture.In another research conducted in Nigeria, the techno-economic feasibility of a solar-based plant for powering a mobile base transceiver station was analyzed using HOMER.The results suggested a 10-kW PV, a 5.5-kW DG, and 64 units of Trojan battery.It also ended up saving 16.4 tons of CO 2 per year in comparison with an autonomous DG plant [49].Brenna et al. [50] utilized HOMER to assess the different integrations of PV, WT, hydro, DG, battery, and H 2 for electrifying a rural area in Ethiopia.Isa et al. [51] proposed an on-grid PV/FC/ battery system for co-supplying electricity and heat to a hospital in Malaysia with a TNPC of $106,551, an LCOE of 0.091 $/kWh, and an emission reduction of 25,873 kg/yr.Singh and Baredar [52] studied the electricity generation potential via an off-grid system consisting of PV, FC, BMG, and battery.Singh et al. [53] employed HOMER to verify the results of optimal sizing of an HRE system which was obtained by swarm-based artificial bee colony and particle swarm optimization.The comparison proved that these three methods had close results as to the sizing of the components.Das et al. [54] explored the techno-economic feasibility of an HRE system using HOMER.The software was employed to compare PV/FC/battery and PV/battery systems with a benchmark of the current DG plant [55].Khemariya et al. [56] used HOMER to evaluate a PV/FC/battery/electrolyzer system to electrify a village in India.Similarly, HOMER optimized a PV/BMG/DG/battery system considering different peak loads, energy demands, and grid availability [57].Shahzad et al. [58] designed an integrated PV/BMG/battery system for the irrigation and residential applications in Pakistan.Furthermore, HOMER assessed the techno-economic feasibility of a grid-tied PV/WT/BMG system for electrifying a village in Pakistan [59].Duman and Guler [60] explored the utilization of an autonomous PV/WT/FC plant to supply electricity for vacation homes in Turkey.HOMER revealed that the displacement of FC with batteries would turn the system more cost-effective.Moreover, stand-alone HRE systems were scrutinized using HOMER for a nursing home in Turkey [61].Apart from those above-mentioned studies, these papers  also applied HOMER for different purposes, which are helpful in learning about the software and its practicality.

Geographical Specifications
Mesr village, also locally called Mesr desert, is located in the province of Isfahan and is surrounded by the Great Salt Desert.This remote area is situated at coordinates 34 • 04 N and 54 • 48 E, with an altitude of rough 800 m above the sea level [89].According to the latest population and housing census conducted in 2016, the village has 41 households with 183 population [90].By virtue of Mesr's location, many national and international tourists visit this place every year.For this, in addition to the number of personal houses, there are a dozen traditional hotels to accommodate visitors, a healthcare center, a mosque, and a small school.After collecting the information about electricity demand from Isfahan Province Electricity Distribution Company, the daily average electricity load would be approximately 366 kWh.The reason why this village consumes this amount of energy is that the number of tourists visiting this touristy rustic area is high in several months of the year.Analyzing the figures and the common electricity-consuming devices, almost 163 kWh/d of electricity would be used for heating equipment, such as space heating, water heating, and cooking.The rest, which is 203 kWh/d, would be needed by other equipment and appliances, such as air conditioning, fan, bulbs, ironing, television, washing machine, dish washer, refrigerator, etc.Thus, two types of load, electric load and thermal load, are assumed in the proposed models.

Materials and Methods
HOMER Pro software was developed by the National Renewable Energy Laboratory (NREL) in the USA and has the ability to evaluate HRE systems in both modes of off-grid and on-grid [91].This software has been extensively utilized because of its suitable accuracy and high speed for optimization [92] and is comparable to meta-heuristic optimization algorithms [53].Hence, HOMER is applied to serve the purpose of this study.

Economic Analysis
The key economic criterion via which HOMER proposes optimal sizing of HRE systems and prioritizes different configurations is TNPC.The formula to compute TNPC reads as Equation (1) [93].
where C total,an is designated to the total annualized cost of all equipment to be purchased and installed in the system.CRF, i, and n, respectively, represent the capital recovery factor, the real yearly discount rate, and the project lifetime.Equation ( 2) can evaluate the term expressed as the denominator of the above-mentioned relation [94].HOMER also uses Equation ( 3) in order to project the real yearly discount rate in which i and f are, respectively, the nominal discount rate and yearly inflation [95].
The other prominent economic factor is the average price of generating one kWh of electricity via the designed system and is called LCOE, which can be obtained using Equation ( 4) [51].
Here, E Prim,AC and E Prim,DC signify the total amount of alternating current (AC) primary load and that of direct current (DC) primary load met by the system within a year.Additionally, E grid,Ex is the total amount of electricity which is exported to the local/national grid.
It is worth mentioning that another economic factor relatively important to the investors is the scrap value of the apparatus at the end of the project lifetime.This salvage value is directly contingent upon the remaining life of the nominated component holding a linear relation with a corresponding replacement cost.The relations are expressed by Equation ( 5) [51].
where C rep refers to the replacement cost of the component when its lifetime, n com , has finished.
The other economic output variable represented by HOMER is operating cost, equating to the annualized value of all costs and revenues except for the initial capital costs.Equation ( 6) is the relation to project operating cost [96].
where C total,ini is designated to the total initial capital costs, which equates to the summation of costs of all equipment in HRE systems.

PV Modeling
Owing to the geographical location of Iran, the country sees approximately 300 sunny days within a year, due to which the total annual solar radiation can reach 2200 kWh [83].This high potential of solar energy is ready and available to be exploited by PV systems in most parts of the country.To model PV arrays into renewable power generation systems, the output of PV is of high importance and needs to be calculated.HOMER software utilizes Equation (7) to project power output generated by PV [97].
where Y PV and f PV , respectively, denote the nominal capacity of the PV generator and the derating or reduction factor available on the technical sheet of the product.The latter is applied by HOMER to simulate the real-world operating conditions and then to acquire the most reliable and precise results.G T signifies the amount of solar radiation reaching the surface of the PV panels, and G T,STC represents the amount under the standard test environment.α P is designated to show the temperature coefficient of power.T c and T c,STC are, respectively, the temperature of the PV cell and that at standard test environment.In order to predict the cell temperature, Equations ( 8) and ( 9) are used [97].
where τ and α are, respectively, transmittance of PV cover and the solar absorption of the PV panels.U L means the coefficient via which heat is transferred to the surrounding environment, and T α refers to the amount of the ambient temperature.Furthermore, η c is the electrical conversion efficiency of the PV panels.HOMER uses Equation (10) in order to measure the cell temperature [85].It should be noted that, due to the difficulty of calculating the value of τα U L , Equation ( 11) is utilized to substitute the aforementioned term with the measurable one.Moreover, the software presumes that the value of τα in Equation ( 9) is 0.9 [96].
where T c,NOCT is the temperature usually reported by the manufacturers and called the nominal operating cell temperature.This value is obtained under the conditions of no-load operating, η c = 0, an ambient temperature of 20 °C, T α,NOCT = 20 °C, and an incident radiation of 800 W/m 2 , G T,NOCT = 800 W/m 2 .

WT Modeling
To capture the kinetic energy of air in motion, WT units are installed and then turn mechanical energy available in wind into electricity via a rotor with two or more blades coupled to an electric generator [70].Wind speed varies as the height changes [98]; hence, it is crucial to calculate wind velocity at the hub height.For this, HOMER applies the logarithmic law, as shown in Equation ( 12), to guess wind speed at the desired elevation [67].
Here, U hub and U anem are wind speed at the hub height, h hub , and wind speed collected at the height of anemometer, h anem , respectively.Furthermore, h 0 refers to the surface roughness length.Finally, the power produced by the WT generator is estimated using Equation ( 13) [97].
where ρ and ρ 0 , respectively, denote the actual density of air and the density of air at standard test conditions with certain pressure and temperature.P WT,STC is the WT output under standard test conditions which is to be computed by HOMER using the WT's power curve.Two other important values reported on technical sheet of the turbines are cut-off and cut-out wind speeds; no power will be generated if wind speed at the hub height is less than the former or more than the latter [99].

Electrolyzer Modeling
H 2 is the most favorable, futuristic energy carrier as it releases virtually zero emissions when being applied.There are no pure H 2 molecules available in the surrounding environment.However, it can be obtained via several chemical processes.To this end, water electrolysis is the promising method in which an electrolyzer utilizes electricity to make the decompaction of water into H 2 and O 2 molecules happen [99].The essential factor connected with modeling and sizing of an electrolyzer is the H 2 production rate which can be projected via Equation ( 14) [79].
Here, n c denotes the number of cells in series, F is the faraday constant, and I Ele is designated to represent the electrolyzer current.HOMER utilizes Equation (15) to evaluate the required input energy for electrolyzer [100].
Processes 2021, 9, 1468 9 of 25 where B Ele and A Ele are known as curve consumption coefficients.Q N and Q represent the nominal mass flowrate of H 2 and the mass flowrate of H 2 , respectively.

FC Modeling
Renewable-based power generation systems usually suffer from a major weakness derived from the intermittent nature of renewable resources during hours when the sunshine is not available or wind in not blowing.To solve this interrupting issue, one or more types of energy storage or backup systems should be added to the plant.Among which, FC has grasped the attention of system designers especially for rural electrification purposes, as it can be fed by H 2 generated via the aid of the excess electricity [16,88].In this study, it is postulated that the type of FC is proton exchange membranes, owing to its benefits such as low-cost maintenance, high energy conversion efficiency, and low operating temperature [101].To evaluate the output voltage of a FC, Equation ( 16) can be used [102].
where the four variables of E, V act , V ohm , V conc refer to open circuit voltage, activation FC overvoltage, ohmic FC overvoltage, and concentration FC overvoltage, respectively.It is worth mentioning that, in HOMER software, to add a FC to the proposed system, it should be selected from the category of generators.HOMER has provided the option of modifying the characteristics of a generator to simulate and imitate the specifications of a FC [103].Equation ( 17) is deployed to predict efficiency of a FC [85].
Here, V c , n c,FC , and I FC represent the average voltage of a cell in FC, the total number of cells in FC, and the FC current, respectively.HHV hyd is the gross calorific value existing in H 2 which is considered usually between 120 and 142 MJ/kg [104].

BMG Modeling
Living in the countryside and rural areas usually associates with livestock activities, thus the waste of cattle and other animals is amply available in these particular regions.However, this invaluable resource of energy is mostly being discarded without being properly exploited.In the case study area, most villagers have cattle, including cows, goats, sheep, chickens, camels, mules, and horses, to provide organic products or recreational activities for tourists.Therefore, animal manure is assumed to be the biomass source of the generator simulated and modeled in one of the proposed systems.To turn the energy prevailing in the bio waste into electricity, a gasifier is required for transforming solid bio-residue into a gaseous fuel.Equation ( 18) is given to project the yearly output electricity of a BMG [53].
where P bio and CUF signify the rating of BMG and the capacity utilization factor, respectively.Equation (19) evaluates the maximum rating of BMG which is deployed in the proposed system [105].
Here, T bio represents the total amount of biomass available to be used in the gasifier.C bio and η bio , respectively, denote the calorific value of the type of biomass to be applied and the efficiency of converting biomass to electricity.O bio refers to the hours of operating BMG.

Converter Modeling
One of the indispensable components of an HRE system is the converter in which an inverter and a rectifier are essential parts in order to make the transition from AC to DC, and vice versa, occur [68].The selection of inverter depends on AC loads and its efficiency can be calculated using Equation (20) [106].
where P inv,out and P inv,in are the output power and the input power of inverter, respectively.

Thermal Load Controller (TLC) Modeling
The existence of TLC is of high importance for the heat production system and so it is added to the proposed system in order to transfer the surplus electricity generated via the system to the thermal load.The main benefit brought to the designed system by TLC is improvement of renewable fraction, resulting in fewer emissions [97].To simulate the role of TLC in the system made by HOMER, it is just required to determine the size and the cost of TLC.

H 2 Tank Modeling
To store the produced H 2 in excess of the amount required by FC or boiler, a storage should be incorporated in the proposed system.The most vital feature of this component is its ability to endure holding the high-pressure H 2 [92].Equation ( 21) can calculate the pressure of H 2 stored in the tank [107].
Here, C H2 , θ H2 , and T represent the compressibility rate of H 2 , the gas constant, with a H 2 value at 4124.18 Nm/kg•K, and temperature, respectively.V H2 is also the specific volume.Similar to TLC, HOMER software just allows the users to set the volume of H 2 tank and its price.

Boiler Modeling
In order to add a boiler into the models, the type and price of the fuel, emission factors, and efficiency should be set for the software to simulate the unit.In this regard, efficiency means the fraction of the fuel's energy which is being turned into heat [96].

Technical Characteristics, Cost of Equipment and Assumptions
To simulate the electricity and heat co-generation system and optimize the size of each component, HOMER needs some technical features and the price of involved equipment.For this, Table 1 lists the technical specifications and costs of all components employed in the system.A nominal discount rate with a mean value of 17.5% and inflation with that of 18% are deemed the benchmark case for further judgment.Then, to assess the impact of economic uncertainty on optimal sizing, the nominal discount rate and inflation are postulated to fluctuate in the range of 15-20% and 10-26%, respectively, based on the figures of a period of 20 years [23].This consideration seems to be essential in order to obtain the most prospective results bounded by the worst-case and best-case scenarios.As mentioned, for the case when the boiler consumes NG, it is assumed that the price of fuel would equal $0.3/m 3 [92].Finally, it is presumed that the project lifespan would be 25 years.
As discussed, Mesr village is a touristy area attracting many visitors throughout the year.Therefore, electrical demand is relatively high compared to its number of residents.According to the load profile reported by Isfahan Province Electricity Distribution Company, electricity consumption for purposes other than heating possesses an average of 203.35 kWh/d with a peak of 27.44 kW.On the other side, daily average electricity requirement for thermal usage would be 163.23 kWh with a peak of 42.91 kW.It should be clarified that the figures associated with the electrical load revealed that electricity served for heating applications, such as space heating, water heating, and cooking, would be almost 80% of that served for other usages, such as lighting, cooling, kitchen appliances, etc. Figures 3 and 4 depict these two profiles on an hourly, daily, and monthly basis.

H2 Tank tank (H2Tank)
Initial tank level: 0 600 ($/kg) 600 ($/kg) 10 ($/yr) [109] A nominal discount rate with a mean value of 17.5% and inflation with that of 18% are deemed the benchmark case for further judgment.Then, to assess the impact of economic uncertainty on optimal sizing, the nominal discount rate and inflation are postulated to fluctuate in the range of 15-20% and 10-26%, respectively, based on the figures of a period of 20 years [23].This consideration seems to be essential in order to obtain the most prospective results bounded by the worst-case and best-case scenarios.As mentioned, for the case when the boiler consumes NG, it is assumed that the price of fuel would equal $0.3/m 3 [92].Finally, it is presumed that the project lifespan would be 25 years.
As discussed, Mesr village is a touristy area attracting many visitors throughout the year.Therefore, electrical demand is relatively high compared to its number of residents.According to the load profile reported by Isfahan Province Electricity Distribution Company, electricity consumption for purposes other than heating possesses an average of 203.35 kWh/d with a peak of 27.44 kW.On the other side, daily average electricity requirement for thermal usage would be 163.23 kWh with a peak of 42.91 kW.It should be clarified that the figures associated with the electrical load revealed that electricity served for heating applications, such as space heating, water heating, and cooking, would be almost 80% of that served for other usages, such as lighting, cooling, kitchen appliances, etc. Figures 3 and 4 depict these two profiles on an hourly, daily, and monthly basis.In this study, the three following configurations are analyzed as illustrated in Figures 5-7

The Benchmark Case (𝑖 = 17.5% and 𝑓 = 18%)
In the first model simulated and analyzed by HOMER, it was assumed that FC and boiler backups consume renewable H2 yielded by electrolyzer.In the second mode boiler, it would use NG as fuel and just FC would consume renewable H2.In the third model, BMG would be utilized instead of FC and just boiler would run on renewable H2 LCOE and  of the first model in its benchmark case would equate to 0.33 $/kWh and $647,708, respectively, which are higher than those of the two other configurations  In the first model simulated and analyzed by HOMER, it was assumed that FC and boiler backups consume renewable H 2 yielded by electrolyzer.In the second model boiler, it would use NG as fuel and just FC would consume renewable H 2 .In the third model, BMG would be utilized instead of FC and just boiler would run on renewable H 2 .LCOE and TNPC of the first model in its benchmark case would equate to 0.33 $/kWh and $647,708, respectively, which are higher than those of the two other configurations.However, at the end of the project lifetime, the salvage value of components utilized in the first model would be greater than that of the other models.Table 2 provides the results of optimal sizing and economic assessment of the models.The numbers with * mean that HOMER does not optimize them.So, they should be determined by the user of the software.For this, the size of electrolyzer is set 20 kW in the second model in which just FC consumes H 2 , while in the other models is 30 and 40 kW.
In relation to the technical results of these systems under the benchmark case, Table 3 demonstrates total electricity and thermal energy generation along with other findings.With regard to the unmet electric load, the first and second proposed systems showed promising figures, which would breed high reliability.Even though the third system had less TNPC compared to the two other models, its unmet electric load revealed that it could not be relied upon.Comparing the first and second models in terms of their TNPC, the latter would be more cost-efficient.Irrespective of TNPC, the second model would roughly emit 11.5 tons/year of CO 2 which is its main disadvantage, as it is not being eco-friendly.Based on the findings represented in Table 3, the first model would perfectly meet all electric and thermal energy demands, and would also be the most environmentally-friendly one.Thus, in the following sub-section, this model is analyzed considering economic uncertainty.

Analysis of the First Model under Economic Uncertainty
As mentioned, predicting the future value of money in underdeveloped countries may be onerous due to the high fluctuations in the discount rate and inflation.Hence, to obtain the most likely results as to the economics of PV/WT/electrolyzer/H 2 -based FC/H 2 -based boiler system, the least and the highest prospective amounts of these rates were incorporated into the projections.To this end, the nominal discount rate and inflation were considered to vary in the range of 15-20% and 10-26%, respectively, based on the 20-year average from mid-2000 to mid-2020.
The results proved that, if the nominal discount rate, i, rose, and the inflation rate, f , declined, then LCOE would increase.For instance, when i = 15% and f = 26%, the amount of LCOE would equate to 0.102 $/kWh, which is less than one-third of the benchmark value = 0.33 $/kWh.On the other side, when i = 20% and f = 10%, LCOE would augment significantly over twofold of the benchmark value to 0.662 $/kWh.Figure 8 depicts the surface plot of the variations in TNPC, considering the fluctuations of these two rates.

Analysis of the First Model under Economic Uncertainty
As mentioned, predicting the future value of money in underdeveloped countries may be onerous due to the high fluctuations in the discount rate and inflation.Hence, to obtain the most likely results as to the economics of PV/WT/electrolyzer/H2-based FC/H2-based boiler system, the least and the highest prospective amounts of these rates were incorporated into the projections.To this end, the nominal discount rate and inflation were considered to vary in the range of 15-20% and 10-26%, respectively, based on the 20-year average from mid-2000 to mid-2020.
The results proved that, if the nominal discount rate, , rose, and the inflation rate, , declined, then LCOE would increase.For instance, when  = 15% and  = 26%, the amount of LCOE would equate to 0.102 $/kWh, which is less than one-third of the benchmark value = 0.33 $/kWh.On the other side, when  = 20% and  = 10%, LCOE would augment significantly over twofold of the benchmark value to 0.662 $/kWh.Figure 8 depicts the surface plot of the variations in TNPC, considering the fluctuations of these two rates.Figure 9 demonstrates the surface plot concerned with annual electric production as a consequence of the variations in the rates.It can be inferred that, for the case when  lowered and  grew, the electric production would increase and LCOE would dimin- Figure 9 demonstrates the surface plot concerned with annual electric production as a consequence of the variations in the rates.It can be inferred that, for the case when i lowered and f grew, the electric production would increase and LCOE would diminish.Furthermore, Figure 10 provides the status of annual thermal energy generation considering different values of i and f .Figure 9 demonstrates the surface plot concerned with annual electric production as a consequence of the variations in the rates.It can be inferred that, for the case when  lowered and  grew, the electric production would increase and LCOE would diminish.Furthermore, Figure 10 provides the status of annual thermal energy generation considering different values of  and .Finally, Figures 11-13, respectively, show electricity generation via FC, via PV, and via WT units, when  and  vary.Finally, Figures 11-13, respectively, show electricity generation via FC, via PV, and via WT units, when i and f vary.Finally, Figures 11-13, respectively, show electricity generation via FC, via PV, and via WT units, when  and  vary.In Figure 13, the central part of the chart shows the values related to the benchmark case.

Discussion and Suggestions for Implementation
One of the most vital criteria that would hinder the exploitation of renewable energies in underdeveloped nations is economic uncertainty.Owing to some unforeseeable factors impinging upon the value of currency in these countries, speculation on economics of projects seems to be difficult.For example, Iran has experienced many ups and downs in its currency value due to lifting and imposing international sanctions.However, this issue can be dealt with via considering the fluctuations in the discount rate and inflation.To this end, in this study, the mean values of these rates were first incorporated into the calculations to set a benchmark case.Then, the most possible range of changes for these two rates were taken into account.
As the results of economic uncertainty analysis indicated, LCOE of the nominated system, PV/WT/electrolyzer/H2-based FC/H2-based boiler, would not be cost-competitive to conventional methods of electricity generation under the majority of likely values of the discount rate and inflation.Here, the role of the government and its underwriting policies are extremely decisive in the implementation of the proposed system in the nominated remote area.To achieve energy security, to decarbonize and decentralize energy sector, to meet the Paris Agreement emission targets, and to sustain the environment, governments should underpin such schemes by the following.
(I) Furnishing the investors or private companies with zero percent or low-rate loans.(II) Introducing carbon tax to encourage the generation and use of renewable electricity, (III) Setting strict rules and regulations against carbon-intensive means of generating electricity.In Figure 13, the central part of the chart shows the values related to the benchmark case.

Discussion and Suggestions for Implementation
One of the most vital criteria that would hinder the exploitation of renewable energies in underdeveloped nations is economic uncertainty.Owing to some unforeseeable factors impinging upon the value of currency in these countries, speculation on economics of projects seems to be difficult.For example, Iran has experienced many ups and downs in its currency value due to lifting and imposing international sanctions.However, this issue can be dealt with via considering the fluctuations in the discount rate and inflation.To this end, in this study, the mean values of these rates were first incorporated into the calculations to set a benchmark case.Then, the most possible range of changes for these two rates were taken into account.
As the results of economic uncertainty analysis indicated, LCOE of the nominated system, PV/WT/electrolyzer/H 2 -based FC/H 2 -based boiler, would not be cost-competitive to conventional methods of electricity generation under the majority of likely values of the discount rate and inflation.Here, the role of the government and its underwriting policies are extremely decisive in the implementation of the proposed system in the nominated remote area.To achieve energy security, to decarbonize and decentralize energy sector, to meet the Paris Agreement emission targets, and to sustain the environment, governments should underpin such schemes by the following.
(I) Furnishing the investors or private companies with zero percent or low-rate loans.
(II) Introducing carbon tax to encourage the generation and use of renewable electricity.(III) Setting strict rules and regulations against carbon-intensive means of generating electricity.(IV) Developing the concept of green tourism to attract as many national and international visitors as possible.The corresponding revenues can cover a substantial proportion of the project's costs.(V) Subsidizing the price of renewable electricity for residents (can be achieved from the resource of funding allocated to operating and maintaining the transmission and distribution network as it would no longer be needed).(VI) Lifting tariffs on importing equipment such as PV, WT, electrolyzer, FC, etc.
It is less doubtful that these carbon neutral remote power generation systems will be eventually valued in the marketplace, owing to their differentiated and distinguished quality of energy, if serious constraints in the energy market are imposed on carbon production.

Conclusions
In the first part of analysis, three models for co-supplying a touristy village with electricity and heat were proposed and analyzed in terms of technical, economic, environmental, and reliability aspects under the conditions of the benchmark case, i = 17.5% and f = 18%, taken from the average value of a 20-year period.The following results were obtained:

•
The first model, the PV/WT/electrolyzer/H 2 -based FC/H 2 -based boiler, had the highest TNPC ($647,708), the lowest unmet electric load, and the highest reliability without any detrimental impact on the environment.

•
The second model, the PV/WT/electrolyzer/H 2 -based FC/NG-based boiler, possessed the second least TNPC ($548,906), and it could meet almost all electric demand.Whereas, utilizing it would end up releasing some 11.5 tons of CO 2 per year.This carbon footprint constitutes a challenging negative point for the second model which may strongly inhibit all the attempts to accomplish the Paris Agreement targets.

•
The techno-economic analysis of the third model, the PV/WT/BMG/electrolyzer/H 2based boiler, showed that it would not be reliable, as 20.5% of total electric load could not be met via this system.However, its TNPC, $488,878, was the least amongst the three analyzed configurations.
In the second part of the paper, the fluctuations in i and f were in accordance with the 20-year average from mid-2000 to mid-2020, respectively, between 15-20% and 10-26%.The results of technical and economic assessment of the first model under economic uncertainty and comparison with the related values of the benchmark case are as follows:

•
The amount of LCOE would vary from 0.102 $/kWh to 0.662 $/kWh, meaning LCOE could be between one-third of the benchmark value and two-fold that (LCOE for the benchmark case = 0.33 $/kWh).Additionally, TNPC would fluctuate between $478,704 and $814,905 from 26% less than the benchmark value up to 26% more than that (TNPC for the benchmark case = $647,708).

•
The optimal size of PV and the number of WT units would change from 25.9 to 52.5 kW and from 11 to 18, respectively.Comparing with the benchmark case (PV size = 33.8kW and number of WT units = 14), the PV size could vary from an amount of 23% less than the benchmark case up to 55% more than that, and corresponding figures for WT would be 21% and 29%, respectively.

•
The amount of renewable H 2 consumed by boiler and FC would be in the ranges of 1815-1962 kg and 559-665 kg, respectively.When comparing with the benchmark (H 2 consumption in boiler = 1922 kg and that in FC = 618), the former would fluctuate from an amount of 6% less than the benchmark value up to an amount of 2% more than that, and related numbers for FC would be 10% and 8%, respectively.

Figure 1
Figure 1 illustrates the waves of the nominal discount rate and inflation in Iran during the above-mentioned period.

Figure 1 .
Figure 1.The fluctuations in (a) the nominal discount rate and (b) inflation mid-2020.

Figure 1 .
Figure 1.The fluctuations in (a) the nominal discount rate and (b) inflation from mid-2000 to mid-2020.

Figure 2 .
Figure 2. The approach of the study.

Figure 3 .
Figure 3.The hourly, daily and monthly profile of electric load.

Figure 3 .
Figure 3.The hourly, daily and monthly profile of electric load.

Figure 4 .
Figure 4.The hourly, daily and monthly profile of thermal load.

Figure 4 .
Figure 4.The hourly, daily and monthly profile of thermal load.

Figure 4 .Figure 5 .
Figure 4.The hourly, daily and monthly profile of thermal load.

Figure 5 . 25 Figure 6 .
Figure 5. Schematic of the first model simulated by HOMER.Processes 2021, 9, x FOR PEER REVIEW 14 of 25

Figure 6 .
Figure 6.Schematic of the second model simulated by HOMER.

Figure 6 .
Figure 6.Schematic of the second model simulated by HOMER.

Figure 7 .
Figure 7. Schematic of the third model simulated by HOMER.

Figure 7 .
Figure 7. Schematic of the third model simulated by HOMER.

Figure 8 .
Figure 8. Changes of TNPC considering the fluctuations of the nominal discount rate and inflation (numbers inside the surface plot refer to LCOE).

Figure 8 .
Figure 8. Changes of TNPC considering the fluctuations of the nominal discount rate and inflation (numbers inside the surface plot refer to LCOE).

Figure 8 .
Figure 8. Changes of TNPC considering the fluctuations of the nominal discount rate and inflation (numbers inside the surface plot refer to LCOE).

Figure 9 . 25 Figure 9 .
Figure 9. Changes of annual electric production considering the fluctuations of the nominal discount rate and inflation (Numbers inside the surface plot refer to LCOE).

Figure 10 .
Figure 10.Changes of annual thermal energy generation considering the fluctuations of the nominal discount rate and inflation (numbers inside the surface plot refer to the amount of H2 consumed by boiler).

Figure 11 .
Figure 11.Changes of electricity generation via FC considering the fluctuations of the nominal discount rate and inflation (numbers inside the surface plot refer to the amount of H2 consumed by FC).

Figure 10 .
Figure 10.Changes of annual thermal energy generation considering the fluctuations of the nominal discount rate and inflation (numbers inside the surface plot refer to the amount of H 2 consumed by boiler).

25 Figure 9 .
Figure 9. Changes of annual electric production considering the fluctuations of the nominal discount rate and inflation (Numbers inside the surface plot refer to LCOE).

Figure 10 .
Figure 10.Changes of annual thermal energy generation considering the fluctuations of the nominal discount rate and inflation (numbers inside the surface plot refer to the amount of H2 consumed by boiler).

Figure 11 .
Figure 11.Changes of electricity generation via FC considering the fluctuations of the nominal discount rate and inflation (numbers inside the surface plot refer to the amount of H2 consumed by FC).

Figure 11 .
Figure 11.Changes of electricity generation via FC considering the fluctuations of the nominal discount rate and inflation (numbers inside the surface plot refer to the amount of H 2 consumed by FC).

Figure 11 .
Figure 11.Changes of electricity generation via FC considering the fluctuations of the nominal discount rate and inflation (numbers inside the surface plot refer to the amount of H2 consumed by FC).

Figure 12 .
Figure 12.Changes of electricity generation via PV considering the fluctuations of the nominal discount rate and inflation (numbers inside the surface plot refer to the optimal size of PV).

Figure 12 . 25 Figure 13 .
Figure 12.Changes of electricity generation via PV considering the fluctuations of the nominal discount rate and inflation (numbers inside the surface plot refer to the optimal size of PV).Processes 2021, 9, x FOR PEER REVIEW 18 of 25

Figure 13 .
Figure 13.Changes of electricity generation via WT units considering the fluctuations of the nominal discount rate and inflation (numbers inside the surface plot refer to the optimal number of the wind turbines).

Table 1 .
Technical specifications and costs of the components.

Table 2 .
Optimization results of simulating the models.

Table 3 .
Technical results of analyzing the models.