Modeling of plasma transport during edge-localized mode in tokamak using a kinetic Vlasov-Poisson code

All'interno delle massicce macchine a forma di ciambella note come tokamak, gli scienziati stanno cercando di ricreare la potenza delle stelle per generare energia pulita. Questi dispositivi mantengono un gas surriscaldato, chiamato plasma, in posizione utilizzando potenti campi magnetici. Tuttavia, questo plasma non è perfettamente calmo; occasionalmente subisce improvvise e violente eruzioni note come modi localizzati al bordo (edge-localized modes). Immaginate il plasma come una pentola che bolle dove la superficie improvvisamente esplode, inviando un'ondata di calore intenso e particelle che si abbatte sulle pareti interne della macchina. Queste pareti sono rivestite di materiali speciali, spesso tungsteno, progettati per resistere al calore, ma esplosioni ripetute possono fonderli o staccarne dei frammenti. Se questi pezzi erosi tornano a fluttuare nel nucleo del plasma, possono avvelenare la reazione e interromperla. Comprendere esattamente come queste eruzioni viaggino dal centro della macchina alle pareti è fondamentale per costruire un reattore a fusione sicuro e duraturo, ma la fisica coinvolta è incredibilmente complessa, poiché comporta particelle che si muovono a velocità diverse e interagiscono in modi che i semplici modelli fluidodinamici non possono catturare appieno.

Per affrontare questa sfida, i ricercatori Ce Wang, Sven Van Loo e Geert Verdoolaege dell'Università di Ghent hanno sviluppato un nuovo strumento di simulazione al computer chiamato KOBRA. Invece di trattare il plasma come un fluido liscio, un approccio comune ma talvolta impreciso per questi eventi ad alta velocità, il loro codice traccia il comportamento di singole particelle mentre si muovono nello spazio e nel tempo. Questo metodo, noto come modellazione cinetica, permette di vedere dettagli fini, come il modo in cui gli elettroni veloci si comportano diversamente dagli ioni più pesanti. Il problema specifico che hanno affrontato era una questione di scala. La distanza che il plasma percorre dal centro della macchina alle pareti è enorme, circa dieci o venti metri, mentre la minuscola distanza su cui agiscono le forze elettriche è solo larga quanto un capello umano. Simulare ogni singolo dettaglio attraverso questa vasta distanza con una griglia uniforme standard richiederebbe un computer capace di eseguire un numero di calcoli così elevato da risultare proibitivo e lento.

Il team ha risolto questo problema utilizzando una tecnica chiamata raffinamento adattivo della mesh (adaptive mesh refinement). Pensate a questo come a una fotocamera digitale che zooma automaticamente per scattare una foto ad alta risoluzione solo dove sta accadendo l'azione, mantenendo lo sfondo a una risoluzione inferiore per risparmiare spazio. Nella loro simulazione, la griglia del computer diventa molto fine e dettagliata solo nelle regioni in cui la distribuzione delle particelle cambia rapidamente, come nel bordo d'attacco dell'esplosione del plasma, e rimane più grossolana nelle regioni più calme e regolari. Questo approccio ha permesso loro di eseguire una simulazione completa del trasporto del plasma dal piano medio della macchina ai target del divertore, ovvero le piastre specifiche progettate per assorbire il calore. Confrontando il loro metodo adattivo con una simulazione tradizionale che utilizzava una griglia ad alta risoluzione ovunque, hanno scoperto che il nuovo metodo produceva risultati quasi identici riguardo al flusso di particelle e calore, con differenze inferiori all'uno percento.

Le simulazioni hanno rivelato la specifica meccanica di come si sviluppano queste eruzioni. Nei primissimi istanti, un piccolo gruppo di elettroni ad alta energia corre davanti al resto del plasma, raggiungendo la parete prima che gli ioni più pesanti possano raggiungerli. Questo arrivo anticipato crea una separazione temporanea di carica, che genera un forte campo elettrico. Questo campo agisce quindi come un freno per gli elettroni veloci e una spinta per gli ioni, rimodellando rapidamente l'intera nube di plasma. I ricercatori hanno osservato che il flusso di particelle alla parete aumenta bruscamente e poi decade lentamente, un modello che corrisponde alle precedenti osservazioni sperimentali. Fondamentalmente, il metodo della mesh adattiva ha catturato questo picco iniziale e la successiva evoluzione con la stessa precisione del metodo a griglia uniforme molto più pesante, dimostrando che la strategia dello "zoom" non sacrifica la verità fisica in favore della velocità.

Oltre all'accuratezza, lo studio ha evidenziato i drammatici guadagni di efficienza di questo nuovo approccio. Sebbene la simulazione adattiva fosse inizialmente più lenta perché doveva riorganizzare costantemente la sua griglia, è diventata molto più veloce man mano che il plasma si stabilizzava e la necessità di un dettaglio estremo diminuiva. Alla fine della simulazione, il metodo adattivo ha completato il compito in circa 68 ore, mentre il metodo a griglia uniforme ha richiesto 113 ore, tagliando efficacemente il tempo di calcolo di quasi la metà. Inoltre, il metodo adattivo ha richiesto molta meno memoria informatica, utilizzando solo circa il 18 percento della memoria necessaria dal modello uniforme alla fine dell'esecuzione. Questa efficienza è vitale perché apre la porta a simulazioni più complesse e ad alta dimensionalità che erano precedentemente impossibili da eseguire. Il lavoro dimostra che, regolando intelligentemente la risoluzione della simulazione per adattarla al comportamento fisico del plasma, gli scienziati possono modellare queste pericolose eruzioni con maggiore velocità e minore consumo di risorse, fornendo una via più chiara per comprendere e mitigare i danni causati da questi eventi nei futuri reattori a fusione.