Since 18 of December 2019 uses Nucleus credentials. Visit our help pages for information on how to Register and Sign-in using Nucleus.
22-27 October 2018
Mahatma Mandir Conference Centre
Asia/Kolkata timezone

Nonlinear 3D simulations of Vertical Displacement Events in tokamaks

26 Oct 2018, 14:00
4h 45m
Mahatma Mandir Conference Centre

Mahatma Mandir Conference Centre

Gandhinagar (nearest Airport: Ahmedabad), India
Poster P8 Posters


Dr Nathaniel M. Ferraro (UsPPPL)


Vertical displacement events (VDEs) where the plasma moves rapidly towards the wall can cause large electromagnetic forces on the vessel structures with possible damaging effects for large tokamaks. Non-axisymmetric modes developing during the VDE can lead to asymmetric, sometimes rotating forces on the vessel which can be even more severe. Large-scale 3D simulations play a crucial role on the path towards assessing and preventing the damaging effects of VDEs on vessel components in future large tokamaks like ITER.

We use the high-order finite element code M3D-C^1 [1] to perform 2D and 3D nonlinear MHD simulations of VDEs in tokamaks including a resistive wall model [2]. In order to develop predictive capabilities, the simulation results are benchmarked with other codes as well as validated against existing experimental measurements.

2D and 3D nonlinear MHD simulations of VDEs are based on and validated against discharges in NSTX [3] as well as DIII-D. The results of a set of axisymmetric VDE calculations based on NSTX discharge #132859 show the sensitivities of the early VDE evolution to different parameters, in particular the halo resistivity. 3D simulations show how non-axisymmetric modes arise in the late VDE phase and lead to a stochastization of the magnetic field lines which allows for an efficient release of thermal energy into the wall. The thermal quench is followed by a fast decay of the plasma current and rise of the wall current.

A detailed benchmarking activity between the M3D-C^1 code and the 3D nonlinear MHD code NIMROD [4] based on an NSTX discharge is being performed. The comparison of axisymmetric VDE calculations is focused on the early VDE growth and the wall forces. We plan to extent this benchmark to 3D simulations with a 2D wall. In addition, an axisymmetric benchmark between the M3D-C^1 code and the CarMaONL code [5] based on a standard ITER scenario using a simplified 2D model of the ITER first wall is in progress.

[1] S. C. Jardin, et al., Comput. Sci. Discov. 5, 014002 (2012).

[2] N. M. Ferraro, et al., Phys. Plasmas 23, 056114 (2016).

[3] D. Pfefferle, et al., Phys. Plasmas (2018). Submittted.

[4] C. R. Sovinec, et al., J. Comput. Phys. 195, 355 (2004).

[5] F. Villone, et al, Plasma Phys. Controlled Fusion 55, 095008 (2013).

Paper Number TH/P8-10
Country or International Organization United States of America

Primary author

Dr Isabel Krebs (UsPPPL)


Prof. Carl Sovinec (UsUWisMad) Dr David Pfefferle (The University of Western Australia, Perth, WA 6009, Australia) Dr Ellen Seegyoung Seol (Rensselaer Polytechnic Institute, Troy, NY 12180, USA) Prof. Fabio Villone (CREATE, DIETI, Universita degli Studi di Napoli Federico II, Napoli, Italy) Prof. Guglielmo Rubinacci (CREATE, DIETI, Universita degli Studi di Napoli Federico II, Napoli, Italy) Dr Henry Strauss (HRS Fusion, West Orange, NJ 07052, USA) Dr Jin Chen (UsPPPL) Dr Kyle Bunkers (UsUWisMad) Dr Lang Lao (UsGA) Prof. Mark S. Shephard (Rensselaer Polytechnic Institute, Troy, NY 12180, USA) Dr Nathaniel M. Ferraro (UsPPPL) Dr Nicholas Eidietis (UsGA) Prof. Stephen C. Jardin (UsPPPL)

Presentation Materials