Published on *EME 810: Solar Resource Assessment and Economics* (https://www.e-education.psu.edu/eme810)

- J.R. Brownson,
*Solar Energy Conversion Systems*(SECS),**Chapter 8**: Measure & Estimation of the Solar Resource (Isotropic and Anisotropic Sky Models.) - [repeated] C. A. Gueymard (2008) From Global Horizontal to Global Tilted Irradiance: How accurate are solar energy engineering predictions in practice? [1] Solar 2008 Conference, San Diego, CA, American Solar Energy Society

Please make sure you read all of Ch 8 in SECS for this lesson, still related to "Empirical Correlation for Components," but paying attention to the isotropic and anisotropic sky models and this page content. This is the third page for which we have included review of the Gueymard paper, so you should be familiar with the findings by now. The Perez et al. paper will be useful to you in Learning Activity 4.2 and for your Lesson Quiz.

Earlier, we discussed the different components of light, beam and diffuse on a horizontal surface. Now, we will discuss how these components can be estimated for tilted surfaces through isotropic or anisotropic diffuse sky/ground models of light source components.

We shall see that we do not need to measure every component of light (scattered and unscattered) to make estimations on the contributions of each component to the total irradiation incident on the aperture of interest. We can rely somewhat on decades of historical observation and empirical correlation by solar scientists and engineers for hourly, daily, and monthly average day data.

The main tools we need are the equations for **hourly** and **daily extraterrestrial irradiance** (Air Mass Zero, or **AM0**) and the **integrated energy density** (irradiation: $MJ/m^2$) gathered from a **horizontally mounted pyranometer**, which you learned of in the last section. We shall also find that we can infer more than just the components of light from the ratios of measured irradiation to AM0 calculated irradiation--we can **describe the fractions of days in a given month where lighting conditions will be clear or overcast/cloudy.**

Following our step to break apart the beam horizontal component from the diffuse horizontal components, we then estimate the components on a tilted surface.

For a tilted plane of array,

Total Radiation = beam + diffuse, *sky* + diffuse, *ground*

$${G}_{t}={G}_{b,t}+{G}_{b,d}+{G}_{g,t}$$

A simple calculation of the beam component can be achieved using

$${G}_{b,t}={G}_{b,n}.cos\theta $$

Radiation on a sloped surface can be calculated for the **beam component** of irradiation by the geometric scaling factor of

$${R}_{b}=\frac{cos(\theta )}{cos({\theta}_{z})}$$

In order to estimate the** diffuse component**, we use alternate models that become increasingly better fits with the empirical data. We can integrate any of these equations over an hour or a day (irradiation). I prefer to offer the irradiance version as a bit easier to read. Note: * all* of these estimation models use irradiation values that were measured from a

The isotropic sky model was developed in the 1960s to estimate the diffuse sky on a tilted surface, complemented by an estimate for diffuse light from the ground. This model assumes that the sky is uniform in composition across the sky dome.

The following expression gives the total solar irradiance incident on a tilted surface as

$$G={\text{beam+diffuse}}_{\text{sky}}{\text{+diffuse}}_{\text{ground}}$$ $${G}_{total}={G}_{b}{R}_{b}+{G}_{d}({F}_{surface-sky})+G\rho ({F}_{surface-ground})$$

where,

$${F}_{surface-sky}=\frac{1+cos\beta}{2}$$ $${F}_{surface-ground}=\frac{1-cos\beta}{2}$$ $${G}_{d,t}={G}_{d}.\frac{1+cos\beta}{2}$$ $${G}_{g,t}={\rho}_{g}({G}_{b}+{G}_{d}).\frac{1-cos\beta}{2}$$

The fraction proportional to the collector tilt is called the diffuse sky irradiance tilt factor for an isotropic sky model, and the reflectance of the ground is called the albedo (a fraction between 0 and 1), and is multiplied by the GHI and the diffuse ground irradiance tilt factor for an isotropic sky model.

Note: "Surface": the aperture. $\rho $ : is the collective reflectivity of the ground (the **albedo**). $\rho $ reduces the irradiance G by a value between 0--1. On an inclined surface, *G _{d,ground}* increases, relative to a horizontal collector.

This model incorporates isotropic diffuse, circumsolar radiation and horizontal brightening. It also employs an anisotropic index A defined mathematically as

$$A=\frac{{G}_{N,I}}{G}$$

The total irradiance on a tilted surface is then calculated by using

$${G}_{T}=({G}_{b}+{G}_{d}A){R}_{B}+{G}_{d}(1-A)\frac{1+\mathrm{cos}\beta}{2}\times \left[1+\sqrt{\frac{{G}_{b}}{G}}{\mathrm{sin}}^{3}\frac{\beta}{2}\right]+G\rho \frac{1-\mathrm{cos}\beta}{2}$$

Go to Kalogirou (Solar Energy Engineering) Ch 2 (pdf from Library) this will be labeled the "Reindl model"

This is an **anisotropic diffuse sky model **that takes into consideration the real observations of subcomponents of diffuse light. The Perez model adds the **circumsolar diffuse component** and the **horizon diffuse component** to the diffuse$_{sky}$ component of the isotropic model. Notice how the beam component is not mentioned here--it doesn't change.

*Sidenote:* Richard Perez is a Senior Research Associate in the Atmospheric Sciences Research Center in SUNY Albany. He has a great website [2].

$${G}_{d}={G}_{iso}+{G}_{cir}+{G}_{hor}+diffus{e}_{ground}$$ $${G}_{d}=\left[{G}_{d}(1-{F}_{1})\frac{1+cos\beta}{2}\right]+\left[{G}_{d}({F}_{1})\frac{cos\theta}{cos{\theta}_{z}}\right]+\left[Gd({F}_{2})sin\beta \right]+\left[{G}_{d}\rho \frac{1-cos\beta}{2}\right]$$ The shape factors (F) in this model can be reviewed in the original article by Perez et. al (1990). However, we can inspect the equations and observe in the equation that $F_{surface-sky}$ is *reduced *by a proportion of $F_1$ (circumsolar radiance), and $F_2$ can either increase or decrease the contribution of horizon radiance.

- You can open up SAM [3] at this point and click the button "Create a New File."
- Click on "Photovoltaics (Detailed)" on the left of the pop up window ("Choose a performance model").
- Click on "Residential (Distributed)" on the right of the pop up window ("then choose from the available financial models").
- Hit "OK."
- You will now have a default residential PV project, based in Phoenix, AZ. There will likely be a SunPower (California company) PV module, and an SMA America dc-ac power inverter listed on the left column of shortcut tabs. They're not really "tabs" like the tab at the top of the window, but I'm going to use that term for convenience in the class.
- Click on the shortcut tab called "Location and Resource" (will list "Location," latitude ($\phi$), longitude ($\lambda$), and elevation above sea level in the shortcut.
- At the bottom of the screen on the right is a large "+" sign labeled "PV Albedo and Radiation". Please expand the section. I want to draw your attention to the
*little section*called "Diffuse Sky Model (Advanced)." Yep, it's so advanced that you don't even realize this is where a powerful data transformation sits!

You can see that one may select "Irradiance Components used for Calculation": this is specifying the type of horizontal irradiation components that you will use in your tilted model. In a data set called the Typical Meteorological Year, the data for the beam is often not actually a measured value.

You can also see that one may select three diffuse sky/ground transposition models to transform the "Irradiance Components" (horizontal) to tilted values. The default is the Perez model that we describe below (and in your supplemental reading). The isotropic model is not used in practice, but it contains the basis for the other anisotropic models of Hay-Davies-Klutcher-Reindl (HDKR) and Perez et al. 1990.

- Tiny box hidden within a plus sign, whole lotta power. Keep that in mind.
- Also, when in doubt: Read the Fine Manual. SAM has a detailed Help manual with links to the literature for every aspect of the program.

- Perez, Ineichen, and Seals (1990) Modeling Daylight Availability and Irradiance Components from Direct and Global Irradiance [4]. Solar Energy J. 44(5), 271-289.

Liu, B.Y.H., Jordan, R.C., 1960. The interrelationship and characteristic distribution of direct, diffuse, and total solar radiation. Solar Energy 4(3),1-19

Perez, R., Ineichen, P., Seals, R., Michalsky, J., Stewart, R., 1990. Modeling daylight and Irradiance components from direct and global irradiance. Solar Energy 44(5), 271-289

**Links**

[1] https://www.researchgate.net/publication/236314649_From_global_horizontal_to_global_tilted_irradiance_How_accurate_are_solar_energy_engineering_predictions_in_practice

[2] http://www.asrc.cestm.albany.edu/perez/

[3] https://sam.nrel.gov/

[4] http://www.sciencedirect.com/science/article/pii/0038092X9090055H