Benchmark of gyrokinetic, kinetic MHD and gyrofluid codes for the linear calculation of fast particle driven TAE dynamics

Könies, A.; Briguglio, S.; Gorelenkov, N. N.; Fehér, Tamás; Isaev, M. Yu.; Lauber, P.; Mishchenko, A.; Spong, D. A.; Todo, Y.; Cooper, W. A.; Hatzky, R.; Kleiber, R.; Borchardt, M.; Vlad, G.; Biancalani, A.; Bottino, A.; Tg, Itpa Ep

doi:10.1088/1741-4326/aae4e6

Cited by 62 publications

(81 citation statements)

References 24 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The solver modifications needed for the pullback scheme implementation have been described. The new scheme has been verified using the ITPA-TAE benchmark [19,20] both in the linear and nonlinear regimes. A considerable improvement of the code efficiency has been observed.…”

Section: Discussionmentioning

confidence: 99%

“…Whereas the control variate mitigation is capable reproducing only the first three points, becoming numerically unstable for larger time steps, the pullback works robustly even for time steps two orders of magnitude larger than the time steps typically used for the electromagnetic simulations with the control variate only. In this simulation, the fast-ion density has been ten times larger (corresponding to 3% particle content) comparing to the standard "ITPA-TAE" benchmark [19,20]. This makes the simulations faster and more robust.…”

Section: Toroidal Alfvén Eigenmodementioning

confidence: 93%

See 1 more Smart Citation

Pullback scheme implementation in ORB5

Mishchenko

Bottino

Biancalani

et al. 2019

Computer Physics Communications

Self Cite

View full text Add to dashboard Cite

The pullback scheme is implemented in the global gyrokinetic particle-incell code ORB5 [S. Jolliet et al, Comp. Phys. Comm., 177, 409 (2007)] to mitigate the cancellation problem in electromagnetic simulations. The equations and the discretisation used by the code are described. Numerical simulations of the Toroidal Alfvén Eigenmodes are performed in linear and nonlinear regimes to verify the scheme. A considerable improvement in the code efficiency is observed. For the internal kink mode, it is shown that the pullback mitigation efficiently cures a numerical instability which would make the simulation more costly otherwise.

show abstract

Section: Discussionmentioning

confidence: 99%

Section: Toroidal Alfvén Eigenmodementioning

confidence: 93%

Pullback scheme implementation in ORB5

Mishchenko

Bottino

Biancalani

et al. 2019

Computer Physics Communications

Self Cite

View full text Add to dashboard Cite

show abstract

“…The equations of motion for each computational particle are solved using a fourth-order Runge-Kutta method. MEGA code participated in the code benchmark of the Energetic Particle Physics Topical Group of the International Tokamak Physics Activity [43]. Good agreements were found in the spatial profile, frequency, and growth rate of a TAE among the seven codes without the energetic ion FLR effects and among the six codes with the FLR effects.…”

Section: Simulation Model 21 Hybrid Simulation Model For Energetic mentioning

confidence: 97%

Multi-phase hybrid simulation of energetic particle driven magnetohydrodynamic instabilities in tokamak plasmas

Todo¹

2016

New J. Phys.

View full text Add to dashboard Cite

Magnetohydrodynamic (MHD) instabilities driven by energetic particles in tokamak plasmas and the energetic particle distribution formed with the instabilities, neutral beam injection, and collisions are investigated with hybrid simulations for energetic particles and an MHD fluid. The multi-phase simulation, which is a combination of classical simulation and hybrid simulation, is applied to examine the distribution formation process in the collisional slowing-down time scale of energetic ions for various beam deposition power (P NBI ) and slowing-down time (t s ). The physical parameters other than P NBI and t s are similar to those of a Tokamak Fusion Test Reactor (TFTR) experiment (Wong et al 1991 Phys. Rev. Lett. 66 1874). For P NBI = 10 MW and t s = 100 ms, which is similar to the TFTR experiment, the bursts of toroidal Alfvén eigenmodes take place with a time interval 2 ms, which is close to that observed in the experiment. The maximum radial velocity amplitude (v r ) of the dominant TAE at the bursts in the simulation is~´v v 3 10 r A 3 where v A is the Alfvén velocity at the plasma center. For P NBI = 5 MW and t s = 20 ms, the amplitude of the dominant TAE is kept at a constant level~´v v 4 10 r A 4 . The intermittency of TAE rises with increasing P NBI and increasing t s (= decreasing collision frequency). With increasing volume-averaged classical energetic ion pressure, which is well proportional to t P NBI s , the energetic ion confinement degrades monotonically due to the transport by the instabilities. The volume-averaged energetic ion pressure depends only on the volume-averaged classical energetic ion pressure, not independently on P NBI or t s . The energetic ion pressure profile resiliency, where the increase in energetic ion pressure profile is saturated, is found for the cases with the highest t P NBI s where the TAE bursts take place.

show abstract

“…When considered, the fast particles have a density profile peaked on axis (see Fig.3 on the left). The magnetic equilibrium and profiles are those of the ITPA-TAE international benchmark case [6] and the safety factor profile is shown in Fig.3 on the right. The main results that will be displayed in this section, have been obtained considering heavier electrons: m e = m H /200.…”

Section: Landau Dampingmentioning

confidence: 99%

“…In the simulations in Section 4, a small but finite inverse aspect ratio will be considered, using the equilibrium profiles of the International Tokamak Physics Activity (ITPA, see Ref. [6]). In Section 5 the studies on the linear and nonlinear growth rate and frequency spectra conducted considering experimental profiles from the NLED-AUG case [7] will be presented.…”

Section: Introductionmentioning

confidence: 99%

Gyrokinetic investigation of the damping channels of Alfvén modes in ASDEX Upgrade

et al. 2020

Self Cite

View full text Add to dashboard Cite

The linear destabilization and nonlinear saturation of energetic-particle driven Alfvénic instabilities in tokamaks strongly depend on the damping channels. In this work, the collisionless damping mechanisms of Alfvénic modes are investigated within a gyrokinetic framework, by means of global simulations with the particle-in-cell code ORB5, and compared with the eigenvalue code LIGKA and reduced models. In particular, the continuum damping and the Landau damping (of ions and electrons) are considered. The electron Landau damping is found to be dominant on the ion Landau damping for experimentally relevant cases. As an application, the linear and nonlinear dynamics of toroidicity induced Alfvén eigenmodes and energetic-particle driven modes in ASDEX Upgrade is investigated theoretically and compared with experimental measurements.

show abstract

Benchmark of gyrokinetic, kinetic MHD and gyrofluid codes for the linear calculation of fast particle driven TAE dynamics

Cited by 62 publications

References 24 publications

Pullback scheme implementation in ORB5

Pullback scheme implementation in ORB5

Multi-phase hybrid simulation of energetic particle driven magnetohydrodynamic instabilities in tokamak plasmas

Gyrokinetic investigation of the damping channels of Alfvén modes in ASDEX Upgrade

Contact Info

Product

Resources

About