Since 18 of December 2019 uses Nucleus credentials. Visit our help pages for information on how to Register and Sign-in using Nucleus.

May 10 – 15, 2021
Virtual Event
Europe/Vienna timezone
The Conference will be held virtually from 10-15 May 2021

Nonlinear MHD modeling of divertor striations in DIII-D RMP ELM suppressed discharges

May 11, 2021, 8:30 AM
Virtual Event

Virtual Event

Regular Poster Magnetic Fusion Experiments P1 Posters 1


Dmitri Orlov (University of California San Diego)


In this work we show that nonlinear MHD plasma response simulations are essential for understanding and predicting accurate heat and particle flux striations in DIII-D during ELM suppression. Understanding the nature of heat and particle distributions on the divertor plates due to splitting of the separatrix by 3D magnetic perturbation in RMP ELM suppressed discharges is an important issue for preventing the premature degradation of divertor components in ITER. Our simulations now match the experimental observations in DIII-D. We further show that in cases with no observable heat flux splitting during ITER-like shaped DIII-D discharges with RMP ELM suppression the loss of heat flux splitting may be due to an increase in carbon radiation above the surface. Since carbon will not be present in ITER, these results raise a question concerning whether a similar effect will be observed in ITER during divertor detachment operations with medium Z impurities such as neon or nitrogen.

(a) HFS divertor footprints in JOREK simulations without plasma response (top), with plasma response without neutrals (middle) and with neutrals (bottom). Dashed line indicates location of visual imaging camera. (b) Comparison of measured HFS divertor particle flux striation peaks (blue bars) with vacuum (magenta) and nonlinear plasma response with neutrals (black) JOREK MHD simulations in DIII-D discharge 166439 at 3300 ms. Field lines with connection lengths of 75 m and longer correspond to particle flux lobes on the 45-degree tile and the center post.

In DIII-D, strike point splitting is routinely observed in the divertor particle flux during operation with RMPs (1). We use small modulations of the in-vessel I-coil n=3 RMP current amplitudes or toroidal rotation of n=2 RMP fields in order to modify the poloidal spectrum of applied perturbation fields. These modulations affect the position and the size of the divertor footprints, and provide useful information to validate different numerical models for non-axisymmetric footprint lobes by comparing to high-precision visible imaging and IR tomography experimental measurement. The observed splitting is consistent with the toroidal mode number n of the applied perturbation but the measured separation of the divertor particle flux footprint lobes exceeds predictions of vacuum (TRIP3D) and linear plasma response (MARS) models by factors of 3-5 (2). The plasma response to the RMP in ITER-like conditions using linear, resistive MHD simulations (M3D-C1, NIMROD) with both single-fluid and two-fluid models is dominantly a screening response that reduces the divertor footprint splitting below the vacuum model predictions. The nonlinear MHD simulations using JOREK, a fluid model for the main plasma and the neutrals, show a much better match to the measured separation of the lobes (Figure 1). The time evolution of the neutrals in JOREK is described by a diffusion model combined with a boundary condition that reflects outgoing ions as incoming neutrals (3). In this work we also examine the possibility of near-SOL field lines affecting the formation of the outer-most lobes in the measured particle flux footprints.

(a) Evolution of the ISP heat flux profile from temperature rise on divertor tile as measured by IRTV in ISS discharge 147170 with n=3 RMP ELM suppression. (b) Evolution of ISP surface power calculated from IRTV measurements and radiated power in lower divertor region from bolometer measurements. Vertical lines correspond to times in panel (a). 2D tomographic reconstruction of CIII emission in the lower divertor in shot 147170 at 2104 ms (c) and 4509 ms (d) assuming toroidal symmetry.

At the same time, the heat flux to the divertor often does not show significant splitting in DIII-D ITER-like RMP ELM-suppressed discharge, as displayed in Figure 2, which is potentially good for ITER. One hypothesis suggested by observations is that the lack of splitting in heat flux may be related to the C-III volumetric radiation immediately above the surface that obscures the lobe structure. We performed fully 3D plasma-fluid and kinetic edge neutral transport Monte-Carlo EMC3-EIRENE simulations using vacuum and linear plasma response MHD solutions, and the results highlight the effect of different levels of carbon impurity radiation near the strike point on the divertor footprints. We also examined a possibility of heat flux lobe smearing partially due to the ion grad-B drifts as MAFOT simulations suggest. MAFOT results also indicate minor effects on heat flux footprints due to ExB fields and plasma sheath near the divertor surfaces.

(1) Moyer et al., RSI 89, 10E106 (2018)
(2) Orlov et al., IAEA FEC 2018 (Ahmedabad, India)
(3) Hujismans et al., 46th EPS, P1.1059 (2018)

This work is supported by the US Department of Energy under DE-FG02-07ER54917, DE-FG02-05ER54809, DE-FC02-04ER54698, DE-SC0012706, DE-AC52-07NA27344, DE-NA0003525, DE-AC02-09CH11466 and DE-AC04-94AL85000.

Affiliation UC San Diego
Country or International Organization United States

Primary authors

Dmitri Orlov (University of California San Diego) Todd E. Evans (General Atomics) Richard Moyer (University of California San Diego) Dr Igor Bykov (University of California San Diego) Max Fenstermacher (LLNL @ DIII-D) Charles Lasnier (Laurence Livermore National Laboratory) Dr Huiqian Wang (General Atomics, San Diego, CA 92186, USA) Stefano Munaretto (General Atomics) Dr Florian Effenberg (Princeton Plasma Physics Laboratory, Princeton, NJ, 08543 USA) Marina Becoulet (IRFM/CEA) Dr Eric Howell (Tech-X) Jacob King (Tech-X Corporation) Dr Andreas Wingen (ORNL) Dr Guido Hujismans (CEA, France)

Presentation materials