Uitleg

Bericht van: Robert (Sint Philipsland) , 13-07-2022 16:09 

De gemeten globale straling G is de som van de directe straling S en de diffuse straling D:

G = S + D

Door de diffuse straling van de globale straling af te trekken krijg je dus de directe straling:

S = G - D

Volgens vergelijking (2.8) uit het technisch rapport kan tijdens wisselend bewolkt weer de diffuse straling berekend worden uit het minimum van de globale straling in een tijdsinterval (wanneer er dus een wolk voor de zon schuift en de directe straling zo goed als wegvalt):

D = (1+a)*Gmin

De theoretische directe straling I kan berekend worden als (zie vergelijking (2.4) en (2.6) uit het technisch rapport, let op dat er in (2.6) een paar haakjes vergeten zijn)

I = I0 * sin(h) * exp(-TL/(0.9 + 9.4*sin(h)))

waarbij I0 de zonneconstante is (ca. 1350 W/m^2), h de hoogte van de zon boven de horizon (in graden of radialen) en TL de Linke troebelheidsfactor (typisch 2-4). De laatste factor corrigeert voor de verzwakking van zonlicht door kleine deeltjes. Deze verzwakking neemt toe naarmate het zonlicht een langere afstand door de atmosfeer moet afleggen. De fractie zonneschijnduur is dan de geschatte directe straling gedeeld door de theoretische directe straling:

Fr = S/I = (G - D)/I = (G - (1+a)*Gmin) / (I0 * sin(h) * exp(-TL/(0.9 + 9.4*sin(h)))

Bovenstaande vergelijking is identiek aan (2.11) uit het technisch rapport, behalve dat ze daar wederom een paar haakjes vergeten zijn. In de overzichten van de algoritmes van Slob en Schipper staat het trouwens wel goed. Voor G wordt de gemiddelde globale straling in een tijdsinterval gebruikt. De waarden van de parameters die zij gebruiken zijn:

a = 0.2 en TL = 4 (Slob)

of

a = 0.27 en TL = 13.03 (Schipper).

Wat kan hier fout gaan? Stel dat de gemiddelde G bijna gelijk is aan de minimale G (weinig variatie in bewolking i.p.v. wisselend bewolkt). De geschatte diffuse straling D wordt dan 1.2 (Slob) of 1.27 (Schipper) keer zo groot als de totale straling G! Dat zou betekenen dat de hoeveelheid directe straling negatief zou zijn. De zonneschijnduur wordt dan ook negatief en dat kan natuurlijk niet. Daarom kan Fr niet lager worden dan 0. Maar eigenlijk past deze situatie beter bij geheel onbewolkt (dun laagje sluierbewolking) of geheel bewolkt.

Een andere situatie waarin het fout gaat is wanneer de geschatte directe straling veel groter is dan de theoretische directe straling. Dit is het geval wanneer het verschil tussen de gemiddelde G en de minimale G juist heel groot is. Vooral de methode van Schipper heeft hier last van, omdat in deze methode een erg hoge troebelheidsfactor gebruikt wordt waardoor de theoretische directe straling vermoedelijk te laag ingeschat wordt. Mogelijk is dit om het effect van sluierbewolking na te bootsen. Hoe dan ook, de zonneschijnfractie kan niet hoger zijn dan 1.

Volgens de inleiding van het technisch rapport is er een nauwkeurigere methode om de zonneschijnduur te bepalen door niet alleen de globale straling, maar ook de directe straling te meten. Voor deze methode was echter een dure zonnevolger nodig die gedeeltelijk handmatig bediend moest worden dus vandaar dat alleen de globale straling gemeten wordt op de automatische stations. Maar het zou goed kunnen dat het met de huidige stand van de techniek wel haalbaar is om ook de directe zonnestraling te meten.

Bericht laatst bijgewerkt: 13-07-2022 16:37

Zonneschijnduur berekenen algoritme Schipper   ( 585)
Luc (Kessel Lb) ( 25m) -- 13-07-2022 11:42
Re: Zonneschijnduur berekenen algoritme Schipper   ( 266)
Peter (Wiltz -Luxemburg) ( 381m) -- 13-07-2022 12:19
Re: Zonneschijnduur berekenen algoritme Schipper   ( 241)
Luc (Kessel Lb) ( 25m) -- 13-07-2022 12:33
Re: Zonneschijnduur berekenen algoritme Schipper   ( 304)
Mark (Kaulille, B Limburg) ( 47m) -- 13-07-2022 12:47
Re: Zonneschijnduur berekenen algoritme Schipper   ( 292)
Luc (Kessel Lb) ( 25m) -- 13-07-2022 13:14
Antwoord uit Handboek Waarnemingen, hfdstk. 8   ( 305)
Ben (Lelystad) ( 13m) -- 13-07-2022 14:04
Top! Bedankt voor de info.   ( 208)
Luc (Kessel Lb) ( 25m) -- 13-07-2022 14:19
Als je data wilt van zowel globale staling als (gemeten) zonuren...   ( 234)
Gerhard (Vaassen) ( 15m) -- 13-07-2022 16:42
Re: Zonneschijnduur berekenen algoritme Schipper   ( 236)
Mark (Kaulille, B Limburg) ( 47m) -- 13-07-2022 14:20
Re: Zonneschijnduur berekenen algoritme Schipper   ( 216)
Hendrik (Londerzeel-B) ( 7m) -- 13-07-2022 14:38
Re: Zonneschijnduur berekenen algoritme Schipper   ( 184)
Mark (Kaulille, B Limburg) ( 47m) -- 13-07-2022 14:58
Niet helemaal   ( 200)
Luc (Kessel Lb) ( 25m) -- 13-07-2022 15:24
Re: Niet helemaal   ( 171)
Mark (Kaulille, B Limburg) ( 47m) -- 13-07-2022 15:31
Re: Niet helemaal   ( 173)
Luc (Kessel Lb) ( 25m) -- 13-07-2022 15:58
Het zou hier maar donker zijn   ( 169)
Hendrik (Londerzeel-B) ( 7m) -- 13-07-2022 18:19
Ik ben geen held in berekeningen   ( 219)
Hendrik (Londerzeel-B) ( 7m) -- 13-07-2022 13:07
Uitleg   ( 282)
Robert (Sint Philipsland) -- 13-07-2022 16:09
Re: Uitleg   ( 256)
Rutger (Meppel) -- 13-07-2022 16:14
Interessante data en artikel Rutger   ( 191)
Robert (Sint Philipsland) -- 13-07-2022 17:49
Re: Zonneschijnduur berekenen algoritme Schipper   ( 223)
Luc (Kessel Lb) ( 25m) -- 14-07-2022 11:18
Re: Zonneschijnduur berekenen algoritme Schipper   ( 179)
Mark (Kaulille, B Limburg) ( 47m) -- 14-07-2022 11:37
Re: Zonneschijnduur berekenen algoritme Schipper   ( 197)
Luc (Kessel Lb) ( 25m) -- 14-07-2022 13:30
Re: Zonneschijnduur berekenen algoritme Schipper   ( 167)
Mark (Kaulille, B Limburg) ( 47m) -- 14-07-2022 13:50
Bedankt!   ( 168)
Luc (Kessel Lb) ( 25m) -- 14-07-2022 15:22
Re: Bedankt!   ( 163)
Mark (Kaulille, B Limburg) ( 47m) -- 14-07-2022 15:48
Even meerekenen...   ( 254)
Ben (Lelystad) ( 13m) -- 14-07-2022 14:49
Re: Even meerekenen...   ( 223)
Luc (Kessel Lb) ( 25m) -- 14-07-2022 15:22