The onset of abrupt magnetic reconnection events, observed in the nonlinear evolution of double tearing modes (DTM), is investigated via reduced resistive magnetohydrodynamic simulations. We have identified the critical threshold for the parameters characterizing the linear DTM stability leading to the bifurcation to the explosive dynamics. A new type of secondary instability is discovered that is excited once the magnetic islands on each rational surface reach a critical structure characterized here by the width and the angle rating their triangularization. This new instability is an island structure-driven nonlinear instability, identified as the trigger of the subsequent nonlinear dynamics which couples flow and flux perturbations. This instability only weakly depends on resistivity.

Turbulent transport near the critical gradient in toroidal plasmas is studied based on global Landau-fluid simulations and an extended predator-prey theoretical model of ion temperature gradient turbulence. A new type of intermittent transport associated with the emission and propagation of a geodesic acoustic mode (GAM) is found near the critical gradient regime, which is referred to as GAM intermittency. The intermittency is characterized by new time scales of trigger, damping, and recursion due to GAM damping. During the recursion of intermittent bursts, stationary zonal flow increases with a slow time scale due to the accumulation of undamped residues and eventually quenches the turbulence, suggesting that a nonlinear upshift of the critical gradient, i.e., Dimits shift, is established through such a dynamical process.

Shear flows induced by nonlinear evolution of double tearing modes are investigated in a resistive magnetohydrodynamic model with slab geometry. It is found that intensive and thin poloidal shear flow layers are generated in the magnetic island region driven by coupled reconnection process at both rational surfaces. The structure of the flow layers keeps evolving after the merging of magnetic separatrices and forms a few narrow vortices along the open field lines in the final stage of magnetic reconnection. The effects of the distance between both rational surfaces and the initial magnetic shear on the nonlinear evolution of the plasma flows are also taken into consideration and the relevant mechanism is discussed.

The radial profiles of perpendicular flows in the presence of the magnetic island were firstly measured in the HL-2A tokamak by hopping the work frequency of the Doppler backward scattering reflectometer system along with a two-dimensional electron cyclotron emission imaging diagnostic identifying the island locations. It has been observed that across the O-point cut the perpendicular flow is quite small at the center of the island and strongly enhanced around the boundary of the island, resulting in a large increase of the flow shear in the outer half island, while across the X-point cut the flow is almost flat in the whole island region. Meanwhile it was found that the density fluctuations are generally weakened inside the island. The results indicate that both the perpendicular flow and the density fluctuation level are modulated by the naturally rotating tearing mode near the island boundary. The cross-correlation between the perpendicular flows and the oscillating electron temperature further reveals that the modulation of the perpendicular flow occurs mainly inside and in the vicinity of the island.

In order to better understand the feature of turbulent transport in tokamak plasmas dominated by zonal flows and nonlinearly excited tertiary (Kelvin–Helmholtz-like) waves, the radial dependency of the turbulent structure in the back of steady-state zonal flows was analyzed based on gyrofluid simulation of the electron temperature gradient turbulence in a sheared-slab configuration. The zonal flows change the characteristics of plasma turbulence from homogeneous to inhomogeneous. Thus, it is found that turbulent vortices are regulated by a radially alternative potential pattern where the spatiotemporal characteristic scales with respect to local fluctuations change radially. From the cross spectrum analysis between the poloidal electric field and the pressure perturbation, it was found that, in plasmas dominated by zonal flows, the restriction of heat flux is conducted by two mechanisms, i.e. the reduction of coherence and the phase synchronization.