Author: Nathaniel M Ferraro
Requested Type: Poster
Submitted: 2018-02-28 14:40:16
Co-authors: S.C. Jardin, I. Krebs, L.L. Lao, B.C. Lyons, H.R. Strauss
Princeton Plasma Physics Laboratory
100 Stellarator Road
Princeton, NJ 08540
We demonstrate new capabilities for integrated modeling of tokamak disruptions with M3D-C1. Among these new capabilities is the ability to model impurity species. Ionization, recombination, and radiation from impurities are evolved using the KPRAD model, and impurity charge state densities are advected with the ExB velocity. Separate equations for electron and ion temperature are included. A simple model for runaway electron generation is also now implemented. M3D-C1 includes a fully compressible extended-MHD model for the plasma and open field-line region, and includes a model for a resistive wall of arbitrary thickness that allows spatially resolved halo currents and eddy currents. The halo width is determined self-consistently by the anisotropy of the thermal diffusivity in the temperature evolution equations. Together, these capabilities enable M3D-C1 to model the precursor, thermal quench, and current quench phases of the disruption self-consistently and in great detail. M3D-C1 is presently being used to model disruptions in NSTX, NSTX-U, DIII-D, and JET. These integrated simulations are being validated using synthetic flux loop and magnetics diagnostics. Ultimately, we seek to provide a validated predictive model of wall forces and thermal loads due to disruptions in tokamaks. This work is supported by the US Department of Energy through the Center for Tokamak Transients Simulation (CTTS) SciDAC program and contract DE-AC02-09-CH11466.