For many locations in the world the average daily irradiation on a horizontal surface can be looked up. For a more accurate indication of how much energy is falling on a solar panel titled at an angle from the horizontal, a series of calculations must be done. The following sections contain two methods: section 4.1 contains a simplified calculation to find the total daily irradiation falling on a tilted surface that faces towards the equator; section 4.2 describes a more complex calculation that calculates the irradiation on a titled plane hour by hour through the day. Both methods follow the same basic path.

The measured irradiation data that can be referenced is normally in the form of the monthly average daily (MAD) irradiation on a horizontal plane; that is the total irradiation measured in each month divided by the number of days in the month. This figure is made up of the beam and diffused components described in section 2.5 and it is assumed that no reflected radiation is incident on a horizontal surface. To find the MAD irradiation on a tilted surface at the same location the measured data is broken down into its beam and diffused component parts, these are then adjusted for the chosen panel angle and a reflected component added. Because the beam irradiation is a directional and the diffused radiation is not the adjustment for these two values is not the same. Finally the three components are added together to give the total MAD irradiation on the titled panel.

In general:

(W/m^{2})(4.1)

Where:

I(β) = terrestrial insolation intensity on a plane tilted at β degrees to the horizontal surface (W/m^{2});

I_{Bh} = terrestrial beam insolation intensity on a horizontal plane (W/m^{2});

θ_{i} = angle of incidence of the radiation on the tilted plane;

I_{D} = diffused irradiation component on a horizontal plane (W/m^{2});

R_{D} = a convention factor to adjust I_{D} to a value for a titled plane;

I_{Th} = the total terrestrial insolation intensity on a horizontal plane (W/m^{2});

R_{A} = a convention factor used to find the reflected (or albedo) component for a titled plane.

Note that if θ_{i }> 90° the titled plane will be shading itself. Note also that θ_{i }changes throughout the day as do R_{D} and R_{A}. I(β) can then be used to find the total irradiation incident on a titled plane (H(β)) during a whole day by integrating equation 3.14 with respect to time using the sunrise and sunset angles as limits. Thus:

(J/m^{2})(4.2)

Equation 3.15 can be used to take account of a titled surface shading itself by using values of ω_{0} described in section 3.2.

The equations can be used for monthly average values or values for single days.

### 4.1. Simplified Method For The Surfaces Tilted Towards The Equator

** ***Work In Radians, Be Careful To Work Only In Wh/m2, kWh/m2, J/m2**

#### 4.1.1. Calculating The Energy Falling On An Extraterrestrial Horizontal Surface

Calculate the insolation intensity on a plane perpendicular to sun rays (I_{0}) on the edge of the atmosphere by correcting I_{SC} for Earth’s elliptical orbit.

(W/m^{2})(4.3)

Calculate the extraterrestrial insolation on a plane horizontal to the Earth’s surface (I_{0h}) at the site’s latitude.

(W/m^{2})(4.4)

Where:

And:

(rads)(4.5)

Equations 4.3 and 4.5 are now combined to give:

(W/m^{2})(4.6)

To find the extraterrestrial irradiation energy falling on a plane horizontal to the Earth’s surface throughout a whole day, integrate equation 4.6 with respect to time between sunrise (ω = -ω_{S}) and sunset (ω = ω_{S}). The resulting equation is:

(J/m^{2})(4.7)

Note that 86 400 = number of seconds in 24 hours and equation 4.7 gives H_{0h} in J/m^{2}. Since there are 3600 seconds in an hour and 86 000/3 600 = 24, the following equation gives H_{0h} in Wh/m^{2}.

(Wh/m^{2})(4.8)

Where:

(4.9)

If n is set to the 15^{th} of each month (i.e. n = 15 for January and n = 46 for February etc), equation 4.8 can be used to calculate monthly average daily H_{0h} values. When MAD irradiation values are used the symbol replaces H_{0h} to indicate average values.

#### 4.1.2. The Clearness Index

Next the clearness index (K_{T}) is defined as the ratio of the total irradiation reaching a horizontal plane at the location on the Earth’s surface and the extraterrestrial irradiation on a horizontal plane above the location. Thus K_{T} is an indication of how much of the Sun’s radiation is lost to scattering and absorption in the atmosphere. If MAD values are used:

(4.10)

Where:

= the monthly average daily irradiation on a horizontal plane at the Earth’s surface. This is the value that can be referenced from weather station data.

= the monthly average daily value of extraterrestrial radiation energy falling on a horizontal plane. Calculated from equation 4.7 or 4.8.

K_{T }values calculated in this way are valid in the range 0.3 ≤ K_{T} ≤ 0.8 between v = ±55°. K_{T} can be used to find the ratio of diffused irradiation to total irradiation for a horizontal surface:

Therefore:

(J/m^{2} or Wh/m^{2})(4.11)

Where is the monthly average daily value for the diffused irradiation incident on a horizontal surface.

#### 4.1.3. Calculation Of The Beam Irradiation On A Tilted Surface

It is assumed that no irradiation reflected from the ground makes it onto a horizontal surface. Therefore:

(4.12)

Where is the monthly average daily value for the beam irradiation on a horizontal plane. Since is measured and has been calculated from equation 4.10, the beam irradiation component incident on a horizontal surface on the average day for each month is given by:

(4.13)

For a titled plane at angle β to horizontal we can define the ratio R_{B}:

(4.14)

It is assumed that the ratio of extraterrestrial irradiation on a surface titled at angle β to the extraterrestrial irradiation on a horizontal plane also equals R_{B} so that:

(4.15)

In equations 4.6 and 4.7 ω_{S} is used for a horizontal plane, however for titled plane ω_{0} must be used as described in section 3.2, also for a tilted plane φ is replaced by (φ – β) (from equation 3.13). So that:

(4.16)

where ω_{0} is either ω_{S} or ω_{S}‘, depending on which is the smaller value and:

(4.17)

Thus, from figure 4.6 and 4.7:

(4.18)

(J/m^{2} or Wh/m^{2})(4.19)

Values of ωS and ω0 should be kept in radians in equation 4.19.

#### 4.1.4. Calculation Of The Diffused Component On A Tilted Surface

The diffused irradiation on a horizontal surface was calculated with equation 4.11, this value is in turn used to calculate the diffused irradiation on a surface tilted from the horizontal. To this end a convention factor is defined as the ratio:

(4.20)

Thus:

(J/m^{2} or Wh/m^{2})(4.21)

#### 4.1.5. Calculation Of The Reflected Component On A Tilted Surface

It is assumed that no reflected (albedo) radiation is incident on a horizontal plane so that to calculate the monthly average daily reflected radiation on a tilted surface () is found by multiplying by a conversion factor:

(4.22)

(J/m^{2} or Wh/m^{2})(4.33)

where ρ is the reflectivity of the surrounding ground. Values of ρ are given in table 4.1.

Ground Cover | Reflectivity | Ground Cover | Reflectivity |
---|---|---|---|

Dry bare ground | 0.2 | Pale soil | 0.3 |

Dry grassland | 0.3 | Dark soil | 0.1 |

Desert sand | 0.4 | Water | 0.1 |

Snow | 0.5-0.8 | Vegetation | 0.2 |

#### 4.1.6. Calculation Of The Total Irradiation On A Tilted Surface

The beam, diffused and reflected components are now added together to give the monthly average daily total irradiation on the titled surface. Thus:

(J/m^{2} or Wh/m^{2})(4.34)

Often written as:

(J/m^{2} or Wh/m^{2})(4.35)

Note that if β = 0 (surface horizontal) then cos β = 1 and reflected irradiation does not contribute to the total irradiation.

The method detailed in this section can be used for an average day for each month and for any value of β. A spreadsheet can be used to calculate values for all possible angles and therefore the optimum angle can be found for each month.

### 4.2. Hourly Method For Surfaces Tilted In Any Direction

#### 4.2.1. Calculation Of Hourly Irradiation Values

Firstly, equation 4.7 is used to calculate for each month. Then is calculated from equation 4.10 (0.3≤K_{T}≤0.8 still applies). For this method equation 4.11 is replaced with the following relationships:

when ω_{S}<81.4° (4.36)

when ω_{S}>81.4° (4.37)

Where ω_{S} is calculated from equation 3.13.

The average daily irradiation values for each month are broken down into hourly values using the ratio r_{t}:

Note that to calculate the hour angle (ω) the hours should be centred on the half hour (e.g. 11:00-12:00).

(4.38)

where:

Where ω is the hour angle in the middle of each hour and is quoted in radians. ω_{ S} is calculated from equation 3.13 and is also in radians. ω is negative in the morning, zero at solar noon and positive in the afternoon.

The total irradiation for each hour ( ) of the ‘average’ day of each month is then calculated:

(4.39)

The diffused component of the total irradiation for each hour is now calculated using the ratio:

(4.40)

So that the hourly diffused irradiation for the ‘average’ day of each month:

(4.41)

The hourly beam component ( ) is calculated from the relationship stated in equations 4.12 and 4.13:

(4.42)

#### 4.2.2. Calculation Of The Total Hourly Irradiation On A Titled Surface

The hourly irradiation on a tilted surface ( ) is calculated using an equation similar to equation 4.35:

In this equation:

Where θ_{i} is the angle of incidence of the Sun’s rays on the tilted plane (figure 3.2) calculated from equation 3.11, and θ_{Z} is the solar zenith angle (figure 1.7)calculated from equation 3.12.

Using this method, values of ω_{S}’ are not used (section 3.2) and since hourly figures of θ_{i} and θ_{Z} are calculated, self-shading of the tilted plane is apparent when θ_{i} is greater than 90° and becomes negative – and therefore is taken as zero. Note the is an hourly figure.

The daily totals for the average day of each month are then found by adding the hourly figures together.