A Comprehensive Technoeconomic and Environmental Evaluation of a Hybrid Renewable Energy System for a Smart Farm in South Korea

. The farming sector like any other industry needs to be decarbonized. Hence, it is essential to meet the energy demands of the farms by adopting energy systems with a low-carbon footprint. Depending on the weather conditions, heating or cooling is needed. Within this context, this study presents a new hybrid renewable decentralized energy system that is designed to satisfy the requirements for heating, cooling, and electricity of a smart farm in South Korea. The under-investigation energy system comprises solar PV arrays, heat pumps, thermal energy storage tanks, and a wood pellet boiler. This study is the ﬁ rst to conduct an inclusive techno-enviroeconomic assessment of such a hybrid energy system by utilizing actual meteorological data on an hourly basis. This enables the model to be dynamic and facilitate accurate and reliable assessments. The modelling e ﬀ orts have been performed in Aspen Plus and MATLAB to investigate the thermodynamic behaviour of the system. The investigation shows that the proposed system has achieved a daily average temperature of around 23.9 ° C inside the farm throughout the year with a standard deviation of 2.16 ° C. For the economic assessment, the levelized cost of energy has been selected as the main economic indicator, and this has been estimated at $0.218/kWh. It is found that the PV panels and the biomass boiler dominate the capital expenditures, and the biomass feedstock is the major contributor to the operating expenditures. Further, the proposed energy system reduces CO 2 emissions, by up to 88.94%, when compared to conventional fossil-based energy systems. The outcomes of this study represent a holistic evaluation for such a low-carbon hybrid energy system when applied to greenhouses in Korea and in similar locations.


Introduction
The globally unprecedented pace of climate change has introduced a new dimension to the generation and consumption of energy, and this has been exacerbated by the existing energy crisis. Therefore, most countries have attempted to accelerate their rate of implementing renewable and sustainable energy sources, since this is a very important and vital approach for all the world to reduce the reliance on fossil fuels. One pathway that is being explored to achieve this transition is encouraging, namely, the development of hybrid energy systems, especially from the use of the available local renewable sources, such as biomass and solar [1]. In addition, most governments have implemented several interesting policies that could incentivise the deployment of clean energy solutions. Some of these policies are presented in the form of investment grants, feed-in tariffs, feed-in premiums, and green energy credits [2]. One of the most interesting initiatives was the introduction of the concept of "prosumers of energy," which has opened the door for the utilization of different clean and renewable energy sources [3,4]. This can be in the form of hybrid energy systems, which introduce several ingenious strategies to improve the energy efficiency, such as the cogeneration of heating and electricity.
Recently, the agricultural sector in the Republic of Korea has been paying more attention to the mitigation of the GHG emissions, and this is believed to be possible to be attained through the employment of smart farms as a replacement to the traditional farms [5]. The aim of the smart farms is to achieve the prosumer concept, in which the energy is produced from the renewable sources and consumed on site. In addition, the excess energy is supplied to the grid or stored. However, the smart farms require high investment costs to meet the high energy demands which are mainly required for the heating purposes [6]. About 53.4% of the total energy consumed by the agricultural sector is devoted to the heating and cooling systems, where this demand is mainly covered by fossil fuels [7]. Therefore, the need to develop an efficient and cost-effective energy system based on renewable energy sources is necessary to meet the carbon neutralization targets [8,9]. It is important to note that the agricultural sector in the Republic of Korea targets to achieve 4.8% of the total national targeted emission reduction of CO 2 [6]. This can be achieved by adopting different renewable energy sources in hybrid configurations.
It is important to note that the cogeneration of power and heat from the hybrid renewable energy systems can increase the efficiency and the reliability of such a system while maintaining low-carbon emissions. The photovoltaic power generation has evolved significantly over the recent years, and it has been proved to be suitable and reliable to produce energy over a long period of time [10]. Also, when it comes to the efficient heating equipment, the heat pumps are on top of the list for power-to-heat applications, with a high coefficient of performance (COP) that can reach over 6 as reported in [11]. The heat pumps can be distinguished by one of the two main types, ground-source heat pump (GSHP) and air-source heat pump (ASHP).
Many research studies have been conducted to evaluate the performance of the GSHP for greenhouse heating and cooling purposes. For instance, Chai et al. [12] investigated the performance of a GSHP for greenhouse heating as a replacement of the traditional coal-or gas-fired heating system in Beijing, China. The heat pump has proved its suitability to satisfy the heating load with an average COP of 3.87 at even lower daily heating cost than the gas-fired heater by an average of 10.15%. In addition, even though the GSHP was not cost-competitive against the coal-fired heater, the environmental impact should be taken into consideration. Similarly, Luo et al. [13] compared the thermoeconomic performance of a coal-fired heating system against a GSHP for greenhouse heating in China. The maximum thermal efficiency of the boiler was 68%, while the COP of the heat pump was 3.3 for cooling and 4.1 for heating, and this represents a high level of energy savings by the employment of the heat pump. Furthermore, the considered lifetime of the coal-fired boilers is short, 10 years vs. 20 years for the GSHP. This in turn favours the economic performance of the GSHP over its lifetime, where the average energy price was $0.023/kWh and $0.04/kWh for the cooling and heating, respectively, whereas for the coal-fired boiler, it was $0.052/kWh. The significant difference between the COPs of the heat pumps and the thermal efficiencies of the conventional systems also implies high CO 2 emission reduction (CO 2 ER). For instance, the heat pump has achieved an average reduction of 43.25% when compared to the coal-fired heating [12]. Interestingly, these high COPs favoured the economic and environmental performance, and this still can be observed for combined systems. In a study conducted in Japan by Zhou et al. [14], they showed that combining a GSHP with a fuel oil heater as a hybrid system for greenhouse heating can reduce the energy consumption by 22.8% and the CO 2 emissions by 35.5% when compared to the employment of only an oil heater.
The performance of the heat pump can be further improved by utilizing solar energy, which can be in the form of thermal energy or electrical energy. The solar thermal energy can be employed to heat the water or the air for the evaporator of the heat pump using the solar collectors [7,[15][16][17][18]. In order to evaluate this integration, Hassanien et al. [16] have investigated the feasibility of implementing a solar-assisted heat pump for the heating of a greenhouse in China. The integration was achieved by employing an evacuated tube solar collector with an efficiency of 49% to heat the water for the heat pump. Through this integration, the solar collector, by itself, has supplied more than 35% of the heat demand over one year. This represents the estimated energy savings over one year from this integrated system, which is mainly attributed to the contribution of the solar thermal energy that increases the COP of the heat pump [15].
Recently, in order to reduce the carbon footprint even further, the PV technology has been increasingly adopted for the power generation as well as for heating and cooling purposes through the integration with the heat pumps [10,19]. For heating purposes, Russo et al. [20] compared the environmental impact of a GSHP powered by solar PV panels and the grid to an LPG-based heat generation system in Italy. It was found that the GSHP-PV combination has achieved a 50% emission reduction of CO 2 relative to the LPG-based heating. This can be even higher, up to 80%, when compared to an electricity-based heating system based in the Netherlands [21]. Also, for the PV-based combined heating and power system (CHP), Yildirim and Bilir [22] have evaluated the performance of such a system in a greenhouse by employing solar PV panels along with a GSHP to provide power and heat to the farm with low emissions in Turkey. The PV system has satisfied between 33.2% and 67.2% of the total energy demand during the summer months, but the net annual coverage ratio of the energy load by the solar PV panels varied between 86.6% and 104.5%.
It can be observed from the literature that integrating the solar energy with the heat pump represents a promising approach to mitigate the environmental impact of a greenhouse CHP system. However, there are several limitations on the performance of the solar energy systems (i.e., the collectors and the PV panels), which are mainly associated to 2 International Journal of Energy Research the weather condition. To overcome these limitations, the solar PV panels can be backed up by electricity from the grid to run the heat pump when applicable. In addition, the heat pump is commonly integrated with a boiler as a backup heat source to the solar collector, or as an alternative to the heat pump itself. In this regard, biomass-fuelled boilers are considered as carbon neutral systems, especially when the biomass is grown at a rate equivalent to, or more than, its consumption. The other advantage of biomass-based systems is their ability to run on several locally available biomass fuels [23]. Finally, the utilization of actual meteorological data is of paramount importance for carrying out accurate technoeconomic and environmental assessments [12].
In particular, in the Republic of Korea, a hybrid system of a hydrothermal heat pump and a wood pellet boiler has been investigated for greenhouse heating [24]. The hybrid system has achieved higher thermal efficiency and lower running cost than the pellet boiler alone. Also, when compared to the conventional fossil fuel boiler, the carbon emission has been reduced by up to 71% [24]. However, when the system is integrated with a solar water heater, lower energy efficiency has been estimated, and therefore, it was not cost-effective. This is due to the limited availability of heat during the winter season in Korea. On the other hand, the PV-based hybrid systems can represent a more suitable approach, especially during cold periods. In addition, another advantage of the PV-ASHP energy system is the ability to run the heat pump for the farms in rural areas (off the grid), and in case of the grid-connected PV systems, the surplus electricity can be exported to the grid.
To the best of our knowledge, there is no study on solar PV-ASHP hybrid systems for greenhouse heating in Korea. Also, the previous studies are time-constrained, as they consider only a few operational days of the year, and this represents a shortcoming in the economic and environmental assessment. In view of these, the present study is aimed at investigating the performance of a hybrid renewable energy system comprising PV solar panels coupled with ASHPs and a wood pellet boiler to supply the energy demand of a smart farm in the Republic of Korea. In order to evaluate the feasibility of such a system, a comprehensive technoeconomic and environmental investigation has been conducted. The model utilizes the onsite operating data of the farm along with the hourly based meteorological data that is obtained from the same location. The outcomes of this study can and do provide meaningful insights into the viability of such a hybrid energy system and how it can contribute to the transition of the agricultural sector towards renewable decentralized energy solutions in Korea and in similar locations. Further, due to the modularity of the individual energy components, the proposed energy system has great scalability potential.

System Description
The smart farm that is investigated is located at the Jeollanamdo Agricultural Research and Extension Services, 1508, Senam-ro, Sanpo-myeon, Naju-si, Jeollanam-do, Republic of Korea. Figure 1 shows a schematic overview of the farm in addition to examples of the cultivated crops in 2021. The farm's total area is 1,378 m 2 and the height is 4.2 m.
The under-investigation energy system is an integration of solar PV-powered heat pumps and a backup wood pellet boiler in order to improve the reliability of the system. The system conditions the air in the farm all year round, and hence, it operates in winter and summer modes. Throughout the year, the grid-connected solar PV system is the primary source of power, and it is used to operate the heat pumps, the water circulation circuit, and the fan coil units (FCUs). The PV panels operate throughout the whole year, regardless of how much the heat pump operates, while any deficiency or excessiveness of the solar power is backed up and balanced by the grid. In case of the absence of solar energy, the wood pellet boiler is engaged to satisfy the heat load of the smart farm; this is only necessary for the wintertime. When no heating or cooling is required by the greenhouse farm, all the energy from the solar PV panels is exported to the grid. The actual layout of the piping configuration for the smart farm energy system is presented in Figure 2.
The layout exhibits the configuration of the energy system that is planned to be commissioned. The smart farm is designed to be sectioned to two planting areas (60%, i.e., 1 st smart farm, and 40%, i.e., 2 nd smart farm) to host different agricultural crops, and hence, the FCUs are distributed accordingly to sets of 15 and 10 units, respectively. Also, it shows the piping network as well as the temperature-based control valves that are employed to control the temperature on the heat pumps and the storage tank. For example, when the heat pumps produce hot water (over 50°C) and it is circulated between the water storage tank and the heat pump, the temperature control valve will open to control the water temperature by increasing the amount of the cold return water. The energy system is designed to have 5 heat pumps and a wood pellet boiler, and their specifications, as provided by the suppliers, are presented in Table 1.
The combined system has been modelled in Aspen Plus and MATLAB. Aspen Plus is a Fortran-based process modelling software package that has been extensively deployed to model various processes, including energy systems. It has an inclusive database of chemical, thermal, and physical properties of several compounds and holds a library of commonly used operation blocks. The developed model will inform and dictate the farm's operational strategies.

Process Modelling Approach.
For the solar PV panels, the hourly based simulation of the power generation from the PV module can be obtained by the following equation [25][26][27]: where P STC (W) is the module maximum power at standard test conditions, G h and G STC (W/m 2 ) are the hourly global solar radiation and the solar radiation at the test conditions, 3 International Journal of Energy Research and α (%/°C) is the temperature coefficient. T c and T a (°C) are the cell and ambient temperatures, respectively, and F diss (-) is a loss factor that accounts for the power dissipation due to dirt, wires, module mismatch, and other losses. The site and the required weather data are acquired from a station that is located in the city of Naju, Korea. All the weather data can be found in the Supplementary Material (available here). The cell temperature T c (°C) is obtained from the following expression [25,28]: where NOCT (°C) is the nominal operating cell temperature and T a (°C) is the ambient temperature. The specifications of the adopted PV panels are shown in Table 2. The power produced from the PV array is utilized to operate the heat pumps, which are simulated along with the wood pellet boiler in Aspen Plus. The five heat pumps are simulated as a single heat pump with the overall cooling capacity of 325 kW th . In the heat pump cycle, the R410a refrigerant is utilized as the working fluid and is defined in Aspen Plus by its composition, which is a mixture of 50% difluoromethane (R-32) and 50% pentafluoroethane (R125) [30]. Also, since the wood pellets are not of unique chemical compound, it is defined in the software as a nonconventional solid component. Later in the process, it is recognized by the solver when decomposed into the conventional compounds of the ultimate analysis (i.e., C, H 2 , O 2 , N 2 , and S) [31]. Furthermore, depending on this analysis, along with the proximate analysis, the properties of nonconventional materials comprising the biomass and ash are defined, where DCOA-LIGT and HCOALGEN are models that have been used to determine the density and enthalpy, respectively. The different properties of the working fluids in Aspen Plus can be determined based on the selected property method. In this model, a combination of the selected methods, including REFPROP, Peng-Robinson with the Boston-Mathias alpha function (PR-BM), and STEAM-TA, is utilized to estimate the properties of the refrigerant, air and chemical compounds, and water, respectively.  International Journal of Energy Research Figure 3 shows the process flow diagram of the energy system during the winter mode, which consists of the heat pump, the storage tank, the wood pellet boiler, and the FCUs. In this system configuration, the air-source heat pump works on the vapor compression cycle, and it is utilized to heat the water that is circulated to the storage tank. The refrigerant of the heat pump is compressed in the compressor with a pressure ratio of 3.5, where the expansion pressure is 6.5 bar, and this is within the typical operation ranges of the refrigeration cycle [32,33]. In addition, the isentropic and mechanical efficiencies of the compression work are 75% and 95%, respectively. The specifications of the evaporator and condenser heat exchangers are a vapor fraction of 1 and a subcooling with 3 degrees, respectively. Accordingly, the heat pump cycle with respect to the saturation lines of the refrigerant is plotted on the P-h diagram as shown in Figure 4. For the summer mode, another Aspen file is created in which the two heat exchangers of the heat pump are exchanged, such that the circulating water of the storage tank is connected to the evaporator instead of the condenser.
The sizing of the system is based on the total cooling capacity of the heat pumps as mentioned in Table 1, which is 325 kW th for the evaporator, and therefore, the flow rate of the refrigerant is determined to be 7250 kg/h. Based on this capacity and with the assistance of the utility tool in Aspen Plus, the corresponding flow rates of the air over the evaporator and the water through the condenser are 44.44 kg/s and 5.5 kg/s, respectively. The water is implemented to extract heat from the condenser of the heat pump in winter and to reject the heat to the evaporator in the summer mode. Here, a thermal storage tank has been incorporated into the design to store the energy of the heated or the cooled water, and it is modelled as a crossflow heat exchanger. The heat exchanger is pinched to 5°C at the cold side. This storage tank, in return, heats up or cools down another water stream that circulates to the FCUs in the winter and summer, respectively.
The COP of the heat pump is evaluated as the ratio between the amount of heat transferred to the water in the condenser (wintertime) or from the water in the evaporator (summertime) _ Q water (W) and the compressor power _ W comp (W), namely, This compressor work is provided by the solar PV array in both the cases of heating and cooling; however, this energy source is not available all day, and this is where the  International Journal of Energy Research role of the wood pellet-fired boiler augments the reliability of the system. Furthermore, the boiler is also ready to engage into the system when the ambient air temperature is too low during the wintertime. The wood pellet boiler is simulated firstly as a complete combustion of the pellets by following the equilibrium modelling approach. Secondly, the hot flue gas from the combustor is to be directed to a heat exchanger to boil the water for heating. It is worth noting that the wood pellet is defined as a nonconventional component and, therefore, its heat of combustion is to be defined in Aspen Plus as a pure component parameter. The heating value and the ultimate and proximate analyses, as received from a local supplier in the farm, are reported in Table 3.
The simulation procedure of the biomass combustion starts with the decomposition of the biomass in a RYield reactor (DECOMP) to its conventional components, i.e., carbon, H 2 , N 2 , O 2 , S, and H 2 O, in addition to the ash fraction. This yield is worked out through a Fortran calculator to preserve the mass balance between the biomass and the yielded components. The ash content is to be separated from the main material stream before the combustor, where the combustion reactions take place following the equilibrium approach in the modelling. This equilibrium condition is achieved by attaining the minimum Gibbs free energy in the RGibbs reactor (BURN), where the combustible components along with the air are introduced. The amount of air that is required to completely burn the fuel is estimated based on the stoichiometric air to fuel mass ratio ðA/f Þ stoich which can be calculated as follows [34]: where the values of C, H, S, and O are the dry ash-free mass fractions of the biomass and can be obtained from the ultimate analysis in Table 3. It is considered that the air is introduced with an excess ratio of 1.8, which is within the applicable range of boilers [35,36], and it is calculated in a calculator block (AFR-BURN). Thus, the mass flow rate of the air, in accordance with the mass flow rate of the wood pellets, is estimated and applied to the air stream flow. For the heat integration between the decomposition and combustion of biomass, a heat stream (Q-DECOMP) is connected between the RYield reactor (DECOMP) and the RGibbs reactor (BURN). The heating capacity of the boiler is 454 kW with an efficiency of 82.5%, and the wood feeding rate is 106.5 kg/h. The flue gas is sent to a common separator to remove the ash stream, and then, it is directed to the boiler heat exchanger (BOILER). The specification of this heat exchanger is that the exit temperature of the flue gas is 150°C and hence the boiling water is fed by a flow rate of 0.1416 kg/s to produce saturated water vapor. It is important to point out that the storage tank is set up to serve the heat pumps and hence it is connected to the water circuit that exchanges the heat with the heat pump as can be seen in Figure 3. Since the heat pump is the main component for the cooling and heating operation, the storage tank is not connected to the wood pellet boiler, which operates with a stable and consistent performance, and it is less dependent on the ambient condition. The storage tank is modelled in Aspen as a crossflow heat exchanger, where its specification is set to "Hot outlet cold inlet temperature difference = 5°C" for the winter mode and "Hot inlet cold outlet temp. difference = 5°C" for the summer mode. In the smart farm energy system configuration, the water from the storage tank or the boiler is provided to the FCUs, which have been incorporated into the design to perform the sensible heating or cooling of the air that is circulated for conditioning the farm. Two main Aspen Plus blocks have been utilized to model this air-conditioning system, namely, a fan block (FANS) with a multiplier block as the design is to include 25 fan units and a heat exchanger (SFARM).  The fan block draws the ambient air at a flow rate of 0.444 kg/s and blows it over the hot or cold tubes containing the flowing water from the storage tank or the boiler. Given that the temperature of the boiler water is much higher than the storage tank temperature in winter, the specification of the FCUs' heat exchanger has two configurations in order to avoid a simulation error. When the heat pump is the source of hot water, the specification of the heat exchanger is set to "Hot outlet cold inlet temperature difference = 10°C " and when the boiler is engaged, the specification is chan-ged to "Hot/Cold outlet temperature approach = 5°C". The control of this conditional specification of the heat exchanger and any control of the variables in the combined system are coded in MATLAB.
The operation of the model starts by importing the weather data, which was collected over the year of 2020 on an hourly basis. The power of the PV array is calculated, and then, MATLAB takes control of running the Aspen files. The control of the switch between the heat pump and the boiler is achieved by MATLAB through the stream selector   block "SELECTOR" by connecting the required stream, either the one from the storage tank or the one from the boiler as can be shown in Figure 3. A connection has been established using the MATLAB code to operate and control the Aspen operation through an interchangeable data bridge between the two software platforms using the COM server that is called "actxserver." This server creates an automated object linking and embedding (OLE) on the local system using the program identifier of Aspen Plus, and this is named "Apwn.Document.37.0" [37]. The simulation is carried out for one year on an hourly basis, i.e., a total of 8760 hrs, where MATLAB selects the targeted Aspen file (winter or summer) based on the ambient temperature (T amb ).
The calculation algorithm and the interaction between Aspen Plus and MATLAB is illustrated in Figure 5. As shown, the main condition to operate the system in the heating mode (wintertime) is when the ambient temperature is below or equal to 18°C, while the cooling mode is operated when the ambient temperature is above or equal to 28°C. Following this, the corresponding Aspen file is opened, and the ambient temperature is set to all inlet air streams. In the heating mode of the simulation, the selection of the boiler over the heat pump is based on two conditions, namely, when the solar power ðP solar Þ = 0 kW or T amb ≤ 0°C. It is worth mentioning that when any solar energy is available, the grid interacts with the PV array either by covering the shortage of power or receiving the surplus energy produced.

Economic
Analysis. An economic assessment is carried out based on the simulation data, which, in principle, can be achieved by calculating the economic indicators, such as the levelized cost of energy (LCOE), the annualized total cost, or the net present cost of the system. In this study, the LCOE, which is the unit cost of energy supplied to the system over its life cycle, has been deployed as the economic indicator to present the viability of the plant over its entire life. It is the ratio of the net present cost (NPC) of energy to the total energy demand (E) and is calculated as follows [38]: The NPC comprises the annualized value of the capital expenditures (ACAPEX) and the annual operation expenditures (OPEX); the latter consists of the maintenance and variable costs. This can be expressed as follows:   International Journal of Energy Research CAPEX is the equipment cost with any related installation and acquisition costs, and if distributed over the years of operation under a certain discount rate, the value of ACAPEX would be a one-year fraction of the total CAPEX. However, as the name suggests, NPC considers the ACA-PEX after being converted back to the initial time of purchase. Following a discount rate of (i d ) and a project lifetime of (n) years, the ACAPEX can be evaluated by multiplying the capital recovery factor (CRF) to the CAPEX as follows [39,40]: For this study, the discount rate and the lifetime are assumed to be 10% and 20 years, respectively. In order to evaluate the capital expenditures (CAPEX), the utilization of a bottom-up technique is followed, in which the cost components are evaluated as fractions of the purchase cost of the equipment (PCE) [41,42]. This approach considers the direct costs, such as the installation costs, as well as the indirect costs of the system components as shown in Table 4.
The price of a unit with specific size can be estimated in accordance with a reported price for a given unit size in a certain year. To consider the effect of scale on the purchase cost, as well as the price change over time, the scaling exponent correlation and the chemical engineering plant cost index (CEPCI) are utilized, and this can be expressed in the following form [43]: where C is the actual current cost of the equipment, C 0 is the unit base cost at the reference year, S and S 0 are the sizes of the actual unit and the reference unit, respectively, and f is an empirically determined exponent.  Table 5. For the OPEX, the employed fixed and variable costs of the different equipment are presented in Table 6.

Environmental Analysis.
Another major factor to take into account when introducing a low-carbon energy system is the environmental impact of the hybrid system. This can be achieved by comparing the CO 2 emissions of the system to that of a conventional heating system. Accordingly, the CO 2 emission reduction ratio (CO 2 ER) is presented, and it can be evaluated from the following expression [53]: where X i CO 2 is the emission factor of the conventional heating system (i) and X Q CO 2 and X e CO 2 are the emission factors of the proposed energy system for the heat load, Q CCHP , and for the electricity consumed by the system, E CCHP , respectively. Q Conv is the total heat supply when generated by a conventional system. It should be noted that in our case, the Q CCHP refers only to the boiler heating generation as the heating/cooling produced by the heat pump is incorporated in the energy of the CCHP system (E CCHP ). Figure 6 depicts the ambient temperature, the solar irradiation, and the corresponding power generated from the PV modules over one year (8760 hours). It can be seen that the temperature fluctuates between -7.5°C (in winter) and as high as 31.5°C (in summer), while the global horizontal irradiation in this location reaches 1018 W/m 2 in the summer months. Accordingly, the PV system generates up to 200 kW during winter, while in the early months of the summer, a peak generation of about 339 kW is achieved.

Technical Assessment of the System.
The system of the PV panels is connected to the grid as a backup source of electricity and as a receiver of surplus  9 International Journal of Energy Research power. The simulation results of the system show the power consumption and how much electricity is being exchanged between the PV panels and the grid, and this is exhibited in Figure 7. Most of the electricity that is received from the grid is in the low temperature period of time, while in most of the summertime, the electricity from the PV panels is exported to the grid. The total annual electricity consump-tion by the farm is 367. 16 MWh. The solar panels provide 234.12 MWh in real-time operation while the remaining 133.04 MWh is covered by the grid. However, when considering the annual accumulated electricity generation from PV arrays, the system exports 332.62 MWh to the grid, and hence, the energy system is a net exporter of electricity with 199.58 MWh.

10
International Journal of Energy Research The generated power from the solar arrays, backed by the grid, is directed to run the heat pump as the primary option to meet the demand of the smart farm. The COP of the heat pump is estimated at about 4 for the heating mode and at about 3 for the cooling mode as shown in Figure 8. The heat pump is running at a constant COP value in both the heating and cooling modes, and this is because of the assumed fixed pressure ratio and the fixed design specifications of the two heat exchangers (evaporator and condenser) as discussed in Section 3.1. Nevertheless, Figure 8 shows that in winter, at the relatively intense cold temperatures (<4°C), the COP expectedly drops down to approximately 3.4 as there is not enough heat available from the air to achieve the design specification of the evaporator, i.e., the refrigerant vapor fraction = 1. When the ambient temperature is between 18°C and 28°C, the heat pump is shut down, and hence, no relevant COP values are presented in Figure 8. The average COP over one year, also known as the seasonal performance factor, is calculated as 3.85.
The predicted average temperature and the average heat load of the farm over one year on a daily basis are shown in Figure 9. It can be seen in Figure 9(a) that the temperature is, in the vast majority of the cases, above 20°C and mostly fluctuates around 24°C and these conditions are satisfactory for the farm operation [18]. In Figure 9(b), the average farm heat load is at the highest value in the winter period at about 220 kW when the ambient temperature is low, and this load decreases as it moves towards the summer period. Also, it is noticeable that for about 100 days (from about the 160 th day to about the 260 th day), the heat load is negligible, and this is because the ambient temperature is in the targeted range of the operating temperature. In the summer, the cooling load of the farm has been satisfied by the heat pump for limited operational days, which are represented by the red dots in Figure 9(b).

Economic Evaluation.
Based on the operating conditions and the simulation data, the economic evaluation of the energy system has been carried out. Figure 10 shows how the CAPEX is broken down between the different components of the system. It can be seen that the PV panels and the biomass boiler share most of the CAPEX with about 80% and with almost equal share of this amount. In addition, the heat pump is responsible for 16%, and the balance is attributed to the FCUs and the thermal storage tank. From a CAPEX perspective, the PV-heat pump combination is more expensive than the biomass boiler as modes of operation for heating. Nevertheless, the purchase cost of the equipment (PCE) or CAPEX cannot be the one and only economic criterion of any system, but rather, the OPEX is of the same vitality as an economic evaluating parameter.
As depicted in Figure 11, the fixed and variable costs of the biomass boiler (i.e., the operation and maintenance cost and the wood pellet cost, respectively) are responsible for the majority of the expenditures by a high margin, i.e., about 68% of the total OPEX, whereas the share of both the heat pump and the PV panels is only 25%, and hence, this combination is less OPEX-intensive than the boiler. The rest of the OPEX is 6% for the FCUs and 1% for the thermal storage tank.
The total annual energy demand by the farm is 1075.87 MWh th , which is utilized to estimate the ACAPEX and OPEX of the energy system. Accordingly, the LCOE has been estimated by using equation (5) at $0.218/kWh. The estimated LCOE lies within the range of the reported energy costs ($0.14-0.287/kWh) for greenhouses coupled to   11 International Journal of Energy Research geothermal heat pumps in the IPCC report "Summary for Policymakers of 2012" [54,55]. For comparison purposes, a grid-connected GSHP system in China for greenhouse heating has an estimated LCOE of $0.4/kWh [13].
Further, Figure 12 depicts that the energy system is CAPEX-intensive over the OPEX and the boiler system (due to the biomass cost) is the major cost contributor to   Figure 11: Breakdown of the OPEX between the system components.

12
International Journal of Energy Research the OPEX as the PV followed by the heat pump do not need major operation and maintenance costs. Also, it can be seen that the surplus electricity from the PV arrays over one year (199.58 MWh) can reduce the LCOE by $0.018/kWh. This represents a reduction by 23.1% of the net annual contribution by the PV system (i.e., $0.078/kWh).

Sensitivity
Analysis. The calculation of the LCOE is based on the applied cost for the different system components, but this cost may have a considerable margin of uncertainty. Hence, a sensitivity analysis can give a good indication on the contribution of each factor into the LCOE. Figure 13 illustrates the share of the most important cost components that would affect the LCOE, such as the pur-chase prices of the PV panels, the heat pump, and the boiler. The adopted cost for the PV corresponds to commercial/ industrial scale (i.e., $1700/kW); we have assumed as lower bound for this cost the utility-scale PV cost in Korea, i.e., $1300/kW [44]. A similar increase by $300/kW is assumed as the possible higher price. For the heat pump and the boiler PCE, a variation of ±50% is considered. Also, two main variable costs are considered: the biomass price and the electricity selling price. A wide range has been considered for the feedstock cost to capture the relatively high uncertainty of wood pellet cost. Similarly, a wide range of values is considered for the electricity costs. The boundaries of change for each investigated cost component are shown in Table 7. 13 International Journal of Energy Research As shown in Figure 13 and due to the relatively high boiler purchase cost, the change in its price leads to a noticeable variation in the LCOE of around ±20.8%. The variation in the investment cost of the PV panels and the heat pump results in the LCOE to vary within ±6% and ±7.2%, respectively. Further, it can be observed that the LCOE is highly sensitive to biomass cost, and for the extreme case of $310/ tonne of wood pellets, the LCOE increases to $0.273/kWh (an increase of 25.2%). Hence, the establishment of a longstanding biomass procurement deal with determined quantities and charges should be chased with an agricultural or forest management company. This highlights the importance of establishing the long-term supply of biomass at a fixed price. In addition, the electricity selling price to the grid can contribute to the reduction of the energy price down to $0.2/kWh, and this can be achieved as the electricity exported may be eligible to receive renewable energy tariffs. It can be noticed that the LCOE is more sensitive to the boiler costs, and hence, the combination of the PV panels with the heat pump should have the priority over the boiler to operate when the solar energy is available.

Carbon Footprint.
In order to study the carbon footprint of the smart farm energy system in comparison to other conventional systems, the emission factors for the different heating technologies are tabulated in Table 8. The emission factors of the PV array (0.0826 kgCO 2 -eq/kWh) and the grid electricity in Korea (0.666 kgCO 2 -eq/kWh) are taken from the SimaPro software, while for the boilers, the emission factors have been taken from relevant literatures. In this study, it is assumed that the conventional electric-based heating is to run a heat pump with the same COP as that of the smart farm's heat pump in order to satisfy the total energy demand of the farm. Also, as discussed in Section 4.1, the majority of the energy demand is for heating, and therefore, the other conventional energy systems are different types of boilers with the cooling load neglected. This means that the proposed new energy system can achieve even higher CO 2 reductions, and this is very important. The calculated emission factor of the proposed system is 0.056 kgCO 2 -eq/kWh. This is much less than the hybrid system that has been investigated for another Korean smart farm with a wood pellet boiler and a grid electricity-based hydrothermal heat pump, i.e., 0.424 kgCO 2 -eq/kWh [24]. This highlights the sustainability gains of our proposed design and the necessity to design systems with low reliance on the grid. Also, it can be seen that the proposed smart farm can achieve a significant emission reduction of between 85.42% and 88.94% when replacing a grid-connected heat pump or fossil fuelbased heating boilers. The achieved emission reduction is higher than the value that has been reported by Zhou et al. [14], which was 35.5%. This value was attained when using a GSHP connected with a fuel oil heater in a hybrid system as a compared to the oil heater only. This shows the significant environmental benefits of employing the pellet boiler and the solar PV panels, as in the current study.
Compared to the wood chip boiler, the hybrid energy system has higher emission factor; however, the smart farm still reports a small CO 2 ER because of the heat pump COP being high enough (seasonal performance factor of 3.85). This implies that the overall amount of energy needed, to achieve the same heating load, is less than the boiler. However, when compared to the wood pellet boiler, the smart farm emits more CO 2 . It should be noted that the greenhouse energy system is a net exporter of electricity, and this has not been considered in the current carbon footprint analysis. Therefore, the net emission of the smart farm could have been reduced even further if we were to consider that the green-exported electricity displaces grid electricity. Figure 14 illustrates the amount of the CO 2 emission that is being released from the smart farm in comparison to the fossil fuel-based heating, biomass-based heating, and electricity-based heating. The smart farm emits only 60.2 tCO 2 /y, while the emissions from the fossil fuel-based heating systems can go up to 544.5 tCO 2 /y. In comparison to the biomass-based heating boilers, the green energy system is still comparable, and even the smart farm energy system runs more hours for cooling. Since the smart farm energy system is a combination of a wood pellet boiler and PV panels, the carbon footprint of the smart farm is slightly  higher than the single wood pellet boiler by 11.62 tCO 2 /y. This is because of the relatively higher emission factor of the PV array.

Conclusion
This study presents a techno-enviroeconomic investigation of a decentralized renewable CCHP system for a smart farm located in South Korea. This energy system comprises solar PV panels, ASHP, and a wood pellet boiler in a hybrid configuration for greenhouse farm application. A novel dynamic modelling approach has been developed and applied to simulate the operation of the energy system on an hourly basis. The model implements an interchangeable operation of the PV-heat pump and the wood pellet boiler. The major and important outcomes of this study are as follows: (i) The synchronized operation of the hybrid energy system model has provided reliable dynamic operation to maintain the energy demand of the farm with any surplus electricity being exported to the grid (ii) Technically, the hybrid energy system has achieved an average daily temperature in the farm around 23.9°C and fluctuates with a standard deviation value of 2.16°C, while the seasonal performance factor of the heat pump is estimated at 3.85 over one year (iii) By utilizing the technical thermodynamic data, the LCOE is estimated at $0.218/kWh, and the total annual energy demand of the farm is 1075.87 MWh th . This implies that the system without subsidies cannot compete with conventional energy systems (iv) The priority for prime energy is given to the PVheat pump combination over the wood pellet boiler due to the higher energy efficiencies (v) For the given economic assumptions, as the energy dependency moves from the PV to the boiler, the higher the cost of energy becomes, and this is attributed to the high running costs of the boiler (vi) It is worth noting that the PV panels help in reducing the LCOE by $0.018/kWh through exporting the surplus electricity to the grid over the year (vii) From an environmental perspective, the carbon footprint of the greenhouse energy system is significantly lower than the conventional heating systems, such that the achieved CO 2 ER varies between 85.42% and 88.94% when compared to the electric-based or fossil fuel-based heating systems (viii) The quantity of CO 2 that could be emitted by the greenhouse energy system is 60.2 tCO 2 /y, which is equivalent to an emission factor of 0.056 kgCO 2eq/kWh The current study presents reliable design and assessments that can provide meaningful insights for the transition of the farming sector to renewable decentralized energy solutions. The resilience of such hybrid energy systems can be increased by including waste to energy technologies, such as gasification and biomethanation along with different energy storage technologies.

Data Availability
Modelling data will be available upon request.

Conflicts of Interest
The authors declare that there is no conflict of interest that could be perceived to affect the objectivity or neutrality of the manuscript.