Magnetohydrodynamic properties of the intermediate and late nonlinear phases of peeling-ballooning modes: An analytic and computational approach

Analytic and numerical techniques are employed to describe the nonlinear evolution properties of edge-localized MHD instabilities. It is shown that the nonlinear dynamics of these instabilities depend upon the nature of the linear MHD instability drive. Edge-localized ideal MHD modes for two shifted-circle tokamak equilibria are examined nonlinearly using the extended-MHD code NIMROD. The equilibria have H-mode-like profiles with differing edge currents and pressure gradients to produce linear stability spectra that peak at low-n (peeling-like) and high-n (ballooning- like). Numerical simulations characterize the evolution from the linear to intermediate nonlinear regime, which is characterized by an ideal MHD displacement amplitude comparable to the mode width. For the peeling-like plasma, an edge-localized, n = 0, sheared cross-field plasma flow arises. As the instability enters the intermediate nonlinear regime the flow-shearing rate becomes comparable to the linear growth rate and nonlinearly slows the mode growth. This result is in contrast to the ballooning-like equilibria where a smaller shear rate is generated and the mode grows exponentially with the linear growth rate well into the intermediate nonlinear regime. The numerical results are consistent with an analytic theory derived to describe the nonlinear evolution properties of edge localized ideal MHD modes. A set of nonlinear evolution equations for the intermediate nonlinear phase of peeling mode instability is derived using the ratio of the radial mode width to the minor radius as an asymptotic expansion parameter. This differs from previous work on the nonlinear ballooning instability where a conventional ballooning order is used to derive the nonlinear evolution equations.;Lastly, the formation of filaments is observed in the late nonlinear regime, simulations using temperature dependent Spitzer resistivity profiles show that plasma with dissipation parameters typical of DIII-D discharges evolve ideally, and the addition of dissipation parameters typically required for numerical convergence lead to dramatically different plasma dynamics.;*The research was performed under appointment to the Fusion Energy Sciences Fellowship Program administered by Oak Ridge Institute for Science and Education under a contract between the U.S. Department of Energy and the Oak Ridge Associated Universities.