How solar calculators work
Hub
7 July 2023
There is an elaborate mathematical system for calculating electricity yield from solar radiation. In this article, I will briefly describe the core algorithm behind most solar yield calculators.
Note: The SPAC (Solar Panel Angle Calculator) application incorporates several key improvements over this standard baseline algorithm.
Popular online solar yield estimation tools:
- Google Sunroof
- EnergySage Solar Calculator
- NREL PVWatts (a industry benchmark alongside SPAC)
Pyranometer and Pyrheliometer
Calculating basic theoretical solar energy potential for ideal conditions is relatively straightforward. However, real-world yield depends heavily on local cloud cover, atmospheric humidity, and pollution. Historically—when PV investments carried high capital costs—onsite ground irradiance measurements were mandatory to establish precise feasibility data.
Solar irradiation is primarily measured using two distinct sensor types:
- Pyranometers: Measure total global solar radiation (both direct and diffuse) falling on a flat surface.
- Pyrheliometers: Measure exclusively the direct beam component of solar radiation by tracking the sun continuously.
Pyranometers are the most common sensors in solar engineering. Two primary operational designs dominate the market: Li-Cor silicon photodiode sensors and Thermopile black-and-white sensors.
Pyranometer (Global Solar Irradiance)
Pyrheliometer (Direct Normal Irradiance)
The Working Principle
Black-and-white thermopile pyranometers detect radiation through temperature differentials between alternating black (absorbing) and white (reflecting) segments. While thermopiles capture the entire solar spectrum evenly, silicon photodiode pyranometers (like the popular Li-Cor series) have wavelength-dependent spectral responses—typically measuring up to $1100\text{ nm}$. Because crystalline silicon PV modules share a similar spectral response cutoff, silicon pyranometers became an industry standard for PV site assessments.
To estimate the energy collected on a tilted panel over a specific timeframe (e.g., monthly), solar calculators first decompose total measured horizontal irradiance ($I_h$) into its direct horizontal component ($I_{bh}$) and diffuse component ($I_{dh}$):
Clear Sky Index ($K_{th}$)
Decomposing horizontal radiation requires calculating the clearness index ($K_{th}$). It represents the ratio of measured horizontal insolation on Earth's surface ($I_h$) to theoretical extraterrestrial horizontal insolation ($I_o$) at the same latitude and time:
A higher $K_{th}$ value indicates clear, cloudless skies and low atmospheric turbidity, whereas lower values represent overcast or polluted atmospheric conditions.
Mean Horizontal Extraterrestrial Insolation
Integrating extraterrestrial solar radiation from sunrise to sunset projected onto a flat surface yields daily extraterrestrial horizontal irradiation ($I_o$):
- $SC$ – Solar Constant ($\approx 1367 \text{ W/m}^2$)
- $n$ – Day number of the year ($1$ to $365$)
- $L$ – Geographic Latitude
- $H_{sr}$ – Sunrise hour angle (in radians)
- $\delta$ – Solar declination angle for the given period
Once $K_{th}$ is derived, empirical correlations (such as the Liu-Jordan method) are applied to calculate the diffuse ($I_{dc}$) and ground-reflected ($I_{rc}$) radiation components reaching a tilted surface at slope angle $\Sigma$:
Direct beam radiation hitting the tilted collector ($I_{bc}$) depends heavily on the solar incidence angle ($\theta$), which is governed by panel tilt, azimuth, solar elevation angle, and solar azimuth.
Direct horizontal radiation ($I_{bh}$) relates to direct normal irradiance ($I_b$) via solar altitude angle ($\alpha$):
Thus, direct beam radiation incident on the tilted module surface ($I_{bc}$) is expressed as:
Slope Factor ($R_b$)
Combining the direct radiation equations yields the relationship between direct horizontal radiation and direct radiation on the tilted module surface:
The geometric ratio $R_b$ is known as the slope factor. For short time steps (e.g., 10-minute intervals), $R_b$ is calculated directly from instantaneous sun position geometry.
For daily or monthly averages based on summary irradiance data ($I_h$), an integrated mean slope factor ($\overline{R}_b$) must be calculated:
Implementation example of the solar geometric calculations
PV Module Modeling & Thermal Losses
Next, electrical output is modeled based on specific module parameters (STC rating, NOCT, temperature coefficients) and local ambient conditions.
Because cell efficiency decreases as operating temperature rises above Standard Test Conditions ($25^\circ\text{C}$), average operating cell temperature ($T_{cell}$) is calculated using Nominal Operating Cell Temperature (NOCT):
Finally, DC array output is converted to net AC power output by applying systemic derating factors (inverter efficiency, mismatch, wiring resistance, and module soiling):
Final Energy Calculation
Assuming an MPPT-enabled inverter, total expected energy output ($W$) over a given timeframe is derived as:
This represents the standard calculation workflow used by most traditional solar web tools. However, to eliminate excessive monthly approximations and account for exact incident angle modifier dynamics, I developed the SPAC (Solar Panel Angle Calculator) algorithm to perform step-by-step vector calculations for superior accuracy.
Related posts
Off Grid DC Voltage
Selecting the correct DC system voltage for an off-grid photovoltaic (PV)...
Read moreSPAC 1.1 UI
The Solar Panel Angle Calculator (SPAC) is a web application designed to compute optimal tilt...
Read moreSPAC 1.27 UI
The Solar Panel Angle Calculator (SPAC) has evolved from a simple tilt-angle computation tool...
Read more