FR 3067616 A1 20181221 FR 1755596 A 20170620 Arrière-plan de l'invention L'invention concerne le domaine des systèmes de production d'hydrogène gazeux, et plus particulièrement des systèmes de production d'hydrogène gazeux de type « multi-pulse à accumulation ». De façon connue, les systèmes de production d'hydrogène gazeux de type « multi-pulse à accumulation » sont aptes à générer de l'hydrogène gazeux à partir d'une pluralité de chargements unitaires indépendants disposés dans une enceinte, et initiés (i.e. allumés) successivement en fonction de la quantité d'hydrogène gazeux nécessaire. L'hydrogène gazeux produit est stocké dans le système et délivré sur demande à un ou plusieurs équipements reliés au système, comme par exemple à une pile à combustible. Chaque chargement unitaire, lorsqu'il est initié, dégage du gaz. La réaction chimique qui se produit étant exothermique, on assiste à un échauffement du chargement unitaire initié le temps de la libération du gaz. Tout le matériau constituant le chargement unitaire ne se transforme pas toutefois en gaz lors de cette réaction chimique, et il subsiste des résidus, qui sont eux-mêmes chauds. La chaleur accumulée dans ces résidus vient alors, par conduction thermique, chauffer les chargements unitaires voisins qui n'ont pas encore été initiés. Il s'ensuit un risque d'augmentation trop importante de la température des chargements unitaires voisins non encore initiés, qui peut conduire à une initiation involontaire de ces chargements. Pour réduire ce risque, il est connu de disposer les chargements unitaires du système de production d'hydrogène gazeux dans un matériau isolant thermiquement. L'épaisseur de matériau isolant entre deux chargements unitaires permet de limiter la conduction thermique. Toutefois, le choix de l'épaisseur de matériau isolant à prévoir autour de chaque chargement unitaire s'avère difficile. En effet, si l'épaisseur de matériau isolant est choisie trop petite, le risque d'avoir une initiation involontaire d'un ou de plusieurs 35 chargements unitaires subsiste. 2 Inversement, une épaisseur de matériau isolant trop importante autour de chaque chargement unitaire pénalise le système de production d'hydrogène gazeux en termes d'encombrement et de masse, et diminue significativement son rendement massique (défini comme le ratio de la 5 masse de gaz libéré par le système sur la masse totale du système) et donc ses performances. Objet et résumé de l'invention La présente invention permet de pallier notamment les 10 inconvénients précités en proposant un procédé d'identification d'une séquence d'initiation d'une pluralité de chargements indépendants notés CH1,CH2...,CHN, N désignant un entier supérieur à 1, chaque chargement indépendant comprenant un matériau apte à générer de l'hydrogène gazeux sous sollicitation thermique, les chargements indépendants 15 CH1,CH2...,CHN étant maintenus espacés les uns des autres dans une enceinte par un bloc de maintien en un matériau isolant thermiquement. Ce procédé est remarquable en ce qu'il comprend : - une étape de détermination d'une pluralité de séquences d'initiation des chargements CH1,CH2...,CHN, chaque séquence d'initiation 20 fournissant un ordre d'initiation différent des chargements CH1,CH2...,CHN ; - pour chaque séquence d'initiation déterminée : o une étape de calcul, pour chaque chargement CHn, n désignant un entier compris entre 1 et N, d'un échauffement de ce 25 chargement CHn dû à une initiation avant l'initiation du chargement CHn d'au moins un autre chargement CHk adjacent au chargement CHn, k désignant un entier différent de n; o une étape d'évaluation, pour ladite séquence d'initiation, de la valeur maximale des échauffements calculés pour les 30 chargements CH1,CH2,...,CHN ; et une étape de sélection d'une séquence d'initiation parmi ladite pluralité de séquences d'initiation pour laquelle la valeur maximale évaluée est égale à la plus petite des valeurs maximales évaluées pour ladite pluralité de séquences d'initiation ou est inférieure ou égale à un seuil d'échauffement prédéterminé. 3 Corrélativement, elle vise également un dispositif d'identification d'une séquence d'initiation d'une pluralité de chargements indépendants notés CH1,CH2...,CHN, N désignant un entier supérieur à 1, chaque chargement indépendant comprenant un matériau apte à générer de 5 l'hydrogène gazeux sous sollicitation thermique, les chargements indépendants CH1,CH2...,CHN étant maintenus espacés les uns des autres dans une enceinte par un bloc de maintien en un matériau isolant thernniquement. Ce dispositif est remarquable en ce qu'il comprend : un module de détermination configuré pour déterminer une pluralité de 10 séquences d'initiation des chargements CH1,CH2...,CHN, chaque séquence d'initiation fournissant un ordre d'initiation différent des chargements CH1,CH2...,CHN; des modules, activés pour chaque séquence d'initiation déterminée, et comprenant : 15 o un module d'estimation, configuré pour estimer pour chaque chargement CHn, n désignant un entier compris entre 1 et N, un échauffement de ce chargement CHn dû à une initiation avant l'initiation du chargement CHn d'au moins un autre chargement CHk adjacent au chargement CHn, k désignant un entier 20 différent de n; o un module d'évaluation, configuré pour évaluer pour ladite séquence d'initiation, la valeur maximale des échauffements estimés pour les chargements CH1,CH2,...,CHN ; et un module de sélection, configuré pour sélectionner une séquence 25 d'initiation parmi ladite pluralité de séquences d'initiation pour laquelle la valeur maximale évaluée est égale à la plus petite des valeurs maximales évaluées pour la pluralité de séquences d'initiation ou est inférieure ou égale à un seuil d'échauffement unitaire prédéterminé. L'invention propose ainsi une alternative particulièrement 30 avantageuse à l'état de la technique et au dimensionnement délicat de l'épaisseur d'isolant thermique nécessaire autour de chaque chargement indépendant du système de production d'hydrogène gazeux, qui consiste à déterminer une séquence d'initiation « optimale » des chargements indépendants du système qui permet de minimiser la valeur maximale 35 atteinte par les échauffements des chargements, ou tout du moins à s'assurer que cette valeur maximale reste inférieure à un seuil 4 prédéterminé. Cette séquence d'initiation « optimale » permet de réduire considérablement les risques d'une initiation involontaire d'un ou de plusieurs chargements sans pour autant sacrifier le rendement massique du système de production d'hydrogène gazeux et donc ses performances. 5 Les inventeurs ont en effet constaté que pour un même ensemble de chargements indépendants, toutes les séquences d'initiation ne sont pas équivalentes en matière d'échauffement des chargements du système de production d'hydrogène gazeux. Certaines séquences d'initiation mal choisies peuvent conduire à un échauffement très élevé 10 d'un ou de plusieurs chargements comparativement à d'autres séquences d'initiation. L'invention propose donc de minimiser les risques d'initiation involontaire de chargements indépendants dans un système de production d'hydrogène gazeux en optimisant le choix de la séquence d'initiation de 15 ces chargements plutôt que de chercher à optimiser l'épaisseur de matériau isolant présente autour de chacun des chargements. A cet effet, l'invention propose de considérer comme critère l'échauffement maximal subi par les chargements CH1,..., CHN et de sélectionner une séquence qui permet de minimiser ce critère. L'échauffement maximal subi par les 20 chargements CH1,..., CHN, s'il est trop élevé, induit un risque important d'initiation involontaire. En utilisant cet échauffement maximal comme critère de sélection, l'invention permet d'identifier de façon simple et rapide une séquence limitant un tel risque. Les risques d'initiation involontaire peuvent ainsi être facilement 25 réduits sans dégrader le rendement massique du système de production d'hydrogène gazeux ni ses performances. Dans un mode particulier de réalisation, lors de l'étape d'estimation, l'échauffement de chaque chargement CHn est estimé à partir : 30 - des instants d'initiation des chargements CHk initiés avant l'initiation du chargement CHn et adjacents au chargement CHn, lesdits instants d'initiation étant exprimé à partir de multiples d'une durée élémentaire At séparant deux initiations de chargements consécutives ; et - d'un échauffement unitaire noté Au désignant l'échauffement par 35 conduction d'un chargement résultant d'un autre chargement adjacent initié à un instant At plus tôt. 5 Ce mode de réalisation propose une méthode de calcul très simple de l'échauffement de chaque chargement indépendant pour la séquence d'initiation considérée. Cette méthode consiste à lister pour chaque chargement indépendant, les chargements adjacents à celui-ci 5 initiés avant lui et à définir les instants d'initiation de ces chargements adjacents comme des multiples d'une quantité déterminée At séparant deux initiations de chargements consécutives. En s'appuyant sur l'hypothèse qu'un chargement initié à un instant m.At plus tôt qu'un chargement donné (m désignant un nombre entier supérieur ou égal à 1) 10 conduit à un échauffement du chargement donné égal à m.Au, il est facile de calculer l'échauffement total d'un chargement résultant des initiations des chargements adjacents à ce chargement en fonction de l'échauffement unitaire Au, et de comparer l'impact des différentes séquences d'initiation sur l'échauffement de ce chargement, ou plus 15 généralement sur le système de production d'hydrogène gazeux en s'intéressant à la valeur maximale des échauffements obtenus pour les différents chargements du système. La pluralité de séquences d'initiation sur laquelle sont appliquées les étapes d'estimation et d'évaluation pour identifier une 20 séquence d'initiation optimale conformément à l'invention peut être déterminée de différentes façons. Ainsi, l'étape de détermination des séquences d'initiation peut par exemple comprendre une détermination exhaustive de toutes les séquences d'initiation possibles des N chargements. 25 Une telle détermination exhaustive conduit à considérer toutes les séquences d'initiation possibles des N chargements pour déterminer la séquence d'initiation optimale : il en existe N! (i.e. factorielle N). De cette sorte, on est sûr que lors de l'étape de sélection, la séquence d'initiation qui est sélectionnée est bien une séquence optimale conduisant à une 30 valeur minimale des échauffements des chargements ou tout du moins, à une valeur inférieure à un seuil que l'on s'est fixé préalablement. Dans une variante de réalisation, l'étape de détermination tient compte pour déterminer la pluralité de séquences d'initiation d'au moins une symétrie présentée par la disposition des chargements dans 35 l'enceinte. 6 La prise en compte de symétrie(s) permet de réduire le nombre de séquences d'initiation sur lesquelles les étapes d'estimation et de d'évaluation sont appliquées sans pour autant sacrifier l'optirnalité de la séquence d'initiation sélectionnée. Cette variante peut permettre de 5 réduire d'un facteur non négligeable la complexité de mise en oeuvre de l'invention, ce qui peut être particulièrement avantageux lorsque le nombre N de chargements considérés dans le système de production d'hydrogène gazeux est important. Dans une autre variante de réalisation, les étapes de 10 détermination, d'estimation, d'évaluation et une partie au moins de l'étape de sélection sont mises en oeuvre en utilisant un algorithme de recuit simulé. De façon connue, l'algorithme de recuit simulé est un algorithme itératif qui permet de déterminer de façon rapide et globale les extrema 15 d'une fonction. En l'espèce ici, il est utilisé pour déterminer une séquence d'initiation optimale permettant de minimiser la valeur maximale des échauffements des chargements résultant de la séquence d'initiation. L'utilisation de l'algorithme de recuit simulé permet donc de déterminer très rapidement la séquence d'initiation optimale d'un ensemble de 20 chargements donné. Elle permet en outre de réduire encore davantage la complexité de mise en oeuvre de l'invention par rapport aux deux variantes précitées en limitant drastiquement les séquences d'initiation sur lesquelles sont appliquées les étapes d'estimation et de calcul. Dans cette variante de réalisation basée sur l'utilisation d'un 25 algorithme de recuit simulé, celui-ci peut être mis en oeuvre une pluralité de fois en étant initialisé à chaque fois au moyen d'une séquence initiale d'initiation des chargements distincte et choisie aléatoirement. En exécutant l'algorithme de recuit simulé plusieurs fois en changeant à chaque fois ses conditions initiales, on évite de sélectionner 30 une séquence d'initiation correspondant à un minimum local de la fonction que l'algorithme cherche à minimiser. L'algorithme de recuit simulé peut par exemple être exécuté une centaine de fois de sorte à s'assurer qu'un minimum local est bien écarté. Dans un mode particulier de réalisation, lorsqu'il existe, parmi la 35 pluralité de séquences d'initiation, plusieurs séquences d'initiation ayant une valeur maximale évaluée égale à la plus petite des valeurs maximales 7 évaluées ou égale au seuil d'échauffement prédéterminé, la séquence d'initiation sélectionnée lors de l'étape de sélection correspond à la séquence d'initiation ayant la valeur moyenne des échauffements des chargements la plus faible. 5 Ce critère permet, lorsque plusieurs séquences d'initiation conduisent à une même valeur maximale, de sélectionner une séquence d'initiation qui présente en moyenne l'impact le plus réduit en matière d'échauffement sur ses chargements. D'autres critères en sus ou en variante des critères précités 10 peuvent être appliqués pour sélectionner la séquence d'initiation lors de l'étape de sélection. Ainsi par exemple, si plusieurs séquences d'initiation correspondent à une même valeur moyenne des échauffements des chargements (typiquement la plus faible), on peut sélectionner une séquence d'initiation qui comprend au plus un nombre prédéterminé de 15 chargements ayant une valeur supérieure à cette valeur moyenne. Dans un mode particulier de réalisation, les différentes étapes du procédé d'identification sont déterminées par des instructions de programmes d'ordinateurs. En conséquence, l'invention vise aussi un programme 20 d'ordinateur sur un support d'informations, ce programme étant susceptible d'être mis en oeuvre dans un dispositif d'identification ou plus généralement dans un ordinateur, ce programme comportant des instructions adaptées à la mise en oeuvre des étapes d'un procédé d'identification tel que décrit ci-dessus. 25 Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de code intermédiaire entre code source et code objet, tel que dans une forme partiellement compilée, ou dans n'importe quelle autre forme souhaitable. 30 L'invention vise aussi un support d'informations lisible par un ordinateur, et comportant des instructions d'un programme d'ordinateur tel que mentionné ci-dessus. Le support d'informations peut être n'importe quelle entité ou dispositif capable de stocker le programme. Par exemple, le support peut 35 comporter un moyen de stockage, tel qu'une ROM, par exemple un CD ROM ou une ROM de circuit microélectronique, ou encore un moyen 8 d'enregistrement magnétique, par exemple une disquette (floppy disc) ou un disque dur. D'autre part, le support d'informations peut être un support transmissible tel qu'un signal électrique ou optique, qui peut être 5 acheminé via un câble électrique ou optique, par radio ou par d'autres moyens. Le programme selon l'invention peut être en particulier téléchargé sur un réseau de type Internet. Alternativement, le support d'informations peut être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté 10 pour exécuter ou pour être utilisé dans l'exécution du procédé en question. Comme mentionné précédemment, grâce au procédé (et au dispositif) d'identification selon l'invention, il est possible, pour un système de production d'hydrogène gazeux comprenant un ensemble de 15 chargements indépendants donné, de déterminer une séquence d'initiation des chargements de ce système permettant de limiter, à chaque initiation, l'échauffement des chargements non encore initiés du système, et par là-même, le risque d'initiation involontaire de ces chargements non encore initiés. 20 De façon très avantageuse, le procédé d'identification proposé par l'invention peut, en complément, être utilisé pour concevoir un système de production d'hydrogène gazeux ayant des performances optimisées, et ce de différentes manières. Ainsi, l'invention a également pour objet un premier procédé de 25 conception d'un système de production d'hydrogène gazeux comprenant : - une étape de détermination d'une pluralité de chargements indépendants CH1,CH2...,CHN pour le système de production d'hydrogène gazeux, N désignant un nombre entier supérieur à 1, chaque chargement indépendant comprenant un matériau apte à 30 générer de l'hydrogène gazeux sous sollicitation thermique, les chargements indépendants étant maintenus espacés les uns des autres dans une enceinte par un bloc de maintien en un matériau isolant thermiquement ; - une étape d'identification d'une séquence d'initiation des chargements 35 CH1,CH2...,CHN en utilisant un procédé d'identification selon l'invention ; 9 - une étape de sélection d'au moins un chargement parmi la pluralité de chargements CH1,...,CHN, ledit au moins un chargement sélectionné présentant l'échauffement le plus élevé pour la séquence d'initiation identifiée parmi ladite pluralité de chargements CH1,...,CHN ; et 5 - une étape de renforcement de l'isolation thermique dudit au moins un chargement sélectionné dans l'enceinte. Autrement dit, le procédé d'identification est avantageusement utilisé pour déterminer, dans un système de production d'hydrogène gazeux donné, un ou plusieurs chargements pour lesquels la séquence 10 d'initiation optimale conduit à un échauffement plus élevé que les autres chargements, et pour compenser cet échauffement plus élevé par un renforcement de l'isolation thermique du ou des chargements concernés. De cette sorte, on abaisse le niveau d'échauffement maximal subi par le système de production d'hydrogène gazeux tout en ayant une action 15 ciblée sur un nombre limité de chargements. Le renforcement local de l'isolation thermique d'un nombre réduit de chargements judicieusement choisis permet de diminuer encore davantage le risque d'initiation involontaire de chargements, sans augmenter de façon drastique la masse du système de production d'hydrogène gazeux. Les conséquences sur le 20 rendement massique du système de production d'hydrogène gazeux sont donc limitées. On note qu'aucune limitation n'est attachée à la façon dont l'isolation thermique du ou des chargements identifiés est réalisée. L'étape de renforcement de l'isolation thermique peut 25 comprendre par exemple une augmentation de l'épaisseur du matériau isolant thermiquement du bloc de maintien autour dudit au moins un chargement sélectionné. En variante, l'étape de renforcement de l'isolation thermique peut comprendre l'utilisation pour ledit au moins un chargement 30 sélectionné d'un matériau isolant thermiquement ayant de meilleures capacités d'isolation thermique que le matériau isolant utilisé pour les autres chargements. Dans une autre variante encore, on peut envisager de combiner ces deux modes de renforcement de l'isolation thermique du ou des 35 chargements identifiés. 10 Par ailleurs, le procédé d'identification selon l'invention peut être exploité d'une autre manière pour concevoir un système de production d'hydrogène gazeux. Ainsi, l'invention vise un second procédé de conception d'un 5 système de production d'hydrogène gazeux comprenant : - une étape de détermination d'une pluralité de chargements indépendants CH1,CH2...,CHN pour le système de production d'hydrogène gazeux, N désignant un nombre entier supérieur à 1, chaque chargement indépendant comprenant un matériau apte à 10 générer de l'hydrogène gazeux sous sollicitation thermique, les chargements indépendants étant maintenus espacés les uns des autres dans une enceinte par un bloc de maintien en un matériau isolant thermiquement ; - une étape d'identification d'une séquence d'initiation des chargements 15 CH1,CH2...,CHN en utilisant un procédé d'identification selon l'invention; - une étape de sélection d'au moins un chargement parmi la pluralité de chargements CH1,CH2...,CHN, ledit au moins un chargement sélectionné présentant l'échauffement le plus élevé pour la séquence 20 d'initiation identifiée parmi la pluralité de chargements CH1,...,CHN ; et - une étape de retrait de l'enceinte dudit au moins un chargement sélectionné. Selon ce procédé, on part d'une configuration initiale du système de production d'hydrogène gazeux comprenant un nombre N de 25 chargements indépendants et on retire de cette configuration initiale un ou plusieurs chargements subissant des échauffements élevés. En retirant ces chargements de la configuration initiale, il est possible d'abaisser significativement l'échauffement maximal subi par le système de production d'hydrogène gazeux sans avoir besoin de renforcer l'isolation 30 thermique des chargements du système. Le nombre de chargements indépendants à retirer peut être paramétré en fonction de l'encombrement par exemple souhaité pour le système de production d'hydrogène gazeux que l'on cherche à concevoir. En variante, le retrait des chargements identifiés peut être 35 accompagné d'une diminution de l'isolation thermique (par exemple en réduisant l'épaisseur de cette isolation) de tout ou partie des chargements 11 du système (par exemple en maintenant constant l'échauffement maximal). Une telle variante permet d'obtenir un gain de masse. Dans un mode particulier de réalisation, le deuxième procédé de conception peut comprendre, suite à l'étape de retrait, une étape 5 d'identification d'une nouvelle séquence d'initiation de la pluralité de chargements obtenue après retrait dudit au moins un chargement sélectionné en utilisant un procédé d'identification selon l'invention. Il convient de noter que suivant les contraintes imposées lors de la conception du système de production d'hydrogène gazeux ou des 10 critères recherchés pour ce système, notamment en matière de masse, d'encombrement, de durée At entre deux initiations successives, etc., l'un ou l'autre des procédés de conception proposés peut être utilisé pour concevoir un système performant comprenant un nombre déterminé de chargements. 15 Brève description des dessins D'autres caractéristiques et avantages de la présente invention ressortiront de la description faite ci-dessous, en référence aux dessins annexés qui en illustrent un exemple de réalisation dépourvu de tout 20 caractère limitatif. Sur les figures : - la figure 1. représente dans son environnement un dispositif d'identification d'une séquence d'initiation d'une pluralité de chargements d'un système de production d'hydrogène gazeux conforme à l'invention ; 25 - la figure 2 illustre l'architecture matérielle du dispositif d'identification de la figure 1, dans un mode particulier de réalisation ; - la figure 3 représente, sous forme d'ordinogramme, les principales étapes d'un procédé d'identification selon l'invention, tel qu'il est mis en oeuvre par le dispositif d'identification de la figure 1 dans un 30 premier mode de réalisation ; - la figure 4 représente, sous forme d'ordinogramme, les principales étapes d'un procédé d'identification selon l'invention, tel qu'il est mis en oeuvre par le dispositif d'identification de la figure 1 dans un deuxième mode de réalisation ; 35 - la figure 5 représente, sous forme d'ordinogramme, les principales étapes d'un premier procédé de conception d'un système de 12 production d'hydrogène gazeux selon l'invention, utilisant un procédé d'identification selon l'invention ; et la figure 6 représente, sous forme d'ordinogramme, les principales étapes d'un deuxième procédé de conception d'un système de 5 production d'hydrogène gazeux selon l'invention, utilisant un procédé d'identification selon l'invention. Description détaillée de l'invention La figure 1 représente, dans son environnement, un dispositif 1 10 d'identification conforme à l'invention, apte à identifier une séquence d'initiation (i.e. d'allumage) d'une pluralité de chargements indépendants d'un système 2 de production d'hydrogène gazeux. Dans l'exemple envisagé à la figure 1, le système 2 de production d'hydrogène gazeux est un système de type « multi-pulse à 15 accumulation », apte à générer de l'hydrogène gazeux à partir d'une pluralité de chargements unitaires indépendants CH1, CH2,..., CHN, N désignant un nombre entier supérieur à 1 (par exemple N=14 sur la figure 1). Ces chargements unitaires CH1, CH2,..., CHN sont disposés axialement dans une enceinte 20, et maintenus espacés les uns des autres dans des 20 logements prévus à cet effet dans un bloc de maintien 21 cylindrique, réalisé en un matériau isolant thernniquement. Chaque chargement unitaire CHn, n=1,...,N renferme ici une pluralité de segments annulaires empilés les uns sur les autres (non représentés), et constitués d'un matériau apte à générer de l'hydrogène 25 gazeux sous sollicitation thermique, comme par exemple un matériau à base d'ammoniaborane NH3BH3. Chaque chargement unitaire CHn est muni d'un initiateur relié à un système de commande 3 permettant d'allumer les segments du chargement unitaire et de provoquer leur décomposition thermique. 30 En variante, chaque chargement unitaire peut renfermer un seul bloc monolithique de matériau apte à générer de l'hydrogène gazeux sous sollicitation thermique et dont la décomposition est initiée au moyen d'un initiateur. Comme indiqué précédemment, le système de commande 3 est 35 relié à chaque initiateur des chargements unitaires CH1, CH2,..., CHN afin de commander l'allumage séquentiel de ceux-ci, typiquement en fonction 13 de la quantité d'hydrogène gazeux nécessaire. L'hydrogène gazeux produit est stocké dans le système 2 de production d'hydrogène gazeux et délivré sur demande à un ou plusieurs équipements reliés au système (non représentés), comme par exemple à une pile à combustible. 5 Chaque chargement unitaire CHn, n=1,...,N, lorsqu'il est initié par le système de commande 3, dégage du gaz. La réaction chimique qui se produit étant exothermique, on assiste à un échauffement du chargement unitaire initié le temps de la libération du gaz. Tout le matériau constituant le chargement unitaire ne se transformant toutefois 10 pas en gaz lors de cette réaction chimique, il subsiste des résidus de matériau dans le chargement unitaire eux-mêmes chauds. La chaleur accumulée dans ces résidus vient alors, par conduction thermique, chauffer les chargements unitaires voisins (i.e. adjacents) qui n'ont pas encore été initiés. Il s'ensuit une augmentation de la température (i.e. un 15 échauffement) des chargements unitaires voisins non encore initiés, qui peut conduire à une initiation involontaire de ces chargements. Le dispositif 1 d'identification est configuré, conformément à l'invention, pour déterminer une séquence d'initiation SEQopt des chargements unitaires CH1, CH2,..., CHN (c'est-à-dire un ordre d'initiation des chargements 20 unitaires) qui permet de réduire l'échauffement maximal susceptible d'être subi par les chargements du système 2 de production d'hydrogène gazeux. Dans le mode de réalisation décrit ici, le dispositif 1 d'identification a l'architecture matérielle d'un ordinateur, telle qu'illustrée 25 à la figure 2. Il comprend notamment un processeur 4, une mémoire vive 5, une mémoire morte 6, une mémoire flash non volatile 7, des interfaces d'entrée/sortie 8 (comme par exemple un écran, un clavier, etc.), ainsi que des moyens de communication 9. Ces moyens de communication 9 30 permettent au dispositif 1 d'identification de communiquer par exemple via un réseau de télécommunications avec le système de commande 3 pour lui transmettre la séquence d'initiation SEQopt qu'il a identifiée pour les chargements CH1, CH2,..., CHN du système 2 de production d'hydrogène gazeux. 14 En variante, le système de commande 3 est paramétré par un opérateur à partir de la séquence d'initiation SEQopt fournie par le dispositif 1 d'identification. La mémoire morte 6 du dispositif 1 d'identification constitue un 5 support d'enregistrement conforme à l'invention, lisible par le processeur 4 et sur lequel est enregistré ici un programme d'ordinateur PROG conforme à l'invention. Le programme d'ordinateur PROG définit des modules fonctionnels (et logiciels ici), configurés pour mettre en oeuvre les étapes 10 du procédé d'identification selon l'invention. Les modules fonctionnels définis par le programme PROG s'appuient sur et/ou commandent les éléments matériels 4-9 du dispositif 1 d'identification mentionnés précédemment. Ils comprennent notamment ici, comme illustré sur la figure 1 : 15 - un module de détermination 1A, configuré pour déterminer une pluralité de séquences d'initiation des chargements indépendants CH1,CH2...,CHN. Chaque séquence d'initiation notée SEQi, i=1,...,I où I désigne un entier supérieur à 1, fournit un ordre d'initiation différent des chargements indépendants CH1,CH2...,CHN ; 20 - des modules, activés pour chaque séquence d'initiation SEQi, i=1,...,I déterminée par le module de détermination 1A, et comprenant : o un module de calcul 1B, configuré pour calculer pour chaque chargement indépendant CHn, un échauffement noté AU(CHn,SEQ0 de ce chargement indépendant CHn dû à une 25 initiation avant l'initiation du chargement indépendant CHn d'au moins un autre chargement indépendant CHk adjacent au chargement indépendant CHn, k désignant un entier différent de n; o un module d'évaluation 1C, configuré pour évaluer pour la 30 séquence d'initiation SEQi considérée, la valeur maximale AUmax(SEQD des échauffements calculés pour les chargements indépendants CH1,CH2,...,CHN ; et un module de sélection 1D, configuré pour sélectionner une séquence d'initiation SEQopt parmi ladite pluralité de séquences d'initiation 35 SEQ1,...,SEQI pour laquelle la valeur maximale déterminée AUmax(SEQi) est égale à la plus petite des valeurs maximales 15 déterminées pour ladite pluralité de séquences d'initiation ou qui est inférieure à un seuil d'échauffement THR prédéterminé. Les différentes fonctions des modules 1A-1D précités sont décrites maintenant plus en détail en référence aux étapes du procédé 5 d'identification selon l'invention, dans plusieurs modes de réalisation. La figure 3 représente les différentes étapes du procédé d'identification selon l'invention, telles qu'elles sont mises en oeuvre par le dispositif 1 d'identification pour identifier une séquence d'initiation des chargements CH1,..., CHN, dans un premier mode de réalisation. 10 Conformément à l'invention, le dispositif 1 d'identification, par le biais de son module de détermination 1A, détermine une pluralité de séquences d'initiation SEQ1,...,SEQI, des chargements CH1,..., CHN, I désignant un entier supérieur à 1 (étape E10). Chaque séquence d'initiation SEQi, i=1,...,I définit un ordre différent d'initiation des 15 chargements CH1,..., CHN. Dans le premier mode de réalisation décrit ici, par souci de simplification, on suppose que deux initiations successives de chargements définies dans les séquences d'initiation SEQi sont séparées par une durée élémentaire notée At, .8.t désignant un nombre réel strictement positif. 20 A titre illustratif, pour la configuration à 14 chargements illustrée à la figure 1, deux exemples de séquences d'initiation possibles peuvent être - une première séquence SEQ1={CH1,CH2,CH3,...,CH14}, où les chargements sont initiés dans l'ordre croissant de leur indice tous les 25 At, autrement dit le chargement CH1 est initié en premier à un instant tO>0, puis le chargement CH2 à tO+At, puis le chargement CH3 à t0+2& et ainsi de suite jusqu'au chargement CH14 initié à At+13At ; - une deuxième séquence SEQ2={CH1, CH2, CH4, CH3, CH5, CH7, CH8, CH6, CH9, CH1, CH10, CH12, CH13, CH14}, selon laquelle le 30 chargement CH1 est initié en premier, suivi du chargement CH2, du chargement CH4, du chargement CH3, et ainsi de suite jusqu'au chargement CH14; et - une troisième séquence SEQ3={CH1, CH3, CH2, CH4, CH7, CH5, CH6, CH8, CH1, CH9, CH10, CH12, CH13, CH14} ; 35 - etc. Bien entendu ces exemples ne sont donnés qu'à titre illustratif. 16 Différentes stratégies peuvent être adoptées par le module de détermination 1A pour déterminer la pluralité de séquences d'initiation SEQ1,..., SEQI. Ainsi, si le nombre de chargements N compris dans le système 2 5 de production d'hydrogène gazeux le permet, le module de détermination lA peut déterminer de façon exhaustive toutes les séquences d'initiation possibles des N chargements. Il en existe I=N! (i.e. factorielle N, soit Nx(N-1)x(N-2)x...x2x1), soit dans l'exemple illustré à la figure 1 où N=14, 1=87 178 291 200 séquences. 10 Dans une variante de réalisation, pour réduire le nombre de séquences d'initiation SEQ1,..., SEQI à considérer, le module de détermination 1A peut tenir compte de symétrie(s) présentée(s) par la disposition des chargements CH1, ..., CHN au sein de l'enceinte 20. Ainsi, la configuration de chargements illustrée à la figure 1 présente diverses 15 symétries qui rendent équivalentes, en matière d'échauffement provoqué sur les chargements du système 2, les séquences d'initiation commençant par le chargement CH1 et les séquences d'initiation commençant par les chargements CH2, CH13 ou encore CH14. Plus généralement au regard de cette configuration et des symétries qu'elle présente, il est suffisant pour 20 le module de détermination 1A de considérer les séquences d'initiation commençant respectivement par les chargements CH1, CH3, CH4, CH6 et CH7 pour en déduire l'effet en terme d'échauffement de toutes les séquences d'initiation commençant par les autres chargements du système 2 de production d'hydrogène gazeux. Ceci permet de réduire le 25 nombre de séquences d'initiation I de N!=14! à I=(5/14).14! soit I=5.13!=31 135 104 000 séquences. Dans un deuxième mode de réalisation décrit plus en détail ultérieurement en référence à la figure 4, les séquences d'initiation déterminées par le module de détermination 1A peuvent répondre à une 30 logique de détermination appliquée par un algorithme d'optimisation, tel que par exemple un algorithme de recuit simulé. Dans ce cas, le nombre I de séquences d'initiation déterminées par le module d'initiation 1 peut être soit prédéterminé (et par exemple lié à un nombre maximal d'itérations mis en oeuvre par l'algorithme d'optimisation), soit dépendre de la 35 convergence de l'algorithme. 17 Pour chacune des séquences d'initiation SEQi, i=1,...,I ainsi déterminées par le module de détermination 1A (boucle représentée par les étapes E20 et E80), le dispositif 1 d'identification, par le biais de son module d'estimation 1B, estime pour chaque chargement unitaire CHn, 5 n=1,...,N, (étapes d'initialisation E40, de test E60 et d'itération E70) l'échauffement AU(CHn,SEQi) de ce chargement dû le cas échéant à une initiation antérieure des chargements CHk, k*n, qui lui sont adjacents (étape E50). Autrement dit, pour chaque chargement unitaire CHn, n=1,...,N, 10 le module d'estimation 1B détermine si un ou plusieurs chargements adjacents (voisin) à CHk a ou ont été initié(s), et le cas échéant, estime l'échauffement du chargement CHn résultant de cette ou ces initiations antérieures. Pour estimer l'échauffement AU(CHn,SEQi) du chargement CHn 15 pour la séquence d'initiation SEQi, le module d'estimation 1B procède, dans le premier mode de réalisation décrit ici, de la façon suivante. Il détermine tout d'abord, au regard de la séquence d'initiation SEQi considérée, tous les chargements adjacents au chargement CHn initiés avant le chargement CHn. A titre illustratif, dans l'exemple de la 20 configuration de chargements représentés à la figure 1 et pour la séquence d'initiation SEQ2={CH1, CH2, CH4, CH3, CH5, CH7, CH8, CH6, CH9, CH11, CH10, CH12, CH13, CH14}, le module d'estimation 1B identifie pour le chargement CH7, les chargements voisins CH3 et CH4 initiés avant le chargement CH7. Les chargements CH1, CH2 et CH5 ont également été 25 initiés avant le chargement CH7, mais ils ne sont pas adjacents au chargement CH7 dans la configuration illustrée à la figure 1. Les chargements CH6, CH8, CH10 et CH11 sont adjacents au chargement C7, mais ils sont initiés selon la séquence d'initiation SEQ2 après le chargement C7. Ils ne sont pas pris en compte par le module d'estimation 30 1B. Chaque chargement adjacent au chargement CHn et initiés avant le chargement CHn ainsi identifié par le module d'estimation 1B est associé à un instant d'initiation de ce chargement adjacent, qui peut être exprimé sous la forme d'un multiple de la durée élémentaire At séparant 35 deux initiations de chargements consécutives. Dans le premier mode de réalisation décrit ici, l'instant d'initiation du chargement adjacent considéré 18 est défini de façon relative par rapport à l'initiation du chargement CHn, autrement dit on définit cet instant d'initiation par rapport à l'instant d'initiation du chargement CHn. Ainsi, dans l'exemple illustratif considéré ci-avant, le chargement adjacent CH3 a été initié 2.At avant le chargement 5 CH7, de sorte que le module d'estimation 1B lui associe un instant d'initiation relatif égal à 2.1t. De façon similaire, le chargement CH4 a été initié 3.At avant le chargement CH7. Le module d'estimation 1B déduit de ces instants d'initiation relatifs des chargements CH3 et CH4 adjacents au chargement CH7, 10 l'échauffement résultant sur le chargement CH7 dû aux initiations antérieures des chargements CH3 et CH4 en s'appuyant sur l'hypothèse qu'un chargement initié à un instant mat plus tôt que le chargement CH7 provoque, m désignant un entier supérieur ou égal à 1, par conduction thermique, sur le chargement CH7 un échauffement égal à m.AU où AU 15 désigne un échauffement unitaire désignant l'échauffement par conduction d'un chargement unitaire résultant d'un autre chargement unitaire adjacent initié à un instant At plus tôt. En d'autres mots, dans l'exemple illustratif envisagé ci-avant : - le chargement CH3 initié 2.At avant le chargement CH7 provoque sur 20 le chargement CH7 un échauffement égal à 2.AU ; et - le chargement CH4 initié 3.At avant le chargement CH7 provoque sur le chargement CH7 un échauffement égal à 3.AU. Il en résulte un échauffement total pour le chargement CH7 et la séquence d'initiation SEQ2 égal à : 25 AU(CH7,SEQ2)=5.1U De façon générale, pour chaque chargement CHn, le module d'estimation 1B estime l'échauffement de ce chargement CHn pour la séquence SEQi considérée selon l'équation (Eq. 1) suivante : AU(CHn, SEQi) - Eike n,CHk adjacent à CHn et initié mkAt avant CHn selon SEQi Mk- AU 30 (Eq. 1) Le module d'évaluation 1C du dispositif 1 d'identification évalue à partir des échauffements AU(CHn,SEQi), n=1,...,N estimés par le module d'estimation 1B la valeur maximale AUmax(SEQi) de ces échauffements donnée par l'équation (Eq. 2) (étape E80) : 35 AUmax(SEQ0 = inaXn=1,...,N A Umax(CHn, SEQi) (Eq. 2) 19 Il réitère les étapes de détermination E30-E60 et d'évaluation E80 pour chaque séquence d'initiation SEQi, i=1,...,I (étapes d'itération E90 et de test E30). Une fois toutes les séquences d'initiation SEQi considérées 5 (réponse non à l'étape test E30), le module de sélection 1D du dispositif 1 d'identification sélectionne parmi les séquences SEQ1,...,SEQI une séquence d'initiation SEQopt des chargements CH1,...,CHN (étape E100). Conformément à l'invention, il sélectionne une séquence d'initiation parmi ladite pluralité de séquences d'initiation SEQ1,...,SEQI 10 pour laquelle la valeur maximale AUmax(SEQi) déterminée est égale à la plus petite des valeurs maximales déterminées pour la pluralité de séquences d'initiation soit : SEQopt = argmini=1,.. JAUmax(SEQO (Eq. 3) En variante, le module de sélection 1D sélectionne une 15 séquence SEQopt dont la valeur maximale AUmax(SEQi) est inférieure (ou égale) à un seuil d'échauffement prédéterminé THR. Par exemple, s'il est connu qu'une configuration de système de production d'hydrogène gazeux comprenant N=14 cartouches ne peut pas conduire à un échauffement maximal inférieur à 7.AU quelle que soit la séquence d'initiation 20 considérée, un seuil THR pouvant être utilisé par le module de sélection 1D est 7.AU. Si plusieurs séquences SEQi parmi la pluralité de séquences SEQ1,...,SEQI répondent au critère de sélection appliqué par le module de sélection ID pour sélectionner la séquence d'initiation SEQopt, le module 25 de sélection 1D peut appliquer en complément d'autres critères. Par exemple, il peut sélectionner parmi les séquences SEQi présélectionnées (autrement dit répondant au critère de valeur minimale des valeurs maximales ou de valeur maximale inférieure au seuil THR), une séquence d'initiation correspondant à la valeur moyenne des échauffements AU 30 (CHn,SEQi) des chargements CH1,...,CHN la plus faible. Bien entendu d'autres critères complémentaires peuvent être appliqués en variante. La séquence d'initiation SEQopt ainsi sélectionnée peut ensuite être appliquée par le système de commande 3 aux chargements CH1, ..., 35 CHN du système 2 de production d'hydrogène gazeux. 20 Dans le premier mode de réalisation décrit ici, le dispositif 1 d'identification considère de façon exhaustive, pour sélectionner la séquence SEQopt, toutes les séquences d'initiation SEQ1,...,SEQI possibles des chargements CH1,...,CHN, ou un nombre réduit de séquences obtenu 5 en tenant compte de symétries existant dans la disposition des chargements au sein de l'enceinte 20. On note que les estimations et évaluations mises en oeuvre sur les différentes séquences d'initiation envisagées peuvent être réalisées séquentiellement ou en parallèle. Dans un deuxième mode de réalisation décrit maintenant en 10 référence à la figure 4, le dispositif 1 d'identification utilise, pour mettre en oeuvre le procédé d'identification selon l'invention, un algorithme d'optimisation, et plus particulièrement un algorithme de recuit simulé. Cet algorithme permet de réduire drastiquement le nombre de séquences d'initiation considérées pour identifier la séquence « optimale » SEQopt et 15 donc d'accélérer le procédé d'identification selon l'invention. De façon connue, l'algorithme de recuit simulé est inspiré d'un processus utilisé en métallurgie dans lequel on alterne des cycles de refroidissement lent et de réchauffage (« recuit ») pour minimiser l'énergie d'un matériau ou d'un système et atteindre une configuration 20 stable des atomes le constituant. Le principe général appliqué par l'algorithme de recuit simulé est le suivant. Par analogie avec le processus utilisé en métallurgie, l'algorithme cherche à minimiser une fonction E déterminée assimilée à l'énergie du système. Pour ce faire, l'algorithme comprend une pluralité 25 d'itérations j=1,...,J, J désignant un entier supérieur à 1, au cours desquelles il part d'un état s du système, modifie cet état en lui appliquant une modification élémentaire, puis détermine si le nouvel état obtenu sj améliore ou pas la fonction qu'il cherche à minimiser (en calculant E(sj)) : si le nouvel état obtenu sj résulte en une variation négative AE de la 30 fonction (i.e. la fonction diminue avec le nouvel état), il est conservé et considéré comme état de départ s=sj à l'itération suivante j+1. Sinon, le nouvel état sj est conservé et utilisé comme état de départ à l'itération suivante j+1 avec une probabilité exp(-AE/T) où T désigne un paramètre fictif assimilé à la température du système. 35 Suivant l'approche retenue, le paramètre T peut être gardé constant (égal à une température initiale TO) sur l'ensemble des itérations 21 ou varier en fonction des itérations (ex. baisser de façon continue selon une loi de décroissance prédéterminée). On note qu'accepter un nouvel état améliorant la fonction revient à chercher l'optimum dans le voisinage de l'état précédent. L'acceptation d'un « mauvais » état permet d'explorer 5 une plus grande partie de l'espace des états et tend à éviter la recherche d'un optimum local. Le nombre d'itérations 3 mises en oeuvre par l'algorithme peut être prédéterminé, ou dépendre d'un critère de convergence testé à chaque itération. 10 Dans le deuxième mode de réalisation décrit ici, le dispositif 1 d'identification applique un tel algorithme de recuit simulé pour minimiser la fonction AUmax(SEQ0. Cette fonction est calculée, comme dans le premier mode de réalisation décrit précédemment, au moyen des équations (Eq. 1) et (Eq. 2) définissant respectivement le modèle 15 d'échauffement de chaque chargement CHn (cf. Eq. 1), et l'échauffement maximal affectant les chargements CH1,..., CHN pour la séquence d'initiation SEQi considérée (cf. Eq. 2). Dans le deuxième mode de réalisation décrit ici, l'algorithme de recuit simulé est appliqué par le dispositif 1 d'identification une pluralité de 20 fois (i.e. Q fois, Q désignant un entier supérieur à 1, par exemple Q=100) au cours desquelles certaines des conditions initiales de l'algorithme sont modifiées (étape d'initialisation F10, de test F50 et d'incrémentation F60). Cette application une pluralité de fois vise à éviter à l'algorithme de recuit simulé de converger vers un minimum local. En changeant les conditions 25 initiales de l'algorithme, on peut obtenir à chaque application une séquence d'initiation « optimale » puis comparer ces séquences d'initiations « optimales » afin de sélectionner celle qui conduit à l'échauffement maximal des chargements le plus faible, comme décrit plus en détail ultérieurement. 30 Dans le deuxième mode de réalisation décrit ici, les conditions initiales de l'algorithme de recuit simulé appliqué par le dispositif 1 d'identification comprennent : - une séquence d'initiation initiale SEQinit (utilisé comme état initial du système à la première itération de l'algorithme) ; 35 - un nombre d'itérations 3 mis en oeuvre à chaque application de l'algorithme ; et 22 - un intervalle de variation du paramètre de température [Tmax ;Trnin] (on envisage par exemple une loi de décroissance de la température à chaque itération de la forme logarithmique Tj=Tmax exp(-(j/J) In(Tmax/Tmin) où Tmax et Tmin désignent respectivement 5 une température maximale et une température minimale données par l'intervalle de variation du paramètre de température). Dans l'exemple envisagé ici, la séquence d'initiation initiale est modifiée à chaque application de l'algorithme (étape F20). Elle est choisie de façon aléatoire à chaque application. Par ailleurs, on définit comme 10 intervalle de variation de la température T, l'intervalle [150°C; 0.05°C] (autrement dit T0=150°C), un pas de décroissance pour la température X=0.3.10-6 et le nombre d'itérations J est fixé à 5000. Bien entendu, ces conditions initiales ne sont données qu'à titre illustratif et d'autres conditions initiales peuvent être envisagées. Par 15 exemple, J peut ne pas être fixé et dépendre d'un critère de convergence de l'algorithme de recuit simulé évalué à chaque itération pour déterminer si l'algorithme est stoppé ou si une nouvelle itération est mise en oeuvre. Pour chaque valeur de l'indice q (q=1,...,Q), l'algorithme de recuit simulé est appliqué par le dispositif 1 d'identification et exécute lors 20 de cette application J=5000 itérations (étape F30). Au cours de chaque itération indexée par j (étape F31 d'initialisation de l'indice j, F38 de test que j REVENDICATIONS1. Procédé d'identification d'une séquence d'initiation (SEQopt) d'une pluralité de chargements indépendants notés CH1,CH2...,CHN, N 5 désignant un entier supérieur à 1, chaque chargement indépendant comprenant un matériau apte à générer de l'hydrogène gazeux sous sollicitation thermique, les chargements indépendants CH1,CH2...,CHN étant maintenus espacés les uns des autres dans une enceinte (20) par un bloc de maintien en un matériau isolant thermiquement, ledit procédé 10 étant caractérisé en ce qu'il comprend : - une étape de détermination (E10,F32) d'une pluralité de séquences d'initiation des chargements CH1,CH2...,CHN, chaque séquence d'initiation fournissant un ordre d'initiation différent des chargements CH1,CH2...,CHN; 15 - pour chaque séquence d'initiation déterminée : o une étape d'estimation (E50,F33), pour chaque chargement CHn, n désignant un entier compris entre 1 et N, d'un échauffement de ce chargement CHn dû à une initiation avant l'initiation du chargement CHn d'au moins un autre chargement 20 CHk adjacent au chargement CHn, k désignant un entier différent de n; o une étape d'évaluation (E80,F34), pour ladite séquence d'initiation, de la valeur maximale des échauffements estimés pour les chargements CH1,CH2,...,CHN ; et 25 - une étape de sélection (E100, F70) d'une séquence d'initiation parmi ladite pluralité de séquences d'initiation pour laquelle la valeur maximale évaluée est égale à la plus petite des valeurs maximales évaluées pour ladite pluralité de séquences d'initiation ou qui est inférieure ou égale à un seuil d'échauffement prédéterminé. 30 2. Procédé d'identification selon la revendication 1 dans lequel lors de l'étape d'estimation, l'échauffement de chaque chargement CHn est estimé à partir : - des instants d'initiation des chargements CHk initiés avant l'initiation 35 du chargement CHn et adjacents au chargement CHn, lesdits instants 30 d'initiation étant des multiples d'une durée élémentaire At séparant deux initiations de chargements consécutives ; et d'un échauffement unitaire noté Au désignant l'échauffement par conduction d'un chargement résultant d'un autre chargement adjacent initié à un instant At plus tôt. 3. Procédé d'identification selon la revendication 1 ou 2 dans lequel l'étape de détermination (E10) comprend une détermination exhaustive de toutes les séquences d'initiation possibles des N 10 chargements. 4. Procédé d'identification selon la revendication 1 ou 2 dans lequel l'étape de détermination (E10) tient compte, pour déterminer la pluralité de séquences d'initiation, d'au moins une symétrie présentée par 15 la disposition des chargements dans l'enceinte. 5. Procédé d'identification selon la revendication 1 ou 2 dans lequel les étapes de détermination (F32), d'estimation (F33), d'évaluation (F34) et une partie au moins de l'étape de sélection (F40) sont mises en 20 oeuvre en utilisant un algorithme de recuit simulé. 6. Procédé d'identification selon la revendication 5 dans lequel l'algorithme de recuit simulé est mis en oeuvre une pluralité de fois, ledit algorithme de recuit simulé étant initialisé à chaque fois au moyen d'une 25 séquence initiale d'initiation des chargements distincte choisie aléatoirement. 7. Procédé d'identification selon l'une quelconque des revendications 1 à 6 dans lequel, lorsqu'il existe, parmi la pluralité de 30 séquences d'initiation, plusieurs séquences d'initiation ayant une valeur maximale estimée égale à la plus petite des valeurs maximales estimées ou égale au seuil d'échauffement prédéterminé, la séquence d'initiation sélectionnée lors de l'étape de sélection correspond à la séquence d'initiation ayant la valeur moyenne des échauffements des chargements 35 la plus faible. 31 8. Programme d'ordinateur (PROG) comportant des instructions pour l'exécution des étapes du procédé d'identification selon l'une quelconque des revendications 1 à 7 lorsque ledit programme est exécuté par un ordinateur. 5 9. Support d'enregistrement (6) lisible par un ordinateur sur lequel est enregistré un programme d'ordinateur (PROG) selon la revendication 8. 10. Procédé de conception d'un système (2) de production d'hydrogène gazeux, ledit procédé comprenant : - une étape de détermination (G10) d'une pluralité de chargements indépendants CH1,CH2...,CHN pour le système de production d'hydrogène gazeux, N désignant un nombre entier supérieur à 1, chaque chargement indépendant comprenant un matériau apte à générer de l'hydrogène gazeux sous sollicitation thermique, les chargements indépendants étant maintenus espacés les uns des autres dans une enceinte par un bloc de maintien en un matériau isolant thermiquement ; - une étape d'identification (G20) d'une séquence d'initiation des chargements CH1,CH2...,CHN en utilisant un procédé d'identification selon l'une quelconque des revendications 1 à 7; - une étape de sélection (G30) d'au moins un chargement parmi la pluralité de chargements CH1,...,CHN, ledit au moins un chargement sélectionné présentant l'échauffement le plus élevé pour la séquence d'initiation identifiée parmi ladite pluralité de chargements CH1,...,CHN; et - une étape de renforcement (G40) de l'isolation thermique dudit au moins un chargement sélectionné dans l'enceinte. 11. Procédé de conception selon la revendication 10 dans lequel l'étape de renforcement de l'isolation thermique comprend une augmentation de l'épaisseur du matériau isolant thermiquement du bloc de maintien autour dudit au moins un chargement sélectionné. 32 12. Procédé de conception selon la revendication 10 dans lequel l'étape de renforcement de l'isolation thermique comprend l'utilisation pour ledit au moins un chargement sélectionné d'un matériau isolant thermiquement ayant de meilleures capacités d'isolation thermique que le 5 matériau isolant utilisé pour les autres chargements. 13. Procédé de conception d'un système (2) de production d'hydrogène gazeux, ledit procédé comprenant : - une étape de détermination (H10) d'une pluralité de chargements 10 indépendants CH1,CH2...,CHN pour le système de production d'hydrogène gazeux, N désignant un nombre entier supérieur à 1, chaque chargement indépendant comprenant un matériau apte à générer de l'hydrogène gazeux sous sollicitation thermique, les chargements indépendants étant maintenus espacés les uns des autres 15 dans une enceinte par un bloc de maintien en un matériau isolant thermiquement ; une étape d'identification (H20) d'une séquence d'initiation des chargements CH1,CH2...,CHN en utilisant un procédé d'identification selon l'une quelconque des revendications 1 à 7; 20 une étape de sélection (H30) d'au moins un chargement parmi la pluralité chargements CH1,CH2,...,CHN, ledit au moins un chargement sélectionné présentant l'échauffement le plus élevé pour la séquence d'initiation identifiée parmi ladite pluralité de chargements CH1,...,CHN; et 25 - une étape de retrait (H40) de l'enceinte dudit au moins un chargement sélectionné. 14. Procédé de conception d'un système de production d'hydrogène gazeux selon la revendication 13 comprenant, suite à l'étape 30 de retrait, une étape d'identification (H50) d'une nouvelle séquence d'initiation de la pluralité de chargements obtenue après retrait dudit au moins un chargement sélectionné en utilisant un procédé d'identification selon l'une quelconque des revendications 1 à 7. 35 15. Dispositif d'identification (1) d'une séquence d'initiation d'une pluralité de chargements indépendants notés CH1,CH2...,CHN, N 33 désignant un entier supérieur à 1, chaque chargement indépendant comprenant un matériau apte à générer de l'hydrogène gazeux sous sollicitation thermique, les chargements indépendants CH1,CH2...,CHN étant maintenus espacés les uns des autres dans une enceinte par un bloc 5 de maintien en un matériau isolant thermiquement, ledit dispositif étant caractérisé en ce qu'il comprend : - un module de détermination (1A) configuré pour déterminer une pluralité de séquences d'initiation des chargements CH1,CH2...,CHN, chaque séquence d'initiation fournissant un ordre d'initiation différent 10 des chargements CH1,CH2...,CHN; - des modules, activés pour chaque séquence d'initiation déterminée, et comprenant : o un module d'estimation (1B), configuré pour estimer pour chaque chargement CHn, n désignant un entier compris entre 1 15 et N, un échauffement de ce chargement CHn dû à une initiation avant l'initiation du chargement CHn d'au moins un autre chargement CHk adjacent au chargement CHn, k désignant un entier différent de n; o un module d'évaluation (1C), configuré pour évaluer pour ladite 20 séquence d'initiation, la valeur maximale des échauffements estimés pour les chargements CH1,CH2,...,CHN ; et - un module de sélection (1D), configuré pour sélectionner une séquence d'initiation parmi ladite pluralité de séquences d'initiation pour laquelle la valeur maximale évaluée est égale à la plus petite des 25 valeurs maximales évaluées pour ladite pluralité de séquences d'initiation ou qui est inférieure ou égale à un seuil d'échauffement prédéterminé.