Next Article in Journal
Power System Restoration Planning Strategy Based on Optimal Energizing Time of Sectionalizing Islands
Next Article in Special Issue
The Influence of Local Environmental, Economic and Social Variables on the Spatial Distribution of Photovoltaic Applications across China’s Urban Areas
Previous Article in Journal
Active-Reactive Additional Damping Control of a Doubly-Fed Induction Generator Based on Active Disturbance Rejection Control
Previous Article in Special Issue
Simulation and Experimental Study on the Optical Performance of a Fixed-Focus Fresnel Lens Solar Concentrator Using Polar-Axis Tracking
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comprehensive Approach for Modelling Horizontal Diffuse Radiation, Direct Normal Irradiance and Total Tilted Solar Radiation Based on Global Radiation under Danish Climate Conditions

Department of Civil Engineering, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
*
Author to whom correspondence should be addressed.
Energies 2018, 11(5), 1315; https://doi.org/10.3390/en11051315
Submission received: 17 April 2018 / Revised: 3 May 2018 / Accepted: 5 May 2018 / Published: 22 May 2018
(This article belongs to the Special Issue Solar Thermal Energy Utilization Technologies in Buildings)

Abstract

:
A novel combined solar heating plant with flat plate collectors (FPC) and parabolic trough collectors (PTC) was constructed and put into operation in Taars, 30 km north of Aalborg, Denmark in August 2015. To assess the thermal performance of the solar heating plant, global radiation, direct normal irradiance (DNI) and total radiation on the tilted collector plane of the flat plate collector field were measured. To determine the accuracy of the measurements, the calculated solar radiations, including horizontal diffuse radiation, DNI and total tilted solar radiation with seven empirical models, were compared each month based on an hourly time step. In addition, the split of measured global radiation into diffuse and beam radiation based on a model developed by DTU (Technical University of Denmark) and the Reduced Reindl correlation model was investigated. A new method of combining empirical models, only based on measured global radiation, was proposed for estimating hourly total radiation on tilted surfaces. The results showed that the DTU model could be used to calculate diffuse radiation on the horizontal surface, and that the anisotropic models (Perez I and Perez II) were the most accurate for calculation of total radiation on tilted collector surfaces based only on global radiation under Danish climate conditions. The proposed method was used to determine reliable horizontal diffuse radiation, DNI and total tilted radiation with only the measurement of global radiation. Only a small difference compared to measured data, was found. The proposed method was cost-effective and needed fewer measurements to obtain reliable DNI and total radiation on the tilted plane. This method may be extended to other Nordic areas that have similar weather.

1. Introduction

Energy consumption in the building sector accounts for about 40% of society’s energy consumption in developed countries. Using renewable energy, especially solar energy, for heating and cooling in the building sector is a promising way to reduce the fossil energy consumption of buildings [1,2]. Solar thermal energy is one of the most commercial renewable energies in the building sector [3,4]. Large solar heating plants connected to district heating networks have been of great success in Europe, especially in Denmark. Most large scale solar heating plants in Europe, even worldwide, are constructed in Denmark. Denmark is the first and the only country with a mature commercial market for solar district heating plants. By the end of 2016, more than 1.3 million m2 solar heating plants were in operation in Denmark [5]. Real-time solar radiation data is widely used as a basic input to control large-scale solar collector fields across thousands of square meters during their lifetime. Furthermore, accurate solar radiation data are very important for designing solar heating systems and estimating the thermal performance of solar district heating plants. Compared to global irradiance, the direct beam component shows much more variability in space and time. Global radiation split into beams and diffuse radiation on the collector plane is important for evaluation of the performance of different collector types and collector field designs as well. In the past, in most cases, inexpensive and inaccurate solar radiation sensors have been used to measure solar radiation on collector planes in solar district heating plants in Denmark. Few technicians onsite in solar heating plants have paid much attention to the accuracy of measurements about solar radiation. Poor solar radiation may result in the wrong control strategies for such large-scale solar collector fields, which can influence the cost-performance of solar heating plants significantly. A simple model for accurately modelling solar radiation is needed to double-check solar radiation measurements onsite, in a cost-effective and fast way.

1.1. State of the Art

Generally, climate stations measure global radiation and, only in rare cases, accurate DNI or diffuse solar radiation on the horizontal surface. Therefore, total irradiation on tilted surfaces and DNI in most cases is calculated using measured global irradiation by means of empirical models for general use. Shukla et al. [6] carried out a comparative study of isotropic and anisotropic sky models to estimate solar radiation incidence on tilted surfaces in India. Demain et al. [7] evaluated 14 empirical models to predict global radiation on inclined surfaces. A hybrid model from the coupling of three models under different sky conditions have been developed for Belgium. Khorasanizadeh et al. set up a new diffuse solar radiation model to determine the optimum tilt angle of surfaces in Tabass, Iran [8]. Marques Filho et al. carried out observational charaterisation and empirical modelling of global radiation, diffuse and direct solar radiation on surfaces in the city of Rio de Janeiro [9]. El Mghouchi et al. evaluated four empirical models to predict daily direct diffuse and global radiation in Tutuan city, north of Morocco [10]. Jakhrani et al. investigated the accuracy of different empirical models for calculating total solar radiation on tilted surfaces [11]. It was found that the isotopic model (Liu and Jordan model) was better for the prediction of solar energy radiation in cloudy weather conditions and could be used to calculate available solar radiation on tilted surfaces in overcast skies under Malaysian climate conditions. El-Sebaii et al. also calculated diffuse radiation on horizontal surfaces and total solar radiation on tilted surfaces using empirical models [12]. They [12] also found that the isotropic model (Liu and Jordan model) could be used to calculate total radiation on tilted surfaces with good accuracy in Jeddah, Saudi Arabia. Gopinathan K.K. investigated solar radiation on variously oriented sloping surfaces in Lesotho, South Africa, with the isotropic model [13]. Li et al. carried out estimation of daily global solar radiation in China [14]. Alyahya et al. analysed the new solar radiation Atlas for Saudi Arabia [15]. Bird et al. developed a simple solar spectral model for direct and diffuse irradiance on horizontal and tilted planes on the earth’s surface for cloudless atmospheres [16]. There have also been several studies on the prediction of solar radiation using machine learning and multivariable regression methods [17,18]. Despotovic et al. [19] investigated the accuracy of different empirical models in predicting total tilted solar radiation and diffuse horizontal solar radiation, respectively. Ineichen P. concluded that the Perez model is slightly better (in terms of RMSD) than other models in any case, even with synthetic data [20]. Gueymard et al. [21] carried out a comprehensive evaluation study of the performance of 140 separation models selected from the literature to predict direct normal irradiance from global horizontal irradiance. The evaluation was based on measured, high-quality, 1-min data of global horizontal irradiance and DNI at 54 research-class stations from seven continents. Only two models consistently delivered the best predictions over the arid, temperate and tropical zones and no model performed consistently well over the high-albedo zone. A comparative study of the impact of horizontal-to-tilted solar irradiance conversion in modelling small PV array performance was presented in [22]. A neural network model was employed to predict daily direct solar radiation in [23]. Frydrychowicz-Jastrzębska et al. compared selected isotropic and anisotropic mathematical models to calculate the distribution of solar radiation on the photovoltaic module plane with any spatial orientation for Poland [24].
Mubarak et al. [25] compared five empirical models for PV applications. The authors concluded that the models of Hay and Davies and Reindl are recommended to estimate tilted irradiance for south-facing modules in regions with mainly cloudy conditions and when albedo measurements are not available. The Hay and Davies model is useful for vertical surfaces (e.g., facades and glazing), whereas the Perez model is recommended for sunny sites and when albedo measurements are available. Lee et al. [26] investigated solar radiation models to estimate direct normal irradiance for Korea. The Reindl-2 model was selected as the best among the evaluated ten existing models for Korea. Different conclusions can be drawn for different locations. Using previous empirical models to convert global solar radiation data for general use in high latitude areas, such as Denmark, may not give highly accurate results. Furthermore, limited literature was found on the analysis and prediction of total tilted solar radiation and DNI at high latitudes. A novel combined solar heating plant with a 4039 m2 parabolic trough collector field and a 5960 m2 flat plate collector field in Taars was put into operation in August 2015 [27]. To evaluate the thermal performance of the plant and the accuracy of the calculated solar radiation, total tilted and horizontal solar radiation was measured in the collector field. In addition, a weather station was in operation close to the solar collector fields to ensure that the pyranometers in the plant had correct values to reduce systemic errors and to measure direct normal irradiance (DNI).
Diffuse radiation influences the thermal performance of the flat plate collector field. In this study, diffuse horizontal radiation was estimated using the RR model (Reduced Reindl correlation model) [28] and the DTU (Technical University of Denmark) model [29]. The DTU model was developed based on measurements from 2006–2010 at a climate station at DTU [29] and was used in this paper to calculate diffuse radiation on the horizontal surface with only global radiation as an input. The RR model was developed by Reindl in 1990 for general use to calculate diffuse radiation on the horizontal surface with only global radiation as an input [30]. These two models were compared to the measured data from the Taars plant.
When diffuse radiation on the horizontal surface has been calculated, direct radiation on the same surface can be derived by the subtraction of beam radiation from global radiation. DNI can then be determined by dividing by the cosine of the zenith angle indirectly. The last two steps for direct radiation are exact numerical conversions without calculation error.

1.2. Scope and Objective

The novel contributions of this paper are as follows: (1) A solar radiation model developed for the Danish climate conditions in the Nordic area was validated; (2) The measured data were from a pilot solar heating plant, not a laboratory, which is more practical and the whole chain of calculations, up to the long term performance of the solar collectors, can be validated; and (3) DNI, diffuse radiation on the horizontal surface and total radiation on the tilted surface during the whole year were analysed.
This paper validated the performance of the DTU model for the derivation of horizontal diffuse irradiance and beam radiation based on more widely available global horizontal irradiance data under Danish climate conditions. One isotropic model and four anisotropic models for the calculation of total titled radiation were also investigated. The difference between measured solar radiation and modelled solar radiation estimated by the empirical formulas under Danish climate conditions, including DNI, diffuse horizontal radiation and total titled solar radiation, were shown. MBE, RMSE, MAPE and RPE were used to assess the feasibility of the investigated empirical models under Danish climate conditions. Calculated total tilted radiation only based on global horizontal radiation and on both global horizontal radiation and beam radiation were discussed, which could provide a new method to calculate total tilted radiation with less measurements under Danish climate conditions and may be extended to other Nordic areas that have similar weather. The combined method to calculate total tilted solar radiation could be a useful tool for design large-scale solar district heating plants, which have been of great success in Denmark.
The aim of this article was to develop a solar radiation model to predict DNI and total tilted solar radiation accurately for solar thermal systems under Danish climate conditions.
The structure of the article is summarised as follows: Section 1 is the introduction; Section 2 is the introduction of the measurements; Section 3 and Section 4 present the method and detailed empirical models used in this study; Section 5 shows the validation of the empirical models; Section 6 shows the predicted DNI and total tilted solar radiation of the selected models; Section 7 is the discussion; and Section 8 is the conclusion.

2. Data Collection and Location Description

As is shown in Figure 1, Denmark has six solar radiation zones with different yearly global radiation, around 1000–1200 kWh/m2. The Taars plant is located in the first solar radiation zone, in the northern part of Denmark [31,32,33]. Figure 2 illustrates the locations of the weather station and the pyranometers in the flat plate collector field. The weather station is next to a solar heating plant. There are several pyranometers to measure global solar radiation and total radiation on the tilted plane of the flat plate collectors in the middle of flat plate collector field (Figure 1). The latitude of Taars is 57.39°N and the longitude is 10.11°E.
As is shown in Figure 2 and Figure 3, four south facing pyranometers with a tilt 50° were installed on the top of a flat plate collector plane in the middle of the flat plate collector field [35]. Three Apogee Pyranometer SP-110 are used as backup sensors to double check the measured total radiation, as is shown in the Figure 3 left. Two of the pyranometers to measure solar radiation on the horizontal surface and solar radiation on the titled collector plane were Kipp & Zonen SMP11 (see Figure 3), which are used in this study. DNI was measured with a PMO6-CC pyrheliometer with the sun tracking platform Sunscanner SC1 in the weather station next to the solar heating plant (see Figure 2 and Figure 4). Table 1 and Table 2 show the technical specifications of the Kipp & Zonen SMP11 pyranometer and PMO6-CC pyrheliometer [36,37]. The measurements, including global radiation, total tilted solar radiation and DNI, were recorded in 2 minutes intervals from the middle of August 2015.

3. Methodology

Figure 5a gives a schematic illustration of this study. Firstly, the DTU model and RR models were used to calculate horizontal diffuse radiation based on measured global radiation. Five calculation models for total radiation on tilted surfaces for general use were investigated: one isotropic model and four anisotropic models. Circumsolar diffuse and horizon-brightening components on the tilted surfaces were taken into consideration in the anisotropic models, but not in the isotropic model. Calculated total tilted solar radiation of the empirical models based on measured global radiation and DNI was derived. Measured total tilted radiation was used to evaluate the suitability of the empirical models for measuring total tilted radiation in Danish conditions (see the validation cycle in Figure 5a), as elaborated on in Section 5.
Then, the selected empirical models based on calculated diffuse radiation and beam radiation were employed to calculate total tilted radiation (Figure 5b), as described further in Section 6. Calculated total tilted radiation using the DTU model and the investigated empirical models (Perez models) only based on global radiation showed good agreement with the measured values from September 2015 to August 2016 (Figure 5b). DNI calculated by the DTU model also had good agreement with measured DNI. In summary, the proposed method to calculate total tilted solar radiation only based on measured global horizontal radiation (red flow chart) is a new, simple and cost-effective approach to obtain accurate total tilted solar radiation for Danish conditions, as measuredglobal radiation data is always available from climate stations. Furthermore, DNI and diffuse radiation measurements are relatively costly both in terms of equipment and manpower. Accurate long-term data for these variables are seldom available in most cases. Therefore, accurate calculated DNI, diffuse radiation and total tilted radiation based only on measured global radiation using the method proposed in this paper is very valuable.

4. Empirical Models

DNI, global radiation and total tilted solar radiation on the top of a 50° tilted, south facing collector were measured with a high time resolution of 2 min. Hourly mean values were calculated based on the measured values. Calculated solar radiation in this study was based on the mean data of a 1 h time step. Both the DTU model and the RR model were used to calculate diffuse radiation on the horizontal surface. Five other empirical models (one isotropic model and four anisotropic models) were used to calculate total solar radiation on the tilted surface. Ground reflectance or albedo was assumed to be 0.1. This value was a reasonable estimation of ground reflectance when shadows between collector rows in the solar heating plant were considered.

4.1. Measured Horizontal Diffuse Radiation

Diffuse radiation on the horizontal surface was not measured directly in the Taars plant. However, diffuse radiation can be derived accurately as the difference between total radiation and beam radiation. Measured beam radiation was calculated by measured DNI and solar zenith angle using Equation (1). Measured diffuse radiation on the horizontal surface was determined as the difference between the measured global radiation and beam radiation components, indirectly, using Equation (2).
G b = DNI × cos θ z
G d = G G b

4.2. Modelled Horizontal Diffuse Radiation

(1)
DTU model
Dragsted et al. measured and analysed solar radiation from a climate station at the Technical University of Denmark from 2006 to 2010, and developed an empirical model to calculate horizontal diffuse radiation from global radiation on the horizontal surface for Danish climate conditions [29]. The empirical model is as follows in Equations (3)–(7):
K T = G / G 0
G d / G = 0.60921 K T 3 + 1.9982 K T 2 0.2787 K T + 1 ,   0.00 K T < 0.29
G d / G = 3.99 K T 3 7.1469 K T 2 + 2.3996 K T + 0.746 ,   0.29 K T < 0.72
G d / G = 288.63 K T 4 625.26 K T 3 + 448.06 K T 2 105.84 K T ,   0.72 K T < 0.80
G d / G = 65.89 K T 4 210.69 K T 3 + 222.91 K T 2 77.203 K T ,   0.80 K T < 1.20
(2)
Reduced Reindl correlation model
The Reduced Reindl correlation model is based on the relationships developed by Reindl et al. [30]. The Reduced Reindl model uses clearness index and solar altitude angle to estimate diffuse radiation on the horizontal surface. The correlation is given by Equations (8)–(10):
G d / G = 1.020 0.254 K T + 0.0123 sin α , 0 K T 0.3 , G d / G 1.0
G d / G = 1.400 1.794 K T + 0.177 sin α , 0.3 < K T < 0.78 , 0.1 G d / G 0.97
G d / G = 0.486 K T 0.182 sin α , 0.78 K T , 0.1 G d / G

4.3. Modelled Total Tilted Solar Radiation

(1)
Isotropic model
The typical isotropic model was developed by Liu and Jordan (Liu–Jordan model; Equations (11) and (12)) [38] and has been used widely in recent decades. The isotropic model assumes that diffuse radiation is uniformly distributed over the complete sky dome and that reflection on the ground is diffuse.
R b = cos θ cos θ z
G T = G b R b + G d ( 1 + cos β 2 ) + G ρ g ( 1 cos β 2 )
(2)
Anisotropic model
(a)
Hay and Davies model (HD model)
The Hay and Davies model (Equations (13) and (14)) accounts for both circumsolar and isotropic diffuse radiation [39,40]. Horizon brightening is not taken into account. There is an increased intensity of diffuse radiation in the area around the sun (circumsolar diffuse radiation). An anisotropy index Ai is introduced in the HD model to weight the amount of circumsolar diffuse radiation. The anisotropy index is used to quantify a portion of the diffuse radiation treated as circumsolar, with the remaining portion of diffuse radiation assumed isotropic. The circumsolar component is assumed to be from the sun’s position.
A i = G b / G 0
G T = ( G b + G d A i ) R b + G d ( 1 A i ) ( 1 + cos β 2 ) + G ρ g ( 1 cos β 2 )
  • (b)
    Hay, Davies, Klucher and Reindl model (HDKR model)
A horizon brightening diffuse term was added to the HD model by Reindl et al. in the HDKR model [39]. Horizon brightening is combined with the isotropic diffuse term and the magnitude is named by a modulating factor G b / G , as is shown in the Equation (15).
G T = ( G b + G d A i ) R b + G d ( 1 A i ) ( 1 + cos β 2 ) ( 1 + G b G sin 3 ( β 2 ) ) + G ρ g ( 1 cos β 2 )
  • (c)
    Perez I model
Compared to the other models described, the Perez model is more computationally intensive and represents a more detailed analysis of isotropic diffuse, circumsolar and horizon brightening radiation by using empirically derived coefficients [41]. Perez et al. developed the model accounting for circumsolar, horizon brightening and isotropic diffuse radiation with an empirically derived “reduced brightness coefficient” [42] in 1988. This is called the Perez I model, and is given by Equations (16)–(21). The coefficients of the Perez I model are listed in Table 3.
G T = G b R b + G d [ ( 1 F 1 ) ( 1 + cos β 2 ) + F 1 ( a c ) + F 2 sin β ] + G ρ g ( 1 cos β 2 )
a / c = Max [ 0 , cos θ ] Max [ cos 85 , cos θ z ]
ε = [ 1 + G N G d + 1.041 θ z 3 ] [ 1 + 1.041 θ z 3 ]
Δ = G d G 0
F 1 = f 11 ( ε ) + f 12 ( ε ) Δ + f 13 ( ε ) θ z
F 2 = f 21 ( ε ) + f 22 ( ε ) Δ + f 23 ( ε ) θ z
  • (d)
    Perez II model
The Perez II model has the same formulation as the Perez I model [43]. Both models differ only in the F1 and F2 coefficients. The method for calculating the detailed parameters a, c, F1 and F2 in the Perez I and Perez II models can be found in Equations (17)–(21). The coefficients of the Perez II model are shown in Table 4.

5. Diffuse Radiation and Total Tilted Radiation

5.1. Diffuse Radiation on the Horizontal Surface

Measured diffuse radiation and calculated diffuse radiation based on the DTU and RR models only using measured global radiation is shown in Figure 6. Yearly measured diffuse radiation and calculated diffuse radiation according to the DTU model was 524 and 510 kWh/m2, respectively. Yearly calculated diffuse radiation according to the RR model was 494 kWh/m2. Monthly calculated results according to the RR model were 6% lower than the measured values on average in Taars. Diffuse radiation calculated by the DTU model was closer to the measured values than the RR model. The difference between measured and simulated diffuse radiation according to the DTU model was about 3% on average. It may be concluded that the DTU model is more suitable for Danish climate conditions compared to the other universal and classic empirical model.

5.2. Total Radiation on the Tilted Surface Based on Global Radiation and Beam Radiation

Calculated total radiation on the tilted surface by use of the isotropic and anisotropic models based on measured total horizontal radiation and measured beam radiation from September 2015 to August 2016 is shown in Figure 7 together with measured values. The surface is facing south with a tilt of 50°.
The calculated monthly total tilted radiation levels according to the isotropic model were much lower than the measured values in Figure 7 compared to the anisotropic models. The measured yearly total radiation was 1170 kWh/m2. The monthly total tilted solar radiation in November 2015, December 2015 and January 2016, was around 20 kWh/m2. Contrary to the conclusions derived under Saudi Arabian and Malaysian weather conditions in past studies [11,12], in the present study, the anisotropic models were better than the isotropic model under Danish climate conditions. For the four anisotropic models, the calculated total tilted radiation according to the Perez II model and the Perez I model gave results closest to the sum of measured values, with only average differences of 1–2%, which is similar to the results reported by Andersen et al. [44].

5.3. Comparison of the Different Models

Measured data from the Taars solar heating plant were used to evaluate the models. Four statistical error parameters were introduced to evaluate the monthly results from September 2015 to August 2016 in Figure 6 and Figure 7 to determine the accuracy of the models for Danish climate conditions.
(1)
MBE, mean bias error
MBE = 1 k 1 k ( G C a l c u l a t e d i G M e a s u r e e d i )
(2)
RMSE, root mean square error
RMSE = ( 1 k 1 k ( G C a l c u l a t e d i G M e a s u r e e d i ) 2 ) 1 / 2
(3)
MAPE, mean absolute percentage error
MAPE = 1 k 1 k | G C a l c u l a t e d i G M e a s u r e e d i G M e a s u r e e d i |
(4)
RPE, relative percentage error
RPE = 1 k ( G C a l c u l a t e d i G M e a s u r e e d i ) 1 k G M e a s u r e e d i
Comparisons between measured values and calculated values of diffuse radiation on the horizontal surface and total radiation on the tilted surface are shown in Table 5 and Table 6, respectively. The lower the MBE and RMSE are, the better the agreement between the measured and calculated values. For MBE, a positive value means an overestimation of the calculated values and a negative value means an underestimation of the calculated values. A drawback of MBE is that one positive value in one calculation step may cancel a negative value in another calculation step. RMSE is always positive. MAPE is positive and a low MAPE means the model is accurate. A negative RPE means the proposed model slightly underestimates radiation. Table 5 shows that the DTU model is more accurate than the RR model for Danish conditions based on the four investigated criteria. The RMSE of the DTU model and the RR model were 2 and 3 kWh/m2, respectively. The RMSE and MAPE of the Perez models were much lower than other models, as shown in Table 6. It can be concluded that the anisotropic models (Perez II model and Perez I model) were the most accurate among the investigated empirical models and the most suitable for calculations of total tilted radiation under Danish climate conditions.

6. DNI and Total Tilted Radiation Only Based on Global Radiation

6.1. Measured and Calculated DNI Only Based on Global Radiation

Global radiation data is available from climate stations at the Danish Meteorological Institute (DMI). DNI is not measured at climate stations and only seldom in solar heating plants in Denmark. Moreover, DNI is a very important design parameter for concentrating solar collectors, such as the parabolic trough collectors in Taars. As was shown in Section 5.1 and Section 5.3, diffuse radiation calculated by the DTU model was more accurate than the RR empirical model under Danish conditions. Using Equations (1) and (2), diffuse radiation calculated by the DTU model was used to calculate DNI or beam radiation. Therefore, the DTU model was used as the optimal model in this section to predict DNI. Figure 8 shows monthly calculated DNI (DTU) and measured DNI from September 2015 to August 2016. The calculated yearly DNI (997 kWh/m2) was about 1% larger than measured yearly DNI (990 kWh/m2) for the period from September 2015 to August 2016, which was within measuring accuracy. Monthly DNI can be higher than 100 kWh/m2 in the summer season, as is shown in Figure 8.

6.2. Measured and Calculated Total Tilted Radiation Only Based on Global Radiation

As mentioned, normally global radiation from the Danish Meteorological Institute is available. Total radiation on collector surfaces is measured at most solar heating plants but with a poor accuracy in Denmark. Using the DTU model and Equations (1) and (2), calculated diffuse radiation and beam radiation could be obtained based on measured global radiation on the horizontal surface. In addition, because the isotropic model could be used easily and widely and the anisotropic models (Perez II model and Perez I model) were closest to the measured values, as described in Section 5.2 and Section 5.3, the isotropic model and the anisotropic models (Perez II and Perez I) were selected to calculate total radiation on the tilted surface based on calculated diffuse radiation and calculated beam radiation from the DTU model, which was calculated only from measured global radiation. The calculated total radiation on the tilted surface using the isotropic model (1070 kWh/m2) was 8% lower than the measured one (1170 kWh/m2) from September 2015 to August 2016 (Figure 9a). The calculated total tilted radiation by the Perez I model (1160 kWh/m2) and Perez II model (1169 kWh/m2) was less than 1% different than the yearly measured total radiation (Figure 9b,c). Both of the Perez models had the best agreement with the measurements of the investigated three empirical models.
Figure 10 shows the daily measured total tilted solar radiation as a function of the modelled total tilted solar radiation from the DTU and Perez II model. The trend in Figure 10 demonstrates good agreement between the daily measured and modelled data. It was also found that maximum daily total tilted solar radiation could be higher than 8 kWh/m2. These results are in good agreement with conclusions presented elsewhere [44]. From the above results, it was found that the DTU model together with the Perez II and I models could be used to predict total radiation on tilted surfaces based only on measured global radiation under Danish conditions very accurately. Furthermore, the proposed models could be employed to check measured total radiation on tilted flat plate collector planes in solar district heating plants in Denmark. The proposed method could also be used to derive solar radiation data for planning solar collector fields based on available horizontal global radiation measurements.

7. Discussion

In general, the anisotropic sky models (Hay and Davies, Reindl, and Perez) provide comparable estimates of the total radiation on a tilted surface and are recommended for general use [43]. The Hay and Davies and the Reindl models are computationally simple when compared to the Perez model. The isotropic sky model under-predicts total radiation on a tilted surface and is not recommended for general use. The HD and Reindl models were recommended in the mentioned references [25,26]. The HD model has also been selected to predict total tilted solar radiation in Greece [45]. The Perez models were the best models under Danish climate conditions in this study, which aligns with other past research [46,47].

8. Conclusions

Measured and calculated monthly horizontal diffuse solar radiation and total tilted solar radiation from September 2015 to August 2016 (a full year) in a demonstration solar district heating plant in Denmark were analysed in this study using an hourly time step. The DTU model, developed for the calculation of horizontal diffuse radiation in Danish climate conditions, was evaluated and validated using the measured data. Calculated monthly DNI based on the DTU model with only measured global radiation as an input was also investigated with good agreement with measurements. Furthermore, one isotropic model and four anisotropic models for general use were investigated for the calculation of total monthly radiation on the tilted surface under Danish climate conditions. From these results, the following conclusions can be drawn:
(1)
It was found that the DTU model could be used for the calculation of diffuse radiation on the horizontal surface or DNI in Denmark with better accuracy than the other classic empirical model.
(2)
Anisotropic models could be used to calculate total radiation on tilted surfaces with better accuracy than the isotropic model under Danish conditions.
(3)
The Perez models together with the DTU model could be a suitable new method to determine total radiation on tilted surfaces and double-check real-time measured solar radiation for Danish solar heating plants. The only input for this method was global radiation measurement.
(4)
Yearly global radiation and DNI was around 1000 kWh/m2 and total tilted solar radiation was around 1200 kWh/m2 in this study.
The proposed method was simple, cost-effective and gave relatively accurate measurements of total tilted radiation under Danish conditions.

Author Contributions

Zhiyong Tian set up the model, carried out the calculations and wrote the original article; Bengt Perers has contribution on the literature search and data collection. Simon Furbo has contribution to the writing and data collection; Jianhua Fan has contribution to the writing and data collection; Jie Deng provided assistance to set up the model; Janne Dragsted contributed to the modelling

Acknowledgments

Special thanks are expressed to AalborgCSP A/S (Andreas Zourellis, Jan Holst Rothman, Steffen Rovsing Møller and Per Aasted) for the information provided. The first author acknowledges financial support from the China Scholarship Council for PhD study at the Technical University of Denmark (No.201506120074). This work was also a part of an EUDP project financed by the Danish Energy Agency and the IEA-SHC Task 55 “Towards the Integration of Large SHC Systems into DHC Networks”.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

DNIDirect normal irradiance, W/m2
MBEMean bias error, kWh/m2
RMSERoot mean square error, kWh/m2
MAPEMean absolute percentage error
RPERelative percentage error
DTUTechnical University of Denmark
RR modelReduced Reindl correlation model
PTCparabolic trough collector
FPCflat plate collector
Nomenclature
RbThe ratio of beam radiation on the tilted surface to that on a horizontal surface at any time
AiAnisotropy index
kNumber of calculated values
iEvery calculated value
GMean total radiation on the horizontal surface, W/m2
GdMean diffuse radiation on the horizontal surface, W/m2
G0Mean extraterrestrial radiation on the horizontal surface, W/m2
GTMean total radiation on the tilted surface, W/m2
GbMean beam radiation on the horizontal surface, W/m2
GNMean direct normal beam radiation, W/m2
KTClearness index
a/cWeighted circumsolar solid angle
GCalculatedCalculated solar radiation, kWh/m2
GMeasuredMeasured solar radiation, kWh/m2
F1Reduced brightness coefficient (circumsolar)
F2Reduced brightness coefficient (horizon brightening)
Greek Letters
β Slope
θ z Zenith angle
θ Incident angle
α Solar altitude angle
ρ g Diffuse reflectance for the total solar radiation
ε Sky clearness parameter
Δ Sky brightness parameter

References

  1. Liu, Z.; Wu, D.; Yu, H.; Ma, W.; Jin, G. Field measurement and numerical simulation of combined solar heating operation modes for domestic buildings based on the Qinghai–Tibetan plateau case. Energy Build. 2018, 167, 312–321. [Google Scholar] [CrossRef]
  2. Yu, B.; Jiang, Q.; He, W.; Liu, S.; Zhou, F.; Ji, J.; Xu, G.; Chen, H. Performance study on a novel hybrid solar gradient utilization system for combined photocatalytic oxidation technology and photovoltaic/thermal technology. Appl. Energy 2018, 215, 699–716. [Google Scholar] [CrossRef]
  3. Liu, Z.; Li, H.; Liu, K.; Yu, H.; Cheng, K. Design of high-performance water-in-glass evacuated tube solar water heaters by a high-throughput screening based on machine learning: A combined modeling and experimental study. Sol. Energy 2017, 142, 61–67. [Google Scholar] [CrossRef]
  4. Li, H.; Liu, Z.; Liu, K.; Zhang, Z. Predictive Power of Machine Learning for Optimizing Solar Water Heater Performance: The Potential Application of High-Throughput Screening. Int. J. Photoenergy 2017, 2017, 1–10. [Google Scholar] [CrossRef]
  5. Weiss, W.; Monika Spörk-Dür, F.M. Solar Heat Worldwide-Global Market Development and Trends in 2016-Detailed Market Figures 2015 (2017 Version). 2017. Available online: http://www.iea-shc.org/solar-heat-worldwide (accessed on 7 May 2018).
  6. Shukla, K.N.; Rangnekar, S.; Sudhakar, K. Comparative study of isotropic and anisotropic sky models to estimate solar radiation incident on tilted surface: A case study for Bhopal, India. Energy Rep. 2015, 1, 96–103. [Google Scholar] [CrossRef]
  7. Demain, C.; Journée, M.; Bertrand, C. Evaluation of different models to estimate the global solar radiation on inclined surfaces. Renew. Energy 2013, 50, 710–721. [Google Scholar] [CrossRef]
  8. Khorasanizadeh, H.; Mohammadi, K.; Mostafaeipour, A. Establishing a diffuse solar radiation model for determining the optimum tilt angle of solar surfaces in Tabass, Iran. Energy Convers. Manag. 2014, 78, 805–814. [Google Scholar] [CrossRef]
  9. Marques Filho, E.P.; Oliveira, A.P.; Vita, W.A.; Mesquita, F.L.L.; Codato, G.; Escobedo, J.F.; Cassol, M.; França, J.R.A. Global, diffuse and direct solar radiation at the surface in the city of Rio de Janeiro: Observational characterization and empirical modeling. Renew. Energy 2016, 91, 64–74. [Google Scholar] [CrossRef]
  10. El Mghouchi, Y.; El Bouardi, A.; Choulli, Z.; Ajzoul, T. Models for obtaining the daily direct, diffuse and global solar radiations. Renew. Sustain. Energy Rev. 2016, 56, 87–99. [Google Scholar] [CrossRef]
  11. Jakhrani, A.Q.; Othman, A.; Rigit, A.R.H.; Samo, S.R.; Ahmed, S. Estimation of Incident Solar Radiation on Tilted Surface by Different Empirical Models. Int. J. Sci. Res. Publ. 2012, 2, 15–20. [Google Scholar]
  12. El-Sebaii, A.A.; Al-Hazmi, F.S.; Al-Ghamdi, A.A.; Yaghmour, S.J. Global, direct and diffuse solar radiation on horizontal and tilted surfaces in Jeddah, Saudi Arabia. Appl. Energy 2010, 87, 568–576. [Google Scholar] [CrossRef]
  13. Gopinathan, K.K. Solar radiation on variously oriented sloping surfaces. Sol. Energy 1991, 47, 173–179. [Google Scholar] [CrossRef]
  14. Li, H.; Ma, W.; Lian, Y.; Wang, X. Estimating daily global solar radiation by day of year in China. Appl. Energy 2010, 87, 3011–3017. [Google Scholar] [CrossRef]
  15. AlYahya, S.; Irfan, M.A. Analysis from the new solar radiation Atlas for Saudi Arabia. Sol. Energy 2016, 130, 116–127. [Google Scholar] [CrossRef]
  16. Bird, R.E.; Riordan, C. Simple Solar Spectral Model for Direct and Diffuse Irradiance on Horizontal and Tilted Planes at the Earth’s Surface for Cloudless Atmospheres. J. Clim. Appl. Meteorol. 1986, 25, 87–97. [Google Scholar] [CrossRef]
  17. Lou, S.; Li, D.H.W.; Lam, J.C.; Chan, W.W.H. Prediction of diffuse solar irradiance using machine learning and multivariable regression. Appl. Energy 2016, 181, 367–374. [Google Scholar] [CrossRef]
  18. Kashyap, Y.; Bansal, A.; Sao, A.K. Solar radiation forecasting with multiple parameters neural networks. Renew. Sustain. Energy Rev. 2015, 49, 825–835. [Google Scholar] [CrossRef]
  19. Despotovic, M.; Nedic, V.; Despotovic, D.; Cvetanovic, S. Evaluation of empirical models for predicting monthly mean horizontal diffuse solar radiation. Renew. Sustain. Energy Rev. 2016, 56, 246–260. [Google Scholar] [CrossRef]
  20. Ineichen, P. Global Irradiance on Tilted and Oriented Planes: Model Validations; Technical Report; University Geneva: Geneva, Switzerland, October 2011. [Google Scholar]
  21. 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]
  22. Polo, J.; Garcia-Bouhaben, S.; Alonso-García, M.C. A comparative study of the impact of horizontal-to-tilted solar irradiance conversion in modelling small PV array performance. J. Renew. Sustain. Energy 2016, 8. [Google Scholar] [CrossRef]
  23. Boussaada, Z.; Curea, O.; Remaci, A.; Camblong, H.; Mrabet Bellaaj, N. A Nonlinear Autoregressive Exogenous (NARX) Neural Network Model for the Prediction of the Daily Direct Solar Radiation. Energies 2018, 11, 620. [Google Scholar] [CrossRef]
  24. Frydrychowicz-Jastrzębska, G.; Bugała, A. Modeling the Distribution of Solar Radiation on a Two-Axis Tracking Plane for Photovoltaic Conversion. Energies 2015, 8, 1025–1041. [Google Scholar] [CrossRef]
  25. Mubarak, R.; Hofmann, M.; Riechelmann, S.; Seckmeyer, G. Comparison of modelled and measured tilted solar irradiance for photovoltaic applications. Energies 2017, 10, 1688. [Google Scholar] [CrossRef]
  26. Lee, H.J.; Kim, S.Y.; Yun, C.Y. Comparison of solar radiation models to estimate direct normal irradiance for Korea. Energies 2017, 10, 594. [Google Scholar] [CrossRef]
  27. Aalborg CSP A/S. Available online: http://www.aalborgcsp.com/ (accessed on 7 May 2018).
  28. TRNSYS Website TRNSYS 17—A TRaNsient System Simulation program-Standard Component Library Overview and Mathematical Reference. Available online: http://www.trnsys.com/assets/docs/03-ComponentLibraryOverview.pdf (accessed on 28 March 2017).
  29. Dragsted, J.; Furbo, S. Solar Radiation and Thermal Performance of Solar Collectors for Denmark; DTU Civil Engineering Reports, No. R-275; Department of Civil Engineering, Technical University of Denmark: Copenhagen, Denmark, 2012. [Google Scholar]
  30. Reindl, D.T.; Beckman, W.A.; Duffie, J.A. Diffuse fraction correlations. Sol. Energy 1990, 45, 1–7. [Google Scholar] [CrossRef]
  31. Tian, Z.; Perers, B.; Furbo, S.; Fan, J. Analysis and validation of a quasi-dynamic model for a solar collector field with flat plate collectors and parabolic trough collectors in series for district heating. Energy 2018, 142, 130–138. [Google Scholar] [CrossRef]
  32. Tian, Z.; Perers, B.; Furbo, S.; Fan, J. Thermo-economic optimization of a hybrid solar district heating plant with flat plate collectors and parabolic trough collectors in series. Energy Convers. Manag. 2018, 165, 92–101. [Google Scholar] [CrossRef]
  33. Tian, Z.; Perers, B.; Furbo, S.; Fan, J. Annual measured and simulated thermal performance analysis of a hybrid solar district heating plant with flat plate collectors and parabolic trough collectors in series. Appl. Energy 2017, 205, 417–427. [Google Scholar] [CrossRef]
  34. DMI Danish Meteorological Institute. 2013. Available online: http://www.dmi.dk/fileadmin/user_upload/Rapporter/TR/2013/TR13-19.pdf (accessed on 7 May 2018).
  35. Tian, Z.; Perers, B.; Furbo, S.; Fan, J. Analysis of measured and modeled solar radiation at the taars solar heating plant in denmark. In Proceeding of the 11th Ises Eurosun 2016 confrenence, Palma de Mallorca, Spain, 11–14 October 2016. [Google Scholar]
  36. SMP11, K. Kipp&Zonen SMP11. 2017. Available online: http://www.kippzonen.com/Product/202/SMP11Pyranometer#.V35QMvl95mM (accessed on 7 May 2018).
  37. PMO6-CC PMO6-CC Pyrheliometer. 2016. Available online: https://www.pmodwrc.ch/pmod.php?topic=pmo6 (accessed on 7 May 2018).
  38. 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]
  39. Duffie, J.A.; Beckman, W.A. Solar Engineering of Thermal Process; John Wiley & Sons Inc.: Hoboken, NJ, USA, 2013. [Google Scholar]
  40. Hay, J.E. Calculations of the solar radiation incident on an inclined surface. In Proceedings of the First Canadian Solar Radiation Data Workshop, Toronto, ON, Canada, 17–19 April 1978; pp. 59–72. [Google Scholar]
  41. Loutzenhiser, P.G.; Manz, H.; Felsmann, C.; Strachan, P.A.; Frank, T.; Maxwell, G.M. Empirical validation of models to compute solar irradiance on inclined surfaces for building energy simulation. Sol. Energy 2007, 81, 254–267. [Google Scholar] [CrossRef] [Green Version]
  42. Perez, R.; Stewart, R.; Seals, R.; Guertin, T. The Development and Verification of the Perez Diffuse Radiation Model; Atmospheric SciencesResearch Center: Albany, NY, USA, 1988. [Google Scholar]
  43. Thermal Energy System Specialists (TRNSYS). TRNSYS 17—Volume 7 Programmer’s Guide; TRNSYS: Madison, WI, USA, 2014. [Google Scholar]
  44. Andersen, E.; Hans, L.; Furbo, S. The influence of the solar radiation model on the calculated solar radiation from a horizontal surface to a tilted surface. In Proceedings of the EuroSun Congress, Freiburg, Germany, 20–23 June 2014. [Google Scholar]
  45. Raptis, P.I.; Kazadzis, S.; Psiloglou, B.; Kouremeti, N.; Kosmopoulos, P.; Kazantzidis, A. Measurements and model simulations of solar radiation at tilted planes, towards the maximization of energy capture. Energy 2017, 130, 570–580. [Google Scholar] [CrossRef]
  46. De Simón-Martín, M.; Alonso-Tristán, C.; Díez-Mediavilla, M. Diffuse solar irradiance estimation on building’s façades: Review, classification and benchmarking of 30 models under all sky conditions. Renew. Sustain. Energy Rev. 2017, 77, 783–802. [Google Scholar] [CrossRef]
  47. Yang, D. Solar radiation on inclined surfaces: Corrections and benchmarks. Sol. Energy 2016, 136, 288–302. [Google Scholar] [CrossRef]
Figure 1. Location of Taars in Denmark [34].
Figure 1. Location of Taars in Denmark [34].
Energies 11 01315 g001
Figure 2. Location of the weather station and pyranometers (PTC: parabolic trough collector, FPC: flat plate collector) [27].
Figure 2. Location of the weather station and pyranometers (PTC: parabolic trough collector, FPC: flat plate collector) [27].
Energies 11 01315 g002
Figure 3. The pyranometers in the middle of the flat plate collector field.
Figure 3. The pyranometers in the middle of the flat plate collector field.
Energies 11 01315 g003
Figure 4. Weather station used and the pyrheliometer of the Taars solar heating plant. (a) Weather station; (b) the weather sensors including pyrheliometer; and (c) used pyrhelionmeter and sunscanner.
Figure 4. Weather station used and the pyrheliometer of the Taars solar heating plant. (a) Weather station; (b) the weather sensors including pyrheliometer; and (c) used pyrhelionmeter and sunscanner.
Energies 11 01315 g004
Figure 5. Schematic illustration of the methodology. (a) Validation cycle; (b) Calculation cycle.
Figure 5. Schematic illustration of the methodology. (a) Validation cycle; (b) Calculation cycle.
Energies 11 01315 g005
Figure 6. Calculated and measured diffuse radiation on the horizontal surface (September 2015–August 2016).
Figure 6. Calculated and measured diffuse radiation on the horizontal surface (September 2015–August 2016).
Energies 11 01315 g006
Figure 7. Calculated and measured total radiation on the 50° tilted south facing surface (September 2015–August 2016).
Figure 7. Calculated and measured total radiation on the 50° tilted south facing surface (September 2015–August 2016).
Energies 11 01315 g007
Figure 8. Calculated DNI (DTU) and measured DNI in the Taars solar heating plant.
Figure 8. Calculated DNI (DTU) and measured DNI in the Taars solar heating plant.
Energies 11 01315 g008
Figure 9. Measured monthly tilted total radiation and calculated tilted total radiation based on calculated diffuse radiation and beam radiation (September 2015–August 2016: (a) DTU and Isotropic model; (b) DTU and Perez I model; (c) DTU and Perez II model).
Figure 9. Measured monthly tilted total radiation and calculated tilted total radiation based on calculated diffuse radiation and beam radiation (September 2015–August 2016: (a) DTU and Isotropic model; (b) DTU and Perez I model; (c) DTU and Perez II model).
Energies 11 01315 g009aEnergies 11 01315 g009b
Figure 10. Daily measured total tilted solar radiation as a function of modelled total tilted solar radiation from the DTU and Perez II model.
Figure 10. Daily measured total tilted solar radiation as a function of modelled total tilted solar radiation from the DTU and Perez II model.
Energies 11 01315 g010
Table 1. Specifications of Kipp & Zonen SMP11 pyranometer.
Table 1. Specifications of Kipp & Zonen SMP11 pyranometer.
ParameterValues
Spectral range (50% points)285 to 2800 nm
Response time (63%)<0.7 s
Response time (95%)<2 s
Zero offset A<7 W/m2
Zero offset B<2 W/m2
Directional response (up to 80° with 1000 W/m² beam)<10 W/m2
Temperature dependence of sensitivity (−20 °C to +50 °C)<1%
Analogue output (−V version)0 to 1 V
Analogue output (−A version)4 to 20 mA
Table 2. Specifications of PMO6-CC pyrheliometer.
Table 2. Specifications of PMO6-CC pyrheliometer.
ParameterValues
Dimension80 × 80 × 230 mm
Mass2.15 kg
Field of view (full angle)
Slope angle
Rangeup to 1400 W/m2 (or custom design available)
Traceability to WRR<0.1%
Operating temperature range−25 °C to +50 °C
Table 3. The coefficients of the Perez I model
Table 3. The coefficients of the Perez I model
ε BinUpper Limit for ε Cases (%)f11f12f13f21f22f23
11.06513.6−0.1961.084−0.006−0.1140.18−0.019
21.235.60.2360.519−0.18−0.0110.2−0.038
31.57.520.4540.321−0.2550.072−0.098−0.046
41.958.870.866−0.381−0.3750.203−0.403−0.049
52.813.171.026−0.711−0.4260.273−0.602−0.061
64.521.450.978−0.986−0.350.28−0.915−0.024
76.216.060.748−0.913−0.2360.173−1.0450.065
8-13.730.318−0.7570.1030.062−1.6980.236
Table 4. The coefficients of the Perez II model.
Table 4. The coefficients of the Perez II model.
ε Binf11f12f13f21f22f23
1−0.008310.58773−0.06206−0.059600.07212−0.02202
20.129990.68260−0.15138−0.018930.06597−0.02887
30.329700.48687−0.22110−0.055414−0.06396−0.02605
40.568210.18745−0.29513−0.10886−0.15192−0.01398
50.87303−0.39204−0.361620.22556−0.46204−0.00124
61.13261−1.23673−0.411850.28778−0.823040.05587
71.06016−1.59999−0.358920.26421−1.127230.13107
80.67775−0.32726−0.250430.26421−1.376500.25062
Table 5. Measured and calculated MBE (kWh/m2), RMSE (kWh/m2), MAPE (%) and RPE (%) for diffuse horizontal radiation.
Table 5. Measured and calculated MBE (kWh/m2), RMSE (kWh/m2), MAPE (%) and RPE (%) for diffuse horizontal radiation.
ItemsDTU ModelRR Model
MBE−1.3−2.5
RMSE2.03.0
MAPE3.5%8.1%
RPE−2.9%−5.7%
Table 6. Measured and calculated MBE (kWh/m2), RMSE (kWh/m2), MAPE (%) and RPE (%) for monthly total tilted radiation.
Table 6. Measured and calculated MBE (kWh/m2), RMSE (kWh/m2), MAPE (%) and RPE (%) for monthly total tilted radiation.
ItemsPerez II ModelPerez I ModelHDKR ModelHD ModelIsotropic Model
MBE−2.4−3.4−8.4−10.0−18.6
RMSE2.02.64.95.810.0
MAPE2.1%2.8%5.7%5.9%12.0%
RPE−1.2%−1.8%−4.3%−5.2%−9.7%

Share and Cite

MDPI and ACS Style

Tian, Z.; Perers, B.; Furbo, S.; Fan, J.; Deng, J.; Dragsted, J. A Comprehensive Approach for Modelling Horizontal Diffuse Radiation, Direct Normal Irradiance and Total Tilted Solar Radiation Based on Global Radiation under Danish Climate Conditions. Energies 2018, 11, 1315. https://doi.org/10.3390/en11051315

AMA Style

Tian Z, Perers B, Furbo S, Fan J, Deng J, Dragsted J. A Comprehensive Approach for Modelling Horizontal Diffuse Radiation, Direct Normal Irradiance and Total Tilted Solar Radiation Based on Global Radiation under Danish Climate Conditions. Energies. 2018; 11(5):1315. https://doi.org/10.3390/en11051315

Chicago/Turabian Style

Tian, Zhiyong, Bengt Perers, Simon Furbo, Jianhua Fan, Jie Deng, and Janne Dragsted. 2018. "A Comprehensive Approach for Modelling Horizontal Diffuse Radiation, Direct Normal Irradiance and Total Tilted Solar Radiation Based on Global Radiation under Danish Climate Conditions" Energies 11, no. 5: 1315. https://doi.org/10.3390/en11051315

APA Style

Tian, Z., Perers, B., Furbo, S., Fan, J., Deng, J., & Dragsted, J. (2018). A Comprehensive Approach for Modelling Horizontal Diffuse Radiation, Direct Normal Irradiance and Total Tilted Solar Radiation Based on Global Radiation under Danish Climate Conditions. Energies, 11(5), 1315. https://doi.org/10.3390/en11051315

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