Next Article in Journal
Experimental Study on Performance of Liquid–Gas Jet Pump with Square Nozzle
Previous Article in Journal
A Physics-Informed Neural Network Approach for Surrogating a Numerical Simulation of Fractured Horizontal Well Production Prediction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of Circumsolar Irradiance and Its Impact on CSP Systems under Clear Skies

by
Edgar F.M. Abreu
1,*,
Paulo Canhoto
1,2 and
Maria João Costa
1,3,4
1
Instituto de Ciências da Terra, Universidade de Évora, Rua Romão Ramalho 59, 7000-671 Évora, Portugal
2
Departamento de Engenharia Mecatrónica, Escola de Ciências e Tecnologia, Universidade de Évora, Rua Romão Ramalho 59, 7000-671 Évora, Portugal
3
Departamento de Física, Escola de Ciências e Tecnologia, Universidade de Évora, Rua Romão Ramalho 59, 7000-671 Évora, Portugal
4
Earth Remote Sensing Laboratory—EaRSLab, Universidade de Évora, Rua Romão Ramalho 59, 7000-671 Évora, Portugal
*
Author to whom correspondence should be addressed.
Energies 2023, 16(24), 7950; https://doi.org/10.3390/en16247950
Submission received: 14 October 2023 / Revised: 27 November 2023 / Accepted: 2 December 2023 / Published: 7 December 2023
(This article belongs to the Section A2: Solar Energy and Photovoltaic Systems)

Abstract

:
In this work, a model to estimate circumsolar normal irradiance (CSNI) for several half-opening angles under clear skies was developed. This approach used a look-up table to determine the model parameters and estimate CSNI for half-opening angles between 0.5° and 5°. To develop and validate the proposed model, data from five locations worldwide were used. It was found that the proposed model performs better at the locations under study than the models available in the literature, with relative mean bias error ranging from −13.94% to 0.70%. The impact of CSNI for these different half-opening angles on concentrating solar power (CSP) systems was also studied. It was found that neglecting CSNI could lead to up to a 7% difference between the direct normal irradiance (DNI) measured by a field pyrheliometer and the DNI that is captured by CSP systems. Additionally, a case study for parabolic trough concentrators was performed as a way to estimate the impact of higher circumsolar ratios (CSR) on the decrease of the intercept factor for these systems. It was also concluded that if parabolic trough designers aim to reduce the impact of CSNI variation on the intercept factor, then parabolic troughs with higher rim angles are preferred.

1. Introduction

Concentrating solar power (CSP) systems highly depend on the intensity of direct normal irradiance (DNI) in their aperture angle [1]. Therefore, it is crucial that designers have access to accurate DNI data to avoid compromising the proper operation and economic viability of such systems. The most reliable way of obtaining accurate DNI data is to use a ground-based station equipped with a pyrheliometer in a sun-tracker system. This instrument needs to track the apparent movement of the sun along its path in the sky to measure the direct beam. However, as this tracking cannot be done with perfect accuracy, these instruments are designed in such a way that they receive radiation not only from the solar disk but also from the surrounding area, known as the circumsolar region [2]; that is, the aperture half-angle of the instrument is larger than the apparent sun disk radius in order to accommodate the tracking misalignments and errors. On the other hand, this radiosity in the circumsolar region around the sun disk, known as circumsolar irradiance, is not negligible when designing CSP systems, because these systems usually have opening half-angles lower than those of the pyrheliometers, depending on the concentration factor.
The World Meteorological Organization (WMO) recommends that pyrheliometers have an opening half-angle of 2.5° to enable a fair comparison between the DNI measured at different locations around the globe [3]. However, the upper limit of the opening half-angles of common CSP technologies such as parabolic trough, linear Fresnel, dish-Stirling and solar tower are 0.8°, 1.0°, 1.6°, and 1.8°, respectively [2,4], that is, lower than the standard aperture half-angles of pyrheliometers. In other words, CSP systems will not receive the same levels of irradiance as measured by collocated standard pyrheliometers. This difference will increase when the circumsolar radiosity also increases as a result of solar radiation scattering in the atmosphere due to a higher concentration of aerosols. The quantification of this difference between the measured DNI and the DNI that is reflected by the mirrors of a given CSP system is one of the main motivations of this study.
The air molecules, aerosols, and cirrus clouds scatter the sun rays from the sun disk region to the circumsolar region [5]. When this irradiance is measured in a surface normal to the sun line, it is known as circumsolar normal irradiance (CSNI). In that sense, while experimental values of DNI obtained from pyrheliometer measurements are related only to the total scattering effect and transmittance of the atmosphere, the CSNI is also related to the angular distance from the center of the sun disk, i.e., the scattering angle [6,7]. On the other hand, measurement of CSNI is not straightforward because of the sharp decrease in intensity between the edge of the sun disk and the limit of the field of view of the instrument [2]. However, some attempts to measure CSNI have been made and a review on these can be found in Abreu et al. [6].
Regarding CSNI modelling, it is usually carried out using radiative transfer models (RTM) such as libRadtran [8]. However, because RTM models often need high-quality atmospheric data such as aerosol optical depth and precipitable water vapor, some authors have proposed simpler models to predict CSNI. In these models, CSNI is usually estimated as a function of more commonly available solar irradiance data and presented as a circumsolar ratio (CSR), i.e., the ratio between CSNI and DNI [5,6,9,10]. Such models are defined for a fixed aperture angle except the model from Eissa et al. [5] in which, in addition, the authors presented a way of varying the aperture angle. This takes a similar form of the fixed aperture angle model, but the model coefficients are determined using a sixth-degree polynomial of the aperture angle and allows the CSNI to be estimated for any defined aperture angle in the range of 0.4°–5.0°. Although the model of Abreu et al. [6] was initially developed for a fixed aperture angle, the way it was constructed also allows its generalization for any aperture angle, with the advantage that it only needs widely available solar irradiance data to predict CSNI; this being one of the main objectives of this study.
Furthermore, information on CSNI is crucial for the design and operation of CSP systems. In CSP plants, the solar irradiance is focused onto an absorber using mirrors [11]. However, due to the difference between the opening half-angles of the pyrheliometers and the CSP systems mentioned above, the use of less accurate DNI data (without correctly accounting for the differences of the opening half-angle, that is, without considering CSNI) leads to incorrect assumptions about the solar irradiance that is focused by the CSP mirrors onto the receiver. This will increase the uncertainty of the energy generation predictions of CSP systems, which ultimately increases the difficulty of its operation.
The impact of CSNI on CSP systems can be assessed using a variety of tools, namely, ray tracing tools, analytical optical performance models, and models that use look-up tables or parametrizations of the solar position relative to the CSP mirrors [2]. Ray tracing models describe the solar irradiance as a multitude of solar rays that originate from the sun, reach the concentrator, and finally the receiver. Examples of ray tracing models can be STRAL [12] and MIRVAL [13], to name a few. Analytical models use logical equations that can describe the ray’s path through the optical system. Examples of analytical models are the Bendt-Rabl model [14] and the HFLCAL method [15], to name a few. Lastly, models based on look-up tables use parametrizations or look-up tables to assess the optical performance of a CSP collector according to the solar position. SAM [16] and Greenius [17,18] are examples of models that use a look-up table. The approach followed in this study was to use an analytical model to determine how the intercept factor (the ratio of the irradiance reaching the receiver over the incident irradiance) is affected according to CSR.
In this work, a simple and fast model to derive CSNI is proposed for clear-sky conditions at any desired aperture angle. In other words, the proposed model is able to determine the CSNI that reaches the receiver of any CSP system, such as parabolic troughs, tower systems, or dish systems, to name a few. The developed model also aims to identify the difference between using ground-based DNI measurements (for the aperture half-angle of the field pyrheliometer) and DNI modelled data according to the aperture half-angles of typical CSP technologies (with lower aperture angles) at several locations around the globe. This allows quantification of the errors associated with using directly measured DNI data on the simulation of CSP systems. To the authors’ best knowledge, the model proposed in this work is the only fast and simple model capable of estimating CSNI at several locations in more than one climate zone, which constitutes another novelty here. Furthermore, the impact of CSNI on CSP systems, specifically in parabolic trough systems, was also studied using the analytical optical performance model from Bendt and Rabl [14].
This paper is organized as follows: Section 2 describes the data used; in Section 3, the model development and performance assessment are presented; in Section 4, prediction of CSNI and its impact on CSP systems, namely, in a parabolic trough, is assessed; and in Section 5, conclusions are drawn.

2. Data

2.1. Experimental Solar Radiation Data

In this work, CSNI data corresponding to several half-opening angles necessary to develop the proposed model were obtained using the version 1.7 of the libRadtran RTM [8] simulations based on AERONET [19] input data for five locations worldwide. AERONET is a ground-based aerosol network that has provided aerosol optical data and microphysical and radiative properties for aerosol research and characterization for more than 25 years. The values retrieved from the AERONET database were the following: aerosol optical depth, aerosol single scattering albedo, aerosol phase function, surface albedo, and precipitable water vapor. The AERONET data were retrieved at the wavelengths of 440 nm, 675 nm, 870 nm, and 1020 nm. However, the simulations were carried out for the wavelength interval between 200 nm and 5000 nm. This wavelength interval covers not only the spectral response of a conventional pyrheliometer (300–4000 nm) but also the spectral response of windowless absolute cavity radiometers (the 5000 nm upper limit), against which field pyrheliometers are commonly calibrated. To extrapolate the AERONET aerosol data to 200 nm and 5000 nm, Ångström’s law was used. For further details on this process, the reader can consult a previous work by Abreu et al. [6]. Ground-based DNI values from the BSRN database were also used in the validation process of the CSNI data for the same locations. BSRN is a network of ground-based radiometric stations with the aim of detecting changes in the solar radiation field at the Earth’s surface [20].
The five locations analyzed in this work are scattered around the globe and have collocated AERONET and BSRN stations, thus allowing measurements to be assumed from both stations at close timestamps. The time difference between AERONET and the radiometric stations data is not higher than 1 min and therefore it is assumed that, under clear skies, they are representative of the same atmospheric conditions. More information about the locations is presented in Table 1.
Regarding data quality control, AERONET data Level 1.5 Version 3 was used here for the following two main reasons: to use single scattering albedo measurements instead of mean values, and because the Version 3 processing algorithm marks a significant improvement in the quality control of the sun photometer AOD measurements, eliminating the need for manual quality control and cloud screening by an analyst [21,22]. Concerning solar radiation data, the quality filters from BSRN [20] were used to ensure that extremely rare or physically impossible values were discarded.
The quality control procedure used here ensures the following: (i) that the simulations to determine DNI and CSNI are as accurate as possible; (ii) that the model developed here was created and validated using the most accurate solar irradiance ground-based data available. Moreover, the data used in this work were screened and classified as clear sky by the AERONET algorithm alongside solar radiation data from collocated radiometric stations, in order to ensure that both AERONET and BSRN/Evora data are representative of the same atmospheric conditions.

2.2. Modelled DNI and CSNI Data

In this work, modelled DNI and CSNI data from the libRadtran RTM are used to develop the proposed CSNI model. The generation of these kind of data is complex because of (i) the availability and processing of the required inputs, (ii) the computation time required to simulate extensive datasets such as the one used here. Therefore, the DNI and CSNI data were generated and validated in a previous work from the same authors [6]. However, it is worth mentioning that the data processing in this work is not the same as in the previously mentioned study [6].
In the previous study mentioned above [6], the libRadtran radiative transfer model was used alongside AERONET to generate DNI and the sky radiance. The latter was then used to generate CSNI data through the integration of sky radiance from the RTM model output in the circumsolar region. The quality of the modelled DNI and CSNI data was then assessed using ground-based DNI data. This validation process was indirect, i.e., the DNI modelled data was compared against the ground-based DNI data, with and without its CSNI counterpart. Since the comparison between modelled and ground-based data was better when the modelled CSNI counterpart was used, it was then concluded that the modelled CSNI data exhibited an acceptable accuracy. Further information of the modelled data, the procedure to generate it, and the validation procedure can be found in [6].

3. Model Development and Assessment

3.1. Model Development

Each pyrheliometer model has its own penumbra function because it has its own geometrical characteristics. This function accounts for the pyrheliometer’s response as a function of the scattering angle due to the effect of the opening window edge, which attenuates the intensity of radiation in a transitions range between illuminated and non-illuminated areas [6]. However, assuming a simpler model, the penumbra function can be assumed as a rectangular function [2], i.e., there is no transition range, and an ideal DNI for the opening half-angle α can be defined as follows:
B n i d e a l α 2 π 0 α L ξ sin ( ξ ) d ξ ,
where L is the sky radiance and ξ is the scattering angle (the angular distance from the center of the sun). In this way, the direct normal irradiance from the sun disk ( B n s u n ), can be defined replacing α by the half-angle of the solar disk ( δ s ) in Equation (1). In a similar way, the ideal CSNI, C S n i d e a l , can be defined as follows:
C S n i d e a l δ S , α 2 π δ s α L ξ sin ξ d ξ ,
where the constraint α δ S is required. In ideal conditions, following a fundamental closure relationship, we can define the ideal DNI as follows:
B n i d e a l α = C S n i d e a l δ s , α + B n s u n ,
The circumsolar ratio (CSR) for this ideal case can also be defined as follows:
C S R α = C S n i d e a l δ s , α B n i d e a l α ,
The model proposed in this work is based on these fundamental relationships regarding the simpler penumbra model of the instruments and, therefore, CSNI and CSR can be determined using only α. To that end, starting from the results of the libRadtran simulations [6], CSNI and CSR were calculated for half-opening angles ranging from 0.5° to 5° with steps of 0.1°, resulting in 46 data points for each record timestamp in each location. Then, the adjustable parameters a, b, and c from the CSR model for a given half-opening angle developed by Abreu et al. [6] (described below by Equation (5) and collocated text) were fitted to the data points for each value of α, for all five locations. Finally, the model parameters a, b, and c (46 data points for each parameter) were modelled by fitting a polynomial relation with α according to Equation (6), where f(α) represents the required model parameter. The referred CSR model [6] is given by the following:
C S R 1 K d = a b + ( K t K b / K d ) c ,
where K t is the sky clearness index, i.e., the ratio between global horizontal irradiance at the Earth’s surface and its counterpart at the top of the atmosphere [23], K b is the beam clearness index, i.e., the ratio between DNI at the Earth’s surface and DNI at the top of atmosphere [24], and K d is the diffuse fraction, i.e., the ratio between diffuse horizontal and global horizontal irradiance [25]. The general polynomial form to model these parameters as a function of the half-opening angle is in the following form:
f α = C 0 + C 1 α + C 2 α 2 ,
Regarding the polynomial fitting to determine the new model parameters, it was verified that the parameter a can be obtained accurately through a first-degree polynomial, while parameters b and c are better adjusted using a second-degree polynomial. The polynomial coefficients to obtain the model parameters are presented in Table 2 according to each location as well as their respective coefficients of determination. The locations are grouped according to climatic zone. The model parameters were obtained using a nonlinear least squares method through the Matlab® fit function. It is worth noting that training and validation datasets comprised of approximately two thirds and one third of the original datasets shown in Table 1, respectively, were created, in order to avoid overfitting.
The polynomial coefficients derived here differ from station to station and from climate zone to climate zone. This can be an indication of the differences in the mean aerosol values and other relevant local meteorological or surface conditions at various stations and climate zones. For example, the atmospheric conditions at EVR and SMS should be very different even though these two stations belong to the same climate zone. Whilst the EVR station is located at a semi-rural small city, the SMS station is located in a rural area, causing these two sites to have different aerosol profiles. Whilst at SMS the predominant aerosol should be the rural aerosol, at EVR a combination of rural and urban aerosols is shown. This may be the reason why a global model (using all of the data) would not produce acceptable results and, therefore, is not shown here.

3.2. Performance Assessment

To assess the performance of the proposed procedure based on the model of Abreu et al. [6], including the determination of model parameters through Equation (6), the CSNI values obtained here and the values from the model of Eissa et al. [5] were compared against the libRadtran CSNI values, for the following values of α: 0.8°, 1.0°, 1.6°, and 1.8°. These α values were selected because they represent the upper bound acceptance half-angles of the following CSP technologies: parabolic through, linear Fresnel, dish-Stirling, and solar tower, respectively [2,4]. The performance assessment was carried out using the following statistical indicators: relative mean bias error (rMBE), relative root mean square error (rRMSE), fractional bias (FB), fractional gross error (FGE), and correlation coefficient (R), defined as follows, respectively.
r M B E = 1 n i = 1 n p i o i o ¯ ,
r R M S E = 1 n i = 1 n ( p i o i ) 2 o ¯
F B = 2 n i = 1 n p i o i p i + o i
F G E = 2 n i = 1 n p i o i p i + o i
R = i = 1 n ( p i p ¯ ) ( o i o ¯ ) i = 1 n ( p i p ¯ ) 2 i = 1 n ( o i o ¯ ) 2
where p i are the model prediction, o i are the corresponding libRadtran simulations, p ¯ is the average of the model predictions, o ¯ is the average of the libRadtran simulations, and n is the total number of data points.
Since there are no CSNI measured data at the required half-opening angles and locations, the proposed CSNI model was compared against the models from Eissa et al. [5]. In the work of Eissa et al. [5], the authors developed three models: one model for Sollar Village (SV in their work, SOV here), one model for Tamanrasset (TAM), and a third model that resulted in a combination of the data from these two stations. It is worth mentioning that these models were derived for a desert environment (i.e., arid climate zone) and its performance is being assessed here at different climate zones. However, to the best knowledge of the authors, there are no other CSNI models available in the literature derived for the remaining climate zones. In this way, the statistical analysis results of the models from Eissa et al. [5] and the model proposed here are presented in Table 3 as an example for the half-opening angles of 0.8°. Values in bold identify the model with the best results according to each statistical indicator. The statistical analysis for the remaining half-opening angles of 1.0°, 1.6°, and 1.8° can be found in Appendix A.
This statistical analysis shows that the model proposed in this work generates more accurate results than the other models available in the literature for the datasets used in this study, regardless of the value of the half-opening angle α, for all climate zones (see also tables in Appendix A).
Looking more closely at the results, the proposed CSNI model performs better at DAR (tropical climate zone) the higher that α is, according to all statistical indicators except rMBE. The same overall trend is shown for the stations in the temperate climate zone (EVR and SMS stations). Regarding the arid climate zone, all statistical indicators of the present model show better results with the increase of α at GOB. However, the same is not readily shown at TAM, where the only statistical parameters that increase performance with the increase of α are R and FGE.
The aerosol characteristics and meteorological conditions of the different climate zones have a direct impact on the results presented here. The mean statistical indicators for each climate zone show that the tropical climate zone is where the proposed model performs best according to rMBE, R, FB, and FGE. Regarding the R statistical indicator, the climate zone in which the proposed model provides the best results is the arid climate zone.

4. Prediction of CSNI and of Its Impact on the Energy Capture of CSP Systems

4.1. Impact on the Energy That Reaches the Aperture of CSP Systems

In the work of Abreu et al. [6], modelled DNI values with and without CSNI were validated against ground-based DNI measurements for an aperture half-angle equal to the typical field pyrheliometer. The modelled DNI and CSNI values were then obtained through the integration of sky radiance from libRadtran simulations on the aperture solid angle of the pyrheliometer, with an overall uncertainty of 4.44%. This overall accuracy value is for the entire data set used in this work that results from the AERONET clear-sky screening algorithm. To establish the dependence of the simulation accuracy with the apparent position of the sun in the sky and the impact on the proposed model performance, a detailed analysis procedure is further needed in the future which is beyond the scope of the present work. In the present work, the validated DNI data are compared against the modelled DNI obtained in the same way but for the typical half-opening values of the CSP technologies: 0.8°, 1.0°, 1.6°, and 1.8°, for parabolic through, linear Fresnel, dish-Stirling, and solar tower, respectively [2,4]. This comparison was performed using the mean average difference (in percentage) between the DNI corresponding to the half-opening angle of the pyrheliometer installed at each location (DNI(αpyr)) and the modelled DNI corresponding to each of the half-opening angles mentioned above (DNI(α)), and is a measure of the error or bias between the energy assessment using the pyrheliometer values and the real value of irradiance that is captured by the CSP systems, as shown in Figure 1. In the stations DAR, SMS, and TAM, DNI is measured using an Eppley NIP (α = 2.9°) whilst at EVR and GOB, DNI is measured using a Kipp&Zonnen CHP1 (α = 2.5°).
The average difference between D N I ( α p y r ) and D N I ( α ) is higher for lower half-opening angles as expected, and for locations with higher CSR. The locations with higher average difference are TAM and DAR, whilst locations such as GOB and SMS show low average differences between the DNI of the pyrheliometer and the DNI for lower half-opening angles. This is related to the CSR magnitude across different locations, and in turn, with the composition of the atmosphere, in particular with the type and concentration of aerosols. It is worth mentioning that not only the magnitude of the average difference across locations depends on the composition of the atmosphere, but also the slope of the lines shown in Figure 1. Locations with higher aerosol concentration have higher slopes (e.g., TAM) as can be seen in Figure 2, where boxplots for aerosol optical depth (AOD) and single scattering albedo (SSA) at 675 nm are shown. Figure 1 and Figure 2 highlight the need for the consideration of both the half-opening aperture angle and the CSNI when designing and simulating CSP systems, especially in locations with higher aerosol concentration. Differences in the incoming DNI of around 7% and 2% for locations such as TAM and DAR, respectively, can lead to erroneous energy generation estimates, which in turn can jeopardize the economic viability of the CSP projects.

4.2. Impact on the Intercept Factor of CSP Systems: The Case of Parabolic trough Concentrators

The model from Bendt et al. [14] was used here to quantify in a simple way the impact of CSR in the optical efficiency of a parabolic trough CSP system with cylindrical receiver.
In that work [14], an analytical approach was proposed to perform the optical analysis of a solar concentrator instead of using ray-tracing software. The former approach is simpler and faster than the later, hence it was used here. By assuming a Gaussian distribution for the sun shape, Bendt et al. [14] were able to determine the intercept factor (the ratio of the irradiance reaching the receiver over the incident irradiance) as a function of the group σ C , where σ stands for total optical error and C for concentration ratio.
To achieve this, the model by Bendt et al. [14] determines an effective source that accounts for the shape of the sun, and for all optical errors from the parabolic trough system. In this way, it is possible to determine the intercept factor.
The total optical error ( σ ) is defined as the root mean square of angular spread caused by all optical errors of the concentrator ( σ o p t ) and the angular width of the sun shape in line focus geometry ( σ s u n ), through the following equation:
σ = σ o p t 2 + σ s u n 2 ,
To study the impact of CSNI (or CSR) variation on CSP systems, namely, parabolic trough systems, a relationship between σ s u n and CSR was established using data from Table 4-1 in Bendt et al. [14]. This table contains data from the 16 standard Lawrence Berkeley Laboratory (LBL) circumsolar irradiance scans [26], namely, σ s u n and CSR, and its relationship is shown in Figure 3.
Then, a fit to the data shown in Figure 4–1a of Bendt et al. [14] was performed, enabling the determination of the intercept factor according to CSR using the procedure described below. Firstly, σ o p t was fixed at 10 mrad and σ s u n varied from 2.5 mrad to 10 mrad. In this way, it was possible to determine the intercept factor for various combinations of concentration factors (C) and rim angles, as shown in Figure 4 for a given parabolic trough system.
From Figure 4, it is clear that higher concentration ratios imply a larger effect of the CSR on the intercept factor of a parabolic trough system, regardless of the selected rim angle. Higher CSR values stand for higher dispersion of reflected solar rays from the parabolic through the concentrator to the receiver, resulting in a lower intercept factor.
For lower rim angles, the CSR seems to have a larger impact on the intercept factor than for higher rim angles, for the same concentration ratio. This has to do with the aperture angle of the reflector, because lower rim angles correspond to lower acceptance angles of the reflector. Therefore, lower reflector aperture angles result in less irradiance reflected when the circumsolar region is larger (i.e., higher CSR).
Regarding the mitigation of CSR effects on the intercept factor of parabolic trough systems, it appears that using higher rim angles is the best approach based on the results from Figure 4. It is expected that the same conclusions could be draw for other CSP systems. For instance, Rabl and Bendt [27] found that the intercept factor of a parabolic dish strongly depended on the rim angle. They also stated that concentrators with a rim angle of 30º are two times more sensitive to CSNI variation than those with a rim angle of 60º. However, those systems are not addressed here for conciseness.

5. Conclusions

In this work, a model to predict circumsolar normal irradiance (CSNI) for several half-opening angles was developed. The model is based on the CSR model developed by Abreu et al. [6] and uses a polynomial fitting of the model parameters to estimate CSNI for half-opening angles of interest for CSP systems, such as parabolic trough and parabolic disk. It was found that the proposed model performs significantly better at the locations under study here than the other models available in the literature. It was also found that the local aerosol regime and atmospheric conditions have a higher impact on model fitting and performance than the overall climate zone of the locations under study.
Regarding the impact of CSNI in CSP systems, it was found that discarding an accurate CSNI estimate could lead to up to a 7% difference between the DNI that is measured by a pyrheliometer and the DNI that is effectively captured by the system. It was also found that these differences are greater for lower half-opening angles, as expected, since the difference between the half-opening angle of the pyrheliometer and the half-opening angle of the CSP system is greater.
The impact of CSNI in the operation of a CSP system, namely with parabolic trough concentrators, was also studied in this work. The authors found that higher CSR values lead to lower intercept factors (the ratio between the irradiance reaching the receiver over the incident irradiance) for these systems. It was also found that if parabolic trough designers aim to reduce the impact of CSNI on the intercept factor, then parabolic troughs with higher rim angles are preferred.
This study also found that further study on the direct impact of the different aerosol and atmospheric characteristics on the development of local or global CSNI models, as well as on the evaluation of the impact of CSNI in other CSP systems, is needed and should be addressed in future work.

Author Contributions

Conceptualization, E.F.M.A., P.C. and M.J.C.; Methodology, E.F.M.A. and P.C.; validation, E.F.M.A.; writing—original draft preparation, E.F.M.A.; writing—review and editing, P.C. and M.J.C.; supervision, P.C. and M.J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by FCT (the Portuguese Science and Technology Foundation), grant number SFRH/BD/136433/2018, co-funded by National funds through FCT—Fundação para a Ciência e Tecnologia, I.P., projects UIDB/04683/2020 and UIDP/04683/2020, and it was also funded by the European Union through the European Regional Development Fund, included in the COMPETE 2020 (Operational Program Competitiveness and Internationalization) through the project DNI-Alentejo with reference ALT20-03-0145-FEDER-000011. The APC was funded by the Institute of Earth Sciences.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

Table A1. Statistical analysis of circumsolar models for a half-opening value of 1.0°. Legend: TR—Tropical; TM—Temperate; AR—Arid. Values in bold indicate the best performing model according to each statistical indicator.
Table A1. Statistical analysis of circumsolar models for a half-opening value of 1.0°. Legend: TR—Tropical; TM—Temperate; AR—Arid. Values in bold indicate the best performing model according to each statistical indicator.
StationStatistical IndicatorModels
Eissa et al. [5] SOVEissa et al. [5] TAMEissa et al. [5] CombinedThis Work
DAR (TR)rMBE (%)−8.28−16.53−12.94−2.25
rRMSE (%)74.3476.3875.2364.85
R0.46160.45160.45850.6277
FB0.0647−0.01370.01800.0694
FGE0.39660.39790.39460.3282
EVR (TM)rMBE (%)−17.47−24.30−21.61−4.13
rRMSE (%)68.0672.6070.2848.07
R0.75150.74750.75060.8641
FB0.0549−0.01160.00840.0634
FGE0.45750.46600.45770.3657
SMS (TM)rMBE (%)51.2341.2344.25−7.88
rRMSE (%)94.9187.5390.0359.32
R0.59440.59520.59630.7761
FB0.59350.53990.5512−0.0536
FGE0.66570.63210.63520.4101
GOB (AR)rMBE (%)−9.66−16.54−13.970.40
rRMSE (%)76.1678.5877.2963.26
R0.59630.59050.59500.7356
FB0.12610.06260.08120.1404
FGE0.45740.45420.45060.4007
TAM (AR)rMBE (%)−46.42−51.60−49.26−11.08
rRMSE (%)64.1968.5266.4835.98
R0.65870.65910.65960.8105
FB−0.6053−0.6832−0.6521−0.0800
FGE0.61790.68820.65970.2596
Table A2. Statistical analysis of circumsolar models for a half-opening value of 1.6°. Legend: TR—Tropical; TM—Temperate; AR—Arid. Values in bold indicate the best performing model according to each statistical indicator.
Table A2. Statistical analysis of circumsolar models for a half-opening value of 1.6°. Legend: TR—Tropical; TM—Temperate; AR—Arid. Values in bold indicate the best performing model according to each statistical indicator.
StationStatistical IndicatorModels
Eissa et al. [5] SOVEissa et al. [5] TAMEissa et al. [5] CombinedThis Work
DAR (TR)rMBE (%)−7.32−14.46−11.38−2.97
rRMSE (%)71.5273.2472.2661.41
R0.48840.47920.48560.6585
FB0.0677−0.00040.02620.0495
FGE0.38610.38650.38360.3096
EVR (TM)rMBE (%)−18.22−24.22−21.90−6.00
rRMSE (%)69.0073.0770.9647.49
R0.76510.76150.76440.8773
FB0.0501−0.01260.00500.0379
FGE0.44790.45380.44700.3495
SMS (TM)rMBE (%)55.7146.3048.97−8.79
rRMSE (%)97.6590.2592.7558.16
R0.60260.60600.60580.7872
FB0.61280.55990.5702−0.0604
FGE0.67480.63590.64090.3856
GOB (AR)rMBE (%)−9.47−15.53−13.340.35
rRMSE (%)74.0176.0074.9360.11
R0.60990.60580.60950.7565
FB0.12080.06140.07760.1306
FGE0.44570.43890.43670.3819
TAM (AR)rMBE (%)−47.77−52.19−50.16−13.12
rRMSE (%)65.6169.3867.5636.73
R0.67360.67440.67480.8230
FB−0.6253−0.6952−0.6679−0.1024
FGE0.63340.69810.67280.2593
Table A3. Statistical analysis of circumsolar models for a half-opening value of 1.8°. Legend: TR—Tropical; TM—Temperate; AR—Arid. Values in bold indicate the best performing model according to each statistical indicator.
Table A3. Statistical analysis of circumsolar models for a half-opening value of 1.8°. Legend: TR—Tropical; TM—Temperate; AR—Arid. Values in bold indicate the best performing model according to each statistical indicator.
StationStatistical IndicatorModels
Eissa et al. [5] SOVEissa et al. [5] TAMEissa et al. [5] CombinedThis Work
DAR (TR)rMBE (%)−11.28−17.81−15.01−3.11
rRMSE (%)69.4771.4970.3958.15
R0.50530.49690.50290.6813
FB0.0232−0.0431−0.01760.0445
FGE0.38090.38560.38150.3012
EVR (TM)rMBE (%)−22.25−27.84−25.70−6.46
rRMSE (%)70.1174.1972.1046.61
R0.77470.77170.77440.8822
FB0.0009−0.0634−0.04560.0261
FGE0.44340.45180.44450.3420
SMS (TM)rMBE (%)50.2541.0643.62−8.91
rRMSE (%)92.8486.0188.3156.76
R0.60650.61180.61060.7919
FB0.57660.51980.5312−0.0668
FGE0.64540.60360.60940.3789
GOB (AR)rMBE (%)−13.74−19.40−17.390.01
rRMSE (%)74.1076.2475.1459.07
R0.61380.61100.61400.7617
FB0.07190.01090.02730.1230
FGE0.43100.42700.42440.3730
TAM (AR)rMBE (%)−50.66−54.66−52.82−13.94
rRMSE (%)68.0471.5669.8537.16
R0.67780.67870.67920.8243
FB−0.6735−0.7412−0.7149−0.1121
FGE0.67720.74260.71710.2589

References

  1. Polo, J.; Fernández-Peruchena, C.; Gastón, M. Analysis on the Long-Term Relationship between DNI and CSP Yield Production for Different Technologies. Sol. Energy 2017, 155, 1121–1129. [Google Scholar] [CrossRef]
  2. Blanc, P.; Espinar, B.; Geuder, N.; Gueymard, C.; Meyer, R.; Pitz-Paal, R.; Reinhardt, B.; Renné, D.; Sengupta, M.; Wald, L.; et al. Direct Normal Irradiance Related Definitions and Applications: The Circumsolar Issue. Sol. Energy 2014, 110, 561–577. [Google Scholar] [CrossRef]
  3. World Meteorological Organization (WMO). Guide to Instruments and Methods of Observation. Volume I—Measurement of Meteorological Variables; World Meteorological Organization (WMO): Geneva, Switzerland, 2021. [Google Scholar]
  4. Rabl, A. Comparison of Solar Concentrators. Sol. Energy 1976, 18, 93–111. [Google Scholar] [CrossRef]
  5. Eissa, Y.; Blanc, P.; Ghedira, H.; Oumbe, A.; Wald, L. A Fast and Simple Model to Estimate the Contribution of the Circumsolar Irradiance to Measured Broadband Beam Irradiance under Cloud-Free Conditions in Desert Environment. Sol. Energy 2018, 163, 497–509. [Google Scholar] [CrossRef]
  6. Abreu, E.F.M.; Canhoto, P.; Costa, M.J. Development of a Clear-Sky Model to Determine Circumsolar Irradiance Using Widely Available Solar Radiation Data. Sol. Energy 2020, 205, 88–101. [Google Scholar] [CrossRef]
  7. Abreu, E.F.M.; Canhoto, P.; Costa, M.J. Circumsolar Irradiance Modelling Using libRadtran and AERONET Data. AIP Conf. Proc. 2019, 2126, 190001. [Google Scholar]
  8. Mayer, B.; Kylling, A. Technical Note: The libRadtran Software Package for Radiative Transfer Calculations—Description and Examples of Use. Atmos. Chem. Phys. 2005, 23, 1855–1877. [Google Scholar] [CrossRef]
  9. Ivanova, S.M. Modelling of Circumsolar Direct Irradiance as a Function of GHI, DHI and DNI. In Proceedings of the 8th International Conference on Solar Radiation and Daylighting Solaris-2017, London, UK, 27–28 July 2017. [Google Scholar]
  10. Neumann, A.; von der Au, B.; Heller, P. Measurements of Circumsolar Radiation at the Plataforma Solar (Spain) and at DLR (Germany); ASME International Solar Energy Conference: Albuquerque, NM, USA, 1998; pp. 429–438. [Google Scholar]
  11. Forman, P.; Penkert, S.; Kämper, C.; Stallmann, T.; Mark, P.; Schnell, J. A Survey of Solar Concrete Shell Collectors for Parabolic Troughs. Renew. Sustain. Energy Rev. 2020, 134, 110331. [Google Scholar] [CrossRef]
  12. Belhomme, B.; Pitz-Paal, R.; Schwarzbözl, P.; Ulmer, S. A New Fast Ray Tracing Tool for High-Precision Simulation of Heliostat Fields. J. Sol. Energy Eng. 2009, 131, 031002. [Google Scholar] [CrossRef]
  13. Leary, P.L.; Hankins, J.D. User’s Guide for MIRVAL: A Computer Code for Comparing Designs of Heliostat-Receiver Optics for Central Receiver Solar Power Plants; Sandia National Lab. (SNL-CA): Livermore, CA, USA, 1979. [Google Scholar]
  14. Bendt, P.; Rabl, A.; Gaul, H.W.; Reed, K.A. Optical Analysis and Optimization of Line Focus Solar Collectors; Solar Energy Research Institution: Golden, CO, USA, 1979. [Google Scholar]
  15. Schwarzbözl, P.; Pitz-Paal, R.; Schmitz, M. Visual HFLCAL—A Software Tool for Layout and Optimisation of Heliostat Fields. In Proceedings of the SolarPACES 2009, Berlin, Germany, 15–18 September 2009. [Google Scholar]
  16. Gilman, P.; Blair, N.; Mehos, M.; Christensen, C.; Janzou, S.; Cameron, C. Solar Advisor Model User Guide for Version 2.0; National Renewable Energy Laboratory: Colorado, CO, USA, 2008; NREL/TP-670-43704, 937349. [Google Scholar]
  17. Quaschning, V.; Ortmanns, W.; Kistner, R.; Geyer, M. Greenius: A New Simulation Environment for Technical and Economical Analysis of Renewable Independent Power Projects. In Proceedings of the SED2001, Solar Engineering 2001: (FORUM 2001: Solar Energy—The Power to Choose), Washington, DC, USA, 21 April 2001; pp. 413–417. [Google Scholar]
  18. Dersch, J.; Schwarzbözl, P.; Richert, T. Annual Yield Analysis of Solar Tower Power Plants With GREENIUS. J. Sol. Energy Eng. 2011, 133, 031017. [Google Scholar] [CrossRef]
  19. Holben, B.N.; Eck, T.F.; Slutsker, I.; Tanré, D.; Buis, J.P.; Setzer, A.; Vermote, E.; Reagan, J.A.; Kaufman, Y.J.; Nakajima, T.; et al. AERONET—A Federated Instrument Network and Data Archive for Aerosol Characterization. Remote Sens. Environ. 1998, 66, 1–16. [Google Scholar] [CrossRef]
  20. Driemel, A.; Augustine, J.; Behrens, K.; Colle, S.; Cox, C.; Cuevas-Agulló, E.; Denn, F.M.; Duprat, T.; Fukuda, M.; Grobe, H.; et al. Baseline Surface Radiation Network (BSRN): Structure and Data description (1992–2017). Earth Syst. Sci. Data 2018, 10, 1491–1501. [Google Scholar] [CrossRef]
  21. Giles, D.M.; Sinyuk, A.; Sorokin, M.G.; Schafer, J.S.; Smirnov, A.; Slutsker, I.; Eck, T.F.; Holben, B.N.; Lewis, J.R.; Campbell, J.R.; et al. Advancements in the Aerosol Robotic Network (AERONET) Version 3 Database—Automated near-Real-Time Quality Control Algorithm with Improved Cloud Screening for Sun Photometer Aerosol Optical Depth (AOD) Measurements. Atmos. Meas. Tech. 2019, 12, 169–209. [Google Scholar] [CrossRef]
  22. Sinyuk, A.; Holben, B.N.; Eck, T.F.; Giles, D.M.; Slutsker, I.; Korkin, S.; Schafer, J.S.; Smirnov, A.; Sorokin, M.; Lyapustin, A. The AERONET Version 3 Aerosol Retrieval Algorithm, Associated Uncertainties and Comparisons to Version 2. Atmos. Meas. Tech. 2020, 13, 3375–3411. [Google Scholar] [CrossRef]
  23. Liu, B.Y.H.; Jordan, R.C. The Interrelationship and Characteristic Distribution of Direct, Diffuse and Total Solar Radiation. Sol. Energy 1960, 4, 1–19. [Google Scholar] [CrossRef]
  24. Gueymard, C.A.; Ruiz-Arias, J.A. Extensive Worldwide Validation and Climate Sensitivity Analysis of Direct Irradiance Predictions from 1-Min Global Irradiance. Sol. Energy 2016, 128, 1–30. [Google Scholar] [CrossRef]
  25. Gueymard, C.A.; Bright, J.M.; Lingfors, D.; Habte, A.; Sengupta, M. A Posteriori Clear-Sky Identification Methods in Solar Irradiance Time Series: Review and Preliminary Validation Using Sky Imagers. Renew. Sustain. Energy Rev. 2019, 109, 412–427. [Google Scholar] [CrossRef]
  26. Noring, J.E.; Grether, D.F.; Hunt, A.J. Circumsolar Radiation Data: The Lawrence Berkeley Laboratory Reduced Data Base. NASA STI/Recon Tech. Rep. N 1991, 92, 19628. [Google Scholar]
  27. Rabl, A.; Bendt, P. Effect of Circumsolar Radiation on Performance of Focusing Collectors. J. Sol. Energy Eng. 1982, 104, 237–250. [Google Scholar] [CrossRef]
Figure 1. Average difference between D N I ( α p y r ) and D N I ( α ) .
Figure 1. Average difference between D N I ( α p y r ) and D N I ( α ) .
Energies 16 07950 g001
Figure 2. Boxplots for aerosol optical depth (AOD) and single scattering albedo (SSA) for the stations analyzed in this work.
Figure 2. Boxplots for aerosol optical depth (AOD) and single scattering albedo (SSA) for the stations analyzed in this work.
Energies 16 07950 g002
Figure 3. Relationship between CSR and σ s u n .
Figure 3. Relationship between CSR and σ s u n .
Energies 16 07950 g003
Figure 4. Variation of intercept factor according to CSR for a parabolic trough system with rim angle of 45° (top), 75° (middle), and 90° (bottom).
Figure 4. Variation of intercept factor according to CSR for a parabolic trough system with rim angle of 45° (top), 75° (middle), and 90° (bottom).
Energies 16 07950 g004
Table 1. Information on the AERONET and radiometric stations and data period used in this work. Legend: Lat.—Latitude, Long.—Longitude, Alt.—Altitude, TR—Tropical, TM—Temperate, AR—Arid.
Table 1. Information on the AERONET and radiometric stations and data period used in this work. Legend: Lat.—Latitude, Long.—Longitude, Alt.—Altitude, TR—Tropical, TM—Temperate, AR—Arid.
LocationCodeLat.
(°N)
Long.
(°E)
Alt.
(m)
Climate
Zone
PeriodData Points
Darwin, AUSDAR−12.425130.83130TR2012–2014974
Évora, PRTEVR38.568−7.912293TM2015–20171163
Gobabeb, NAMGOB−23.56115.042407AR2015–20171833
S. M. da Serra, BRASMS−29.443−53.823489TM2014–2016614
Tamanrasset, DZATAM22.7905.5291385AR2014–20161060
Table 2. Polynomial coefficients and respective coefficients of determination of model parameters a, b , and c fitting as a function of alpha. Legend: TR—Tropical; TM—Temperate; AR—Arid.
Table 2. Polynomial coefficients and respective coefficients of determination of model parameters a, b , and c fitting as a function of alpha. Legend: TR—Tropical; TM—Temperate; AR—Arid.
StationParameterPolynomial CoefficientsR2
C 0 C 1 C 3
DAR (TR)a−0.00940.013700.9978
b0.7667−0.07270.01520.4325
c1.4199−0.00630.00450.3446
EVR (TM)a−0.01520.020600.9962
b0.74030.00330.00560.8524
c1.67700.0343−0.00170.8817
SMS (TM)a−0.01580.021300.9931
b2.04060.1906−0.02270.3278
c2.19550.0439−0.01210.2684
GOB (AR)a−0.02380.030300.9938
b1.25460.02860.01280.9250
c1.74570.01610.00160.6406
TAM (AR)a−0.06980.074000.9801
b1.48900.09440.04270.9919
c1.40510.07390.00240.9359
Table 3. Statistical analysis of circumsolar models for a half-opening value of 0.8°. Legend: TR—Tropical; TM—Temperate; AR—Arid. Values in bold indicate the best performing model according to each statistical indicator.
Table 3. Statistical analysis of circumsolar models for a half-opening value of 0.8°. Legend: TR—Tropical; TM—Temperate; AR—Arid. Values in bold indicate the best performing model according to each statistical indicator.
StationStatistical IndicatorModels
Eissa et al. [5] SOVEissa et al. [5] TAMEissa et al. [5] CombinedThis Work
DAR (TR)rMBE (%)−3.69−12.79−8.82−1.96
rRMSE (%)73.6375.1874.2165.51
R0.45580.44620.45290.6098
FB0.11110.02800.06190.0783
FGE0.40590.39860.39840.3373
EVR (TM)rMBE (%)−14.56−21.99−19.06−5.38
rRMSE (%)69.7674.1471.8651.53
R0.73720.73330.73630.8503
FB0.09600.02560.04730.0501
FGE0.47200.47710.47000.3787
SMS (TM)rMBE (%)54.8643.9247.29−7.90
rRMSE (%)99.2891.2993.9862.40
R0.58860.58870.59020.7664
FB0.62480.56830.5809−0.0267
FGE0.69400.65700.66160.4119
GOB (AR)rMBE (%)−5.53−13.16−10.290.70
rRMSE (%)77.7879.8978.6966.07
R0.58570.58000.58450.7205
FB0.17360.10590.12640.1519
FGE0.47390.46570.46370.4149
TAM (AR)rMBE (%)−43.20−48.86−46.32−10.16
rRMSE (%)63.0167.4565.3637.98
R0.63650.63600.63700.7870
FB−0.5482−0.6302−0.5975−0.0651
FGE0.57740.64400.61710.2760
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Abreu, E.F.M.; Canhoto, P.; Costa, M.J. Prediction of Circumsolar Irradiance and Its Impact on CSP Systems under Clear Skies. Energies 2023, 16, 7950. https://doi.org/10.3390/en16247950

AMA Style

Abreu EFM, Canhoto P, Costa MJ. Prediction of Circumsolar Irradiance and Its Impact on CSP Systems under Clear Skies. Energies. 2023; 16(24):7950. https://doi.org/10.3390/en16247950

Chicago/Turabian Style

Abreu, Edgar F.M., Paulo Canhoto, and Maria João Costa. 2023. "Prediction of Circumsolar Irradiance and Its Impact on CSP Systems under Clear Skies" Energies 16, no. 24: 7950. https://doi.org/10.3390/en16247950

APA Style

Abreu, E. F. M., Canhoto, P., & Costa, M. J. (2023). Prediction of Circumsolar Irradiance and Its Impact on CSP Systems under Clear Skies. Energies, 16(24), 7950. https://doi.org/10.3390/en16247950

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop