Theory and simulation of quasilinear transport from external magnetic field perturbations in a DIII-D plasma

The linear response profiles for the 3D perturbed magnetic fields, currents, ion velocities, plasma density, pressures, and electric potential from low-n external resonant magnetic field perturbations (RMPs) are obtained from the collisional two-fluid M3D-C¹ code [N. M. Ferraro and S. C. Jardin, J. Comput. Phys. 228, 7742 (2009)]. A newly developed post-processing RMPtran code computes the resulting quasilinear E×B and magnetic (J×B) radial transport flows with respect to the unperturbed flux surfaces in all channels. RMPtran simulations focus on ion (center of mass) particle and transient non-ambipolar current flows, as well as the toroidal angular momentum flow. The paper attempts to delineate the RMP transport mechanisms that might be responsible for the RMP density pump-out seen in DIII-D [M. A. Mahdavi and J. L. Luxon, Fusion Sci. Technol. 48, 2 (2005)]. Experimentally, the starting high toroidal rotation does not brake to a significantly lower rotation after the pump-out suggesting that convective and E×B transport mechanisms dominate. The direct J×B torque from the transient non-ambipolar radial current expected to accelerate plasma rotation is shown to cancel much of the Maxwell stress J×B torque expected to brake the plasma rotation. The dominant E×B Reynolds stress accelerates rotation at the top of the pedestal while braking rotation further down the pedestal.