Comparative Analysis and Design of a Solar-Based Parabolic Trough–ORC Cogeneration Plant for a Commercial Center

: This paper performs technical, economic and environmental feasibility analyses of two di ﬀ erent solar cogeneration plants, consisting of a solar system (a parabolic trough collector ﬁeld coupled with thermal energy storage), an Organic Rankine Cycle (ORC), and mechanical chillers, that should cover the electrical and cooling demands of a commercial center located in Zaragoza (Spain). System A is hybridized with an auxiliary biomass boiler that complements the solar system’s thermal production, providing a constant heat supply to the ORC, which operates at full load during the operating hours of the solar system. In contrast, system B is not hybridized with biomass, so the ORC is fully driven by the solar system, operating at partial load according to the solar resource availability. Both systems are connected to the electrical grid, allowing electricity purchases and sales when needed. The design procedure involves the sizing of the equipment as well as the modelling of the hourly behavior of each system throughout the year. The physical analysis is complemented by an economic assessment, which considers investment and variable costs, as well as an estimate of the signiﬁcant environmental beneﬁts of the proposed plants. The solar plants are compared to a conventional system in which all the electrical consumption is covered with electricity purchased from the grid. The costs of the electricity produced by systems A and B are estimated at 0.2030 EUR / kWh and 0.1458 EUR / kWh, which are about 49% and 7% higher than the electricity purchase price in Spain (0.1363 EUR / kWh). These results indicate that while none of the solar plants are presently competitive with the conventional system, system B (without biomass hybridization) is actually closer to economic feasibility in the short and medium term than system A (with biomass hybridization).


Introduction
The COVID-19 pandemic puts the climate crisis into perspective as governments around the world rapidly act to safeguard the health, safety, and wellbeing of their citizens in what represents one of the greatest challenges to face society in generations. This public health crisis must be considered in the larger context of climate change and global warming, which also requires urgent and strong measures. In the wake of the pandemic, the world must confront the task of kickstarting the economy while also consolidating climate change commitments. The International Energy Agency (IEA) [1] and the International Renewable Energy Agency (IRENA) [2] regard energy efficiency, renewable energy, and energy storage as key pillars of the energy transition, as they contribute to environmental benefits while also fostering economic growth, the creation of jobs, human wellbeing, the security of the energy supply, and reductions in energy costs. To reach these objectives, it is crucial to design new energy supply systems that are carbon neutral, resilient, and less polluting. refrigeration cycle [30], and a reversible ORC/heat pump [43]. Moreover, different operation modes are proposed, such as the ORC plant operating in cogeneration mode in the winter and in trigeneration mode in the summer [41], and producing power and heating in the winter and power and cooling in the summer [42].
Solar-driven ORC plants with and without biomass hybridization have been extensively analyzed in the literature. Zhao et al. [38] explored the configuration effects of CCHP systems employing PTCs, single-effect absorption chillers for cooling production, and heat exchangers for heating production. The analyzed configurations are classified into sequential (cascaded PTC field, ORC, cooling, and heating cycles) and parallel (ORC, cooling, and heating cycles supplied by the PTC field), and the results demonstrate that the sequential configuration is more compact and presents better thermodynamic and economic performances. Roumpedakis et al. [44] proposed a design optimization and performance assessment of small-scale solar-driven ORC plants in the southeast Mediterranean region (Cyprus, Greece, and Turkey), considering, in addition, different working fluids and solar thermal collectors (flat plate, evacuated tube, and PTC). A multi-objective optimization was developed using a genetic algorithm to maximize exergy efficiency and minimize the payback period. While the exergy efficiencies of the collectors were about the same, the PTCs resulted in better economic performances, with the lowest payback periods. Pantaleo et al. [45] carried out a thermodynamic and economic analysis of a novel hybrid solar-biomass CHP system composed of a topping PTC-biomass externally fired gas turbine and a bottoming ORC. Different system configurations were considered, including power-only and cogeneration modes, under the Italian energy policy scenario. The results indicated the low economic profitability of PTC integration relative to the biomass-only plant. Wu et al. [37] assessed the economic and environmental performances of three solar-biomass hybrid CCHP systems, namely CCHP-ORC (the ORC is driven by waste heat from a gas engine), solar-ORC (the ORC is driven by a PTC field), and biomass-ORC (the ORC is driven by a biomass boiler). A multi-objective optimization model based on mixed-integer linear programming was developed to determine the optimal system configuration and operating strategies for a set of hotel and office buildings in Shangai, China. The results obtained suggested that the solar-ORC had the best economic performance, whereas the biomass-ORC enjoyed the best environmental benefit. Petrollese and Cocco [32] proposed a novel robust optimization approach for the preliminary design of solar-based ORC systems that included the effects of heat source and heat sink variations, especially regarding the availability of solar radiation and the variations in ambient temperature. The proposed methodology was tested by considering an existing medium-sized solar ORC facility in Italy. Bellos et al. [40] investigated a solar and biomass polygeneration system for the climate conditions of Athens (Greece) from energy, exergy, and financial viewpoints. The system includes a PTC, a storage tank, a biomass boiler, an ORC, and a vapor compression chiller, and produces electricity and cooling, as well as heating at two temperature levels.
The interest of solar and biomass ORC plants has also been highlighted by the commissioning of various pioneering facilities. In this regard, Denmark already has cogeneration systems based on ORC integrated with CSP and biomass in operation, demonstrating their techno-economic feasibility at a district level. Examples include (i) the Brønderslev facility, with 16.6 MW t of installed capacity and CSP aperture area of 26,929 m 2 [46] and (ii) the Taars facility, with 2.5 MW t of installed capacity and the support of flat-plate solar thermal collectors [47]. Recently, a pilot micro-trigeneration plant with an ORC, biomass and a linear Fresnel solar field started operating in Spain, commissioned by the European project Tribar [48]. The plant produces 50 kW of electricity, 400 kW of heat, and 35 kW of cooling.
Our literature survey indicates the great interest in hybrid solar PTC and biomass ORC-based trigeneration systems. These systems represent an attractive solution that not only promotes the combined production of energy services (electricity, heating, and cooling), but also combines two energy sources that complement each other perfectly. In order to achieve the potential benefits of these systems, it is necessary to perform an adequate design procedure. Nevertheless, many reported studies Energies 2020, 13, 4807 5 of 29 make major simplifications or fail to incorporate the dynamic conditions of the system into the analysis, such as the partial load operation of the ORC, and the daily and hourly variability of solar thermal production throughout the year. Additionally, most studies propose cogeneration (power and heat) or trigeneration (power, heat, and cooling) systems, but few focus on the combined production of power and cooling only, with very low or no heat demand, which is a common situation in some kinds of buildings, such as commercial centers, depending, among other factors, on their geographical location.
In this direction, the present paper investigates the technical, economic, and environmental feasibility of two different solar PTC-ORC plants that cover the electricity and cooling demands of a commercial center in Zaragoza (Spain). Both systems include a solar system, composed of a PTC field integrated with TES tanks, an ORC, and vapor compression mechanical chillers. One of the systems is hybridized with a biomass boiler, while the other is not. In this work, the preliminary design procedure and the operation of the system include: (i) a detailed physical simulation of the solar system (PTC field integrated with TES tanks), using the System Advisor Model (SAM) software [49], which captures its dynamic behavior throughout the hours of the day and the days of the year, such as the variations in ambient temperature, solar radiation, and thermal energy production, among others; and (ii) the modelling of the different pieces of equipment based on technical data from commercially available devices, such as the ORC's part-load performances, obtained from the manufacturer, Turboden [50], among others. Solar ORC plants are proposed with the aim of significantly contributing to the decarbonization of energy production, achieving a very high renewable fraction. In this regard, the plants' main energy resources are renewable energies (i.e., solar and biomass), and the auxiliary resource is the grid electricity, which contains a non-renewable fraction, whereas studies in the literature commonly consider the renewable resource as a support to partially displace non-renewable conventional fuel.
Even though the case study results are unique to the problem at hand, they will give an insight into the interest and the level of economic maturity of this type of system in the short and medium term in the Spanish context. Furthermore, the methodology described herein can be applied to different consumer centers (e.g., office, hotel and hospital buildings) in different geographical locations, and the ORC plant can be modified accordingly, for instance, to produce heat as a final product or to incorporate absorption chillers for cooling production.

Energy Demands of the Commercial Centre
The proposed energy systems must attend to the electricity and cooling demands of a commercial center located in Zaragoza, Spain. Cooling is required for air conditioning throughout the year, even during winter months, to maintain comfortable indoor temperature conditions. The annual thermal cooling demand amounts to 6412 MWh t and is satisfied by chilled water from 7 to 12 • C, produced by electricity-driven vapor compression refrigeration chillers. The annual thermal cooling demand is distributed throughout the year, considering the monthly thermal cooling demand data simulated for a similar consumer center and the hourly occupation rates on working days and holidays [51]. The monthly distribution of thermal cooling is shown in Figure 1. The occupancy of the consumer center is 18 h per day, corresponding to 12 h where it is open to the public (from hours 09 to 21) plus 6 h for services and supply (from hours 06 to 09 and 21 to 00). The occupancy profiles for working days and holidays are shown in Figure 2.  The hourly thermal cooling demand throughout the year is shown in Figure 3 for working days and holidays. Thermal cooling loads are notably higher on holidays, when the maximum occupancy is reached, than on working days. Moreover, they are significantly higher in the summer than in the winter months.   The hourly thermal cooling demand throughout the year is shown in Figure 3 for working days and holidays. Thermal cooling loads are notably higher on holidays, when the maximum occupancy is reached, than on working days. Moreover, they are significantly higher in the summer than in the winter months. The hourly thermal cooling demand throughout the year is shown in Figure 3 for working days and holidays. Thermal cooling loads are notably higher on holidays, when the maximum occupancy is reached, than on working days. Moreover, they are significantly higher in the summer than in the winter months. The annual electrical demand required for lighting, services, etc. (excluding the electricity consumption for chilled water production in the vapor compression refrigeration chillers) amounts to 1370 MWhe. The daily electricity demands on a working day and on a holiday are 3465 kWhe and 4264 kWhe, respectively, and are the same throughout the year. The hourly distributions for each type of day are estimated according to the load profiles shown in Figure 4.   The annual electrical demand required for lighting, services, etc. (excluding the electricity consumption for chilled water production in the vapor compression refrigeration chillers) amounts to 1370 MWh e . The daily electricity demands on a working day and on a holiday are 3465 kWh e and 4264 kWh e , respectively, and are the same throughout the year. The hourly distributions for each type of day are estimated according to the load profiles shown in Figure 4.  The annual electrical demand required for lighting, services, etc. (excluding the electricity consumption for chilled water production in the vapor compression refrigeration chillers) amounts to 1370 MWhe. The daily electricity demands on a working day and on a holiday are 3465 kWhe and 4264 kWhe, respectively, and are the same throughout the year. The hourly distributions for each type of day are estimated according to the load profiles shown in Figure 4.

Description of the Analyzed Systems
A simplified diagram of the solar-based cogeneration systems analyzed herein is depicted in Figure 5. System A consists of a solar system, composed of a solar parabolic trough collector field with thermal energy storage, integrated with an ORC and hybridized with a biomass boiler. It is assumed that the ORC only operates at full load and the biomass boiler complements the thermal output of the solar system whenever the solar resources are low, guaranteeing a constant thermal energy supply to the ORC, both in quantity and temperature. The system is in operation whenever there are enough solar resources available to operate the ORC, that is, the biomass boiler only operates during the operating hours of the solar system. In contrast, system B consists of a solar system integrated with an ORC without the biomass boiler hybridization. In this case, the ORC is exclusively driven by solar thermal energy and is allowed to operate at partial load whenever the solar resources are not enough to allow for full load operation.

Description of the Analyzed Systems
A simplified diagram of the solar-based cogeneration systems analyzed herein is depicted in Figure 5. System A consists of a solar system, composed of a solar parabolic trough collector field with thermal energy storage, integrated with an ORC and hybridized with a biomass boiler. It is assumed that the ORC only operates at full load and the biomass boiler complements the thermal output of the solar system whenever the solar resources are low, guaranteeing a constant thermal energy supply to the ORC, both in quantity and temperature. The system is in operation whenever there are enough solar resources available to operate the ORC, that is, the biomass boiler only operates during the operating hours of the solar system. In contrast, system B consists of a solar system integrated with an ORC without the biomass boiler hybridization. In this case, the ORC is exclusively driven by solar thermal energy and is allowed to operate at partial load whenever the solar resources are not enough to allow for full load operation. Solar radiation, more specifically solar Direct Normal Irradiation (DNI), is harnessed by the parabolic trough collector field and transferred to the heat transfer fluid (HTF) as thermal energy. The solar collector field is connected to a TES in order to prolong the operation time and reduce the production variability, partially mitigating the non-simultaneity between the time period of solar resource availability and consumption. The assembly is connected to a heat exchanger (the ORC Figure 5. Simplified diagram of the analyzed parabolic trough-ORC-biomass cogeneration system. Adapted from Wagner and Gilman [52]. Solar radiation, more specifically solar Direct Normal Irradiation (DNI), is harnessed by the parabolic trough collector field and transferred to the heat transfer fluid (HTF) as thermal energy. The solar collector field is connected to a TES in order to prolong the operation time and reduce the production variability, partially mitigating the non-simultaneity between the time period of solar resource availability and consumption. The assembly is connected to a heat exchanger (the ORC evaporator) where the thermal energy of the HTF is transferred to the ORC's working fluid, allowing for the conversion of thermal energy into mechanical energy in the ORC. An alternator directly connected to the ORC turbine shaft transforms the mechanical energy into electricity, which is used to cover the electricity demand required for lighting and services, and to drive the vapor compression refrigeration chillers that produce the low-temperature thermal energy required for the cooling demand. The system is connected to the electrical grid, which allows electricity purchase when the demand exceeds the production, and electricity sale when there is an excess of electricity production. Thus, purchasing electricity from the grid ensures an appropriate electricity supply, and selling electricity to the grid generates revenue to improve the economic return of the investment.
Both systems A and B have been sized so that the annual electrical energy production is similar to that of a conventional system that covers the energy services with electricity purchased from the grid. The procedure for sizing the plant components of both systems is explained in detail in previous works developed by the authors [53,54]. The first step was the selection of the vapor compression refrigeration chillers that must cover the thermal cooling demands of the commercial center. In this way, the total annual electrical demand of the commercial center (electricity required for cooling and for other services) was determined. With this information, the ORCs for systems A and B that should produce the annual electricity consumed by the commercial center were selected and the heat input required by each ORC was thus determined. This information was used to size and design both solar systems, as well as, in the case of system A, the auxiliary biomass boiler, which complements the solar field's thermal energy production in order to guarantee a constant thermal energy supply to the ORC. The appropriate sizing of the solar field area and TES capacity was determined for each system in order to minimize heat loss and maximize the solar thermal energy production with the quality (temperature level) required by each ORC. An iterative procedure was applied to properly match the thermal energy produced in each solar system and the heat input required by each ORC. The solar systems were assumed to be installed in an area next to the consumer center. The available surface area was sufficiently large so that it did not constitute an active restriction to the problem.
The technical data of the plant components selected are shown in Table 1. The following subsections describe the different plant components and explain the physical model developed to design and study the behavior and performance of both systems A and B throughout the year.

Solar System
The solar system was modelled and simulated using the System Advisor Model (SAM) [49], which is an open-source model oriented towards the feasibility assessment of renewable energy projects. A brief description of the mathematical model used for the simulation of the solar system is presented in the following paragraphs. The SAM simulation engine, integrated with TRNSYS [55], provides, among other things, a detailed hourly simulation performance of concentrating parabolic trough solar collectors, piping and thermal energy storage based on the physical principles of thermodynamics and heat transfer, considering transient effects related to variable solar radiation during the day and throughout the year, the thermal capacity of the HTF in the field piping, the headers, and the balance of the whole system [56]. The solar field consists of solar collector assemblies (SCA) that are combined in one or more parallel loops ( Figure 5). Each SCA is composed of a parabolic collector and its receiver, as depicted in Figure 6. SCAs are connected in series in each loop in order to incrementally increase the HTF temperature to the designed outlet loop HTF temperature value. In this model, the SCA serves as the lowest level of discretization and they are considered as a single unit when calculating the model. Thus, the absorbed energy, losses, temperature, pressure drop, and other performance values are calculated independently for each SCA [49,56]. The solar field consists of solar collector assemblies (SCA) that are combined in one or more parallel loops ( Figure 5). Each SCA is composed of a parabolic collector and its receiver, as depicted in Figure 6. SCAs are connected in series in each loop in order to incrementally increase the HTF temperature to the designed outlet loop HTF temperature value. In this model, the SCA serves as the lowest level of discretization and they are considered as a single unit when calculating the model. Thus, the absorbed energy, losses, temperature, pressure drop, and other performance values are calculated independently for each SCA [49,56]. The energy balance of a control volume encompassing the HTF within the absorber tubes and piping for a single SCA is where the internal energy, U, term represents the change of energy of the control volume with time, m is the mass of the HTF in the control volume of the SCA, chtf is the specific heat of the HTF, L is the The energy balance of a control volume encompassing the HTF within the absorber tubes and piping for a single SCA is being ∂U ∂t = m·c ht f + (mc) bal,sca ·L ∂T ∂t (2) where the internal energy, U, term represents the change of energy of the control volume with time, m is the mass of the HTF in the control volume of the SCA, c htf is the specific heat of the HTF, L is the length of a single SCA, (mc) bal,sca is the thermal capacitance per unit of length of joints, piping, insulation and other SCA pieces of equipment that also experience thermal inertia effects, T in and T out are the inlet and outlet temperatures of the HTF to each SCA, . m ht f is the HTF mass flow rate, which is constant across the boundary of the SCA and . q abs is the rate of absorbed solar thermal energy in the receiver. By applying the previous energy balance, which is a first order differential equation, successively to each SCA, all the temperature values of the HTF in the solar field are obtained, considering that the HTF inlet temperature to each SCA is the outlet temperature of the previous SCA. In order to meet the design value of the outlet loop HTF temperature, the mass flow rate of the HTF is varied. If the calculated HTF mass flow rate in the loop is lower than the minimum accepted value, then the software resets the mass flow rate to the minimum allowable and the outlet loop HTF temperature will be lower than the design value. When the solar field absorbs more energy than required to reach the design value of the outlet loop HTF temperature, even with the maximum HTF mass flow rate accepted value, the system controller defocuses some collectors in order to reduce the harnessed solar energy and to avoid overheating the HTF [52,56].
The collector model determines the magnitude of incident solar flux on the receiver, considering constant and variable optical losses. Constant optical losses encompass tracking errors, imperfections in the collector geometry, reflectance and soiling of the mirrors. Variable optical losses depend on the solar position. The most important variable losses are the so-called cosine losses, proportional to the cosine of the incidence angle θ of the solar irradiation, which is the angle between the solar irradiation and the normal to the aperture plane, as shown in Figure 6. Other position losses considered are losses of reflected radiation at the end of the collector row (spillage losses), shadowing caused by other collector rows, available solar energy not being harnessed early in the morning and late in the evening, and an incident angle modifier accounting for the remaining variable optical losses [52,56].
The receiver is modelled as a one-dimensional energy flow, following the model developed by Forristall [59]. Only temperature gradient in the radial direction is considered significant. Figure 7a presents a diagram of a quarter of a cross section of the receiver. Temperatures T 1 to T 5 are calculated by applying thermal energy balance and a thermal resistance heat transfer model ( Figure 7b). The geometry of the receiver is specified as well as the glass envelope absorptance and transmittance. Concentrated solar irradiation flux passes through the glass envelope, where a small fraction is absorbed, and unabsorbed irradiation reaches the absorber tube. Thermal energy is transferred through the absorption tube wall to the HTF. Thermal losses occur in the absorber tube and glass envelope surfaces due to convection and radiation mechanisms [52,56].
The SAM model also calculates the pressure drop occurring in several field components, e.g., piping from and to the solar field headers, receiver tube piping, elbows and valves. The total HTF volume and surface area of the piping are considered, with the exception of the surface area of the receiver absorber. The model applies a specific heat loss coefficient to the piping surface area, which depends on its temperature, to calculate the thermal energy loss from pipes [52,56]. by applying thermal energy balance and a thermal resistance heat transfer model (Figure 7b). The geometry of the receiver is specified as well as the glass envelope absorptance and transmittance. Concentrated solar irradiation flux passes through the glass envelope, where a small fraction is absorbed, and unabsorbed irradiation reaches the absorber tube. Thermal energy is transferred through the absorption tube wall to the HTF. Thermal losses occur in the absorber tube and glass envelope surfaces due to convection and radiation mechanisms [52,56]. The SAM model also calculates the pressure drop occurring in several field components, e.g., piping from and to the solar field headers, receiver tube piping, elbows and valves. The total HTF volume and surface area of the piping are considered, with the exception of the surface area of the receiver absorber. The model applies a specific heat loss coefficient to the piping surface area, which depends on its temperature, to calculate the thermal energy loss from pipes [52,56].
As shown in Figure 5, we also consider the integration of TES into the system, in order to increase the stability of the solar field production. Both cold and hot tanks are designed to store the entire HTF volume. The tanks are assumed to be fully thermally mixed without stratification. The capacity of the tank is expressed in equivalent hours of the solar system operation at full load, that is, the number of hours that the storage system can supply energy at a design point. In terms of thermal energy, the total TES thermal capacity, Etes, is equal to the design point thermal requirement that the solar system should provide times the total number of desired storage hours. The volume of HTF, Vtes, required to match the TES energy content is calculated by Equation (3), considering average material properties for the hot and cold tank design temperatures [52,56]. As shown in Figure 5, we also consider the integration of TES into the system, in order to increase the stability of the solar field production. Both cold and hot tanks are designed to store the entire HTF volume. The tanks are assumed to be fully thermally mixed without stratification. The capacity of the tank is expressed in equivalent hours of the solar system operation at full load, that is, the number of hours that the storage system can supply energy at a design point. In terms of thermal energy, the total TES thermal capacity, E tes , is equal to the design point thermal requirement that the solar system should provide times the total number of desired storage hours. The volume of HTF, V tes , required to match the TES energy content is calculated by Equation (3), considering average material properties for the hot and cold tank design temperatures [52,56].
where ρ tes,ave and c tes,ave are, respectively, the average density and specific heat of the HTF, and the design temperature difference is equal to the hot solar field outlet temperature T sf,out minus the cold solar field return temperature T sf,in multiplied by the heat exchanger factor f hx , which is the ratio of the temperature difference on the TES side of the heat exchanger to the solar field side temperature difference. The hot and cold tanks are modelled separately using the variable tank model (type 39) in TRNSYS [55]. The heat exchanger connecting the solar field and the TES (see Figure 5) is modelled with the effectiveness-NTU method, considering the counterflow arrangement. SAM calculates the size of the heat exchanger from the solar system design values [52,56]. The input parameters for the simulation of the solar systems with SAM are described as follows. For both systems A and B, the solar collector model Siemens SunField 6 and the receiver Siemens UVAC 2010 were selected. All required data about geometry, materials and required technical specifications are available in SAM's libraries [49]. The selected HTF that flows through the receiver tube is Therminol 66. It is able to cover an extended operating range, from −85 • C to 400 • C, and its common applications are heating and cooling processes, exploiting the fluid as a thermal vector for transporting energy [60]. The solar field outlet temperature T sf,out and solar field return temperature T sf,in values, as well as the thermal power provided by the solar system, were selected in an appropriate range to operate the selected ORCs. The chosen collectors' orientation was the north-south direction, which collects more energy in summer than in winter. This alternative has been considered as the best choice because the application of this system aims to cover the cooling demand, which is higher in summer. In order to design the collector parabolic field, the model uses the design point at the DNI value (950 W/m 2 ) with the sun position at noon on the summer solstice (21 June north of the equator, and 21 December south of the equator). Weather conditions at the plant location (Zaragoza, Spain) were taken from the Meteonorm database [61].
The appropriate match between the solar field area and the TES capacity was also taken into account. By keeping the TES capacity constant, increasing the solar field area enabled a higher annual solar thermal energy production, but this effect was gradually reduced because, with an oversized solar field, a greater amount of harnessed solar energy could not be stored and was not used. Furthermore, for the same solar field area, increasing the TES capacity promoted a rise in the annual solar thermal energy production, but when the TES capacity was oversized the increase in the solar thermal production was reduced. A more detailed analysis is provided in [53,54].

Biomass Boiler
The hybridization of the solar system with a biomass boiler is only considered in the case of system A. The biomass boiler complements the solar field production whenever the solar resources are not able to supply the thermal power required by the ORC to run at full load, thus allowing for a more stable operation of the power production system. Therefore, the use of the biomass boiler, which is installed in parallel with the solar system, as shown in Figure 5, is only allowed during the operating hours of the solar system. The considered biomass boiler is a moving grate furnace with a nominal thermal power of 2.8 MW t and energy efficiency (LHV basis) of 0.82, fed with commercial biomass pellets [57].

Organic Rankine Cycle
The purpose of the ORC is to produce electricity from the heat supplied by the solar system hybridized with the biomass boiler, in the case of system A, and by the solar system alone, in the case of system B. Covering the annual electricity consumption of the commercial center was the design criterion of both systems. Therefore, following the procedure described in [53,54], and considering the technical features and data of Turboden 6/7 HR and Turboden 10/14 HR [50], the selected ORCs were: (i) 500 kW e of nominal electrical power, with a nominal energy efficiency of 18% in the case of system A, in which the ORC only operates at full load; and (ii) 978 kW e of nominal electrical power, with the same nominal energy efficiency of 18% in the case of system B, in which the ORC is allowed to operate at partial load. As a result, the thermal power input required by each ORC at full load operation corresponds to 2.78 MW t and 5.43 MW t for systems A and B, respectively. The inlet and outlet temperatures of the HTF that feeds each ORC are 214 • C and 121 • C (Table 1), respectively. These temperature values are in agreement with the technical data of selected ORCs [50].
The selected ORCs are single-pressure cycles with regenerative preheating of the organic fluid [50]. The appropriate integration of each ORC with the heat source in all the analyzed cases was verified through a thermodynamic model of the Rankine Cycle, built with the Engineering Equation Solver (EES) software [62], considering R245fa as the working fluid due to its favorable behavior [63]. Taking into account the technical characteristics and operating conditions of the ORCs considered herein, the thermodynamic cycles were replicated, assuming reasonable values for the technical performance parameters (e.g., isentropic efficiencies in turbine, pump, regenerative heat exchanger effectiveness, minimum temperature differences in evaporator and condenser) of the different pieces of equipment of the ORCs [28]. The results were in agreement with commercial equipment and performance data provided by Turboden [50]. Based on Turboden technical data, the ORC can provide part-load operation down to 10% of nominal load and maintain 90% of the cycle efficiency down to a 50% load, as shown in Figure 8. It was confirmed that, even at the partial load operation of the ORC integrated with the solar system in the case of system B, fundamental physical laws, particularly the Second Law of thermodynamics, were fulfilled. of equipment of the ORCs [28]. The results were in agreement with commercial equipment and performance data provided by Turboden [50]. Based on Turboden technical data, the ORC can provide part-load operation down to 10% of nominal load and maintain 90% of the cycle efficiency down to a 50% load, as shown in Figure 8. It was confirmed that, even at the partial load operation of the ORC integrated with the solar system in the case of system B, fundamental physical laws, particularly the Second Law of thermodynamics, were fulfilled.

Mechanical Chiller
Vapor compression refrigeration chillers (mechanical chillers), based on a reverse Rankine Cycle, cover the annual thermal cooling demand (6455 MWht) of the considered commercial center, consuming electricity produced by the ORC and/or electricity purchased from the electric grid, and producing chilled water at 7 °C. The commercial chiller selected was COBALT W 153.3 [58], whose technical data are shown in Table 1. Monthly and daily electrical demand profiles for cooling throughout the year have been evaluated and their shapes are similar to the thermal cooling demand profiles shown in Figure 3, but corrected with the chiller's Energy Efficiency of Refrigeration (EER = 4.68). In order to guarantee the cooling demand, three mechanical chillers were considered to be installed in both systems A and B, with a total nominal cooling thermal power capacity of 4581 kWt and 978 kWe nominal electrical power.

Operation of the Analyzed Systems
This section describes the annual, monthly, and hourly operation of systems A and B. It is worth keeping in mind the following: (i) there is a cooling demand for air conditioning throughout the year (see Figures 1 and 3); (ii) the ORC is in operation only during the operating hours of the solar

Mechanical Chiller
Vapor compression refrigeration chillers (mechanical chillers), based on a reverse Rankine Cycle, cover the annual thermal cooling demand (6455 MWh t ) of the considered commercial center, consuming electricity produced by the ORC and/or electricity purchased from the electric grid, and producing chilled water at 7 • C. The commercial chiller selected was COBALT W 153.3 [58], whose technical data are shown in Table 1. Monthly and daily electrical demand profiles for cooling throughout the year have been evaluated and their shapes are similar to the thermal cooling demand profiles shown in Figure 3, but corrected with the chiller's Energy Efficiency of Refrigeration (EER = 4.68). In order to guarantee the cooling demand, three mechanical chillers were considered to be installed in both systems A and B, with a total nominal cooling thermal power capacity of 4581 kW t and 978 kW e nominal electrical power.

Operation of the Analyzed Systems
This section describes the annual, monthly, and hourly operation of systems A and B. It is worth keeping in mind the following: (i) there is a cooling demand for air conditioning throughout the year (see Figures 1 and 3); (ii) the ORC is in operation only during the operating hours of the solar system-for the rest of the time, the electricity required must be purchased from the grid; (iii) in system A, the ORC is operating at full load with the support of the biomass boiler, which complements the thermal power produced by the solar system; (iv) in system B, there is not an auxiliary biomass boiler and the ORC can modulate operating at partial load according to the thermal power supplied by the solar system; and (v) electricity purchase and sale are allowed when there is deficit or excess of production, respectively.

System A: Hybrid Solar Parabolic Trough-ORC-Biomass System
It has been considered that system A starts its operation when the harnessed solar power is higher than 0.5 MW t . Therefore, the total annual hours of operation of the system reach 5117 h/year. The annual energy balance of the whole system is shown in Figure 9.
On an annual basis, the solar system (#2) produces 72% of the thermal energy (#4) and electrical energy (#5) of the cogeneration plant, while the biomass boiler contributes the remaining 28% (#3). The total electricity produced by the ORC (#5) amounts to 2560 MWh e . This would be enough to cover almost all (93%) of the electrical requirements (#8) for cooling (#6) and for electricity demands (#7). However, due to the irregular availability of solar resources, even with the support of the thermal energy storage and the biomass boiler, only 73% of the electricity produced (#5) is actually consumed by the system (#9), while the remaining 27% is sold to the grid (#11), and the system must purchase 867 MWh e of electricity from the grid (#10). The obtained renewable (solar + biomass) and solar fractions in the electricity consumed by the system (#8) are about 68% (#9/#8) and 49% (#2/#4·#9/#8), respectively. Noteworthily, the electricity required for cooling (#6) is virtually the same as the electricity demand of the commercial center (#7).
complements the thermal power produced by the solar system; (iv) in system B, there is not an auxiliary biomass boiler and the ORC can modulate operating at partial load according to the thermal power supplied by the solar system; and (v) electricity purchase and sale are allowed when there is deficit or excess of production, respectively.

System A: Hybrid Solar Parabolic Trough-ORC-Biomass System
It has been considered that system A starts its operation when the harnessed solar power is higher than 0.5 MWt. Therefore, the total annual hours of operation of the system reach 5117 h/year. The annual energy balance of the whole system is shown in Figure 9.  On an annual basis, the solar system (#2) produces 72% of the thermal energy (#4) and electrical energy (#5) of the cogeneration plant, while the biomass boiler contributes the remaining 28% (#3). The total electricity produced by the ORC (#5) amounts to 2560 MWhe. This would be enough to cover almost all (93%) of the electrical requirements (#8) for cooling (#6) and for electricity demands (#7). However, due to the irregular availability of solar resources, even with the support of the thermal energy storage and the biomass boiler, only 73% of the electricity produced (#5) is actually consumed by the system (#9), while the remaining 27% is sold to the grid (#11), and the system must purchase 867 MWhe of electricity from the grid (#10). The obtained renewable (solar + biomass) and solar fractions in the electricity consumed by the system (#8) are about 68% (#9/#8) and 49% (#2/#4·#9/#8), respectively. Noteworthily, the electricity required for cooling (#6) is virtually the same as the electricity demand of the commercial center (#7). Table 2 shows the monthly energy balance of system A. Due to the higher availability of solar resources in the summer, with more direct sunlight hours than in the winter: (i) the contribution of the solar system (#2) increases in the summer, reaching about 93% of the thermal production (#4) in July, compared to only 37% in December; (ii) the solar system operates for more hours in the summer, so that electricity production in the ORC also increases, being almost three times higher in July than in December. These effects are explained further for the hourly operation of the system.  Table 2 shows the monthly energy balance of system A. Due to the higher availability of solar resources in the summer, with more direct sunlight hours than in the winter: (i) the contribution of the solar system (#2) increases in the summer, reaching about 93% of the thermal production (#4) in July, compared to only 37% in December; (ii) the solar system operates for more hours in the summer, so that electricity production in the ORC also increases, being almost three times higher in July than in December. These effects are explained further for the hourly operation of the system.  1  651  847  1630  1943  2493  2666  3200  2567  1955  1117  767  567  2  280  423  863  1011  1272  1314  1519  1331  1019  596  352  226  3  410  356  430  407  366  271  119  134  399  438  315  377  4  689  778  1292  1417  1637  1584  1637  1465  1417  1034  667  603  5  124  140  232  255  294  285  294  263  255  186  120  109  6  44  42  53  80  182  229  225  199  129  85  54  46  7  116  105  116  113  116  113  116  116  113  116  113  116  8  161  147  169  193  299  341  342  315  242  201  167  162  9  83  97  147  171  217  216  222  212  194  151  88  72  10  78  50  22  21  82  125  119  104  48  50  79  90  11  41  43  85  84  78  69  72  52  61  35  32  36  12  208  196  249  374  854  1070  1055  932  606  397  254  214  13  500  434  524  496  446  331  145  164  486  535  385  460 The dynamic hourly behavior throughout the year of the solar system (solar field integrated with TES) is obtained with SAM [49]. Figure 10 shows the hourly profiles on a typical day in selected months of the year: (i) the total incident thermal power (MW t ) to the solar field considering the deviation with respect to the normal solar irradiance; (ii) the thermal power (MW t ) harnessed by the HTF in the receiver; (iii) the thermal power (MW t ) released to the ORC (heat sink); (iv) the thermal power (MW t ) discharged from the TES; and (v) the thermal power (MW t ) charged to the TES. Note that hourly values of these magnitudes for the whole year have been calculated, obtaining detailed information about the hourly available solar resources and the hourly solar thermal energy released to the ORC, as well as detailed information about each piece of equipment in the solar system, including the hourly behavior of the TES throughout the year. The monthly thermal energy provided by the solar system and biomass boiler is shown on the left side of Figure 11. On the right side are the hourly thermal productions during a day in January (winter) and a day in July (summer). As previously explained, the operation of the solar system starts after the thermal output reaches the threshold of 0.5 MWt, so that the operating hours of the solar system vary significantly throughout the year according to the solar resource availability, as in January (8 h/day) and in July (19 h/day). Moreover, since the biomass boiler only works during the operating hours of the solar system, to complement the thermal output necessary to reach the ORC's required thermal input nominal power: (i) the sum of the solar and biomass thermal energy in each hour is always equal to 2.78 MWt; and (ii) the hourly contribution of the biomass boiler is higher in the winter months, when there are less available solar resources, as in January the biomass boiler covers about 60% of the daily thermal production, with a well-distributed load throughout the day, while in July it covers only 7%, being much more prominent at hours 01 and 07. The monthly thermal energy provided by the solar system and biomass boiler is shown on the left side of Figure 11. On the right side are the hourly thermal productions during a day in January (winter) and a day in July (summer). As previously explained, the operation of the solar system starts after the thermal output reaches the threshold of 0.5 MW t , so that the operating hours of the solar system vary significantly throughout the year according to the solar resource availability, as in January (8 h/day) and in July (19 h/day). Moreover, since the biomass boiler only works during the operating hours of the solar system, to complement the thermal output necessary to reach the ORC's required thermal input nominal power: (i) the sum of the solar and biomass thermal energy in each hour is always equal to 2.78 MW t ; and (ii) the hourly contribution of the biomass boiler is higher in the winter months, when there are less available solar resources, as in January the biomass boiler covers about 60% of the daily thermal production, with a well-distributed load throughout the day, while in July it covers only 7%, being much more prominent at hours 01 and 07.
While the thermal energy production of the solar system and biomass boiler is the same on each day of the month, the electricity requirements (electricity demand and electricity for cooling demands) of the commercial center presented two different profiles, namely working days and holidays, as explained in Section 2. Figure 12 shows the distribution of the electricity produced by the ORC, bought from and sold to the electric grid on working days; on the left side, the monthly distribution is shown and on the right side the hourly values for January (winter) and July (summer) are shown. As can be seen, electricity production in the ORC follows the thermal energy availability, being higher in the summer months and lower in winter. Regarding the hourly profiles, for a working day in January, the system generates surplus electricity, which is sold to the grid from hours 10 to 17, and there is no need for electricity purchase (during operating hours of the solar system). On the other hand, for a working day in July, the high electricity requirements lead the system to purchase electricity from hours 11 to 20. While the thermal energy production of the solar system and biomass boiler is the same on each day of the month, the electricity requirements (electricity demand and electricity for cooling demands) of the commercial center presented two different profiles, namely working days and holidays, as explained in Section 2. Figure 12 shows the distribution of the electricity produced by the ORC, bought from and sold to the electric grid on working days; on the left side, the monthly distribution is shown and on the right side the hourly values for January (winter) and July (summer) are shown. As can be seen, electricity production in the ORC follows the thermal energy availability, being higher in the summer months and lower in winter. Regarding the hourly profiles, for a working day in January, the system generates surplus electricity, which is sold to the grid from hours 10 to 17, and there is no need for electricity purchase (during operating hours of the solar system). On the other hand, for a working day in July, the high electricity requirements lead the system to purchase electricity from hours 11 to 20. Similarly, Figure 13 presents the monthly and hourly distributions of the electricity produced by the ORC, bought from and sold to the electric grid on holidays by system A. It should be noted that the monthly values shown in Figure 13 are lower than the ones shown in Figure 12 because the number of holidays per month is about half of the number of working days. In fact, as explained in Section 2, taken individually, holidays presented higher electricity requirements (electricity demand and electricity for cooling demand) than working days. On a holiday in January (winter), while the  While the thermal energy production of the solar system and biomass boiler is the same on each day of the month, the electricity requirements (electricity demand and electricity for cooling demands) of the commercial center presented two different profiles, namely working days and holidays, as explained in Section 2. Figure 12 shows the distribution of the electricity produced by the ORC, bought from and sold to the electric grid on working days; on the left side, the monthly distribution is shown and on the right side the hourly values for January (winter) and July (summer) are shown. As can be seen, electricity production in the ORC follows the thermal energy availability, being higher in the summer months and lower in winter. Regarding the hourly profiles, for a working day in January, the system generates surplus electricity, which is sold to the grid from hours 10 to 17, and there is no need for electricity purchase (during operating hours of the solar system). On the other hand, for a working day in July, the high electricity requirements lead the system to purchase electricity from hours 11 to 20. Similarly, Figure 13 presents the monthly and hourly distributions of the electricity produced by the ORC, bought from and sold to the electric grid on holidays by system A. It should be noted that the monthly values shown in Figure 13 are lower than the ones shown in Figure 12 because the number of holidays per month is about half of the number of working days. In fact, as explained in Section 2, taken individually, holidays presented higher electricity requirements (electricity demand and electricity for cooling demand) than working days. On a holiday in January (winter), while the Similarly, Figure 13 presents the monthly and hourly distributions of the electricity produced by the ORC, bought from and sold to the electric grid on holidays by system A. It should be noted that the monthly values shown in Figure 13 are lower than the ones shown in Figure 12 because the number of holidays per month is about half of the number of working days. In fact, as explained in Section 2, taken individually, holidays presented higher electricity requirements (electricity demand and electricity for cooling demand) than working days. On a holiday in January (winter), while the system is in operation, it is not necessary to purchase electricity from the grid, and even a small part of the electricity produced is sold. However, on a holiday in July (summer), the higher electricity requirements lead the system to purchase electricity from the grid from hours 10 to 20, limiting the electricity sale.
Energies 2020, 13, x FOR PEER REVIEW 18 of 31 system is in operation, it is not necessary to purchase electricity from the grid, and even a small part of the electricity produced is sold. However, on a holiday in July (summer), the higher electricity requirements lead the system to purchase electricity from the grid from hours 10 to 20, limiting the electricity sale.

System B: Parabolic Trough-ORC System
In the previous section, the ORC of the proposed hybridized system was always operating at full load (500 kWe), thanks to the support of the biomass auxiliary boiler, which complemented the thermal power produced during the operating hours of the solar collectors. This section analyzes the behavior of the cogeneration plant without the support of the biomass boiler and allowing partial load operation of the ORC, according to the technical data of the manufacturer of the selected ORC shown in Figure 8 [50]. A simplified diagram of system B is shown in Figure 14. The main components (i.e., solar system and ORC) have been resized following the procedure explained in [53,54] to fulfil the criterion of supplying a similar amount of annual electrical energy to that of a conventional system based on electricity purchased from the grid of the commercial center. The technical specifications of system B are shown in Table 1. The annual energy balance of system B for a year of operation is shown in Figure 14. In this case, the solar system produces 2979 MWhe (#5), which is slightly higher than the annual electrical demand, 2740 MWhe (#8), of the commercial center for cooling (#6) and for other electrical demands (#7). However, due to the irregular availability of solar resources, even with the support of thermal energy storage, only about 69% of the electricity produced by the ORC is actually consumed (#9), while the

System B: Parabolic Trough-ORC System
In the previous section, the ORC of the proposed hybridized system was always operating at full load (500 kW e ), thanks to the support of the biomass auxiliary boiler, which complemented the thermal power produced during the operating hours of the solar collectors. This section analyzes the behavior of the cogeneration plant without the support of the biomass boiler and allowing partial load operation of the ORC, according to the technical data of the manufacturer of the selected ORC shown in Figure 8 [50]. A simplified diagram of system B is shown in Figure 14. The main components (i.e., solar system and ORC) have been resized following the procedure explained in [53,54] to fulfil the criterion of supplying a similar amount of annual electrical energy to that of a conventional system based on electricity purchased from the grid of the commercial center. The technical specifications of system B are shown in Table 1.
Energies 2020, 13, x FOR PEER REVIEW 18 of 31 system is in operation, it is not necessary to purchase electricity from the grid, and even a small part of the electricity produced is sold. However, on a holiday in July (summer), the higher electricity requirements lead the system to purchase electricity from the grid from hours 10 to 20, limiting the electricity sale.

System B: Parabolic Trough-ORC System
In the previous section, the ORC of the proposed hybridized system was always operating at full load (500 kWe), thanks to the support of the biomass auxiliary boiler, which complemented the thermal power produced during the operating hours of the solar collectors. This section analyzes the behavior of the cogeneration plant without the support of the biomass boiler and allowing partial load operation of the ORC, according to the technical data of the manufacturer of the selected ORC shown in Figure 8 [50]. A simplified diagram of system B is shown in Figure 14. The main components (i.e., solar system and ORC) have been resized following the procedure explained in [53,54] to fulfil the criterion of supplying a similar amount of annual electrical energy to that of a conventional system based on electricity purchased from the grid of the commercial center. The technical specifications of system B are shown in Table 1. The annual energy balance of system B for a year of operation is shown in Figure 14. In this case, the solar system produces 2979 MWhe (#5), which is slightly higher than the annual electrical demand, 2740 MWhe (#8), of the commercial center for cooling (#6) and for other electrical demands (#7). However, due to the irregular availability of solar resources, even with the support of thermal energy storage, only about 69% of the electricity produced by the ORC is actually consumed (#9), while the The annual energy balance of system B for a year of operation is shown in Figure 14. In this case, the solar system produces 2979 MWh e (#5), which is slightly higher than the annual electrical demand, 2740 MWh e (#8), of the commercial center for cooling (#6) and for other electrical demands (#7). However, due to the irregular availability of solar resources, even with the support of thermal energy storage, only about 69% of the electricity produced by the ORC is actually consumed (#9), while the remaining 31% is sold to the electric grid (#11), and the system is required to purchase 690 MWh e from the electric grid (#10). The resulting solar fraction in the electricity consumed by the system (#8) is thus approximately 75% (#9/#8). Table 3 shows the monthly energy balance of system B. This shows how the solar system is able to produce (#5) the electricity required for cooling (#6) during all the months of the year, except in December, and is able to produce more electricity (#5) than demanded (#8) for seven months, from March to September. Nevertheless, it should be noted that even in summer months, when the electricity production is higher because there is a higher availability of solar resources, it is also necessary to purchase electricity in order to fulfil the electrical demands of the commercial center during the night, e.g., in July, the electricity produced (#5) is close to 150% of the electricity demanded (#8), but the system must still purchase about 26 MWh e from the grid (#10).  Table 4 shows the hourly values of thermal power provided by the solar system, obtained from the dynamic simulation of the solar system in SAM [49]. These values represent the inputted thermal energy of the ORC for the production of electricity. According to the partial load operation chart of the ORC shown in Figure 8, the minimum thermal input required for the operation of the ORC is about 1 MW t and the maximum at full load is 5.43 MW t . Therefore, the ORC could be in operation for about 4874 h per year. The nominal electrical efficiency value of the ORC at full load is 18% (see Table 1), and during the annual operation of system B, the annual average electrical efficiency of the ORC reaches 16.7% (Figure 14). According to Figure 8, the energy efficiency of the ORC remains at its maximum value (18%) for operating loads between 80% and 100%. That is, when the thermal input is higher than 4 MW t , the energy efficiency of the ORC is maximum, and for a thermal input over 3 MW t (50% load), the energy efficiency of the ORC is about 16%. The values presented in Table 4 show that the ORC was operating at loads higher than 50% for most of the year (8 months), being especially high in the summer months (May to September, higher than 80% and reaching full load most of the time), while falling below 50% in the winter (November to February).
The monthly values of thermal energy production of the solar system of system B are shown on the left side of Figure 15. As expected, in summer months, the thermal energy production is five times higher than in winter months, which is a wider variation than in the case of system A, supported with a biomass boiler. The hourly thermal productions of the hybrid system on a day in January and in July are shown on the right side of Figure 15. In January (winter), the system is in operation for 8 h/day, during which the thermal power provided to the ORC is always lower than 3 MW t . This means that the ORC operates below a 50% load at all the operation times, even being at its technical limit of a 10% load at hour 17, which is an unfavorable operation condition in terms of the technical feasibility and reliability of the system. In July (summer), the system is in operation for 19 h/day, and its operation is close to 100% load for about 12 h. The rest of the time, the system is operating above an 80% load and modulates its operation between 30 and 50% loads at the times at which the system starts up and shuts down. The monthly values of thermal energy production of the solar system of system B are shown on the left side of Figure 15. As expected, in summer months, the thermal energy production is five times higher than in winter months, which is a wider variation than in the case of system A, supported with a biomass boiler. The hourly thermal productions of the hybrid system on a day in January and in July are shown on the right side of Figure 15. In January (winter), the system is in operation for 8 h/day, during which the thermal power provided to the ORC is always lower than 3 MWt. This means that the ORC operates below a 50% load at all the operation times, even being at its technical limit of a 10% load at hour 17, which is an unfavorable operation condition in terms of the technical feasibility and reliability of the system. In July (summer), the system is in operation for 19 h/day, and its operation is close to 100% load for about 12 h. The rest of the time, the system is operating above an 80% load and modulates its operation between 30 and 50% loads at the times at which the system starts up and shuts down. For system B, the monthly electricity produced, bought and sold is presented for working days in Figure 16 and for holidays in Figure 17. As can be seen, electricity purchases on working days are more restricted to colder months (from October to February), being particularly low during the rest of the year, when electricity sale is more pronounced. On the other hand, on holidays the system must purchase electricity all year round, except for the months of March and April, when it is able to sell relatively less electricity to the grid. As explained for system A, it should be noted that the For system B, the monthly electricity produced, bought and sold is presented for working days in Figure 16 and for holidays in Figure 17. As can be seen, electricity purchases on working days are more restricted to colder months (from October to February), being particularly low during the rest of the year, when electricity sale is more pronounced. On the other hand, on holidays the system must purchase electricity all year round, except for the months of March and April, when it is able to sell relatively less electricity to the grid. As explained for system A, it should be noted that the monthly values shown in Figure 17 are lower than the ones shown in Figure 16 because the number of holidays per month is about half of the number of working days.  Figure 17 are lower than the ones shown in Figure 16 because the number of holidays per month is about half of the number of working days.  The hourly electricity produced, bought from and sold to the electric grid on working days and holidays in January (winter) and July (summer) is shown on the right side of Figures 16 and 17. It is shown that, on a working day in January, the system is able to cover virtually all of the electricity requirements of the commercial center for almost 4 h, with it even being possible to sell a very reduced quantity of electricity at hours 11 and 15. In contrast, on a holiday in January, the system only reaches electricity self-sufficiency for about 2 h, and no electricity sale is possible. In July, when the electricity requirements (electricity for cooling and electricity demand) are considerably higher, the electricity produced on a working day is sufficiently high to cover the electricity demand and sell to the grid during all operating hours. In the case of a holiday in July, electricity purchase is required for 8 h.

Economic and Environmental Analysis
Once the technical feasibility of the proposed systems (A and B) has been proven, an estimate of their economic cost is presented, as well as of their potential environmental benefit with respect to a conventional system in which all the electrical requirements are covered with electricity purchased  monthly values shown in Figure 17 are lower than the ones shown in Figure 16 because the number of holidays per month is about half of the number of working days.  The hourly electricity produced, bought from and sold to the electric grid on working days and holidays in January (winter) and July (summer) is shown on the right side of Figures 16 and 17. It is shown that, on a working day in January, the system is able to cover virtually all of the electricity requirements of the commercial center for almost 4 h, with it even being possible to sell a very reduced quantity of electricity at hours 11 and 15. In contrast, on a holiday in January, the system only reaches electricity self-sufficiency for about 2 h, and no electricity sale is possible. In July, when the electricity requirements (electricity for cooling and electricity demand) are considerably higher, the electricity produced on a working day is sufficiently high to cover the electricity demand and sell to the grid during all operating hours. In the case of a holiday in July, electricity purchase is required for 8 h.

Economic and Environmental Analysis
Once the technical feasibility of the proposed systems (A and B) has been proven, an estimate of their economic cost is presented, as well as of their potential environmental benefit with respect to a conventional system in which all the electrical requirements are covered with electricity purchased The hourly electricity produced, bought from and sold to the electric grid on working days and holidays in January (winter) and July (summer) is shown on the right side of Figures 16 and 17. It is shown that, on a working day in January, the system is able to cover virtually all of the electricity requirements of the commercial center for almost 4 h, with it even being possible to sell a very reduced quantity of electricity at hours 11 and 15. In contrast, on a holiday in January, the system only reaches electricity self-sufficiency for about 2 h, and no electricity sale is possible. In July, when the electricity requirements (electricity for cooling and electricity demand) are considerably higher, the electricity produced on a working day is sufficiently high to cover the electricity demand and sell to the grid during all operating hours. In the case of a holiday in July, electricity purchase is required for 8 h.

Economic and Environmental Analysis
Once the technical feasibility of the proposed systems (A and B) has been proven, an estimate of their economic cost is presented, as well as of their potential environmental benefit with respect to a conventional system in which all the electrical requirements are covered with electricity purchased from the grid, including the electricity required to produce cooling in the mechanical chillers. A summary of these analyses is presented in Table 5.

Economic Analysis
With respect to the economic analysis, a rough estimate is presented, considering 20 years of useful life and a linear amortization rate of 0.05 year −1 . Operation and maintenance costs and other indirect costs have not been considered in the analysis. The investment costs of the various pieces of equipment shown in Table 5 have been obtained from scientific specialized literature. The price of the electricity purchased from the grid corresponds to the average value of the electricity price of the year 2019 in Spain for industrial consumption [68]. The price of the electricity sold has been estimated at about 50% of the purchase price. The price of biomass pellets has been taken from the technical documents of the Spanish Government, by collecting the evolution of the price of biomass pellets for thermal use in Spain [69].
The total investment cost of the conventional system, consisting of the mechanical chillers, was EUR 824,580 and its total annual cost was estimated at 414,691 EUR/year, consisting of a chiller annual cost = 41,229 EUR/year and an electricity cost = 373,462 EUR/year. Relative to the conventional system, the total annual costs of system A (hybridized with biomass) 632,350 EUR/year and system B (non-hybridized solar cogeneration), 506,135 EUR/year, are 52% and 22% higher, respectively. Therefore, neither of the analyzed systems is presently profitable.
In the case of system A, the additional investment cost required with respect to the conventional system (6,865,880 − 824,580 = EUR 6,041,300) is not recovered with the annual energy cost savings in the energy bill (373,462 − 289,056 = 84,406 EUR/year) during the plant's lifetime of 20 years. In fact, the simple payback period (annual energy cost savings divided by the additional investment cost) obtained is higher than 70 years and the corresponding net present value NPV A = −5,177,760 EUR is a negative value. The cost of the electricity produced by the solar system hybridized with a biomass boiler has been estimated by dividing the sum of the annual investment cost of the solar field, TES, biomass boiler, and ORC, plus the annual energy cost of the biomass pellets (519,803 EUR/year), by the annual electricity produced in the ORC (2560 MWh e /year), obtaining the value of 0.2030 EUR/kWh, which is about 1.5 times greater than the cost of the electricity purchased from the grid (0.1363 EUR/kWh).
In the case of system B, which does not include the biomass boiler, and the ORC is fully driven by solar thermal energy operating at a partial load according to the solar resource availability; the installed capacities of the solar field and ORC are higher than those of system A in order to guarantee a similar annual electrical energy production (electricity for cooling and electricity demand). As a result, the investment cost of system B amounts to EUR 9,510,280, which is significantly higher than the investment cost of system A (EUR 6,865,880) because of the larger installed capacities of the solar system and ORC, even without the biomass boiler. On the other hand, the resulting total annual cost of system B (506,135 EUR/year) is actually lower than that of system A (632,350EUR/year), due to important total annual energy cost savings (289,056 − 30,621 = 258,435 EUR/year) from electricity sale revenues and the non-consumption of biomass pellets. Relative to the conventional system, the total annual energy cost savings (373,462 − 30,621 = 342,841 EUR/year) are still not enough to recover the additional investment cost (9,510,280 − 824,580 = EUR 8,685,700) during the 20-year plant lifetime. The resulting simple payback period is about 25 years and the corresponding net present value, NPV B = −2,653,460 EUR, is also a negative value. For system B, the cost of the electricity produced by the solar system is 0.1458 EUR/kWh, which is only 7% higher than the electricity purchase price (0.1363 EUR/kWh).

Environmental Analysis
The environmental analysis (Table 5) was performed in terms of the total annual CO 2 emissions and total annual non-renewable primary energy consumption. It should be noted that only the environmental burden associated to the operation of the system has been considered (i.e., environmental loads associated to the electricity from the grid and the pellets consumed), while the materials of the pieces of equipment and the construction of the installation were not taken into account. Regarding the conventional system, its total annual CO 2 emissions and total annual non-renewable primary energy consumption were estimated at 906,940 kgCO 2 /year and 5,353,960 kWh/year, respectively. As can be seen from Table 5, the total annual CO 2 emissions generated by systems A and B have been significantly reduced with respect to the conventional system.
In system A, the reduction in CO 2 emissions and non-renewable primary energy consumption reached about 85%. The CO 2 emissions associated with the electricity produced in system A were estimated at 0.034 kgCO 2 /kWh e , corresponding to the CO 2 emissions generated from biomass burning (88,272 kgCO 2 /year) divided by the annual electricity produced in the ORC (2560 MWh e /year). This value is one order of magnitude lower than that of the electricity from the Spanish grid (0.331 kgCO 2 /kWh e ).
In system B, the annual CO 2 emissions and non-renewable primary energy consumption take on negative values. This means that the proposed system is not only carbon neutral, but also offsets part of the greenhouse gas emissions and non-renewable primary energy consumption of the electrical network. This is obtained thanks to the renewable electricity produced in the solar system and sold to the electric grid, which displaces grid electricity, avoiding the emission of 0.331 kgCO 2 /kWh e and the consumption of 1.954 kWh/kWh e of non-renewable primary energy. In contrast to system A, in which the system generates some CO 2 emissions due to biomass consumption, system B only consumes solar energy and, therefore, the CO 2 emissions associated with the electricity produced by the system are zero. Note that, in this analysis, only environmental loads associated to the operation of the system (i.e., electricity consumed from the grid and biomass pellets) were considered.

Sensitivity Analysis
Even though the economic analysis performed in Section 5.1 showed that systems A and B are presently not feasible from an economic viewpoint, the costs of the produced electricity are relatively low. Thus, in this section, we will analyze the factors influencing the costs of both systems.
As can be seen in Table 5, for system A, the investment cost of the solar field of parabolic trough collectors (153,690 EUR/year) represents almost 45% of the total annual investment cost (343,294 EUR/year) and 24% of the total annual cost (632,350 EUR/year). Regarding the operation costs, the energy cost of biomass is also quite high (217,738 EUR/year) and represents about 30% of the total annual cost. In the case of system B, the investment cost of the solar field of parabolic trough collectors (230,535 EUR/year) represents more than 45% of both the total annual investment costs and total annual costs. Note that, in system B, the total annual energy costs are very low, first, because there is no biomass consumption, and, second, because the electricity sale offsets about two thirds of the electricity purchase costs.
Since the investment costs (systems A and B) and the cost of biomass (system A only) accounted for significant shares of the total annual costs, we decided to assess the influence of the unit investment cost and the biomass purchase price on the cost of the electricity produced by the solar systems. The results presented in Figure 18 were obtained by modifying said costs independently and/or simultaneously. The line "Investment cost reduction" in system A was obtained by reducing the unit investment cost of the solar collectors (235 EUR/m 2 ), the TES (65 EUR/kWh t ), the biomass boiler (315 EUR/kW t ) and the ORC (2000 EUR/kW e ) in proportion to the cost factor (a cost factor equal to one means that the cost remains unchanged, and a cost factor equal to 0.8 means that the cost of the considered component has experienced a 20% cost reduction, and so on), while keeping unchanged the cost of the rest of the plant components (chillers), energy resources and considered parameters (annual amortization factor and plant lifetime). Similarly, the line "Biomass price reduction" was obtained by reducing only the price of the biomass (0.0444 EUR/kWh). In contrast, the line "Both biomass and investment cost reduction" was obtained by simultaneously reducing the unit investment cost of the solar collectors, TES, biomass boiler, ORC and biomass price by applying the same cost factor. Clearly, only the effect of the unit investment cost of equipment (i.e., solar collectors, TES and ORC) was analyzed for system B, since there is no biomass boiler and no consumption of pellets. For system A, the results show the significant influence of the biomass price. Only producing electricity with biomass is expensive due to the low energy efficiency of the ORC (about 20%), which means that the cost of each kWh of electricity is more than six times higher than the cost of each kWh of biomass consumed in the boiler. Note that the cost of the electricity produced is higher than the price of pellets divided by the efficiencies of the biomass boiler and the ORC. For this reason, reducing the biomass purchase cost has a significant effect on the cost of the electricity produced in system A. Nevertheless, a notable reduction in the biomass price (about 80%) is required to reach the electricity cost of the electricity produced in system A, which is similar to the price of electricity (0.1363 EUR/kWh) purchased from the grid. Reducing the investment cost of equipment (concentrating on the parabolic trough solar collectors, TES, biomass boiler and ORC) has a higher influence. Thus, in system A, it is possible to reach an electricity production cost lower than the electricity purchase price by reducing about 60% the investment cost. Combining both factors, with a reduction of 35% in the investment cost of equipment (concentrating parabolic trough solar collectors, TES, biomass boiler and ORC) and in the price of biomass, the cost of electricity produced would be lower than the price of electricity purchased.
For system B, in which the electricity cost is close to the electricity purchase price, with less than a 10% reduction in the investment cost of equipment (concentrating parabolic trough solar collectors, TES and ORC), it would already be possible to produce electricity at a lower cost than the electricity purchased from the grid.
Taking into account only the main parameters with the biggest influences, the obtained results show that, if properly optimized, a solar-based parabolic trough integrated with TES and ORC cogeneration plants, with possible hybridization using biomass, could become competitive in the medium term, provided that the significant cost reduction trend of this type of technology is maintained, which is expected for the next few years [20].

Conclusions
The technical feasibility of two novel solar-based PTC-ORC cogeneration systems, producing power and cooling only, and designed to cover the annual energy demands of a commercial center located in Zaragoza (Spain), has been proven. The commercial center considered herein, despite its high energy consumption, had a very low or null heat demand, which is a common situation in some kinds of buildings, such as commercial centers, depending, among other things, on their geographical location.
The analyzed cogeneration systems consist of solar parabolic trough collectors integrated with thermal energy storage, an Organic Rankine Cycle, and vapor compression mechanical chillers. System A is hybridized with a biomass boiler that supports solar thermal production and allows full For system A, the results show the significant influence of the biomass price. Only producing electricity with biomass is expensive due to the low energy efficiency of the ORC (about 20%), which means that the cost of each kWh of electricity is more than six times higher than the cost of each kWh of biomass consumed in the boiler. Note that the cost of the electricity produced is higher than the price of pellets divided by the efficiencies of the biomass boiler and the ORC. For this reason, reducing the biomass purchase cost has a significant effect on the cost of the electricity produced in system A. Nevertheless, a notable reduction in the biomass price (about 80%) is required to reach the electricity cost of the electricity produced in system A, which is similar to the price of electricity (0.1363 EUR/kWh) purchased from the grid. Reducing the investment cost of equipment (concentrating on the parabolic trough solar collectors, TES, biomass boiler and ORC) has a higher influence. Thus, in system A, it is possible to reach an electricity production cost lower than the electricity purchase price by reducing about 60% the investment cost. Combining both factors, with a reduction of 35% in the investment cost of equipment (concentrating parabolic trough solar collectors, TES, biomass boiler and ORC) and in the price of biomass, the cost of electricity produced would be lower than the price of electricity purchased.
For system B, in which the electricity cost is close to the electricity purchase price, with less than a 10% reduction in the investment cost of equipment (concentrating parabolic trough solar collectors, TES and ORC), it would already be possible to produce electricity at a lower cost than the electricity purchased from the grid.
Taking into account only the main parameters with the biggest influences, the obtained results show that, if properly optimized, a solar-based parabolic trough integrated with TES and ORC cogeneration plants, with possible hybridization using biomass, could become competitive in the medium term, provided that the significant cost reduction trend of this type of technology is maintained, which is expected for the next few years [20].

Conclusions
The technical feasibility of two novel solar-based PTC-ORC cogeneration systems, producing power and cooling only, and designed to cover the annual energy demands of a commercial center located in Zaragoza (Spain), has been proven. The commercial center considered herein, despite its high energy consumption, had a very low or null heat demand, which is a common situation in some kinds of buildings, such as commercial centers, depending, among other things, on their geographical location. The analyzed cogeneration systems consist of solar parabolic trough collectors integrated with thermal energy storage, an Organic Rankine Cycle, and vapor compression mechanical chillers. System A is hybridized with a biomass boiler that supports solar thermal production and allows full load operation of the ORC throughout the solar system's operating time. System B, on the other hand, is not hybridized with biomass. For this reason, the installed capacities of the solar system and ORC are larger, and it is considered that the ORC can modulate, operating at partial load, in order to produce the electrical requirements of the commercial center (electricity for cooling production and electricity demand). Both systems were connected to the electrical grid, regarded as an auxiliary system, thereby selling the surplus electricity, and purchasing the required electricity when needed.
The present work was developed by applying a systematic procedure to obtain the preliminary sizing of the plant components considering commercially available pieces of equipment. Moreover, a detailed physical simulation of the operation of the system throughout the day and the year was also performed, capturing its hourly, daily, and monthly dynamic behavior, considering: (i) the variations in ambient temperature, solar radiation, and thermal energy production, among others; and (ii) the modelling of the different pieces of equipment based on technical data and the investment costs of commercially available devices.
The technical feasibility analysis was complemented by simplified economic and environmental analyses. From an economic viewpoint, the proposed solar-based cogeneration systems were not competitive with conventional systems, in which all required electricity is purchased from the grid. The configuration without biomass hybridization (system B) was more economically interesting than the one with a biomass boiler (system A), producing electricity at 0.1458 EUR/kWh, which is only 7% more expensive than the electricity purchase price in Spain. This was mainly due to the elimination of the purchase of biomass pellets, at the expense of the operational reliability provided by a dispatchable renewable energy source (i.e., biomass) and of a simpler operating strategy due to the partial load operation of the ORC in system B. All things considered, the results are promising and reveal the systems' potential interest and competitiveness in the short to medium term, particularly in a context in which external environmental costs could also be considered. From an environmental viewpoint, the proposed PTC-ORC systems present considerable benefits, contributing significantly to the decarbonization of energy production, and achieving a very high renewable fraction. Note that the plants' main energy resources were renewable energies (i.e., solar and biomass), and their auxiliary resource was grid electricity. Nevertheless, it would be advisable to develop this analysis by including the environmental burden associated with the equipment construction.
Future research efforts are required in order to develop optimized cogeneration and trigenerations systems that are mainly driven by renewable energies in terms of: (i) sizing; (ii) more appropriate energy integration schemes, considering the possibility of covering all energy services of a building (heating, cooling and electricity), implementing absorption chillers as well as thermal energy storage of the heating and cooling produced; and (iii) reducing the investment costs of these systems.