pvperfsolar performance modeling

Reference / Perez settings

Model option

PVsyst 8 settings

1361 W/m², the direct normal irradiance computed from the global and diffuse values, a 0.9 cap and the Kasten 1966 air mass

Inputs used: Direct normal irradiance, Station pressure, Altitude

Method

The Perez equations run with the inputs and constants below. The solar constant is 1361 W/m². The direct normal irradiance is computed from the file's global and diffuse horizontal values; the file's direct normal column is not used. The direct sunlight and the circumsolar part are both zero while the sun is less than 2° above the horizon. Before it enters the sky clearness, the direct normal irradiance is capped at 0.9 of the extraterrestrial irradiance. The sky brightness uses the Kasten 1966 air mass and an altitude factor that decreases with the site elevation. The circumsolar coefficient has no floor: under some overcast skies it is negative, which makes the circumsolar part negative and raises the uniform part.

Settings

Setting Value
Solar constant SS 1361 W/m²
Eccentricity factor ete_t 1+0.033cos⁡2πt87601 + 0.033\cos\frac{2\pi t}{8760}, with tt in hours
Direct normal irradiance DNIused\text{DNI}_{used} (Gh−D)/sin⁡α(G_h - D)/\sin\alpha from the file's global and diffuse horizontal values; zero when α<2°\alpha < 2°
Cap on the direct part of the clearness cεc_\varepsilon, as a fraction of E0E_0 0.9
Air mass AM\text{AM} Kasten 1966
Altitude factor fzf_z e−1.30×10−4 ze^{-1.30 \times 10^{-4}\,z}, with zz the site elevation in m
Floor of the circumsolar coefficient F1,minF_{1,min} none; negative values are kept
Circumsolar cutoff elevation αc\alpha_c 2°
Coefficient table Perez 1990, all-sites composite, 8 clearness bins

The math

  1. Extraterrestrial irradiance
  2. Direct sunlight
  3. Sky clearness
  4. Sky brightness
  5. Circumsolar and horizon-band coefficients
  6. Uniform part
  7. Circumsolar part
  8. Horizon-band part
  9. Diffuse sky light on the plane
E0=1361 (1+0.033cos⁡2πt8760)(1)E_0 = 1361\,\left(1 + 0.033\cos\frac{2\pi t}{8760}\right) \tag{1} DNIused={Gh−Dsin⁡αα≥2°0α<2°,B=DNIused max⁡(0, cos⁡θ)(2)\text{DNI}_{used} = \begin{cases} \dfrac{G_h - D}{\sin\alpha} & \alpha \ge 2° \\ 0 & \alpha < 2° \end{cases},\qquad B = \text{DNI}_{used}\,\max(0,\ \cos\theta) \tag{2} ε=(D+DNIε)/D+1.041 θz31+1.041 θz3,DNIε=min⁡ ⁣(DNIused, 0.9 E0)(3)\varepsilon = \frac{(D + \text{DNI}_\varepsilon)/D + 1.041\,\theta_z^3}{1 + 1.041\,\theta_z^3},\qquad \text{DNI}_\varepsilon = \min\!\left(\text{DNI}_{used},\ 0.9\,E_0\right) \tag{3} Δ=D AM fzE0,AM=1sin⁡α+0.15 (α+3.885)−1.253,fz=e−1.30×10−4 z(4)\Delta = \frac{D\,\text{AM}\,f_z}{E_0},\qquad \text{AM} = \frac{1}{\sin\alpha + 0.15\,(\alpha + 3.885)^{-1.253}},\qquad f_z = e^{-1.30 \times 10^{-4}\,z} \tag{4} F1=f11+f12 Δ+f13 θz,F2=f21+f22 Δ+f23 θz(5)F_1 = f_{11} + f_{12}\,\Delta + f_{13}\,\theta_z,\qquad F_2 = f_{21} + f_{22}\,\Delta + f_{23}\,\theta_z \tag{5} Diso=D (1−F1) 1+cos⁡β2(6)D_{iso} = D\,(1 - F_1)\,\frac{1 + \cos\beta}{2} \tag{6} Dcirc={D F1 max⁡(0, cos⁡θ)max⁡(sin⁡5°, sin⁡α)α≥2°0α<2°(7)D_{circ} = \begin{cases} D\,F_1\,\dfrac{\max(0,\ \cos\theta)}{\max(\sin 5°,\ \sin\alpha)} & \alpha \ge 2° \\ 0 & \alpha < 2° \end{cases} \tag{7} Dhb=D F2sin⁡β(8)D_{hb} = D\,F_2\sin\beta \tag{8} Dplane=Diso+Dcirc+Dhb(9)D_{plane} = D_{iso} + D_{circ} + D_{hb} \tag{9}

where

Symbol Meaning Unit
E0E_0 the extraterrestrial normal irradiance at time tt W/m²
tt the time from 1 January 00:00 hours
DNIused\text{DNI}_{used} the direct normal irradiance the model uses W/m²
GhG_h the global horizontal irradiance W/m²
DD the diffuse horizontal irradiance of the weather file W/m²
α\alpha the elevation of the sun above the horizon degrees
BB the direct sunlight on the module plane W/m²
θ\theta the angle between the sun and the perpendicular to the module plane degrees
ε\varepsilon the sky clearness, which selects the row of the coefficient table —
DNIε\text{DNI}_\varepsilon the direct normal irradiance that enters the sky clearness W/m²
θz\theta_z the zenith angle of the sun, 90° minus α\alpha radians
Δ\Delta the sky brightness —
AM\text{AM} the relative air mass of Kasten 1966 —
fzf_z the altitude factor of the sky brightness —
zz the site elevation m
F1F_1 the circumsolar coefficient —
F2F_2 the horizon-band coefficient —
f11…f23f_{11} \ldots f_{23} the six coefficients of the clearness row, from the published table —
DisoD_{iso} the uniform part of the diffuse light on the plane W/m²
β\beta the tilt of the module plane from the horizontal degrees
DcircD_{circ} the circumsolar part of the diffuse light on the plane W/m²
DhbD_{hb} the horizon-band part of the diffuse light on the plane W/m²
DplaneD_{plane} the diffuse sky light on the module plane W/m²

References

  1. Perez, Ineichen, Seals, Michalsky & Stewart 1990, Modeling daylight availability and irradiance components from direct and global irradiance, Solar Energy 44(5) 271–289
  2. Kasten 1966, A new table and approximation formula for the relative optical air mass, Archiv für Meteorologie, Geophysik und Bioklimatologie B14 206–223
  3. PVsyst 8 help, Transposition model
  4. Analysis of PVsyst 8.1.5 outputs shows this convention in its results.

API

array.perez_settings = pvsyst_8

Symbols follow the symbol registry.