Development and Validation of a Thermo-Economic Model for Design Optimisation and Off-Design Performance Evaluation of a Pure Solar Microturbine

The aim of this paper is to present a thermo-economic model of a microturbine for solar dish applications, which demonstrates the applicability and accuracy of the model for off-design performance evaluation and techno-economic optimisation purposes. The model is built using an object-oriented programming approach. Each component is represented using a class made of functions that perform a one-dimensional physical design, off-design performance analysis and the component cost evaluation. Compressor, recuperator, receiver and turbine models are presented and validated against experimental data available in literature, and each demonstrated good accuracy for a wide range of operating conditions. A 7-kWe microturbine and solar irradiation data available for Rome between 2004 and 2005 were considered as a case study, and the thermo-economic analysis of the plant was performed to estimate the levelised cost of electricity based on the annual performance of the plant. The overall energy produced by the plant is 10,682 kWh, the capital cost has been estimated to be EUR 27,051 and, consequently, the specific cost of the plant, defined as the ratio between the cost of components and output power in design condition, has been estimated to be around EUR 3980/kWe. Results from the levelised cost of electricity (LCOE) analysis demonstrate a levelised cost of electricity of EUR 22.81/kWh considering a plant lifetime of 25 years. The results of the present case study have been compared with the results from IPSEpro 7 where the same component characteristic maps and operational strategy were considered. This comparison was aimed to verify the component matching procedure adopted for the present model. A plant sizing optimisation was then performed to determine the plant size which minimises the levelised cost of electricity. The design space of the optimisation variable is limited to the values 0.07–0.16 kg/s. Results of the optimisation demonstrate a minimum LCOE of 21.5 [EUR/kWh] for a design point mass flow rate of about 0.11 kg/s. This corresponds to an overall cost of the plant of around EUR 32,600, with a dish diameter of 9.4 m and an annual electricity production of 13,700 [kWh].


Introduction
The rapid growth of the world population, as well as the expansion of the industrial sector, has led to a continuously growing demand for electrical power. This growth, coupled with the need for reduction of greenhouse gas emissions, has created serious global challenges, e.g., achieving energy demand without violating increasingly stricter environmental regulations. This challenge has led to a CSP systems are generally designed and analysed based on the nominal design point; given the intermittence of the solar radiation, they operate for most of the time in off-design conditions. Modelling plant performance is then an important aspect for detailed analysis and design and operation optimisation. The accuracy of the models and their applicability are key aspects to consider in software development. Mazzoni et al. [7] developed a simulation tool for a pure solar dish-micro gas turbine application tailored for the OMSoP project. The set-up of the CSP plant simulator is built on a modular arrangement at the component level as presented by Cerri [8]. A library of suitable component models for the plant was realised. The simulator has been adopted to perform steady state simulations for predicting the performance under variable direct normal irradiation. The biggest advantage of the modular approach is the prospect to generate any possible plant layout and integrate new components. This characteristic creates the possibility to include the desired number of compression and expansion stages and the option to integrate any number of heat exchangers for intercooled recuperative cycles. Moreover, modularity of the tool permits an easy integration of new components, such as a combustion chamber for hybrid operations and reheating cycles or thermal storage units.
On the other hand, for optimisation purposes it is important to evaluate different plant designs considering dissimilar main cycle design parameters. To evaluate the correct size of the plant is then necessary to design and evaluate off-design performances of a wide range of possible component assembly solutions. In this situation new components must be designed. In most of the commercial software available on the market for power plant analysis, the user is often provided with a library of components to choose from, resulting in less or limited flexibility, lower applicability in terms of component selection and poor adaptability to optimisation processes.
Gavagnin et al. [6] presented a paper aimed to find the optimum design and performance of a solar micro-turbine powered by dish collectors using an innovative methodology which integrates the design and off-design models of the system, tailored over existing components of the project OMSoP. This was in contrast to the common practice of assigning an estimated efficiency to the engine turbomachinery. The design model was run for two different sets of turbine inlet temperature and recuperator effectiveness (TIT and ε): 800 °C, 85% and 900 °C, 90%. Different systems were designed for each case, yielding different off-design performance maps that link the net power output of the system to the boundary conditions of the site (ambient temperature and direct normal irradiation or DNI). With this information, annual simulations were performed in three different locations (Beijing, Seville and San Diego) based on standard performance maps for a reference DNI and on designs tailored to their particular DNI.
Ghavami [9] developed a computational model for pure solar dish-MGT systems. The model combines the cycle analysis of the MGT with component models to perform design point performance simulation, generate component performance maps and perform off-design CSP systems are generally designed and analysed based on the nominal design point; given the intermittence of the solar radiation, they operate for most of the time in off-design conditions. Modelling plant performance is then an important aspect for detailed analysis and design and operation optimisation. The accuracy of the models and their applicability are key aspects to consider in software development. Mazzoni et al. [7] developed a simulation tool for a pure solar dish-micro gas turbine application tailored for the OMSoP project. The set-up of the CSP plant simulator is built on a modular arrangement at the component level as presented by Cerri [8]. A library of suitable component models for the plant was realised. The simulator has been adopted to perform steady state simulations for predicting the performance under variable direct normal irradiation. The biggest advantage of the modular approach is the prospect to generate any possible plant layout and integrate new components. This characteristic creates the possibility to include the desired number of compression and expansion stages and the option to integrate any number of heat exchangers for intercooled recuperative cycles. Moreover, modularity of the tool permits an easy integration of new components, such as a combustion chamber for hybrid operations and reheating cycles or thermal storage units.
On the other hand, for optimisation purposes it is important to evaluate different plant designs considering dissimilar main cycle design parameters. To evaluate the correct size of the plant is then necessary to design and evaluate off-design performances of a wide range of possible component assembly solutions. In this situation new components must be designed. In most of the commercial software available on the market for power plant analysis, the user is often provided with a library of components to choose from, resulting in less or limited flexibility, lower applicability in terms of component selection and poor adaptability to optimisation processes.
Gavagnin et al. [6] presented a paper aimed to find the optimum design and performance of a solar micro-turbine powered by dish collectors using an innovative methodology which integrates the design and off-design models of the system, tailored over existing components of the project OMSoP. This was in contrast to the common practice of assigning an estimated efficiency to the engine turbomachinery. The design model was run for two different sets of turbine inlet temperature and recuperator effectiveness (TIT and ε): 800 • C, 85% and 900 • C, 90%. Different systems were designed for each case, yielding different off-design performance maps that link the net power output of the system to the boundary conditions of the site (ambient temperature and direct normal irradiation or DNI). With this information, annual simulations were performed in three different locations (Beijing, Seville and San Diego) based on standard performance maps for a reference DNI and on designs tailored to their particular DNI.
Ghavami [9] developed a computational model for pure solar dish-MGT systems. The model combines the cycle analysis of the MGT with component models to perform design point performance simulation, generate component performance maps and perform off-design performance evaluation. The computational model has been coupled with an up-to-date economic model which was specifically developed through the OMSoP project for dish-MGT systems.
The approach to dish technologies has been mostly technical in the past, although for a complete assessment of a solar dish MGT power plant an economic analysis is required. More recently, Semprini et al. [10] presented a thermo-economic evaluation of a dish-MGT solar power plant during hybrid and pure solar operation. The study included a comprehensive cost analysis covering manufacturing, transportation and installation for a dish and MGT solar power generator, considering pure solar and hybrid operations. The study demonstrated that for the pure solar system, the average system purchase equipment cost, considering a production volume of 1000 units/year, was EUR 3262/kW e for the solar-only configuration and EUR 3302/kW e for the hybrid version. Regarding the dependence on MGT size, the minimum specific cost seems to be obtained for an air mass flow rate in the order of 200 g/s even if further, more refined optimisation is still needed to find the true optimal combination of size, technical specifications, layouts and configuration.
Ghavami [6] also performed a thermo-economic optimisation of a pure solar MGT. The model already described in the previous paragraph was coupled to an optimisation tool to find system designs which led to optimal thermo-economic performance for a 5-kW e system. Then, the optimisation was extended over the rated power of 5-30 kWe to find the power rating which results in the minimum cost of generated electricity by the dish-MGT systems. To generate new performance maps of the components the model adopted scaling techniques.
Despite years of development, solar micro turbines are still in an early stage of development and needs further technological progress in order to achieve competitiveness and to be ready for a global scale commercialisation. This is especially true for the MGT dish, for which development is still at a demonstration level and, in general, the following flaws have been identified in the cited previous works in the topic:

•
Low aptitude for optimisation. Most existing tools and studies are based on existing plants, resulting in reduced capability of the tool to design and evaluate accurately off-design performance of different components. • Scarce flexibility. Most existing tools do not allow the possibility to integrate new components and consequently compare dissimilar layouts. • Few thermo-economic assessments and optimisations concerning solar dish-MGT, and none considering integrated hybridization (fuel back-up) and/or thermal storage.
In this paper, a validated model of a microturbine for solar dish applications is presented to be used for techno-economic optimisation purposes. The aim of this study is to create a flexible and modular tool for design and off-design performance evaluation, which can be adopted for different plant layouts and operational strategies. The developed tool is not intended to be used only for optimisation purposes but can also be adopted to evaluate the performance of existing plants and components. In the first part of the paper the logic of the tool, together with the components' models, are presented and validated. In the following results section the model is adopted for a thermo-economic assessment of a 7-kW e power plant. Finally, a plant size optimisation is presented to find the optimum size of the plant in respect of the annual power output and the component's cost.

Materials and Methods
In this section, an overview of the model is presented. First, the logic approach of the software is introduced together with the components' model description and validation. An object-oriented programming approach was used to guarantee modularity of the components in the code and to enable the possibility to integrate new modules in the future, such as combustion chambers and thermal storage. Each module represents a plant's component, performing design and off-design calculations. A case study follows to demonstrate the capability of the software for performance prediction and cost evaluation. Finally, a plant-size optimisation is performed to evaluate the best size of the plant that can guarantee the right trade-off between plant cost and annual power generation.
To evaluate the performance of the plant in different operational conditions, a model for each component of the plant has been developed using a modular approach. The flowchart of the model is shown in Figure 2 where the generic Jth module is characterised by two functions. The first is the 'design function' which calculates the design performance and the geometrical characteristics of the component based on the desired set of design input variables. The second is the 'off-design function' which evaluates the performance of the plant at off-design conditions based on the geometrical characteristics. The latter is calculated by the design function, the required set of input variables at off-design conditions and the operational strategy, for instance constant TIT. The model has been programmed using the object-oriented language C++.
The biggest advantage of the modular approach is the possibility to generate any desired plant layout. This characteristic allows the possibility of including the desired number of compression and expansion stages, the possibility to integrate any number of heat exchangers for intercooled recuperative cycles, or the possibility to integrate a combustion chamber for hybrid operations and reheating cycles. To evaluate the performance of the plant in different operational conditions, a model for each component of the plant has been developed using a modular approach. The flowchart of the model is shown in Figure 2 where the generic jth module is characterised by two functions. The first is the 'design function' which calculates the design performance and the geometrical characteristics of the component based on the desired set of design input variables. The second is the 'off-design function' which evaluates the performance of the plant at off-design conditions based on the geometrical characteristics. The latter is calculated by the design function, the required set of input variables at off-design conditions and the operational strategy, for instance constant TIT. The model has been programmed using the object-oriented language C++.
The biggest advantage of the modular approach is the possibility to generate any desired plant layout. This characteristic allows the possibility of including the desired number of compression and expansion stages, the possibility to integrate any number of heat exchangers for intercooled recuperative cycles, or the possibility to integrate a combustion chamber for hybrid operations and reheating cycles. In this study, a recuperative cycle is presented. Once the components are selected, the first calculation carried out by the software is to design the plant and its main components using design point input parameters. After the design procedure is completed, the program evaluates the off-design performance of the components based on the desired operational strategy. The off-design procedure starts by assuming a first attempt mass flow rate and iterates until all the components match each other. During iteration the matching guess is continuously updated until component matching is achieved. Further constraints can be included in the software according to the operational strategy. In this case another iteration is required and the system is solved via serial nested loops. The loops matching guesses and constraints are solved in a nested sequence until convergence.

Fluid Properties Module
The fluid properties are characterised using a specific class, which can be accessed by any components' classes. Through this module, each working fluid can be described as an object with its composition and thermodynamic properties. The specific heat of the working fluid as a function of temperature is evaluated using the first order polynomial regression given in Equation (1).
(1) In this study, a recuperative cycle is presented. Once the components are selected, the first calculation carried out by the software is to design the plant and its main components using design point input parameters. After the design procedure is completed, the program evaluates the off-design performance of the components based on the desired operational strategy. The off-design procedure starts by assuming a first attempt mass flow rate and iterates until all the components match each other. During iteration the matching guess is continuously updated until component matching is achieved. Further constraints can be included in the software according to the operational strategy. In this case another iteration is required and the system is solved via serial nested loops. The loops matching guesses and constraints are solved in a nested sequence until convergence.

Fluid Properties Module
The fluid properties are characterised using a specific class, which can be accessed by any components' classes. Through this module, each working fluid can be described as an object with its composition and thermodynamic properties. The specific heat of the working fluid C P as a function of temperature T is evaluated using the first order polynomial regression given in Equation (1).
where A and B are two coefficients characteristic of the single component k and µ is the mass fraction of the component k. N is the total number of the constituents of the fluid. To evaluate thermal conductivity and viscosity, for the purpose of this study, the class was calibrated to consider dry air as the working fluid, using a second order interpolation of the known properties of air as a function of the temperature [11]. The effect of pressure is negligible and hence was not considered.

Microturbine
The MGT is modelled by its three main components: the compressor, the recuperator and the turbine. Typically, for MGT applications, centrifugal compressors are adopted. This is mainly due to their robustness for the range of mass flow rates and pressures in which MGTs operate. Additionally, a previous study conducted by the authors demonstrated that an optimised design that utilises a vaneless diffuser can lead to a wider range of operation with relatively high efficiency compared to vaned diffuser designs [12]. Hence, a centrifugal compressor design with stationary casing, impeller, stationary vaneless diffuser and volute has been adopted for this study. The design and off-design performance evaluation have been performed using a one-dimensional model based on mass, momentum and energy conservation principles as given in Equations (2) and (3).
In Equation (2), R is the radius, b is the blade width, ρ is the density and C a is the axial component of the absolute velocity. The subscripts t and h refer to the inlet hub and tip. In Equation (3), T 0 and p 0 are the total temperature and pressure, η pc is the polytrophic compressor efficiency and γ the specific heats ratio. In both equations, the suffixes in and out refer to the inlet and outlet, respectively.
The efficiency of the compressor was evaluated using a set of standard empirical relationships available in [13]. In the adopted loss model, two categories of losses have been modeled: external (parasitic) losses and internal losses. Internal losses include all the losses that cause a pressure drop across the main channel, such as incidence losses (∆h i ), blade loading (∆h bl ), skin friction (∆h s f ), mixing losses (∆h ml ), tip clearance losses (∆h cl ) and vaneless diffuser losses (∆h cd ). The external losses are all those losses characterised by an enthalpy loss but not a pressure loss across the main channel. The disk friction losses (∆h d f ), the leakage losses (∆h lk ) and the recirculation losses (∆h rc ) belong to this category. The isentropic compressor efficiency then can be calculated using Equation (4), where ∆h s is the isentropic compression enthalpy change.
The recuperator chosen for this study is a cross flow plate-fin heat exchanger with off-set strip fins. This arrangement has one of the highest heat transfer to volume ratios of all heat exchangers available in the market [11]. Design and off-design performances of the heat exchanger have been evaluated using the ε-NTU method which provides a model for the effectiveness and the overall heat transfer coefficient as given in Equation (5). where, ε is the recuperator effectiveness, and C min and C max are the minimum and maximum heat capacity rates between the hot and cold fluids defined as the product of the mass flow rate and the heat capacity at constant pressure. NTU is the number of transfer units defined as U A/C min and χ is the ratio between C min and C max . U and A are the overall heat transfer coefficient and the total heat transfer rate, respectively. The heat transfer mode is mainly convective in the recuperator and the heat transfer coefficients (h) are predicted using Chilton-Colburn J-factor analogy [11] as given in Equation (6) as a function of the Colburn number, (J), which is a function of the Reynolds number and the geometry.
where G is the volumetric flow rate, C p is the specific heat at constant pressure and Pr is the Prandtl number. Pressure losses due to friction have been evaluated using the semi-empirical correlations proposed by Manglik et al. [14] which relates the Colburn number and the friction factor to the Reynolds number and the geometry of the channels of the heat exchanger. The considered expander is a radial inflow turbine. Using the same approach adopted for the compressor, the design performance evaluation has been performed using a one-dimensional model based on conservation of mass, momentum and energy for the flow in the turbine.
The efficiency of the turbine has been evaluated using a set of empirical relationships that can be found in the literature [15]. The losses considered in this model are: incidence loss, passage loss, skin friction, blade loading, mixing losses, losses due to secondary flows, tip clearance loss, windage loss, and scroll loss. Incidence loss is caused by the deviation of the in-flow angle from the blade optimum design angle. Passage loss includes all the losses occurring in the main impeller passage. Tip clearance loss is a consequence of leakage from the pressure to suction sides occurring at the tip of the blade and can be evaluated by modelling the leakage flow and considering a simplified geometry of the shroud with a clear separation between axial and radial components. Windage loss is caused by the flow leaked onto the back face of the turbine disk and nozzle and scroll loss is due to the skin friction occurring in the scroll. Considering all the losses discussed so far, it is possible to evaluate the turbine efficiency using Equation (9), where ∆h s is the isentropic expansion enthalpy change.
For off-design calculations, a relationship between the mass flow rate and pressure ratio is required. The mass flow rate can be expressed as a function of the outlet static pressure, rotational speed and inlet temperature. The relationship between these three parameters can be represented in a Cartesian coordinate system with the law of the ellipse or Stodola cone. This relationship has been widely applied for steam turbines and is acceptably accurate for a high number of stages, but is not particularly suitable for radial in-flow turbine application. Moreover, this method does not consider the dependency on the rotational speed [16]. For the purpose of this study, the model was modified in order to take into account the dependency on the rotational speed and the turbine polytrophic efficiency. where ϑ is a function of the rotational speed, µ is the corrected mass flow, m is the expansion polytrophic exponent of the thermodynamic transformation and ϑ is defined as:

Receiver and Dish
The receiver considered for the purpose of this study is a cylindrical air tube cavity receiver (Qiu, 2015). As shown in Figure 3, the receiver is composed of a cylindrical cavity surrounded by insulation. At the bottom of the cavity an optical splitter is used to ensure effective distribution of the solar beam over the inner surface of the receiver and, at the top of the cavity, a quartz glass is presented to minimise the convective heat loss between the receiver and the ambient. The performance of an air tube cavity receiver is strongly influenced by the losses due to natural convection and re-radiation to the surroundings [17] with predominant natural convection losses. The losses caused by the wind forced convection are negligible compared to the natural convection losses. The losses, and thus the efficiency of the receiver can be evaluated using Equations (12)- (15).
where the adduction heat transfer coefficient h conv is calculated though the Nusselt number described by Equation (13).
where, G r is the Grashoff number, T w is the wall temperature, T amb is the ambient temperature and θ i is the inclination angle of the receiver to the vertical direction. Properties of the fluid are calculated at the film temperature, the arithmetic average between the ambient and the wall temperature. The other contribution to losses is due to irradiation phenomena and can be decomposed into two components: the re-irradiation, which is the characteristic of high temperature materials (Equation (14)) and the radiation absorbed by the glass and reflected by the pipe surface (Equation (15)).
where ε p is the pipe emissivity, σ b is the Stefan-Boltzmann constant, α p and α g are the pipe material and the glass absorptivity, respectively. Q sun is the heat provided by the sun for the given dish diameter and dish optical efficiency. The absorbed heat by the receiver is passed to air, assuming perfect insulation, as given by (Equation (16)). where: In the previous equations DNIis the direct normal irradiation provided by the sun expressed in [W/m 2 ], A dish is the dish area and η dish is the dish optical efficiency. In these equations, h p and A p are the pipe heat transfer coefficient and area respectively, T in and T out are the air inlet and outlet temperature. The dish optical efficiency characterises the capacity of the dish to concentrate the incident sunrays at the receiver focal spot. The main factors reducing efficiency are the surface finishing of the mirror and the deviation of the dish surface from parabolic shape which causes the beam to spread over a wider area than the receiver window.
forced convection are negligible compared to the natural convection losses. The losses, and thus the efficiency of the receiver can be evaluated using Equations (12)- (15).
where the adduction heat transfer coefficient ℎ is calculated though the Nusselt number described by Equation (13). Figure 3. Breakthrough of the adopted cavity receiver.
Where, Gr is the Grashoff number, is the wall temperature, is the ambient temperature and is the inclination angle of the receiver to the vertical direction. Properties of the fluid are calculated at the film temperature, the arithmetic average between the ambient and the wall temperature. The other contribution to losses is due to irradiation phenomena and can be decomposed into two components: the re-irradiation, which is the characteristic of high temperature materials (Equation (14)) and the radiation absorbed by the glass and reflected by the pipe surface (Equation (15)).
where is the pipe emissivity, is the Stefan-Boltzmann constant, and are the pipe material and the glass absorptivity, respectively. is the heat provided by the sun for the given dish diameter and dish optical efficiency. The absorbed heat by the receiver is passed to air, assuming perfect insulation, as given by (Equation (16)).
where: The heat transfer and the pressure losses within the receiver pipe are evaluated considering a helically coiled tube. Many validated correlations to describe the heat transfer characteristics of a helicoidally coiled pipe are available in the literature. The model adopted in this study is the one proposed by Xin and Ebadian [18].

Cost Analysis
A detailed techno-economic analysis was carried out to evaluate the cost and the levelised cost of electricity of the plant in different arrangements. The analysis is based on the components' models and a series of cost functions to evaluate individual cost of each components, and other elements of cost of the system such as installation and maintenance.
Unlike large-scale gas turbines, only a few cost evaluation studies are available in literature for microturbines. For this study, the approach used by Galanti [19] has been used and the cost functions for the different components are described by the following equations where all the costs are given in euros.
where PR c and PR t are the compressor and turbine pressure ratios, p in is the inlet pressure expressed in bar, ∆p is the recuperator percentage pressure loss, . m is the mass flow rate in kg/s, ε is the recuperator effectiveness and P is the net power output in kW. The cost of the generator also includes the cost of power electronics and has been evaluated using Equation (21). The receiver cost function was interpolated from a study conducted by Pioneer Engineering [20] for a pressurised air solar receiver, adaptable for a solar dish and microturbine system. The study considered direct labour costs, direct material costs and manufacturing costs. Results from this study were interpolated and Equation (22) was obtained which expresses the cost of the receiver as a function of absorbed heat. In this equation, Q Rec is the receiver input heat expressed in watts. The dish cost by far has the biggest contribution to the plant cost [7] and it needs careful consideration for adequate prediction of the cost. In this study, the dish cost was evaluated using a regression model as described by Equation (23) based on a study conducted by [21]. The study evaluated the specific prices of collectors in the range from 67 m 2 to 230 m 2 for different production volumes. According to [7] the exponent n Dish does not vary significantly from unity for a wide range of dish diameters, and, for this reason, a linear dependency of the dish cost on the area was considered.
where A Dish is the dish surface area in [m 2 ]. The overall cost can be then calculated using Equation (24).
In the previous equation, C inst is the installation cost, which includes the cost of installation and other costs not related directly to components. On the base of the findings of a study conducted by Gavagnin et al. [7], this has been approximated as 17% of the total components' costs. Using the Chemical Engineering Plant Cost Index (CEPCI) each calculated cost was discounted to update the costs at the current year as described in (Equation (25)).
where C j is the cost at year j, CEPCI j is the Chemical Engineering Plant Cost Index for year, j, C ref is a known cost of the components for a particular year and CEPCI ref is taken for the same reference year. The levelised cost of electricity therefore can be calculated using Equation (26).
The capital recovery factor α has been evaluated considering typical values of i = 0.07 for the actual interest rate and n = 25 years for the lifetime of the plant. Operation and maintenance costs have been estimated to be 5% of the total cost as suggested by [22]. Costs of operation and maintenance has been considered as 5% of total costs, which includes capital and installation costs. Scale economy related to different production rates were not considered in this study.

Model Validation
In this section the model validation is presented. Given the lack of experimental data regarding solar dish-MGT power plant, the validation has been performed for each component. One or more validation case studies for similar designs to those investigated in this paper were considered for each component and results output from the present model have been compared with experimental data available in literature. Data about the models validation can be found in the supplementary materials provided by the authors.

Compressor Off-Design Model Validation
The compressor model has been validated against experimental data for three different compressor sizes [23]. Figure 4 shows the main geometrical characteristics of the impeller for these compressors. The results obtained from the model for the compressor efficiency and pressure ratio for different mass flow rates and rotational speeds are given in Figures 5-7. For the Eckardt impeller O, the validation demonstrates good accuracy for a wide range of operating conditions showing a maximum relative error of 2.5% for the efficiency and 3% for the pressure ratio. Moreover, the curves demonstrate comparable behaviour and the error is located only in short sections of the curves. The same observations are valid for the Eckardt impeller A. In this case the model shows a maximum relative error of 4% for the efficiency and 4.3% for the pressure ratio. Despite the higher maximum relative error, the overall accuracy of the model is higher when compared with the Eckardt impeller A, especially regarding efficiency. On the other hand, the pressure ratio shows a more accentuated slope of the curve when getting close to the choking limit. Nevertheless, the accuracy of the model was considered acceptable for the purpose of this study.
For the case of KIMM compressor geometry, the efficiency chart was not available and the comparison has been performed only on the pressure ratio. Results show better agreement with the experimental results when compared to the other two compressor geometries showing a maximum relative error of 1.5%. This demonstrates the suitability of the present model for smaller size compressors.

Recuperator Off-Design Model Validation
The model for the recuperator has been validated against available experimental data for a range of fin geometries as proposed by [24]. The fin geometries and arrangements used for the validation, namely, 1/10-27.03 and 1/8-16.00(D) fin geometries, are presented in Figure 8 and more details of the geometry are given in [24]. Figure 9 shows the modelling results of the Colburn factor (J) and friction factor (f) for both 1/10-27.03 and 1/8-16.00(D) fin geometries. The modelling for the first geometry presents a maximum relative error of 3% for both Colburn factor and friction factor. Nevertheless, there is some inaccuracy in the results of the second fin geometry, with a maximum relative error around 18% for the Colburn number and 6% for the friction factor. Given the higher accuracy of the model in predicting the performance of the 1/10-27.03 fin geometry, this layout will be adopted by the model for the recuperator design.

Recuperator Off-Design Model Validation
The model for the recuperator has been validated against available experimental data for a range of fin geometries as proposed by [24]. The fin geometries and arrangements used for the validation, namely, 1/10-27.03 and 1/8-16.00(D) fin geometries, are presented in Figure 8 and more details of the geometry are given in [24]. Figure 9 shows the modelling results of the Colburn factor (J) and friction factor (f) for both 1/10-27.03 and 1/8-16.00(D) fin geometries. The modelling for the first geometry presents a maximum relative error of 3% for both Colburn factor and friction factor. Nevertheless, there is some inaccuracy in the results of the second fin geometry, with a maximum relative error around 18% for the Colburn number and 6% for the friction factor. Given the higher accuracy of the model in predicting the performance of the 1/10-27.03 fin geometry, this layout will be adopted by the model for the recuperator design. relative error of 3% for both Colburn factor and friction factor. Nevertheless, there is some inaccuracy in the results of the second fin geometry, with a maximum relative error around 18% for the Colburn number and 6% for the friction factor. Given the higher accuracy of the model in predicting the performance of the 1/10-27.03 fin geometry, this layout will be adopted by the model for the recuperator design.

Receiver off-Design Model Validation
The receiver model was validated by comparing the results obtained from the one-dimensional model with experimental results available in literature [25]. The receiver consists of 15-turn helically coiled copper tube with inner diameter of 12 mm and thickness of 1 mm. The inner diameter and height of the cavity are 140 mm and 200 mm, respectively.
Results shown in Figure 10 demonstrate good agreement concerning outlet temperature for a wide range of operating points with a maximum relative error below 2.5%, which is considered adequate for the purpose of this study.

Receiver off-Design Model Validation
The receiver model was validated by comparing the results obtained from the one-dimensional model with experimental results available in literature [25]. The receiver consists of 15-turn helically coiled copper tube with inner diameter of 12 mm and thickness of 1 mm. The inner diameter and height of the cavity are 140 mm and 200 mm, respectively.
Results shown in Figure 10 demonstrate good agreement concerning outlet temperature for a wide range of operating points with a maximum relative error below 2.5%, which is considered adequate for the purpose of this study. Results shown in Figure 10 demonstrate good agreement concerning outlet temperature for a wide range of operating points with a maximum relative error below 2.5%, which is considered adequate for the purpose of this study.

Turbine Off-Design Model Validation
The turbine model was validated against experimental data available in literature for a similar design [26] in terms of size and expansion ratio. Table 1 presents the main geometrical parameters of the turbine adopted for validation. Results of turbine model validation demonstrate good accuracy for a wide range of operational conditions, which also covers the actual operating condition of turbines of interest for this study ( Figure 11). The model reproduces satisfactorily the behavior of the turbine, especially near choking condition, while for a given rotational speed, it seems to lose its accuracy towards lower mass flow rates. Nevertheless, the maximum relative error is low, apart from in a few working conditions where it reaches 13.1%. Although the coefficient of determination, calculated using Equation (28), decreases with the rotational speed it stays above 0.994. This ensures sufficient accuracy when predicting the corrected mass flow rate and the efficiency as functions of rotational speed and pressure ratio. Moreover, the inaccuracy of the model is associated within working points not representing typical working conditions and may not affect the accuracy of the cycle simulation.
In Equation (28), y i is the predicted value,ŷ i is the observed value and y i is the mean value of the observed data. working conditions and may not affect the accuracy of the cycle simulation.
In Equation (28), is the predicted value, � is the observed value and � is the mean value of the observed data.

Results
In this section a case study for a 7-kW e MGT is discussed and the techno-economic performance of the plant is evaluated using the presented models. For verification of the components' matching a model developed in the commercial software IPSEpro, was used. The benchmark results for a pure-solar MGT-based CSP power plant are presented below. Finally, a plant-size optimisation has been performed to identify feasible economic solutions for a given set of conditions such as local solar irradiation pattern.

Case Study
In this section a case study considering a 7-kW e pure-solar plant powered by microturbine is presented. The plant consists of a compressor, recuperator, solar receiver, turbine and a high-speed generator to convert the available mechanical energy to electricity. Preliminary design characteristics have been chosen based on the findings of the OMSoP project [6] and a previous study conducted by the authors [27]. Table 2 reports the preliminary design characteristics of the plant. The study has been performed considering data available in literature for solar irradiation in Rome for the period between 2004 and 2005. Figure 12 shows the DNI frequency diagram and the moving average of the DNI for a constant variation of 100 W/m 2 in the DNI. As can be seen in this figure, the highest average frequency levels are for the DNI are around the value of 800 W/m 2 . For this reason, this value has been chosen as design DNI while DNIs lower than 300 W/m 2 have not been investigated due to negligible power that can be produced in such solar irradiation level.
The study has been performed considering data available in literature for solar irradiation in Rome for the period between 2004 and 2005. Figure 12 shows the DNI frequency diagram and the moving average of the DNI for a constant variation of 100 W/m 2 in the DNI. As can be seen in this figure, the highest average frequency levels are for the DNI are around the value of 800 W/m 2 . For this reason, this value has been chosen as design DNI while DNIs lower than 300 W/m 2 have not been investigated due to negligible power that can be produced in such solar irradiation level.

Design Results
The results from the design procedure are presented in Table 3 based on the design requirements presented in Table 2 and the chosen DNI. The net power output has been evaluated considering mechanical and electrical efficiencies of 85% and 90%, respectively.
Using the models discussed previously, each component of the system was designed based on the performance parameters given in Table 3. The compressor is composed of an impeller with 12 blades and a vaneless diffuser of 60 mm radius. The impeller has a tip exit radius of 36 mm and width of 1.5 mm; the isentropic efficiency of the compressor is 0.80. The fin geometry 1/10-27.03 has been adopted for the recuperator and it has an overall heat transfer area of about 17 m 2 and a heat transfer to volume ratio of around 2900 m −1 . Its effectiveness was set to 0.85 and the heat transfer coefficient was evaluated at 163 W/m 2 K and 119 W/m 2 K on the hot and cold sides, respectively. The pressure drop is estimated to be 1% on the hot side and 3% on the cold side. The receiver is characterised by a cavity diameter of 19 cm and length of 33.6 cm. The overall length of the pipe, which has a diameter of 4 cm, is 3.748 m, and is formed of a helical coiled pipe with 8 turns. The efficiency of the receiver is 0.76 while the overall pressure loss is around 3%. Based on the receiver performance, the dish aperture diameter required to produce the desired electrical power of 7 kWe is 8.37 m, assuming dish optical efficiency of 0.8. The turbine impeller, with isentropic efficiency of 0.799, consists of 11 blades with an inlet radius of 38.5 mm, exit hub radius of 14.9 mm and exit tip radius of 28.6 mm. Appendix A summarises the design parameters of the plant components.

IPSEpro Simulation Set-Up
In addition to modelling the plant using the developed code, the microturbine was also modelled by IPSEpro, a commercial heat and mass balance software package for process and cycle analysis. The layout of the process units of the model is shown in Figure 1. Characteristics maps for the turbomachinery and recuperator have been generated based on the validated models of components as discussed in the preceding section and are implemented in the IPSEpro model for off-design calculations. Input parameters are rotational speed, DNI and ambient conditions, i.e., air relative humidity, ambient pressure and ambient temperature. The TIT is considered as constraint with maximum allowable limit of 1073 ± 1 K. When the limit of TIT is reached, the rotational speed, thus the air mass flow rate, and turbine outlet temperature are regulated to avoid exceeding the limit. The model outputs are pressure, temperature, mass flow rate and enthalpy downstream and upstream of the components, by which performance indicators such as power and electrical efficiency are calculated. The physical properties of gas components are calculated with polynomials derived from the JANAF (Joint Army Navy Air Force) Thermodynamic Tables [28]. In IPSEpro, all calculations were performed assuming that the gas components being ideal gases.

Off-Design Results
For a typical pure-solar application, the MGT operates at a constant TIT (at the highest allowable level) to maximise the overall efficiency of the plant. It should be noted that here, the TIT has been limited due to material limitation of the receiver rather than the turbine material. In case of OMSoP plant this limit was estimated to 1073 K. To achieve a constant TIT, it is necessary to operate at variable mass flow rate and consequently variable rotational speed. For lower DNIs it would not be possible to achieve the highest allowable TIT, which means that the MGT will operate at lower efficiency. Moreover, the power output for these operating points is relatively low, as is their occurrence frequency as shown in Figure 13. As such, a minimum output power of the plant was fixed to 2000 W and points with lower power output were not considered. It should be emphasised that considering any power below this level would neither change the annual produced electricity significantly nor the levelised cost of electricity.

Off-Design Results
For a typical pure-solar application, the MGT operates at a constant TIT (at the highest allowable level) to maximise the overall efficiency of the plant. It should be noted that here, the TIT has been limited due to material limitation of the receiver rather than the turbine material. In case of OMSoP plant this limit was estimated to 1073 K. To achieve a constant TIT, it is necessary to operate at variable mass flow rate and consequently variable rotational speed. For lower DNIs it would not be possible to achieve the highest allowable TIT, which means that the MGT will operate at lower efficiency. Moreover, the power output for these operating points is relatively low, as is their occurrence frequency as shown in Figure 13. As such, a minimum output power of the plant was fixed to 2000 W and points with lower power output were not considered. It should be emphasised that considering any power below this level would neither change the annual produced electricity significantly nor the levelised cost of electricity. The results of the present model have been compared with the results from IPSEpro where the same component characteristic maps and operational strategy were considered. This comparison was aimed to verify the component matching procedure adopted for the present model. The results in Figure 13 demonstrate good accuracy for the whole range of operation. The discrepancy in the results is due to the setup in IPSEpro where the TIT is limited not to a single value, 1073 K in this case, but The results of the present model have been compared with the results from IPSEpro where the same component characteristic maps and operational strategy were considered. This comparison was aimed to verify the component matching procedure adopted for the present model. The results in Figure 13 demonstrate good accuracy for the whole range of operation. The discrepancy in the results is due to the setup in IPSEpro where the TIT is limited not to a single value, 1073 K in this case, but to a range with a minimum tolerance of 1 K. Thus, while the desirable TIT is set to 1073 K in the present model this parameter was set to 1073 ± 1 K in the IPSEpro model.

Cost Analysis
The annual electricity generation and the levelised cost of electricity have been calculated using the components' characteristic maps and the plant overall performance parameters. The analysis has been carried out considering the DNI frequencies in Figure 13 based on hourly average values of Rome in the period between 2004 and 2005. The overall energy produced by the plant is 10,682 kWh. Figures 14 and 15 shows the breakdown of the overall energy produced along the year by each DNI level and, not surprisingly, the electricity is mostly produced during the high sun irradiance period in the range between 650 and 850 W/m 2 . to a range with a minimum tolerance of 1 K. Thus, while the desirable TIT is set to 1073 K in the present model this parameter was set to 1073 ± 1 K in the IPSEpro model.

Cost Analysis
The annual electricity generation and the levelised cost of electricity have been calculated using the components' characteristic maps and the plant overall performance parameters. The analysis has been carried out considering the DNI frequencies in Figure 13 based on hourly average values of Rome in the period between 2004 and 2005. The overall energy produced by the plant is 10,682 kWh. Figures 14 and 15 shows the breakdown of the overall energy produced along the year by each DNI level and, not surprisingly, the electricity is mostly produced during the high sun irradiance period in the range between 650 and 850 W/m 2 . The cost functions described previously have been used to estimate the cost of the plant and cost of producing electricity. Figure 15 shows the plant cost breakdown for different components of the plant. The capital cost has been estimated to be EUR 27,051 and, consequently, the specific cost of the plant, defined as the ratio between the cost of components and output power in design condition, has been estimated to be around EUR 3980/kWel. As can be seen in Figure 14

Plant Size Optimsation
The analysis discussed above was repeated for different plant design parameters to optimise the plant size in order to achieve minimum LCOE while maximising the annual electrical production. Rome was selected as location site, characterised by DNI levels and frequencies already shown in Figure 12. For each given DNI an off-design performance analysis was carried out considering constant TIT as operational strategy. The overall energy produced during the year was then The cost functions described previously have been used to estimate the cost of the plant and cost of producing electricity. Figure 15 shows the plant cost breakdown for different components of the plant. The capital cost has been estimated to be EUR 27,051 and, consequently, the specific cost of the plant, defined as the ratio between the cost of components and output power in design condition, has been estimated to be around EUR 3980/kWel. As can be seen in Figure 14, the most costly component is the dish, representing 63.5% of the overall cost, while the MGT represents 18.3% of the overall cost. Results from the LCOE analysis demonstrate a levelised cost of electricity of EUR 22.81/kWh considering a lifetime of the plant of 25 years.

Plant Size Optimsation
The analysis discussed above was repeated for different plant design parameters to optimise the plant size in order to achieve minimum LCOE while maximising the annual electrical production. Rome was selected as location site, characterised by DNI levels and frequencies already shown in Figure 12. For each given DNI an off-design performance analysis was carried out considering constant TIT as operational strategy. The overall energy produced during the year was then calculated together with the levelised cost of electricity. As discussed previously, the most costly component of the plant is the dish and its size is directly proportional to the mass flow rate for a given TIT. On the other hand, the mass flow rate is the main parameter influencing the power output of the plant and, subsequently, the annual electricity generation. For this reason, the mass flow rate was chosen as the optimisation variable, while pressure ratio and TIT were kept constant and equal to 3 and 1073 K, respectively. The design space of the optimisation variable is limited to the values 0.07-0.16 kg/s. These limits have been evaluated to keep the compressor efficiency as high as possible by limiting the compressor specific speed to its optimum value. At both ends of this range, the compressor efficiency drops significantly. Figure 16 shows the results of the optimisation. The performances of the plant considering dissimilar plant design mass flow rates have been compared in term of levelised cost of energy and annual energy produced. As shown in Figure 16a, results of the optimisation demonstrate a minimum LCOE of 21.5 [EUR/kWh] for a design point mass flow rate of about 0.11 kg/s. This corresponds to an overall cost of the plant around EUR 32,600, with a dish diameter of 9.4 m and an annual electricity production of 13,700 [kWh]. As expected, the annual energy production increases with the mass flow rate as has been shown in Figure 16b. As shown in Figure 17, the range of DNIs, in which the plant produces more than 2 kWe, is getting wider with an increase of the mass flow rate. Nevertheless, the dish size and, consequently, the overall cost is proportional to the mass flow rate; the right trade-off between annual energy production and levelised cost of electricity should be decided based on particular power demand for each application. Figure 17 highlights that for the whole range of mass flow rates the majority of the annual energy produced still occurs in the range between 650 and 850 [W/m 2 K].

Discussions and Conclusions
In this paper a techno-economic model for design optimisation and off-design performance evaluation of a pure-solar MGT-based plant was presented. The model adopted a modular approach using object-orientated programming where, for each component of the plant, a model was developed and validated against experimental data available in the literature. Results of the validation demonstrate good accuracy of the proposed models in a wide range of operating conditions.
A commercial software program, IPSEpro 7, for power cycle and process analysis, was used to verify the component matching procedure and cycle calculations for off-design conditions. The performance of a 7-kWe solar dish microturbine power plant was then analysed with the present platform and in the IPSEpro 7 environment. The components characteristic maps were generated based on one-dimensional modelling of the components and integrated into the models. These maps were also integrated in commercial software in order to validate the component matching routine.

Discussions and Conclusions
In this paper a techno-economic model for design optimisation and off-design performance evaluation of a pure-solar MGT-based plant was presented. The model adopted a modular approach using object-orientated programming where, for each component of the plant, a model was developed and validated against experimental data available in the literature. Results of the validation demonstrate good accuracy of the proposed models in a wide range of operating conditions. A commercial software program, IPSEpro 7, for power cycle and process analysis, was used to verify the component matching procedure and cycle calculations for off-design conditions. The performance of a 7-kWe solar dish microturbine power plant was then analysed with the present platform and in the IPSEpro 7 environment. The components characteristic maps were generated based on one-dimensional modelling of the components and integrated into the models. These maps were also integrated in commercial software in order to validate the component matching routine. Results of the validation demonstrate good accuracy of the model in any working condition of the plant. From an economic point of view, the plant presents levelised cost of electricity of EUR 22.06/kWh, a capital cost of EUR 27,051 and a specific cost of the plant of about EUR 3980/kW e . These results are comparable with those obtained by Gavagnin [7] considering 0.1 kg/s as design mass flow rate and 500 unit/year production rate. The author estimates a specific cost of the plant to be around EUR 3700/kWel. The cost breakdown also demonstrates comparable results. Gavagnin [7], based on real component costs provided by the OMSoP project supplier, estimates MGT representing 17% of the overall cost, the dish representing 62% and a receiver representing 4%. The present study estimates the most costly component as the dish with 63.5% of the overall cost while the MGT represents 18.3% of the overall cost and the receiver 3%. This demonstrates the applicability of the cost functions adopted in this study to estimate the overall cost of the plant. Nevertheless, comparison with previous studies [6,7], highlights the necessity of including the production rate as a further variable in the cost analysis.
A plant-size optimisation was then carried out to identify the plant size which minimises levelised cost of electricity while keeping the annual energy production as high as possible. Results of the optimisation demonstrates a minimum LCOE of 21.5 [EUR/kWh] for a design point mass flow rate of about 0.11 kg/s. This corresponds to an overall plant cost of EUR 32,600, with a dish diameter of 9.4 m and an annual electricity generation of 13,700 [kWh]. Results from optimisation are comparable with those obtained by Ghavami [6]. The author performed an optimisation considering dissimilar plant size and sets of turbine inlet temperatures and recuperator effectiveness. The study highlights minimum levelised cost of energy around EUR 17/kWh, when the TIT is increased to 950 • C and the recuperator effectiveness is improved up to 90%. A parametric analysis considering dissimilar TIT and recuperator effectiveness, with due regard for components' material limitation and cost, will be considered in future studies. Future studies are also aimed to investigate the possibility to include two new modules in the software, such as combustion chamber and thermal energy storage unit. Integrated hybridization and/or thermal storage improves the dispatchability of the plant with consequent benefit for annual power output. Nevertheless, the right trade-off between component costs must be found. Furthermore, the software can be easily adopted to generate a database, which can be used to train machine learning algorithm. These algorithms can be easily integrated with a genetic algorithm for multi-objective optimisation purposes.

Conflicts of Interest:
The authors declare no conflict of interest.    Figure A2 shows an example of the recuperator maps. The map has been calculated considering variable mass flow rate and constant inlet temperature for both fluids at the design point: Thot,in = 869 K and Tcold,in = 436 K. Similar charts can be developed considering variable cold and hot inlet temperature.  Figure A2 shows an example of the recuperator maps. The map has been calculated considering variable mass flow rate and constant inlet temperature for both fluids at the design point: T hot,in = 869 K and T cold,in = 436 K. Similar charts can be developed considering variable cold and hot inlet temperature. Figure A2 shows an example of the recuperator maps. The map has been calculated considering variable mass flow rate and constant inlet temperature for both fluids at the design point: Thot,in = 869 K and Tcold,in = 436 K. Similar charts can be developed considering variable cold and hot inlet temperature. The receiver maps have been generated by varying the DNI under variable mass flow rate at a given inlet temperature ( Figure A3). The proposed one-dimensional model was used to calculate the outlet temperature, the receiver efficiency and the pressure drop. The results suggest that the efficiency of the receiver is directly proportional to the mass flow rate, but the DNI has the most influence on the efficiency and the maximum achievable efficiency is dictated by the DNI. The calculations have been performed for ambient temperature of 298 K, inlet temperature of 800 K and inlet pressure of 308,600 Pa. The receiver maps have been generated by varying the DNI under variable mass flow rate at a given inlet temperature ( Figure A3). The proposed one-dimensional model was used to calculate the outlet temperature, the receiver efficiency and the pressure drop. The results suggest that the efficiency of the receiver is directly proportional to the mass flow rate, but the DNI has the most influence on the efficiency and the maximum achievable efficiency is dictated by the DNI. The calculations have been performed for ambient temperature of 298 K, inlet temperature of 800 K and inlet pressure of 308,600 Pa.  Figure A4 shows the turbine maps obtained for different rotational speeds, in the same operating range of speed as the compressor and, similar to the compressor maps, each speed line is limited by the efficiency drop.  Figure A4 shows the turbine maps obtained for different rotational speeds, in the same operating range of speed as the compressor and, similar to the compressor maps, each speed line is limited by the efficiency drop.

Abbreviations
(c) Figure A3. Receiver characteristic maps: variation of outlet temperature (a), efficiency (b) and pressure drop (c) with air corrected mass flow rate. Figure A4 shows the turbine maps obtained for different rotational speeds, in the same operating range of speed as the compressor and, similar to the compressor maps, each speed line is limited by the efficiency drop.