FR 3138181 A1 20240126 FR 2207428 A 20220720 Procédé (100) de remplissage d’un réservoir (1) de gaz sous pression, mis en œuvre par un dispositif (10) de remplissage comprenant une source (2) de gaz, une conduite (3) de remplissage reliant la source (2) au réservoir (1), une vanne (4) de régulation de débit et/ou de pression, et un organe (5) de commande électronique configuré pour arrêter le remplissage lorsqu’une température (ou une masse volumique) estimée du gaz présent dans le réservoir (1) atteint une valeur limite de température (ou respectivement une valeur limite de masse volumique), le procédé comprenant, avant le remplissage du réservoir (1), une étape (S1) de détermination de la température ambiante (Tamb) au niveau du dispositif (10) de remplissage, une étape (S2) de détermination de la pression (Pini) du gaz présent dans le réservoir (1) et une étape (S3) préliminaire d’estimation de la température initiale (Tini) du gaz présent dans le réservoir (1), la température initiale (Tini) du gaz présent dans le réservoir étant une valeur estimée en fonction de la température ambiante (Tamb) et en fonction de la pression (P) du gaz présent dans le réservoir avant remplissage, la température initiale (Tini) du gaz présent dans le réservoir étant supérieure ou égale ou inférieure ou égale à la température ambiante (Tamb) Figure de l’abrégé : Fig. 14 Dispositif et procédé de remplissage de réservoir de gaz sous pression L’invention concerne un dispositif et un procédé de remplissage de réservoir de gaz sous pression. L’invention s’applique en particulier au remplissage de réservoirs d’hydrogène. L’invention concerne plus particulièrement un dispositif de remplissage d’un réservoir de gaz sous pression, comprenant une source de gaz, une conduite de remplissage reliant la source au réservoir, une vanne de régulation de débit et/ou de pression dans la conduite de remplissage, un capteur configuré pour mesurer la pression dans le réservoir et/ou la température ambiante au niveau du dispositif de remplissage, et un organe de commande électronique configuré pour: - arrêter le remplissage lorsqu’une température (respectivement une masse volumique) estimée du gaz présent dans le réservoir atteint une valeur limite de température (respectivement une valeur limite de masse volumique), - estimer avant le remplissage du réservoir, une température initiale du gaz présent dans le réservoir. Dans un dispositif tel que décrit ci-dessus, la température du gaz présent dans le réservoir, et en particulier la température avant remplissage, reste difficilement accessible. De façon connue, la température initiale retenue dans les protocoles de remplissage existants est issue de préconisations de la SAE (“Society of Automotive Engineering”) relatives au remplissage de véhicules à pile à combustible. La température initiale est donnée en fonction de la température ambiante selon une courbe prédéfinie pour le cas chaud et une courbe prédéfinie pour le cas froid. Or, la pratique montre que la température initiale du réservoir à remplir n’évolue pas toujours suivant ces courbes préconisées par la SAE. Ainsi, avec une température initiale du réservoir située en dehors du corridor formé par les courbes préconisées par la SAE, le réservoir se trouve exposé notamment à un risque de surchauffe ou de sur-remplissage. Le risque de surchauffe est accru en cas de remplissages successifs. En effet, alors qu’un remplissage entraîne généralement sur le gaz présent dans le réservoir une élévation de température, les protocoles actuels considèrent que pour un remplissage ultérieur la température initiale du réservoir reste dans le corridor défini par la SAE. Ainsi, en raison d’une température initiale du gaz sous-estimée lors du remplissage ultérieur, la limite de température peut être rapidement franchie, en contradiction avec les prédictions des modèles déployés dans les protocoles de remplissage actuels. 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 dispositif 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 ce que la température initiale est une valeur estimée en fonction de la température ambiante et en fonction de la pression du gaz présent dans le réservoir avant remplissage, la température initiale du gaz dans le réservoir étant supérieure ou égale ou inférieure ou égale à la température ambiante. Ainsi, l’invention améliore la précision avec laquelle est estimée la température initiale du gaz présent dans le réservoir, limitant ainsi le risque de surchauffe ou de sur-remplissage sur le réservoir. L’invention concerne également un procédé de remplissage d’un réservoir de gaz sous pression, mis en œuvre par un dispositif de remplissage comprenant une source de gaz, une conduite de remplissage reliant la source au réservoir, une vanne de régulation de débit et/ou de pression dans la conduite de remplissage, et un organe de commande électronique configuré pour arrêter le remplissage lorsqu’une température (ou une masse volumique) estimée du gaz présent dans le réservoir atteint une valeur limite de température (ou respectivement une valeur limite de masse volumique). Le procédé comprend, avant le remplissage du réservoir, une étape de détermination de la température ambiante au niveau du dispositif de remplissage, une étape de détermination de la pression du gaz présent dans le réservoir et un étape préliminaire d’estimation de la température initiale du gaz présent dans le réservoir. La température initiale du gaz dans le réservoir est une valeur estimée en fonction de la température ambiante et en fonction de la pression du gaz présent dans le réservoir avant remplissage, ladite température initiale du gaz dans le réservoir étant supérieure ou égale ou inférieure ou égale à la température ambiante. Par ailleurs, des modes de réalisation du procédé selon l’invention peuvent comporter l'une ou plusieurs des caractéristiques suivantes : - la température initiale du gaz dans le réservoir est l’une au moins parmi: une première température initiale calculée correspondant à un état du réservoir considéré comme récemment rempli à une première masse volumique initiale, une seconde température initiale calculée correspondant à un état du réservoir considéré comme récemment soutiré jusqu’à une seconde masse volumique initiale; - la première température initiale du gaz dans le réservoir est comprise dans une plage de températures hautes ayant une borne inférieure déterminée supérieure ou égale à la température ambiante et une borne supérieure correspondant à un maximum de limite de température déterminé, par exemple égal à 85°C; - la première température initiale est déterminée à partir d’une première courbe prédéterminée prédictive “TPgas_hot” fournie par un modèle physique prédéterminé simulant un remplissage de référence du réservoir sur la plage de températures hautes; - la seconde température initiale du gaz dans le réservoir est comprise dans une plage de températures basses ayant une borne supérieure déterminée inférieure ou égale à la température ambiante et une borne inférieure correspondant à un minimum de limite de température déterminé, par exemple compris entre zéro et -5°C; - la seconde température initiale est déterminée à partir d’une seconde courbe prédéterminée prédictive “TPgas_cold” fournie par un modèle physique prédéterminé simulant une vidange de référence du réservoir sur la plage de températures basses; - le procédé comprend une étape consistant à modéliser, lors du remplissage, une première courbe de variation de la température (ou respectivement une première courbe de variation de la masse volumique) du gaz présent dans le réservoir, ayant comme condition de départ la première température initiale (respectivement la première masse volumique initiale) du gaz présent dans le réservoir, et/ou une deuxième courbe de variation de la température (ou respectivement une deuxième courbe de variation de la masse volumique) du gaz présent dans le réservoir, ayant comme condition de départ la seconde température initiale (respectivement la seconde masse volumique initiale) du gaz présent dans le réservoir; - le procédé comprend une étape consistant à estimer, lors du remplissage, une courbe de variation de la température (ou respectivement une courbe de variation de la masse volumique) du gaz présent dans le réservoir, ladite courbe étant comprise entre la première courbe de variation de la température (respectivement la première courbe de variation de la masse volumique) et la deuxième courbe de variation de la température (respectivement la deuxième courbe de variation de la masse volumique); - la première température initiale du réservoir et/ou la seconde température initiale est recalculée en cours de remplissage en fonction du débit et en fonction de la température du gaz présent dans la conduite de remplissage, ledit débit et ladite température étant déterminés par calcul et/ou par des capteurs au niveau du dispositif de remplissage; - la valeur limite de température (respectivement la valeur limite de masse volumique) est une valeur fixe déterminée, par exemple 85°C (respectivement, par exemple 24,1 kg/m3 ou 40,2 kg/m3) ou une valeur donnée par une courbe “TPgas_max” de température de référence (respectivement par une courbe “RHO_cold” de masse volumique de référence ), ladite courbe étant fournie par un modèle physique prédéterminé qui simule les conditions thermodynamiques du gaz lors d’un remplissage de référence du réservoir; - le modèle physique est basé sur un système d’équations comprenant au moins l’une parmi: une équation de bilan d’énergie interne appliqué au gaz présent dans le réservoir, une équation de bilan de masse appliqué au gaz présent dans le réservoir, une équation de conservation d’énergie dans une paroi du réservoir, une équation de continuité de flux de chaleur entre le gaz présent dans le réservoir et la paroi du réservoir, une équation de continuité de flux de chaleur entre la paroi du réservoir et l’air ambiant, et une équation de débit reliant un débit massique du dispositif de remplissage à une différence de pression entre le dispositif de remplissage et le réservoir; - la première température initiale du réservoir (respectivement la seconde température initiale) du réservoir) et la pression initiale du gaz présent dans le réservoir sont obtenues par résolution dudit système d’équations, par exemple à l’aide d’un algorithme de type matriciel tridiagonal; - la température ambiante du dispositif de remplissage et la pression initiale du gaz présent dans le réservoir sont déterminées par calcul et/ou mesurées par des capteurs au niveau du dispositif de remplissage; - l’organe de commande électronique est configuré pour contrôler la vanne de régulation de débit et/ou de pression pour réaliser une courbe ou une rampe de pression prédéterminée lors du remplissage; - l’organe de commande électronique est configuré pour simuler et estimer la courbe de variation de la masse volumique (et/ou la courbe de variation de la température) du gaz présent dans le réservoir de façon dynamique pendant le remplissage du réservoir et/ou de façon anticipée, c’est-à-dire avant le remplissage. L’invention peut concerner également tout dispositif ou procédé alternatif comprenant toute combinaison des caractéristiques ci-dessus ou ci-dessous dans le cadre des revendications. D’autres particularités et avantages apparaîtront à la lecture de la description ci-après, faite en référence aux figures suivantes dans lesquelles : illustre schématiquement et de façon partielle un dispositif de remplissage selon l’invention; illustre schématiquement des étapes d’un procédé de remplissage au moyen du dispositif de remplissage; illustre une courbe de température de référence donnant, pour une température ambiante donnée, une évolution de la limite de température du gaz en fonction de sa pression; illustre une évolution par rapport à la température ambiante d’une borne inférieure (respectivement d’une borne supérieure) associée à une plage de températures hautes (respectivement une plage de températures basses) dans laquelle intervient un remplissage de référence; illustre des courbes donnant, pour un remplissage de référence simulé à une température ambiante donnée, la variation de la température du réservoir considéré comme récemment rempli (respectivement considéré comme récemment soutiré), la variation de la température, de la masse volumique et de la pression du gaz présent dans ledit réservoir, la variation de la température et de la pression du gaz en sortie du dispositif de remplissage; illustre une courbe surface correspondant à plusieurs courbes de température de référence établies pour différentes températures ambiantes ; illustre la courbe de température de référence de la avec une variable « température » réduite entre 0 (zéro) et 1 (un); illustre des courbes de température de référence et des courbes prédictives de la première température initiale du réservoir ou du gaz présent dans ledit réservoir; illustre des courbes représentant une évolution de la température du réservoir et de la température du gaz présent dans ledit réservoir, ladite évolution étant liée à une diffusion thermique postérieure à un remplissage de référence; illustre les courbes de la après réduction de la variable « température » ; illustre des courbes prédictives de la deuxième température initiale du réservoir ou du gaz présent dans ledit réservoir; illustre les courbes de la après réduction de la variable « température » et réduction de la variable « pression » entre zéro (0) et un (1); illustre un exemple de mise en œuvre d’une courbe de variation de la température du gaz pour piloter le remplissage du réservoir et prévenir sa surchauffe; illustre un exemple de mise en œuvre d’une courbe de variation de la température du gaz en combinaison avec la courbe de température de référence pour piloter le remplissage du réservoir et prévenir sa surchauffe; représente un exemple de discrétisation radiale d’une paroi du réservoir ; représente des détails d’une maille obtenue après la discrétisation de la . Le dispositif 10 de remplissage de réservoirs 1 de gaz sous pression est par exemple une station de remplissage de réservoirs d’hydrogène sous pression. Le dispositif 10 comprend une source 2 de gaz, une conduite 3 de remplissage reliant la source 2 au réservoir 1, une vanne 4 de régulation de débit et/ou de pression dans la conduite 3 de remplissage, un ensemble de capteur(s) 6a, 6b configuré pour mesurer la pression dans le réservoir 1 et/ou la température ambiante Tamb au niveau du dispositif 10 de remplissage. En outre, le dispositif 10 comprend un organe 5 de commande électronique configuré pour contrôler le remplissage et notamment pour arrêter le remplissage lorsqu’une température (respectivement une masse volumique) estimée du gaz présent dans le réservoir 1 atteint une valeur limite de température (respectivement une valeur limite de masse volumique). L’organe 5 de commande est également configuré pour estimer avant le remplissage du réservoir 1, une température initiale du gaz présent dans le réservoir 1. L’organe 5 de commande électronique comprend par exemple un microprocesseur, un ordinateur ou tout contrôleur électronique approprié. Selon l’invention, la température initiale Tini du gaz présent dans le réservoir 1 est une valeur estimée en fonction de la température ambiante Tamb et en fonction de la pression Pini du gaz présent dans le réservoir 1 avant remplissage, ladite température initiale Tini étant supérieure ou égale ou inférieure ou égale à la température ambiante Tamb. Ainsi, pour une température ambiante Tamb donnée, la température initiale Tini du gaz présent dans le réservoir n’est plus une grandeur fixe comme le prévoit l’art antérieur, mais plutôt une grandeur qui dépend de la pression initiale Pini du gaz présent dans le réservoir 1. De façon avantageuse, la température initiale Tini du gaz dans le réservoir 1 est l’une au moins parmi: une première température initiale Tini_1 calculée correspondant à un état du réservoir 1 considéré comme récemment rempli à une première masse volumique ρini_1, et une seconde température initiale Tini_2 calculée correspondant à un état du réservoir 1 considéré comme récemment soutiré jusqu’à une seconde masse volumique initiale ρini_2. La température ambiante Tamb et la pression initiale Pini sont déterminées respectivement au cours d’une étape S1 et une étape S2 d’un procédé de remplissage mettant en œuvre le dispositif de remplissage et illustré de façon schématique à la . Les étapes S1 et S2 peuvent être exécutées simultanément ou non. Quant à la température initiale Tini, qui se décline en une première température initiale Tini_1 et une deuxième température initiale Tini_2, elle est déterminée au cours d’une étape S3 également illustrée à la . Dans ce qui suit, le réservoir 1 considéré comme récemment rempli sera également assimilé à un réservoir récemment exposé en plein soleil ou dans un milieu relativement chaud ou encore à un réservoir en cours d’un remplissage à chaud. Le réservoir considéré comme récemment soutiré sera également assimilé à un réservoir récemment exposé dans un milieu relativement froid ou un réservoir en cours d’un remplissage à froid. Le réservoir considéré comme récemment rempli et le réservoir considéré comme récemment soutiré sont des réservoirs virtuels modélisés qui peuvent être identiques ou non. Dans le second cas, le réservoir considéré comme récemment soutiré peut notamment présenter un rapport surface/volume supérieur à celui du réservoir considéré comme récemment rempli, ce qui le rend plus apte à évacuer la chaleur et à se refroidir. De façon avantageuse, la première température initiale Tini_1 du réservoir est comprise dans une plage de températures hautes ayant une borne inférieure déterminée TLgas_min_1 supérieure ou égale à la température ambiante et une borne supérieure TLgas_max_1 correspondant à un maximum de limite de température déterminé. L’évolution de la borne inférieure TLgas_min_1 en fonction de la température ambiante du dispositif de remplissage est donnée par une courbe “Thotsoak” (voir ) issue par exemple de préconisations de la SAE relatives au remplissage d’un réservoir chaud. La borne supérieure TLgasmax_1 est fixée par exemple à 85°C. Cette valeur tient compte des propriétés thermophysiques du réservoir. De façon avantageuse, la seconde température initiale Tini_2 du réservoir est comprise dans une plage de températures basses ayant une borne supérieure déterminée TLgas_max_2 inférieure ou égale à la température ambiante Tamb et une borne inférieure TLgas_min_2 correspondant à un minimum de limite de température déterminé, par exemple comprise entre zéro et -5°C. En particulier, la borne supérieure TLgas_max_2 est fonction de la température ambiante et suit une courbe “Tcoldsoak” (voir ) issue par exemple de préconisations de la SAE relatives à un remplissage d’un réservoir froid. De façon avantageuse, la valeur limite de température (respectivement la valeur limite de masse volumique) est une valeur fixe déterminée ou une valeur donnée par une courbe “TPgasmax” de température de référence (respectivement une courbe de masse volumique de référence). En particulier, la courbe “TPgasmax” de température de référence (respectivement la courbe de masse volumique de référence) permet, pour une pression donnée du gaz dans le réservoir, d’estimer la valeur limite de température (respectivement la valeur limite de masse volumique). Un exemple de courbe “TPgas_max” de température de référence est illustré à la . Pour une pression maximale du gaz présent dans le réservoir 1, la valeur limite de température est de 85°C et coïncide avec la borne supérieure TLgasmax_1 de la plage de températures hautes. Pour une pression minimale du gaz dans le réservoir, la valeur limite de température est d’environ 20°C et coïncide avec la borne inférieure TLgasmin_1 de la plage de températures hautes. La valeur limite de masse volumique dépend de la pression de service du réservoir. Par exemple, pour une pression de service de 350 barg, la valeur limite de masse volumique est fixée par exemple à 24.1 kg/m3. Pour une pression de service de 700 barg, la valeur limite de masse volumique est fixée par exemple à 40.2 kg/m3. De manière générale, la valeur limite de masse volumique peut être définie par un pourcentage d’état de remplissage “SOC” (pour “State of Charge”, en anglais) qui définit une quantité maximale d’hydrogène préconisée dans un réservoir. Détermination de la courbe “TPgasmax” de température de référence et de la courbe de masse volumique de référence De façon avantageuse, la courbe “TPgasmax” de température de référence (respectivement la courbe de masse volumique de référence) est fournie par un premier modèle physique qui simule un remplissage de référence du réservoir sur la plage de températures hautes (respectivement sur la plage de températures basses) avec un pré-refroidissement au niveau de la source. Le remplissage de référence est simulé depuis une pression minimale qui est ici fixée à 5 barg jusqu’à une pression maximale correspondant à la pression de service du réservoir. En particulier, la courbe “TPgasmax” de température de référence (respectivement la courbe de masse volumique de référence) est obtenue à partir d’une courbe “Pgashot” donnant l’évolution de la pression du gaz présent dans le réservoir en fonction du temps et à partir d’une courbe “Tgasmax” (respectivement une courbe “RHO_cold”) donnant l’évolution en fonction du temps de la température (respectivement l’évolution de la masse volumique) du gaz contenu dans le réservoir. Les courbes “Pgashot”, “Tgasmax” (respectivement “RHO_cold”) sont également fournies par le premier modèle physique prédéterminé. Elles sont illustrées à la . Pour différentes températures ambiantes, le premier modèle physique prédéterminé prévoit autant de courbes “TPgasmax” de température de référence (respectivement autant de courbes de masse volumique de référence) qui forment une courbe surface de température de référence (respectivement une courbe surface de masse volumique de référence). La illustre un exemple d’une courbe surface de température de référence. Afin de faciliter l’exploitation de la courbe “TPgasmax” de température de référence (respectivement la courbe de masse volumique de référence), la variable “température” est réduite entre la valeur zéro (0) correspondant à la borne inférieure TLgasmin_1 et la valeur un (1) correspondant à la borne supérieure TLgasmax_1. La illustre une courbe ”Reduc_TPgasmax” obtenue après réduction de la variable “température” sur la courbe “TPgasmax” de température de référence. Détermination de la première température initiale Tini_1 du gaz La première température initiale Tini_1 du gaz présent dans le réservoir est donnée par une deuxième courbe prédictive prédéterminée “TPgashot” fonction de la pression initiale Pini du gaz (voir ). En d’autres termes, pour une pression initiale Pini donnée du gaz présent dans le réservoir, la courbe prédictive “TPgashot” permet de déterminer la première température initiale Tini_1. La courbe prédictive “TPgashot” est tirée d’un deuxième modèle physique qui simule un remplissage de référence du réservoir sur la plage de températures hautes avec éventuellement un pré-refroidissement du gaz fourni par la source. De préférence, le deuxième modèle physique simule en plus une diffusion thermique au niveau du réservoir et/ou au niveau du gaz présent dans le réservoir. Ainsi, le deuxième modèle physique peut différer du premier modèle physique par la prise en compte de la diffusion thermique. La diffusion thermique est postérieure au remplissage de référence du réservoir et au pré-refroidissement de la source. De façon avantageuse, la diffusion thermique est simulée sur une période de 2 minutes correspondant au temps minimum requis entre deux remplissages successifs du réservoir sur la station de remplissage. Encore de façon avantageuse, au cours de la période de diffusion thermique, la pression du gaz dans le réservoir est supposée constante. Pour une température ambiante donnée, la diffusion thermique induit une chute de température au niveau du gaz présent dans le réservoir. Une telle chute est représentée par une courbe “STgasmax” illustrée à la . Il est à noter que la courbe “TPgashot” est obtenue à partir de la courbe “Pgashot” et de la courbe “Tgashot” établies toutes les deux pour une même température ambiante. En particulier, la courbe “Tgashot” résulte d’une soustraction opérée entre la courbe “Tgasmax” et la courbe “STgasmax”. De façon avantageuse, afin de faciliter l’exploitation de la courbe “TPgashot”, la variable “température” peut être réduite entre la valeur zéro (0) correspondant à la borne inférieure TLgasmin_1 et la valeur un (1) correspondant à la borne supérieure TLgasmax_1. La courbe “Reduc_TPgashot” ainsi obtenue est illustrée à la . Détermination de la deuxième température initiale Tini_2 du gaz La deuxième température initiale Tini_2 du gaz est obtenue à partir d’une troisième courbe “TPgascold” fonction de la pression initiale Pini du gaz (voir ). En d’autres termes, sur la base de la pression initiale Pini du gaz présent dans le réservoir, la courbe “TPgascold” permet de déterminer la deuxième température initiale Tini_2 du gaz présent dans le réservoir. La courbe “TPgascold” peut être fournie par un troisième modèle physique prédéterminé (décrit plus loin) qui simule une vidange de référence du réservoir sur la plage de températures basses, c’est-à-dire entre la borne supérieure TLgasmax_2 et la borne inférieure TLgasmin_2. La borne supérieure TLgasmax_2 correspond au début de la vidange de référence où le réservoir est dans un état rempli au moins en partie, et de préférence à 100%. A cette borne supérieure TLgasmax_2 correspond une pression maximale Pgasmax_2, dite pression de service (“Nominal Working Pressure” en anglais), et une masse volumique maximale ρgasmax_2. Il est à noter que la pression maximale Pgasmax_2 peut s’exprimer en fonction de la borne supérieure TLgasmax_2 à travers l’expression mathématique suivante tirée de la littérature: Pgasmax_2 [barg] = A * TLgasmax_2 [°C] + B avec A compris entre 1 et 2, et B compris entre 300 et 400. La borne inférieure TLgasmin_2 correspond à la fin de vidange de référence où le réservoir présente une pression minimale Pgasmin_2, ici fixée à 5 barg, et une masse volumique minimale ρgasmin_2, dite seconde masse volumique initiale ρini_2. En outre, la borne inférieure TLgasmin_2 peut être obtenue à partir de l’expression mathématique ci-dessous tirée du troisième modèle physique: Tgasmin_2 = A * Tamb + B où Tamb est la température ambiante au niveau du dispositif de remplissage; A est compris entre 0,5 et 1,5; B est compris entre -30 et -15. Pour différentes températures ambiantes, le troisième modèle physique fournit autant de courbes “TPgascold” qui forment une courbe surface (non illustrée). Afin de faciliter l’exploitation de ces courbes “TPgascold”, la variable “température” sur chacune d’elles est réduite entre la valeur zéro (0) correspondant à la borne inférieure TLgasmin_2 et la valeur un (1) correspondant à la borne supérieure TLgasmax_2. De même, la variable “pression” peut être réduite entre la valeur zéro (0) correspondant à la pression minimale Pgasmin_2 et la valeur un (1) correspondant à la pression maximale Pgasmax_2 avant vidange du réservoir. La courbe “Reduc_TPgascold” illustrée à la est le résultat d’une telle réduction opérée sur la courbe “TPgascold”. En définitive, pour une pression initiale Pini mesurée dans le réservoir et une température ambiante Tamb relevée au niveau du dispositif de remplissage, au moins un modèle physique prédéterminé considère que cette pression initiale Pini est la pression finale d’un remplissage précédent (“cas chaud”) ou d’un soutirage précédent (“cas froid”). Ces deux hypothèses définissent respectivement une valeur pour la première température initiale Tini_1 (correspondant au cas chaud) et une valeur pour la deuxième température initiale Tini_2 (correspondant au cas froid). La valeur de la première température initiale Tini_1 et la valeur de la deuxième température initiale Tini_2 ainsi définies peuvent être différentes de la valeur de température initiale déterminée selon l’art antérieur, par exemple s’écartant de 10 à 20°C degrés de la température ambiante. Estimation de la courbe de variation de la température du gaz entre une première courbe de variation de la température du gaz (cas chaud) et une seconde courbe de variation de la température du gaz (cas froid). De façon avantageuse, l’organe 5 de commande peut être configuré pour estimer, lors du remplissage, une courbe de variation de la température du gaz présent dans le réservoir en fonction de la température initiale du gaz présent dans le réservoir, c’est-à-dire avec comme point de départ (ou condition de départ), la température initiale du gaz présent dans le réservoir. Plus spécifiquement, l’organe 5 de commande peut être configuré pour estimer une telle courbe entre une première courbe de variation de la température du gaz (cas chaud) et une deuxième courbe de variation de la température du gaz (cas froid). Ainsi, pour une pression donnée du gaz présent dans le réservoir, la température correspondante est estimée entre une première température située sur la première courbe de variation de la température du gaz (cas chaud) et une seconde température située sur la deuxième courbe de variation de la température du gaz (cas froid). La première courbe (respectivement la deuxième courbe) de variation de la température du gaz a pour point de départ (ou condition de départ) la première température initiale Tini_1 (respectivement la deuxième température initiale Tini_2) du gaz. En outre, la première courbe (respectivement la deuxième courbe) de variation de la température du gaz est fournie par un modèle physique prédéterminé qui simule un remplissage de référence du réservoir depuis la première température initiale Tini_1 (respectivement depuis la deuxième température initiale Tini_2) du gaz, ledit remplissage étant accompagné d’un pré-refroidissement au niveau de la source et suivi d’une diffusion thermique. Il est à noter que le modèle physique à l’origine de la courbe “TPgashot” prédictive de la première température initiale Tini_1 du gaz est de préférence le même que celui à l’origine de la première courbe de variation de la température du gaz. Les deux courbes se distinguent par leurs points de départ qui sont respectivement la borne inférieure TLgasmin_1 (pour la courbe “TPgashot” prédictive de la première température initiale Tini_1) et la première température initiale Tini_1 (pour la première courbe de variation de la température du gaz). Estimation de la masse volumique du gaz présent dans le réservoir entre une première courbe de variation de la masse volumique (cas chaud) et une seconde courbe de variation de la masse volumique (cas froid) De façon avantageuse, l’organe 5 de commande peut être configuré pour estimer (par exemple calculer), lors du remplissage, une courbe de variation de la masse volumique du gaz présent dans le réservoir en fonction de la masse volumique initiale du gaz présent dans le réservoir. Cela signifie que la courbe de variation de la masse volumique a pour point de départ la masse volumique initiale du gaz présent dans le réservoir. Il est à noter que la masse volumique initiale peut être calculée par l’équation des gaz parfaits ou réels en connaissant le volume du réservoir (connu ou estimé, par exemple via la technique connue basée sur une impulsion “pulse” de pression avant remplissage) , la pression et la température du gaz. Plus spécifiquement, l’organe 5 de commande peut être configuré pour estimer la courbe de variation de la masse volumique du gaz entre une première courbe de variation de la masse volumique du gaz (cas chaud) et une deuxième courbe de variation de la masse volumique du gaz (cas froid). Ainsi, pour une pression donnée du gaz présent dans le réservoir, la masse volumique correspondante est estimée entre une première masse volumique située sur la première courbe de variation de la masse volumique (cas chaud) et une seconde masse volumique située sur la seconde courbe de variation de la masse volumique (cas froid). La première courbe (respectivement la deuxième courbe) de variation de la masse volumique a pour point de départ la première masse volumique initiale (respectivement la seconde masse volumique initiale) du gaz présent dans le réservoir. En outre, la première courbe (respectivement la deuxième courbe) de variation de la masse volumique du gaz est fournie par un modèle physique qui simule un remplissage de référence du réservoir depuis la première masse volumique initiale (respectivement depuis la deuxième masse volumique initiale) du gaz avec un pré-refroidissement au niveau de la source. Détermination de la courbe de température de référence du réservoir et des courbes de variation de la température du réservoir De façon avantageuse, la condition d’arrêt de remplissage du réservoir portant sur une valeur limite de la température du gaz présent dans le réservoir peut être remplacée par une condition d’arrêt portant sur une valeur limite de la température du réservoir. La valeur limite de la température du réservoir peut être donnée par une courbe de température de référence du réservoir, ladite courbe étant fournie par le premier modèle physique. En particulier, pour un réservoir formé de plusieurs couches, par exemple une couche interne dite “liner” en contact avec le gaz et une couche externe dite “composite” en contact avec l’air ambiant, le premier modèle physique fournit trois courbes de température de référence illustrées à la , à savoir une première courbe “TPgasliner_max” pour la paroi interne, une deuxième courbe “TPlinercompo_max” pour l’interface couche interne/couche externe, et une troisième courbe “TPexternal_max” pour la paroi externe. Les courbes “TPgasliner_max”, “TPexternal_max”, “TPlinercompo_max” peuvent être obtenues à partir de la courbe “Pgashot” en combinaison respectivement avec des courbes “Tgasliner_hot”, “Texternal_hot”, “Tlinercompo_hot” donnant l’évolution de la température en fonction du temps pour la paroi interne, l’interface couche interne/couche externe, et la paroi externe. Les courbes “Tgasliner_hot”, “Texternal_hot” sont illustrées à la . En outre, pour être facilement exploitables, les courbes “TPgasliner_max”, “TPexternal_max”, “TPlinercompo_max” peuvent être transformées respectivement en une courbe “Reduc_TPgasliner_max”, une courbe “Reduc_TPexternal_max”, et une courbe “Reduc_TPlinercompo_max” définies chacune entre une valeur zéro (1) correspondant à la borne inférieure TLgasmin_1 de température et une valeur un (1) correspondant à la borne supérieure TLgasmax_1 de température. Estimation de la température initiale du réservoir entre une première température initiale (cas chaud) et une deuxième température initiale du réservoir (cas froid) L’organe 5 de commande peut être configuré pour estimer la température initiale du réservoir entre une première température initiale du réservoir correspondant à un état du réservoir considéré comme récemment rempli à la première masse volumique initiale (cas chaud), et une deuxième température initiale du réservoir correspondant à un état du réservoir considéré comme récemment soutiré jusqu’à une seconde masse volumique initiale (cas froid). La première température initiale Tini_1 du réservoir peut être donnée par une courbe prédictive “TPgasliner_hot”, “TPlinercompo_hot”, “TPexternal_hot” (voir ) fonction de la pression initiale Pini du gaz. En d’autres termes, pour une pression initiale Pini donnée du gaz présent dans le réservoir, la courbe prédictive “TPgasliner_hot”, “TPlinercompo_hot”, “TPexternal_hot” permet de déterminer la première température initiale Tini_1 du réservoir. Les courbes “TPgasliner_hot”, “TPexternal_hot”, “TPlinercompo_hot” peuvent être fournies par le deuxième modèle physique. Plus spécifiquement, ces courbes sont obtenues à partir de la courbe “Pgashot” en combinaison respectivement avec des courbes “Tgasliner_hot_bis”, “Texternal_hot_bis” “Tlinercompo_hot_bis ” (non illustrées) donnant l’évolution de la température en fonction du temps pour respectivement la paroi interne, la paroi externe et l’interface couche interne/couche externe du réservoir. La courbe “Tgasliner_hot_bis” est obtenue à partir de la courbe “Tgasliner_hot” modulée par une courbe “STgasliner”. La courbe “Texternal_hot_bis” est obtenue à partir de la courbe “Texternal_hot” modulée par une courbe “STexternal”. La courbe “Tlinercompo_hot_bis” est obtenue à partir de la courbe “Tlinercompo_hot” modulé par une courbe “STlinercompo”. Les courbes “STgasliner”, “STexternal”, “STlinercompo” sont illustrées à la . Elles représentent une diffusion thermique postérieure au remplissage de référence et qui intervient au niveau de la paroi interne, la paroi externe et l’interface couche interne/couche externe du réservoir. Dans l’exemple illustré à la , la diffusion thermique correspond à une chute de température au cours du temps pour la paroi interne du réservoir. En revanche, pour la paroi externe et l’interface couche interne/couche externe du réservoir, la diffusion thermique correspond à une augmentation de température. Il est à noter que les courbes “TPgasliner_hot”, “TPexternal_hot”, “TPlinercompo_hot” peuvent être réduites respectivement en une courbe “Reduc_TPgasliner_hot”, une courbe “Reduc_TPexternal_hot”, et une courbe “Reduc_TPlinercompo_hot”, définies chacune entre une valeur zéro (1) correspondant à la borne inférieure TLgas_min_1 de température et une valeur un (1) correspondant à la borne supérieure TLgas_max_1 de température (voir ). La deuxième température initiale Tini_2 du réservoir est donnée par une courbe prédictive “TPgasliner_cold”, “TPexternal_cold”, “TPlinercompo_cold” fonction de la pression initiale Pini du gaz (voir ). En d’autres termes, pour une pression initiale Pini donnée du gaz, la courbe prédictive “TPgasliner_cold”, “TPexternal_cold”, “TPlinercompo_cold” permet de déterminer la deuxième température initiale Tini_2 du réservoir. Il est à noter que la courbe prédictive “TPgasliner_cold”, “TPexternal_cold”, “TPlinercompo_cold” peut être fournie par le troisième modèle physique qui simule une vidange de référence du réservoir. A la fin de cette vidange de référence, la température du réservoir atteint une valeur minimale qui varie à travers l’épaisseur du réservoir. Pour un réservoir formé de deux couches, la température minimale Tmin au niveau de la paroi interne, la paroi externe ou l’interface couche interne/couche externe peut être donnée par exemple par une loi linéaire de type Tmin = A*Tamb +B, où Tamb est la température ambiante; A est compris entre 0,5 et 1,5; B est compris entre -20 et -10. Les courbes prédictives “TPgasliner_cold”, “TPexternal_cold”, “TPlinercompo_cold” peuvent être obtenues à partir à de la courbe “Pgascold” (donnant l’évolution de la pression du gaz en fonction du temps) en combinaison respectivement avec les courbes “Tgasliner_cold”, “Texternal_cold”, “Tlinercompo_cold” donnant l’évolution de la température du réservoir en fonction du temps (cas froid). Les courbes “Tgasliner_cold”, ”Texternal_cold” et “Pgas_cold” sont illustrées à la . Estimation de la courbe de variation de la température du réservoir entre une première courbe de variation de la température du réservoir (cas chaud) et une seconde courbe de variation de la température du réservoir (cas froid). L’organe 5 de commande peut être configuré pour estimer, lors du remplissage, une courbe de variation de la température du réservoir (voir étape S5 du procédé illustré à la ). Cette estimation est faite entre une première courbe de variation de la température du réservoir (cas chaud) et une deuxième courbe de variation de la température du réservoir (cas froid). Ladite première courbe et la dite deuxième courbe sont déterminées lors d’une étape S4 du procédé illustré à la . La première courbe de variation a pour point de départ la première température initiale Tini_1 du réservoir. La deuxième courbe de variation a pour point de départ la deuxième température initiale Tini_2 du réservoir. La première courbe (respectivement la deuxième courbe) de variation de la température du réservoir peut être obtenue à partir d’un modèle qui simule un remplissage de référence du réservoir depuis la première température initiale Tini_1 (respectivement depuis la deuxième température initiale Tini_2) du réservoir, avec un pré-refroidissement du gaz au niveau de la source, et une diffusion thermique postérieure au remplissage de référence. Effet de la diffusion thermique sur les courbes prédictives de la première température initiale Tini_1 du réservoir ou du gaz présent dans le réservoir Une analyse de la donnant les courbes de température de référence du réservoir - ou du gaz présent dans le réservoir - (c’est-à-dire les courbes “TPgasliner_max”, “TPlinercompo_max”, “TPexternal_max”, “TPgas_max”), et les courbes prédictives de la première température initiale Tini_1 du réservoir - ou du gaz dans le réservoir - (c’est-à-dire les courbes “TPgasliner_hot”, “TPlinercompo_hot”, “TPexternal_hot”, “TPgas_hot”) montre que le gaz et la paroi interne du réservoir voient leurs températures diminuer à la suite de la diffusion thermique, alors que la paroi externe et l’interface couche interne/couche externe du réservoir voient leurs températures augmenter à la suite de la diffusion thermique. La prise en compte de cette évolution singulière de la température du réservoir à la suite de la diffusion thermique permet d’améliorer la précision de la courbe de température de référence du réservoir et/ou la courbe de température de référence du gaz présent dans le réservoir. Ainsi, l’invention permet d’améliorer la sécurité lors du remplissage en limitant le risque de franchissement de la valeur limite de température. Exemples de mise en œuvre de la courbe “Tgas_max” de température de référence et de la première courbe “Tgas_hot_bis” de variation de la température du gaz Lors du remplissage du réservoir, afin de prévenir toute surchauffe, la première courbe “Tgas_hot_bis” de variation de la température du gaz présent dans le réservoir peut être exploitée seule (voir exemple de la ) ou en combinaison avec la courbe “Tgasmax” de température de référence du gaz présent dans le réservoir (voir exemple de la ). Comme le montre chacun de ces exemples, après raccordement du véhicule à l’extrémité de la conduite de remplissage (étape 1) et initialisation du remplissage par le dispositif de remplissage (étape 2), une première surchauffe est détectée au niveau du gaz présent dans le réservoir (étape 3). Dans l’exemple de la , la surchauffe est détectée au moment où la courbe “Tgas_hot_bis” atteint la valeur limite de température représentée ici par une ligne horizontale à 85°C. Dans l’exemple de la , la première surchauffe est détectée en dessous de la valeur limite de température (85°C). Cette surchauffe correspond au moment où la courbe “Tgas_hot_bis” « rattrape » la courbe “Tgas_max”. A la suite de cette première surchauffe, l’organe de commande arrête le remplissage alors que le réservoir n’est pas rempli à son maximum. L’arrêt se prolonge (étape 4) pendant une certaine durée requise avant un nouveau cycle de remplissage. Les étapes 5, 6, 7 correspondent à ce nouveau cycle de remplissage pour lequel la température initiale du gaz présent dans le réservoir est déterminée en fonction de la température atteinte par le même gaz à la fin du cycle précédent, mais aussi en fonction de la pression de ce même gaz au début du nouveau cycle. Grâce à la prise en compte de la température atteinte lors d’un cycle de remplissage antérieur, l’invention limite lors d’un cycle de remplissage ultérieur le risque de surchauffe, et ouvre la possibilité d’atteindre un maximum de charge dans le réservoir. Bien entendu, les courbes “Tgas_max” et “Tgas_hot_bis” des exemples ci-dessus peuvent être remplacées par les courbes “TPgasmax” et “TPgas_hot” pour piloter le remplissage du réservoir. Lorsque la condition d’arrêt porte sur une valeur limite de la température du réservoir, le remplissage peut être piloté à l’aide de l’une des courbes (cas chaud) donnant la variation de la température du réservoir (par exemple, la courbe “TPgasliner_hot”). Cette courbe de variation peut être exploitée seule ou en combinaison avec la courbe de température de référence associée. (par exemple, la courbe “TPgasliner_max”). Présentation des modèles physiques Les modèles physiques évoqués ci-dessus à l’origine des différentes courbes prédictives présentées dans cette description s’appuient sur un système d’équations de bilan d'énergie appliqué au réservoir et au gaz présent dans le réservoir. Ces équations sont présentées ci-dessous. L’équation de bilan d’énergie interne du gaz dans les réservoirs s’écrit de la façon suivante : (1) Ce bilan fait intervenir la masse du gaz dans les réservoirs , l’énergie interne massique du gaz dans ce réservoir , l’échange de chaleur entre le gaz dans les réservoirs et la paroi interne de celui-ci avec le coefficient d’échange et la surface interne de ces réservoirs , le débit massique total d’hydrogène ainsi que l’enthalpie massique entrante du gaz dans les réservoirs. et sont respectivement les températures moyennes spatiales du gaz dans le volume de ces réservoirs et de la paroi interne de celui-ci. Un des deux critères de sécurité à ne pas dépasser fait intervenir la densité dans le modèle. Comme le volume total des réservoirs du véhicule (Compressed Hydrogen Storage System) peut être connu par l’automate, l’équation (1) peut s’écrire, en utilisant : (2) Pour résoudre cette équation (2), il est nécessaire de la discrétiser en remplaçant les dérivées temporelles par des différences finies. On obtient ainsi l’équation, en divisant par : (3) où est la valeur d’un pas de temps (environ une seconde). L’équation de bilan de matière dans les réservoirs permet d’écrire : (4) soit, en divisant par : (5) L’équation (5) permet d’obtenir la valeur de . L’équation (3) permet donc d’obtenir la valeur de l’énergie interne : (6) L’équation de conservation d’énergie dans la paroi s’écrit : (7) Ce bilan fait intervenir la masse volumique de la paroi ρ, sa capacité calorifique c p et sa conductivité thermique λ, ainsi que le rayon r à partir du centre du réservoir. À l’interface entre le gaz et le liner, l’équation de continuité du flux s’écrit : (8) À l’interface entre le composite et le milieu ambiant, l’équation de continuité du flux s’écrit : (9) Où est le coefficient d’échange de chaleur entre le milieu ambiant et la paroi externe du réservoir, est la surface externe de ces réservoirs et est la température ambiante. Enfin, à l’interface entre les deux matériaux liner et composite à l’intérieur de la paroi, l’équation de continuité du flux s’écrit : (10) Discrétisation des équations dans la paroi Pour résoudre les équations (7), (8), (9) et (10), il est nécessaire de discrétiser radialement la paroi en différentes mailles. Le nombre de mailles dans le liner n_liner et dans le composite n_comp peut varier. Les températures de paroi sont calculées à chaque nœud, et il y a n_mesh = n_liner + n_comp + 1 nœuds. La taille des mailles est donc dans le liner et dans le composite, comme le montre la Figure 13. Les mailles aux interfaces ont une épaisseur de . Dans l’exemple illustré, 3 mailles sont prévues dans le liner et 3 autres mailles sont prévues dans le composite. L’équation (7), pour être résolue, est discrétisée et appliquée pour chacune des différentes mailles constituant la paroi. La Figure 14 montre la discrétisation radiale utilisée pour résoudre l’équation. Sur cette figure, une maille, centrée sur le point P, est entourée d’un point W à l’ouest et d’un point E à l’est, l’intérieur du réservoir se trouvant à l’ouest. Entre les points W et P est définie une face fictive w, et entre les points P et E est définie une deuxième face fictive e. Les distances entre les points W et P d’une part, P et E de l’autre sont notées respectivement et . Les faces w et e sont distantes de . Les points W, P et E sont situés respectivement à , et du centre du réservoir. L’équation (7) est intégrée sur un pas de temps et spatialement d’ouest en est, pour obtenir : (11) Le terme de gauche est appelé A, tandis que le terme de droite est appelé B. En faisant l’hypothèse que la densité ρ et la capacité thermique changent peu entre les instants t et t+Δt, le terme de gauche s’écrit : (12) L’hypothèse est faite que la différence de température reste constante entre les points w et e et est égale à . Le terme A s’écrit : (13) En écrivant : , le terme A s’écrit finalement : (14) Le terme B, terme de droite de l’équation (11), peut s’écrire : (15) En faisant l’hypothèse que et , et en supposant que la conductivité thermique λ change peu entre les instants t et t+Δt, le terme B s’écrit : (16) On note avec un exposant 0 les valeurs prises à l’instant t, et avec un exposant 1 les valeurs prises à l’instant t+Δt. On note avec un indice W , P et E les valeurs prises aux points W, P et E respectivement. Par exemple, la valeur est notée . Un schéma implicite est utilisé, donc l’équation (11) devient : (17) À l’interface entre le gaz et le liner, l’équation de continuité du flux (8) est discrétisée de la manière suivante : (18) À l’interface entre le composite et le milieu ambiant, l’équation de continuité du flux (9) est discrétisée de la manière suivante : (19) Enfin, à l’interface entre les deux matériaux liner et composite à l’intérieur de la paroi, l’équation de continuité du flux (10) est discrétisée de la manière suivante : (20) Le système d’équations (17) pour chaque point P qui ne se trouve pas à une interface et (18), (19) et (20) peuvent être résolues par l’algorithme de Thomas, aussi appelé TDMA (pour Tri-Diagonal Matrix Algorithm ). Propriétés physiques et coefficients Tg et Pg à partir de ρ et u À chaque pas de temps, le modèle calcule la température et la pression du gaz à partir de la densité ρ et l’énergie interne u du gaz : (21) (22) Ces fonctions f 1 et f 2 sont ajustées pour correspondre au mieux aux propriétés thermophysiques du gaz. Coefficients de transferts thermiques ki et ke Le modèle a également besoin de valeurs pour les coefficients de transferts thermiques k i et k e . Pour le coefficient de transfert interne k i , une fonction est utilisée, calibrée à partir de la corrélation présentée par Bourgeois, T. et al. (voir “Bourgeois T, Ammouri F, Weber M, Knapik C. Evaluating the temperature inside a tank during a filling with highly-pressurized gas. International Journal of Hydrogen Energy 2015;40:11748-55”) (23) Pour le coefficient de transfert externe k e , une fonction est utilisée, calibrée à partir de la corrélation présentée par Massard, F. (voir Aide-mémoire du thermicien. ELSEVIER; 1997) : (24.1) (24.2) (24.3) (24.4) (24.5) De façon avantageuse, les modèles physiques décrits ci-dessus peuvent être associés à un modèle de transfert thermique entre le dispositif de remplissage et le ou les réservoirs. Ce modèle de transfert est décrit ci-dessous. Le terme d’enthalpie massique entrante du gaz de l’équation (1) est obtenu à partir de l’enthalpie massique du gaz au niveau du dispenser grâce à un modèle de transfert thermique. Celui-ci a pour but d’estimer les échanges thermiques du gaz entre le dispenser du dispositif de remplissage et l’entrée dans les réservoirs du véhicule, tout au long du flexible et des autres éléments reliant le dispenser aux réservoirs. L’équation de conservation de l’énergie appliquée au gaz à l’intérieur de la tuyauterie, dénommée piping, est écrite dans l’équation (25), pour chaque longueur élémentaire dx le long du piping. On fait ici l’hypothèse que le piping a une température homogène sur toute sa longueur : (25) où est le débit massique de gaz, est la capacité thermique massique du gaz, est l’abscisse le long du piping, la température du gaz dans le piping dépendant de et , est le coefficient d’échange interne dans le piping, est le diamètre interne moyen du piping et est la température du piping dépendant uniquement de . L’équation (25) peut être intégrée entre 0 (le dispenser) et la position x le long du piping pour obtenir l’expression de . Pour plus de clarté, nous notons une longueur caractéristique : (26.1) (26.2) (26.3) L’équation de conservation de l’énergie appliquée au piping est écrite dans l’équation (27) : ) (27) où est la capacité thermique du piping, est l’énergie transférée du gaz vers le piping, est le coefficient d’échange externe du piping et est la surface externe du piping. Le terme peut s’écrire, grâce à l’équation (26.3) : (28.1) (28.2) (28.3) où est la longueur du piping. L’équation (27) est discrétisée : (29) Procédé (100) de remplissage d’un réservoir (1) de gaz sous pression, mis en œuvre par un dispositif (10) de remplissage comprenant une source (2) de gaz, une conduite (3) de remplissage reliant la source (2) au réservoir (1), une vanne (4) de régulation de débit et/ou de pression dans la conduite de remplissage, et un organe (5) de commande électronique configuré pour arrêter le remplissage lorsqu’une température (ou une masse volumique) estimée du gaz présent dans le réservoir (1) atteint une valeur limite de température (ou respectivement une valeur limite de masse volumique), le procédé comprenant, avant le remplissage du réservoir (1), une étape (S1) de détermination de la température ambiante (Tamb) au niveau du dispositif (10) de remplissage, une étape (S2) de détermination de la pression (Pini) du gaz présent dans le réservoir (1) et une étape (S3) préliminaire d’estimation de la température initiale (Tini) du gaz présent dans le réservoir (1), la température initiale (Tini) du gaz présent dans le réservoir (1) étant une valeur estimée en fonction de la température ambiante (Tamb) et en fonction de la pression (Pini) du gaz présent dans le réservoir (1) avant remplissage, ladite température initiale (Tini) du gaz présent dans le réservoir (1) étant supérieure ou égale ou inférieure ou égale à la température ambiante (Tamb). Procédé selon la revendication 1, caractérisé en ce que la température initiale (Tini) du gaz présent dans le réservoir (1) est l’une au moins parmi: une première température initiale calculée (Tini_1) correspondant à un état du réservoir (1) considéré comme récemment rempli à une première masse volumique initiale (ρini_1), une seconde température initiale calculée (Tini_2) correspondant à un état du réservoir (1) considéré comme récemment soutiré jusqu’à une seconde masse volumique initiale (ρini_2). Procédé selon la revendication 2, caractérisé en ce que la première température initiale (Tini_1) du gaz présent dans le réservoir (1) est comprise dans une plage de températures hautes ayant une borne inférieure (TLgas_max_1) déterminée supérieure ou égale à la température ambiante (Tamb) et une borne supérieure (TLgas_max_2) correspondant à un maximum de limite de température déterminé, par exemple égal à 85°C. Procédé selon l'une quelconque des revendications 2 ou 3, caractérisé en ce que la première température initiale (Tini_1) est déterminée à partir d’une première courbe prédéterminée prédictive “TPgas_hot” fournie par un modèle physique prédéterminé simulant un remplissage de référence du réservoir (1) sur la plage de températures hautes. Procédé selon l’une quelconque des revendications 2 à 4, caractérisé en ce que la seconde température initiale (Tini_2) du gaz présent dans le réservoir (1) est comprise dans une plage de températures basses ayant une borne supérieure (TLgas_max_2) déterminée inférieure ou égale à la température ambiante (Tamb) et une borne inférieure (TLgas_min_2) correspondant à un minimum de limite de température déterminé, par exemple compris entre zéro et -5°C. Procédé selon l’une quelconque des revendications 2 à 5, caractérisé en ce que la seconde température initiale (Tini_2) est déterminée à partir d’une seconde courbe prédéterminée prédictive “TPgas_cold” fournie par un modèle physique prédéterminé simulant une vidange de référence du réservoir (1) sur la plage de températures basses. Procédé selon l’une quelconque des revendications 2 à 6, comprenant une étape (S4) consistant à modéliser, lors du remplissage, une première courbe de variation de la température (ou respectivement une première courbe de variation de la masse volumique) du gaz présent dans le réservoir (1), ayant comme condition de départ la première température initiale (Tini_1) (respectivement la première masse volumique initiale (ρini_1 )) du gaz présent dans le réservoir (1), et/ou une deuxième courbe de variation de la température (respectivement une deuxième courbe de variation de la masse volumique) du gaz présent dans le réservoir (1), ayant comme condition de départ la seconde température initiale ( Tini_2) (respectivement la seconde masse volumique initiale (ρini_2)) du gaz présent dans le réservoir (1). Procédé selon la revendication 7, comprenant une étape (S5) consistant à estimer, lors du remplissage, une courbe de variation de la température (ou respectivement une courbe de variation de la masse volumique) du gaz présent dans le réservoir (1), ladite courbe étant comprise entre la première courbe de variation de la température (respectivement la première courbe de variation de la masse volumique) et la deuxième courbe de variation de la température (respectivement la deuxième courbe de variation de la masse volumique). Procédé selon l’une quelconque des revendications 2 à 8, caractérisé en ce que la première température initiale (Tini_1) et/ou la seconde température initiale (Tini_2) est recalculée en cours de remplissage en fonction du débit et en fonction de la température du gaz présent dans la conduite de remplissage, ledit débit et ladite température étant déterminés par calcul et/ou par des capteurs au niveau du dispositif (10) de remplissage. Procédé selon l’une quelconque des revendications 1 à 9, caractérisé en ce que la valeur limite de température (respectivement la valeur limite de masse volumique) est une valeur fixe déterminée, par exemple 85°C (respectivement, par exemple 24,1 kg/m3 ou 40,2 kg/m3) ou une valeur donnée par une courbe “TPgas_max” de température de référence (respectivement par une courbe “RHO_cold” de masse volumique de référence), ladite courbe étant fournie par un modèle physique prédéterminé qui simule les conditions thermodynamiques du gaz lors d’un remplissage de référence du réservoir (1). Procédé selon l’une quelconque des revendications 4 à 10, caractérisé en ce que le modèle physique est basé sur un système d’équations comprenant au moins l’une parmi: une équation de bilan d’énergie interne appliqué au gaz présent dans le réservoir (1), une équation de bilan de masse appliqué au gaz présent dans le réservoir (1), une équation de conservation d’énergie dans une paroi du réservoir (1), une équation de continuité de flux de chaleur entre le gaz présent dans le réservoir (1) et la paroi du réservoir (1), une équation de continuité de flux de chaleur entre la paroi du réservoir (1) et l’air ambiant, et une équation de débit reliant un débit massique du dispositif (10) de remplissage à une différence de pression entre le dispositif (10) de remplissage et le réservoir (1). Procédé selon la revendication 11, caractérisé en ce que la première température initiale (Tini_1) (respectivement la seconde température initiale (Tini_2)) et la pression initiale (Pini) du gaz présent dans le réservoir (1) sont obtenues par résolution dudit système d’équations, par exemple à l’aide d’un algorithme de type matriciel tridiagonal. Procédé selon l’une quelconque des revendications 1 à 12, caractérisé en ce que la température ambiante (Tamb) du dispositif (10) de remplissage et la pression initiale (Pini) du gaz présent dans le réservoir (1) sont déterminées par calcul et/ou mesurées par des capteurs au niveau du dispositif (10) de remplissage. Procédé selon l’une quelconque des revendications 1 à 13, caractérisé en ce que l’organe (5) de commande électronique est configuré pour contrôler la vanne (4) de régulation de débit et/ou de pression pour réaliser une courbe ou une rampe de pression prédéterminée lors du remplissage. Procédé selon l’une quelconque des revendications 7 à 14, caractérisé en ce que l’organe (5) de commande électronique est configuré pour simuler et estimer la courbe de variation de la température et/ou la courbe de variation de la masse volumique du gaz présent dans le réservoir (1) de façon dynamique pendant le remplissage du réservoir (1) et/ou de façon anticipée, c’est-à-dire avant le remplissage. Dispositif (10) de remplissage d’un réservoir (1) de gaz sous pression, comprenant une source (2) de gaz, une conduite (3) de remplissage reliant la source (2) au réservoir (1), une vanne (3) de régulation de débit et/ou de pression dans la conduite (3) de remplissage, un ensemble de capteur(s) (6a, 6b) configuré pour mesurer la pression dans le réservoir (1) et/ou la température ambiante au niveau du dispositif (10) de remplissage, et un organe (5) de commande électronique configuré pour: -arrêter le remplissage lorsqu’une température (respectivement une masse volumique) estimée du gaz présent dans le réservoir (1) atteint une valeur limite de température (respectivement une valeur limite de masse volumique), - estimer avant le remplissage du réservoir (1), une température initiale (Tini) du gaz présent dans le réservoir (1), caractérisé en ce que la température initiale (Tini) du gaz présent dans le réservoir (1) est une valeur estimée en fonction de la température ambiante (Tamb) et en fonction de la pression (Pini) du gaz présent dans le réservoir (1) avant remplissage, la température initiale (Tini) du gaz présent dans le réservoir (1) étant supérieure ou égale ou inférieure ou égale à la température ambiante (Tamb).