Imprimer cette page

Calcul des trajectoires du soleil et de l'insolation

Ce logiciel freeware calcule les trajectoires du soleil dans le ciel, en propose plusieurs représentations (3D et projections usuelles), et calcule les insolations instantanées avec possibilité d'exportation des résultats (journalier, mensuel, annuel) sous forme de tableau au format html.

Il fonctionne sans installation, il suffit d'exécuter le fichier téléchargé : AZprocede_solaire.exe. Ce fichier peut par exemple être stocké dans un répertoire quelconque de votre PC (Mes documents, C:\AZprocede), ou sur un support de stockage mobile (clef USB) pour être exécuté n'importe ou.

Il calcule, en temps réel ou en accéleré, pour n'importe quel endroit du globe terrestre, la course du soleil dans le ciel, l'insolation qui en résulte, au sol ou sur un capteur solaire plan d'orientation et d'inclinaison donné. Les équations permettant de calculer toutes les données sont fournies dans la partie "Modèle Numérique". Différents paramètres sont à la main de l'utilisateur permettant d'effectuer toutes sortes de simulation. Plusieurs représentations graphiques sont accessibles pour visualiser au mieux la course du soleil dans le ciel, ou représenter les heures de lever et de coucher du soleil, l'énergie reçue par le capteur, etc... sur une période d'une année.

Les paramètres réglables sont:

  • Modification de la vitesse 1X, 60X (1s=1mn), 600X (1s=10mn) ou 3600X (1s=1h),
  • Latitude et longitude du point d'observation, et longitude de référence pour fuseau horaire (par pas de 15°),
  • Décalage heure d'été et heure d'hiver en manuel ou en automatique (dernier dimanche de mars à 3h, dernier dimanche d'octobre à 3h),
  • Temps solaire ou temps légal, numéro du jour de l'année, année (au démarrage, la date réglée est celle du jour),
  • Pour le capteur solaire: inclinaison i et orientation o réglées manuellement, ou orientation automatique sur l'azimut a, ou inclinaison automatique sur la hauteur h du soleil, ou orientation et inclinaison automatiques,
  • Pour l'atmosphère: température, altitude, facteur de trouble de Linke, humidité relative (%),
  • Temps solaire vrai TSV=12+(AH/15), AH angle horaire du soleil en degrés, 
  • ...

Outils graphiques d’analyse du procédé

Représentation de la course du soleil dans le ciel en temps réel ou accéléré (jusqu'à 3600X - 1s=1h, et V=16X):

  • course du soleil en temps réel dans le ciel (représentation et affichage hauteur et azimut en °)
  • projection cylindrique, projection orthogonale, ou projection sétréographique (projection sur le plan équatorial en passant par le nadir)
  • vue en trois dimensions de la course du soleil (projection axonométrique avec angle réglable et rotation d'azimut),
  • courses du soleil pour les dates particulières suivantes: solstices d'hiver (21 décembre, jour le plus court n° 355 en bleu), solstice d'été (21 juin, jour le plus long n° 172 en rouge), équinoxe de printemps (21 mars, jour n°80 en vert, confondue avec l'équinoxe d'automne le 21 septembre)
  • courbes solaires annuelles avec heures de lever et de coucher de soleil, durée d'insolation, heure de passage au méridien et hauteur maximale pour chaque jour,
  • course du soleil à heure fixe chaque jour de l'année (tracé de l'analemme),
  • ...

Modèle numérique

  • Energie solaire moyenne à l'entrée dans l'atmosphère: E0=1367 W.m-2,
  • Energie solaire en fonction du jour de l'année j: ESol=1367×(1+0.0334×cos(360×(j-2.7206)/365.25)) W.m-2,
  • Pression atmosphérique Patm=101325×(1-2.26.e-5×z)5.26, en Pa, avec z altitude en m,
  • Pression de vapeur saturante P°eau=2.165×(1.098+T/100)8.02, en mmHg, avec T température en °C,
  • Pression partielle de vapeur d'eau Peau=ε×P°eau, en mmHg, avec ε humidité relative en %,
  • Masse d'air optique mair=Patm/(101325×sin(h)+15198.75/(3.885+h)1.253, avec h hauteur du soleil en °,
  • Epaisseur optique de Rayleigh ER=1/(0.9×m+9.4), atténuation par diffusion,
  • Facteur de trouble de Linke TL=2.4+14.6×B+0.4×(1+2×B)×ln(Peau),
  • B varie de 0.02 (montagne) à 0.05 (campagne), 0.10 (ville) et à 0.20 (atmosphère polluée),,
  • Rayonnement solaire direct: I*=ESol×exp(-ER×m×sin(h)), en W.m-2,
  • ou formule simplifiée I*=ESol×exp(-TL/(0.9+9.4×sin(h))), en W.m-2,
  • Orientation du plan récepteur o suivant la rose des vents (0° pour le sud, -90° pour l'est, 90° pour l'ouest, 180° pour le sud),
  • Inclinaison du capteur i 0° pour le plan horizontal, 90° pour le plan vertical,
  • Coefficient d'incidence: CI=sin(i)×cos(h)×cos(orientation-azimut)+cos((i)*sin(h),
  • Rayonnement solaire diffus: D*=125×sin(h)0.4×((1+cos(i)/2)+211.86×sin(h)1.22×((1-cos(i)/2), en W.m-2,
  • Rayonnement solaire global: G*=S*+D*=I*×CI+D*
  • ...
...lire...

texte déroulé