Dachamatics

Direct solar radiation and shading

Duration of direct solar illumination of a bed over a day, and the share of the daily direct-beam flux that reaches it in the presence of shading objects.

Quantities

SymbolQuantityUnit
φ, λlatitude and longitude of the plotdegrees
hsolar elevation above the horizondegrees
zzenith distance, z = 90° − hdegrees
mrelative optical air massdimensionless
Ishare of the direct beam reaching a horizontal surfacedimensionless
s(t)share of the bed's area not shaded at time t0 to 1
Δtintegration step0.25 h
Tduration of direct illuminationh/day
Eshare of the daily flux reaching the beddimensionless

Formulas

Optical air mass (Kasten & Young, 1989)

m = 1 / (cos z + 0.50572 · (96.07995 − z)^(−1.6364))

Attenuation of the direct beam (Meinel & Meinel, 1976)

I = 0.7^(m^0.678) · cos z

The cos z factor accounts for the inclination of the beam on a horizontal surface. I is normalised to the flux at z = 0 outside the atmosphere and is not expressed in W/m².

Duration of illumination

T = Σ Δt · s(t), summed over the moments t for which h > 0

Share of the daily flux

E = Σ I(t) · Δt · s(t)  ÷  Σ I(t) · Δt

The denominator is evaluated for the same bed on the same day with s(t) ≡ 1. E = 1 denotes the absence of shading and is attainable on any day of the year.

Solar position

Elevation and azimuth follow the NOAA Solar Calculator algorithm: solar declination, equation of time, hour angle. The equation of time shifts true noon relative to mean noon by up to ±16 min depending on the date; omitting it displaces the shadow by the corresponding amount.

Shadow geometry

A tree crown is approximated by a sphere centred at height H − R, where H is the tree height and R the crown radius; its shadow on a horizontal plane is an ellipse elongated by 1/sin h. A building is approximated by a prism up to the eaves; when a ridge height is given, a gable roof is added above it. The building's shadow is its own footprint displaced away from the Sun, together with the roof's shadow: each footprint edge yields a quadrilateral joining that edge's shadow cast from eaves height to the shadow of its projection onto the ridge cast from ridge height. For a rectangular footprint with the ridge along the long axis this gives two slopes and two vertical gable ends. The transmittance of crown and wall is taken as zero.

Illuminated share of the area

The bed polygon is covered by a regular 5 × 5 grid of points. At each point the transmittances of all covering shadows are multiplied; s(t) is the mean over the points. Shading is evaluated for h ≥ 5°: at lower elevations shadow length grows without bound and the direct component becomes negligible.

Applicability and limitations

Sources

← Back to contents