Saturday 15 September 2012

Energy balance

Figure 1shows the energy sources that are considered when estimating the energy balance at the pavement surface. The energy that enters (or leaves) the ground comes from different sources such as the sun, the objects that surround the pavement and the air. This energy balance can be expressed as (Dempsey, 1985):
Qg = Qs + Ql + Qc                                                                                                                                            (1)
where
Qg is the energy absorbed by the ground [W m-2]
Qs is the net short wave radiation [W m-2]
Ql is the net long wave radiation [W m-2]
Qc is the convective heat transfer [W m-2]

clip_image002[4]

Figure 1: Energy exchange between pavement surface and environment
It has to be noted that this approach does not include the effects of transpiration, condensation, evaporation and sublimation.
Qs – net short wave radiation
The short wave radiation can be identified with the high frequency energy radiated by the sun that reaches the pavement surface and part of which is visible to the human eye. It is, therefore, a function of the extraterrestrial radiation incident on the atmosphere at a particular latitude, day and time and of the filtering and reflecting effect of the clouds that may be present at that location (ME-PDG, 2004).
 clip_image005                                                                                                                       (2)
where
as is the surface short wave absorptivity (0.8 – 0.98)
R is the potential extraterrestrial radiation [W m-2]
Sc is the percentage of sunshine
A and B are constants that account for diffuse scattering and absorption by the atmosphere. The ME-PDG suggests values of 0.202 and 0.539 respectively but they will change with location.
The potential extraterrestrial radiation R is a function of latitude Ф [rad], declination δ [rad] and time of day. Its average value during one day, Ra, can be calculated by:
 clip_image008                             (3)
where
Ф is the latitude expressed in radians
δ is the declination expressed in radians and calculated by:
 clip_image011                                                                                             (4)
where day is the day of the year expressed as an integer.
The average value of Ra does not serve the purpose of estimating temperatures on an hourly basis. Nonetheless, assuming a parabolic distribution of R during daytime that starts from sunrise and finishes at sunset (Dempsey, 1985), it is possible to use Ra to estimate R at any given time. The total energy radiated during daytime must be equal to Ra multiplied by 24h, i.e. the areas of Figure 2(a) and (b) must be the same.
clip_image012 clip_image015

Figure 2: Calculating R from Ra (a) (b)
The time of sunrise and sunset for one particular day and latitude can be calculated as follows (University of Oregon, 2009):
Sunrise_angle = clip_image018 [rad]                                                                               (5)
Sunset_angle = clip_image021 [rad]                                                                        (6)
Sunrise_time = clip_image024 [s]                                                              (7)
Sunset_time = clip_image027 [s]                                                               (8)
From what given above, the parabola’s coefficients can be calculated:
clip_image030                                                                                       (9)
clip_image033                                                                                         (10)
clip_image036 [s]                                                                                    (11)
clip_image039 for clip_image001                  (12)
R=0 for  clip_image005[4] or clip_image008[22](13)
Ql – net long wave radiation
Any object that is at a temperature higher than the absolute zero radiates energy. Since we are dealing with objects at low temperatures (certainly not comparable to the sun’s temperature), they radiate at low frequencies / long wave lengths, well outside the visible spectrum.
The general formula to calculate this energy is given by:
clip_image050 [W m-2]                                                                                                                         (14)
where:
ε is the object’s emissivity
σ is the Stefan-Boltzman constant [W m-2 K-4]
T is the object’s temperature [K]
The net long wave radiation that contributes to the energy balance at the pavement’s surface is calculated considering the incoming (positive) energy radiated by the air and the outgoing (negative) energy radiated by the pavement itself. Therefore:
clip_image053 [W m-2]                                                                                                          (15)
where:
Rldc is the clear sky downward radiation:
clip_image056                                                                                                                                         (16)
where the air emissivity εa can be calculated as a function of vapour pressure and Ta is the air temperature (Jacobs, 2004).
Qx is the outgoing radiation:
clip_image059                                                                                                                                           (17)
where εs is the surface emissivity and Ts is the surface temperature.
Ccc is the cloud cover correction factor:
clip_image062                                                                                                                                      (18)
where N is the cloud base factor (0.9-0.8 depending on height) and W is the percentage of clouds.
Qc – convective heat transfer
The pavement’s surface exchanges heat with the surrounding air by convection. This energy transfer is function of air and surface temperature and wind speed. The ME-PDG suggests the following equation:
clip_image065                                                                                    (19)
where:
Tm = average of Ta and Ts [K]
U = average daily wind speed [m s-1]

0 comments:

Post a Comment