Hydrodynamic simulations of expanded warm dense foil heated by pulsed-power

Ce document présente un cadre de modélisation robuste qui couple des simulations de circuits électriques à puissance pulsée avec des codes hydrodynamiques unidimensionnels afin de concevoir et d'optimiser avec précision des expériences pour la génération de matière chaude dense expansée dans des feuilles métalliques minces confinées dans des cellules de saphir.

L'essentiel

La vue d'ensemble : Qu'est-ce que la « matière dense chaude » ?

Imaginez un matériau qui n'est ni tout à fait un solide, ni tout à fait un liquide, ni tout à fait un gaz. Il se trouve dans un entre-deux étrange et désordonné appelé matière dense chaude (MDH). Voyez cela comme une piste de danse bondée où tout le monde bouge rapidement (chaud) mais continue de s'entrechoquer (dense).

Les scientifiques doivent comprendre cet état de la matière pour étudier des choses comme l'intérieur des planètes géantes ou pour construire de meilleurs réacteurs à fusion nucléaire. Mais c'est difficile à étudier car c'est complexe à créer en laboratoire et encore plus difficile à prédire avec les mathématiques.

L'expérience : La « crêpe électrique »

Les chercheurs ont mis en place une expérience pour créer cet état.

• La configuration : Ils ont pris une feuille de métal très mince (comme une feuille d'aluminium microscopique) et l'ont prise en sandwich entre deux plaques épaisses et dures en saphir (comme le verre d'une montre).
• L'action : Ils ont bombardé ce sandwich métallique avec une décharge électrique massive et ultra-rapide (puissance pulsée).
• Le résultat : L'électricité chauffe le métal si rapidement (en moins d'un millionième de seconde) qu'il fond, bout et se transforme en un plasma chaud et en expansion. Comme les plaques de saphir le retiennent, le métal ne peut s'étendre que dans une seule direction, comme une crêpe qui gonfle.

Le problème : La « boîte noire »

Le défi est que lorsque vous bombardez le métal, deux choses se produisent simultanément :

1. Le circuit électrique : L'électricité circule à travers les fils, l'interrupteur et le métal. À mesure que le métal chauffe et change de forme, sa capacité à conduire l'électricité change, ce qui modifie le flux de courant.
2. Le mouvement physique : Le métal chauffe, se dilate et se déplace. En se déplaçant, il change la forme du circuit, ce qui modifie à nouveau l'électricité.

C'est une boucle de rétroaction. Si vous essayez de calculer l'électricité sans savoir comment le métal bouge, vous vous trompez. Si vous essayez de calculer le mouvement du métal sans connaître l'électricité, vous vous trompez aussi.

La solution : Le modèle du « vélo tandem »

Les auteurs ont construit un programme informatique qui agit comme un vélo tandem.

• Le cycliste 1 (Le modèle électrique) : Cette partie simule le générateur de puissance, l'interrupteur et les fils. Elle calcule la quantité de courant qui circule.
• Le cycliste 2 (Le modèle hydrodynamique) : Cette partie simule la feuille de métal. Elle calcule comment le métal chauffe, se dilate et change de densité.

Ces deux cyclistes sont liés ensemble. Chaque fraction de seconde, ils communiquent :

• « Hé, le métal vient de devenir plus chaud et plus fin », dit le Cycliste 2.
• « D'accord, je vais ajuster le flux de courant car le métal est maintenant un moins bon conducteur », dit le Cycliste 1.
• « D'accord, je vais mettre à jour la chaleur et la pression en fonction de ce nouveau courant », dit le Cycliste 2.

Comment ils l'ont testé

Pour s'assurer que leur « vélo tandem » fonctionne, ils l'ont testé de trois manières différentes, comme on vérifierait le moteur d'une voiture à différents niveaux :

1. Le test du « courant connu » : Ils ont injecté dans l'ordinateur les mesures de l'électricité réelle de l'expérience réelle et ont demandé : « Pouvez-vous prédire comment le métal bouge ? »

   • Résultat : Oui, très bien. L'ordinateur a prédit la vitesse et l'expansion du métal presque parfaitement. Cela leur a permis de savoir quelles « règles » mathématiques (équations d'état) décrivent le mieux le comportement du métal.

2. Le test de la « conductivité connue » : Ils ont injecté dans l'ordinateur la conductivité électrique réelle du métal (sa capacité à conduire) et ont demandé : « Pouvez-vous prédire l'électricité et le mouvement ? »

   • Résultat : Oui. L'ordinateur a prédit avec succès la tension et le courant, correspondant à l'expérience réelle. Cela a prouvé que les deux parties du modèle communiquent correctement entre elles.

3. Le test de la « prédiction pure » : C'était le plus difficile. Ils n'ont donné à l'ordinateur aucune donnée de l'expérience réelle. Ils lui ont simplement donné les lois de la physique et ont demandé : « Pouvez-vous prédire toute l'expérience à partir de zéro ? »

   • Résultat : C'était très proche. L'ordinateur a prédit la vitesse, le courant et la tension avec une bonne précision. Il y avait de petites différences (comme une erreur de 10 % sur la tension à la toute fin), mais l'image globale était correcte.

Pourquoi cela importe

L'article conclut que ce modèle informatique est un outil robuste et efficace.

Au lieu de simplement deviner comment configurer les futures expériences, les scientifiques peuvent désormais utiliser ce modèle de « vélo tandem » pour les concevoir. Ils peuvent simuler différents scénarios sur ordinateur pour voir ce qui se passera avant même d'allumer la machine réelle. Cela aide à comprendre la physique de la matière dense chaude sans avoir à compter uniquement sur des expériences coûteuses et difficiles.

En bref : Ils ont construit un jumeau numérique d'une explosion électrique à haute vitesse. Ils ont prouvé qu'il fonctionne en le comparant à de vraies explosions, et ils peuvent maintenant l'utiliser pour planifier de futures expériences en toute confiance.

Résumé technique

Résumé Technique : Simulations Hydrodynamiques de Feuilles Métalliques en Régime de Chauffage par Puissance Pulsée en État de Dilatation

Énoncé du Problème
La Matière Dense Chaude (MDC/WDM), le régime intermédiaire reliant la matière condensée et le plasma, est cruciale pour des domaines allant de la science planétaire à la fusion par confinement inertiel (FCI). Cependant, la modélisation de ce régime est complexe en raison de la coexistence des effets de corrélation ionique et de la dégénérescence électronique, que ni la physique des plasmas classiques ni les théories conventionnelles de la matière condensée ne décrivent avec précision. Bien que diverses techniques expérimentales existent pour produire de la MDC, les données concernant le régime de dilatation (densité inférieure à la densité nominale) sont particulièrement précieuses pour valider les équations d'état (EOS) et les modèles de conductivité pertinents pour les configurations de hohlraums de la FCI.

Le défi spécifique abordé dans ce travail est la conception et la simulation d'une expérience de puissance pulsée où de fines feuilles métalliques, confinées dans une cellule de saphir, sont chauffées par effet Joule pour atteindre le régime de la MDC dilatée. La conception de telles expériences est difficile car elle nécessite de prédire simultanément la réponse électrique du pilote de puissance pulsée et l'évolution hydrodynamique du matériau chauffé. De plus, une interprétation précise de la transition du solide vers un plasma partiellement ionisé nécessite un cadre numérique entièrement couplé qui lie le dépôt d'énergie électrique à l'évolution thermodynamique.

Méthodologie
Les auteurs ont développé un cadre de modélisation qui couple une description électrique d'un système de puissance pulsée avec un code hydrodynamique unidimensionnel (1D), ESTHER. L'approche repose sur une géométrie planaire 1D, ce qui simplifie le problème en supposant que la dilatation se produit dans une seule direction préférentielle en raison du confinement par haute impédance acoustique (saphir), permettant ainsi de négliger les effets de champ magnétique et de réduire le coût de calcul tout en conservant la physique essentielle.

La méthodologie se déroule en trois étapes de modélisation distinctes :

1. Modélisation du Circuit Électrique : Le pilote de puissance pulsée, l'interrupteur à éclateur (spark-gap) et la charge sont modélisés comme un circuit RLC équivalent. L'évolution du courant est résolue numériquement à l'aide d'un schéma de différence finie centrée de second ordre pour gérer la résistance et l'inductance dépendantes du temps.

   • Pilote : Caractérisé par des paramètres intrinsèques fixes ($R_f, L_f, C$).
   • Interrupteur : Modélisé à l'aide d'un modèle d'arc de Braginskii modifié, où la résistance et l'inductance de l'éclateur évoluent en fonction de l'intégrale du courant et du nombre de canaux parallèles ($N$).
   • Validation : Ce modèle électrique a été validé par des mesures en court-circuit de quatre générateurs différents (EPP1, EPP2, EOLE, GEPI2) couvrant des courants allant d'environ 140 kA à environ 5 MA.

2. Simulation Hydrodynamique : Le code ESTHER, un solveur hydrodynamique lagrangien 1D, a été utilisé pour simuler l'évolution thermodynamique de la feuille. Le code résout les équations de conservation de la masse, de la quantité de mouvement et de l'énergie, en incorporant une EOS et une conduction thermique. Le terme source de chauffage Joule est directement intégré dans l'équation de l'énergie.

3. Stratégies de Couplage : Trois approches ont été testées pour lier les domaines électrique et hydrodynamique :

   • Approche A (Puissance Expérimentale) : Le terme source est la puissance mesurée expérimentalement ($P = I_{exp} \times U_{exp}$). Cela permet d'isoler la performance de différents modèles d'EOS (SESAME, BLF, H´ebert et al.).
   • Approche B (Conductivité Expérimentale) : Le terme source est dérivé de la conductivité électrique expérimentale ($\sigma_{exp}$) combinée au solveur de circuit. Cela permet de valider le couplage du solveur de courant avec l'hydrodynamique.
   • Approche C (Entièrement Numérique) : Une approche auto-cohérente où la conductivité électrique est calculée à partir d'un modèle tabulé (Lee-More/Desjarlais modifié) basé sur l'état thermodynamique local (densité et température) de chaque cellule hydrodynamique. Cela élimine la dépendance vis-à-vis des signaux d'entrée expérimentaux.

Résultats Clés

• Validation Électrique : Le modèle de circuit proposé, incluant la dynamique de l'éclateur, a montré un excellent accord avec les mesures de courant en court-circuit pour une large gamme de générateurs et de conditions de fonctionnement (50 kA à 5 MA).
• Sélection de l'EOS : En utilisant la puissance expérimentale comme source, les modèles d'EOS BLF et H´ebert et al. ont reproduit nettement mieux les vitesses d'expansion, les densités et les énergies internes que le modèle SESAME, particulièrement concernant la transition solide-liquide et le plateau de vitesse subséquent.
• Validation du Couplage : En utilisant la conductivité expérimentale comme terme source, la simulation couplée a reproduit avec précision les signaux de courant et de tension expérimentaux. Les modèles BLF et H´ebert et al. ont fourni le meilleur accord pour l'évolution de la tension, bien que de légères divergences dans les calculs d'épaisseur de la feuille aient été notées.
• Simulation Entièrement Numérique : L'approche entièrement auto-cohérente (Approche C) a réussi à reproduire les tendances générales de la vitesse d'expansion, du courant et de la tension. Bien que de légers décalages temporels aient été observés (attribués à l'interaction entre les températures prédites par l'EOS et le modèle de conductivité), les simulations ont efficacement encadré le comportement expérimental. L'étude a souligné que les divergences dans l'approche entièrement numérique proviennent probablement de l'influence combinée de l'EOS (affectant les prédictions de température) et du modèle de conductivité, plutôt que d'une déficience unique.

Signification et Revendications
L'article affirme que cette approche numérique constitue une « méthode robuste et efficace » pour concevoir et optimiser les futures expériences de Matière Dense Chaude utilisant des installations de puissance pulsée. En couplant avec succès la réponse électrique du pilote avec l'évolution hydrodynamique de la charge, le cadre permet :

1. Optimisation de la Conception : Prédire la trajectoire thermodynamique des cibles avant l'expérimentation.
2. Test de Modèles : Fournir une méthodologie directe pour tester différents modèles d'EOS et de conductivité électrique sous des conditions pertinentes.
3. Accès aux Grandeurs Incommensurables : Donner accès à des propriétés telles que la température et l'énergie interne, qui sont difficiles ou impossibles à mesurer directement dans ces expériences.

Les auteurs soulignent que, bien que les simulations entièrement numériques présentent de légères divergences par rapport aux données expérimentales, l'accord global valide le cadre comme un outil puissant pour interpréter les expériences existantes et guider les conceptions futures, sans nécessiter l'invention de nouvelles applications au-delà de la portée des expériences de feuilles par puissance pulsée décrites ici.