FR 3036159 A1 20161118 FR 1554207 A 20150512 1 La présente invention concerne un procédé de remplissage ou de soutirage d'un réservoir de gaz sous pression. L'invention concerne plus un procédé de remplissage ou de soutirage d'un réservoir de gaz sous pression, notamment de carburant tel que de l'hydrogène gazeux, le réservoir étant délimité par une paroi de forme générale cylindrique ayant des dimensions et des propriétés thermo-physiques déterminées et connues, le procédé comprenant une régulation du débit du flux de gaz introduit ou, respectivement du gaz soutiré, et/ou de la température dudit gaz introduit, pour éviter que le réservoir n'atteigne un seuil haut déterminé de température ou un seuil bas déterminé de température. L'invention concerne en particulier le remplissage ou le soutirage de réservoirs d'hydrogène sous pression en particulier des réservoirs composites stockant du gaz à des pressions élevée, notamment 200bar, 700bar ou 1000bar. Ces réservoirs composites (notamment de type IV) ont une résistance aux températures extrêmes plus limitée que les réservoirs métalliques. Par exemple, les températures limites d'utilisation de tels réservoirs doivent être comprises entre -40°C et +85°C. Ainsi, les processus de remplissage (produisant un échauffement) et le processus de soutirage (produisant un refroidissement) doivent être contrôlés. Ce contrôle peut être réalisé notamment en contrôlant le débit de gaz et/ou sa température. Un standard (SAE) de remplissage de véhicules automobiles a été développé pour établir des recommandations de refroidissement et des vitesses de remplissage en fonction de conditions initiales (température ambiante, pression initiale dans le réservoir...). Ce modèle se base sur la température du gaz dans le réservoir et non pas la température du réservoir (qui est généralement inférieure à celle du gaz). La mesure de la température maximale atteinte par le réservoir est cependant difficilement mesurable. Dans ces conditions, les recommandations des conditions de remplissage établies par ce standard sont considérées trop contraignantes et donc en partie inutiles et générant des surcoûts importants. 3036159 2 Un but de la présente invention est de pallier tout ou partie des inconvénients de l'art antérieur relevés ci-dessus. A cette fin, le procédé selon l'invention, par ailleurs conforme à la définition générique qu'en donne le préambule ci-dessus, est essentiellement caractérisé en 5 ce qu'il comprend une étape d'estimation en temps réel par calcul d'au moins une température du réservoir parmi : la température moyenne de la paroi du réservoir Twall,average(r,t) en fonction du temps (t), la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps, la température minimale atteinte par la paroi du réservoir Twall,min(t) en fonction du temps et en ce que, la 10 régulation du débit du flux de gaz ou de la température dudit gaz est réalisée en fonction de ladite température du réservoir calculée. Selon l'invention il est ainsi possible d'estimer, via un modèle mathématique basé sur des mesures de pression dans le réservoir et de la température ambiante, la température maximale ou minimale atteint par la paroi du réservoir. 15 Les valeurs de pression du gaz et de température ambiante sont des données externes mesurées. C'est-à-dire qu'il n'est pas nécessaire de mesures complexes dans le réservoir rempli ou soutiré. Un modèle mathématique dynamique à trois dimensions peut être utilisé pour calculer la température dans le réservoir, y compris au niveau de sa paroi. 20 Cependant, de préférence l'invention utilise un modèle simplifié. Ce modèle simplifié peut être basé notamment sur les données d'entrée que sont : la température du gaz dans la ligne de transfert vers le réservoir, sa pression, les caractéristiques géométriques du réservoir, la température ambiante. D'autres paramètres thermo-physiques du réservoir peuvent être utilisés 25 comme par exemple la conductivité thermique, la capacité calorifique et la masse volumique. Comme illustré à la figure 3, la masse initiale m(0) de gaz dans le réservoir 1 peut également une donnée d'entrée connue ou estimée. De même pour le débit Q(t) de gaz transféré au temps t. 30 De préférence, le modèle calcule sur cette base et en temps réel la température moyenne du gaz dans le réservoir et son état de remplissage (« SOC= pourcentage de remplissage dans lequel 100% correspond à la quantité massique de gaz du réservoir rempli à sa pression de service et une température 3036159 3 homogène déterminée, notamment 15°C). Cet état de charge SOC est directement lié à la température moyenne du gaz dans le réservoir et à sa pression via une équation d'état du gaz. De préférence, le modèle calcule sur cette base en temps réel également la 5 température moyenne de la paroi du réservoir selon deux dimensions Tw,average(r,t) dans laquelle r est la coordonnée spatiale (notamment le rayon r par rapport à l'axe central longitudinal du réservoir cylindrique) et t le temps. De préférence, le modèle calcule sur cette base en temps réel également la température maximale de la paroi du réservoir Tw,max(t) dans laquelle t est le 10 temps. Le modèle peut calculer également la température minimale de la paroi du réservoir Tw,min(t) en fonction du temps t. Ces données peuvent être obtenues le cas échéant à partir de la température Tgas,average(t) du gaz (calculée). Ce modèle mathématique et thermodynamique peut être intégré à un 15 contrôleur logique programmable ("PLC") permettant d'estimer en temps réel la température extrême durant un remplissage ou un soutirage. La température moyenne du gaz dans le réservoir peut également être calculée. A partir de cette donnée de la température extrême de la paroi, le contrôleur peut ajuster la vitesse de transfert de gaz (remplissage ou soutirage) 20 et/ou le niveau de refroidissement du gaz (remplissage) pour rester dans les températures de sécurité. Cette solution permet d'éviter des régulations de débit ou des refroidissements surdimensionnés par rapport au besoin réel. Selon d'autres particularités possibles : 25 - lors d'un remplissage, lorsque la température du réservoir calculée atteint un seuil haut (SH) déterminé, le flux et notamment le débit de gaz est diminué et/ou la température du gaz fourni au réservoir est diminuée par échange thermique avec une source de froid, - lors d'un remplissage, lorsque la température du réservoir calculée est 30 inférieure au seuil haut d'une valeur déterminée, le flux de gaz et notamment son débit est augmenté et/ou la température dudit gaz et/ou du réservoir est augmentée par échange thermique avec une source de chaleur, 3036159 4 - le procédé comprend une étape de calcul du Nombre de Richardson (Ri) pour le gaz dans le réservoir en fonction du temps, une étape de comparaison du nombre de Nombre de Richardson calculé avec une valeur de référence (Vr) déterminée comprise entre 0,05 et 1,5 et de préférence entre 0,05 et 0,15 et, 5 lorsque nombre de Nombre de Richardson calculé est inférieur à valeur de référence déterminée le température du gaz dans le réservoir est considérée homogène c'est-à-dire que la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps est égale à la température moyenne de la paroi du réservoir Twall,average(r,t) en contact avec le gaz en fonction du 10 temps (t) : Twall,max(t)= Twall,average(r=r_liner,t). - lorsque le nombre de Nombre de Richardson calculé est supérieure à valeur de référence déterminée, la température du gaz dans le réservoir est considérée hétérogène c'est-à-dire que la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps n'est pas égale à la 15 température moyenne de la paroi du réservoir en contact avec le gaz en fonction du temps (t) c'est-à-dire au niveau de ce qui forme le « liner » du réservoir, Twall,average(r= rayon du liner du réservoir, t), et dans ces conditions, le procédé comporte une étape d'augmentation du débit de gaz fourni au réservoir pour diminuer la valeur du Nombre de Richardson calculé au-dessous de la valeur de 20 référence (Vr) déterminée et rendre ainsi le gaz homogène en température, - lors d'un remplissage, lorsque le Nombre de Richardson est supérieur à la valeur de référence, le flux de gaz et notamment son débit est augmenté, - le procédé comprend, avant remplissage ou, respectivement, avant soutirage, une étape de détermination ou de détection par capteur(s), de la 25 température initiale T(0) du gaz dans le réservoir, de la pression initiale P(0) du gaz dans le réservoir, de la température moyenne initiale de la paroi du réservoir Tw, average(0) et une étape de détermination de la masse de gaz initiale dans le réservoir m(0) puis, lors d'un remplissage lorsque le nombre de Nombre de Richardson calculé est inférieure à valeur de référence déterminée le température 30 du gaz dans le réservoir est considérée homogène, c'est-à-dire que la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps est égale à la température moyenne de la paroi du réservoir Twall,average(r= r=rayon du liner du réservoir,t) en contact avec le gaz en fonction 3036159 5 du temps (t), et dans ces conditions, le procédé comporte, au cours du remplissage, une étape de calcul de la température Tgas,average(t) moyenne du gaz dans le réservoir en temps réel en fonction du temps et de la température moyenne de la paroi Twall,average(r,t) en temps réel en fonction du temps (t) à 5 partir d'un bilan de masse et d'enthalpie appliqué au gaz dans le réservoir et à partir également d'un bilan d'énergie dans la paroi du réservoir, de l'équation d'état du gaz (notamment équation des gaz réels), et d'un bilan des échanges thermiques entre le gaz et la paroi, et entre la paroi du réservoir est l'extérieur, - le procédé comprend, avant remplissage ou, respectivement, avant 10 soutirage, une étape de détermination ou de détection par capteur(s), de la température initiale T(0) du gaz dans le réservoir, de la pression initiale P(0) du gaz dans le réservoir, de la température moyenne initiale de la paroi du réservoir Tw, average(0) et une étape de détermination de la masse de gaz initiale dans le réservoir m(0) puis, lors d'un remplissage lorsque le nombre de Nombre de 15 Richardson calculé est supérieure à valeur de référence déterminée la température du gaz dans le réservoir est considérée hétérogène c'est-à-dire que la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps n'est pas égale à la température moyenne de la paroi du réservoir Twall,average(r)=r_liner,t) en contact avec le gaz en fonction du temps (t), et dans 20 ces conditions, le procédé comporte, au cours du remplissage, une étape de calcul de la température Tgas,average(t) moyenne du gaz dans le réservoir en temps réel en fonction du temps (t) et de la température moyenne de la paroi Twall,average(r,t) en temps réel en fonction du temps à partir d'un bilan de masse et d'enthalpie appliqué au gaz dans le réservoir et d'un bilan d'énergie dans la 25 paroi du réservoir, de l'équation d'état du gaz (équation des gaz parfaits ou réels), et d'un bilan des échanges thermiques entre la paroi du réservoir est l'extérieur et en ce que le procédé comprend une étape de calcul de la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps, cette température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction 30 du temps étant obtenue par corrélation à partir de la température Tgas,average(t) moyenne du gaz dans le réservoir calculée en temps réel en fonction du temps et en fonction de la température moyenne de la paroi Twall,average(r,t) en temps réel en fonction du temps (t), 3036159 6 - le procédé comprend une étape de calcul de l'enthalpie hin(t) du gaz entrant ou sortant du réservoir en fonction du temps, une étape de mesure ou de calcul de la masse de gaz m(t) introduite ou soutiré du réservoir en fonction du temps ou, respectivement, une étape de détermination de la pression P(t) dans le 5 réservoir en fonction du temps, le procédé comprenant une étape de détermination de la température moyenne du gaz Tgas,average(t) à l'instant t dans le réservoir en degré K, cette température moyenne Tgas,average(t) étant exprimée comme une fonction du premier degré de la température moyenne du gaz T(t-1) à l'instant précédent (t-1) et d'un coefficient d'échange chaleur convectif 10 entre le gaz et la paroi interne du réservoir (1) à l'instant (t-1) en W.m-2.K1, dans lequel le coefficient d'échange chaleur kg(t-1) est donné par la relation kg=(Àg/Dint).Nuint dans laquelle Àg est la conductivité thermique du gaz dans le réservoir en W.m-1.K1, Dint est le diamètre interne du réservoir (1) en mètre et NuDint le nombre de Nusselt du gaz dans le réservoir (1) (sans dimension), et 15 dans lequel le nombre de Nusselt du gaz est exprimé en fonction du nombre de Reynolds (Redin) (sans dimension) relatif à la convection forcée dans le réservoir et du nombre de Rayleigh (RaDint) (sans dimension) relatif à la convection naturelle interne dans le réservoir selon une formule NuDint= a.RaDintb + c.Redind dans laquelle a et c sont des coefficients sans dimension fonction du rapport 20 (Lint/Dint) entre la longueur interne du réservoir Lint en mètre et le diamètre interne du réservoir Dint en mètre et du rapport (Dint/di) entre le diamètre interne du réservoir Dint en mètre et le diamètre de l'injecteur di en mètre, a, b, c et d étant des nombres réels positifs sans dimension, a étant compris entre 0 et 1, b étant compris entre 0,2 et 0,5, c étant compris entre 0 et 1 et d étant compris entre 25 0.5 et 0.9, - le procédé comprend une étape d'estimation de la température maximale atteinte dans l'épaisseur du réservoir Twall,max(t) en fonction du temps, ou respectivement la température minimale atteinte dans l'épaisseur du réservoir Twall,min(t) en fonction du temps, à partir de la température Tgas,average(t) 30 moyenne du gaz dans le réservoir calculée en temps réel en fonction du temps (t), ladite température maximale ou minimale étant obtenue à partir de tableau(x) de correspondance obtenu(s) par des expérimentations, et/ou à partir de calcul et/ou de simulations de sorte qu'à une température Tgas,average(t) moyenne du gaz 3036159 7 dans le réservoir calculée en temps réel en fonction du temps (t) correspond, selon les conditions connues de débit de gaz, du diamètre d'injecteur de gaz dans le réservoir et des dimensions et caractéristiques du réservoir, une température maximale ou, respectivement minimale, dans l'épaisseur du réservoir en fonction 5 du temps, - le Nombre de Richardson (Ri) pour le gaz dans le réservoir en fonction du temps est calculé à partir du nombre de Grashoff (Gr) pour le gaz et du nombre de Reynolds (Re) pour le gaz au temps (t) selon la formule suivante : Ri= Gr/Re2, dans lequel les nombres de Grashof (Gr) et de Reynolds (Re) sont des données 10 connues ou calculées à partir de la valeur de pression du gaz mesurée ou de la masse de gaz dans le réservoir et de la température du gaz Tgas,average(t) moyenne du gaz dans le réservoir, - lors d'un soutirage, lorsque la température du réservoir calculée atteint un seuil bas déterminé, le flux de gaz est diminué, 15 - le procédé comprend une étape de mesure ou de calcul en continu de l'un au moins parmi : la pression Pin(t) au temps (t) du gaz introduit dans ou soutiré du réservoir, la température Tin(t) au temps (t) du gaz introduit dans ou soutiré du réservoir, la P(t) pression dans le réservoir au temps (t), la température ambiante (Tamb(t) au temps (t), la masse de gaz m(t) dans le réservoir au temps (t), 20 - le procédé est mis en oeuvre par une station de remplissage de réservoirs de gaz hydrogène comprenant au moins une source d'hydrogène à haute pression, au moins une conduite de transfert de reliant sélectivement la source à un réservoir, et un organe électronique d'acquisition, de stockage et de traitement de données tel qu'un ordinateur ou microprocesseur, ledit organe pilotant le 25 transfert de gaz entre la source et le réservoir, la logique (4) électronique étant programmée pour calculer en temps réel au moins une température du réservoir parmi : la température moyenne de la paroi du réservoir Twall,average(r,t) en fonction du temps (t), la température maximale atteinte dans l'épaisseur du réservoir Twall,max(t) en fonction du temps, la température minimale atteinte dans 30 l'épaisseur du réservoir Twall,min(t) en fonction du temps et pour réguler le débit du flux de gaz et/ou de la température dudit gaz en fonction de ladite température du réservoir calculée, 3036159 8 - lors d'un soutirage, lorsque la température du réservoir calculée est supérieure au seuil bas (SB), le flux de gaz est augmenté, - la température initiale T(0) du gaz dans le réservoir, et la température moyenne initiale de la paroi du réservoir T,(0) sont approximée comme étant 5 égales à la température ambiante mesurée ou renseignée avant remplissage, - les dimensions géométriques du réservoir sont connues, par exemple communiquées avant remplissage notamment : la longueur interne (Lint) du réservoir, le diamètre interne (Dint) du réservoir, - le diamètre (di) de l'injecteur en mètre est également un paramètre connu 10 avant remplissage, - la masse de gaz initiale m(0) dans le réservoir (en kg) est calculée par une équation d'état du gaz à partir des valeurs connues de pression initiale P(0) (en Pa) et de température initiale T(0) du gaz (en K), - en cas de soutirage le nombre NuDint de Nusselt (sans dimension) du gaz 15 dans le réservoir (1) est donné par la formule NuDint= a.RaDintb' - les propriétés thermo-physiques déterminées et connues du réservoir comprennent au moins l'un parmi : la masse volumique, la capacité calorifique spécifique, la conductivité thermique de l'enveloppe intérieure du réservoir (liner), la masse volumique, la capacité calorifique spécifique, la conductivité thermique 20 du composite, le volume du réservoir, les longueurs et diamètres internes et externes du réservoirs, l'épaisseur du liner et du composite, la surface interne et externe du réservoir, la masse totale du liner, la masse totale du composite, le diamètre d'injecteur dans le réservoir, le coefficient de perte de charge à l'entrée de la bouteille, l'émissivité de la surface externe du réservoir, 25 - la température Tgas,average(t) moyenne du gaz dans le réservoir est la température moyennée pour le volume de stockage du réservoir (en K), - le réservoir est de forme générale cylindrique, la température moyenne de la paroi du réservoir Twall,average(r, t) en fonction d'une variable d'espace (r) et du temps (t), est la température dans l'épaisseur de la paroi du réservoir selon une 30 direction transverse à la paroi, la variable d'espace (r) étant la position dans l'épaisseur de la paroi selon une direction parallèle à un rayon du cylindre à partir de l'axe longitudinal central du cylindre (en m), 3036 159 9 - la température maximale atteinte dans l'épaisseur du réservoir Twall,max(t) en fonction du temps est la température maximale (en K) atteinte par la paroi au temps (t), quel que soit la portion de la paroi, c'est-à-dire une température maximale au sens spatial du terme, 5 - la température initiale T(0) du gaz dans le réservoir est mesurée par un capteur de température, par exemple un capteur de température ambiante, - la pression initiale P(0) du gaz dans le réservoir (en Pa) est mesurée par un capteur de pression à l'entrée ou dans le réservoir, - la température moyenne initiale de la paroi du réservoir T,(0) (en K) est 10 mesurée par un capteur de température ou est approximée à la température initiale du gaz dans le réservoir T(0), - la détermination de la masse de gaz initiale dans le réservoir m(0) est calculée via l'équation des gaz parfaits ou réels à partir des données suivantes connues : la pression initiale P(0) du gaz dans le réservoir, le volume V du 15 réservoir, la température initiale du gaz dans le réservoir T(0) et du facteur de compressibilité Z du gaz réel (sans dimension) ou de la constante R des gaz parfaits, - le seuil haut (SH) de température est compris entre 65 et 105°C et de préférence égal à 85°C, 20 - le flux de gaz est diminué lors d'un remplissage en réduisant la vitesse d'augmentation de pression dans le réservoir, - le seuil bas de température est compris entre -50°C et -30°C et de préférence égal à -40°C, - la température et la pression du gaz dans le réservoir sont considérées 25 homogènes dans les estimations ou calculs, - la vitesse de remplissage est contrôlée au moyen d'une vanne ou un régulateur de pression de type proportionnel intégral, - le réservoir est un réservoir composite de type IV ou III, - la paroi du réservoir composite comprend un liner plastique ou métallique 30 et une couche de matériau composite, - la température moyenne du gaz Tgas,average(t) est déterminée en temps réel par la résolution numérique de deux bilans enthalpiques : un premier bilan enthalpique appliqué au gaz injecté dans le réservoir et un second bilan 3036159 10 enthalpique appliqué à la paroi du réservoir; de la résolution numérique de l'équation de la chaleur dans la paroi ;en utilisant l'équation d'état des gaz réels appliquée au gaz dans le réservoir, les propriétés thermophysiques du gaz tel que le facteur de compressibilité Z, la capacité thermique massique cp et l'enthalpie 5 massique du gaz étant connues, - la température moyenne (moyenne en deux dimension 2D) de la paroi du réservoir Tw(r,t) est déterminée en temps réel par la résolution numérique de deux bilans enthalpiques : un premier bilan enthalpique appliqué au gaz injecté dans le réservoir et un second bilan enthalpique appliqué à la paroi du réservoir; de la 10 résolution numérique de l'équation de la chaleur dans la paroi ; en utilisant l'équation d'état des gaz réels appliquée au gaz dans le réservoir, les propriétés thermophysiques du gaz tel que le facteur de compressibilité Z, la capacité thermique massique Cp et l'enthalpie massique h du gaz étant connues, - le nombre de Richardson (sans dimension) est calculé comme le rapport 15 entre le nombre de Grashoff (sans dimension) et le nombre de Reynolds (sans dimension) au carré, - le nombre de Reynolds peut être égal au nombre de Reynolds du gaz en entrée de réservoir (masse volumique en entrée que multiplie la vitesse en entrée que multiplie le diamètre de 1"injecteur de gaz et divisé par la viscosité dynamique 20 à l'entrée), - le nombre de Reynolds peut être égal au nombre de Reynolds au fond de la bouteille (masse volumique du gaz dans la bouteille que multiplie la vitesse du gaz au fond de la bouteille que multiplie le diamètre de la bouteille et divisé par la viscosité dynamique du gaz, 25 - la vitesse au fond de la bouteille (vf) (en m.s-1) peut être calculée à partir de la vitesse du gaz en entrée de bouteille (ve), du diamètre de l'injecteur (di) (en mètre), de la longueur interne de la bouteille (Lint) (en mètre) via la formule suivante : vf = (ve.6,2.di)/Lint - le nombre de Grashof Gr peut être calculé par la formule : constante 30 d'accélération du gaz * coefficient de dilatation thermique isobare * (T (t) - Tw(r = rliner, t))*(Dint)3/(viscosité cinématique du gaz)2, 3036159 11 - le réservoir a un rapport (Lint/Dint) entre la longueur interne du réservoir Lint en mètre et le diamètre interne du réservoir Dint en mètre compris entre 1 et 7 et de préférence compris entre 1,8 et 6,6, et de préférence inférieur à 4,5, - le rapport (Dint/di) entre le diamètre interne du réservoir Dint en mètre et 5 le diamètre di de l'injecteur en mètre est compris entre 30 et 100 et de préférence compris entre 35,0 et 72,3, - le nombre de Reynolds (Red) (sans dimension) relatif à la convection forcée dans le réservoir est compris entre 5,6.104 et 2,0.106 - le nombre de Rayleigh (Radint) (sans dimension) relatif à la convection 10 naturelle interne dans le réservoir (1) est compris entre 8,0.108 et 1,0.1012 - la masse mt(t) de gaz dans le réservoir est calculée en temps réel via l'équation d'état du gaz (équation des gaz réels de préférence) à partir de la valeur de température moyenne du gaz Tgas,average(t) dans le réservoir calculée (1) et de la pression P(t) du gaz dans le réservoir mesurée en temps réel 15 L'invention peut concerner également un dispositif de remplissage ou de soutirage de gaz d'un réservoir comprenant une conduite de transfert comprenant une vanne et raccordable à un réservoir, le dispositif comprenant un organe électronique d'acquisition, de stockage et de traitement de données tel qu'un 20 ordinateur ou microprocesseur, ledit organe pilotant le transfert de gaz entre la source et le réservoir, la logique électronique étant programmée pour calculer au moins une température du réservoir parmi : la température moyenne de la paroi du réservoir Twall,average(r,t) selon son épaisseur (rayon r) et en fonction du temps (t), la température maximale atteinte dans l'épaisseur du réservoir 25 Twall,max(t) en fonction du temps, la température minimale atteinte dans l'épaisseur du réservoir Twall,min(t) en fonction du temps et pour réguler le du débit du flux de gaz en fonction de ladite température du réservoir courante calculée.. L'invention peut concerner également tout dispositif ou procédé alternatif 30 comprenant toute combinaison des caractéristiques ci-dessus ou ci-dessous. D'autres particularités et avantages apparaîtront à la lecture de la description ci-après, faite en référence aux figures dans lesquelles : 3036159 12 - la figure 1 représente une vue schématique et partielle illustrant un exemple d'installation de remplissage de réservoir pouvant mettre en ouvre l'invention, - les figures 2 à 6 illustrent de façon schématique et partielle différents 5 logigrammes illustrant des exemples d'étapes pouvant être mises en oeuvre par le procédé selon l'invention. La figure 1 illustre schématiquement une station de remplissage (ou de soutirage) d'un réservoir 1 de gaz sous pression, notamment de carburant tel que de l'hydrogène gazeux. 10 Le réservoir 1, par exemple un réservoir composite de type IV, est délimité par une paroi 1 de forme générale cylindrique ayant des dimensions et des propriétés thermo-physiques déterminées et connues. La station peut comprendre au moins une source 10 d'hydrogène à haute pression, au moins une conduite 2 de transfert reliant sélectivement la source 1 au 15 réservoir 1 et un organe 4 électronique d'acquisition, de stockage et de traitement de données tel qu'un ordinateur ou microprocesseur. L'organe 4 électronique pilote le transfert de gaz entre la source 10 et le réservoir 1 et peut être programmée pour calculer en temps réel au moins une température du réservoir (1) parmi : la température moyenne de la paroi du 20 réservoir Twall,average(r,t) en fonction du temps (t), la température maximale atteinte dans l'épaisseur du réservoir Twall,max(t) en fonction du temps, la température minimale atteinte dans l'épaisseur du réservoir Twall,min(t) en fonction du temps et pour réguler le débit du flux de gaz et/ou de la température dudit gaz en fonction de ladite température du réservoir calculée (cf. figure 3). 25 Bien entendu, classiquement, en plus du contrôle de la température extrême atteinte par la paroi du réservoir (minimale et/ou maximale), la masse de gaz dans le réservoir est de préférence contrôlée également (ou tout autre paramètre reflétant la quantité de gaz dans le réservoir). Cette masse de gaz peut être calculée classiquement à partir de la température du gaz calculée et de la 30 pression du gaz mesurée. Les paramètres d'entrée connus pour ce ou ces calculs d'estimations comprennent par exemple : 3036159 13 - les propriétés thermodynamiques du gaz (du réservoir notamment de son liner et de structure composite), la géométrie GE du réservoir (longueur, diamètre, ... ). Ces données sont des constantes connues. - les conditions connues de Pression P(0), température T(0) et 5 température de la paroi Twall(0) à l'instant initial t=0. Ces conditions peuvent être mesurées ou déterminées ou approximées, - les conditions en temps réel de Pression du gaz entrant au temps t Pin(t), de température du gaz entrant Tin(t) au temps t, la pression du gaz dans le réservoir P(t) au temps t (mesurée par exemple dans la conduite reliée à 10 l'entrée/sortie du réservoir 1), la température ambiante Tamb(t) au temps t. - des coefficients a, b, c, d d'échange de chaleur et de corrélation (explicités ci-après). Les coefficients a, b, c, d peuvent être obtenus par des essais expérimentaux pour chaque type de réservoir à partir de tests de montée et de 15 chute de pression dans le réservoir. Ces coefficients peuvent le cas échéant être corrélés à des dimensions ou rapports de dimensions du réservoir. A partir de ces données d'entrée connues, l'organe 4 électronique peut être configuré pour calculer en temps réel les données de sortie suivantes : -la température moyenne du gaz dans le réservoir T(t)= Tgas,average(t) 20 au temps t en fonction du temps t, - la masse m(t) de gaz dans le réservoir au temps, - la température moyenne de la paroi du réservoir Twall,average(r,t) en fonction du temps t, r étant la coordonnée de rayon à partir de l'axe longitudinal du réservoir, 25 - la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps, - le nombre de Richardson Ri(t) pour le gaz dans le réservoir au temps t (cf. ci-après). La température moyenne de la paroi du réservoir Twall,average(r,t) 30 représente la moyenne selon deux dimensions (2D), c'est-à-dire qu'elle représente la température de la couche de paroi à la coordonnée r prise à partir de l'axe longitudinal du réservoir. Cette température est homogène selon deux 3036159 14 dimensions mais peut varier selon le rayon r. Cette température moyenne est calculée en résolvant par exemple l'équation de la chaleur dans la paroi. La température maximale (respectivement minimale) atteinte par la paroi du réservoir Twall,max(t) en fonction du temps peut être la température de la paroi 5 au temps t au niveau de l'interface entre le gaz et la paroi. Ci-après sera décrit un exemple d'utilisation d'un tel modèle pour une station de remplissage. Supposons un remplissage réalisé à une température de gaz d'entrée constante. A chaque intervalle de temps le modèle (mis en oeuvre par l'organe 4 10 électronique qui pilote le remplissage/soutirage) estime en temps réel la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps. Si cette température maximale devient proche de la limite admissible (85°C par exemple), dans ce cas l'organe 4 de contrôle peut réduire la montée en 15 pression en agissant par exemple sur la vanne 3 de contrôle. Cette diminution de la rampe de pression (augmentation de pression par unité de temps) réduit ou supprime l'augmentation de la température. Si la température maximale est en-dessous de la valeur admissible, la vitesse d'augmentation de pression peut être augmentée. 20 Il s'agit donc d'un procédé de contrôle du débit de remplissage en fonction de la température maximale de la paroi du réservoir estimée/calculée. Un exemple d'application est illustré schématiquement à la figure 5. Ainsi, à partir des paramètres d'entrée (étapes 11, 12, 13, 14,15) les températures moyenne du gaz Tgas,average(t) en fonction du temps t et la température 25 maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps sont calculées (étape 16). Sur cette base, l'état de charge SOC est calculé (étape 17) puis cet état de charge est comparé avec l'état de charge cible (100%). Si l'état de charge cible est atteint le remplissage est interrompu (étape 19) sinon il est poursuivi et la température maximale du réservoir Twall,max(t) est comparée au 30 seuil limite 85°C (moins un facteur de sécurité TS) cf. étape 20. Si cette température reste dans les limites admissibles (étape 21) la vitesse de remplissage peut être maintenue ou augmentée et le processus retourne à l'étape de mesure de la pression de (l'étape 15). Si cette température n'est pas dans les 3036159 15 limites admissibles (étape 22) la vitesse de remplissage peut être réduite et le processus retourne à l'étape de mesure de la pression P(t) de (l'étape 15). Alternativement ou cumulativement au contrôle du débit, la température du gaz peut être contrôlée (le gaz est refroidi ou son refroidissement est augmenté 5 si la température maximale s'approche de la limite admissible). Un tel exemple est illustré à la figure 6. Le processus de la figure 6 se distingue de celui de la figure 5 uniquement en ce que, à l'issue de l'étape 17 de calcul de l'état de charge SOC, le processus comprend une étape 27 durant laquelle le nombre de Richardson Ri du gaz est comparé à une valeur de référence Vr (ici Vr=1). Si le nombre de 10 Richardson R1 excède la valeur de référence (« Y, » étape 21) la vitesse de remplissage peut être maintenue ou augmentée et le processus retourne à l'étape 16 de calcul de la température moyenne du gaz Tgas,average en fonction du temps t et de la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps. Si le nombre de Richardson Ri est inférieur à la 15 valeur de référence (« N », étape 21) la température maximale du réservoir Twall,max(t) est comparée au seuil limite (moins un facteur de sécurité TS) cf. étape 20. Selon que cette température maximale du réservoir Twall,max(t) atteint le seuil limite (moins un facteur de sécurité TS (« Y étape 122), une nouvelle température de gaz d'entrée est calculée (c'est-à-dire que la température du gaz 20 de remplissage est abaissée). Le processus revient à l'étape 11 dans laquelle la température du gaz Tin(t) est fournie au modèle. Si en revanche cette température maximale du réservoir Twall,max(t) reste inférieure au seuil limite (moins un facteur de sécurité TS (« Y étape 121), une nouvelle température de gaz d'entrée est calculée (c'est-à-dire que la 25 température du gaz de remplissage est augmentée). Le processus revient à l'étape 11 dans laquelle la température du gaz entrant Tin(t) est fournie au modèle. Comme illustré schématiquement à la figure 2, lors d'un remplissage, lorsque la température du réservoir (1) calculée atteint un seuil haut (SH) 30 déterminé, le flux et notamment le débit de gaz est diminué et/ou la température du gaz fourni au réservoir (1) est diminuée par échange thermique avec une source de froid. De même, lors d'un soutirage, lorsque la température du réservoir calculée atteint un seuil bas déterminé (SB), le flux de gaz peut être diminué. 3036159 16 La figure 4 illustre plus en détail une prise en compte possible du nombre de Richardson Ri. Le nombre de Richardson Ri pour le gaz dans le réservoir 1 en fonction du temps est calculé. Le procédé comprend une étape 27 de comparaison du 5 nombre de nombre de Richardson Ri calculé avec une valeur de référence Vr déterminée comprise entre 0,05 et 1,5 et de préférence entre 0,05 et 0,15 et notamment égale à 0,1. Lorsque nombre de Nombre de Richardson Ri calculé est inférieur à valeur de référence Vr (étape 127) déterminée, la température du gaz dans le réservoir 1 est considérée homogène c'est-à-dire que la température 10 maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps est considérée égale à la température moyenne de la paroi du réservoir Twall,average(r,t) en contact avec le gaz en fonction du temps (t) : Twall,max(t)= Twall,average(r,t), r étant le rayon à partir de l'axe longitudinal de symétrie du réservoir cylindrique. Au niveau de l'interface en contact avec le gaz (c'est-à-dire 15 r=rayon au niveau de ce qui forme le « liner » du réservoir), r=rayon du liner. Lorsque le nombre de Nombre de Richardson Ri calculé est supérieure à valeur de référence Vr déterminée, la température du gaz dans le réservoir 1 est considérée hétérogène c'est-à-dire que la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps n'est pas égale à la 20 température moyenne de la paroi du réservoir Twall,average(r=rayon du liner,t) en contact avec le gaz en fonction du temps (t), et dans ces conditions, le procédé peut comporter une étape d'augmentation du débit fourni au réservoir 1 pour diminuer la valeur du Nombre de Richardson (Ri) calculé au-dessous de la valeur de référence Vr déterminée (étape 227 figure 4). 25 Un exemple non limitatif de modèle de calcul de la température moyenne du gaz dans le réservoir Tgas,average(t) et de la température extrême Twall,max(t) (maximale ou minimale) de la paroi du réservoir va être décrit à présent. Ledit modèle peut être basé sur : 30 - un bilan de masse et d'énergie du gaz dans le réservoir 1, - une équation d'état du gaz, 3036159 17 des corrélations avec des coefficients sans dimension modélisant les échanges de chaleur entre le gaz et la paroi du réservoir, et entre la paroi du réservoir et l'environnement extérieur une équation de chaleur à une dimension dans la paroi du réservoir, 5 une corrélation entre la température maximale de la paroi du réservoir en fonction du temps et la température moyenne du gaz et de la paroi en contact avec le gaz en fonction du temps. Cette corrélation peut être obtenue par des essais et/ou des simulations. On pourra se référer par exemple au document W02013014346A1 (ou 10 l'article « Evaluating the temperature inside a tank during a filling with highly- pressurized gas », publié en 2014, auteurs: Thomas Bourgeois and àl. Seoul (Korea) : Proceedings of the 20th World Hydrogen Energy Conference, 2014. Des détails de calcul dudit modèle vont être décrits ci-dessous. Par soucis de simplification, le réservoir est considéré rempli par du gaz. Cependant, 15 l'adaptation du modèle au cas d'un soutirage ou-d'une stabilisation (ni remplissage ni soutirage) sera décrite ap.rès. Le modèle combine les bilans de masse et d'énergié du gaz et de l'équation d'état du gaz. Pour un remplissage, la température et la pression du gaz dans le réservoir 20 sont considérées homogènes. Le gaz entrant dans le réservoir 1 possède une enthalpie hl, et le gaz dans le réservoir est considéré échanger de la chaleur avec la paroi via un coefficient d'échange thermique kg. La masse de gaz dans le réservoir est considérée varier directement selon le débit de gaz entrant. Dans ce cas, dm 25 dt film Dans la suite, la variation de masse en fonction du temps sera appelée dm_dt 111 La première équation d'énergie du modèle et le bilan d'enthalpie sont appliqués au système ouvert de l'intérieur du réservoir 1. On néglige l'énergie 30 cinétique et les variations d'énergie gravitationnelle. En connaissant le volume V du réservoir, la surface interne Shit du réservoir et l'enthalpie spécifique du gaz on a l'expression : 3036159 18 dh dP mdt = V(--F kgSintCrwi - di(hin - h) (expression 1) Avec ces hypothèses, les variations de l'enthalpie sont dues à trois facteurs : la compression du gaz, l'enthalpie entrant, et les échanges de chaleur avec la paroi. La seconde équation d'énergie est la définition de la variation 5 d'enthalpie d'un gaz réel : dP dh = cp dT + (1 - 13T) -p (expression 2) En. combinant les expressions 1 et 2 on obtient une équation décrivant l'évolution de la température du gaz dans le réservoir en fonction de l'augmentation de pression, de la température de la paroi, du débit de gaz et de 10 l'enthalpie du gaz entrant. dT dP ni CP dt = VfIT - kg Si (Twi T) rii(he - h) (expression 3) Pour compléter le modèle, l'équation suivante d'état d'un gaz réel peut être utilisée : P V= ri R Z(T, r(expression 4) 15 Pour estimer l'évolution de la température du gaz le système d'équations peut être discrétisé en considérant que certaines dérivées peuvent être calculée comme des variations et que certaines variables à l'instant t sont proches des valeurs à l'instant (t-1) précédent. Des paramètres thermodynamiques Cp, p, h et Z peuvent être estimés pour 20 chaque pression et température en utilisant les tables du standard NIST. A ce stade deux discrétisations peuvent être faites selon le choix des paramètres d'entrée : débit massique Q(t), ou pression du gaz P(t). Le terme « paramètre d'entrée » désigne une variable est connue soit par sa mesure (exemple la pression mesurée dans le conduit de remplissage/soutirage) soit 25 connue car renseignée dans le modèle (par exemple une augmentation de pression de 0,2bar par seconde). Les paragraphes suivants concernent la discrétisation des équations dans le cas de la donnée d'entrée-débit massique ou pression. En combinant et discrétisant les expressions 3 et 4 on obtient : 3036159 19 mt-i 1 - At kgt_iSint(T'it_i - Tt_i) + (mt - mt-i)(htntt-i m Tt = t-i 1 - At-1 mt At-i) avec, mt = mt-i + QtAt avec, . At-i 9t-im Pt-1Tt-1R ht_l ) 5 A partir der ces équations, en connaissant l'état au temps précédent (t-1), en plus du signal de débit Q(t) au temps t, il est possible de déterminer la température du gaz au temps t. A partir de l'équation d'état du gaz, connaissant la pression P au temps (t-1) 10 la température du gaz T au temps t et la masse m(t) au temps t, la pression P(t) au temps t peut être calculée. Avec les notations suivantes : VPt-iTt-i(Pt--Pt-1) mt-1 A - kg t_iStnewit-i- Tt-i) At (htutt-i ht-i) At 15 B = RZ t_i, P At Il est possible de discrétiser les équations précédentes pour obtenir l'expression de la température T à l'instant t. Apt-1T teint-i It + -t-1 A t-i 20 Ici à nouveau, en connaissant le système au temps précédent (t-1) et la valeur de pression au temps t il est possible de calculer la température du gaz au temps t. A partir de l'équation_ d'état, en connaissant Tgas,average(t) et P(t) il est possible de calculer la masse de gaz dans le réservoir au temps t. 25 Dans le cas d'un soutirage, la température du gaz dans le réservoir peut être calculée. La seule différence avec les équations précédentes réside dans le fait que l'enthalpie à l'entrée est à présent l'enthalpie à la sortie et est considérée P tV M Tt 3036159 20 égale à celle du gaz dans le réservoir. Ainsi le terme (hin-h) de l'expression est nul. La donnée d'entrée dans ce cas peut être soit la pression dans le réservoir ou le débit massique sortant. La modélisation des échanges de chaleur est un paramètre important du modèle. Au contraire de modélisations complexes via les équations de Navier Stokes, les êchanges thermiques entre le gaz et la paroi peuvent être modélisés via des corrélations basées sur des nombres sans dimension. Dans le cas d'un remplissage de réservoir à haute pression 200bar, 700bar ou 1000bar, on préconise une corrélation basée sur les nombres de Nusselt, 10 Rayleigh et Prandtl NUDinty Rame and Rech.) par exemple selon la méthode décrite dans le document EP2824378A1. L'expression est par exemple : NuDint = a RaDintb + c Redind Le nombre de Nusselt (NuDint) est basé sur le diamètre interne du réservoir 15 et représente les 'échanges de chaleur convectifs entre le gaz et la paroi. La corrélation est basée sur deux termes. Un premier terme représente la convection naturelle (basée sur le nombre de Rayleigh) tandis que le second terme représente les échanges de convection forcée et dépend du nombre de Reynolds. Les coefficients a, b, c 'et d de corrélation sont supposés constants et 20 dépendent uniquement de la géométrie du réservoir et la nature de l'écoulement de gaz en son sein. Ce modèle peut déterminer cette expression (et donc les coefficients) par des essais. Ils sont donc supposés connus et figés pour différentes conditions de remplissage. 25 Dans le cas d'un soutirage l'expression peut être NUDint = a Rgnint b C'est-à-dire que les échanges de chaleur sont dus uniquement à la convection naturelle basée sur le nombre de Rayleigh. Ce type de corrélation est bien connu dans la littérature. 3036159 21 Lorsque du vent souffle autour du réservoir, les échanges de chaleur externes entre le réservoir et son environnement peuvent être modélisés avec une équation de convection forcée entre l'air et un cylindre selon une formule du type : kaDext 0,5 2 Nua = = ((O. 4 Reije, + O. 06 Re ) airl14) A-air Dext 5 Si le vent est nul le nombre de Reynolds peut être considéré nul. Il est possible de choisir une corrélation de convection libre. La modélisation du gaz, des échanges de chaleur entre le gaz et la paroi ont ainsi- été explicités, le principe du modèle calculant l'évolution de la température de la paroi va à présent,être décrit. 10 Pour résoudre l'équation de chaleur dans la paroi, la paroi va être modélisée selon une dimension. Un bilan de chaleur est réalisé dans un élément de volume élémentaire dV compris entre les portions de rayon r et r+dr (par rapport l'axe longitudinal du réservoir). 15 ' Ce volume élémentaire est un cylindre clos supposé homogène en température (T(r,t)). Ce volume élémentaire a une épaisseur homogène dr et un diamètre interne r. sa longueur interne est Lint+2(r - rint) avec Lint et ont étant respectivement la longueur interne et le rayon d'un cylindre équivalent à zéro dimension (« OD »). 20 Les échanges de chaleur avec le volume élémentaire et à travers un flux de chaleur jr peut être exprimé selon la loi de Fourier j = -À gradT(r, t). Le bilan de chaleur dans le volume élémentaire peut être exprimé comme suit : OT ai- 8T pcpdV -à7 = -i1r Or 1 Sr + 2i-r+dr nr r+ci I Sr+dr r ur 25 A l'interface entre le gaz et la paroi interne (liner) on considère une continuité de la puissance au temps t + hit. Par conséquent : 3036159 22 P = kg(t Lit) Aint (7' gcts(t + - T(runent + lit)) = -2-(runer,t + At) gradT . dA Aint OT = -2.(runer,t + Lit) -dr dA A int Au niveau de l'interface entre le composite et le liner on considère une égalité des flux (indice « w » pour le flux dans un sens, par exemple ouest, et indice « e » pour le sens opposé). t + At ) (-0T) = -2(r e, t + At ) HaT) 5 Or w ar e On considère un unique point de température pour cette interface entre le composite et le liner exprimée par T(rLC, t). Au niveau de l'interface avec l'environnement on a l'expression : t + At ) (-ai) ka(t + At)(T(D't, t + At) - Tamb) ECT (nD ext, t Atr - T amb4) ôr 10 Le terme de- gauché représente le flux dans le réservoir tandis que le terme de droite est le flux en dehors du réservoir, à la surface du réservoir. La discrétisation des équations de paroi à une dimension a été décrite précédemment avec la discrétisation à zéro dimension de la masse de gaz et le bilan d'énergie. De plus la corrélation précédemment décrite des flux de chaleur 15 avec la paroi permet de déterminer (calculer) les valeurs de température moyenne du gaz en fonction du temps Tgas,average(t) et de température moyenne de la paroi du réservoir Twall,average(r,t). A partir ce là, le dispositif doit déterminer la température maximale de la paroi Twall,max(t). 20 Pour cela le procédé utilise une corrélation entre la température moyenne de la paroi du réservoir Twall,average(r,t) et la température maximale de la paroi Twall,max(t). Les hétérogénéités de température dans le réservoir et la paroi sont considérées dépendre essentiellement des débits et vitesses de gaz. Une 25 corrélation possible a la forme suivante : L D ..) Tweil,max(0 = f (Twall,average(r. liner ,t),T gas,average(0, -D , --d ,m 3036159 23 Pour déterminer cette corrélation (fonction f) des expérimentations (remplissages, soutirages) peuvent être réalisées en mesurant les températures de gaz et les températures obtenues au niveau de la paroi. Des simulations à deux ou trois dimensions peuvent permettre de calculer les points chauds/froids 5 pendant le remplissage/soutirage. A partir de mesures expérimentales il a été observé que le nombre de Richardson est très utile pour déterminer les conditions d'homogénéité ou d'hétérogénéité de la température dans le réservoir. Le nombre de Richardson est donné par Gr 10 Ri = Re2 Gr = v2 Avec g Le nombre de Reynolds peut être écrit de la façon suivante : P inv indin Redira = m Avec vin = 4 Ilitdrn Pfiv bottomp int 15 Ou bien ReDint= Fig Avec d inV in vbottom = 6, 2 X , I" int efté-expression est dépendante du temps t et peut être calculée à chaque étape par le modèle. Le nombre de Richardson calculé peut être comparé à une valeur de référence Vr. Ceci permet d'indiquer le niveau d'homogénéité pendant le processus de remplissage/soutirage. Pour des réservoirs de forme générale cylindrique avec un ratio L/D (Longueur L sur diamètre D) inférieur à 4,5 (L/D 1 indique des conditions non homogènes. eigi Tg - Twi alt 20 25 3036159 24 Pendant le transfert de gaz les conditions de transfert peuvent être adaptées pour maintenir les conditions d'homogénéité. Dans ces conditions d'homogénéité, la température maximale atteinte par la paroi du réservoir Twall,max(t) est la température moyenne de la paroi en contact avec le gaz. 5 Aucune corrélation n'est nécessaire dans ce cas entre les température moyenne de paroi en contact avec le gaz et maximale de paroi. Pour des réservoirs dans lesquels le ratio L/D>4,5 les conditions d'hétérogénéité peuvent être considérées toujours présentes. Une corrélation est dans ce cas nécessaire. 10 Ce procédé (calcul) peut être appliqué pendant le transfert de gaz dans une station de remplissage/soutirage. Bien entendu, ces calculs peuvent être réalisés a priori pour chaque type de réservoir pour préétablir les conditions optimales de transfert de gaz. Des simulations peuvent notamment être faites pour déterminer des 15 vitesses de remplissages différentes et les profils de température obtenus. De cette façon, il est possible de déterminer les conditions optimales de remplissage au préalable (vitesse, débit, refroidissement). L'invention s'applique aussi bien au remplissage pour contrôler l'échauffement du réservoir qu'au soutirage pour contrôler le refroidissement du 20 réservoir. 3036159 25 Nomenclature et unités des termes utilisés surface (m2) A a Conductivité thermique (m2.s-1) cp Chaleur massique du gaz à pression constante (J.kg-1.K1) cv Chaleur massique du gaz à volume constant (J.kg-1.K1) Diamètre de l'injecteur de gaz dans le réservoir (m) ci D Diamètre du reservoir (m) Épaisseur (m) g Accélération due à la gravité (m.s-2) h Enthalpie spécifique (J.kg-1) Coefficient d'échange de chaleur (W.m-2.K-1) L Longueur du reservoir (m) m Masse dans le réservoir (kg) Débit massique dans le réservoir (également « Q') (kg.s-1) M Masse molaire du gaz utilisé (kg.mo1-1) P Pression du gaz (Pa) Po Pression atmosphérique (Pa) Flux de chaleur (J.s-1.m-2) r Rayon, coordonnée dans la paroi du réservoir à partir de l'axe longitudinal (m) R Constante des gaz parfaits (J.mo1-1.K-1) S Surface du reservoir (m2) t Temps (s) T Température (K) Tfe Température de l'air au contact de la paroi extérieure du réservoir (K) Tfi Température du gaz en contact avec la paroi interne du réservoir (K) Twe Température de la paroi en contact avec l'air (K) Tw, Température de la paroi en contact avec le gaz (K) / Volume du reservoir (m3) Volume molaire (m3.mo1-1) V Volume spécifique (m3.kg-1) Z Facteur de compressibilité du gaz utilise dans l'équation des gaz réels (sans dimension) Coefficient d'expansion thermique (K1) Y Rapport de chaleurs massiques r/ r Emissivité de la surface externe du réservoir (sans dimension) A Conductivité thermique (W.m-1.K-1) p Viscosité dynamique (Pa.$) Coefficient Joule-Thomson du gaz (K.Pa-1) p 1-4JT Densité (kg.m-3) Vitesse du gaz (m.s -1) 3036159 26 Signification des indices air Air amb Ambiant comp Propriété relative au composite ext Paroi externe f Final g Gaz in Propriété ou nature au niveau de l'entrée du réservoir int Interne liner Propriété de l'enveloppe interne (liner) du reservoir 0 Initial w Propriété de la paroi Paramètre sans dimension Nug Nombre de Nusselt du gaz, Nug = koint Nliair RaDint 5 RaDext RedIn RaDext 10 Prair kaDext Nombre de Nusselt de l'air, Nuair = , air Nombre de Rayleigh du gaz Nombre de Rayleigh de l'ambiant Nombre de Reynolds du flux Nombre de Reynolds de l'air Nombre de Prandtl de l'air, Prair gair cpair REVENDICATIONS1. Procédé de remplissage ou de soutirage d'un réservoir (1) de gaz sous pression, notamment de carburant tel que de l'hydrogène gazeux, le réservoir étant délimité par une paroi de forme générale cylindrique ayant des dimensions et des propriétés thermo-physiques déterminées et connues, le procédé comprenant une régulation du débit du flux de gaz introduit ou, respectivement du gaz soutiré, et/ou de la température dudit gaz introduit, pour éviter que le réservoir n'atteigne un seuil haut déterminé de température ou un seuil bas déterminé de température, le procédé étant caractérisé en ce qu'il comprend une étape d'estimation en temps réel par calcul d'au moins une température du réservoir parmi : la température moyenne de la paroi du réservoir Twall,average(r,t) en fonction du temps (t), la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps, la température minimale atteinte par la paroi du réservoir Twall,min(t) en fonction du temps et en ce que, la régulation du débit du flux de gaz ou de la température dudit gaz est réalisée en fonction de ladite température du réservoir calculée. 2. Procédé selon la revendication 1, caractérisé en ce que, lors d'un remplissage, lorsque la température du réservoir (1) calculée atteint un seuil haut (SH) déterminé, le flux et notamment le débit de gaz est diminué et/ou la température du gaz fourni au réservoir (1) est diminuée par échange thermique avec une source de froid. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que, lors d'un remplissage, lorsque la température du réservoir calculée est inférieure au seuil haut (SH) d'une valeur déterminée, le flux de gaz et notamment son débit est augmenté et/ou la température dudit gaz et/ou du réservoir est augmentée par échange thermique avec une source de chaleur. 4. Procédé selon l'une quelconque des revendications 1 à 3 caractérisé en ce qu'il comprend une étape de calcul du Nombre de Richardson (Ri) pour le gaz dans le réservoir (1) en fonction du temps, une étape de comparaison du nombre de Nombre de Richardson (Ri) calculé avec une valeur de référence (Vr) déterminée comprise entre 0,05 et 1,5 et de préférence entre 3036159 28 0,05 et 0,15 et, lorsque nombre de Nombre de Richardson (Ri) calculé est inférieur à valeur de référence (Vr) déterminée le température du gaz dans le réservoir (1) est considérée homogène c'est-à-dire que la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps 5 est égale à la température moyenne de la paroi du réservoir Twall,average(r,t) en contact avec le gaz en fonction du temps (t) : Twall,max(t)--= Twall,average(r=r_liner,t). 5. Procédé selon la revendication 4 caractérisé en ce que lorsque le nombre de Nombre de Richardson (Ri) calculé est supérieure à valeur de 10 référence (Vr) déterminée, la température du gaz dans le réservoir .(1) est considérée hétérogène c'est-à-dire que la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps n'est pas égale à la température moyenne de la paroi du réservoir en contact avec le gaz en fonction du temps (t) c'est-à-dire au niveau de ce qui forme le « liner » du 15 réservoir, Twall,average(r= rayon du liner du réservoir, t), et dans ces conditions, le procédé comporte une étape d'augmentation du débit de gaz fourni au réservoir (1) pour diminuer la valeur du Nombre de Richardson (Ri) calculé au-dessous de la valeur de référence (Vr) déterminée et rendre ainsi le gaz homogène en température. 20 6. Procédé selon la revendication 4 ou 5, caractérisé en ce que, lors d'un remplissage, lorsque le Nombre de Richardson (Ri) est supérieur à la valeur de référence (Vr), le flux de gaz et notamment son débit est augmenté. 7. Procédé selon l'une quelconque des revendications 4 à 6, caractérisé en ce qu'il comprend, avant remplissage ou, respectivement, 25 avant soutirage, une étape de détermination ou de détection par capteur(s), de la température initiale T(0) du gaz dans le réservoir (1), de la pression initiale P(0) du gaz dans te réservoir (1), de la température moyenne initiale de la paroi du réservoir Tw, average(0) et une étape de détermination de la masse de gaz initiale dans le réservoir m(0) puis, lors d'un remplissage 30 lorsque le nombre de Nombre de Richardson (Ri) calculé est inférieure à valeur de référence (Vr) déterminée le température du gaz dans le réservoir (1) est considérée homogène, c'est-à-dire que la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps est égale 3036159 29 à la température moyenne de la paroi du réservoir Twall,average(r= r=rayon du liner du réservoir,t) en contact avec le gaz en fonction du temps (t), et dans ces conditions, le procédé comporte, au cours du remplissage, une étape de calcul de la température Tgas,average(t) moyenne du gaz dans le 5 réservoir (1) en temps réel en fonction du temps et de la température moyenne de la paroi Twall,average(r,t) en temps réel en fonction du temps (t) à partir d'un bilan de masse et d'enthalpie appliqué au gaz dans le réservoir (1) et à partir également d'un bilan d'énergie dans la paroi du réservoir, de l'équation d'état du gaz (notamment équation des gaz réels), et d'un bilan 10 des échanges thermiques entre le gaz et la paroi, et entre la paroi du réservoir est l'extérieur. 8. Procédé selon l'une quelconque des revendications 4 à 7, caractérisé en ce qu'il comprend, avant remplissage ou, respectivement, avant soutirage, une étape de détermination ou de détection par capteur(s), 15 de la température initiale T(0) du gaz dans le réservoir (1), de la pression initiale P(0) du gaz dans le réservoir (1), de la température moyenne initiale de la paroi du réservoir Tw, average(0) et une étape de détermination de la masse de gaz initiale dans le réservoir m(0) puis, lors d'un remplissage lorsque le nombre de Nombre de Richardson (Ri) calculé est supérieure à 20 valeur de référence (Vr) déterminée la température du gaz dans le réservoir (1) est considérée hétérogène c'est-à-dire que la température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps n'est pas égale à la température moyenne de la paroi du réservoir Twall,average(r)=r_liner,t) en contact avec le gaz en fonction du temps (t), et 25 dans ces conditions, le procédé comporte, au cours du remplissage, une étape de calcul de la température Tgas,average(t) moyenne du gaz dans le réservoir (1) an temps réel en fonction du temps (t) et de la température moyenne de la paroi Twall,average(r,t) en temps réel en fonction du temps à partir d'un bilan de masse et d'enthalpie appliqué au gaz dans le réservoir 30 (1) et d'un bilan d'énergie dans la paroi du réservoir, de l'équation d'état du gaz (équation des gaz parfaits ou réels), et d'un bilan des échanges thermiques entre la paroi du réservoir est l'extérieur et en ce que le procédé comprend une étape de calcul de la température maximale atteinte par la 3036159 30 paroi du réservoir Twall,max(t) en fonction du temps, cette température maximale atteinte par la paroi du réservoir Twall,max(t) en fonction du temps étant obtenue par corrélation à partir de la température Tgas,average(t) moyenne du gaz dans le réservoir (1) calculée en temps réel en fonction du 5 temps et en fonction de la température moyenne de la paroi Twall,average(r,t) en temps réel en fonction du temps (t). 9. Procédé selon la revendication 7 ou 8, caractérisé en ce qu'il comprend une étape de calcul de l'enthalpie hin(t) du gaz entrant ou sortant du réservoir (1) en fonction du temps, une étape de mesure ou de calcul de 10 la masse de gaz m(t) introduite ou soutiré du réservoir en fonction du temps ou, respectivement, une étape de détermination de la pression P(t) dans le réservoir (1) en fonction du temps, le procédé comprenant une étape de détermination de la température moyenne du gaz Tgas,average(t) à l'instant t dans le réservoir (1) en degré K, cette température moyenne Tgas,average(t) 15 étant exprimée comme une fonction du premier degré de la température moyenne du gaz T(t-1) à l'instant précédent (t-1) et d'un coefficient d'échange chaleur convectif entre le gaz et la paroi interne du réservoir (1) à l'instant (t-1) en W.rn2.K-1, dans lequel le coefficient d'échange chaleur kg(t-1) est donné par la relation kg=(Ag/Dint).Nuint dans laquelle Ag est la 20 conductivité thermique du gaz dans le réservoir en W.m-1.1 10.Procédé selon l'une quelconque des revendications 7 à 9, caractérisé en ce qu'il comprend une étape d'estimation de la température maximale atteinte dans l'épaisseur du réservoir Twall,max(t) en fonction du temps, ou respectivement la température minimale atteinte dans l'épaisseur 5 du réservoir Twall,min(t) en fonction du temps, à partir de la température Tgas,average(t) moyenne du gaz dans le réservoir (1) calculée en temps réel en fonction du temps (t), ladite température maximale ou minimale étant obtenue à partir de tableau(x) de correspondance obtenu(s) par des expérimentations, et/ou à partir de calcul et/ou de simulations de sorte qu'à 10 une température Tgas,average(t) moyenne du gaz dans le réservoir (1) calculée en temps réel en fonction du temps (t) correspond, selon les conditions connues de débit de gaz, du diamètre d'injecteur de gaz dans le réservoir et des dimensions et caractéristiques du réservoir (1), une température maximale ou, respectivement minimale, dans l'épaisseur du 15 réservoir en fonction du temps. 11.Procédé selon l'une quelconque des revendications 4 à 9, caractérisé en ce que le Nombre de Richardson (Ri) pour le gaz dans le réservoir (1) en fonction du temps est calculé à partir du nombre de Grashoff (Gr) pour le gaz et du nombre de Reynolds (Re) pour le gaz au temps (t) 20 selon la formule suivante : Ri= Gr/Re2, dans lequel les nombres de Grashof (Gr) et de Reynolds (Re) sont des données connues ou calculées à partir de la valeur de pression du gaz mesurée ou de la masse de gaz dans le réservoir et de la température du gaz Tgas,average(t) moyenne du gaz dans le réservoir (1). 25 12.Procédé selon l'une quelconque des revendications 1 à 10, caractérisé en ce que, lors d'un soutirage, lorsque la température du réservoir calculée atteint un seuil bas déterminé (SB), le flux de gaz est diminué. 13.Procédé selon l'une quelconque des revendications 1 à 12, 30 caractérisé en ce qu'il comprend une étape de mesure ou de calcul en continu de l'un au moins parmi : la pression Pin(t) au temps (t) du gaz introduit dans ou soutiré du réservoir, la température Tin(t) au temps (t) du gaz introduit dans ou soutiré du réservoir, la P(t) pression dans le réservoir 3036159 32 au temps (t), la température ambiante (Tamb(t) au temps (t), la masse de gaz m(t) dans le réservoir (1) au temps (t). 14. Procédé selon l'une quelconque des revendications 1 à 13, caractérisé en ce qu'il est mis en oeuvre par une station de remplissage de 5 réservoirs de gaz hydrogène comprenant au moins une source (10) d'hydrogène à haute pression, au moins une conduite (2) de transfert de gaz reliant sélectivement la source (1) à un réservoir (1), et un organe (4) électronique d'acquisition, de stockage et de traitement de données tel qu'un ordinateur ou microprocesseur, ledit organe (4) pilotant le transfert de gaz 10 entre la source (10) et le réservoir (1), caractérisé en ce que la logique (4) électronique est programmée pour calculer en temps réel au moins une température du réservoir (1) parmi : la température moyenne de la paroi du réservoir Twall,average(r,t) en fonction du temps (t), la température maximale atteinte dans l'épaisseur du réservoir Twall,max(t) en fonction du temps, la 15 température minimale atteinte dans l'épaisseur du réservoir Twall,min(t) en fonction du temps et pour réguler le débit du flux de gaz et/ou de la température dudit gaz en fonction de ladite température du réservoir calculée. 15.Dispositif de remplissage ou de soutirage de gaz d'un réservoir (1) 20 comprenant une conduite (2) de transfert comprenant une vanne (3) et raccordable à un réservoir (1), le dispositif comprenant un organe (4) électronique d'acquisition, de stockage et de traitement de données tel qu'un ordinateur ou microprocesseur, ledit organe (4) pilotant le transfert de gaz entre la source (10) et le réservoir (1), caractérisé en ce que la logique (4) 25 électronique est programmée pour calculer au moins une température du réservoir (1) parmi : la température moyenne de la paroi du réservoir Twall,average(r,t) selon son épaisseur (rayon r) et en fonction du temps (t), la température maximale atteinte dans l'épaisseur du réservoir Twall,max(t) en fonction du temps, la température minimale atteinte dans l'épaisseur du 30 réservoir Twall,min(t) en fonction du temps et pour réguler le débit du flux de gaz en fonction de ladite température du réservoir courante calculée.