<p>After a brief introduction on present day simulation and modeling of burning plasmas in magnetically confined devices, the results of a non-linear benchmark are presented, undertaken among several state-of-the-art codes available to study the self-consistent interaction of an Energetic Particle (EP) population with shear Alfvén waves: HYMAGYC, MEGA, ORB5 and XTOR-K. The first two codes, HYMAGYC and MEGA, are hybrid codes: HYMAGYC is a MHD-Gyrokinetic code (the bulk plasma is represented by MHD equations, while the EP species is treated using the gyrokinetic formalism), while MEGA is an MHD-Drift-Kinetic code (the bulk plasma is represented by MHD equations, while the EP species is treated using the drift-kinetic formalism, with the possibility of an ad-hoc gyroaveraging); ORB5 is a global electromagnetic gyrokinetic code (both bulk and EP species are treated using the gyrokinetic formalism); XTOR-K is a non-linear kinetic-MHD code (the bulk plasma is described by a set of non-linear resistive two-fluid MHD equations, extended to include kinetic effects of multiple ion species with a fully kinetic PIC module). The equilibrium of the so-called NLED-AUG reference case, in its version with peaked off-axis EP density profile, has been used, while considering a <InlineEquation ID="IEq1"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="41614_2025_199_Article_IEq1.gif" Format="GIF" Height="19" Rendition="HTML" Resolution="72" Type="Linedraw" Width="52" /> </InlineMediaObject> <EquationSource Format="TEX">$|n|=1$</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mo stretchy="false">|</mo> <mi>n</mi> <mo stretchy="false">|</mo> <mo>=</mo> <mn>1</mn> </mrow> </math></EquationSource> </InlineEquation> perturbation. This non-linear benchmark is the natural continuation of the linear benchmark already considered in the recent past, and represent a first-ever code comparison in the deep non-linear stage. In the present study fluid non-linearities are omitted, and the focus will be on comparing wave-particle interactions effects across codes. Characteristics of the non-linear saturation of the mode, self-consistent modification to the EP density profile and other features are compared among the considered codes. This brief review presents the state of the art of gyrokinetic and hybrid MHD-kinetic codes emerging as tools for studying reactor-relevant burning plasmas in realistic conditions.</p>

State of the art of gyrokinetic and hybrid MHD-kinetic codes through non-linear benchmarking to study reactor relevant burning plasmas

期刊论文
出版模式：

开放获取
发表日期：

2025年9月9日

查看全文

Gregorio Vlad,
Xin Wang,
Francesco Vannini,
Hinrich Lütjens,
Sergio Briguglio,
Nakia Carlevaro,
Matteo V. Falessi,
Giuliana Fogaccia,
Valeria Fusco,
Fulvio Zonca,
Alessandro Biancalani,
Alberto Bottino,
Thomas Hayward-Schneider,
Philipp Lauber,
Baruch Rofman,
Laurent Villard

摘要

After a brief introduction on present day simulation and modeling of burning plasmas in magnetically confined devices, the results of a non-linear benchmark are presented, undertaken among several state-of-the-art codes available to study the self-consistent interaction of an Energetic Particle (EP) population with shear Alfvén waves: HYMAGYC, MEGA, ORB5 and XTOR-K. The first two codes, HYMAGYC and MEGA, are hybrid codes: HYMAGYC is a MHD-Gyrokinetic code (the bulk plasma is represented by MHD equations, while the EP species is treated using the gyrokinetic formalism), while MEGA is an MHD-Drift-Kinetic code (the bulk plasma is represented by MHD equations, while the EP species is treated using the drift-kinetic formalism, with the possibility of an ad-hoc gyroaveraging); ORB5 is a global electromagnetic gyrokinetic code (both bulk and EP species are treated using the gyrokinetic formalism); XTOR-K is a non-linear kinetic-MHD code (the bulk plasma is described by a set of non-linear resistive two-fluid MHD equations, extended to include kinetic effects of multiple ion species with a fully kinetic PIC module). The equilibrium of the so-called NLED-AUG reference case, in its version with peaked off-axis EP density profile, has been used, while considering a $|n|=1$ $| n | = 1$ perturbation. This non-linear benchmark is the natural continuation of the linear benchmark already considered in the recent past, and represent a first-ever code comparison in the deep non-linear stage. In the present study fluid non-linearities are omitted, and the focus will be on comparing wave-particle interactions effects across codes. Characteristics of the non-linear saturation of the mode, self-consistent modification to the EP density profile and other features are compared among the considered codes. This brief review presents the state of the art of gyrokinetic and hybrid MHD-kinetic codes emerging as tools for studying reactor-relevant burning plasmas in realistic conditions.