In this paper we present a coupled Finite Element Method - Boundary Element Method (FEM-BEM) approach for the solution of the free-boundary axi-symmetric plasma equilibrium problem. The proposed method, obtained from an improvement of the Hagenow-Lackner coupling method, allows to efficiently model the equilibrium problem in unbounded domains by discretizing only the plasma region; the external conductors can be modelled either as 2D or 3Dmodels, according to the problem of interest. The paper explores different iterative methods for the solution of the nonlinear Grad-Shafranov equation, such as Picard, Newton-Raphson and Newton-Krylov, in order to provide a robust and reliable tool, able to handle large-scale problems (e.g. high resolution equilibria). This method has been implemented in the FRIDA code (FReeboundary Integro-Differential Axisimmetric - https://github.om/matteobonotto/ FRIDA), together with a suitable Adaptive Integration Technique (AIT) for the computation of the source term. FRIDA has been successfully tested and validated against experimental data from RFX-mod device, and numerical equilibria of an ITER-like device.

A Coupled FEM-BEM Approach for the Solution of the Free-Boundary Axi-Symmetric Plasma Equilibrium Problem

Bettini P.;
2021

Abstract

In this paper we present a coupled Finite Element Method - Boundary Element Method (FEM-BEM) approach for the solution of the free-boundary axi-symmetric plasma equilibrium problem. The proposed method, obtained from an improvement of the Hagenow-Lackner coupling method, allows to efficiently model the equilibrium problem in unbounded domains by discretizing only the plasma region; the external conductors can be modelled either as 2D or 3Dmodels, according to the problem of interest. The paper explores different iterative methods for the solution of the nonlinear Grad-Shafranov equation, such as Picard, Newton-Raphson and Newton-Krylov, in order to provide a robust and reliable tool, able to handle large-scale problems (e.g. high resolution equilibria). This method has been implemented in the FRIDA code (FReeboundary Integro-Differential Axisimmetric - https://github.om/matteobonotto/ FRIDA), together with a suitable Adaptive Integration Technique (AIT) for the computation of the source term. FRIDA has been successfully tested and validated against experimental data from RFX-mod device, and numerical equilibria of an ITER-like device.
File in questo prodotto:
Non ci sono file associati a questo prodotto.
Pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11577/3445679
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 6
  • ???jsp.display-item.citation.isi??? 6
  • OpenAlex ND
social impact