Brevet FR3104241A1 — publié le 20210611

PROCEDE DE DETERMINATION D’UN PARAMETRE D’EFFICACITE D’UN CAPTEUR SOLAIRE THERMIQUE AYANT L’ENROBE D’UNE ROUTE POUR ABSORBEUR Un aspect de l’invention concerne un procédé de détermination d’un paramètre d’efficacité d’un capteur (200) solaire thermique comportant au moins un tube (201) véhiculant un fluide caloporteur situé dans un absorbeur configuré pour capter l’énergie solaire et transmettre l’énergie solaire captée au fluide caloporteur sous forme de chaleur, l’absorbeur étant un enrobé (202) d’une route et comportant les étapes suivantes : déterminer des conditions extérieures s’exerçant sur le capteur (200) ; calculer au moins une température du fluide caloporteur en sortie (2012) du tube (201) et une température moyenne de l’enrobé (202) à l’aide d’un modèle numérique du capteur (200) à partir des conditions extérieures déterminées et de paramètres de fonctionnement du capteur (200) fixés ; calculer le paramètre d’efficacité à partir de la température du fluide caloporteur en sortie du tube (201) et/ou de la température moyenne de l’enrobé (202). Figure à publier avec l’abrégé : Figure 3A PROCEDE DE DETERMINATION D’UN PARAMETRE D’EFFICACITE D’UN CAPTEUR SOLAIRE THERMIQUE AYANT L’ENROBE D’UNE ROUTE POUR ABSORBEUR DOMAINE TECHNIQUE DE L’INVENTION Le domaine technique de l’invention est celui des capteurs solaires thermiques et plus particulièrement celui des capteurs solaires thermiques ayant pour absorbeur l’enrobé d’une route. La présente invention concerne un procédé de détermination d’au moins un paramètre d’efficacité d’un capteur solaire thermique et en particulier un procédé de détermination d’au moins un paramètre d’efficacité d’un capteur solaire thermique ayant l’enrobé d’une route pour absorbeur. ARRIERE-PLAN TECHNOLOGIQUE DE L’INVENTION De nombreux modèles numériques reproduisant le comportement thermique des capteurs solaires thermiques cohabitent sous différents environnements de simulation. Ces modèles s’appuyant principalement sur la norme internationale NF EN ISO 9806, permettent de calculer la température moyenne du fluide caloporteur circulant dans le capteur et s’échauffant au contact de l’absorbeur captant l’énergie solaire, pour ainsi évaluer la capacité du capteur à convertir l’énergie solaire en énergie thermique. Ces modèles reposent sur la principale hypothèse que la température moyenne du fluide caloporteur est à tous moments identique à la température de l’absorbeur. Or, cette hypothèse n’est pas vérifiée dans le cas d’un capteur solaire thermique particulier ayant pour absorbeur l’enrobé d’une route, l’écart entre la température de l’absorbeur et la température du fluide caloporteur pouvant atteindre 30°C. Ces modèles ne permettent donc pas de calculer de manière précise la température moyenne du fluide caloporteur et ne permettent pas de calculer la température de l’absorbeur du capteur solaire thermique particulier. Il existe donc un besoin d’évaluer de manière précise la capacité d’un capteur thermique solaire ayant pour absorbeur l’enrobé d’une route, à convertir l’énergie solaire en énergie thermique. L’invention offre une solution aux problèmes évoqués précédemment, en permettant de calculer un paramètre d’efficacité d’un capteur solaire thermique ayant pour absorbeur l’enrobé d’une route. Un premier aspect de l’invention concerne un procédé de détermination par un système informatique d’au moins un paramètre d’efficacité d’un capteur solaire thermique comportant au moins un tube véhiculant un fluide caloporteur d’une entrée du tube vers une sortie du tube, le tube étant situé dans un absorbeur configuré pour capter l’énergie solaire et transmettre l’énergie solaire captée au fluide caloporteur sous forme de chaleur, le paramètre d’efficacité étant représentatif de la capacité du capteur à convertir l’énergie solaire en énergie thermique, l’absorbeur étant un enrobé d’une route et comportant les étapes suivantes : déterminer des conditions extérieures s’exerçant sur le capteur ; calculer au moins une température du fluide caloporteur en sortie du tube et une température moyenne de l’enrobé à l’aide d’un modèle numérique du capteur configuré pour calculer au moins la température moyenne de l’enrobé et la température du fluide caloporteur en sortie du tube à partir des conditions extérieures déterminées et de paramètres de fonctionnement du capteur fixés ; calculer le paramètre d’efficacité à partir de la température du fluide caloporteur en sortie du tube et/ou de la température moyenne de l’enrobé. Grâce à l’invention, il est possible d’évaluer précisément la capacité d’un capteur solaire thermique ayant pour absorbeur l’enrobé d’une route, à convertir l’énergie solaire en énergie thermique en utilisant un modèle numérique propre à ce type de capteurs, possédant deux inconnues, la température du fluide caloporteur en sortie du tube et la température moyenne de l’enrobé, ces températures n’étant pas identiques la plupart du temps. Le paramètre d’efficacité calculé à partir des températures obtenues par le modèle numérique pourra par exemple permettre de prendre une décision concernant le déploiement ou non d’un tel capteur à un endroit donné. Le modèle numérique permet également de faire varier les valeurs des paramètres de fonctionnement, par exemple pour déterminer les paramètres de fonctionnement optimaux permettant de maximiser le paramètre d’efficacité. Outre les caractéristiques qui viennent d’être évoquées dans le paragraphe précédent, le procédé selon un premier aspect de l’invention peut présenter une ou plusieurs caractéristiques complémentaires parmi les suivantes, considérées individuellement ou selon toutes les combinaisons techniquement possibles. Selon une variante de réalisation, le paramètre d’efficacité est calculé sur une période de simulation, les étapes du procédé étant réalisées pour chaque pas de temps de la période de simulation. Ainsi, il est possible d’évaluer la capacité du capteur solaire thermique à convertir l’énergie solaire en énergie thermique sur une période donnée, par exemple une semaine, un mois ou une année. Selon une sous-variante de réalisation de la variante de réalisation précédente, le pas de temps est inférieur à une heure. Ainsi, il est possible de prendre en compte les variations dynamiques des conditions extérieures entre deux pas de temps successifs. Selon une sous-variante de réalisation de la variante de réalisation précédente compatible avec la sous-variante de réalisation précédente, le paramètre d’efficacité peut être choisi parmi au moins le groupe formé des paramètres suivants : rendement énergétique du capteur sur la période de simulation, puissance ou énergie thermique produite par le capteur sur la période de simulation, différence de températures maximale entre l’entrée du tube et la sortie du tube sur la période de simulation, différence de températures moyenne entre l’entrée du tube et la sortie du tube sur la période de simulation, différence de températures minimale entre l’entrée du tube et la sortie du tube sur la période de simulation, température maximale de l’enrobé sur la période de simulation, température moyenne de l’enrobé sur la période de simulation, température minimale de l’enrobé sur la période de simulation, température maximale du fluide caloporteur en sortie du tube sur la période de simulation, température moyenne du fluide caloporteur en sortie du tube sur la période de simulation, température minimale du fluide caloporteur en sortie du tube sur la période de simulation ou nombre de jours par période de simulation pour lesquels le rendement énergétique est supérieur à un seuil de rendement. Ainsi, le paramètre d’efficacité est représentatif de la capacité du capteur à convertir l’énergie solaire en énergie thermique. Selon une variante de réalisation compatible avec la variante de réalisation précédente, les conditions extérieures comportent au moins l’irradiance solaire, l’irradiance thermique par rayonnement provenant de la voûte céleste, la vitesse du vent et la température de l’air ambiant autour du capteur. Ainsi, le modèle thermique peut prendre en compte l’énergie solaire stockée par le capteur et au moins les pertes du capteur par rayonnement et par convection. Selon une variante de réalisation compatible avec les variantes de réalisation précédentes, les paramètres de fonctionnement comportent au moins la température du fluide caloporteur à l’entrée du tube et le débit massique du fluide caloporteur dans le tube. Selon une sous-variante de réalisation de la variante de réalisation précédente, la surface du capteur est divisée en une pluralité de surfaces élémentaires identiques comportant chacune une portion du tube et le modèle numérique est composé d’une pluralité de modèles numériques élémentaires identiques correspondant chacun à une surface élémentaire, le modèle numérique élémentaire étant configuré pour calculer au moins la température moyenne de l’enrobé de la surface élémentaire et la température du fluide caloporteur en sortie de la portion de tube comprise dans la surface élémentaire en utilisant la température du fluide caloporteur en sortie de la portion de tube précédente calculée par le modèle élémentaire précédent pour température du fluide caloporteur en entrée de la portion du tube comprise dans la surface élémentaire. Ainsi, la surface étendue caractéristique de ce type de capteurs est découpée en surfaces élémentaires et le modèle numérique est composée d’une pluralité de modèles numériques élémentaires identiques en série s’appliquant chacun à une surface élémentaire. Il est alors possible d’obtenir la température moyenne de l’enrobé et la température en sortie de la portion du tube, et donc le paramètre d’efficacité pour chaque surface élémentaire, et la surface du capteur pourra donc être choisie par exemple pour maximiser le paramètre d’efficacité. Selon une première alternative de réalisation de la sous-variante de réalisation précédente, le procédé comporte une étape préalable de calibration du modèle numérique élémentaire en utilisant au moins une température moyenne de l’enrobé mesurée. Ainsi, le modèle numérique n’a pas besoin d’être calibré sur la surface totale du capteur, mais seulement sur une surface élémentaire de petite dimension, ce qui facilite la mise en œuvre de la calibration du modèle numérique. Selon une deuxième alternative de réalisation de la sous-variante de réalisation précédente, le procédé comporte une étape préalable de calibration du modèle numérique en utilisant au moins une température moyenne de l’enrobé mesurée. Ainsi, le modèle numérique est calibré sur la surface totale, par exemple à partir de données mesurées sur un démonstrateur, et la précision du modèle numérique est améliorée. Selon une sous-alternative de réalisation de l’alternative de réalisation précédente, l’étape de calibration permet de calibrer au moins les paramètres suivants : le coefficient d’absorption de l’enrobé, le coefficient de transfert thermique surfacique entre le fluide caloporteur et l’enrobé, le coefficient d’échange convectif du capteur, l’émissivité globale en rayonnement de la surface du capteur et la capacité thermique de l’enrobé. Ainsi, le modèle thermique peut prendre en compte l’énergie solaire stockée par le capteur et au moins les pertes du capteur par rayonnement, par convection et par inertie et les échanges entre l’enrobé et le fluide caloporteur. Selon une variante de réalisation compatible avec les variantes de réalisation précédentes, l’étape de détermination des conditions extérieures s’exerçant sur le capteur comporte une sous-étape de récupération de conditions extérieures moyennes s’exerçant dans une localité dans laquelle est situé le capteur, à partir d’une base de données météorologiques et une sous-étape de correction des conditions extérieures moyennes en fonction de la localisation du capteur dans la localité. Ainsi, les conditions extérieures s’exerçant sur le capteur sont précises et obtenues sans avoir à déployer des capteurs à la position envisagée pour le capteur. Selon une variante de réalisation compatible avec les variantes de réalisation précédentes, le modèle numérique du capteur comporte les équations suivantes : Avec : E sol l’énergie solaire globale incidente sur le capteur ; E r le flux énergétique net échangé par rayonnement en grandes longueurs d’ondes entre la voûte céleste et le capteur ; E conv les pertes énergétiques par convection E conv du capteur ; E in l’énergie stockée ou déstockée par inertie ; E p l’énergie thermique produite par le capteur ; C e la capacité thermique de l’enrobé ; T e la température moyenne de l’enrobé ; α e le coefficient d’absorption de l’enrobé ; A g la surface d’échange effective du capteur ; I sol l’irradiance globale solaire reçue par le capteur ; ε r l’émissivité globale en rayonnement de la surface du capteur ; I r l’irradiance thermique provenant de la voûte céleste reçue par le capteur; σ la constante de Stefan-Bolztmann ; K(u) un coefficient d’échange convectif dépendant de la vitesse du vent u ; T a la température de l’air ambiant autour du capteur ; q le débit massique de fluide caloporteur dans le tube ; c pf la chaleur spécifique du fluide caloporteur ; T out la température du fluide caloporteur en sortie du tube ; T in la température du fluide caloporteur à l’entrée du tube ; h t le coefficient de transfert thermique surfacique entre le fluide caloporteur et l’enrobé; S t la surface intérieure du tube. Un deuxième aspect de l’invention concerne un produit-programme d’ordinateur comprenant des instructions qui, lorsque le programme est exécuté par un ordinateur, conduisent celui-ci à mettre en œuvre les étapes du procédé selon le premier aspect de l’invention. Un troisième aspect de l’invention concerne un support d’enregistrement lisible par ordinateur comprenant des instructions qui, lorsqu’elles sont exécutées par un ordinateur, conduisent celui-ci à mettre en œuvre les étapes du procédé selon le premier aspect de l’invention. L’invention et ses différentes applications seront mieux comprises à la lecture de la description qui suit et à l’examen des figures qui l’accompagnent. BREVE DESCRIPTION DES FIGURES Les figures sont présentées à titre indicatif et nullement limitatif de l’invention. La figure 1 est un schéma synoptique illustrant les étapes du procédé de l’invention. La figure 2 est une représentation schématique en coupe du capteur thermique solaire et de ses échanges thermiques avec son environnement. La figure 3A est une représentation schématique vue de dessus du capteur thermique solaire. La figure 3B est une représentation schématique vue de dessus d’une surface élémentaire du capteur thermique solaire. La figure 4 est un graphe illustrant la température mesurée du fluide caloporteur en sortie du tube et la température du fluide caloporteur en sortie du tube calculée par un modèle numérique de l’art antérieur, en fonction du temps. La figure 5 est un graphe illustrant la température mesurée du fluide caloporteur en sortie du tube et la température du fluide caloporteur en sortie du tube calculée par le modèle numérique de l’invention et la température mesurée de l’enrobé et la température de l’enrobé calculée par le modèle numérique de l’invention, en fonction du temps. La figure 6 est un graphe illustrant la température mesurée du fluide caloporteur en sortie du tube et la température du fluide caloporteur en sortie du tube calculée par le modèle numérique de l’invention, en fonction du temps. DESCRIPTION DETAILLEE Sauf précision contraire, un même élément apparaissant sur des figures différentes présente une référence unique. Un premier aspect de l’invention concerne un procédé de détermination d’au moins un paramètre d’efficacité d’un capteur solaire thermique. On entend par « capteur solaire thermique », un capteur configuré pour capter l’énergie solaire transmise par rayonnement et la transférer à un fluide caloporteur sous forme d’énergie thermique. Dans la suite de la description, on utilisera indifféremment les termes « capteur solaire thermique » et « capteur ». On entend par « paramètre d’efficacité d’un capteur solaire thermique », un paramètre représentatif de la capacité du capteur solaire thermique à convertir l’énergie solaire en énergie thermique. Le paramètre d’efficacité est par exemple l’énergie E p ou la puissance thermique que peut fournir le capteur. La figure 3A montre une représentation schématique vue de dessus du capteur 200. Le capteur 200 comporte au moins un tube 201 ayant une entrée 2011 et une sortie 2012. Le tube 201 est configuré pour transporter un fluide caloporteur de son entrée 2011 vers sa sortie 2012. Le fluide caloporteur est par exemple de l’eau glycolée. Sur la figure 3A, le capteur 200 comporte un unique tube 201 mais le capteur 200 pourrait comporter une pluralité de tubes 201. Le capteur 200 comporte également un absorbeur configuré pour capter l’énergie solaire et la transmettre au fluide caloporteur circulant dans le tube 201. L’absorbeur du capteur 200 est un enrobé 202 d’une route, c’est-à-dire le revêtement de la route. L’enrobé 202 est par exemple un mélange de gravier, de sable et d’un liant hydrocarboné comme du bitume. Le capteur 200 pourrait comporter une pluralité de tubes 201 situés à des profondeurs différentes dans l’enrobé 202. L’énergie thermique peut ainsi être stockée dans la route puis restituée à des bâtiments et infrastructures environnants, via un système de pompes à chaleur, par exemple pour chauffer ces bâtiments et infrastructures. Les capteurs 200 ont pour caractéristiques principales une très forte inertie, par exemple de l’ordre de 2,4 MJ/(m 3 .K), une faible diffusivité thermique par exemple de l’ordre de 10 -6 m².s et une surface étendue, par exemple entre plusieurs centaines à plusieurs milliers de m². On entend par « inertie thermique d’un capteur », la prédisposition du capteur à garder longtemps sa température initiale lorsqu’intervient une perturbation de son équilibre thermique. Ainsi, le capteur 200 soumis à une perturbation l’amenant vers une nouvelle température d’équilibre, atteint lentement cette nouvelle température d’équilibre. On entend par « diffusivité thermique d’un capteur », la faculté d’un capteur à transporter une onde thermique entre deux points du capteur. La surface du capteur 200 peut être divisée en une pluralité de surfaces élémentaires 203 identiques. Une surface élémentaire 203 correspond par exemple à une surface de 2 m² du capteur 200. La figure 3B montre une représentation schématique vue de dessus d’une surface élémentaire 203 du capteur 200. Une surface élémentaire 203 comporte une portion 2030 du tube 201 ayant une entrée 2031 et une sortie 2032. Le tube 201 est ainsi découpé artificiellement en une pluralité de portions 2030 de tubes 201 identiques, la première portion 2030 de tube 201 ayant pour entrée 2031 l’entrée 2011 du tube 201 et la dernière portion 2030 du tube 201 ayant pour sortie 2032 la sortie 2012 du tube 201. Comme illustré sur la figure 3A, pour chaque portion 2030 de tube 201 différente de la première portion 2030 du tube 201, à savoir la portion 2030 du tube 201 ayant pour entrée 2031 l’entrée 2011 du tube 201, l’entrée 2031 de la portion 2030 du tube 201 est la sortie 2032 de la portion 2030 du tube 201 précédente. Autrement dit, pour chaque portion 2030 de tube 201 différente de la dernière portion 2030 du tube 201, à savoir la portion 2030 du tube 201 ayant pour sortie 2032 la sortie 2012 du tube 201, la sortie 2032 de la portion 2030 du tube 201 est l’entrée 2032 de la portion 2030 du tube 201 suivante. La figure 1 est un schéma synoptique illustrant les étapes du procédé de l’invention. Une première étape 101 du procédé 100 consiste à déterminer les conditions extérieures s’exerçant sur le capteur 200. Les conditions extérieures comportent au moins l’irradiance globale solaire I sol , l’irradiance thermique par rayonnement I r provenant de la voûte céleste, la vitesse du vent u et la température T a de l’air ambiant autour du capteur 200. Les conditions extérieures sont par exemple obtenues par consultation d’une base de données météorologiques permettant d’obtenir des conditions extérieures moyennes s’exerçant dans la localité dans lequel se trouve le capteur 200, puis par correction des conditions extérieures moyennes en fonction de la position du capteur 200 dans la localité. Par exemple, si le capteur solaire thermique 200 est à une altitude supérieure par rapport au thermomètre ayant permis de relever la température moyenne dans la localité, la température s’exerçant sur le capteur 200 est inférieure à la température moyenne. La correction consiste alors par exemple à diminuer la température de 1°C tous les 100 mètres d’altitude en plus. Une localité correspond par exemple à une surface de 500 m². Les conditions extérieures peuvent être une combinaison de conditions extérieures moyennes obtenues via la base de données météorologiques. Par exemple, l’irradiance thermique par rayonnement I r provenant de la voûte céleste n’est pratiquement jamais fournie par les bases de données météorologiques et est donc reconstruite, par exemple à partir de la température extérieure, de l’hygrométrie extérieure et de la quantité de précipitations Une deuxième étape 102 du procédé 100 consiste à calculer au moins la température moyenne T e de l’enrobé 202 et la température T out du fluide caloporteur en sortie 2012 du tube 201 en utilisant un modèle numérique du capteur 200 à partir des conditions extérieures déterminées à la première étape 101 et de paramètres de fonctionnement du capteur 200 fixés. Le modèle numérique du capteur 200 permet de calculer au moins la température moyenne T e de l’enrobé 202 et la température T out du fluide caloporteur en sortie 2012 du tube 201 à partir des conditions extérieures s’exerçant sur le capteur 200 et de paramètres de fonctionnement du capteur 200 fixés. Les paramètres de fonctionnement comportent par exemple au moins la température T in du fluide caloporteur à l’entrée 2011 du tube 201 et le débit massique q du fluide caloporteur dans le tube 201. Le modèle numérique est par exemple basé sur un bilan thermique du capteur 200 prenant en compte ses échanges thermiques avec son environnement. La figure 2 montre une représentation schématique en coupe du capteur 200 et de ses échanges thermiques avec son environnement. Le capteur 200, et en particulier l’enrobé 202, capte l’énergie solaire et la transforme en énergie thermique et perd notamment de l’énergie par rayonnement, par convection, notamment à cause du vent et par inertie à cause de la très forte inertie du capteur 200. Le modèle numérique prend ainsi en compte l’énergie solaire globale incidente E sol sur le capteur 200, le flux énergétique net échangé par rayonnement en grandes longueurs d’ondes E r entre la voûte céleste et le capteur 200, les pertes énergétiques par convection E conv du capteur 200, l’énergie stockée ou déstockée par inertie E in et l’énergie thermique E p produite par le capteur 200. On obtient ainsi le bilan suivant : Avec α e le coefficient d’absorption de l’enrobé 202. On a par exemple : Avec : C e la capacité thermique de l’enrobé 202 ; T e la température moyenne de l’enrobé 202 ; α e le coefficient d’absorption de l’enrobé 202. A g la surface d’échange effective du capteur 200 ; I sol l’irradiance globale solaire reçue par le capteur 200 ; e r l’émissivité globale en rayonnement de la surface du capteur 200 ; I r l’irradiance thermique provenant de la voûte céleste reçue par le capteur 200 ; σ la constante de Stefan-Bolztmann ; K(u) un coefficient d’échange convectif dépendant de la vitesse du vent u ; T a la température de l’air ambiant autour du capteur 200 ; q le débit massique de fluide caloporteur dans le tube 201 ; c pf la chaleur spécifique du fluide caloporteur ; T out la température du fluide caloporteur en sortie 2012 du tube 201 ; T in la température du fluide caloporteur en entrée 2011 du tube 201. Le modèle numérique prend également en compte les échanges thermiques entre l’enrobé 202 et le fluide caloporteur dans le tube 201. On a par exemple : Avec : q le débit massique de fluide caloporteur dans le tube 201 ; c pf la chaleur spécifique du fluide caloporteur ; T out la température du fluide caloporteur en sortie 2012 du tube 201; T in la température du fluide caloporteur à l’entrée 2011 du tube 201; h t le coefficient de transfert thermique surfacique entre le fluide caloporteur et l’enrobé 202 ; S t la surface intérieure du tube 201 ; T e la température moyenne de l’enrobé 202. Ainsi, en résolvant le système de deux équations du modèle numérique, on peut obtenir la température moyenne T e de l’enrobé 202 et la température T out du fluide caloporteur en sortie 2012 du tube 201. Le modèle numérique est par exemple composé d’une pluralité de modèles numériques élémentaires identiques correspondant chacun à une surface élémentaire 203 de la surface du capteur 200. On fait ainsi l’hypothèse que les échanges thermiques entre le tube 201 et l’enrobé 202 ne sont pas influencés par la proximité des portions 2030 de tube 201 successives. Le modèle numérique élémentaire est configuré pour calculer au moins la température moyenne T e de l’enrobé 202 de la surface élémentaire 203 et la température T out du fluide caloporteur en sortie 2032 de la portion 2030 de tube 201 comprise dans la surface élémentaire 203. Le modèle numérique élémentaire correspond au modèle numérique précédemment décrit appliqué à une surface élémentaire 203. Ainsi, il suffit d’adapter les équations précédemment décrites à une surface élémentaire 203. En reprenant l’exemple de modèle numérique précédent, on a par exemple T out la température du fluide caloporteur en sortie 2032 de la portion 2030 du tube 201 comprise dans la surface élémentaire 203 et T in la température du fluide caloporteur à l’entrée 2031 de la portion 2030 du tube 201 comprise dans la surface élémentaire 203. La température T in du fluide caloporteur en entrée 2031 de la portion 2030 du tube 201 comprise dans la surface élémentaire 203 est la température T out du fluide caloporteur en sortie 2032 de la portion 2030 de tube 201 précédant la portion 2030 du tube 201 de la surface élémentaire 203. Le procédé 100 comporte alors par exemple une étape préalable de calibration du modèle numérique élémentaire en utilisant au moins une température moyenne T e de l’enrobé 202 mesurée. Seule la température T out du fluide caloporteur en sortie 2032 de la portion 2030 de tube 201 varie pendant la calibration. Le procédé 100 peut également comporter une étape préalable de calibration du modèle numérique correspondant à la surface totale du capteur 200 en utilisant au moins une température moyenne T e de l’enrobé 202 mesurée. Seule la température T out du fluide caloporteur en sortie 2012 du tube 201 varie pendant la calibration. L’étape de calibration permet de calibrer au moins les paramètres suivants : le coefficient d’absorption α e de l’enrobé 202, le coefficient de transfert thermique surfacique h t entre le fluide caloporteur et l’enrobé 202, le coefficient d’échange convectif K(u) du capteur 200, l’émissivité globale en rayonnement e r de la surface du capteur 200, et la capacité thermique C e de l’enrobé 202. On fait l’hypothèse que les paramètres calibrés sont les mêmes pour toutes les surfaces élémentaires 203 du capteur 200. Une troisième étape 103 du procédé 100 consiste à calculer le paramètre d’efficacité à partir de la température T out du fluide caloporteur en sortie 2012 du tube 201 et/ou de la température moyenne T e de l’enrobé 202. Le paramètre d’efficacité peut être calculé sur une période de simulation, c’est-à-dire que les étapes 101 à 103 du procédé 100 peuvent être réalisées pour chaque pas de temps de la période de simulation. La période de simulation est par exemple une semaine et le pas de temps, 30 minutes. Le pas de temps est par exemple inférieur à une heure. Le paramètre d’efficacité est par exemple le rendement énergétique du capteur 200 sur la période de simulation, la puissance ou l’énergie thermique produite par le capteur 200 sur la période de simulation, la différence de températures maximale entre l’entrée 2011 du tube 201 et la sortie 2012 du tube 201 sur la période de simulation, la différence de températures moyenne entre l’entrée 2011 du tube 201 et la sortie 2012 du tube 201 sur la période de simulation, la différence de températures minimale entre l’entrée 2011 du tube 201 et la sortie 2012 du tube 201 sur la période de simulation, la température maximale de l’enrobé 202 sur la période de simulation, la température moyenne de l’enrobé 202 sur la période de simulation, la température minimale de l’enrobé 202 sur la période de simulation, la température maximale du fluide caloporteur en sortie 2012 du tube 201 sur la période de simulation, la température moyenne du fluide caloporteur en sortie 2012 du tube 201 sur la période de simulation, la température minimale du fluide caloporteur en sortie 2012 du tube 201 sur la période de simulation ou encore le nombre de jours par période de simulation pour lesquels le rendement énergétique est supérieur à un seuil de rendement. On définit le rendement énergétique η sur une période de simulation comme : Avec : E sol l’énergie solaire globale incidente sur le capteur 200 pendant la période de simulation ; E p l’énergie thermique produite par le capteur 200 pendant la période de simulation. Par exemple, pour calculer la puissance produite par le capteur sur un pas de temps donné, il suffit de diviser l’énergie thermique E p produite par le capteur 200 par le pas de temps correspondant. La figure 4 est un graphe illustrant la température Tout_ref mesurée du fluide caloporteur en sortie 2012 du tube 201 et la température Tout du fluide caloporteur en sortie 2012 du tube 201 calculée par un modèle numérique s’appuyant sur la norme internationale NF EN ISO 9806 de l’art antérieur en fonction du temps. Pour rappel, le modèle numérique de l’art antérieur ne calcule pas de température moyenne T e de l’enrobé 202 car fait l’hypothèse que celle-ci est égale à la température T out du fluide caloporteur en sortie 2012 du tube 201. La figure 5 est un graphe illustrant la température Tout_ref mesurée du fluide caloporteur en sortie 2012 du tube 201 et la température Tout du fluide caloporteur en sortie 2012 du tube 201 calculée par le modèle numérique de l’invention et la température moyenne Te_ref mesurée de l’enrobé 202 et la température moyenne Te de l’enrobé 202 calculée par le modèle numérique de l’invention en fonction du temps. La figure 6 est un graphe illustrant la température Tout_ref mesurée du fluide caloporteur en sortie 2012 du tube 201 et la température Tout du fluide caloporteur en sortie 2012 du tube 201 calculée par le modèle numérique de l’invention en fonction du temps. La figure 6 correspond à la figure 5 à l’échelle de la figure 4 sans les courbes Te et Te_ref correspondant aux températures de l’enrobé 202. Les graphes des figures 4, 5 et 6 correspondent à des mêmes conditions extérieures s’exerçant sur le capteur 200. On constate que l’écart entre la température T out_ref mesurée du fluide caloporteur en sortie 2012 du tube 201 et la température T out du fluide caloporteur en sortie 2012 du tube 201 calculée par le modèle numérique de l’art antérieur visible sur la figure 4 est supérieur à l’écart entre la température T out_ref mesurée du fluide caloporteur en sortie 2012 du tube 201 et la température T out du fluide caloporteur en sortie 2012 du tube 201 calculée par le modèle numérique de l’invention visible sur la figure 6. Le modèle de l’invention est donc plus précis que le modèle de l’art antérieur pour le calcul de la température T out du fluide caloporteur en sortie 2012 du tube 201. De plus, l’écart entre la température moyenne T e_ref mesurée de l’enrobé 202 et la température T out_ref mesurée du fluide caloporteur en sortie 2012 du tube 201 peut être supérieur à 15°C, ce qui montre que l’hypothèse utilisée par la norme est erronée pour le capteur 200. On constate que l’écart entre la température moyenne T e_ref mesurée de l’enrobé 202 et la température moyenne T e de l’enrobé 202 calculée par le modèle numérique de l’invention visible sur la figure 5 est faible. Le modèle de l’invention permet donc de calculer précisément la température moyenne T e de l’enrobé 202. Le procédé 100 est mis en œuvre par un système informatique, c’est-à-dire par un ordinateur comportant un processeur. Procédé (100) de détermination par un système informatique d’au moins un paramètre d’efficacité d’un capteur (200) solaire thermique comportant au moins un tube (201) véhiculant un fluide caloporteur d’une entrée (2011) du tube (201) vers une sortie (2012) du tube (201), le tube (201) étant situé dans un absorbeur (202) configuré pour capter l’énergie solaire et transmettre l’énergie solaire captée au fluide caloporteur sous forme de chaleur, le paramètre d’efficacité étant représentatif de la capacité du capteur (200) à convertir l’énergie solaire en énergie thermique, le procédé (100) étant caractérisé en ce que l’absorbeur (202) est un enrobé (202) d’une route et comporte les étapes suivantes : déterminer des conditions extérieures s’exerçant sur le capteur (200, 101) ; calculer au moins une température (T out ) du fluide caloporteur en sortie (2012) du tube (201) et une température moyenne (T e ) de l’enrobé (202) à l’aide d’un modèle numérique du capteur (200) configuré pour calculer au moins la température moyenne (T e ) de l’enrobé (202) et la température (T out ) du fluide caloporteur en sortie (2012) du tube (201) à partir des conditions extérieures déterminées et de paramètres de fonctionnement du capteur (200) fixés (102) ; calculer le paramètre d’efficacité à partir de la température (T out ) du fluide caloporteur en sortie (2012) du tube (201, 103) et/ou de la température moyenne (T e ) de l’enrobé (202). Procédé (100) selon la revendication 1, caractérisé en ce que le paramètre d’efficacité est calculé sur une période de simulation, les étapes (101, 102, 103) du procédé (100) étant réalisées pour chaque pas de temps de la période de simulation. Procédé (100) selon la revendication 2, caractérisé en ce que le pas de temps est inférieur à une heure. Procédé (100) selon l’une quelconque des revendications 2 ou 3, caractérisé en ce que le paramètre d’efficacité peut être choisi parmi le groupe formé des paramètres suivants : rendement énergétique du capteur (200) sur la période de simulation, puissance ou énergie thermique produite par le capteur (200) sur la période de simulation, différence de températures maximale entre l’entrée (2011) du tube (201) et la sortie (2012) du tube (201) sur la période de simulation, différence de températures moyenne entre l’entrée (2011) du tube (201) et la sortie (2012) du tube (201) sur la période de simulation, différence de températures minimale entre l’entrée (2011) du tube (201) et la sortie (2012) du tube (201) sur la période de simulation, température maximale de l’enrobé (202) sur la période de simulation, température moyenne de l’enrobé (202) sur la période de simulation, température minimale de l’enrobé (202) sur la période de simulation, température maximale du fluide caloporteur en sortie (2012) du tube (201) sur la période de simulation, température moyenne du fluide caloporteur en sortie (2012) du tube (201) sur la période de simulation, température minimale du fluide caloporteur en sortie (2012) du tube (201) sur la période de simulation ou nombre de jours par période de simulation pour lesquels le rendement énergétique est supérieur à un seuil de rendement. Procédé (100) selon l’une quelconque des revendications précédentes, caractérisé en ce que les conditions extérieures comportent l’irradiance solaire, l’irradiance thermique par rayonnement provenant de la voûte céleste, la vitesse du vent et la température de l’air ambiant autour du capteur (200). Procédé (100) selon l’une quelconque des revendications précédentes, caractérisé en ce que les paramètres de fonctionnement comportent la température (T in ) du fluide caloporteur à l’entrée (2011) du tube (201) et le débit massique du fluide caloporteur dans le tube (201). Procédé (100) selon la revendication 6, caractérisé en ce que la surface du capteur (200) est divisée en une pluralité de surfaces élémentaires (203) identiques comportant chacune une portion (2030) du tube (201) et le modèle numérique est composé d’une pluralité de modèles numériques élémentaires identiques correspondant chacun à une surface élémentaire (203), le modèle numérique élémentaire étant configuré pour calculer la température moyenne (T e ) de l’enrobé (202) de la surface élémentaire (203) et la température (T out ) du fluide caloporteur en sortie (2032) de la portion (2030) de tube (201) comprise dans la surface élémentaire (203) en utilisant la température (T out ) du fluide caloporteur en sortie (2032) de la portion (2030) de tube (201) précédente calculée par le modèle élémentaire précédent pour température (T in ) du fluide caloporteur en entrée (2031) de la portion (2030) du tube (201) comprise dans la surface élémentaire (203). Procédé (100) selon la revendication 7, caractérisé en ce qu’il comporte une étape préalable de calibration du modèle numérique élémentaire en utilisant au moins une température moyenne (T e ) de l’enrobé (202) mesurée ou une étape préalable de calibration du modèle numérique en utilisant au moins une température moyenne (T e ) de l’enrobé (202). Procédé (100) selon la revendication 8, caractérisé en ce que l’étape de calibration permet de calibrer les paramètres suivants : le coefficient d’absorption de l’enrobé (202), le coefficient de transfert thermique surfacique entre le fluide caloporteur et l’enrobé (202), le coefficient d’échange convectif du capteur (200), l’émissivité globale en rayonnement de la surface du capteur (200) et la capacité thermique de l’enrobé (202). Procédé (100) selon l’une quelconque des revendications précédentes, caractérisé en ce que l’étape (101) de détermination des conditions extérieures s’exerçant sur le capteur (200) comporte une sous-étape de récupération de conditions extérieures moyennes s’exerçant dans une localité dans laquelle est situé le capteur (200), à partir d’une base de données météorologiques et une sous-étape de correction des conditions extérieures moyennes en fonction de la localisation du capteur (200) dans la localité. Procédé (100) selon l’une quelconque des revendications précédentes, caractérisé en ce que le modèle numérique du capteur (200) comporte les équations suivantes : Avec : E sol l’énergie solaire globale incidente sur le capteur (200) ; E r le flux énergétique net échangé par rayonnement en grandes longueurs d’ondes entre la voûte céleste et le capteur (200) et son environnement ; E conv les pertes énergétiques par convection E conv du capteur (200) ; E in l’énergie stockée ou déstockée par inertie ; E p l’énergie thermique produite par le capteur (200) ; C e la capacité thermique de l’enrobé (202) ; T e la température moyenne de l’enrobé (202) ; α e le coefficient d’absorption de l’enrobé (202) ; A g la surface d’échange effective du capteur (200) ; I sol l’irradiance globale solaire reçue par le capteur (200) ; e r l’émissivité globale en rayonnement de la surface du capteur (200) ; I r l’irradiance thermique provenant de la voûte céleste reçue par le capteur (200) ; σ la constante de Stefan-Bolztmann ; K(u) un coefficient d’échange convectif dépendant de la vitesse du vent u ; T a la température de l’air ambiant autour du capteur (200) ; q le débit massique de fluide caloporteur dans le tube (201) ; c pf la chaleur spécifique du fluide caloporteur ; T out la température du fluide caloporteur en sortie (2012) du tube (201) ; T in la température du fluide caloporteur à l’entrée (2011) du tube (201) ; h t le coefficient de transfert thermique surfacique entre le fluide caloporteur et l’enrobé (202) ; S t la surface intérieure du tube (201). Produit-programme d’ordinateur comprenant des instructions qui, lorsque le programme est exécuté par un ordinateur, conduisent celui-ci à mettre en œuvre les étapes du procédé selon l’une quelconque des revendications 1 à 11.