Gyrokinetic simulation of magnetic-island-induced electric potential vortex mode

Ion temperature gradient(ITG)-driven turbulence with embedded static magnetic islands is simulated by utilizing a gyrokinetic theory-based global turbulence transport code(GKNET)in this work.Different from the traditional equilibrium circular magnetic-surface average(EMSA)method,an advanced algorithm that calculates the perturbed magnetic-surface average(PMSA)of the electric potential has been developed to precisely deal with the zonal flow component in a non-circular magnetic surface perturbed by magnetic islands.Simulations show that the electric potential vortex structure inside islands induced by the magnetic islands is usually of odd parity when using the EMSA method.It is found that the odd symmetry vortex can transfer into an even one after a steep zonal flow gradient,i.e.the flow shear has been built in the vicinity of the magnetic islands by adopting the PMSA algorithm.The phase of the potential vortex in the poloidal cross section is coupled with the zonal flow shear.Such an electric potential vortex mode may be of essential importance in wide topics,such as the turbulence spreading across magnetic islands,neoclassical tearing mode physics,and also the interaction dynamics between the micro-turbulence and MHD activities.

A TVD-WAF-based hybrid finite volume and finite difference scheme for nonlinearly dispersive wave equations

A total variation diminishing-weighted average flux （TVD-WAF）-based hybrid numerical scheme for the enhanced version of nonlinearly dispersive Boussinesq-type equations was developed. The one-dimensional governing equations were rewritten in the conservative form and then discretized on a uniform grid. The finite volume method was used to discretize the flux term while the remaining terms were approximated with the finite difference method. The second-order TVD-WAF method was employed in conjunction with the Harten-Lax-van Leer （HLL） Riemann solver to calculate the numerical flux, and the variables at the cell interface for the local Riemann problem were reconstructed via the fourth- order monotone upstream-centered scheme for conservation laws （MUSCL）. The time marching scheme based on the third-order TVD Runge- Kutta method was used to obtain numerical solutions. The model was validated through a series of numerical tests, in which wave breaking and a moving shoreline were treated. The good agreement between the computed results, documented analytical solutions, and experimental data demonstrates the correct discretization of the governing equations and high accuracy of the proposed scheme, and also conforms the advantages of the proposed shock-capturing scheme for the enhanced version of the Boussinesq model, including the convenience in the treatment of wave breaking and moving shorelines and without the need for a numerical filter.

THE INTERACTION BETWEEN EDDY FLUX AND ZONAL AVERAGED FLOW IN“JUNE SUDDEN CHANGE”

In this paper,Eliassen-Palm theory and associated diagnostic method are used to discuss the dynamical mechanism of the interaction between eddy flux and zonal averaged flow in the sudden change of the general circulation in the Northern Hemisphere in early summer of 1982,which indicates that the northward jump of subtropic jet stream is closely related to the enhanced transport of mid-latitude eddy energy to the subtropic tro- popause.However,the direction of the transport is conditioned by the structure of zonal averaged flow.It is evident that the adjustment of zonal averaged flow responses quickly to the eddy transport during this episode. As regard to the suddenness of the change of circulation,the critical role is played by the eddy forcing,