Access the full text.
Sign up today, get DeepDyve free for 14 days.
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
Solar systems based on the coupling of parabolic concentrating collectors and thermal engines (i.e. dish-Stirling systems) are among the most efficient generators of solar power currently available. This study focuses on the modelling of functioning data from a 32 kWe dish-Stirling solar plant installed at a facility test site on the University of Palermo campus, in Southern Italy. The proposed model, based on real monitored data, the energy balance of the collector and the partial load efficiency of the Stirling engine, can be used easily to simulate the annual energy production of such systems, making use of the solar radiation database, with the aim of encouraging a greater commercialisation of this technology. Introducing further simplifying assumptions based on our experimental data, the model can be linearised providing a new analytical expression of the parameters that characterise the widely used Stine empirical model. The model was calibrated against data corresponding to the collector with clean mirrors and used to predict the net electric production of the dish-Stirling accurately. A numerical method for assessing the daily level of mirror soiling without the use of direct reflectivity measures was also defined. The proposed methodology was used to evaluate the history of mirror soiling for the observation period, which shows a strong correlation with the recorded sequence of rains and dust depositions. The results of this study emphasise how desert dust transport events, frequent occurrences in parts of the Mediterranean, can have a dramatic impact on the electric power generation of dish-Stirling plants. ☼ Published version in Ap p lied Energy : http s://doi.org/10.1016/j.ap energy .2019.114378 available in http s://www.sciencedirect.com/science/article/p ii/S0306261919320653 © <2020>. This manuscrip t version is made available under the under the CC -BY-NC-ND 4.0 license http ://creativecommons.org/licenses/by -nc-nd/4.0/ This p ap er is accep ted for p ublication in Ap p lied Energy (2020). Keywords: Stirling efficiency; solar concentrator; energy analysis; numerical model; mirror soiling; 1. Introduction The rise in carbon emissions, to a large extent due to the combustion of fossil fuels to generate electricity, has been well-documented over the past few decades, with many studies indicating the need for innovation in terms of renewable energy generation and storage . Indeed, the rising global demand for electricity, which is predicted to continue, represents an enormous environmental risk, and how we choose to generate this electricity will be a key factor in determining the future energy policy of the European Union . Concentrated Solar Power (CSP) systems are among the most modern renewable energy systems currently used to harness solar energy , with near zero carbon emissions . These CSP systems operate by using a number of mirrors to concentrate the sun’s rays on a relatively small receiver. The resulting heat can be used either to produce steam to power a steam turbine for electricity generation or directly to power industrial processes . There are four CSP technologies already being used globally in the heat and power industry: parabolic trough concentrators, linear Fresnel reflectors, solar towers, and parabolic dish concentrators , . Compared to other CSP systems, parabolic dish technology has demonstrated the highest levels of efficiency in the conversion of solar to thermal energy , , ,  using a parabolic concentrator which tracks the sun and concentrates the direct component of the solar radiation on a receiver fixed at its focal point. The high-temperature heat in the receiver is transferred to a working fluid in the Stirling engine, which then produces mechanical energy , . The mechanical energy is subsequently transformed into electricity by a generator. Existing plants use Stirling cycle engines powered by air or hydrogen , , while recent prototype plants have been realised with the use of air microturbines. Some advantages of the dish-Stirling system are that it has high cycle thermal efficiency, fewer moving parts than other systems of similar dimensions, and a smooth-running, quiet engine . The greatest developments in parabolic dish systems have taken place since the 70s through continuous research aimed at optimising their components , , . This is demonstrated in a number of works dedicated to: the optimisation of the absorber temperature , , ; the influence of the solar collector design parameters on engine efficiency ; energy and exergetic analyses ,  and the multi-objective optimisation ,  of dish-Stirling systems. Moreover, the continuous advances in the design methods of these CSP systems, such those illustrated in ,  and , have resulted in the realisation of numerous new demonstration plants all over the world . Despite the continuous development of these systems, there have been several barriers which have not favoured their successful entry into the market of renewable electricity generation systems . Some of these are typical of all CSP systems, i.e. the large initial capital investments and the subsequent uncertainty about the operation and maintenance costs , . Others, are specific to dish-Stirling systems such as the difficulty in being integrated with storage systems compared to parabolic trough and linear Fresnel collectors . On the other hand, several recent studies have shown that dish- Stirling systems could have numerous areas of application , such as: micro- cogeneration , potable water production , off-grid electrification and water pumping. Other interesting applications of dish-Stirling systems can certainly be developed in the future considering that this solar technology is very suitable for hybridisation with other renewable energy sources, according to a more general trend that began a few decades ago in this technological sector . Another interesting study, shows that dish-Stirling systems have an environmental impact similar to that of photovoltaic systems, their main competitor, and that hybridisation with other renewable energy sources, such as bio-gas, could be the key factor in their commercialization . Therefore, from a more general perspective of decarbonisation of electricity generation systems in the EU, national and international incentive policies will certainly be needed to encourage greater commercial penetration of this interesting typology of CSP system . In order to implement these policies, it will be necessary to develop guidelines  which will have to be based on simple and realistic analysis tools to allow the evaluation of monthly and annual energy performance of these systems taking advantage, for example, of the hourly-based solar radiation databases usually available for each geographic location. Among several numerical models used to evaluate the energy production of dish- Stirling systems based on real performance data from operating dish-Stirling plants , , , , the Stine model and its subsequent evolutions are the most widely used , . These empirical models are based on the observation that there is usually a linear correlation between Direct Normal Irradiation (DNI) and the electrical power generated by dish-Stirling plants when the collector mirrors are clean , , . The coefficients of the correlation can be calibrated using function test data on a full day with clear skies. These constants can be further corrected to take into account the effect of air temperature changes and degradation of optical efficiency due to soiling of the collector mirrors , , . Above all, the effect of mirror soiling is a very important aspect in assessing the performance of a plant because, as experimental data show, it can lead to a significant reduction in the efficiency of the system. . In this paper we show how data collected during the monitoring program of a demo dish-Stirling plant, recently built at a facility test site at Palermo University (Sicily), was used to assess experimentally the main factors affecting the performance of these systems. The plant site is characterised by a Mediterranean climate . The tests performed on this plant are particularly interesting because there are few demonstrative dish-Stirling plants installed in the Central Mediterranean region and there is therefore no evidence of the effects of the particular micro-climate of this region on the performance and maintenance operations of this type of plant. This geographical area is, in fact, characterised by rather high levels of DNI  particularly during the warmer periods of the year. However, cloud decks associated with typical mid-latitude depressions and frontal passages are evident during cold periods of the year. In addition, desert dust transport events are frequent in this area. Dust transport occurs almost every day in different parts of the Mediterranean, especially during spring and summer , , , . Both the frequent passage of clouds in the winter season and the dust deposition on collector mirrors in the spring and summer seasons can reduce the performance of dish-Stirling systems, yet these effects are often difficult to predict and manage. Using operating data measured on days when the concentrator mirrors were clean, we have elaborated a physical-numerical model based on a simplified energy balance of the collector (according to a well-established technique for modelling the CSP system performance ), taking into account the part-load efficiency curve of the Stirling engine . Assuming further reasonable simplifying hypotheses, based on the experimental evidence, it has been possible to linearise the model obtaining a physical interpretation of the empirical coefficients of the model proposed by Stine. In addition, we compared the performance measured for the dish-Stirling system with soiled mirrors with that expected from the model with clean mirrors, for the same conditions of DNI and air temperature, using a technique proposed in the literature to assess the effect of module soiling on the energy production of photovoltaic systems . Using this method, it was possible to define an average daily cleanliness index of the system for each day during the monitoring period. In a similar way to that proposed for photovoltaic systems , we compared the evolution of this index with the history of rains and dust deposition events recorded in the same period for the area where the plant is installed. 2. Experimental setup and performance data analysis 2.1 The dish-Stirling demo plant in Palermo On 14 December 2017, the installation of a dish-Stirling plant at a facility test site at Department of Engineering of Palermo University was completed (see Fig. 1). The CSP unit, which was manufactured by the Swedish company Ripasso Energy and constructed by the Italian companies Elettrocostruzioni S.r.l and HorizonFirm S.r.l., started operating in January 2018. During the following months, the demo plant underwent an intense test and monitoring program which was intended to assess the performance of the system and its possible commercialisation for the CSP Italian market. In the preceding years, Ripasso Energy had already installed a series of similar dish-Stirling units on a test site at Upington in the Republic of South Africa, which is a location characterised by very high levels of DNI. Figure 1. dish-Stirling unit at the Department of Engineering, University of Palermo, Italy. The electric output of the Ripasso CSP unit is typically linearly proportional to the direct solar radiation with a net electric peak power of 31.5 kW at a DNI value of 960 W/m . Nevertheless, in November 2012, the Ripasso South African plants set the current world record of energy conversion efficiency from solar-to-electric at 32%, making the Ripasso dish-Stirling a state-of-the-art system in the CSP sector . The most relevant technical data on the CSP unit installed at the facility site at the University of Palermo has been summarised in Tab. 1 and is briefly described below. Table 1: Technical data of the Ripasso Energy dish-Stirling unit. Concentrator Dish diameter 12 [m] Focal length 7.45 [m] Effective aperture area 101 [m ] Geometric concentration ratio 3217 [ ˗ ] Reflectivity of clean mirrors 95 [%] Total height during operation 14 [m] Total weight 8000 [kg] Occupation of soil area 500 [m ] Power Conversion Unit Type of Stirling engine 4 cylinders double acting [ ˗ ] Displaced volume 4 × 95 [cm ] Typical output at DNI=960 W/m 31.5 @ 2300 rpm [kW ] Receiver temperature 720 [°C] Working gas Hydrogen [ ˗ ] Max gas pressure 200 [bar] Power control Pressure - speed control [ ˗ ] Weight 700 [kg] The Ripasso Energy dish-Stirling unit consists of a paraboloidal concentrator, a Power Conversion Unit (PCU) and a tracking system. The paraboloidal dish reflector, with a diameter of about 12 m and an effective aperture area of 101 m , is made up of a lightweight, steel dish structure, mirror facets and three truss beams supporting the PCU. The facets, mounted on the steel structure, consist of sandwich panels with a low-iron, thin, back-silvered glass mirror surface characterised by high levels of reflectivity (95% for clean surfaces). The high precision dual-axis tracking system, on the other hand, is composed of the tracker structure, the actuators, the control cabinet and an annular concrete foundation occupying an area of about 500 m and supporting the weight of the whole plant (about 8.7 tons). The tracker steel structure, in turn, is made up of two main parts: a carousel-type structure mounted on four wheels for azimuthal rotations and a cradle-type structure for elevation rotations. Finally, the PCU consists of a heat receiver, with an internal cavity, a Stirling engine, an alternator and a cooling system. The tracking system, continuously following the daily path of the sun, allows the parabolic dish to collect, reflect off and concentrate the direct solar radiation onto the receiver of the PCU located at paraboloid’s focal point (focal length of 7.45 m). The high precision tracking and optics systems permit a uniform distribution of the thermal energy on the receiver surface, which is adjacent to the hot side of the Stirling engine. Thus, the thermal energy absorbed by the receiver heats the working gas of the Stirling engine (hydrogen) to a temperature of about 720 °C. Part of this heat is firstly converted into mechanical energy by the engine and then transformed into three-phase electricity by the alternator connected to the crankshaft of the engine. The remaining heat is, instead, released into the environment by the fan cooling system mounted on the rear of the tracking structure and connected by a pipe and a heat exchanger to the coldest part of the Stirling engine. The cooling system, where water circulates, allows the maintenance of a high temperature difference between the hot and cold sides of the Stirling engine (720 °C and 70 °C, respectively). This is fundamental to guarantee high conversion efficiency from thermal to mechanical energy. Although the Ripasso CSP unit is characterised by numerous innovative and highly efficient technological components, especially the tracking system and the receiver, the core part of the whole system remains the Stirling engine. The highly efficient Ripasso Stirling engine is an updated version of the original model (USAB 4-95) licensed from Kockums AB in 2008. Some data concerning the performance of the original USAB 4-95 Stirling engine, working with a temperature of 720 °C, can be found in . In the past, this engine was installed on different dish-Stirling models, among which are: the MDAC systems located at different test facilities in the USA (1984-1985)  and the SES MPP systems (2008), one of which is located at the Sandia National Laboratories test site in New Mexico, USA . These plants represent precursors to the evolved versions that Ripasso Energy has realised in recent years at the Ripasso test site at Upington, South Africa and at the facility site of the Department of Engineering of Palermo University, Italy. 2.2 Measurement devices and data processing In the period between 11 January and 2 July 2018, an extensive measurement campaign was carried out to collect data related to the operation of the CSP plant located at the facility test site in Palermo. During this period, although the full operation of the solar plant was interrupted several times due to the maintenance and testing stages planned for the demo plant, it was possible to record, for 3,960 hours, the performance of the system under different conditions of solar irradiance, air temperature and cleanliness of the mirrors. The main physical quantities recorded during this stage were: the net electrical power output Ė using an energy meter (ABB - A43 212-100, with an accuracy of the active and reactive energy class of ±1% and ±2% respectively), the solar beam radiation I using a pyrheliometer (Kipp&Zonen – SHP, with an average accuracy of ± 0.5%) and the external ambient temperature T using a weather transmitter (Vaisala - WXT536 with an accuracy air of ±0.3°C at +20°C). The other most relevant physical quantities measured in the same period were the parasitic electric power Ė and the limiting temperatures of the Stirling engine T and T . h c The original measured quantities were collected with a one-second temporal resolution in a MySQL database consisting of 14,256,000 records. The data was preliminarily filtered to avoid points representing mechanical transients. Thus, considering that the sampling frequency f of the signals coming from the dish-Stirling plant is 1 Hz, it was first necessary to define a “moving average filter” procedure useful for both the reduction of the quantity of data and the filtering of the transients. The moving average is a low-pass filter, very common for regulating an array of sampled data, and is calculated according to the following formula : N 1 y i x j (1) j 0 xj yi where is the monitored input signal, is the resulted signal, and N is the considered number of points. In the present case study, i is the index referring to minutes, j is the index referring to the seconds of the considered minute, and N is equal to 60. This kind of filter is optimal for many common problems since it reduces random white noise while keeping the sharpest step response. An example of a moving average filter application is shown in Fig. 2, where the previously described procedure has been applied to a squared signal affected by random white noise. Fig. 2: comparison between an original square signal and a filtered signal expressed by Eqn. (1) over a time series of 300 seconds with 1 Hz of sampling frequency. On average, the amount of noise reduction is equal to the square root of the number of points. In the present case study, a 60 point moving average filter reduces the noise by a factor of 7.7. Thus, after defining a cut-off frequency f for which the attenuation of the output is larger than 3 dB, this frequency is equal to : 0.442947 0.443 f f f 0.00738 Hz (2) c s s N 1 Furthermore, all the data for which the revolutions per minute of the Stirling engine are over 2020 has been filtered out because this is indicative of a transient condition. As a result of the filtering procedures, it was possible to select the data from phases where the system worked continuously and in stationary conditions. This data was used for the analysis presented in this study. 2.3 Analysis and discussion of the experimental results Preliminary to the analysis, the experimental data collected during the testing campaign (occurred between 11 January 2018 and 2 July 2018) was divided into two groups classified in accordance with whether they were recorded on days when the collector mirrors were clean or not. In fact, during the period in which the experimental data was recorded, the dish-Stirling plant was subject to numerous rainfall cycles of different intensity and sirocco events with associated variable levels of dust deposition. For this reason, we assumed that most of the analysed records correspond to conditions in which the collector mirrors were naturally characterised by different (unknown) levels of soiling. On the other hand, (although no reflectometric measurements were carried out during this preliminary phase of the experiment) a total of 13 days were selected from the total data set in which the mirrors were classified as clean also in relation to the special conditions occurring in the preceding days: the first 4 days between 18 January and 2 February immediately after removing the protective films from the mirrors and in a period in which the frequent rains kept them clean; another day on 30 May after a mirror washing had taken place on 23 May; and finally 8 other days between 16 June and 2 July when the mirrors were particularly clean as a result of heavy rains occurring on 14-15 June. On most of the other days, with the exception of these 13, the effect of dust deposition on the surface of mirrors was often observed even by the naked eye (by comparing the mirrors mounted on the collector with a perfectly clean sample mirror), such as on 6 March, as a result of a significant Saharan dust transport event (Sirocco) occurring from 28 February to 5 March. With this in mind, we decided to follow the approach suggested by Sandia National Laboratories for the analysis of the performance of dish-Stirling plants . The data of net electric power, corrected for the effect of the ambient temperature, was correlated with the corresponding values of direct solar irradiation and plotted in Fig. 3. In relation to the correction factor for the air temperature, it should be noted that in the first phase of the testing stage (between January and February) the average outdoor temperature ranged between 15-16 °C while in the last period (between May and July) between 26-28 °C. Fig. 3: Corrected net electric output power  versus DNI for the Palermo dish- Stirling unit (the records have been categorised corresponding to clean or soiled mirrors). The data of air temperature, DNI and electrical power output for a day when the mirrors were clean (20 June 2018) are plotted in Fig. 4 as a function of time. The same data is also highlighted with yellow triangles in Fig. 3. On this day the air temperature was almost constant at about 27 °C, the DNI had a wide variation between 320-800 W/m , while the corresponding net electrical output power ranged between 4 kW and 22 kW. This day, as shown in Fig. 4, was particularly characterised by scattered variations of the DNI levels since a first slow passage of clouds occurred between 12:00 PM and 12:45 PM and a second rapid sequence of cloud transitions started at around 2:00 PM and concluded with a thunderstorm at about 3:00 PM. Fig. 4: A comparison of DNI and net electric output power as a function of time, alongside variation of air temperature, measured on 20 June 2018. Thus, from the analysis of the resulting data plotted in Fig. 3 and Fig. 4 it is possible to observe that: The selected data corresponding to clean mirrors is well-correlated by a linear relation (R =0.96) representing the upper limit of the experimental point cloud all the other data shows that the plant performance was lower compared to that of the first group with similar values of DNI the data recorded on the 6 March (the day when the mirrors were visibly soiled after a Sirocco event) represents the lower limit of the experimental point cloud the data recorded on 20 June is well-aligned with that of clean mirror conditions despite the scattered variations in the DNI levels due to the fast passage of clouds on this day. Therefore, considering both the methods used to select the clean mirror data and the previous observations, it can be deduced that a large majority of the reported production losses during the testing campaign were caused by the soiling of the mirrors and only marginally by other factors such as: sun tracking errors, fast passage of clouds, engine transients, etc. With more specific regard to the clean mirror performance of the dish-Stirling plant in Palermo, the correct net electrical output power data was correlated with the corresponding DNI value s by means of a linear relationship, as suggested by the present experimental data and in accordance with the model proposed by Sandia National Laboratories for the SES system . The equation in question has been reported in Fig. 3. Using this equation with I =960 W/m it is possible to extrapolate a net electric output power equal to =28 kW, which added to the average value of the parasitic electric consumption recorded during the operation of the Palermo system ( =3.6 kW), makes it possible to calculate a gross electric output power equal to =31.6 kW. The latter value is very close to the rated gross electrical power output ( =32.6 kW) that can be calculated from the data declared by the manufacturer (see Tab. 1) for the dish-Stirling system with E E clean mirrors ( =31.5 kW and =1.1 kW). This consistency confirms that the peak performance of the plant is in line with what was expected for this particular system and that the selection of the set of data corresponding to the condition of clean mirrors was apt. On the other hand, the discrepancy between the parasitic electric power consumption values indicates that it will be necessary to work further on the optimisation of the electric consumption of both the sun tracking and cooling systems of the demo plant in Palermo. However, as already illustrated in Fig. 3, the effect of the mirror soiling can produce much greater reductions in performance than those due to parasitic consumptions, especially at higher levels of DNI. In order to quantify these effects on the optical efficiency of the system, we used the method proposed by Stine  through the introduction of a reductive factor to the slope of the linear function of Fig. 3. Therefore, by fitting with a linear expression the data of as a function of DNI relative to the measurements from 6 March, the day when the reflectivity of the mirrors reached its minimum, it was possible to deduce a reduction in the optical efficiency of the collector equal to about 30%, calculated by comparing the slope of the linear expression corresponding to clean mirrors with that of soiled mirrors (see Fig. 3). Finally, Fig. 5 reports the data of the receiver temperature as a function of the DNI variations all of which were recorded during the experimental campaign (with cleaned and soiled collector mirrors). Thus, analysing the complete data set, it is possible to observe that the data varies between a minimum of about T = 710 °C and a maximum of T =728 °C with r r an average of T =719.79 °C. However, in about 85% of the data the temperature is within the narrow range of T = 720±0.5 °C. The remaining values represent conditions in which there have been small transients caused, for example, by the start-up of the system or the rapid passage of cloudy formations, etc. These variations with respect to the average, as can be seen from Fig. 5, are independent of the DNI values. Thus, observing this data it can be deduced that, regardless of the conditions of irradiation and ambient temperature, the receiver temperature remained almost constant at about T =720 °C. This observation, as will be clear later on in the discussion presented in this article, will be fundamental for the development of the numerical model presented in this work. Fig. 5: Temperature of the receiver as a function of the variations of DNI recorded during the testing campaign. The points represent data when the collector mirrors have different levels of soiling. 3. Theoretical analysis 3.1 Energy balance of the dish-Stirling system In dish-Stirling systems (see Fig. 6), the mirrors constituting the primary optics of the collector (usually with a parabolic shape) concentrate the sun’s light on the focal point of the collector, where a cavity-type receiver absorbs and transfers the concentrated solar energy to a Stirling engine driving an electric generator. Unlike other types of solar collectors (photovoltaic and thermal), CSP systems, such as solar dish-Stirling plants, can only transform the direct component of solar radiation into thermal/electric power. Therefore, the rate of solar energy reflected off the mirrors of a collector continuously tracking the sun can be calculated as , , : Q IA (3) sun b n where I is the solar beam radiation and A is the net effective surface of the mirrors b n (projected mirror area) of the dish collector. The latter quantity can be calculated from the total collector aperture area A by: AA (4) n m a where is a factor taking into account all of the causes of the losses of effective reflective area, such as blocking and shading, and the gap between the mirror facets . However, not all of this solar power is effectively reflected off the mirrors, intercepted by the cavity receiver aperture and then absorbed by its surface, because of different factors that affect the optical performances of the collector, the main ones being : reflectance of the mirrors mirror surface slope errors soiling of the mirrors sun tracking system errors receiver positioning errors (often combined with mirror errors) wind-induced vibrations effective absorbance of the cavity receiver water condensation at the mirror surface during the early hours of the day Figure 6. Schematic diagram showing the energy balance of the solar collector. Taking into account a combination of these effects, the concentrated solar power effectively absorbed by the cavity receiver is usually expressed as: QQ (5) r ,in sun o where is the optical efficiency of the concentrator system, defined as: (6) o cle where ρ is the clean mirror reflectance (which also incorporates the losses due to slope mirror errors), is a factor which takes into account the mirror cleanliness, γ is the cle intercept factor of the collector and is the effective cavity receiver absorptance . In other words, optical efficiency can be interpreted as the ratio between solar power absorbed by the receiver and the rate of solar power hitting the collector, while r ,in the intercept factor γ is the ratio between the rate of energy hitting the receiver and the sun effective concentrated solar power reflected off the mirrors of the concentrator collector. The intercept factor γ, in particular, incorporates the combination of errors associated with the control of sun tracking, the erroneous positioning of the receiver and defocusing, due to wind-induced vibrations of the collector structures. It is also necessary to consider the additional convective and radiative thermal losses from the aperture of the cavity receiver, due to the high temperatures reached at its internal surface. Then the total power losses from the receiver can be defined as the following sum: Q Q Q (7) r ,out con rad Q Q where and are the convective and radiative losses respectively. The con rad convective losses, can be expressed as , , : Q h A T T (8) con r r r air where h is the effective convective coefficient of the cavity receiver, A is the r r surface area of the aperture of the receiver and T and T are the internal receiver and r air external ambient temperatures, respectively. The second component of thermal losses, the radiative losses, are generally modelled as , : (9) Q A T T rad r r sky where is the effective emissivity of the cavity receiver , , is the Stefan–Boltzmann constant and T is the effective sky temperature. This temperature sky depends on the atmospheric conditions and is usually related to the air temperature by empirical formulas like the following , , : 1.5 TT 0.0552 (10) sky air It is not trivial to observe that the greater the aperture of the cavity receiver A , the greater the thermal losses and the higher the intercept factor . Thus, the correct r ,out design of the dimension of the receiver aperture should always be based on the trade -off between two opposing tendencies : the reduction of thermal loss from the receiver and the increase in optical efficiency, the latter being proportional to the intercept factor as shown in Eq. 6. From the energy balance of the receiver, it is possible to express the thermal power which is effectively delivered to the Stirling engine by the following difference , , : Q Q Q (11) S ,in r ,in r ,out which, using all of the above equations, can be recast as , , : (12) Q I A A h T T T T S ,in b n o r r r air r sky Thus, dividing Eq. 12 by Eq. 3, it is possible to make explicit the expression describing the thermal efficiency of the solar collector as , , : S ,in 44 (13) h T T T T c o r r air r sky Q IC bg sun where C is the solar collector geometric concentration ratio defined as: C . (14) After calculating the thermal solar input power to the CSP system, the amount of this power that can be transformed into electrical power needs to be determined, but it is first necessary to evaluate the mechanical output power of the Stirling engine by exploiting the energy balance of the PCU (see Fig. 7). Figure 7. Schematic diagram showing the energy balance of the PCU. The mechanical output power of the Stirling engine is related to the input thermal power by the efficiency of the engine itself as: WQ (15) S S S ,in From a theoretical point of view, the Stirling cycle consists of four internally reversible processes: two processes are isothermal compressions at temperature T and T ; h c the other two thermodynamic processes are isochoric (constant volume). During isothermal expansion at temperature T , heat is added to the working fluid, while during isothermal compression at temperature T , heat is rejected to the surroundings. A regenerator with full effectiveness allows the heat rejected during the isochoric rarefaction to be transferred during the isochoric compression. This Stirling cycle exhibits the same Carnot thermal efficiency. The efficiency of the Stirling engine, in turn, could be generally related to the reversible Carnot efficiency of the engine , between the same T and T limiting h c SC, temperatures, as , , : (16) S S,, ex S C where is the exergetic efficiency of the engine. The Carnot cycle efficiency, S, ex which is a function of both the heat input temperature T and the heat rejection temperature T , , : (17) SC , represents the ratio between the maximum theoretical mechanic power output of the ideal Stirling cycle operating between the limiting temperatures T and T , and the actual h c thermal power input The exergetic efficiency , represents the ratio between the S ,in S, ex actual mechanical power output and the maximum mechanical power. The value of the exergetic efficiency takes into account the irreversibilities of the engine  and is found between = 0.55-0.88 , . More generally, it is also necessary to consider that the S, ex performance of these kinds of systems depends on many factors , among which is the effect of the partialisation of the input thermal power due to fluctuations in solar radiation . For solar dish-Stirling engines, it is often assumed = 0.5 , , though there S, ex are relatively few experimental works on the performance of these engines. In addition, using the energy balance of the Stirling engine (see Fig. 7), based on the assumption that the heat delivered to the Stirling engine and the heat dissipated occur S ,out simultaneously , it is also possible to express the thermal power that must be dissipated by the cooling systems as: Q Q W 1 Q (18) S ,out S ,in S s S ,in Based both on the analysis of the real data of partial load operation of the original USAB 4-95 Stirling engine (working with T = 720 °C)  and on the performance of the Stirling engine of the SES plant , we decided to approximate with a linear regression the curve between the thermal input power and the mechanical output power of the engine as: W () a Q a R (19) S 1 S ,in 2 T where a and a are the two positive fitting constants and R is the ratio of a reference 1 2 T temperature T to the current ambient temperature T (both expressed as Kelvin degrees): o air (20) air This last correction factor was introduced, inspired by what was proposed in the Stine model ,  and similar ones , to take into account the effect of variations in outside temperature on the efficiency of the Stirling engine. Finally, combining Eq. 15 with Eq. 19, it follows: aR (21) ST 1 S,in which is an expression of the efficiency of the Stirling engine as a function of the partialisation of thermal input power and the ambient temperature. Using the mechanical output power of the engine, it is then possible to evaluate the gross electric power of PCU by: EW (22) g e S where is the electrical efficiency of the alternator. Finally, the net electric power of the e CSP system can be determined by subtracting the total parasitic electric absorption from the electric gross power: E E E (23) n gp The total parasitic power of the dish-Stirling unit, in turn, is mostly the sum of two contributors: the electric absorption of the control system and tracking motors and pt , those of the pumps and fans of the cooling system , i.e.: pd , E E E (24) p p ,t p ,d Finally, the instantaneous efficiency from solar-to-electric for all the dish-Stirling plants can be calculated as the ratio between the net electrical power produced and the solar power hitting the concentrator mirrors : . (25) DS sun By substituting: the right-hand-side of Eq. 12 into Eq. 19 and, then the right-hand- side of the obtained expression into Eq. 22 and, finally, the resulting expression into Eq. 23, it is possible to obtain, in a compact form, the following equation: (26) E a A R I a Q a R E n e o 1 n T b e 1 r ,out 2 T p where are the thermal losses from the receiver described by Eqs. 7-10. The model r ,out described by Eq. 26, once calibrated with the real operating data, can be applied in order to predict the net electric power output of the dish-Stirling system taking into account: the changes of the direct normal irradiation and external temperature the optical efficiency (by the parameters of Eq. 6) the effect of mirror soiling (by the cleanliness index in Eq. 6) the thermal losses from the receiver (by the parameters of Eqs. 7-10) assuming the receiver temperature is constant during the operation of the plant the expression of the thermal efficiency of the Stirling engine as a function of the partialisation of the thermal input loads (by Eq. 19) and external air temperature (by Eq. 20) the parasitic powers of the tracking and cooling systems of the PCU (by Eq. 24) In Fig. 8, the calculation flowchart describing the calibration procedure of the mode l is presented. The procedure is essentially aimed at the estimation of the two slope a and W Q the intercept a of the linear relationship between and . This algorithm must be S S ,in performed using the dataset corresponding to the condition of clean mirrors of the collector as input. Therefore, in the described procedure it has been assumed that the coefficient cle is always unitary. Finally, Fig. 9, represents the flowchart describing the calculus of the net electric output power from the plant using the model with its parameters already calibrated. Figure 8. Model calibration procedure flowchart using a dataset corresponding to the condition of clean mirrors of the collector. ☼ Published version in Ap p lied Energy : http s://doi.org/10.1016/j.ap energy .2019.114378 available in http s://www.sciencedirect.com/science/article/p ii/S0306261919320653 © <2020>. This manuscrip t version is made available under the under the CC-BY-NC-ND 4.0 license http ://creativecommons.org/licenses/by -nc-nd/4.0/ This p ap er is accep ted for p ublication in Ap p lied Energy (2020). Figure 9. Calculation flowcharts for the determination of the net electric output power. 3.2 A linear model of dish-Stirling electric power generation From the inspection of Eqs. 7-10, using an appropriate set of parameters, it may easily be verified that, if the receiver temperature T is kept almost constant during plant operation (as, for example, observed in the measurements dataset of the dish-Stirling in Palermo), the thermal losses from the receiver are not very sensitive to variations in T . Under this air assumption, it follows from Eq. 26, that the gross electric power output can be expressed as: ave E a A I a Q a R (27) g e o 1 n b e 1 r ,out 2 T ☼ Published version in Ap p lied Energy : http s://doi.org/10.1016/j.ap energy .2019.114378 available in http s://www.sciencedirect.com/science/article/p ii/S0306261919320653 © <2020>. This manuscrip t version is made available under the under the CC -BY-NC-ND 4.0 license http ://creativecommons.org/licenses/by -nc-nd/4.0/ This p ap er is accep ted for p ublication in Ap p lied Energy (2020). ave where is an average value of the thermal losses from the receiver. Considering, r ,out moreover, a constant value of the cleanliness index of the mirrors ( constant from Eq. 6) o and an average value of the air temperature (R constant from Eq. 20), it can be deduced that Eq. 27 simplifies to a linear relationship between the gross electrical power output and the solar beam radiation. Interestingly, this linearity emerges from the analysis of the operational data of existing dish-Stirling plants and is the basis of the empirical model proposed and validated by Stine , . As in the Stine model, Eq. 27 shows that the optical efficiency reduction due to the soiling of mirrors has the practical effect of changing the slope of the linear equation , . Stine in particular introduces this effect through the corrected insolation concept defined as the product of the current beam radiation with a reflectance ratio, the latter defined as the ratio between the measured current reflectance (with soiled mirrors) and that with clean mirrors . On the other hand, unlike the original Stine model, we have defined the correction term R in Eq. 26 in analogy to what is proposed in the Sandia National Laboratories model for SES dish-Stirling plants , i.e., as a ratio between a rated temperature and the current air temperature (see Eq. 20). Moreover, if we further assume a constant value for the parasitic absorptions, equal ave E to the average of the values measured during the operation of the plant , we can obtain from Eq. 27: ave ave (28) E a A R I a Q a R E n e o 1 n T b e 1 r ,out 2 T p which represents a linear expression between the net output power generated by the plant and the normal direct radiation hitting the collector mirrors. Also this linearity has been extensively observed by the experimental performance data of numerous dish-Stirling plants operating with clean mirrors , ,  and is further confirmed by the results from the Palermo plant which are presented in this work. Thus, the analytical expressions represented by Eq. 27 and Eq. 28 can be considered to be the main result of this work as they provide a theoretical interpretation of the empirical models presented by Stine for the SBP dish-Stirling and by Sandia National Laboratories for the SES dish-Stirling. The latter models, calibrated on a limited number of tests carried out on pilot plants, have been used successfully to extrapolate the energy production and economic performance of dish-Stirling plants for potential commercial sites. The main advantage to using the new model presented in this paper, is that it can also be used to carry out technical-economic optimisation analyses of dish-Stirling systems such as those presented, for example, by Sandia National Laboratories for SES plants . Another possible application of the proposed approach is to re-analyse the dish-Stirling performance data already available in the literature , with the purpose of carrying out a comparative analysis of the performance of these systems. For example, it is possible to extrapolate the efficiency curves of the engines with which they are equipped using the system waterfall charts and data on net electric power output as a function of DNI. To this end, we can further simplify Eq. 28 assuming an average value of the ave ave T R ambient temperature and, therefore, an average value of the corrector factor , air T obtaining: E b I b (29) nb12 where the slope is equal to ave 3 b a A R 10 (30) 11 e o n T and the intercept is equal to ave ave ave b a Q a R E (31) 2 e 1 r ,out 2 T p In Eq. 30 and Eq.31 we have considered b expressed in [kW·m /W] and b , a , 1 2 2 ave ave Q E and in [kW]. The parameters b and b in Eq. 29 can be estimated through a r ,out p 1 2 linear regression of the data of as a function of the DNI. Then, using the waterfall charts (elaborated for the peak performance of the collector), it is possible to deduce: the reflectance of the mirror ρ, the intercept factor of the collector γ and the peak value of the max thermal input power to the Stirling engine . Thus, by estimating a realistic value of the S ,in effective absorbance of the receiver α, it is then possible to calculate the average value of the thermal losses from the receiver from the following expression: ave max max Q I Q (32) r ,out b S,in max where I indicates the maximum design input value of DNI. Finally, from Eq. 30, it is possible to calculate: a 10 (33) ave AR e o n T and, after introducing Eq. 33 in Eq. 31: ave 3 ave b E A 10 b Q 2 p o n 1 r ,out a (34) ave AR e o n T where a and a are the slope and the intercept of the Eq. 19, respectively. 1 2 3.3 A method for estimating the daily averaged cleanliness index of mirrors Every solar concentrating device like dish-Stirling system is affected significantly by the presence of dust in two ways: one is the reduction of incoming solar radiation, and the other is the deposition of dust particles on the surface of the reflectors. The accumulation of soiling on the mirrors reduces the reflectivity and can result in an uneven energy flux on the receiver . In other words, progressive soiling leads to an increasing reduction in the optical efficiency of the collector and, therefore, to a decrease in the electrical power generated by the system. When production losses are no longer economically acceptable, it is, therefore, necessary to wash the mirrors. The frequency of washing cycles, however, is clearly dependent on the particular climatic conditions of the site where the plant is installed. In the Mediterranean region, in particular, dust greatly influences radiative transfer. The radiation can be attenuated by a few tens of W/m , and during severe episodes even by a few hundred W/m . The deposition of dust is quite significant. It has been estimated, with experimental and modelling tools, that approximately 108 metric tons of dust are deposited annually over Mediterranean waters, and almost the same amount is transferred towards the landmass of Europe , . The main source of material for soiling are the various types of aerosols in the atmosphere. The aerosols are of natural and/or anthropogenic origin , , . Anthropogenic particles are regularly released by human activities (mainly combustion) or the by-product of conversion of gaseous pollutants to aerosols mainly sulphates and nitrates , , . They are of submicron size and are highly hygroscopic. Deposition of these aerosol particles on the mirrors occurs mainly during high moisture conditions (e.g. during night hours or when drizzling). The main source of natural particles in the Mediterranean region is Saharan dust. Saharan dust transport exhibits seasonal variability with a peak during springtime , , . The dust particle size is from submicron to about 50 μm. The transport is of an episodic type, with each episode lasting from one to three days regularly. These can be predicted relatively easily by utilizing Numerical Weather Prediction (NWP) models with a coupled dust cycle sub-model (e.g. SKIRON/Dust system, ). Such systems have been developed for the Mediterranean region in the 2000s and are regularly in operation by providing 3-5 day predictions of dust concentration and deposition (dry and wet). Other sources of naturally-produced aerosols are sea-salt and pollen from trees. Sea salt is the remnant of seawater droplet evaporation (wave breaking, bubble generation and breaking). Sea salt is highly hygroscopic and sticks to the mirror surfaces relatively easily. Sea-salt regularly affects coastal regions . Pollen is a more hydrophobic type of particle that does not stick easily to the mirrors. The assessment of energy production losses due to soiling of the collector mirrors is an important part of the performance monitoring of a CSP plant. The soiling level is generally estimated by periodically measuring the average reflectivity of the collector mirrors . Experimental data describing the history of reflectivity for a dish-Stirling plant similar to that studied in this work shows that variations in soiling rate after mirror washing are difficult to predict . Usually, the reflectivity measurements are locally conducted at a limited number of points on the facets of the collector; these readings are then averaged. However, since the soiling distribution on the collector is often not uniform, it is difficult to assess how representative this average is. In addition, although the reflectivity can be recorded regularly, it is very difficult to carry out continuous measurements during the ongoing operation of the plant. Therefore, the previously described methodology is difficult to implement in the real-time performance monitoring system of a dish-Stirling plant. Moreover, as described above, the soiling of mirrors not only affects the reflectivity of mirrors but, more generally, reduces the global optical efficiency of the collector. For this reason, in this work, the cleaning index of the mirror has been considered, from a cle conceptual point of view, to be a more correct indicator of the global effect that soiling can induce on the optical performance of the collector compared to the reflectance. Despite the complexity of the problem, after verifying the great impact of mirror soiling on the performance of the Palermo dish-Stirling (see Fig. 3), we decided to develop a simple and practical numerical method for estimating a daily averaged cleanliness index day of mirrors . This index has been defined in analogy to what is proposed in the literature cle to study the soiling level of the modules of photovoltaic systems . The main purpose of the proposed method was to develop a numerical tool for the real-time optimisation of the frequency of washing cycles of the mirrors. The algorithm, described in detail by the calculation flowcharts in Fig. 10, is based on the simple comparison between the measured performance of the system (with soiled mirrors) and that expected using the predictive model previously calibrated using data referring to the condition of clean mirrors of the collector. A similar approach was recently proposed to assess the effect of the soiling of modules on the energy production of large scale photovoltaic plants . Finally, it is important to underline that the method proposed with the algorithm in Fig. 10 should not be intended as an alternative system with respect to traditional in situ measurement techniques to estimate the reflectivity of mirrors. Rather, it is a method that can be used continuously to monitor the daily optical efficiency variations and plan reflectometric measurements of mirrors (whenever the calculated index is below a pre- defined threshold) and, if necessary, arrange for the washing of the mirrors. Figure 10. Calculation flowcharts describing the algorithm for the determination of the daily averaged cleanliness index of the solar collector. 4. Results and Discussion 4.1 Simulation of the dish-Stirling performance with clean mirrors The numerical model described in a compact form by Eq. 26 has been calibrated on the clean mirror data collected for the dish-Stirling collector of Palermo and used to simulate the measured performance of the plant, in order to verify both the predictive capacity of the model and to test the hypotheses underlying it. With this aim, the parameters from Tab. 2 were defined previously using both the technical information made available by Ripasso Energy and the results from specific experimental tests carried out on the dish-Stirling unit located in Palermo. The calibration was done by implementing the algorithm described in Fig. 8 using the parameters from Tab. 2 and the measured data from a single day (20 June 2018). On this day, for which the performance has already been represented in Fig. 4, the mirrors of the collector were clean and there was a large variation in DNI values (300 - 800 W/m ). Through the calibration procedure, the following parameters from Eq. 19 were estimated: a =0.477 and a =3.9 kW. Then, using this 1 2 equation, it was possible to calculate the mechanical power output from the Stirling engine as a function of the thermal power input, estimated through Eq. 12, using both the data from June 20 and that from all the other 12 days when the mirrors of the collector were clean. Table 2: Parameters estimated for the Ripasso Energy dish-Stirling unit. Parameters Value Units Effective aperture area A 101 m Receiver aperture area A 0.0314 m Clean reflectivity of mirrors ρ 0.95 ˗ Collector intercept factor γ 0.973 ˗ Effective receiver absorptance α 0.92 ˗ 0.85 ˗ Clean mirrors optical efficiency Receiver convective coefficient h 10 W/(m ·K) Receiver effective emissivity ε 0.88 ˗ 0.924 ˗ Alternator electrical efficiency Figure 11. Stirling engine mechanical output power versus thermal input power. The comparison shows that the Ripasso engine represents an optimised version of the original USAB 4-95. The results of these calculations have been plotted in Fig. 11 alongside the performance data from another Stirling engine, the USAB 4-95, this engine belonging to the same family as the Ripasso engine, the latter being an improved version. The data depicted in Fig 11, which is elaborated from , refers, in particular, to an advanced model with a rated input power of =84.2 kW. This engine is equipped with anti-friction S ,in bearings and operates with a heater-head temperature of T =720 °C and 1800 rpm. It is interesting to note, as a result of this first calibration stage, that: W Q the values of as a function of calculated using the complete dataset from all S S ,in 13 days (when the collector mirrors were clean) are strongly correlated by a linear relationship that was calibrated using data from a single day real data from Stirling engines operating at partial loads show this similar linear thermal-to-mechanical power relation and the two lines, although independently determined, are close to each other. These observations reinforce the validity of the hypotheses underlying the approach described in this work. Taking the results described above, the relationship between the values of and Q , calculated with Eq. 21 using the previously determined values of Q and R are S ,in S ,in plotted in Fig. 12 alongside the analogous efficiency curve plot of the USAB 4-95 engine. The corresponding values of both the Carnot efficiency (calculated by Eq. 17 using the SC , recorded data of T and T ) and the exergetic efficiency (calculated by Eq. 16, Eq. 17 h c S ,ex and Eq. 21) are also shown in the same figure as a function of thermal input power . S ,in Results from Fig. 12 show that the efficiency of the Ripasso Stirling engine ranges between 0.28 and 0.43, corresponding with a thermal input power rising from 17 kW S ,in up to 74 kW. On the other hand, corresponding with the same variations of the thermal input power , the Carnot efficiency of the engine ranges between 0.64 and 0.71, S ,in SC , with an average value of 0.65. These values have been calculated considering that during operation, the highest limiting engine temperature was almost constant and equal to T =993.15 K while the lowest temperature T ranged between 293.15 K and 360.15 K. Finally, as depicted in Fig. 12, the resulting variations in the exergetic efficiency of the Ripasso engine , calculated as a S ,ex function of the partialisation of the thermal power delivered to the engine, are between 0.45 and 0.67. Figure 12. Stirling engine efficiencies versus thermal input power. The latest results, obtained from real-world operating data of a Stirling engine, are consistent with what is shown in the literature for similar engines , . However, it does not seem justified to assume the common hypothesis of considering a constant value for the exergetic efficiency of the engine equal to =0.5 . S ,ex Once the model was calibrated, it was possible to carry out simulations aimed at predicting the net electricity production of the dish-Stirling plant installed in Palermo. To perform these simulations, the procedure illustrated by the calculation flowchart in Fig. 9 was used. A first test to validate the capabilities of the model was carried out by calculating the rated electrical output power of the plant and comparing this value with that declared by the manufacturer. Thus, assuming that the collector mirrors are clean, and considering both a solar irradiation of DNI=960 W/m and an external temperature T =25 °C, the model predicts an air input power Q =80.72 kW, a mechanical output power W =34.66 kW and a gross electric S ,in S power =32.02 kW. Finally, considering parasitic absorptions only related to the cooling system, approximately equal to =1.1 kW, it is possible to calculate a net electric output power equal to =30.92 kW. The latter quantity, calculated conforming to the definition of net electric power proposed by Ripasso Energy, is very close to the rated electric net power reported in Tab. 1. Afterwards, the numerical model was used to simulate the net electric output power data recorded during the 13 days when the collector mirrors of the dish-Stirling unit were clean. In order to perform these simulations, the measured values of T and DNI were air used. Moreover, the total parasitic absorption power was assumed to equal =3.6 kW which is the average of the record measured when the collector mirrors were clean ( cle =1). A comparison of measured and predicted values is depicted in Fig. 13, showing the ability of the proposed model to predict the electric power output generated by the plant. Figure 13. Measured versus predicted net electric output power for the days on which the collector mirrors were clean. Fig. 14 and 15 show the measured and predicted net electric output power of the system for two clear-sky days, during which the ambient temperature was 15 °C and 25 °C, respectively Figure 14. Measured versus predicted net electric power output on 18 January2018. DNI and air temperature variations are also shown. Figure 15. Measured versus predicted net electric power output on 30 May 2018. DNI and air temperature variations are also shown. Finally, using Eqs. 7-10 with the input values of T from the same dataset, it has air been possible to calculate the thermal losses from the collector receiver. These losses present an approximately linear trend with the variations in external temperature, passing Q 174 Q from kW at T =15 °C down to =173 kW at T =30 °C as depicted in Fig. a a r ,out r ,out 16. Therefore, on the basis of these observations, we noticed that assuming a constant value of the thermal losses from the receiver in Eq. 26, equal to the average of the measurements ave =173.5 kW, no substantial errors were introduced into the results of the simulations. r ,out This important observation allows us to verify the assumption behind the linearised models described by Eqs. 27-29. Figure 16. Thermal losses from the receiver versus changes in air temperature recorded during the testing campaign. 4.2 Example of application of the linearised model. To verify the possibilities offered by the simplified linear model described by Eqs. 29-30, we decided to use it to study the performance data of a different plant from that of Palermo: the Stirling Energy System (SES) equipped with the USAB 4-95 Stirling engine max , . For this system the waterfall chart at peak operation condition ( I =1000 W/m ) and the experimental data of net electrical output power as a function of direct normal irradiation, were used. This information, elaborated from the data available in the literature , is illustrated by Fig. 17 and Fig. 18, respectively. Thus, both using the waterfall chart in Fig. 17 and fitting the data in Fig. 18 with a linear regression, we have estimated the minimum parameters necessary for subsequent processing These parameters have been reported in Tab. 3. Thus, assuming a value of the receiver absorptivity equal to α=0.92 it was possible to estimate =0.81 for the system. Finally, from Eqs. 32-34, we obtained a =0.445 and a =2.54 kW. Figure 17. Waterfall chart for SES dish-Stirling (elaborated from ) Figure 18. Net output power versus DNI of SES dish-Stirling (elaborated from ) Table 3: Parameters estimated for SES dish-Stirling unit Parameter Value Units Effective aperture area A 87.7 m Clean reflectivity of mirrors ρ 0.91 ˗ Collector intercept factor γ 0.97 ˗ 0.95 ˗ Alternator electrical efficiency max 1000 W/m Maximum design DNI ave 1.2 kW Average parasitic absorption ave 4.72 kW Average receiver thermal losses r ,out Slope linear relation b 0.030 (kW·m )/W Intercept linear relation b 5.60 kW After the calibration stage it was possible to calculate the relationship between the net electric output power as a function of the variations of the thermal input power to the Stirling engine. These results, along with those elaborated by a previous study on the same system , have been plotted in Fig. 19. With regards to this last result, it is interesting to note that there is a good agreement between the two analyses, despite the fact that the present study is based on a simplified approach using the limited number of information available in literature for the SES system. Figure 19. net electric output power versus thermal input power for SES dish-Stirling calculated by the method elaborated in this study and compared with the results elaborated by the Sandia National Laboratories Finally, we decided to diagram the partial load efficiency curves of both Ripasso and SES Stirling engines as they have been extrapolated from this study and to compare them with the experimental curve of the USAB 4-95 engine. These curves have been plotted in Fig. 20 in terms of normalised efficiency versus normalised thermal input power. max It is possible to define the normalised efficiency as the ratio of over and the max Q Q normalised thermal input as the ratio of over . As expected, the normalised S ,in S efficiency curves of the USAB 4-95, SES and Ripasso engines are substantially coincident since they all belong to the same family. Figure 20. Normalised efficiencies versus normalised input power for the three different Stirling engines (SES, Ripasso and USAB 4-95) in this study. 4.3 Estimation of the daily averaged cleanliness index of the mirrors The results presented so far have mainly concerned the performance of the dish- Stirling plant with clean collector mirrors. The remaining data regarding the periods during which the mirrors were soiled has been analysed using the methodology defined in the calculation flowchart in Fig. 10 with the purpose of defining the time-series plot of the daily average cleanliness index of the mirrors. An example of the application of the proposed method is shown in Fig. 21. In this figure, the net electric power curve predicted by the model assuming clean mirrors is plotted along with the curve corresponding to the condition of soiled mirrors. The application of the algorithm described in Fig. 10 has made it possible to estimate, for this particular day, a value for the daily average cleanliness index day equal to 0.85 . cle Figure 21. Determination of the daily average cleanliness index of mirrors on 11 June 2018 using the algorithm outlined in Fig. 10. The results of all the elaborations are shown in Fig. 22. This figure shows the variation in the calculated cleanliness index function against the rate of rainfall, expressed as (mm/day), and significant Saharan dust transport events (Sirocco), which are highlighted with yellow bars. The pattern of the calculated variations is similar to measurements of reflectivity history for a dish-Stirling system described in the literature . It is possible to observe in Fig. 22 that during the first few months of monitoring, an intense period of rainfall contributed to the maintenance of clean mirrors; however, from point A to point B, there was a notable drop in the cleanliness index to 0.73, due to a prolonged Sirocco event from 28 February until 5 March (see Fig. 23). Immediately afterwards, a period of rainfall between 7 and 12 March corresponding to points B and C on the graph, was echoed by a rise in the cleanliness coefficient to 0.96. From points C to D, there was a gradual reduction in the cleanliness index and then a subsequent increase due to other rain events (point E). In the period between the points E and F on the graph, there is no data since the plant was off. However, it is plausible that during this period a pattern similar to that between A and C occurred as a consequence of a second Sirocco event from 21 to 24 March (see Fig. 24) and a subsequent period of rain. Gradual soiling continued from point F to G since during this time there were no significant weather events. The next Sirocco event, which was much shorter in duration than the first, was recorded in April from 15 to 17 (see Fig. 25). During this event, the cleanliness index dropped to 0.82, point H, and after the subsequent rains rose again to 0.93, point I. Approximately one month passed before manual cleaning and maintenance of the dish was performed, point J. During this time, the rate of soiling was slow and steady, and indeed no other significant weather events were recorded. A marked difference in the cleanliness index can be seen after the mechanical cleaning of 23 May, with this number rising from 0.82 to 1.0 (point K). There was a final cycle of soiling and washing of the mirrors since a soil excavation, that was going on near the plant (between 31 March and 8 June) produced the last notable drop in the cleanliness index and a subsequent sequence of heavy rains (occurring on 14-16 June) produced a complete washing of the collector mirrors (points L to M). The soil excavation period is highlighted with grey bars on the graph. In conclusion, the application of the proposed method has made it possible to obtain a time-series plot of the daily average cleanliness index, which is congruent with the repeated cycles of rainfall and Saharan dust transport events. Moreover, these results confirm what is already known in the literature: wet deposition is the dominant mechanism and is responsible for about 80% of the deposited dust. However, wet deposition is concentrated on rainy days which are quite few per year (such as between points A and B in Fig. 22). Dry deposition, on the other hand, is less dominant but it occurs almost every day  (such as between points I and J in Fig. 22). Figure 22. Time-series plot of the daily average cleanliness index, rainfall and sirocco events. ☼ Published version in Ap p lied Energy : http s://doi.org/10.1016/j.ap energy .2019.114378 available in http s://www.sciencedirect.com/science/article/p ii/S0306261919320653 © <2020>. This manuscrip t version is made available under the under the CC-BY-NC-ND 4.0 license http ://creativecommons.org/licenses/by -nc-nd/4.0/ This p ap er is accep ted for p ublication in Ap p lied Energy (2020). Figure 23. Total dust load (mrg/m ) predicted with SKIRON/Dust system on 2 March 2018 Figure 24. Total dust load (mrg/m ) predicted with SKIRON/Dust system on 21 March ☼ Published version in Ap p lied Energy : http s://doi.org/10.1016/j.ap energy .2019.114378 available in http s://www.sciencedirect.com/science/article/p ii/S0306261919320653 © <2020>. This manuscrip t version is made available under the under the CC -BY-NC-ND 4.0 license http ://creativecommons.org/licenses/by -nc-nd/4.0/ This p ap er is accep ted for p ublication in Ap p lied Energy (2020). Figure 25. Total dust load (mrg/m ) predicted with SKIRON/Dust system on 15 April 2018 With the aim of having a summary picture of all the elaborations, all the experimental data of output power versus DNI were aggregated in three intervals of cleanliness index and plotted in Fig. 26, along with the lines representing the output of the model at fixed indices. The research activities planned for the plant in the future will provide for the automation of the proposed method, which will be integrated with a campaign of periodic measurements of mirror reflectivity in order to consolidate this approach, which will be ultimately aimed at the economic optimisation of the activities of cyclical washing of the collector mirrors. Figure 26. Net electric power output versus DNI (experimental data aggregated in three different cleaning index intervals and model lines at different fixed values of the same index) 5. Conclusion In this paper, we presented a numerical model for the analysis of energy production of dish-Stirling systems that was validated with experimental data from a 32 kW solar plant recently installed at a facility test site on the University of Palermo campus, in the south of Italy. The model was developed by introducing a simplified expression of the Stirling engine part load efficiency into the set of equations describing the energy balance of the system and was calibrated using data from a single day of plant operation of the Palermo dish-Stirling plant. The results of the simulations carried out with the model show excellent agreement with the experimental data of energy production corresponding to the condition of clean collector mirrors. Moreover, the model can be linearised by introducing further simplifying assumptions based on the experimental evidence. The resulting analytical expression, which is a linear relationship between solar irradiation and electric power output from the system, for the first time suggests a theoretical interpretation of what is extensively observed in literature data and which is the basis for the empirical model proposed by Stine for dish-Stirling plants. The main advantage to using the proposed model in this work compared to similar empirical expression is that it is possible to perform optimisation analyses of dish-Stirling components taking into account the main parameters that govern the overall efficiency of the system, including the operating curve of the Stirling engine. The model was also used to define a numerical procedure for assessing the effect of mirror soiling on the optical efficiency of the solar collector that can be integrated with periodic measurements of mirror reflectivity in order to optimise of the activities of cyclical washing of the collector mirrors. The application of the proposed method to the data collected during the testing campaign of the solar collector of Palermo, has made it possible to obtain a time-series plot of the daily average cleanliness index, which is congruent with the repeated cycles of rainfall and Saharan dust transport events. In future works, the numerical approach proposed in this paper will be used to assess the monthly and annual energy performance of dish-Stirling plants built in the central Mediterranean area, taking into account the particular climatic conditions of this geographical location. Such studies are of fundamental importance for the development of guidelines that can direct national policies in favouring a greater commercial penetration of this interesting type of renewable energy production system. Acknowledgements The authors express gratitude to the companies HorizonFirm S.r.l., Elettrocostruzioni S.r.l. and Ripasso Energy for the support provided without which it would not have been possible to install the dish-Stirling concentrator plant at the University of Palermo campus. Nomenclature A total collector aperture area [m ] A effective reflective surface of the collector (projected mirror area) [m ] A receiver equivalent surface area [m ] a slope of the linear relation between W and Q (dimensionless) S S ,in W Q a intercept of the linear the relationship between and [W] S ,in b slope of the linear relationship between and I [kW·m /W] 1 b b intercept of the linear relationship between and I [W] 2 b C geometric concentration ratio of the collector (dimensionless) gross electric output power [W] net electric output power [W] total parasitic absorption power [W] parasitic absorption power of the dry-cooler and circulation pumps [W] pd , parasitic absorption power of the tracking system [W] pt , f sampling frequency [Hz] f cut-off frequency [Hz] h natural convective coefficient at receiver surface [W/(m ·K)] I solar beam radiation [W/m ] rate of thermal convective losses from the receiver [W] con flux of radiant energy emitted from the receiver [W] rad thermal energy power absorbed by the receiver [W] r ,in thermal losses from the receiver [W] r ,out rate of solar energy incident on the mirrors of the collector [W] sun thermal power delivered to the Stirling engine [W] S ,in Q thermal power rejected from the Stirling engine [W] S ,out R ratio of temperature T to temperature T (dimensionless) T o air W rate of mechanical work of the Stirling engine [W] T air temperature [K] air T reference temperature [K] T rejection temperature of the engine working fluid [K] T temperature of the receiver surface [K] T heat input temperature of the engine working fluid [K] T effective sky temperature [K] sky thermal efficiency of the solar collector (dimensionless) average electric generator efficiency (dimensionless) cleanliness index of the collector mirror (dimensionless) cle factor accounting for the losses of reflective mirror area (dimensionless) effective optical efficiency (dimensionless) efficiency of the Stirling engine (dimensionless) Carnot cycle efficiency of the Stirling engine (dimensionless) SC , exergetic efficiency of the Stirling engine (dimensionless) S ,ex N number of samples in the data record (dimensionless) Superscripts ave averaged value day mean daily value max maximum design value Greek letters effective absorptivity of the receiver (dimensionless) effective emissivity factor of the receiver (dimensionless) clean mirror reflectance (dimensionless) intercept factor of the solar collector (dimensionless) -8 2 4 Stefan–Boltzmann constant, 5.67×10 [W/(m ·K )] References  Bilgili M, Ozbek A, Sahin B, Kahraman A. An overview of renewable electric power capacity and progress in new technologies in the world. Renew Sustain Energy Rev 2015;49:323–34. doi:10.1016/j.rser.2015.04.148.  Joint Research Centre. Energy Technology Reference Indicator projections for 2010 - 2050. 2014. doi:10.2790/057687.  Zhang HL, Baeyens J, Degrève J, Cacères G. Concentrated solar power plants: Review and design methodology. Renew Sustain Energy Rev 2013;22:466–81. doi:10.1016/j.rser.2013.01.032.  Zhu S, Yu G, Ma Y, Cheng Y, Wang Y, Yu S, et al. A free-piston Stirling generator integrated with a parabolic trough collector for thermal-to-electric conversion of solar energy. Appl Energy 2019;242:1248–58. doi:10.1016/j.apenergy.2019.03.169.  Buscemi A, Panno D, Ciulla G, Beccali M, Lo Brano V. Concrete thermal energy storage for linear Fresnel collectors: Exploiting the South Mediterranean’s solar potential for agri-food processes. Energy Convers Manag 2018;166. doi:10.1016/j.enconman.2018.04.075.  Ogunmodimu O, Okoroigwe EC. Concentrating solar power technologies for solar thermal grid electricity in Nigeria: A review. Renew Sustain Energy Rev 2018;90:104–19.  Sharma A, Shukla SK, Rai AK. Finite time thermodynamic analysis and optimization of solar-dish Stirling heat engine with regenerative losses. Therm Sci 2011;15:995–1009. doi:10.2298/TSCI110418101S.  Thombare DG, Verma SK. Technological development in the Stirling cycle engines. Renew Sustain Energy Rev 2008;12:1–38. doi:10.1016/j.rser.2006.07.001.  Mancini T, Heller P, Butler B, Osborn B, Schiel W, Goldberg V, et al. Dish-Stirling Systems: An Overview of Development and Status. J Sol Energy Eng 2003;125:135. doi:10.1115/1.1562634.  Coventry J, Andraka C. Dish systems for CSP. Sol Energy 2017;152:140–70. doi:10.1016/J.SOLENER.2017.02.056.  Ferreira AC, Nunes ML, Teixeira JCF, Martins LASB, Teixeira SFCF. Thermodynamic and economic optimization of a solar-powered Stirling engine for micro-cogeneration purposes. Energy 2016;111:1–17. doi:10.1016/j.energy.2016.05.091.  Reader GT, Hooper C. Stirling Engines. 1983. London E FN Spon n.d.  Al-Dafaie AMA, Dahdolan ME, Al-Nimr MA. Utilizing the heat rejected from a solar dish Stirling engine in potable water production. Sol Energy 2016;136:317–26. doi:10.1016/j.solener.2016.07.007.  Hafez AZ, Soliman A, El-Metwally KA, Ismail IM. Design analysis factors and specifications of solar dish technologies for different systems and applications. Renew Sustain Energy Rev 2017. doi:10.1016/j.rser.2016.09.077.  Singh UR, Kumar A. Review on solar Stirling engine: Development and performance. Therm Sci Eng Prog 2018. doi:10.1016/j.tsep.2018.08.016.  Bǎdescu V. Note concerning the maximal efficiency and the optimal operating temperature of solar converters with or without concentration. Renew Energy 1991;1:131–5. doi:10.1016/0960-1481(91)90114-5.  Chen L, Sun F, Wu C. Optimum collector temperature for solar heat engines. Int J Ambient Energy 1996;17:73–8. doi:10.1080/01430750.1996.9675221.  Kongtragool B, Wongwises S. Optimum absorber temperature of a once-reflecting full conical concentrator of a low temperature differential Stirling engine. Renew Energy 2005;30:1671–87. doi:10.1016/j.renene.2005.01.003.  Ahmadi MH. Investigation of Solar Collector Design Parameters Effect onto Solar Stirling Engine Efficiency. J Appl Mech Eng 2012;01:10–3. doi:10.4172/2168- 9873.1000102.  Reddy VS, Kaushik SC, Tyagi SK. Exergetic analysis and performance evaluation of parabolic trough concentrating solar thermal power plant (PTCSTPP). Energy 2012. doi:10.1016/j.energy.2012.01.023.  Duffie JA, Beckman WA. Solar engineering of thermal processes. John Wiley & Sons; 2013.  Ahmadi MH, Sayyaadi H, Mohammadi AH, Barranco-Jimenez MA. Thermo- economic multi-objective optimization of solar dish-Stirling engine by implementing evolutionary algorithm. Energy Convers Manag 2013;73:370–80. doi:10.1016/j.enconman.2013.05.031.  Carrillo Caballero GE, Mendoza LS, Martinez AM, Silva EE, Melian VR, Venturini OJ, et al. Optimization of a Dish Stirling system working with DIR-type receiver using multi-objective techniques. Appl Energy 2017. doi:10.1016/j.apenergy.2017.07.053.  Beltrán-Chacon R, Leal-Chavez D, Sauceda D, Pellegrini-Cervantes M, Borunda M. Design and analysis of a dead volume control for a solar Stirling engine with induction generator. Energy 2015. doi:10.1016/j.energy.2015.09.046.  Hafez AZ, Soliman A, El-Metwally KA, Ismail IM. Solar parabolic dish Stirling engine system design, simulation, and thermal analysis. Energy Convers Manag 2016. doi:10.1016/j.enconman.2016.07.067.  Mendoza Castellanos LS, Carrillo Caballero GE, Melian Cobas VR, Silva Lora EE, Martinez Reyes AM. Mathematical modeling of the geometrical sizing and thermal performance of a Dish/Stirling system for power generation. Renew Energy 2017. doi:10.1016/j.renene.2017.01.020.  Power CS. Technology Roadmap Concentrating Solar Power. Current 2010. doi:10.1787/9789264088139-en.  Guo S, Liu Q, Sun J, Jin H. A review on the utilization of hybrid renewable energy. Renew Sustain Energy Rev 2018. doi:10.1016/j.rser.2018.04.105.  Bravo Y, Carvalho M, Serra LM, Monné C, Alonso S, Moreno F, et al. Environmental evaluation of dish-Stirling technology for power generation. Sol Energy 2012;86:2811–25. doi:10.1016/j.solener.2012.06.019.  Schiel W, Schweiber A, Stine WB. Evaluation of the 9-kw e dish/stirling system of schlaich bergermann und partner using the proposed iea dish/stirling performance analysis guidelines. Intersoc. Energy Convers. Eng. Conf. 1994, 1994.  Stine WB. Experimentally Validated Long-Term Energy Production Prediction Model for Solar Dish/Stirling Electric Generating Systems. In: D.Y. Goswami, L.D. Kannberg, T.R. Mancini, S. Somasundaram eds., editor. Proc. Intersoc. Energy Convers. Eng. Conf. vol. 2, New York, NY, USA: American Society of Mechanical Engineers; 1995, p. 491–5. doi:AC04-94AL85000.  Igo J, Andraka CE. Solar dish field system model for spacing optimization. Proc. Energy Sustain. Conf. 2007, 2007. doi:10.1115/ES2007-36154.  Nepveu F, Ferriere A, Bataille F. Thermal model of a dish/Stirling systems. Sol Energy 2009. doi:10.1016/j.solener.2008.07.008.  Mendoza Castellanos LS, Galindo Noguera AL, Carrillo Caballero GE, De Souza AL, Melian Cobas VR, Silva Lora EE, et al. Experimental analysis and numerical validation of the solar Dish/Stirling system connected to the electric grid. Rene w Energy 2019. doi:10.1016/j.renene.2018.11.095.  Andraka CE. Cost/performance tradeoffs for reflectors used in solar concentrating dish systems. 2008 Proc. 2nd Int. Conf. Energy Sustain. ES 2008, 2009.  Lopez CW, Stone KW. Performance of the Southern California Edison Company Stirling Dish. NASA STI/Recon Tech Rep N 1993;94.  Panno D, Buscemi A, Beccali M, Chiaruzzi C, Cipriani G, Ciulla G, et al. A solar assisted seasonal borehole thermal energy system for a non-residential building in the Mediterranean area. Sol Energy 2018. doi:10.1016/j.solener.2018.06.014.  Sciortino L, Agnello S, Barbera M, Bonsignore G, Buscemi A, Candia R, et al. Direct sunlight facility for testing and research in HCPV. AIP Conf. Proc., vol. 1616, 2014, p. 158–61.  Kallos G, Astitha M, Katsafados P, Spyrou C. Long-Range Transport of Anthropogenically and Naturally Produced Particulate Matter in the Mediterranean and North Atlantic: Current State of Knowledge. J Appl Meteorol Climatol 2007;46:1230–51. doi:10.1175/JAM2530.1.  Kallos G, Spyrou C, Astitha M, Mitsakou C, Solomos S, Kushta J, et al. Ten-year operational dust forecasting – Recent model development and future plans. IOP Conf. Ser. Earth Environ. Sci., vol. 7, 2009, p. 12012. doi:10.1088/1755- 1307/7/1/012012.  Kallos G, Solomos S, Kushta J, Mitsakou C, Spyrou C, Bartsotas N, et al. Natural and anthropogenic aerosols in the Eastern Mediterranean and Middle East: Possible impacts. Sci Total Environ 2014;488–489:389–97. doi:10.1016/j.scitotenv.2014.02.035.  Spyrou C, Kallos G, Mitsakou C, Athanasiadis P, Kalogeri C. The Effects of Naturally Produced Dust Particles on Radiative Transfer. Adv. Meteorol. Climatol. Atmos. Phys., Springer; 2013, p. 317–23.  Hachicha AA, Yousef BAA, Said Z, Rodríguez I. A review study on the modeling of high-temperature solar thermal collector systems. Renew Sustain Energy Rev 2019. doi:10.1016/j.rser.2019.05.056.  Selcuk MK, Fujita T. A nomographic methodology for use in performance trade-off studies of parabolic dish solar power modules 1984.  Massi Pavan A, Mellit A, De Pieri D. The effect of soiling on energy production for large-scale photovoltaic plants. Sol Energy 2011. doi:10.1016/j.solener.2011.03.006.  Gostein M, Stueve B, Chan M. Accurately measuring PV soiling losses with soiling station employing PV module power measurements. 2017 IEEE 44th Photovolt. Spec. Conf. PVSC 2017, 2017. doi:10.1109/PVSC.2017.8366169.  Smith SW. The scientist and engineer’s guide to digital signal processing. C alifornia Technical Pub; 1997.  Gil R, Monné C, Bernal N, Muñoz M, Moreno F. Thermal model of a dish stirling cavity-receiver. Energies 2015;8:1042–57. doi:10.3390/en8021042.  Eldighidy SM. Optimum outlet temperature of solar collector for maximum work output for an Otto air-standard cycle with ideal regeneration. Sol Energy 1993;51:175–82. doi:10.1016/0038-092X(93)90094-5.  Hogan RE. AEETES-A solar reflux receiver thermal performance numerical model. Sol Energy 1994. doi:10.1016/0038-092X(94)90066-3.  Yaqi L, Yaling H, Weiwei W. Optimization of solar-powered Stirling heat engine with finite-time thermodynamics. Renew Energy 2011;36:421–7. doi:10.1016/j.renene.2010.06.037.  García D, González MA, Prieto JI, Herrero S, López S, Mesonero I, et al. Characterization of the power and efficiency of Stirling engine subsystems. Appl Energy 2014;121:51–63. doi:10.1016/j.apenergy.2014.01.067.  Kallos G, Papadopoulos A, Katsafados P, Nickovic S. Transatlantic Saharan dust transport: Model simulation and results. J Geophys Res 2006;111:D09204. doi:10.1029/2005JD006207.  Kallos G, Katsafados P, Spyrou C, Papadopoulos A. Desert dust deposition over the Mediterranean Sea estimated with the SKIRON/Eta - System validation. 4th EuroGOOS Conf., 2005.  Querol X, Pey J, Pandolfi M, Alastuey A, Cusack M, Pérez N, et al. African dust contributions to mean ambient PM10 mass-levels across the Mediterranean Basin. Atmos Environ 2009;43:4266–77.  Kallos G, Spyrou C, Astitha M, Mitsakou C, Solomos S, Kushta J, et al. Ten-year operational dust forecasting--Recent model development and future plans. IOP Conf. Ser. Earth Environ. Sci., vol. 7, 2009, p. 12012.  Mitsakou C, Kallos G, Papantoniou N, Spyrou C, Solomos S, Astitha M, et al. Saharan dust levels in Greece and received inhalation doses. Atmos Chem Phys 2008;8:7181–92. doi:10.5194/acp-8-7181-2008.  Astitha M, Kallos G. Gas-phase and aerosol chemistry interactions in South Europe and the Mediterranean region. Environ Fluid Mech 2009;9:3–22.  National & Kapodistrian University of Athens School of Physics. SKIRON/Dust system n.d. http://forecast.uoa.gr/dustindx.php?domain=med (accessed October 10, 2019).  Solomos S, Kallos G, Kushta J, Astitha M, Tremback C, Nenes A, et al. An integrated modeling study on the effects of mineral dust and sea salt particles on clouds and precipitation. Atmos Chem Phys 2011;11:873–92. doi:10.5194/acp-11- 873-2011.
Physics – arXiv (Cornell University)
Published: Jan 3, 2020
Access the full text.
Sign up today, get DeepDyve free for 14 days.