Flow Shear Effects in the Onset Physics of Resistive MHD Instabilities in Tokamaks. Final report

The progress in this research centers around the computational analysis of flow shear effects in the onset of a 3/2 mode driven by a 1/1 mode in DIII-D equilibria. The initial idea was to try and calculate, via nonlinear simulations with NIMROD, the effects of rotation shear on driven 3/2 and 2/1 seed island physics, in experimentally relevant DIIID equilibria. The simulations indicated that very small seed islands were directly driven, as shielding between the sawtooth and the surfaces is significant at the high Lundquist numbers of the experiment. Instead, long after the initial crash the difference in linear stability of the 3/2, which remained prevalent despite the flattening of the core profiles from the sawtooth, contributed to a difference in the eventual seed island evolution. Essentially the seed islands grew or decayed long after the sawtooth crash, and not directly from it. Effectively the dominant 1/1 mode was found to be dragging the coupled modes surrounding it at a high rate through the plasma at their surfaces. The 1/1 mode is locked to the local frame of the plasma in the core, where the flow rate is greatest. The resonant perturbations at the surrounding surfaces propagate in the 'high slip regime' in the language of Fitzpatrick. Peaked flux averaged jxb forces (see Figs. 1 and 2) agree with localized flow modifications at the surfaces in analogy with Ebrahimi, PRL 2007. We track the mode into nonlinear saturation and have found oscillatory states in the evolution. During a visit (11/09) to Tulsa by R.J. LaHaye (GA), it became clear that similar oscillatory states are observed in DIII-D for these types of discharges