Brevet FR2963090A1 — publié le 20120127

PROCEDE D'ESTIMATION DE LA CHARGE THERMIQUE IMPOSEE A UN REFRIGERATEUR CRYOGENIQUE, PRODUIT PROGRAMME ASSOCIE ET PROCEDE DE REGULATION DU REFRIGERATEUR. La présente invention se rapporte au domaine de la réfrigération d'installations. La présente invention se rapporte en particulier à la réfrigération cryogénique d'installations. La présente invention se rapporte encore à la réfrigération io cryogénique d'installations susceptibles de fonctionner selon un régime instationnaire. Un régime de fonctionnement instationnaire peut être rencontré dans de nombreuses applications. C'est par exemple le cas des installations mettant en oeuvre 15 des champs magnétiques intenses. Un exemple d'installation mettant en oeuvre des champs magnétiques intenses est le tokamak. Le tokamak est une chambre de confinement magnétique destinée à contrôler un plasma pour étudier la possibilité de produire de l'énergie par fusion nucléaire. 20 Dans les tokamaks, on peut rencontrer un régime de fonctionnement instationnaire pulsé. Dans ce cas, si l'on connaît globalement le profil des charges pulsées imposées au réfrigérateur, on ne sait pas exactement quand une impulsion de charge advient. De plus, d'autres perturbations imprévisibles peuvent advenir, lesquelles sont liées au 25 fonctionnement du tokamak. Pour générer des champs magnétiques intenses, sans provoquer la destruction des électroaimants, on utilise des électroaimants supraconducteurs. Le fonctionnement supraconducteur d'un électroaimant nécessite un maintien de sa température au-dessous de sa température 30 critique. Ceci est obtenu grâce à un réfrigérateur cryogénique, par exemple refroidi par hélium. La figure 1 est un schéma représentatif du fonctionnement d'un tel réfrigérateur cryogénique. Le réfrigérateur cryogénique 1 comprend un compresseur 10 permettant de comprimer un gaz à température ambiante (To = 300K environ), en l'occurrence de l'hélium, depuis la pression atmosphérique Pa à une pression P de 20 bar environ. Plusieurs échangeurs thermiques 20, 21, 22, 23, 24 sont io disposés en parallèle jusqu'à un séparateur de phase 40 comportant un bain d'hélium liquide à la température Tf = 4,5K. Dans le cas d'espèce, cinq échangeurs thermiques sont prévus. Ces échangeurs thermiques, à contre-courant, permettent de diminuer petit à petit la température de l'hélium circulant dans les différents 15 circuits de ces échangeurs thermiques. Par ailleurs, pour chaque échangeur thermique, la pression des deux circuits d'échange de l'échangeur diffère, de sorte que le circuit « amont » est un circuit à haute pression et chaud et le circuit « aval » est un circuit à basse pression, plus froid. Entre le circuit « amont» lié au premier échangeur thermique 20 et le circuit « aval» du 20 dernier échangeur thermique 24, la pression passe ainsi de P = 20 bar en sortie du compresseur à une pression légèrement supérieure à la pression atmosphérique, pour éviter des problèmes de cavitation en basse pression. Le bain 41 d'hélium liquide est donc à la pression atmosphérique. 25 Le réfrigérateur cryogénique 1 comprend également plusieurs moyens pour extraire un travail. Dans le cas présent, ces moyens sont formés par deux turbines 30, 31. La première turbine 30 prélève du travail à la sortie du circuit 30 basse pression du premier échangeur thermique 20 et réinjecte ce travail à l'entrée basse pression du deuxième échangeur thermique 21. La deuxième turbine 31 prélève du travail à la sortie du circuit basse pression du troisième échangeur thermique 22 et réinjecte ce travail à l'entrée basse pression du quatrième échangeur thermique 23. Ces turbines 30, 31 sont complémentaires aux échangeurs 5 thermiques et participent, par le prélèvement d'un travail, au refroidissement de l'hélium. Le réfrigérateur cryogénique 1 comprend enfin une vanne 50 Joule-Thomson, disposée entre la sortie du circuit basse pression du dernier échangeur thermique 24 et le bain 41 d'hélium liquide à 4,5K sous pression io atmosphérique. Cette vanne 50 permet de liquéfier l'hélium gazeux obtenu en sortie du circuit basse pression du dernier échangeur thermique 24 en effectuant une détente s'accompagnant d'une chute de température de l'hélium. Le bain 41 d'hélium liquide permet alors de fournir une 15 puissance frigorifique pour maintenir les électroaimants de l'installation à champs magnétiques intenses, par exemple un tokamak, en fonctionnement supraconducteur, compte tenu de la puissance thermique imposée par l'installation à ce bain 41 d'hélium liquide. La puissance thermique imposée par l'installation au bain 41 20 d'hélium liquide est également appelée charge thermique. Ce type de réfrigérateur est conçu et dimensionné pour un fonctionnement en régime stable ou lentement variable du champ magnétique intense généré par l'installation, c'est-à-dire pour un régime de fonctionnement thermique permanent ou quasi-permanent du réfrigérateur. 25 En effet, ces régimes de fonctionnement de l'installation permettent d'obtenir une charge thermique stable au niveau du réfrigérateur cryogénique. Le fonctionnement en régime thermique permanent ou quasi-permanent du réfrigérateur cryogénique permet d'assurer un refroidissement adéquat des électroaimants, de façon stable et sûre. 30 Cependant, les futures installations visant à étudier la possibilité de production d'énergie par fusion nucléaire activée par des champs magnétiques intenses prévoient de mettre en oeuvre des champs magnétiques pulsés. C'est le cas des projets de tokamaks ITER (France) et JT60_SA (Japon). Les variations du champ magnétique dans le tokamak impliquent alors une variation similaire de la charge thermique imposée au réfrigérateur cryogénique. Un exemple de l'évolution de la charge thermique imposée au réfrigérateur cryogénique sur le futur tokamak japonais JT60_SA, destiné à to fonctionner en régime pulsé, est représenté sur la figure 2. On note que les impulsions sont a priori périodiques, dans le cas d'espèce selon une période de 3000s, la durée d'une impulsion étant environ de 100s. II n'en reste pas moins que des processus aléatoires peuvent altérer cette évolution, comme cela été présenté précédemment. 15 Sur cette figure 2, on peut constater que le régime de fonctionnement pulsé du tokamak imprime un régime de fonctionnement également pulsé (courbe 100) de la charge thermique imposée au réfrigérateur cryogénique. Cette figure 2 représente également (courbe 101) la charge thermique moyenne imposée par le tokamak au réfrigérateur. 20 Or, les réfrigérateurs cryogéniques utilisés sur les installations existantes ne sont pas adaptés pour assurer, dans ce régime pulsé, un refroidissement adapté. En effet, une augmentation de la charge thermique entraîne une augmentation du débit d'hélium refroidi, en retour vers les électroaimants 25 de l'installation. Ceci refroidit l'ensemble du réfrigérateur, car le froid apporté par la charge thermique déséquilibre les échangeurs 20, 21, 22, 23 et 24 entre la partie haute pression et la partie basse pression. De plus, l'évaporation provoquée par la charge sur le bain d'hélium augmente 30 instantanément le débit de retour côté basse pression de chaque échangeur thermique, ce qui déséquilibre tout le réfrigérateur cryogénique. Une augmentation conséquente de la charge thermique peut même conduire à l'arrêt complet du réfrigérateur cryogénique. Pour remédier à ce problème, il a déjà été proposé de lisser l'impact du champ magnétique pulsé sur l'évolution de la charge thermique 5 imposée au réfrigérateur cryogénique. Ce lissage consiste à limiter les variations de la charge thermique réellement imposée au réfrigérateur cryogénique, pour assurer un fonctionnement nominal du réfrigérateur cryogénique et par suite, éviter notamment tout arrêt du réfrigérateur. to A cet effet, il a été proposé d'agir de façon mécanique et/ou de façon thermique, soit en installant des moyens dédiés à l'intérieur du réfrigérateur cryogénique, soit en installant un dispositif additionnel entre l'installation mettant en oeuvre des champs magnétiques intenses et son réfrigérateur cryogénique. 15 Par exemple, le document Dauguet & al., Advances in Cryogenic Engineering : Transactions of the Cryogenic Engineering Conference, CEC Vol. 53, edited by J.G. Weisend II, pp. 564-569 propose de maintenir une charge thermique constante sur le réfrigérateur cryogénique. Pour cela, le réfrigérateur cryogénique comporte de nombreux 20 échangeurs thermiques et vannes cryogéniques activées afin que la puissance frigorifique délivrée par le réfrigérateur reste stable, à une valeur correspondant à une charge thermique moyenne imposée par le tokamak. Ainsi, le réfrigérateur cryogénique fonctionne selon un régime permanent ou quasi-permanent, même en présence d'un régime en fonctionnement pulsé du 25 tokamak. Une solution voisine est proposée dans le document WO 2009/024705. Selon un autre exemple, le document "Design of the ITER- FEAT cryoplant to achieve stable operation over a wide range of experimental 30 parameters and operation scenarios", Claudet & al., Fusion Engineering and Design 58-59 (2001), pp. 205-209 propose d'installer un dispositif intermédiaire entre le tokamak et le réfrigérateur cryogénique. Ce dispositif permet d'effectuer une dérivation, en amont du réfrigérateur, d'une partie du débit d'hélium supplémentaire obtenu lors d'un pic de la charge thermique imposée au réfrigérateur cryogénique. Ainsi, le réfrigérateur cryogénique ne voit pas l'augmentation de débit d'hélium liée au régime de fonctionnement pulsé du tokamak. Les solutions proposées sont basées sur des dispositifs mécaniques et/ou thermiques visant à lisser la variation de charge thermique lo susceptible d'être imposée au réfrigérateur cryogénique. Ces solutions fonctionnent correctement. Toutefois, elles nécessitent un complément d'installation (échangeurs thermiques, vannes,...) qui peut rapidement s'avérer coûteux et, parfois, délicat à mettre en oeuvre. 15 Un objectif de l'invention est de résoudre l'un au moins des inconvénients des réfrigérateurs cryogéniques existants. Un autre objectif est de résoudre l'un au moins des inconvénients des réfrigérateurs cryogéniques existants, lorsque l'installation à réfrigérer fonctionne selon un régime instationnaire. 20 Comme cela a déjà été déjà mentionné précédemment, un régime instationnaire peut être rencontré dans de nombreux domaines. En conséquence, l'invention ne saurait être limitée aux installations générant un champ magnétique intense, pulsé, tel qu'un tokamak, mais peut s'étendre à tous types d'installations prévoyant une réfrigération cryogénique. 25 Pour atteindre l'un au moins de ces objectifs, l'invention propose un procédé d'estimation d'une charge thermique imposée à un réfrigérateur cryogénique comprenant un séparateur de phase comportant un bain de liquide frigorifique, dans lequel on estime cette charge thermique avec un programme mis en oeuvre par ordinateur, ledit programme étant basé sur 30 un bilan de masse effectué sur le séparateur de phase permettant d'exprimer l'évolution de la dérivée temporelle de la hauteur du bain de liquide frigorifique dans le séparateur de phase. Le procédé selon l'invention pourra prévoir d'autres caractéristiques techniques, prises seules ou en combinaison: la charge thermique imposée au réfrigérateur cryogénique est une charge thermique instationnaire ; la charge thermique instationnaire imposée au réfrigérateur cryogénique est pulsée ; il est prévu une étape de détermination des variables pertinentes du io réfrigérateur cryogénique représentant l'évolution de la dérivée temporelle de la hauteur du bain de liquide frigorifique dans le séparateur de phase, préalablement à l'étape d'estimation de la charge thermique avec le programme mis en oeuvre par ordinateur ; le réfrigérateur cryogénique comprenant une vanne à l'entrée du 15 séparateur de phase, lesdites variables pertinentes comprennent au moins le degré d'ouverture de la vanne, la température en amont de la vanne et leurs dérivées temporelles respectives ainsi que la charge thermique ; le programme mis en oeuvre par ordinateur comprend une étape dans 20 laquelle on filtre les données relatives à l'évolution de la dérivée temporelle de la hauteur du bain de liquide frigorifique dans le séparateur de phase. Pour atteindre l'un au moins de ces objectifs, l'invention propose également un procédé de régulation du réfrigérateur cryogénique 25 soumis à une charge thermique, dans lequel on estime la charge thermique imposée au réfrigérateur avec le procédé d'estimation d'une charge thermique imposée à ce réfrigérateur selon l'invention, puis, on régule au moins un paramètre de fonctionnement du réfrigérateur, en fonction de la valeur de la charge thermique préalablement estimée. 30 2963090 s Le procédé de régulation selon l'invention pourra prévoir d'autres caractéristiques techniques, en particulier: une étape dans laquelle la régulation est réalisée en modifiant le degré d'ouverture de la vanne du réfrigérateur, située à l'entrée du séparateur 5 de phase. Pour atteindre l'un au moins de ces objectifs, l'invention propose encore un produit programme d'ordinateur comprenant des instructions de codes de programme pour la mise en oeuvre du procédé d'estimation d'une charge thermique imposée à un réfrigérateur cryogénique io selon l'invention. D'autres caractéristiques, buts et avantages de l'invention seront énoncés dans la description détaillée ci-après faite en référence aux figures suivantes : - la figure 3 est un schéma représentant un bilan de masse 15 réalisé sur le séparateur de phase du réfrigérateur cryogénique, par exemple illustré sur les figures 1 ou 4 ; - la figure 4 est un schéma d'un réfrigérateur cryogénique portant le nom de « station 800W@4.5K » installé au Commissariat à l'Energie Atomique, sur lequel le procédé selon l'invention a été testé et 20 validé ; - la figure 5 est un graphique représentant d'une part, l'évolution de la charge thermique réellement imposée au réfrigérateur cryogénique représenté sur la figure 4 lors d'un essai et, d'autre part, l'évolution de la charge thermique estimée par le procédé conforme à 25 l'invention. Le procédé comprend notamment une étape d'estimation de la charge thermique w imposée au réfrigérateur 1, 1' cryogénique avec un programme mis en oeuvre par ordinateur. L'élaboration de ce programme est détaillée ci-après. Ce programme s'appuie notamment sur un bilan de conservation de la masse effectué sur le séparateur de phase 40, 40' du réfrigérateur 1, 1' cryogénique. Ce bilan de masse effectué sur le séparateur de phase 40, 40' est schématisé sur la figure 3. Typiquement, il s'agit d'un séparateur de phase 40, 40' tel que représenté sur les figures 1 ou 4 respectivement, lequel comporte un bain 41, 41' de liquide frigorifique, par exemple de l'hélium. L'équation de conservation de la masse sur le séparateur de phase 40, 40' s'écrit : i0 dt = m+ - f (u, , P, T) - L (Eq. 1) L, où : 15 h est la hauteur du bain de liquide frigorifique 41, 41' dans le séparateur de phase 40, 40' (% de la hauteur maximale hmax atteignable dans le bain); m+, le débit de gaz entrant dans le séparateur de phase (g/s); m_ , le débit de gaz sortant du séparateur de phase (g/s); w, la charge thermique imposée au séparateur de phase, c'est-à-dire au 20 réfrigérateur, par l'installation (W); ul, le degré d'ouverture de la vanne 50, 50' (%) ; L,,, la chaleur latente de vaporisation (J/g) ; et f(ui, P, T) une fonction dépendant du degré d'ouverture ul de la vanne 50, 50', de la pression P en amont de la vanne et de la température T en amont de 25 cette vanne. La fonction f(ui, P, T) peut prendre différentes formes développées, en fonction de la précision souhaitée sur l'estimation de la i0 dh dérivée temporelle ar de la hauteur h du liquide de refroidissement dans le séparateur de phase 40, 40', laquelle est par la suite notée h. Dans l'exemple développé ci-après, la quantité h dépend des paramètres suivants : ul, T.u1, u, , T et w où u, et T sont les dérivées s temporelles du degré d'ouverture de la vanne 50, 50' et de la température en amont de cette vanne, respectivement. Ces paramètres apparaissent comme représentatifs des paramètres influant principalement sur le bilan de masse réalisé sur le séparateur de phase 40, 40'. Dans cet exemple, l'influence de la pression P dans la io fonction f n'est donc pas prise en considération. Ce choix de paramètres pourrait être différent, en fonction de la précision souhaitée du modèle, notamment en y intégrant des paramètres liés à la pression P. Selon un autre exemple, le modèle prendra en compte, afin 15 de définir la quantité h, les paramètres suivants : u1, T, et leurs dérivées temporelles respectives, ainsi que la charge thermique w. On cherche alors à exprimer l'équation (Eq. 1) sous une forme linéaire, soit sous la forme de (Eq. 2): (Tu. 20 où : al, a2, a3, a4, a5 sont des coefficients à déterminer. L'équation (Eq. 2) peut s'écrire sous la forme de l'équation (Eq. 3): h L q. Il compte tenu des notations suivantes : Pour estimer la charge thermique w imposée au réfrigérateur 1, 1' cryogénique par l'installation, nous nous intéressons au filtrage de la dérivée temporelle h de la hauteur h de liquide de refroidissement dans le séparateur de phase, laquelle est exprimée dans le champ de Laplace. Le filtrage permet d'éliminer tous les « bruits » qui sont io susceptibles d'influer sur la détermination de la variable h. Ainsi, si l'on note H(s), la dérivée temporelle h de cette hauteur h dans le champ de Laplace et, Y(s) la fonction obtenue après filtrage du premier ordre de la quantité H(s), on peut exprimer, a priori, ce lien sous la forme de l'équation (Eq. 4) : s 15 1. + 17°E où rf est la constante de temps de ce filtre du premier ordre, qu'il faut définir. Pour définir la constante de temps, on procède de la façon suivante. On impose un profil d'ouverture à la vanne 50, 50', par exemple un échelon. Puis, on suit l'évolution de la dérivée temporelle de la hauteur du 20 bain de liquide. On est alors en mesure de déterminer la constante de temps, de façon classique pour un homme du métier. L'équation (Eq. 4) correspond alors, dans le champ réel, à l'équation différentielle (Eq. 5) suivante : 25 où : y est la variable réelle associée à la fonction filtrée Y(s) dans le champ de Laplace de l'équation (Eq. 4) et y, la valeur réelle intégrée de la variable y . En insérant la relation fournie par l'équation (Eq. 3) issue du 5 bilan de conservation de la masse effectué sur le séparateur de phase 40, 40' dans l'équation (Eq. 5) précédente, on obtient alors l'équation (Eq. 6) : + L-,r - -- +a On comprend alors que cette équation différentielle traduit le bilan de conservation de la masse effectué sur le séparateur de phase, suite à 10 un filtrage du premier ordre effectué sur la variable h. On introduit alors le vecteur d'état, xW sous la forme de l'égalité (Eq. 7) : Cette notation permet alors d'écrire l'équation (Eq. 6) sous la 15 forme de l'équation (Eq. 8) suivante : Puis, en remplaçant la variable y par la variable 1 [x,,, +LZz], rf on obtient l'équation : 20 ce qui donne finalement la représentation d'état suivante, (Eq. 9) : dont xW est l'état ; q, z et w les entrées et y la sortie. Dans un souci de clarté de présentation, la représentation d'état (Eq. 9) peut encore s'exprimer sous la forme (Eq. 10) suivante : r t; compte tenu des notations de matrices suivantes : G _z.:b La représentation d'état (Eq. 9/Eq. 10) permet finalement d'exprimer la variable y, filtrée, représentative de la hauteur h du liquide de réfrigération dans le séparateur de phase, en fonction de q, z, w ainsi que des coefficients al, a2, a3, a4, a5 et z f . Cette représentation d'état est facilement implémentable dans un automate programmable. Cette représentation d'état peut être étendue pour estimer la charge thermique imposée au réfrigérateur cryogénique. 5 Pour cela, on définit l'état étendu de la représentation d'état (Eq. 9), par la matrice suivante, (Eq. 11) : où le deuxième vecteur d'état s, représente la charge thermique w. On exprime alors la représentation d'état (Eq. 10) sous la forme suivante, nommée (Eq. 12) : 10 Cette représentation suppose que s,,, = 0, c'est-à-dire que la variable s,,, représentative de la charge thermique w imposée au réfrigérateur, présente un profil constant. Cette hypothèse est correcte lorsque la perturbation arrivant sur le bain de liquide se présente sous la forme d'échelons. Cela est notamment le cas pour un tokamak. 15 Bien entendu, d'autres hypothèses pourraient être avancées. Par exemple, on peut supposer que la perturbation arrive sous la forme d'une rampe ou d'une sinusoïde. Cela pourrait être le cas avec des réfrigérateurs cryogéniques, employés pour d'autres installations que le tokamak. 20 15 Si l'on note, (Eq. 13) : .J (T, , alors la représentation d'état de (Eq. 12) s'écrit sous la forme (Eq. 14) : Puis, si l'on note Ad et Bd les matrices discrètes représentatives des matrices A et B respectivement, pour une période d'échantillonnage Te, les relations précédentes s'écrivent : U{T. j . flTr où+ est la forme discrète de la matrice On rappelle que lorsque l'on connaît les matrices continues 15 d'une représentation d'état A, B, C, D, les matrices discrètes associées Ad, Bd, Cd et Dd s'écrivent : T Ad = eAT avec Te la période d'échantillonnage ; Bd = f eAs.B.ds ; Cd = C et Dd 0 =D. 20 Les relations exprimant e et y peuvent s'écrire sous une forme réduite. + [B' U T [ {roc i0 Pour cela, on rappelle l'équation de l'observateur: + A = Ad + Bd U + L où L est une matrice colonne à deux lignes à valeurs réelles, nommé gain de l'observateur. Le gain de l'observateur dépend uniquement du temps de 5 réponse tr choisi pour l'observateur. Puis, on utilise la relation exprimant y pour la variable y . Dans ce cas, y s'écrit : y = [C, + [D] . U(T, u1). En remplaçant cette dernière expression dans l'équation de l'observateur, on aboutit finalement à la forme réduite mentionnée ci-dessus, 10 à savoir l'équation (Eq. 15) : dans laquelle le gain de l'observateur L est choisi pour que les valeurs propres de la matrice Ad - LC soient à parties réelles strictement négatives. 15 Comme l'homme du métier peut le comprendre aisément, les valeurs propres de la matrice Ad - LC correspondent aux pôles de la fonction de transfert correspondante. Ces pôles sont notamment liés à la constante de temps zf , dans la mesure 20 où ce paramètre intervient dans l'expression de cette matrice. Par suite, ces pôles sont liés au temps de réponse tr choisi pour l'observateur. Il faut cependant noter que la variable y n'est pas disponible dans l'équation (Eq. 15). Il faut donc la déterminer par ailleurs. Pour cela, on 25 revient à l'équation différentielle (Eq. 5), à savoir : io Puis, on définit le vecteur : et on réécrit l'équation (Eq. 5) sous la forme suivante : E +h_ que l'on peut exprimer sous la forme de la représentation d'état suivante : Cette représentation d'état peut alors être exprimée de façon discrète, pour une période d'échantillonnage Te sous la forme suivante (Eq. 16): où 77+ est la forme discrète de la variable 77, les matrices [Rd] et [Ed] étant les 15 matrices discrètes de la représentation d'état de l'équation (Eq. 5]. La détermination de ces matrices est réalisée en fonction des matrices continues [R] = [-1tz f ] et E = [-1/1-f ], par les relations : TÎe.E.ds. Rd = e n.Te avec Te la période d'échantillonnage ; Ed = 0 17 + h. En regroupant alors les équations (Eq. 15) et (Eq. 16), on peut écrire un ensemble (Eq. 17) d'équations: + [E -L']. + Tul + L + 1. 15 Cet ensemble d'équations peut être exprimé sous la forme suivante, (Eq. 18) : ( Rd oi,x ., , 2 «3, o, considérant la notation : ou, plus simplement, (Eq. 19) : + = T, y1, h) avec les matrices suivantes : ( Rd 0-1 Ad- Le) u., .x2 01x4 Ed -L D h Cobs = (0,0,1), et L'équation (Eq.19) permet ainsi d'estimer la charge thermique w imposée au réfrigérateur, sous une forme aisément implémentable dans un automate programmable. La charge thermique w ne peut cependant être estimée à partir de (Eq. 19), qu'à partir du moment où les matrices Aobs et Bobs sont parfaitement déterminées, ce qui revient à identifier les coefficients al à a5, déterminer les valeurs de la matrice L et choisir la constante de temps r du filtre de premier ordre choisi. L'identification des coefficients al à a5 s'effectue de la façon suivante. On intègre l'équation (Eq. 2) entre un temps de référence t1 et le temps t, puis on filtre. L'intégration aboutit à l'équation (Eq. 20) : h(t)- h(tl)= M(t) * a où m(t) est la matrice ligne donnée par : 19 T r) ( T(t) - T(t1 et, a le vecteur colonne des coefficients a1 à a5. 5 Le filtrage de l'équation (Eq. 20) est ensuite réalisé avec le filtre suivant : Les coefficients al à a5 sont alors déterminés en minimisant, par la méthode des moindres carrés, les relations formant un système d'équations (Eq. 21) dont les inconnus sont les coefficients de la matrice a: f'..( _11)(t 1I)(t2) .F(h)( * - laquelle permet de résoudre l'équation (Eq. 20) sous forme filtrée, entre le io temps tl initial et le temps final tN, avec N = 5 dans le cas d'espèce. Pour ce faire, il est nécessaire de mesurer les valeurs des coefficients des autres matrices. Ces mesures ont été effectuées sur le réfrigérateur 1' représenté schématiquement sur la figure 4. 15 Ce réfrigérateur 1' est similaire au réfrigérateur 1 décrit à l'appui de la figure 1, notamment au niveau du séparateur de phase 40'. Ce réfrigérateur 1' diffère cependant légèrement du réfrigérateur décrit à l'appui de la figure 1 par différents aspects. En effet, il ne comporte que quatre échangeurs thermiques 20 20', 21', 22', 23' pour baisser la température du fluide de refroidissement, en l'occurrence de l'hélium, de 300K à 4,5K. La pression fournie par le compresseur est de 16 bar. Par ailleurs, une seule turbine 31' est prévue et un échangeur thermique 60' additionnel (azote, 80K) est prévu au niveau du premier échangeur thermique 20'. 25 Les coefficients de la matrice a ont été calculés pour une pression P en sortie du compresseur de 16 bar. Par ailleurs, la température en amont de cette vanne 50' est de T = 7K. Le temps t1 considéré est tl = 1s, à partir duquel les données sont enregistrées. La période considérée pour la prise des mesures est de 3s, si bien que dt = t;.+1 - t; = 3s, pour i variant de 1 à N. Dans ces conditions, la résolution du système d'équations (Eq. 21) permet d'aboutir aux coefficients du tableau 1 suivant : Coefficients de la matrice a valeur al 0,003710768747845 a2 -0,000396506514909 a3 0, 024617305832635 a4 0,286747643082885 a5 -0,000240828045189 Tableau 1 Une fois les coefficients al à a5 obtenus, ceux-ci sont entrés dans le programme, tout comme la constante de temps r f et les coefficients du gain de l'observateur L. Un temps de réponse à 5%, noté tr5%, a été choisi pour l'observation. Dans le cas d'espèce, il vaut tr5% = 60s. Cela signifie que 15 l'observateur est calé afin d'obtenir une estimation à 5% de la charge thermique w, en 60s. La constante de temps z f est déterminée par la relation tr5% =3'r f . Ainsi, la constante de temps r1 vaut z f = 20s. Par ailleurs, pour déterminer les valeurs propres de la matrice 20 Ad - LC , on détermine les pôles z de la fonction de transfert 3.dt correspondante, soit z =e 0'5% = (0.8607 X0.8607, dans la mesure où dt = 3s comme mentionné précédemment. Ces données permettent de définir notamment les matrices Aobs et Bobs. L'estimation de la charge thermique w avec ce programme peut alors être réalisée avec le système d'équations (Eq. 19). La figure 5 représente l'estimation de la charge thermique obtenue par le programme décrit ci-dessus de l'évolution réelle de la charge thermique imposée au réfrigérateur représenté sur la figure 4. Comme on peut le constater sur cette figure 5, bien que la charge thermique imposée au réfrigérateur 1' soit fortement instationnaire, et io même aléatoire, l'estimation est de très bonne qualité. L'estimation suit l'évolution réelle de la charge thermique lors des pics de puissance. Cette estimation est également très bonne en régime permanent, par exemple entre le temps t = Os et le temps t = 3.104 s sur la figure 5. 15 Finalement, il est possible d'estimer correctement, en temps réel, l'évolution de la charge thermique imposée au réfrigérateur cryogénique. Cette estimation est très bonne que le régime de fonctionnement soit permanent ou instationnaire. Ceci est particulièrement intéressant, car il est alors possible 20 d'effectuer une régulation d'un ou de plusieurs paramètres(s) de fonctionnement du réfrigérateur cryogénique. On peut par exemple imaginer de réguler la température en sortie de la turbine 31', la hauteur du bain de liquide dans le séparateur ou d'autres paramètres. 25 Considérant l'exemple ci-dessus, on peut alors physiquement agir sur les variables tels que l'ouverture de la vanne, la température et/ou la pression en amont de cette vanne etc..., afin de réguler le ou chaque paramètre de fonctionnement considéré. En particulier, dans le cas des régimes pulsés d'un tokamak, il 30 est possible d'estimer en temps réel la valeur de la charge thermique imposée au réfrigérateur cryogénique, sans savoir ce qui se passe à l'extérieur du réfrigérateur. Cette régulation permet alors d'assurer un refroidissement adéquat de l'installation, sans risque de provoquer un arrêt du réfrigérateur et 5 par suite, un arrêt de l'installation elle-même. Cette régulation est peu coûteuse, car elle ne nécessite aucun matériel important. Il faut noter que l'identification des coefficients de la matrice a est valable pour le réfrigérateur 1' représenté sur la figure 4, compte tenu des to paramètres choisis dans l'équation (Eq. 2). Une autre identification des valeurs de la matrice a serait nécessaire si d'autres paramètres étaient choisis pour élaborer l'équation (Eq. 2). Par exemple, si l'on souhaite avoir un modèle encore plus précis, il est envisageable de prendre en compte, en outre, la pression P en amont de la 15 vanne 50' ainsi que la dérivée temporelle de cette pression P. Une autre identification des valeurs de la matrice a serait également nécessaire pour ce même réfrigérateur fonctionnant dans des conditions différentes, par exemple si la pression en sortie du compresseur diffère. 20 REVENDICATIONS1. Procédé d'estimation d'une charge thermique imposée à un réfrigérateur (1,1') cryogénique comprenant un séparateur de phase (40, 40') comportant un bain (41, 41') de liquide frigorifique, dans lequel on estime cette charge thermique avec un programme mis en oeuvre par ordinateur, ledit programme étant basé sur un bilan de masse effectué sur le séparateur de phase (40, 40') permettant d'exprimer l'évolution de la dérivée temporelle de la hauteur du bain de liquide frigorifique dans le séparateur de phase. 2. Procédé selon la revendication 1, dans lequel la charge thermique imposée au réfrigérateur (1, 1') cryogénique est une charge thermique 15 instationnaire. 3. Procédé selon l'une des revendications précédentes, dans lequel la charge thermique instationnaire imposée au réfrigérateur (1, 1') cryogénique est pulsée. 20 4. Procédé selon l'une des revendications précédentes, dans lequel il est prévu une étape de détermination des variables pertinentes du réfrigérateur cryogénique (1, 1') représentant l'évolution de la dérivée temporelle de la hauteur du bain (41, 41') de liquide frigorifique dans le 25 séparateur de phase (40, 40'), préalablement à l'étape d'estimation de la charge thermique avec le programme mis en oeuvre par ordinateur. 5. Procédé selon la revendication précédente, dans lequel, le réfrigérateur cryogénique (1, 1') comprenant une vanne (50, 50') à l'entrée du 30 séparateur de phase (40, 40'), lesdites variables pertinentes comprennent au moins le degré d'ouverture de la vanne, la température en amont de la vanne et leurs dérivées temporelles respectives ainsi que la charge thermique. 6. Procédé selon l'une des revendications 4 ou 5, dans lequel le programme mis en oeuvre par ordinateur comprend une étape dans laquelle on filtre les données relatives à l'évolution de la dérivée temporelle de la hauteur du bain (41, 41') de liquide frigorifique dans le séparateur de phase (40, 40'). io 7. Procédé de régulation du réfrigérateur (1, 1') cryogénique soumis à une charge thermique, dans lequel on estime la charge thermique imposée au réfrigérateur avec le procédé selon l'une des revendications précédentes, puis, on régule au moins un paramètre de fonctionnement du réfrigérateur, en fonction de la valeur de la charge thermique 15 préalablement estimée. 8. Procédé selon la revendication précédente, dans lequel l'étape de régulation est réalisée en modifiant le degré d'ouverture de la vanne (50, 50') du réfrigérateur (1, 1'), située à l'entrée du séparateur de 20 phase (40, 40'). 9. Produit programme d'ordinateur comprenant des instructions de codes de programme pour la mise en oeuvre d'un procédé selon l'une des revendications 1 à 6. 25