Brevet FR2982014A1 — publié le 20130503

PROCEDE DE REGULATION DE LA TEMPERATURE D'UN ELEMENT REFROIDI PAR UN CRYOREFRIGERATEUR A FONCTIONNEMENT PERIODIQUE, DISPOSITIF DE MISE EN OEUVRE ASSOCIE ET INSTALLATION CRYOGENIQUE COMPRENANT CE DISPOSITIF. La présente invention se rapporte au domaine de la réfrigération d'un élément, en particulier à la réfrigération cryogénique de cet élément. La présente invention se rapporte plus particulièrement à la réfrigération cryogénique d'un élément au moyen d'un cryoréfrigérateur à fonctionnement périodique. Parmi les cryoréfrigérateurs à fonctionnement périodique connus, on peut citer les tubes à gaz pulsé ou les machines de GiffordMacMahon. Le refroidissement avec un cryoréfrigérateur à fonctionnement périodique est généralement employé pour refroidir à très basse température, typiquement entre 4K et 120K, des installations pour lesquelles les contraintes d'exploitation sont fortes ou pour lesquelles il n'est pas possible d'effectuer une distribution d'hélium liquide. De telles contraintes peuvent être rencontrées dans des applications nucléaires pour lesquelles l'élément à refroidir est disposé dans un milieu radioactif confiné. Dans ce cas, l'utilisation d'un réservoir d'hélium situé à l'extérieur du bâtiment nucléaire est en effet très contraignante, puisqu'il faut transférer l'hélium à l'intérieur du bâtiment, ce qui implique des pertes thermiques importantes. Par ailleurs, la récupération de l'hélium gazeux en sortie peut obliger à prendre des mesures de précautions particulièrement contraignantes. Ces contraintes peuvent également être rencontrées pour des applications spatiales, par exemple pour le refroidissement de certains composants d'un satellite. Un avantage du cryoréfrigérateur à fonctionnement périodique est qu'il fournit du froid sans avoir besoin d'apport d'hélium liquide extérieur par l'intermédiaire d'un réservoir de stockage d'hélium liquide. Cependant, un cryoréfrigérateur à fonctionnement périodique engendre des variations de température dans le temps qui se propagent dans l'installation et par suite, dans l'élément dont on cherche généralement à maintenir la température. Les variations de température subies par l'élément sont atténuées le long de la chaîne de refroidissement de l'installation, mais elles perturbent toujours le fonctionnement de cet élément. Io Un exemple d'application de refroidissement d'un élément par un cryoréfrigérateur à fonctionnement périodique est représenté sur la figure 1. La figure 1 est un schéma représentant une installation cryogénique 1' comportant un dispositif de refroidissement 1 pour un élément 15 10 à refroidir, par exemple la pince cryogénique 10 d'un système permettant d'obtenir une fusion nucléaire inertielle avec un laser de forte puissance, tel que le laser « MégaJoule ». La pince cryogénique 10 forme l'élément dont on cherche à réguler la température. 20 Le dispositif de refroidissement 1 comprend un tube 2 à gaz pulsé, un stabilisateur 3 pour atténuer les variations de température du tube 2 à gaz pulsé et un circuit de refroidissement 4 à l'hélium gazeux. Le circuit de refroidissement 4 comprend un premier échangeur thermique 41 avec le stabilisateur 3, un deuxième échangeur 25 thermique 42 avec la pince 10 et un moyen 43 pour faire circuler l'hélium gazeux dans le circuit de refroidissement 4. L'extrémité de la pince cryogénique 10 supporte une cible 11 pour le faisceau laser. La réussite des expériences de fusion avec le laser sur la 30 cible 11 dépend fortement de la stabilité thermique de cette cible 11. A ce jour, le dispositif de refroidissement 1 permet d'atteindre une variation de température de l'ordre de 10mK au niveau de la cible 11, cette cible 11 présentant une température généralement comprise entre 18K et 23K. La variation relative de la température de la cible 11 apparaît relativement faible, mais elle est encore gênante pour la réussite de certaines expériences de fusion avec ce laser. L'exigence de stabilité de la température au niveau de l'installation à refroidir avec un cryoréfrigérateur à fonctionnement périodique Io intéresse également d'autres applications, notamment dans le domaine spatial. En effet, selon le type de cryoréfrigérateur à fonctionnement périodique utilisé, il a été mesuré que l'amplitude de la variation de température au niveau du bout froid du cryoréfrigérateur pouvait être comprise 15 entre 300mK et 2200mK qui sont difficiles à atténuer malgré la présence d'un stabilisateur au sein du dispositif de refroidissement. Un objectif de l'invention est donc de proposer un procédé de régulation de la température d'une installation refroidie par un cryoréfrigérateur à fonctionnement périodique permettant de diminuer, par 20 rapport aux techniques existantes, les variations de température de cette installation. Un objectif de l'invention est également de proposer un dispositif pour mettre en oeuvre ce procédé. Un autre objectif de l'invention est de proposer une installation 25 cryogénique comportant ce dispositif. Dans le cas particulier d'un laser de forte puissance pour l'étude de la fusion nucléaire inertielle, l'invention vise notamment à obtenir des variations de température de l'ordre du mK au niveau de la cible. Pour atteindre l'un au moins de ces objectifs, l'invention 30 propose un procédé de régulation de la température d'un élément refroidi par un cryoréfrigérateur à fonctionnement périodique, caractérisé en ce qu'il comprend les étapes suivantes : (a) mesurer au moins l'évolution temporelle périodique de la température T de l'élément (10), cette évolution permettant d'obtenir un signal de température représentatif du fonctionnement périodique du cryoréfrigérateur (2) ; (b) construire un signal de flux thermique périodique dit équivalent, en fonction du signal de température obtenu à l'étape (a) ; (c) construire un signal de commande de flux thermique, dépendant du signal de flux thermique équivalent construit à l'étape (b), ledit signal de commande de flux thermique servant également, par retour, à la construction du signal de flux thermique équivalent lors de l'étape (b); (d) injecter un flux thermique dans l'élément (10), sur la base du signal de flux thermique construit à l'étape (c), pour compenser les variations périodiques mesurées de la température T de l'élément. Le procédé selon l'invention pourra comprendre les étapes suivantes, prises seules ou en combinaison : - le signal de commande de flux thermique dépend également d'un signal de commande pour améliorer la compensation des variations périodiques de la température de l'élément aux basses fréquences ; - le flux thermique injecté dans l'élément à l'étape (c) dépend de la différence entre le signal de commande et le signal de flux thermique équivalent construit à l'étape (b), pour limiter la puissance thermique injectée dans l'élément qui est nécessaire à la compensation des variations périodiques de sa température ; - l'étape (b) comprend également une étape additionnelle consistant à retarder le signal de flux thermique équivalent (S,Deq) préalablement construit, par rapport à la phase du signal de température déterminé à l'étape (a) ; il est prévu une étape consistant à estimer la valeur de ce retard temporel ; l'estimation de ce retard temporel comprend les étapes suivantes : (b1) injecter un signal d'identification du retard temporel dans un système thermique comprenant l'élément dont on cherche à réguler la température et un dispositif de mise en oeuvre du procédé, ce signal d'identification présentant une forme et une fréquence proches du signal représentatif du fonctionnement du cryoréfrigérateur; (b2) déterminer le déphasage (1) du signal d'identification en sortie du système thermique, ce déphasage étant associé au retard temporel qu'on cherche à estimer ; (b3) calculer le déphasage (1)0 de ce signal d'identification avec un modèle représentant le comportement thermique du système thermique en l'absence de tout retard; (b4) calculer le déphasage (p, associé au retard temporel par la relation (Pr = (I) Ikt; (b5) estimer le retard temporel, au moins à partir du déphasage (pr calculé à l'étape (b4) et de la fréquence du signal d'identification, ledit retard étant utilisé pour retarder le signal de flux thermique équivalent. Pour atteindre cet objectif, l'invention propose également un dispositif de mise en oeuvre du procédé de régulation de la température d'un élément refroidi par un cryoréfrigérateur à fonctionnement périodique selon l'une des revendications précédentes, caractérisé en ce qu'il comprend : - au moins un moyen pour mesurer l'évolution périodique de la température T de l'élément liée au fonctionnement périodique du cryoréfrigérateur; - au moins un moyen pour déterminer le flux thermique équivalent associé à cette variation périodique de la température T de l'élément, ce moyen ou cet ensemble de moyens comprenant un ordinateur ou un ensemble d'ordinateurs ou un processeur ou un ensemble de processeurs; - au moins un moyen pour injecter un flux thermique dans l'élément, à partir d'un signal de commande de flux thermique dépendant du signal de flux thermique équivalent ; - un moyen pour injecter le signal de commande de flux thermique au niveau de l'entrée dudit au moins un moyen pour déterminer le flux thermique équivalent. Le dispositif selon l'invention pourra comprendre l'une au moins des caractéristiques suivantes, prises seules ou en combinaison : - il est prévu un moyen capable de délivrer un signal de commande io pour améliorer la compensation des variations périodiques de la température de l'élément aux basses fréquences ; - le moyen pour injecter un flux thermique dans l'élément comprend un générateur de champ magnétique variable associé à un bain thermique isotherme, pour aimanter ou désaimanter 15 adiabatiquement un dispositif raccordé à l'élément ; - ledit au moins un moyen pour déterminer l'évolution périodique de la température de l'élément et ledit au moins un moyen pour injecter un flux thermique dans l'élément sont disposés l'un à côté de l'autre. Enfin, l'invention propose également une installation 20 cryogénique pour refroidir un élément avec un cryoréfrigérateur à fonctionnement périodique, ladite installation comprenant un stabilisateur relié au cryoréfrigérateur et un circuit de refroidissement comprenant un premier échangeur thermique relié au stabilisateur et un deuxième échangeur thermique relié à l'élément, caractérisé en ce qu'il comprend, en outre, un 25 dispositif selon l'invention. Les 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 2 est un schéma représentant un mode de 30 réalisation du procédé mis en oeuvre dans le cadre de l'invention; - la figure 3 est un schéma d'un mode de réalisation d'une version améliorée du procédé représenté sur la figure 2, cette version prenant notamment en compte le comportement dynamique de type retard pur, connu, des différents éléments constitutifs du dispositif de mise en oeuvre de l'invention; - la figure 4 est un schéma précisant comment est déterminé un signal de flux thermique, associé à l'évolution périodique de température de l'élément dont on cherche à réguler la température, dans le cadre du procédé mis en oeuvre sur la figure 3 - la figure 5 est un schéma d'un mode de réalisation d'une version améliorée du procédé représenté sur la figure 2, cette version prenant notamment en compte le comportement dynamique de type retard pur, qui n'est pas connu mais estimé, des différents éléments constitutifs du dispositif de mise en oeuvre de l'invention; - la figure 6 est un schéma représentant le fonctionnement d'un filtre non-linéaire susceptible d'être mis en oeuvre dans l'estimation effectuée dans le cadre du procédé représenté sur la figure 5 ; - la figure 7 représente une pince cryogénique selon une vue en perspective ; - la figure 8 est une représentation schématique de la pince de la figure 7 sur laquelle sont précisés l'implantation de différents capteurs de température et de différents actionneurs ; - la figure 9 est un schéma électrique équivalent au schéma pouvant être qualifié de thermique de la figure 8 ; - la figure 10 est un schéma représentant l'évolution de la température de l'élément dont on cherche à réguler la température, en l'occurrence la pince cryogénique du laser « MégaJoule », en fonction du temps, avec un filtre utilisé dans les procédés schématisés sur les figures 3 et 4; - la figure 11 est un schéma représentant l'atténuation de la perturbation de la température de l'élément dont on cherche à réguler la température, en l'occurrence la pince cryogénique du laser « MégaJoule », en fonction du retard généré par différents éléments constitutifs du dispositif de mise en oeuvre des procédés schématisés sur les figures 2 et 3 : - la figure 12 est un schéma d'un dispositif réel servant d'installation expérimentale test pour valider l'intérêt de l'invention ; - la figure 13 est un schéma électrique équivalent à l'installation expérimentale représentée sur la figure 12 ; - la figure 14 fournit l'évolution de la température de l'élément dont on cherche à réguler la température appartenant à l'installation expérimentale de la figure 12. Sur la figure 2, on a représenté le principe général du procédé selon l'invention, ainsi que des moyens du dispositif de mise en oeuvre de ce procédé. Il s'agit d'un procédé de régulation de la température T d'un élément 10 refroidi par un cryoréfrigérateur à fonctionnement périodique. Le procédé comprend les étapes suivantes : (a) mesurer au moins l'évolution temporelle périodique de la température T de l'élément, cette évolution permettant d'obtenir un signal de température représentatif du fonctionnement périodique du cryoréfrigérateur ; (b) construire un signal de flux thermique périodique dit équivalent, en fonction du signal de température obtenu à l'étape (a) ; (c) construire un signal de commande de flux thermique, dépendant du signal de flux thermique équivalent construit à l'étape (b), ledit signal de commande de flux thermique servant également, par retour, à la construction du signal de flux thermique équivalent lors de l'étape (b); (d) injecter un flux thermique dans l'élément, sur la base du signal de flux thermique construit à l'étape (c), pour compenser les variations périodiques mesurées de la température T de l'élément. Le procédé selon l'invention permet ainsi de compenser les 30 effets perturbateurs constatés au niveau de l'élément 10 (variation de température non souhaitée), en raison du fonctionnement périodique du cryoréfrigérateur. Ce procédé est destiné à s'appliquer à un système thermique 100' comprenant ledit élément 10 dont on cherche à réguler la température et l'actionneur 104. La boucle de régulation de la température de cet élément 10 est quant à elle associée à un dispositif 100, ce dernier étant destiné à fournir le signal de commande Sw de flux thermique à l'actionneur 104. Par exemple, l'élément 10 peut être formé par la pince d'un laser de forte puissance destiné à l'étude de la fusion nucléaire inertielle, tel Io que le laser « MégaJoule ». Cette pince cryogénique est en l'occurrence un assemblage de plusieurs pièces à refroidir. La mesure de la température T de l'élément à l'étape (a) peut être effectuée par un capteur de température 101 appartenant au système thermique 100'. 15 Le capteur de température 101 utilisé peut être un capteur résistif. La température est obtenue par calcul à partir de la mesure de la résistance et de la loi en température du capteur 101. Les performances demandées à ce capteur de température 101 sont principalement une grande stabilité, i.e. une faible dérive, et un 20 niveau de bruit aussi faible que possible. En d'autres termes, la reproductibilité de la mesure fournie par le capteur de température est importante. En revanche, la précision de la mesure fournie par ce capteur, c'est-à-dire la différence entre la température réelle et la température 25 mesurée par le capteur 101 est d'importance moindre dans la mesure où il est toujours possible d'effectuer un recalage sur un point fixe, avec des moyens de traitement adaptés connus de l'homme du métier. A cet effet, il est possible d'effectuer un étalonnage du capteur 101 avec un capteur de référence étalonné par un organisme agréé. En variante, il est possible 30 d'étalonner le capteur 101 en visualisant un changement de phase : passage 2 9 82014 liquide /solide (point triple) à l'aide d'un système fonctionnant par technique ombroscopique. Avec un capteur de température 101 de type résistif réalisé en germanium, on peut alors considérer les caractéristiques suivantes AR = 5 0,3m4 pour RN = 200 où RN est la résistance du capteur à la température nominale TN et AR, la variation de cette résistance autour de la température TN. Selon un autre exemple, avec un capteur de température 101 résistif de type CERNOXTM, on a environ AR = 5m0 pour RN = 10000. Le niveau des signaux mesurés ne dépasse pas 3mV à 4mV avec une bande passante de 10 quelques hertz. La température mesurée par le capteur 101 constitue alors un signal ST qui est transmis à un moyen 102, appartenant au dispositif 100 (boucle de régulation), permettant de construire le signal de flux thermique équivalent Scpeq. Ce signal de flux thermique Scpeci est l'équivalent en termes de flux de chaleur, de la température T mesurée par le capteur. Le moyen 102 est en fait un moyen de calcul, tel qu'un ordinateur, un ensemble d'ordinateurs, un processeur ou un ensemble de processeurs, permettant de mettre en oeuvre un modèle thermique représentatif du comportement du système thermique 100' comprenant l'élément 10. L'étape (b) ainsi mise en oeuvre par le moyen 102 permet donc de produire une information associée à la variation de température de l'élément 10 sous la forme d'un flux thermique équivalent. A cet effet, l'étape (b) intègre également le retour du signal Sw au niveau de l'entrée du moyen 102 pour avoir une information sur le flux thermique à injecter ; cf. figures 2 à 4 par exemple. Cette information peut alors être prise en compte par un actionneur 104 permettant d'injecter effectivement un flux thermique W dans l'élément 10 pour compenser les variations de température de l'élément 10. Ce flux thermique équivalent est reconstruit par le biais de l'équation (#11) du modèle qui est ultérieurement à l'appui des relations (#1) à (#11). Un modèle thermique susceptible d'être employé afin d'évaluer le flux thermique équivalent à une variation périodique de la température T de l'élément 10 est décrit ci-après, dans le cas où cet élément est constitué par la pince cryogénique 10 d'un laser à forte puissance pour l'étude de la fusion inertielle. Pour cela, on modélise le comportement thermique du io système thermique 100' comportant la pince 10, l'actionneur 104 et le capteur de température 101. Cette modélisation thermique s'effectue sous la forme d'une représentation d'état discrète. Dans le domaine de l'automatique, une représentation d'état permet de modéliser un système dynamique sous une 15 forme matricielle en utilisant des variables d'état. Cette représentation d'état doit alors rendre compte de l'état du système à n'importe quel instant futur si l'on possède les valeurs initiales. En l'occurrence, cette représentation est discrète pour son implémentation dans un moyen de calcul, tel qu'un processeur ou un ordinateur. 20 Cette représentation d'état discrète s'écrit sous la forme (#1): x(k + 1) = A. x(k) + B. U (k) T (k) = C. x(k) (#1) 25 OÙ : T représente la température de la pince cryogénique 10 qu'on cherche à réguler, U est un vecteur associé au flux thermique reçu par la pince cryogénique 10, 30 x est le vecteur d'état du système thermique 100', k est un incrément de temps discret ; et 2 9 82014 12 A, B et C sont des matrices à déterminer. Le vecteur d'état x du système thermique 100' est un ensemble minimal de variables d'état, c'est-à-dire de grandeurs temporelles, 5 nécessaires et suffisantes pour déterminer l'évolution du système thermique 100' quand on connaît les équations qui décrivent le fonctionnement de ce système thermique 100' et les entrées de ce système 100'. Par la suite, le vecteur x pourra être un vecteur colonne comportant n variables x1, x2,...x, où n, le nombre de composantes du vecteur x, correspond au degré de complexité du modèle thermique. Le nombre n définit ainsi le degré de complexité du modèle thermique modélisant le système thermique 100'. Dans ce cas, la matrice A est une matrice carrée de type nx*nx où nx définit la taille du vecteur d'état x du système thermique 100'. La matrice B est une matrice de type nx*nu où nu est la taille du vecteur U. Par exemple, si l'on considère uniquement un actionneur 104 dans le modèle thermique, alors nu = 1. C'est par exemple ce qui a été considéré dans l'essai réalisé sur la pince du laser MégaJoule, dans le chapitre « résultat d'un essai » présenté ultérieurement. Selon un autre exemple, on peut prendre en considération d'autres sources de flux de chaleur, par exemple en ajoutant celle qui provient de l'échangeur de la pince. C'est ce qui est effectué dans le chapitre « exemple d'application » présenté ci-après. La matrice C est une matrice ligne de type 1*nx. Il s'agit d'une matrice ligne dans la mesure où l'on cherche à réguler une seule température, en l'occurrence la température T de la pince cryogénique 10. Pour le refroidissement de la pince cryogénique 10 de ce type de laser, le cryoréfrigérateur 2 peut être un tube à gaz pulsé dont le fonctionnement périodique peut être représenté par une fonction de forme 30 sinusoïdale. En conséquence, la perturbation en température liée au fonctionnement du tube à gaz pulsé arrivant en entrée de la pince cryogénique 10 a également une forme sinusoïdale. Le signal de flux de chaleur Sp qui induit cette perturbation en température de la pince 10 est donc le signal que le dispositif 102 vise à reconstruire. Autrement dit, le flux de chaleur équivalent Screq provoque le même effet sur la température T de la pince 10 que le signal de flux de chaleur Sp créé par le dispositif de refroidissement de la pince. La représentation d'état discrète du signal sinusoïdal Sp s'écrit Io alors sous la forme (#2) : xsin(k + 1) = A sin. xsin(k) Sp (k) = Csin - X sin(k) (#2) 15 oÙ : A sin = [2. cos (21 fpTs) 11, avec fp la fréquence de fonctionnement du tube à -1 gaz pulsé et Ts la période d'échantillonnage du signal Sp, C sin = (1 0), et xsin(k) est le vecteur d'état associé au signal sinusoïdal Sp à l'instant k; 20 x n(k+1) est le vecteur d'état associé au signal sinusoïdal Sp à l'instant k+1. En combinant les représentations d'état (#1) et (#2), on obtient alors une représentation d'état discrète du comportement du système thermique 100' dans son ensemble, prenant en compte la nature sinusoïdale 25 du signal Sp ramené à un flux de chaleur de perturbation sur l'entrée du système thermique 100', à savoir : x(k +1) xs,n(k +1)_ 0i- A B - C x(k) + B - U (k) T(k) =[CI o A. xsin(k) 0 (#3) co B' Ao xo x(k) xsin(k)_ Les états du modèle thermique de l'équation (#3) sont estimés à l'aide d'un observateur d'états. Pour éliminer l'écart entre les états xp du modèle thermique de l'équation (#3) et leurs équivalents physiques xp , la différence entre la température mesurée par le capteur 101 et la température estimée par l'équation du modèle thermique (#3) est injectée par le vecteur de gain L dans cette équation (#3), conduisant à l'équation matricielle de l'observateur (#4) ci-dessous : x' (k +1) = Ao - xo (k) + Bo - U (k) : prediction T(k)= Co xo(k) erreur(k)=T(k)-Tme,',(k) x'(k +1) = x 0(k +1)+ L* erreur(k) xo(k+1)=x'(k+1)1(k)+L( ) * (Tm .esure .1( z Co xo (k)) : correction prediction x (k +1) = (210 - LC 0)- x '(k) + B ' - U (k) + L-T (k) Dans cette équation (#4), les coefficients des matrices Ao, Bo et Co s'expriment à partir des coefficients des matrices A, B, et C ainsi que des matrices A. et C. qui ont été fournies précédemment. En ce qui concerne les matrices A, B et C, des exemples numériques sont fournis plus loin dans la description. Par ailleurs, la température T est la température estimée par le modèle thermique et la température Tmesure est la température fournie par le capteur 101. Le vecteur de gain L est une matrice colonne, car on a une seule température à réguler. Le vecteur de gain L est déterminé par la méthode des « moindres carrés ». Nous détaillons ci-après la méthode des « moindres carrées » susceptible d'être mise en oeuvre pour déterminer le vecteur de gain L . La méthode des « moindres carrés » permet de minimiser l'erreur e(k), prise à l'instant k, entre les états x'(k) du modèle thermique et leurs équivalents physiques x '(k) à ce même instant k. Cette erreur e(k) est donc un vecteur colonne s'exprimant io sous la forme suivante : e(k) = x p(k)- x 0(k) (#5) L'évolution dans le temps de cette erreur e(k) peut être reliée au vecteur de gain L sous la forme du système d'équations suivant : 15 er (k +1) = AIT eT (k) - C'' v(k) (#6) v(k) -= (k) dans lequel le symbole T mis en exposant correspond à la transposée de la matrice considérée. 20 Compte tenu du système d'équations (#6), le vecteur de gain est alors la solution optimale qui minimise le critère : eT (n) - Q. e(n)+vT (n)- R- y(n)] (#7) 25 OÙ : Q et R sont des matrices de pondération. En utilisant le critère de l'équation (#7), il est possible de pondérer séparément l'erreur d'estimation pour les états xo du modèle thermique de l'équation (#3). En ajustant les valeurs de la matrice de pondération Q, une plus grande importance est donnée aux erreurs liées aux états de la perturbation, associée au signal Sp. La matrice de pondération R permet quant à elle de limiter la norme du vecteur de gain L. Le calcul du vecteur de gain L est effectué au moyen de la commande dlqr() implémentée dans le logiciel Matlab® de la manière io suivante = dlqr(A07 , ,Q,R) (#8) La commande dlqr permet de minimiser le critère de 15 l'équation (#7). Cette commande met en oeuvre la résolution de l'équation matricielle dite de « Riccati » en temps discret, dont la solution est notée S. L'équation matricielle de « Riccati » s'écrit alors dans ce cas : AoSA` -S-(A SC).(C'SCO1)1 -(C' A(1))+Q= 0 (#9) 20 Une fois la solution de l'équation (#9) obtenue, on détermine la transposée du vecteur de gain LT de la façon suivante : =(C'SC7,7 + R)1 SA:,) (#1O) 25 Une fois le vecteur de gain ainsi obtenu, on peut alors l'injecter dans l'équation de l'observateur (#4). Il convient de noter que d'autres méthodes que celles des « moindres carrés » pourraient être employées pour déterminer le vecteur de gain L. Par exemple, on pourrait 30 utiliser la méthode dite du « placement de pôle ». Enfin, pour résoudre l'équation de l'observateur (#4), il faut exprimer le flux thermique W injecté dans le système thermique 100' en fonction du flux thermique équivalent req. Ce flux thermique équivalent (Dec, est obtenu en récupérant l'état correspondant au sinus du vecteur d'état xo(k) de l'équation de l'observateur (#4). Cl)eq = [0 Csinl . x0 (k) (#11) Physiquement, le lien entre le flux thermique W (signal Sw) à to injecter dans la pince cryogénique 10 et le flux thermique équivalent Oeci (signal ScDeq) peut s'effectuer de la façon suivante, cf. figure 2. Pour les fréquences élevées, c'est-à-dire les fréquences typiquement supérieures à 0,01Hz, le signal de flux thermique équivalent Scrieq permet de former le signal Sw, afin de compenser les variations périodiques 15 de la température de la pince cryogénique 10. C'est en particulier le cas à la fréquence de 1,4Hz. En effet, la fréquence de fonctionnement du tube à gaz pulsé est généralement de 1,4Hz. On retrouve donc cette fréquence dans le signal de perturbation Sp et le signal de flux thermique équivalent ScDeq, qui vise à compenser les effets de la 20 perturbation. L'actionneur 104 recevant le signal Sw injecte alors, dans la pince 10, un flux thermique W qui est fonction de la période et de l'amplitude du signal de flux thermique équivalent S(Deq. Toutefois, ceci peut parfois poser des difficultés à basses 25 fréquences, c'est-à-dire les fréquences typiquement inférieures à 0,01Hz. Pour cette raison, le dispositif 100 peut comprendre un moyen 103 pour assurer la régulation du système thermique 100' en basses fréquences afin de garantir le suivi de la consigne. Ce moyen 103 délivre un signal de commande Sc spécifique, qui est amené à être modifié par le signal 30 de flux thermique équivalent Sctleci afin de former le signal de commande de flux thermique Sw arrivant effectivement sur l'actionneur 104. Le signal Sc est une représentation numérique du flux thermique à injecter dans la pince 10 pour assurer la régulation à ces basses fréquences. L'interaction entre les signaux Sc et ScIpeq ne pose pas de difficulté dans la mesure où le moyen 102 prend en compte les actions du moyen 103, par la présence du retour du signal Sw vers l'entrée du moyen 102, cf. figure 2. Le signal de commande Sc peut être une valeur constante. Plus généralement, le signal Sc peut être issu d'un régulateur, par exemple de type PID (pour Proportionnel Intégral Dérivé). L'interaction des signaux Sc et Scbeci commande ainsi la période et l'amplitude du signal Sw de flux thermique parvenant à l'actionneur 104. La présence du signal de commande Sc permet ainsi d'améliorer encore la régulation de température de la pince cryogénique. Cette variante est donc plus avantageuse, bien que la mise en oeuvre du signal de commande ne soit pas obligatoire. Avantageusement, le signal de flux thermique équivalent Seleq déterminé à l'étape (b) est injecté en opposition de phase. Dans ce cas et lorsque le signal de commande Sc est prévu, on peut exprimer la relation 20 entre Sc, Sw et ScIpeq sous la forme suivante : Sw = Sc - Cette expression est alors mise sous forme de matrices de flux discrétisées, dans l'équation matricielle de l'observateur (#4) pour en permettre la résolution. Le fait d'injecter le signal de flux équivalent en opposition de 25 phase permet de minimiser la puissance thermique requise pour l'actionneur 104. Il s'agit donc d'une solution préférentielle. Encore plus préférentiellement, le fait d'injecter un flux de chaleur en opposition de phase au plus proche de la température à réguler permet de minimiser encore plus la puissance thermique devant être injectée 30 dans la pince cryogénique 10 et donc, d'optimiser le dispositif de refroidissement 1. A cet effet, il convient de disposer l'actionneur 104 au plus près de la pince cryogénique 10. Par exemple, l'actionneur 104 pourra avantageusement être disposé à côté du capteur de température de la pince 10. Le dispositif 100 et le système thermique 100' pourront être associés au dispositif de refroidissement 1 décrit précédemment, l'interface entre le dispositif 100 et la pince cryogénique 10 du laser s'effectuant par l'actionneur 104 pour injecter un flux thermique W dans l'élément 10 et le capteur de température 101. L'actionneur 104 peut être un générateur de flux de chaleur capable de délivrer des puissances positive et négative (actionneur thermique). Par exemple, il peut s'agir d'une résistance électrique parcourue par un courant alternatif. Il peut également s'agir d'un actionneur émettant des ondes électromagnétiques dans le domaine infrarouge. En variante, il peut s'agir d'un actionneur 104 comprenant un générateur de champ magnétique variable associé à un bain thermique isotherme, pour aimanter ou désaimanter adiabatiquement un dispositif raccordé à l'élément 10. L'intérêt de cette désaimantation est de contrôler de manière 20 active la température, avec des temps de réponse courts et surtout, sans perte d'énergie. En pratique, il s'agit de mettre un ou plusieurs matériau(x) à l'intérieur d'une bobine supraconductrice, dont les variations de courant sont commandées par le modèle thermique, pour la régulation de la température. Contrairement à l'emploi d'un moyen de chauffage, cette 25 solution évite toute propagation de chaleur fournie à l'élément et permet de s'affranchir d'éventuels problèmes de déphasage entre le flux thermique injecté dans l'élément et le signal de température ST, problèmes qui peuvent être rencontrés avec un moyen de chauffage classique. Le principe de la désaimantation adiabatique est rappelé ci- 30 dessous. 2 982 014 Si l'on considère par exemple un élément tel que la pince cryogénique 10 du laser « MégaJoule » et le dispositif qui lui est raccordé, le matériau formant la pince est constitué d'une pluralité de moments magnétiques assimilables à autant d'aimants de la taille d'un atome ou d'une 5 molécule. Comme des aimants, les moments magnétiques interagissent entre eux au sein du matériau. Si l'on place la pince cryogénique 10 et le dispositif qui lui est raccordé dans un champ magnétique extérieur (au moyen d'un générateur de champ magnétique par exemple), les moments magnétiques du matériau ro s'alignent sur la direction et le sens du champ appliqué. Plus le champ est intense, plus il y a de moments magnétiques alignés. Le désordre est alors réduit et l'entropie diminue. L'alignement des moments magnétiques est exothermique : c'est l'effet magnétocalorique. 15 Ceci permet d'ajouter ou de retirer un flux de chaleur du dispositif raccordé à la pince cryogénique et par suite, de la pince cryogénique 10 elle-même. C'est le principe de l'aimantation adiabatique. En réduisant le champ magnétique, les moments magnétiques retournent à leur orientation aléatoire, ce qui demande de 20 l'énergie. Les moments magnétiques puisent alors dans l'énergie thermique du matériau et la température du matériau baisse : c'est le principe de la désaimantation adiabatique. Le contrôle du champ magnétique permet finalement de contrôler le flux de chaleur 1N injecté dans la pince cryogénique 10. Il faut cependant noter que la présence d'un moyen 101 pour mesurer la température de la pince cryogénique 10 et celle d'un actionneur 104 pour injecter un flux thermique dans la pince cryogénique 10 peuvent impliquer un retard temporel des différents signaux se propageant dans le dispositif 100 et le système thermique 100'. Ce retard temporel se traduit alors par un déphasage entre ces signaux, qui sont sinusoïdaux. Par suite, le signal de commande de flux thermique Sw arrivant à l'actionneur 104 présente une période et une amplitude adaptées pour compenser les variations périodiques de la température T de l'élément, mais une phase qui ne l'est pas. Il peut en résulter une inefficacité du procédé de régulation. Comme cela sera montré à l'aide d'un exemple de réalisation décrit dans la suite de la description, un retard temporel faible peut être négligé dans le modèle thermique et donc dans le procédé conforme à l'invention. Un tel retard n'engendre en effet aucune instabilité du système de io régulation. Avantageusement, on prendra toutefois en compte ce retard temporel dans le modèle thermique. Pour la prise en compte de ce retard temporel dans le modèle thermique (équation de l'observateur #4), deux cas peuvent alors se présenter selon que le retard temporel généré par certains 15 éléments du système thermique 100' est connu ou pas. 1' cas : le retard temporel est connu Nous avons représenté sur la figure 3 le schéma de 20 fonctionnement du procédé selon l'invention en présence d'un retard connu. Celui-ci est à comparer avec le schéma de la figure 2, pour lequel aucune prise en compte de ce retard n'est effectuée. Comme cela peut être constaté en comparant les figures 2 et 3, le système 100' de la figure 3 comprend un moyen 105 représentatif du 25 retard pur connu de l'ensemble des éléments constitutifs du système. Le signal Sw est ainsi retardé de r périodes d'échantillonnage, au niveau du moyen 105, pour former le signal S*w. Le retard pur 105 n'est qu'une représentation du retard susceptible d'être généré par les différents éléments du système thermique 30 100', notamment l'actionneur 104 et le capteur 101 de température. Pour cette raison, le schéma de la figure 3 intègre ce retard pur au système thermique 100' à réguler. Il s'agit d'un artéfact pour modéliser le retard entre l'entrée et la sortie du système thermique 100'. On comprend que le signal Sw s'exprime en fonction du signal S*,Deq issu du moyen 102 et, le cas échéant, du signal de commande Sc comme cela a été expliqué précédemment. On comprend également que le signal de flux thermique équivalent S*,,,eq issu du moyen 102 intègre également ce retard. Pour cela, il convient d'opérer de la façon suivante, comme cela est représenté sur la figure 4. Le moyen 102 comprend en entrée, le signal de commande de flux thermique Sw et le signal de température ST. Le signal Sw est retardé de r périodes d'échantillonnage dans le moyen 102' afin de former le signal S*w déphasé, ce dernier étant injecté dans l'équation de l'observateur (#4) et suivantes. De la résolution de l'équation de l'observateur ressort le signal de flux thermique équivalent ScDeq. Le signal de flux thermique équivalent Scpec, est alors enregistré dans une mémoire tampon (« buffer » selon la terminologie anglo-saxonne ; non représentée) pour former deux signaux S*(Deconf et S*(Decisup. On peut plus précisément exprimer les signaux S*(Deciindk) et S*GI)eqSup(k) sous la forme : SiàeqInf(k) = Soeq(k - [1/(fp.Ts)+ r}) (#12) S;jegsup(k) = Soeq (k - 1.1/(fp- Ts) + r]) (#13) où -11/(fp. Ts) + r 1 caractérise le retard arrondi à l'entier inférieur du signal ScDeci par rapport au signal S*cDeqInf, 11/(fp. Ts) + r I caractérise le retard arrondi à l'entier supérieur du signal ScDeq par rapport au signal S*(DeqSup, avec : 1, l'opérateur mathématique correspondant à l'arrondi à l'entier inférieur. [ 1, l'opérateur mathématique correspondant à l'arrondi à l'entier supérieur ; r, le nombre de périodes d'échantillonnage dont on retarde le signal, lequel est connu. Les deux signaux S*(DeqInf et S*cl)eqSup permettent alors, grâce à une interpolation linéaire, de reconstituer la valeur approchée du signal S *0eq retardé. On s'affranchit ainsi de la limite de précision donnée par la période d'échantillonnage, ce qui permet de rejeter la perturbation périodique. Finalement, on peut représenter le signal S*(Deq(k) de la manière suivante : 4eq (k) = K1 * S;jeqInf(k) + K2 * S",',eqSup(k) (#14) Avec : K2, le reste de la partie entière de (1/fpTs) + r K1 = 1 - K2; fp, la fréquence du signal perturbation Sp, et Ts, la période d'échantillonnage. En sortie du moyen 102, on obtient un signal de flux thermique équivalent S*cpeq déphasé et d'amplitude différente par rapport au signal de température ST. Ceci s'effectue en deux sous-étapes qui sont associées à l'étape (b), à savoir une première sous-étape dans laquelle on reconstruit le signal Scieq (Flux de chaleur) à partir de la mesure de température T de la pince cryogénique 10 et une deuxième sous-étape dans laquelle on retarde le signal ainsi reconstruit. Au cours de l'étape (b), on retarde donc également le signal de flux thermique équivalent préalablement déterminé, par rapport à la phase du signal de température déterminé à l'étape (a). 2ème cas : le retard temporel n'est pas connu Nous avons représenté sur la figure 5 le fonctionnement du io procédé selon l'invention en présence d'un retard qui n'est pas connu. En particulier, le schéma représenté sur la figure 4 est à comparer avec les schémas représentés sur les figures 2 (non prise en compte du retard) et 3 (retard connu). Dans ce cas où le retard lié aux différents éléments 15 constitutifs du système thermique 100' n'est pas connu, celui-ci doit donc être estimé et pris en compte pour reconstruire le flux de chaleur équivalent Sk(Deq. Il existe plusieurs méthodes pour identifier un retard temporel dans des systèmes linéaires. On peut par exemple se référer au document « Survey and Comparison of Time-Delay Estimation Methods in Linear 20 Systems », Svante Beirklund, PhD thesis, 2003 qui présente certaines de ces méthodes. Dans le cas d'espèce, une méthode dite approximative dans le domaine fréquentiel a été choisie. Pour cela, on utilise un signal d'identification du déphasage 25 généré par le système thermique. Ce signal d'identification est fourni par un moyen 106 adapté, tel qu'un générateur de fonctions, et présente une forme proche du signal représentatif du fonctionnement du cryoréfrigérateur, par exemple un tube à gaz pulsé, et donc du signal de perturbation. En l'occurrence, il s'agit donc d'un signal Sm de flux thermique sinusoïdal 30 s'exprimant sous la forme : Swi(k) = Wi. sin (27r. fi. Ts. k) (#15) où: f, est la fréquence du signal d'identification qui est proche de la fréquence de fonctionnement du tube à gaz pulsé, s Ts est la période d'échantillonnage de ce signal, et W, est l'amplitude du signal. Le signal d'identification Sw, est injecté dans le système thermique 100', qui comprend des composants susceptibles de générer un il) retard, à savoir notamment la pince cryogénique 10, l'actionneur 104 et le capteur 101. En sortie de la pince cryogénique 10, à savoir à l'endroit où l'on mesure la température T de la pince 10 qu'on cherche à réguler, le signal d'identification est déphasé d'une phase 4, par rapport au signal provenant du générateur de fonctions. 15 Le déphasage 4) est déterminé dans le moyen 107, en effectuant par exemple une analyse spectrale. Pour cela, le signal d'identification Sw, est injecté depuis le générateur de fonctions vers l'élément 107. En parallèle, le signal ST, qui est récupéré en sortie de la pince 10, est transmis au moyen 107. 20 Pour déterminer ce déphasage 4, le moyen 107 peut effectuer cette analyse spectrale au moyen d'un filtre numérique basé sur l'algorithme de Goertzel. Cet algorithme est connu de l'homme du métier et on pourra obtenir plus d'informations à son sujet dans le document « Goertzel's 25 algorithm », Douglas L. Jones, 2006. La fonction de transfert H entre la sortie du filtre et l'entrée du filtre s'écrit sous la forme : H(Z) = (1-2.cos(2.7r.n.z-l+z-2) 1 (#16) 30 Où: où f est la fréquence que l'on souhaite extraire (analyse spectrale) et z la variable numérique du filtre. Pour déterminer le déphasage phi = (1); on calcule la fonction de transfert du filtre de Goertzel sur N échantillons. Puis, à partir des deux dernières valeurs du filtre (dl, d2), les parties imaginaire (« imag ») et réelle (« real ») sont calculées. On en déduit enfin la phase recherchée, en calculant l'arc tangente (« atan ») du rapport entre la partie imaginaire et la partie réelle. Le programme est le suivant : realW = 2.0*cos(2.0*pi*k/N); imagW = sin(2.0*pi*k/N); dl = 0.0; d2 = 0.0; for (n=0; n En comparant le déphasage (1) déterminé par analyse spectrale au déphasage 4)o déterminé par le modèle thermique présenté à l'appui des relations (#1) à (#11) en l'absence de tout retard, on en déduit que la différence cp, = ti) - (1)0 est un déphasage lié effectivement au retard temporel généré par le système thermique 100'. Il faut ensuite exprimer le déphasage cp, sous la forme d'un retard temporel. Ceci s'effectue en prenant au moins en compte la valeur de ce déphasage cpr et la fréquence fi du signal d'identification. De préférence, on exprime ce retard temporel en un nombre de périodes d'échantillonnage, ce qui fait intervenir la période d'échantillonnage Ts de ce signal. Dans ce cas, à ce déphasage cp, , on peut associer un nombre r de périodes d'échantillonnage de retard, sous la io forme : (Pr (#17) = 27r. f i.Ts Ce nombre de périodes d'échantillonnage est calculé dans le 15 moyen 108. La détermination du nombre r s'effectue en continu et peut présenter une certaine instabilité, du fait que le déphasage déterminé dans le moyen 107 avec l'algorithme de Goertzel n'est pas forcément stable. Pour remédier à cet éventuel inconvénient, le moyen 108 20 effectue avantageusement un filtrage non-linéaire du nombre r . Un exemple de filtrage non-linéaire pouvant être mis en oeuvre est représenté sur la figure 6. La figure 6 représente en abscisse, le nombre r obtenu avec l'équation (#17) et en ordonnée, le nombre r de périodes d'échantillonnage, 25 une fois filtré. On observe que le filtre non-linéaire effectue une correction nommée A. Cette correction A est évaluée de la manière suivante. On sait que le système thermique 100' est robuste à des variations de retard de deux périodes d'échantillonnage. On choisit alors une 30 correction A égale à 1,5 période d'échantillonnage. Il s'agit d'un compromis 2 9 82014 28 entre la réaction du système thermique 100' et la robustesse. Avec cette approche, les valeurs erronées du retard temporel ne seront pas prises en compte. Le nombre r de périodes d'échantillonnage représentatif du 5 retard temporel ainsi estimé est exploité pour le calcul du flux de chaleur équivalent S*(Deq dans le moyen 102. A cet effet, on peut reprendre le descriptif qui a été réalisé précédemment à l'appui du schéma de la figure 4.Cependant, dans ce cas, le moyen 102 comprend en entrée, le signal de commande de flux thermique Sw, le signal de température ST associée à la température mesurée de la pince 10 et ce nombre r de périodes d'échantillonnage, estimé. Il suffit alors de modifier le schéma de la figure 4 en y ajoutant une entrée supplémentaire, connectée au moyen 102'. Bien entendu, si aucun filtrage linéaire conforme à celui de la figure 6 n'est effectué, c'est le nombre T. qui est injecté dans le moyen 102. La figure 5 permet de préciser ce qui se passe dans le moyen 102. On note que la méthodologie employée permet de revenir à une situation analogue à celle du schéma représenté sur la figure 3 (laquelle précise ce qui se passe dans le moyen de calcul 102), où le signal de flux de commande Sw est retardé de r périodes d'échantillonnage avant d'être injecté dans l'équation de l'observateur. Le nombre r de périodes d'échantillonnage est mis à profit pour retarder le signal Sw dans le moyen 102', qui effectue un retard pur pour fournir le signal déphasé S*w. Ce dernier signal (équivalent au flux thermique W) est injecté dans l'équation de l'observateur. Ainsi, en sortie du moyen 102, on obtient un signal de flux thermique équivalent S*(Deq reconstruit déphasé qui peut être injecté en opposition de phase dans l'actionneur 104 et par suite, dans la pince 10. Avantageusement, on forme plus précisément un signal Sw = Sc - S*,Deq qui sera transmis au système thermique 100'. Exemple d'application Afin de mettre en évidence la pertinence de l'invention, une simulation numérique a été effectuée pour la pince du laser « MégaJoule ». Cette pince est représentée sur la figure 7. Par ailleurs, un retard temporel connu, lié au comportement dynamique des différents éléments constitutifs de la pince 10 est pris en considération. Ce retard temporel intègre donc notamment les retards liés io aux capteurs de température et aux actionneurs qui génèrent physiquement le flux thermique. En effet, le modèle considère que trois capteurs de température T1, T2, T3 sont prévus respectivement sur l'échangeur 1030, l'enclume 1020 et l'embase 1010. Par ailleurs, il considère que deux 15 actionneurs sont pris en compte, respectivement au niveau de l'embase 1010 (Wemb) et au niveau de l'échangeur 1030 (Wech). Une représentation montrant l'implantation des capteurs de température et des actionneurs est fournie sur la figure 8. Les conditions choisies sont les suivantes. 20 Les coefficients de la matrice Csin sont connus. Les coefficients de la matrice Asin sont déterminables à partir de la fréquence fp et de la période d'échantillonnage T5. En l'occurrence, la fréquence de fonctionnement fp du tube à gaz pulsé 2 est de fp = 1,4Hz, le signal étant sinusoïdal (cela correspond à une période de l'ordre de 714ms) et 25 la période d'échantillonnage vaut, pour cette simulation, Ts = 80ms. Par ailleurs, le système d'équations (#1) s'établit en effectuant une analogie thermique-électrique sur les différents éléments de la pince 10, en considérant les trois capteurs de température et les deux actionneurs. L'analogie thermique-électrique est bien connue de l'homme 30 du métier. On rappelle cependant qu'une température peut être associée à un potentiel électrique, une capacité thermique (laquelle correspond au produit de la capacité thermique massique d'un objet considéré par sa masse) à une capacité électrique, un flux de chaleur à un courant électrique et une résistance thermique à une résistance électrique. Le schéma électrique équivalent au schéma thermique de la figure 8, est représenté sur la figure 9. Sur cette figure 9, l'indice 1 correspond à l'embase 1010, l'indice 2 à l'enclume 1020 et l'indice 3 à l'échangeur 1030. L'embase 1010 et l'enclume 1020 sont considérées comme étant réalisées en aluminium et l'échangeur est réalisé en cuivre. Par ailleurs, on considère que la masse de l'embase 1010 est de 218g, la masse de l'enclume 1020 de 143g et la masse de l'échangeur 1030 de 692g. En fonction de la nature des matériaux, des hypothèses simplificatrices peuvent être mises en oeuvre. Par ailleurs, il convient de prendre en compte l'aspect dynamique du filtrage actif opéré par l'invention, l'analogie thermique- électrique ne permettant pas à elle seule de définir les matrices A, B et C. Ceci relève du savoir-faire du demandeur. Cependant, aux fins de mise en oeuvre de l'invention, nous précisons numériquement, et pour le cas d'espèce, ces matrices : A=F0,8187 0'0677 0,0008 0 0 0,8714 0,0299 0,1993 0 0,1285 0,9670 0 0 0 0 0,8007 0,1813 B= 0 C = 0 0 0] L'amplitude de la température T de la pince cryogénique 10 en l'absence de la régulation effectuée par le procédé selon l'invention est fixée à 5mK. Cette valeur est réaliste par rapport aux données expérimentales actuellement disponibles sur le dispositif représenté sur la figure 1 pour le laser « MégaJoule ». On comprend donc que l'amplitude de 5mK est celle qui est actuellement obtenue sur la pince du laser « MégaJoule » avec la mise en oeuvre du stabilisateur 3 et du circuit de refroidissement 4, mais sans la régulation effectuée par l'invention. Le retard temporel Tr lié aux différents éléments constitutifs du dispositif 100 est fixé à Tr = 195ms. Dans ce cas, on a K2 = 0,366 (fp= 1,4Hz ; Ts = 80ms ; r = Tr/Ts = 2,4375) et K1 = 0,634. Les résultats obtenus sont représentés sur la figure 10. La figure 10 montre l'évolution de la température T (en ir) Kelvins) de la pince cryogénique 10 en fonction du temps t (en secondes). Entre t = Os et t = 25s, le filtrage selon l'invention n'est pas mis en oeuvre. On observe donc une variation périodique de cette température, dont l'amplitude vaut 5mK (amplitude crête à crête de 10mK) autour d'une valeur de référence de Tréf = 17,98K. 15 Au temps t = 25s, le filtrage actif selon l'invention est mis en oeuvre. On observe alors que la variation de température diminue nettement, l'amplitude de la température T est de l'ordre du mK, autour de cette valeur de référence. Puis, au temps t = 50s, on impose à la pince cryogénique 10 20 une augmentation de la température T de référence, celle-ci étant portée à Tréf = 18K. On observe alors que la variation de température T de la pince cryogénique 10 est de l'ordre du mK autour de cette autre valeur de référence, après un transitoire d'une dizaine de secondes. 25 Comme cela a été précisé précédemment, la mise en oeuvre d'un retard temporel dans le modèle thermique n'est pas obligatoire si celui-ci est « faible ». C'est ce que montre la figure 11. La figure 11 représente l'atténuation des variations de la 30 température T de la pince cryogénique 10, en fonction du retard temporel effectivement généré. Une première courbe Cl montre l'évolution de cette atténuation, en prenant en compte ce retard temporel, connu. Une deuxième courbe C2 montre l'évolution de cette atténuation en l'absence de prise en compte de ce retard temporel dans la régulation effectuée dans le cadre de l'invention. On observe que l'atténuation est comparable à celle de la courbe Cl jusqu'à une valeur de retard temporel de 130ms. Au-delà, il devient nécessaire de prendre en compte le retard temporel dans le modèle thermique. Avantageusement, on prendra donc en compte ce retard temporel dans le modèle thermique et donc dans le procédé selon l'invention, mais, dans certains cas, cette prise en compte peut être négligée. Résultats d'un essai : Afin de confirmer la pertinence de l'invention, un essai sur un système réel a été effectué. Ce système de test est similaire au système comportant la pince du laser « MégaJoule ». Il est représenté sur la figure 12. Contrairement à ce qui est représenté sur la figure 12, ce système de test ne dispose pas de l'embase 1010 destinée à porter la cible 11, et en conséquence il ne dispose pas de la cible 11 elle-même. Il ne dispose pas non plus de l'échangeur 1030, dans la mesure où le moyen 4 est un barreau qui a pour fonction de représenter le circuit de refroidissement de la figure 1. Toutefois, l'embase 1010 et la cible 11 qu'elle porte sont représentées sur la figure 12 pour indiquer au lecteur où elles sont destinées à être installées au sein du système. Par ailleurs, le système est pourvu de deux tubes 2 à gaz pulsé. Dans le cas présent, on cherche donc à réguler la température de l'enclume 1020, laquelle est munie d'un capteur de température. Un actionneur thermique (non représenté) est également prévu sur l'enclume 1020, à proximité du capteur de température. Les conditions choisies sont les suivantes : - les coefficients de la matrice Cs;, sont connus ; s - les coefficients de la matrice As,n sont déterminables à partir de la fréquence fp et de la période d'échantillonnage T. En l'occurrence, la fréquence de fonctionnement fp des tubes à gaz pulsé 2 est de fp = 1,412Hz, le signal étant sinusoïdal (cela correspond à une période de l'ordre de 708,2 ms). La période d'échantillonnage vaut, pour cet essai, Ts 10 = 40ms. Tout comme dans I' « exemple d'application » présenté précédemment, une analogie thermique-électrique est effectuée, en prenant en compte les éléments suivants : stabilisateur 3, barreau 4 et l'enclume 1020. 15 Le stabilisateur 3, le barreau 4 et l'enclume 1020 sont tous réalisés en Cuivre. Le stabilisateur 3 présente une masse de 10,892kg, le barreau 4 de 5,92kg et l'enclume 1020 de 1,74kg. Le schéma électrique équivalent est fourni en figure 13. Sur cette figure 13, l'indice 1 correspond à l'enclume 1020, l'indice 2 au barreau 4 20 et l'indice 3 au stabilisateur 3. Là encore, des hypothèses simplificatrices peuvent êtres prévues. Par ailleurs, il convient de prendre en compte l'aspect dynamique du filtrage actif opéré par l'invention, l'analogie thermique- 25 électrique ne permettant pas à elle seule de définir les matrices A, B et C. Ceci relève du savoir-faire du demandeur. Aux fins de mise en oeuvre de l'invention, les matrices A, B et C sont cependant fournies ci-dessous, dans le cas d'espèce : -9.1560 9.1560 0 A = 2.6911 -4.8776 2.1865 0 1.1884 -1.4601 0.5951 B= 0 0 C = [1 0 0] L'amplitude de la température T de la pince cryogénique 10 en l'absence de la régulation effectuée par le procédé selon l'invention est fixée à environs 10mK. On comprend donc que l'amplitude de 10mK est celle qui est actuellement obtenue sur la pince du laser « MégaJoule » avec la mise en oeuvre du stabilisateur 3 et du remplacement circuit de refroidissement à l'hélium 4 par un barreau de cuivre, mais sans la régulation effectuée par l'invention. Le retard temporel Tr lié aux différents éléments constitutifs du io dispositif 100 est mesuré à environ Tr = 45 ms. Ce retard est donc un retard connu. Dans ce cas, on a K2 = 0,83 (fp= 1,412Hz ; Ts = 40ms ; r = Tr/Ts = 1,125) et K1 = 0,17. Les résultats obtenus sont représentés sur la figure 14. 15 La figure 9 montre l'évolution de la température T (en Kelvins) de la pince cryogénique 10 en fonction du temps t (en secondes). Le régime étant déjà établi, on remarque que les oscillations rapides à 1,4 Hz sont pratiquement rejetées ( REVENDICATIONS1. Procédé de régulation de la température d'un élément (10) refroidi par un cryoréfrigérateur (2) à fonctionnement périodique, caractérisé en ce s qu'il comprend les étapes suivantes : (a) mesurer au moins l'évolution temporelle périodique de la température T de l'élément (10), cette évolution permettant d'obtenir un signal de température représentatif du fonctionnement périodique du cryoréfrigérateur (2) ; 10 (b) construire un signal de flux thermique périodique (Se.eq, S*ceq) dit équivalent, en fonction du signal de température obtenu à l'étape (a) ; (c) construire un signal de commande de flux thermique (Sw), dépendant du signal de flux thermique équivalent construit à l'étape (b), ledit signal de commande de flux thermique (Sw) servant également, 15 par retour, à la construction du signal de flux thermique équivalent (S,Deq, S*,Deq) lors de l'étape (b); (d) injecter un flux thermique (W) dans l'élément (10), sur la base du signal de flux thermique (Sw) construit à l'étape (c), pour compenser les variations périodiques mesurées de la température T de l'élément. 20 2. Procédé selon la revendication 1, dans lequel le signal de commande de flux thermique (Sw) dépend également d'un signal de commande (Sc) pour améliorer la compensation des variations périodiques de la température de l'élément aux basses fréquences. 25 3. Procédé selon la revendication précédente, dans lequel le flux thermique (W) injecté dans l'élément à l'étape (c) dépend de la différence entre le signal de commande (Sc) et le signal de flux thermique équivalent (Scteq, S*(Deq) construit à l'étape (b), pour limiter la 30 puissance thermique injectée dans l'élément qui est nécessaire à la compensation des variations périodiques de sa température. 4. Procédé selon l'une des revendications précédentes, dans lequel l'étape (b) comprend également une étape additionnelle consistant à retarder le signal de flux thermique équivalent (S,Deq) préalablement construit, par rapport à la phase du signal de température déterminé à l'étape (a). 5. Procédé selon la revendication précédente, dans lequel il est prévu une étape consistant à estimer la valeur de ce retard temporel. 6. Procédé selon la revendication précédente, dans lequel l'estimation de ce retard temporel comprend les étapes suivantes : (bi) injecter un signal d'identification du retard temporel dans un système thermique (100') comprenant l'élément dont on cherche à réguler la température et un dispositif (100) de mise en oeuvre du procédé, ce signal d'identification présentant une forme et une fréquence proches du signal représentatif du fonctionnement du cryoréfrigérateur; (b2) déterminer le déphasage (1) du signal d'identification en sortie du système thermique, ce déphasage étant associé au retard temporel qu'on cherche à estimer ; (b3) calculer le déphasage (1)0 de ce signal d'identification avec un modèle représentant le comportement thermique du système thermique (100') en l'absence de tout retard; (b4) calculer le déphasage (p, associé au retard temporel par la relation (Pr_ - (I)o; (b5) estimer le retard temporel, au moins à partir du déphasage (pr calculé à l'étape (b4) et de la fréquence du signal d'identification, ledit retard étant utilisé pour retarder le signal de flux thermique équivalent. 7. Dispositif de mise en oeuvre du procédé de régulation de la température d'un élément (10) refroidi par un cryoréfrigérateur (2) àfonctionnement périodique selon l'une des revendications précédentes, caractérisé en ce qu'il comprend : au moins un moyen (101) pour mesurer l'évolution périodique de la température T de l'élément (10) liée au fonctionnement périodique du cryoréfrigérateur (2); - au moins un moyen (102, 103, 106, 107, 108) pour déterminer le flux thermique équivalent (S(Deq, S*(Deq) associé à cette variation périodique de la température T de l'élément (10), ce moyen ou cet ensemble de moyens comprenant un ordinateur ou un ensemble d'ordinateurs ou un processeur ou un ensemble de processeurs; - au moins un moyen (104) pour injecter un flux thermique (W) dans l'élément (10), à partir d'un signal de commande de flux thermique (Sw) dépendant du signal de flux thermique équivalent (Soeq, S*(peq) ; - un moyen (109) pour injecter le signal de commande de flux thermique (Sw) au niveau de l'entrée dudit au moins un moyen (102) pour déterminer le flux thermique équivalent (Scpeq, S*(p q). 8. Dispositif selon la revendication précédente, dans lequel il est prévu un moyen (103) capable de délivrer un signal de commande (Sc) pour améliorer la compensation des variations périodiques de la température de l'élément aux basses fréquences. 9. Dispositif selon l'une des revendications 7 ou 8, dans lequel le moyen (104) pour injecter un flux thermique dans l'élément (10) comprend un générateur de champ magnétique variable associé à un bain thermique isotherme, pour aimanter ou désaimanter adiabatiquement un dispositif raccordé à l'élément (10). 10. Dispositif selon l'une des revendications 7 à 9, dans lequel ledit au moins un moyen (101) pour déterminer l'évolution périodique de la température de l'élément (10) et ledit au moins un moyen (104) pourinjecter un flux thermique (W) dans l'élément (10) sont disposés l'un à côté de l'autre. 11. Installation cryogénique pour refroidir un élément (10) avec un cryoréfrigérateur à fonctionnement périodique (2), ladite installation comprenant un stabilisateur (3) relié au cryoréfrigérateur (2) et un circuit de refroidissement (4) comprenant un premier échangeur thermique (41) relié au stabilisateur (3) et un deuxième échangeur thermique (42) relié à l'élément (10), caractérisé en ce qu'il comprend, io en outre, un dispositif selon l'une des revendications 7 à 10.