Quantitative 3D non-linear simulations of shattered pellet injection in ASDEX Upgrade using JOREK

Cet article présente des simulations 3D non linéaires quantitatives de l'injection de pellets brisés dans le tokamak ASDEX Upgrade à l'aide du code JOREK, où l'intégration d'une limitation simplifiée du flux de chaleur parallèle permet de résoudre les écarts avec l'expérience et d'obtenir des prédictions fiables pour la validation des modèles de mitigation des disruptions en vue d'ITER.

L'essentiel

🌟 Le Grand Défi : Arrêter l'orage nucléaire

Imaginez que vous essayez de faire bouillir de l'eau dans une casserole, mais au lieu d'un simple feu, vous avez un orage électrique gigantesque à l'intérieur. C'est ce qui se passe dans un tokamak (une machine à fusion nucléaire comme celle d'ITER). Parfois, cet orage devient incontrôlable : c'est une rupture majeure (ou disruption).

Si rien ne se passe, cet orage frappe les murs de la machine comme un marteau-piqueur, risquant de la détruire. Pour éviter cela, les scientifiques ont une solution de secours : l'injection de pellets éclatés (des petits morceaux de glace).

L'analogie de la tempête de neige :
Imaginez que votre four surchauffe. Pour le refroidir, vous lancez des milliers de petits flocons de neige à l'intérieur. Ces flocons fondent, absorbent la chaleur et la transforment en vapeur (rayonnement), ce qui éteint le feu doucement au lieu de le laisser exploser. C'est le principe de l'injection de pellets.

🧊 Le Problème : La "Neige" fondait trop vite

Dans les expériences passées sur la machine ASDEX Upgrade (en Allemagne), les scientifiques ont lancé ces pellets de glace. Mais quand ils ont simulé la scène sur ordinateur avec le logiciel JOREK, quelque chose clochait :

• La réalité : La glace met un certain temps à fondre et à refroidir le four.
• L'ordinateur : Il disait que la glace fondait instantanément et que le four se refroidissait trop vite, comme si la neige se transformait en eau bouillante en une seconde.

Pourquoi ? Parce que le modèle informatique supposait que la chaleur se déplace dans le plasma (le gaz super chaud) comme si c'était de l'eau qui coule dans un tuyau lisse. En réalité, le plasma est plus complexe : la chaleur a du mal à se déplacer le long des lignes magnétiques, un peu comme si le tuyau était bouché ou très rugueux.

🛠️ La Solution : Mettre un "Frein" à la chaleur

Les chercheurs ont compris qu'ils devaient ajouter un frein à leur simulation. Ils ont décidé de réduire la vitesse à laquelle la chaleur voyage le long des lignes magnétiques par un facteur de 10.

L'analogie du trafic routier :

• Avant (Modèle ancien) : C'était comme si tous les voitures (la chaleur) roulaient sur une autoroute vide à 300 km/h. Elles arrivaient trop vite sur le lieu de l'accident (le pellet).
• Après (Nouveau modèle) : Ils ont mis des ralentisseurs et des bouchons. Les voitures roulent maintenant à 30 km/h. La chaleur arrive plus lentement, ce qui correspond exactement à ce que les scientifiques ont vu dans la vraie expérience.

Grâce à ce petit ajustement (diviser la conductivité thermique par 10), la simulation est devenue quantitativement exacte. Elle prédit maintenant la durée du refroidissement et la quantité de lumière émise avec une précision incroyable.

🔍 Ce qu'ils ont appris avec ce nouveau modèle

Une fois le "frein" en place, ils ont pu tester différentes stratégies pour voir laquelle est la meilleure pour protéger la machine :

1. La quantité de "neige" (Néon) :

   • Ils ont testé des pellets avec très peu de néon (trace) et beaucoup de néon (10%).
   • Résultat : Avec peu de néon, le refroidissement est lent et doux (comme une pluie fine). Avec beaucoup de néon, c'est plus rapide (comme un déluge). Les deux fonctionnent, mais le temps de réaction change. Le modèle correspond maintenant parfaitement aux temps observés en laboratoire.

2. La taille des "flocons" (Taille des fragments) :

   • Ils ont comparé des gros morceaux de glace à des milliers de tout petits grains de poussière.
   • Résultat surprenant : On pensait que les petits grains fondraient mieux car ils ont plus de surface totale. Et c'est vrai au début ! Mais les gros morceaux restent dans le cœur de la machine plus longtemps, permettant de refroidir l'ensemble plus efficacement.
   • Le mystère des petits grains : En expérience, les très petits grains semblent parfois ne pas pénétrer assez loin. Les chercheurs pensent que c'est à cause de l'effet fusée : quand un petit grain fond, le gaz qui s'échappe le propulse comme une fusée, le faisant rebondir hors du cœur du plasma. Leur modèle actuel ne le prend pas encore en compte, mais c'est la prochaine étape !

🚀 Pourquoi c'est important pour le futur ?

Cette recherche est cruciale pour ITER, le futur réacteur nucléaire international qui va essayer de produire de l'énergie propre.

• Avant : Les simulations étaient comme des dessins au crayon : on voyait la forme, mais pas les détails précis.
• Maintenant : Grâce à ce "frein à la chaleur", les simulations sont comme des photographies haute définition.

Cela permet aux ingénieurs de dire : "Si on lance ce type de pellet, à telle vitesse, avec telle taille, on sauvera la machine à 100 %." C'est un pas de géant vers la maîtrise de l'énergie des étoiles sur Terre.

En résumé : Les scientifiques ont corrigé un bug dans leur calcul de la chaleur, ce qui a rendu leurs prédictions parfaites. Ils savent maintenant exactement comment lancer la "neige" pour éteindre le "feu" nucléaire sans abîmer la machine.

Résumé technique

Titre de l'étude

Simulations quantitatives 3D non linéaires de l'injection de pellets éclatés (SPI) dans ASDEX Upgrade utilisant JOREK.

1. Problématique

Les disruptions majeures dans les tokamaks de grande taille, tels qu'ITER, posent des défis critiques pour la durée de vie des composants face au plasma et l'intégrité de la machine. L'injection de pellets éclatés (Shattered Pellet Injection - SPI) est la stratégie de mitigation privilégiée pour ITER, visant à injecter massivement du matériau dans le plasma pour rayonner l'énergie uniformément et supprimer la formation d'électrons de fuite (runaway electrons).

Cependant, une modélisation prédictive précise est essentielle pour optimiser les paramètres d'injection. Les travaux précédents de l'équipe (Tang et al., 2025) utilisant le code MHD 3D non linéaire JOREK sur des expériences ASDEX Upgrade (AUG) ont révélé une divergence significative : la phase pré-thermal quench (pré-TQ) simulée était beaucoup plus courte que celle observée expérimentalement. Cette erreur était attribuée à une surestimation du transport de chaleur parallèle aux lignes de champ magnétique, modélisé par la formule classique de Spitzer-Härm, qui ne tient pas compte des effets de limitation du flux de chaleur (heat-flux limiting) dans les régimes de refroidissement rapide induits par le pellet.

2. Méthodologie

Les auteurs ont revisité les simulations SPI d'AUG (décharge #40355 et #40673) en modifiant l'approche de modélisation du transport thermique parallèle ($\chi_\parallel$) :

• Code de simulation : Utilisation du code MHD réduit à deux températures JOREK.
• Configuration : Plasma en mode H, champ toroïdal $B_t = 1.8$ T, courant plasma $I_p = 0.8$ MA.
• Modèle de transport thermique : Au lieu d'utiliser la valeur nominale de Spitzer-Härm ($\chi_{\parallel, SH}$), les auteurs ont appliqué une réduction constante d'un facteur 10 ($\chi_\parallel = 0.1 \cdot \chi_{\parallel, SH}$). Cette approche simplifiée vise à simuler les effets de limitation du flux de chaleur (où le flux est saturé en dessous de la limite de libre parcours moyen, $q_\parallel \sim \alpha n_e T_e v_{th}$).
• Paramètres étudiés :
  • Effet de la réduction de $\chi_\parallel$ sur la dynamique du TQ.
  • Influence de la fraction de néon (0,12 % vs 10 %).
  • Influence de la taille des fragments (53 fragments vs 1105 fragments), en utilisant des distributions de vitesse et de taille basées sur le modèle de fragmentation de Parks.
• Comparaison : Les résultats sont comparés quantitativement aux données expérimentales d'AUG, notamment l'évolution de l'énergie thermique, le courant plasma et la fraction de rayonnement.

3. Contributions Clés

• Validation Quantitative : La principale contribution est la transition d'une comparaison qualitative à une validation quantitative des simulations SPI. La réduction du coefficient de diffusivité thermique parallèle permet de reproduire fidèlement la durée de la phase pré-TQ et la fraction de rayonnement observées expérimentalement.
• Identification des mécanismes physiques : L'étude démontre que la surestimation précédente de $\chi_\parallel$ entraînait un réchauffement trop rapide du nuage de neutres ablaté, accentuant les pertes induites par la dérive du plasmoid (plasmoid-drift) et amplifiant artificiellement la réponse MHD (instabilités de déchirure résistives $2/1$ et $3/2$).
• Optimisation des paramètres SPI : Réévaluation de l'impact de la fraction de néon et de la taille des fragments avec un modèle physique plus réaliste, fournissant des données cruciales pour la conception du système de mitigation d'ITER.

4. Résultats Principaux

A. Impact de la diffusivité thermique réduite

• Durée du pré-TQ : Avec le facteur de réduction de 0,1, la durée du pré-TQ simulée correspond quantitativement aux expériences (environ 1 ms pour 10 % de néon et 3-4 ms pour 0,12 % de néon), contrairement aux simulations précédentes qui étaient trop rapides.
• Dynamique MHD : La réduction de $\chi_\parallel$ atténue la réponse MHD. Les surfaces de flux magnétiques restent plus intactes au cœur du plasma, et la propagation du front de refroidissement est mieux alignée avec la trajectoire des fragments, évitant un effondrement thermique prématuré.
• Fraction de rayonnement : Les simulations avec $\chi_\parallel$ réduit montrent une fraction de rayonnement beaucoup plus élevée et réaliste, indiquant que la perte d'énergie thermique se fait davantage par rayonnement que par conduction/convection excessive.

B. Influence de la fraction de néon

• Le comportement de refroidissement en deux étapes (d'abord conduction/convection, puis rayonnement) persiste.
• Une fraction de néon plus faible (0,12 %) entraîne un processus de TQ plus lent, une durée d'ablation plus longue et une assimilation d'électrons plus élevée, ce qui pourrait être bénéfique pour la suppression des électrons de fuite.

C. Influence de la taille des fragments

• Petits fragments (1105) : Surface totale plus grande, taux d'ablation initial plus élevé, mais durée de pré-TQ plus courte.
• Gros fragments (53) : Durée de pré-TQ plus longue, permettant une assimilation totale légèrement supérieure à la fin du TQ.
• Limitation du modèle : Pour les traces de néon, l'assimilation des petits fragments est sous-estimée par rapport à l'expérience. Les auteurs suggèrent que l'absence du effet fusée (rocket effect) dans le modèle actuel est la cause : les petits fragments, ayant une masse faible, pourraient être expulsés plus rapidement du plasma par cet effet, limitant leur pénétration.

5. Signification et Perspectives

Ce travail est une étape cruciale pour la fiabilité des modèles de mitigation de disruptions pour ITER.

• Confiance accrue : L'intégration d'un traitement simplifié de la limitation du flux de chaleur permet désormais d'utiliser JOREK pour des études prédictives fiables sur les configurations d'injection SPI.
• Optimisation pour ITER : Les résultats aident à définir les paramètres optimaux (taille des fragments, composition, timing) pour minimiser les charges thermiques et les forces électromagnétiques sur les composants d'ITER.
• Travaux futurs : Les auteurs prévoient d'implémenter des modèles plus sophistiqués, notamment un traitement plus réaliste de la conduction parallèle limitée et l'ajout de l'effet fusée pour améliorer la modélisation de la pénétration des petits fragments.

En conclusion, cette étude démontre que la prise en compte correcte de la physique du transport de chaleur parallèle est indispensable pour simuler avec précision les disruptions induites par SPI, transformant JOREK en un outil de validation expérimentale robuste.