Shortwave solar

radiation

1
Calculating equation coefficients

Construction Conservation Equation

Surface Conservation Equation

Fluid Conservation Equation

needs flow estimation

needs radiation and

convection estimation

2
The Sun

 Core temperature 8x106 to 40x106 K.

 Effective black body temperature of 6000 K.
 Solar constant: extraterrestrial flux from the
sun received on a unit area perpendicular to
the direction of propagation – mean
Sun/Earth distance value is 1353 W/m2.
 Actual extraterrestial radiation varies with
time of year as earth-sun distance varies.

3
Energy from the sun
Incoming
Longwave solar energy Reflected shortwave
radiation to space 175 . 1015 W radiation
122.5 . 1015 W 52.5 . 1015 W

Atmospheric boundary

Convection currents (wind and

ocean waves)
368 . 1012 W
Tidal energy
3 . 1012 W
Evaporation of water, heating
of water & ice
Geothermal 40 . 1015 W
energy
32 . 1012 W Photosynthesis on land and sea
98 . 1012 W

Earth’s surface Direct conversion to heat

82 . 1015 W
Formation of fossil fuels
13 . 106 W
4
Atmospheric interactions

 The greater the distance

that the radiation passes
through the atmosphere, the
greater is the frequency
dependent scattering.
Spectra at ground level are
often referred to particular
‘air masses’.
 Air Mass 1 is the thickness
of the atmosphere vertically
above sea level.
 Air Mass 2 is double this
2
thickness (equivalent to 1
Atmosphere
direct solar radiation at an
altitude of 30 degrees). 30°

5
Direct and diffuse radiation

 Solar radiation reaches the Earth

directly from the Sun) and On clear days
diffusely after scattering in the around 90% of
atmosphere and reflected from the total solar
surrounding objects. radiation is direct.

 Only direct radiation can be

focussed.

 The total radiation reaching a On heavily

surface is the summation of the overcast days
direct, sky diffuse and reflected 100% of the solar
components. radiation is
diffuse.

6
Spectral distribution of short-wave solar radiation

NASA/ASTM Standard Spectral Irradiance

Wavelength (μm)

0 - 0.38 0.38 – 0.78 > 0.78

(visible range)
Fraction in range 0.07 0.47 0.46
Energy in range (W/m2) 95 640 618

7
Short-wave radiation
impacts

8
Passive utiulisation

9
Location coordinates

 latitude - angle N or S
above or below equator.

 longitude – angle E or W
from prime meridian
(Greenwich).

 Longitude difference –
angle from location to
local time zone reference
meridian (west –ve).

10
Solar declination

21 December
21 March
summer S
hemisphere

21 June
summer N 21 September 30

hemisphere 20

10
Declination
0
-3 5 65 16 5 2 65 365
-1 0

-2 0

-3 0
D ay o f th e y e ar

11
Solar time

ts – tm = ± Ldiff/15 + (et/60) + ds

where,
ts = solar time
tm = local time
Ldiff = longitude difference
et = equation of time
ds = daylight saving time

12
Solar geometry

 Declination
d = 23.45 sin (280.1 + 0.9863 Y)
where Y = year day number (January 1 =1,
December 31 = 365)

 Altitude
βs = sin-1 [cos L cos d cos θh + sin L sin d ]
where L is site latitude,
θh is hour angle = 15 (12 – ts)

 Azimuth
αs = sin-1 [ cos d sin θh / cos βs ]

 Incidence angle
iβ = cos-1[ sin βs cos (90-βf) + cos βs cos ω sin (90-βf)]
where ω = azimuth angle between sun and surface normal,
βf = surface inclination angle

13
Solar radiation prediction (all W/m2)

Idn - direct normal or “beam” (pyrheliometer)

Idh - direct horizontal Idh = Idnsinβs
known
Ifh - diffuse horizontal (pyranometer with shadow band)
Igh - global horizontal (pyranometer or solarimeter)
rg - ground reflectivity

Idβ - direct radiation on a surface of inclination βf

unknown
Isβ - sky diffuse radiation incident on a surface of inclination βf
Irβ - ground reflected radiation incident on a surface of inclination βf

Igh = Idh +Ifh

= Idn sin βs + Ifh
Solar Altitude, βs
Solar data for simulation:
either: Igh and Ifh or Idn and Ifh

14
Solar radiation measurement

 Pyranometer measures the total solar irradiance

on a planar surface.

 Pyrheliometer measures direct beam solar

radiation by tracking the sun’s position
throughout the day.

15
Solar radiation measurement

 Shaded pyranometer measures diffuse

solar irradiance on a (usually horizontal)
surface.

 The shade blocks direct radiation and

some diffuse radiation (so need to adjust
readings).

 Integrated pyranometer measures both total

and diffuse radiation on a (usually
horizontal) surface.

 Diffuse is calculated based on shading

patterns from internal shades

16
Short-wave flow-paths

A - reflected shortwave flux

B - flux emission by convection and
longwave radiation
C - shortwave flux transmission to cause
opaque surface insolation
D - shortwave transmission to cause
transparent surface insolation
E - shortwave transmission to adjacent zone
F - enclosure reflections
G - shortwave loss
H - solar energy penetration by transient
conduction
I - solar energy absorption prior to
retransmission by the processes of B.

17
Short-wave radiation calculation
iβ - angle between the
incident beam and the
Intensity of direct radiation on surface of inclination β: surface normal vector
Idβ = Idh cos iβ / sin βs ω- surface-solar azimuth
(= |αs − αf|)
Intensity of diffuse radiation on same surface αf, βf - surface azimuth and
ground reflected: Irβ = 0.5 [1- cos (90 – βf)] (Idh + Ifh ) rg inclination respectively
where rg is the ground reflectance αs, βs - solar azimuth and
sky component: Isβ = 0.5 [1+ cos (90 - βf)] Ifh elevation respectively
assuming an isotropic diffuse sky

In practice the sky is not isotropic and so empirically-based models

that correct for circumsolar and horizon brightening are employed:
sky component:  1  cos(90  β f )     I fh 2  
 

3 β f  
I s β  I fh    1  1   2   sin  
 2     I gh    2  
   I 2 fh  
  1  1   2   cos 2 (i β )sin 3 90  β s 
   I  
  gh  

Angle of incidence: i β  cos -1

sin  s cos( 90   f )  cos  s cos  sin( 90   f )  Numerical approach
using 145 sky vault
patches.

18
Surface-solar angles
solar surface
beam normal
N

ψ
βf
cross section
surface inclined at
βs
plan view
angle βf
solar
αs αf
beam

surface 3-D view

normal iβ
ω
solar
beam
S
S βf

19
Solar angle tables (altitude & azimuth)

20
Solar tables (Idv & Idh)

21
PV power output

A simple model:

Example 1 Example 1

Calculate the power output from a PV For the same situation calculate the
panel at 60°C with 840 W/m2 incident power output if the temperature was
solar radiation if the same panel produces 30°C. β is again measured at 0.003 W/K
150 W at STC (1000W/m2 & 25°C). β is
measured at 0.003 W/K

22
Longwave Radiation
Exchange

23
Calculating equation coefficients

Construction Conservation Equation

Surface Conservation Equation

Fluid Conservation Equation

needs flow estimation

needs radiation and

convection estimation

24
Internal long-wave radiation – calculation

25
Internal long-wave radition =ε σA → = A

26
Internal long-wave radiation – numerical method
 Surfaces divided into finite elements and a
unit hemisphere superimposed on each
element.
 Unit hemisphere’s surface divided into patches
representing the radiosity field of the
associated finite element.
 ‘Energy rays’ are formed by connecting the
centre point of the finite element and all
surface patches.
 Each ray is projected to determine an
intersection with another surface.
 At this intersection a surface response model
is invoked to determine the energy absorption
and the number and intensity of exit rays –
these are continually added to the stack of rays
queued for processing.
 Ray processing is discontinued when the
inherent energy level falls below a threshold.
 The energy absorptions for each finite element
are then summated as appropriate to give the
final net longwave radiation exchanges for the
enclosure.
27
External long-wave radiation

28