Thermoeconomic Modeling and Parametric Study of a Photovoltaic-Assisted 1 MW e Combined Cooling , Heating , and Power System

In this study a small-scale, completely autonomous combined cooling, heating, and power (CCHP) system is coupled to a photovoltaic (PV) subsystem, to investigate the possibility of reducing fuel consumption. The CCHP system generates electrical energy with the use of a simple gas turbine cycle, with a rated nominal power output of 1 MWe. The nominal power output of the PV subsystem is examined in a parametric study, ranging from 0 to 600 kWe, to investigate which configuration results in a minimum lifecycle cost (LCC) for a system lifetime of 20 years of service. The load profile considered is applied for a complex of households in Nicosia, Cyprus. The solar data for the PV subsystem are taken on an hourly basis for a whole year. The results suggest that apart from economic benefits, the proposed system also results in high efficiency and reduced CO2 emissions. The parametric study shows that the optimum PV capacity is 300 kWe. The minimum lifecycle cost for the PV-assisted CCHP system is found to be 3.509 million €, as compared to 3.577 million € for a system without a PV subsystem. The total cost for the PV subsystem is 547,445 €, while the total cost for operating the system (fuel) is 731,814 € (compared to 952,201 € for a CCHP system without PVs). Overall the proposed system generates a total energy output of 210,520 kWh (during its whole lifetime), which translates to a unit cost of 17 €/kWh.


Introduction
Cogeneration has gained more attention in recent years in an effort to develop highly efficient systems at all scales from the kW to the MW range [1][2][3][4].It is a method that has the potential of attaining high fuel efficiency, since the exhaust gases leave the system and enter the atmosphere with a very low exergy [5], which also suggests reduction in greenhouse gases (e.g.CO 2 emissions) [6].However, rigorous design considerations have to be made to achieve high system efficiency at minimum lifecycle cost (LCC).Cogeneration can have various applications, which can typically combine the generation of power, heating and cooling.In climates which require both space heating and space cooling in different periods of the year, a promising application of cogeneration is combined cooling, heating, and power (CCHP).Therefore, CCHP systems can be considered as alternatives to large-scale, electricity-only generating power plants, to improve fuel efficiency by recovering the rejected heat from the thermal cycle for district heating and district cooling.Although such a system can be applied in various scales, configurations and capacities, an interesting application is in small-scale, decentralized, completely autonomous systems.These distributed energy systems operate with minimum losses, due to their proximity with the serviced buildings.They typically range from 1 to 10 MW e and can provide power, space heating, domestic hot water and space cooling.Heating, or cooling, can be provided with a simple district energy network connecting the consumption site with the CCHP system.Since heating and cooling loads do not generally coincide, the same district network can be used in the summer to distribute cooling, produced by thermally activated absorption chillers, also located within the CCHP system [5,7].
Photovoltaic technology converts solar radiation directly into electricity and two commercialized types that are currently the most popular ones for practical applications are: (1) polycrystalline silicon photovoltaics (PVs), which is a well-established, tested, and mature technology; and (2) thin film solar cells, which is a relatively new technology, which can operate more efficiently than polycrystalline PVs (~10%) in low-light conditions (e.g., dawn, or cloudy day).In recent years, due to improvements in PV panel efficiency and manufacturing methods, the payback times have fallen to 2-3 years for crystalline silicon PV systems, and to almost 1 year for some thin-film ones, under moderate levels of sunshine [8].However the major problem is the high cost of manufacturing the sheets of semiconductor materials needed for power systems.The cost of PV modules in 2012 was 0.8 USD/W, but it was expected to drop to 0.6 USD/W (thin-film "First Solar") [9].The most important advantages of PV technology are: (a) they can be installed in a wide variety of places; (b) additional modules can be added incrementally; (c) quick installation; (d) little maintenance; (e) power can be produced during the afternoon hours when demand is at its peak; and (f) power generation is often on-site, thereby eliminating transmission losses [10].The economic feasibility of solar energy depends heavily on energy cost and available subsidies.Subsidies for PV systems in Cyprus are generally on the order of 50% of the total price (including purchase and installation costs).In addition, a further step could involve the integration of a hydrogen generating and storage subsystem.This modification could lead to the eventual development of a self-sustainable, hydrocarbon-free CCHP system.This can be achieved by further integration of an electrolyzer-hydrogen storage unit with the CCHP system, where a gas turbine cycle will be able to combust renewable hydrogen when solar energy is unavailable [11].The modular nature of both PVs and electrolyzer stacks is thereby a significant advantage in achieving these targets.
In recent publications, the authors have presented a liquefied natural gas (LNG)-fueled CCHP system, considering different system capacities (1-10 MW e ), to investigate the feasibility of such systems to penetrate the energy market in countries with hot climates requiring both cooling and heating energy [5,7].In this study a PV-assisted 1 MW e CCHP system is considered for application in a hot climatic region, where the system is assumed to be located far from existing central power plants, e.g., to fulfill the energy needs of a remote community.An interesting aspect in terms of fueling could be liquefied natural gas (LNG), which is readily available in different areas of the world.However, it should be noted that the system is not restricted to the use of LNG, and most conventional hydrocarbons could be used to fuel the CCHP system.In this study, a PV subsystem is considered for operation along with the CCHP system, in the vicinity of the latter, to reduce fuel consumption when solar energy is available.The capacity of the integrated units (i.e., available heating and cooling) will vary according to the operating scheme of the gas turbine cycle.The absorption chiller unit is assumed to be of the double-effect LiBr-H 2 O type.Therefore, this research work investigates the above options and modifications to quantify and analyze the thermoeconomic performance of the proposed system.

System Configuration
The system configuration is shown schematically in Figure 1.LNG is transported to the CCHP system location and, before combustion, natural gas is preheated and converted to gasified form.This is achieved by heat exchanger (HEx) 1, which receives atmospheric air.This air is subsequently fed to the air compressor for combustion.The generated flue gas in the combustor is expanded in the gas turbine to produce electricity through the electric generator.The hot exhausted flue gas from the gas turbine can be used in summer operation to produce steam in HEx 2, which is fed to the heat-activated cooling plant comprised of the LiBr-H 2 O absorption chiller.The produced cooling energy is distributed to nearby buildings by means of a district energy network.During winter operation, the heat is recovered by HEx 3 and distributed to the buildings with the district energy network.is distributed to nearby buildings by means of a district energy network.During winter operation, the heat is recovered by HEx 3 and distributed to the buildings with the district energy network.The study assumes 5% heat losses in the heat exchangers and the pipelines of the district energy network.To simplify calculations, natural gas is assumed to consist only of methane.Pressure losses at a rate of 5% are assumed for every component (heat exchangers, pumps, and combustor), and 10% for the district energy network.Power input to the pumps is neglected in the calculations, because it is very small compared to the power generated by the system (less than 1%).The reason is due to the fact that the energy transferred through the network is in a closed circuit, meaning that only minimum amount of power is needed to supplement pressure losses (also the energy network is rather small, with a total length of 2 km).The completely autonomous system follows an electricityled operating strategy to avoid the need of interacting with a central power grid, i.e., import/ export of electricity.The gas turbine cycle does not require a fuel (natural gas) compressor, because the fuel is already at high pressure, specifically at 4.5 MPa [5,7].Finally, the modeling of the proposed system does not include components related to regulating voltage/frequency (i.e., capacitors, batteries, etc.).It is acknowledged that in a complete system, some of these components would have to be included in the system, along with a control system.

Modeling Methodology
The models combine theoretical principles and experimental data, where applicable and available from manufacturers, or the literature.In the following subsections each mathematical model of the proposed system is described in detail and all relationships used are given along with the input values (shown in Table 1).Also details are given for the system energy balance equations, which relate the power, heating, and cooling outputs of the proposed PV-assisted CCHP system to the energy requirements of the serviced buildings (i.e., the load profile).Finally the cost model for each component and the overall system are given in the end of this section.The study assumes 5% heat losses in the heat exchangers and the pipelines of the district energy network.To simplify calculations, natural gas is assumed to consist only of methane.Pressure losses at a rate of 5% are assumed for every component (heat exchangers, pumps, and combustor), and 10% for the district energy network.Power input to the pumps is neglected in the calculations, because it is very small compared to the power generated by the system (less than 1%).The reason is due to the fact that the energy transferred through the network is in a closed circuit, meaning that only minimum amount of power is needed to supplement pressure losses (also the energy network is rather small, with a total length of 2 km).The completely autonomous system follows an electricity-led operating strategy to avoid the need of interacting with a central power grid, i.e., import/ export of electricity.The gas turbine cycle does not require a fuel (natural gas) compressor, because the fuel is already at high pressure, specifically at 4.5 MPa [5,7].Finally, the modeling of the proposed system does not include components related to regulating voltage/frequency (i.e., capacitors, batteries, etc.).It is acknowledged that in a complete system, some of these components would have to be included in the system, along with a control system.

Modeling Methodology
The models combine theoretical principles and experimental data, where applicable and available from manufacturers, or the literature.In the following subsections each mathematical model of the proposed system is described in detail and all relationships used are given along with the input values (shown in Table 1).Also details are given for the system energy balance equations, which relate the power, heating, and cooling outputs of the proposed PV-assisted CCHP system to the energy requirements of the serviced buildings (i.e., the load profile).Finally the cost model for each component and the overall system are given in the end of this section.In summary, the following procedure is carried out to generate the simulation data: 1.
The components/subsystems for the CCHP system are modeled individually from first principles, and then they are combined to produce the total CCHP system model.2.
The CCHP system is simulated at full-load and part-load to generate operating data for the nominal power output and the part-loads, i.e., 0-1 MW e .The power output of the system, along with the isentropic efficiencies for the compressor and the gas turbine of the gas turbine cycle are therefore varied, as input parameters in each simulation.Simulation data are generated for both modes of operation, i.e., summer and winter modes.The output parameters of these simulations are the following: available cooling energy in the district cooling network (DCN) or available heating energy in the district heating network (DHN) (for summer or winter operation, respectively), input fuel energy, net electrical efficiency, primary energy ratio, and mass flow rate of input fuel.

3.
The PV subsystem is modeled and operating data are generated based on Typical Meteorological Year (TMY2) data for Nicosia, Cyprus.The area of the PV array is varied to change the nominal power output of the PV subsystem, as needed.Data are generated for the whole year on an hourly basis, i.e., 8760 time segments.The output data include: the incident solar radiation, the PV array temperature, the maximum power point efficiency of the PV array, and the PV power output.

4.
The generated data from the CCHP system and the PV subsystem are then combined with the load profile to calculate the number of households that can be serviced from the system, cost data and thermodynamic data.

5.
A parametric study is conducted to investigate how different values of nominal power output for the PV subsystem affect the cost.
The system is modeled and simulated in the commercial software Engineering Equation Solver (EES) Professional version, which offers the capability of carrying out parametric studies with different built-in tools, such as parametric tables and thermo-physical data.

Gas Turbine Cycle
The gas turbine cycle includes an air compressor, a combustor and a gas turbine.The cooling energy recovered by the LNG regasification process improves overall system efficiency, since the specific volume of the air is reduced [5,7].The enthalpy at the air compressor at exit is: where h 1 is the enthalpy at inlet and h s,2 is the isentropic enthalpy at exit.The combustion process can be modeled with an energy balance including the reactants (methane, oxygen and nitrogen) and the products (carbon dioxide, water, nitrogen and oxygen).The produced net electrical power output through the gas turbine process, in relation to the consumed fuel, is given as a function of the molar flow rate of methane, enthalpy of the products exiting the combustor per mole of fuel: .
The net electrical efficiency of the system is defined as the ratio of the net electrical power output to the chemical energy of the fuel (based on the lower heating value (LHV)):

Heat Exchangers
The heat exchangers, shown in Figure 1, assume energy balance modeling (cold and hot sides).Therefore the heat transfer rate of the hot and cold sides can be defined as follows: where .
m is the mass flow rate (cold and hot sides) and h is the specific enthalpy.The LNG regasification process is shown to be a simple heat exchanger process in this study.However, it should be noted that this process is much more complicated in an actual system, but since this is not the scope of the study, a simple thermodynamic process is assumed to be adequate.

Absorption Chiller
The absorption chiller is of the LiBr-H 2 O type, since this is an efficient chiller with coefficient of performance (COP) values of 1.2-1.3(double-effect) [12,13].It is suitable for air-conditioning applications, since it generates a cooling output of 5-10 • C in the form of cold water at the evaporator output.A double-effect chiller is chosen over a single-effect one, since the heat source (steam generated in a heat recovery steam generator) is at a suitable temperature to operate a double-effect absorption chiller.The absorption chiller follows the modeling explained in detail by some of the authors in [14].

District Energy Network
As mentioned above, the district energy network consists of two operating modes: summer and winter modes.When operating at summer mode, heat recovered from the gas turbine exhaust's flue gas is used to activate the absorption chiller, generating cold water [13].The cold water is distributed to the building from the CCHP location through the district cooling network's pipelines for space cooling application at a temperature of 22 • C [7].When operating in winter mode, the heat is recovered through a heat exchanger and distributed to the buildings through the district heating network to generate domestic hot water and space heating at a temperature of 50 • C and 22 • C, respectively [7].

Photovoltaic Panels
The PV subsystem is integrated into the CCHP system to reduce fuel consumption.The PV model is based on the theoretical assumptions found in [15].The model calculates the total incident solar radiation on a tilted surface, which considers both the ground-reflected and beam terms, as defined in the results of the HDKR (Hay, Davies, Klucher, and Reindl) model: The temperature of the PV array, T c , is calculated by use of the following relation (the effect of the wind speed is neglected): The maximum power point efficiency of the PV array is given as follows: The power output of the PV array is given by:

Load Profile and System Balance
Since the proposed PV-assisted CCHP system is completely autonomous (no power and/or heating/cooling is imported and/or exported from/to the central power grid/central district energy network, respectively), the selected operating strategy must ensure that it will be able to completely fulfill the load profile [16].For this study the consumption data for the load profile for an average 100 m 2 household in Cyprus are adjusted on an hourly basis.Since power is the most expensive energy type, no surplus power should be produced, and therefore an electricity-led operation must be followed.Since the system will be operated at both full-load and part-load, it means that power-to-heating and power-to-cooling ratios will vary accordingly.Since the available thermal energy output and, in effect, the district energy delivered to the buildings will be inadequate for certain loads, it is necessary to generate additional thermal energy (heating or cooling) with other means.Specifically, thermal energy can be generated with vapor-compression units, i.e., central electric chillers (located at the site of the CCHP system) or with electric heat pumps located inside the buildings.It is assumed in this study that the proposed CCHP system will service buildings that are already equipped with electric heat pumps, and therefore central electric chillers are not considered, as these would add an unnecessary extra cost to the system.However, for new buildings, central chiller units would be a more ideal option, since overall efficiency and performance would improve.The average COP value for the electric heat pumps is assumed to be 3.0.The cooling load is the space cooling energy, whereas the heating load is the combined space heating and domestic hot water load for the buildings.The electrical load is the load needed for lighting and appliances.The electrical energy for the operation of the heat pumps is added in the energy balance, which is shown below.
The electrical load profile (excluding the power input to the electric heat pumps), P cs , is the product of the electrical energy consumption per household, P cs,hs , and the number of households, n: Similarly, the heating load profile, Q h,cs , is the product of the heating energy consumption per household, Q h,cs,hs , and the number of households; while the cooling load profile, Q c,cs , is the product of the cooling energy consumption per household, Q c,cs,hs , and the number of households: The available heating energy output, .
Q DHN , in the district energy network during winter operation is calculated through interpolation of the CCHP system model data for a given input CCHP electrical energy power output, .W gen .Similarly, the fuel input to the system, E f uel,in , net electrical efficiency, η el,net , and fuel mass flow rate, .m LNG , are calculated through interpolation.The available cooling energy output, .Q DCN , is calculated in the same manner.Since a heat-led operating strategy is not followed, the available heating or cooling energy generated by the CCHP system will not match the heating and cooling load profiles.The power input to the electric heat pumps must be provided with additional electrical energy generation from the CCHP system.However, the additional electrical energy generation also translates to additional thermal energy generation.Therefore, an energy balance, which includes all of these parameters, must be formulated to ensure that the system will be able to fully and accurately match the load profile requirements.
The amount of available heating energy in the district energy network which is actually used in the load profile is given as follows: Similarly, the amount of available cooling energy in the district energy network which is actually used in the load profile is given as follows: The additional heating and cooling energy generated with the electric heat pumps, Q hp,h , and Q hp,c , respectively, is used to calculate the required power input to the heat pump, P hp,h , and P hp,c , respectively.
The overall energy balance of the PV-assisted CCHP system can be formulated to include all of the above parameters, and calculate the electrical energy power output of the CCHP-system, P gen , for any given time segment: P gen = max 0, P cs + P hp,h + P hp,c − P pv (18) The last term, P pv , is the power output of the PV subsystem.
To assess the performance of the proposed system, an appraisal criterion parameter, called the "primary energy ratio" (PER), is defined: It should be clarified that this definition is preferred over "total system efficiency", since the latter could be a misleading definition if related to thermodynamic efficiency, since the proposed system combines both thermal and refrigeration cycles.

Cost Model
The cost model, which is partly based on the cost analysis found in Dincer & Zamfirescu [17], includes all cost functions that are needed to evaluate the projected total LCC of each system configuration.As mentioned above, the LCC is calculated in a parametric study, where the PV power output is varied, ranging from 0 to 600 kW e .The cost model is integrated to the PV-assisted CCHP system model, and all the included cost functions are shown in Table 2. Table 3 includes the values for the cost factor input parameters, and Table 4 includes the values for the cost input parameters.These values are taken from [17], unless specified otherwise.The cost of LNG, c LNG , is taken as 3.082 €/GJ [18], while the system lifetime is approximated at 20 years of service [17].The specific cost for the PV subsystem is set at 1.825 €/W, based on an approximation of values found in [9].

Q hcy
Total thermal energy production from combined cooling, heating, and power (CCHP) system (J) Total cost of fan-coil units (€) PA n Uniform series present worth factor (−)

FP n
Compound amount factor (−) Tax deduction on the loan (€) C omi Operation, maintenance and insurance cost (€) The volumetric flow rate of LNG is calculated as follows: It is assumed that each LNG truck has a capacity of truck LNG = 58 m 3 [19].The LNG supply times per year can be calculated as follows: The fan-coil units are assumed to have an average power input, .
Q f c , of 1 kW per unit [17].

Results and Discussion
In this section the results for the proposed PV-assisted CCHP system are analyzed in terms of thermodynamics and cost.Finally, a parametric study is used to evaluate the system configuration resulting in minimum LCC.

Validation of the Model
The CCHP system model is validated with reported data for commercially available gas turbine systems and the results of the validation show good agreement between these data and the modeling.The validation has been shown in detail by the authors in [7].

Simulation of the CCHP System and the PV Subsystem
The CCHP system is simulated at part-load and full-load, i.e., in a power output ranging between 0 and 1 MW e .As shown in Figure 2, the useful energy output of the CCHP system depends on the selected load; at full-load the CCHP system has a heating-to-power ratio of 1.72, and a cooling-to-power ratio of 2.00; at half-load these values are 2.45 and 2.98, respectively.The increase of these ratios is due to the comparatively higher heating available from the gas turbine exhaust, since the gas turbine cycle becomes less efficient (in terms of electrical energy output) as the load decreases.This is shown more clearly in Figure 3, where the PER is shown to remain almost constant at all loads, due to the way it was defined in the previous section, since the electrical energy output decreases; at the same time, the useful heating (or cooling) energy output increases as the load decreases.The PV subsystem integrated into the CCHP system aims to reduce the fuel consumption through the partial production of the power needed to satisfy the load profile.The power output of the PV panels depends on the availability of solar radiation.For example, for a PV subsystem with a nominal power output of 150 kWe, the pattern of power generation is shown in Figure 4, with the electrical energy output (in W-h) shown on an hourly basis for a whole year as applied for the weather conditions in Nicosia, Cyprus.The average power output is 34.1 kWe.The PV subsystem integrated into the CCHP system aims to reduce the fuel consumption through the partial production of the power needed to satisfy the load profile.The power output of the PV panels depends on the availability of solar radiation.For example, for a PV subsystem with a nominal power output of 150 kWe, the pattern of power generation is shown in Figure 4, with the electrical energy output (in W-h) shown on an hourly basis for a whole year as applied for the weather conditions in Nicosia, Cyprus.The average power output is 34.1 kWe.The PV subsystem integrated into the CCHP system aims to reduce the fuel consumption through the partial production of the power needed to satisfy the load profile.The power output of the PV panels depends on the availability of solar radiation.For example, for a PV subsystem with a nominal power output of 150 kW e , the pattern of power generation is shown in Figure 4, with the electrical energy output (in W-h) shown on an hourly basis for a whole year as applied for the weather conditions in Nicosia, Cyprus.The average power output is 34.1 kW e .

Application of the PV-assisted CCHP System to a Selected Load Profile
The load profile consists of the following consumption data: (a) power, i.e., the electrical energy needed to satisfy all building needs (including the heat pumps); (b) heating needed for space heating and domestic hot water (winter operation mode); (c) cooling needed for space cooling (summer operation mode).The summer operation mode assumes that domestic hot water will be available through other means (e.g.solar collectors) and therefore it is neglected in the calculations.The load profile is shown in Figure 5. From an observation, the electrical load profile remains constant throughout the year, while the heating load profile starts in November and ends in April.The cooling load profile starts in May and ends in September.There are some small periods of time between the winter and summer periods, where neither heating nor cooling is required (the so called "shoulder" months).Therefore, from this pattern, it is evident that the CCHP system will need to provide additional electrical energy to operate the heat pumps when heating (or cooling) generation from the CCHP system will be insufficient.The CCHP system can completely fulfill the load profile of 436 households having the characteristics of the load profile shown in Figure 5. Without the integration of a PV subsystem, the CCHP system must generate a total electrical energy of 1959 MWh per year.This figure includes both the electrical load profile and the operation of the electric heat pumps.The latter requires 146 MWh of electrical energy per year.The total fuel input energy is 9825 MWh per year (or 780 tons), generating

Application of the PV-assisted CCHP System to a Selected Load Profile
The load profile consists of the following consumption data: (a) power, i.e., the electrical energy needed to satisfy all building needs (including the heat pumps); (b) heating needed for space heating and domestic hot water (winter operation mode); (c) cooling needed for space cooling (summer operation mode).The summer operation mode assumes that domestic hot water will be available through other means (e.g., solar collectors) and therefore it is neglected in the calculations.The load profile is shown in Figure 5. From an observation, the electrical load profile remains constant throughout the year, while the heating load profile starts in November and ends in April.The cooling load profile starts in May and ends in September.There are some small periods of time between the winter and summer periods, where neither heating nor cooling is required (the so called "shoulder" months).Therefore, from this pattern, it is evident that the CCHP system will need to provide additional electrical energy to operate the heat pumps when heating (or cooling) generation from the CCHP system will be insufficient.

Application of the PV-assisted CCHP System to a Selected Load Profile
The load profile consists of the following consumption data: (a) power, i.e., the electrical energy needed to satisfy all building needs (including the heat pumps); (b) heating needed for space heating and domestic hot water (winter operation mode); (c) cooling needed for space cooling (summer operation mode).The summer operation mode assumes that domestic hot water will be available through other means (e.g.solar collectors) and therefore it is neglected in the calculations.The load profile is shown in Figure 5. From an observation, the electrical load profile remains constant throughout the year, while the heating load profile starts in November and ends in April.The cooling load profile starts in May and ends in September.There are some small periods of time between the winter and summer periods, where neither heating nor cooling is required (the so called "shoulder" months).Therefore, from this pattern, it is evident that the CCHP system will need to provide additional electrical energy to operate the heat pumps when heating (or cooling) generation from the CCHP system will be insufficient.The CCHP system can completely fulfill the load profile of 436 households having the characteristics of the load profile shown in Figure 5. Without the integration of a PV subsystem, the CCHP system must generate a total electrical energy of 1959 MWh per year.This figure includes both the electrical load profile and the operation of the electric heat pumps.The latter requires 146 MWh of electrical energy per year.The total fuel input energy is 9825 MWh per year (or 780 tons), generating The CCHP system can completely fulfill the load profile of 436 households having the characteristics of the load profile shown in Figure 5. Without the integration of a PV subsystem, the CCHP system must generate a total electrical energy of 1959 MWh per year.This figure includes both the electrical load profile and the operation of the electric heat pumps.The latter requires 146 MWh of electrical energy per year.The total fuel input energy is 9825 MWh per year (or 780 tons), generating 2311 tons of CO 2 .The average net electrical efficiency of the system is 19.6%, while the average PER is 39.4%.

Parametric Study
The effect of the integration of the PV subsystem to the CCHP system is shown through a parametric study, where the size of the PV subsystem is varied ranging from 0 to 600 kW e .Figure 6 shows the variation of useful energy output.It is obvious that less heating (or cooling) is generated as the PV capacity increases, since less natural gas is combusted in the gas turbine cycle.As the PV subsystem increases, it is also better for the overall performance of the system, as more heating (or cooling), which is available, is actually used in the load profile.This is shown more clearly in Figure 7, where the average PER increases as the PV nominal power output increases, at the expense, however, of lower average net electrical efficiency.This is inevitable, since the gas turbine cycle is forced to operate at lower loads and, in the case of a 600 kW e PV subsystem, the gas turbine basically operates either at peak loads (i.e., during daytime), or when solar energy is unavailable (i.e., during nighttime).Figure 8 shows the significant reduction in fuel consumption and CO 2 emissions, as the nominal power output of the PV subsystem increases.This has significant benefits in operating cost, and environmentally it can help reduce greenhouse gas emissions, especially if the system is situated in sensitive locations where this is desirable.
Energies 2016, 9, 663 13 of 16 2311 tons of CO2.The average net electrical efficiency of the system is 19.6%, while the average PER is 39.4%.

Parametric Study
The effect of the integration of the PV subsystem to the CCHP system is shown through a parametric study, where the size of the PV subsystem is varied ranging from 0 to 600 kWe. Figure 6 shows the variation of useful energy output.It is obvious that less heating (or cooling) is generated as the PV capacity increases, since less natural gas is combusted in the gas turbine cycle.As the PV subsystem increases, it is also better for the overall performance of the system, as more heating (or cooling), which is available, is actually used in the load profile.This is shown more clearly in Figure 7, where the average PER increases as the PV nominal power output increases, at the expense, however, of lower average net electrical efficiency.This is inevitable, since the gas turbine cycle is forced to operate at lower loads and, in the case of a 600 kWe PV subsystem, the gas turbine basically operates either at peak loads (i.e. during daytime), or when solar energy is unavailable (i.e. during nighttime).Figure 8 shows the significant reduction in fuel consumption and CO2 emissions, as the nominal power output of the PV subsystem increases.This has significant benefits in operating cost, and environmentally it can help reduce greenhouse gas emissions, especially if the system is situated in sensitive locations where this is desirable.2311 tons of CO2.The average net electrical efficiency of the system is 19.6%, while the average PER is 39.4%.

Parametric Study
The effect of the integration of the PV subsystem to the CCHP system is shown through a parametric study, where the size of the PV subsystem is varied ranging from 0 to 600 kWe. Figure 6 shows the variation of useful energy output.It is obvious that less heating (or cooling) is generated as the PV capacity increases, since less natural gas is combusted in the gas turbine cycle.As the PV subsystem increases, it is also better for the overall performance of the system, as more heating (or cooling), which is available, is actually used in the load profile.This is shown more clearly in Figure 7, where the average PER increases as the PV nominal power output increases, at the expense, however, of lower average net electrical efficiency.This is inevitable, since the gas turbine cycle is forced to operate at lower loads and, in the case of a 600 kWe PV subsystem, the gas turbine basically operates either at peak loads (i.e. during daytime), or when solar energy is unavailable (i.e. during nighttime).Figure 8 shows the significant reduction in fuel consumption and CO2 emissions, as the nominal power output of the PV subsystem increases.This has significant benefits in operating cost, and environmentally it can help reduce greenhouse gas emissions, especially if the system is situated in sensitive locations where this is desirable.The cost analysis for the proposed system is shown in Figure 9, where LCC is varied depending on the size of the nominal PV power output.Interestingly, LCC drops as the PV subsystem capacity increases, but only up to 300 kWe.After this point the LCC increases and, at 600 kWe, it reaches its maximum value for this parametric study.The reason is that at 600 kWe the useful power output of the PV subsystem is not enough to account for the higher capital cost of the PV subsystem.In other words, more electrical energy is available for consumption, but since the PV subsystem is oversized, it cannot match the requirements of the load profile.Therefore, in this case, electricity generated from the PV subsystem is wasted.The minimum, and thereby the optimum, value for the LCC is at 300 kWe of nominal power output for the PV subsystem.This pattern could change if the unit cost of PV is decreased in the future.In this case, a higher capacity for the PV subsystem could be desirable.The effect of the specific cost of the PV subsystem on the LCC of the proposed CCHP system is shown in Figure 10, where the specific cost of the PV subsystem is varied from 0.730 to 1.825 €/kWe.Specifically, a PV subsystem with a capacity of 600 kWe results in the lowest LCC in the case of a pv c at 0.730 €/kWe.In this case the LCC is 3.372 million €.The cost analysis for the proposed system is shown in Figure 9, where LCC is varied depending on the size of the nominal PV power output.Interestingly, LCC drops as the PV subsystem capacity increases, but only up to 300 kW e .After this point the LCC increases and, at 600 kW e , it reaches its maximum value for this parametric study.The reason is that at 600 kW e the useful power output of the PV subsystem is not enough to account for the higher capital cost of the PV subsystem.In other words, more electrical energy is available for consumption, but since the PV subsystem is oversized, it cannot match the requirements of the load profile.Therefore, in this case, electricity generated from the PV subsystem is wasted.The minimum, and thereby the optimum, value for the LCC is at 300 kW e of nominal power output for the PV subsystem.This pattern could change if the unit cost of PV is decreased in the future.In this case, a higher capacity for the PV subsystem could be desirable.The effect of the specific cost of the PV subsystem on the LCC of the proposed CCHP system is shown in Figure 10, where the specific cost of the PV subsystem is varied from 0.730 to 1.825 €/kW e .Specifically, a PV subsystem with a capacity of 600 kW e results in the lowest LCC in the case of a c pv at 0.730 €/kW e .In this case the LCC is 3.372 million €.The cost analysis for the proposed system is shown in Figure 9, where LCC is varied depending on the size of the nominal PV power output.Interestingly, LCC drops as the PV subsystem capacity increases, but only up to 300 kWe.After this point the LCC increases and, at 600 kWe, it reaches its maximum value for this parametric study.The reason is that at 600 kWe the useful power output of the PV subsystem is not enough to account for the higher capital cost of the PV subsystem.In other words, more electrical energy is available for consumption, but since the PV subsystem is oversized, it cannot match the requirements of the load profile.Therefore, in this case, electricity generated from the PV subsystem is wasted.The minimum, and thereby the optimum, value for the LCC is at 300 kWe of nominal power output for the PV subsystem.This pattern could change if the unit cost of PV is decreased in the future.In this case, a higher capacity for the PV subsystem could be desirable.The effect of the specific cost of the PV subsystem on the LCC of the proposed CCHP system is shown in Figure 10, where the specific cost of the PV subsystem is varied from 0.730 to 1.825 €/kWe.Specifically, a PV subsystem with a capacity of 600 kWe results in the lowest LCC in the case of a pv c at 0.730 €/kWe.In this case the LCC is 3.372 million €.

Conclusions
In this study the effects of integrating a proposed 1 MWe CCHP system, fueled by LNG, with a PV subsystem are analyzed in terms of thermoeconomic modeling and a parametric study.The study shows that the PV subsystem can significantly reduce the power generation (and thereby fuel consumption and CO2 emissions) needed from the CCHP system.In the parametric study, PV capacity (nominal power output) is varied from 0 to 600 kWe, and the optimum PV capacity is found to be 300 kWe.The minimum LCC for the PV-assisted CCHP system is found to be 3.509 million €, as compared to 3.577 million € for a system without a PV subsystem.The total cost of the PV subsystem is 547,445 €, while the total cost for operating the system (fuel) is 731,814 € (compared to 952,201 € for a CCHP system without PVs).These cost savings could increase in the future as the unit PV cost decreases and the unit cost of natural gas increases.
Overall, the proposed PV-assisted CCHP system generates the following amounts of useful energy (energy that is actually utilized in the load profile throughout the total lifetime (20 years) of the system): 83,140 kWh of electrical energy, 60,960 kWh of heating energy and 66,420 kWh of cooling energy.Therefore, for a total energy output of 210,520 kWh, this results in a unit cost of 17 €/kWh.Additionally, the study shows that if the specific cost of the PV subsystem drops in the future, eventually the gas turbine (or any other conventional prime mover) will switch roles with the PV subsystem.In other words, the gas turbine could act as an auxiliary power generator for periods of time where no solar radiation will be available.Finally, better utilization of the available heating and cooling energy is shown for the CCHP system with a PV subsystem, since the system becomes more electricity-dependent, offering the capability of converting electricity to heating (or cooling) through the electric heat pumps.Thereby, more control and flexibility is provided to adjust the generation of thermal energy as dictated by the varying load profile.

Figure 1 .
Figure 1.Schematic configuration of the proposed photovoltaic (PV)-assisted 1 MWe combined cooling, heating, and power system.

Figure 1 .
Figure 1.Schematic configuration of the proposed photovoltaic (PV)-assisted 1 MW e combined cooling, heating, and power system.

Energies 2016, 9 , 663 11 of 16 Figure 2 .
Figure 2. Simulation of the 1 MWe combined cooling, heating, and power system at part-load and fullload: fuel input vs. useful energy output.

Figure 3 .
Figure 3. Simulation of the 1 MWe combined cooling, heating, and power system at part-load and fullload: variation of efficiency at full-load and part-load operation.

Figure 2 . 16 Figure 2 .
Figure 2. Simulation of the 1 MW e combined cooling, heating, and power system at part-load and full-load: fuel input vs. useful energy output.

Figure 3 .
Figure 3. Simulation of the 1 MWe combined cooling, heating, and power system at part-load and fullload: variation of efficiency at full-load and part-load operation.

Figure 3 .
Figure 3. Simulation of the 1 MW e combined cooling, heating, and power system at part-load and full-load: variation of efficiency at full-load and part-load operation.

Figure 4 .
Figure 4. Generation of electrical energy in the photovoltaic subsystem throughout the year.

Figure 5 .
Figure 5. Annual load profile (power, heating, and cooling) for the average household considered in the study on an hourly basis.

Figure 4 .
Figure 4. Generation of electrical energy in the photovoltaic subsystem throughout the year.

Figure 4 .
Figure 4. Generation of electrical energy in the photovoltaic subsystem throughout the year.

Figure 5 .
Figure 5. Annual load profile (power, heating, and cooling) for the average household considered in the study on an hourly basis.

Figure 5 .
Figure 5. Annual load profile (power, heating, and cooling) for the average household considered in the study on an hourly basis.

Figure 6 .
Figure 6.Parametric study: annual useful energy output for different photovoltaic subsystem capacities ranging from 0 to 600 kWe.

Figure 7 .
Figure 7. Parametric study: annual average net electrical efficiency and average primary energy ratio (PER) for different photovoltaic subsystem capacities ranging from 0 to 600 kWe.

Figure 6 .
Figure 6.Parametric study: annual useful energy output for different photovoltaic subsystem capacities ranging from 0 to 600 kW e .

Figure 6 .
Figure 6.Parametric study: annual useful energy output for different photovoltaic subsystem capacities ranging from 0 to 600 kWe.

Figure 7 .
Figure 7. Parametric study: annual average net electrical efficiency and average primary energy ratio (PER) for different photovoltaic subsystem capacities ranging from 0 to 600 kWe.

Figure 7 .
Figure 7. Parametric study: annual average net electrical efficiency and average primary energy ratio (PER) for different photovoltaic subsystem capacities ranging from 0 to 600 kW e .

Figure 8 .
Figure 8. Parametric study: annual consumption of fuel and CO2 emissions for different photovoltaic subsystem capacities ranging from 0 to 600 kWe.

Figure 9 .
Figure 9. Parametric study: cost analysis for different photovoltaic subsystem capacities ranging from 0 to 600 kWe.

Figure 8 .
Figure 8. Parametric study: annual consumption of fuel and CO 2 emissions for different photovoltaic subsystem capacities ranging from 0 to 600 kW e .

Figure 9 .
Figure 9. Parametric study: cost analysis for different photovoltaic subsystem capacities ranging from 0 to 600 kWe.

Figure 9 .
Figure 9. Parametric study: cost analysis for different photovoltaic subsystem capacities ranging from 0 to 600 kW e .

Figure 10 .
Figure 10.Effect of the specific cost for the photovoltaic subsystem on the lifecycle cost of the proposed system.

Table 1 .
Values of the system input parameters.

Table 2 .
Cost model for the proposed photovoltaic-assisted combined cooling, heating, and power system.Annual cost of fuel (incl.regasification) (€/year) c f y = c f yy + c reg C tr Total cost of LNG transport (€) C tr = c tr,km L tr V truck supply LNG,yr N C cc Total cost of chillers (€)

Table 3 .
Input cost factor values.

Table 4 .
Input values for the cost model.