Cyclokinetic models and simulations for high-frequency turbulence in fusion plasmas

Gyrokinetics is widely applied in plasma physics. However, this framework is limited to weak turbulence levels and low drift-wave frequencies because high-frequency gyro-motion is reduced by the gyro-phase averaging. In order to test where gyrokinetics breaks down, Waltz and Zhao developed a new theory, called cyclokinetics [R. E. Waltz and Zhao Deng, Phys. Plasmas 20, 012507 (2013)]. Cyclokinetics dynamically follows the high-frequency ion gyro-motion which is nonlinearly coupled to the low-frequency drift-waves interrupting and suppressing gyro-averaging. Cyclokinetics is valid in the high-frequency (ion cyclotron frequency) regime or for high turbulence levels. The ratio of the cyclokinetic perturbed distribution function over equilibrium distribution function δf/F can approach 1.This work presents, for the first time, a numerical simulation of nonlinear cyclokinetic theory for ions, and describes the first attempt to completely solve the ion gyro-phase motion in a nonlinear turbulence system. Simulations are performed [Zhao Deng and R. E. Waltz, Phys. Plasmas 22(5), 056101 (2015)] in a local flux-tube geometry with the parallel motion and variation suppressed by using a newly developed code named rCYCLO, which is executed in parallel by using an implicit time-advanced Eulerian (or continuum) scheme [Zhao Deng and R. E. Waltz, Comp. Phys. Comm. 195, 23 (2015)]. A novel numerical treatment of the magnetic moment velocity space derivative operator guarantee saccurate conservation of incremental entropy.By comparing the more fundamental cyclokinetic simulations with the corresponding gyrokinetic simulations, the gyrokinetics breakdown condition is quantitatively tested. Gyrokinetic transport and turbulence level recover those of cyclokinetics at high relative ion cyclotron frequencies and low turbulence levels, as required. Cyclokinetic transport and turbulence level are found to be lower than those of gyrokinetics at high turbulence levels and low-Ω* values with stable ion cyclotron modes. The gyrokinetic approximation is found to break down when the density perturbation exceeds 20%, or when the ratio of nonlinear E×B frequency over ion cyclotron frequency exceeds 20%. This result indicates that the density perturbation of the Tokamak L-mode near-edge is not sufficiently large for breaking the gyro-phase averaging. For cyclokinetic simulations with sufficiently unstable ion cyclotron (IC) modes and sufficiently low Ω* ～10, the high-frequency component of the cyclokinetic transport can exceed that of the gyrokinetic transport. However, the low-frequency component of the cyclokinetic transport does not exceed that of the gyrokinetic transport. For higher and more physically relevant Ω* ?50 values and physically realistic IC driving rates, the low-frequency component of the cyclokinetic transport remains smaller than that of the gyrokinetic transport. In conclusion, the “L-mode near-edge short-fall” phenomenon, observed in some low-frequency gyrokinetic turbulence transport simulations, does not arise owing to the nonlinear coupling of high-frequency ion cyclotron motion to low-frequency drift motion.     

Critically balanced ion temperature gradient turbulence in fusion plasmas

Scaling laws for ion temperature gradient driven turbulence in magnetized toroidal plasmas are derived and compared with direct numerical simulations. Predicted dependences of turbulence fluctuation amplitudes, spatial scales, and resulting heat fluxes on temperature gradient and magnetic field line pitch are found to agree with numerical results in both the driving and inertial ranges. Evidence is provided to support the critical balance conjecture that parallel streaming and nonlinear perpendicular decorrelation times are comparable at all spatial scales, leading to a scaling relationship between parallel and perpendicular spatial scales. This indicates that even strongly magnetized plasma turbulence is intrinsically three?dimensional.     

Nonlinear saturation of trapped electron modes via perpendicular particle diffusion

In magnetized fusion plasmas, trapped electron mode (TEM) turbulence constitutes, together with ion temperature gradient (ITG) turbulence, the dominant source of anomalous transport on ion scales. While ITG modes are known to saturate via nonlinear zonal flow generation, this mechanism is shown to be of little importance for TEM turbulence in the parameter regime explored here. Instead, a careful analysis of the statistical properties of the ExB nonlinearity in the context of gyrokinetic turbulence simulations reveals that perpendicular particle diffusion is the dominant saturation mechanism. These findings allow for the construction of a rather realistic quasilinear model of TEM induced transport.     

Wave turbulent statistics in non-weak wave turbulence

In wave turbulence, which is made by nonlinear interactions among waves, it has been believed that statistical properties are well described by the weak turbulence theory, where separation of linear and nonlinear time scales derived from weak nonlinearity is assumed. However, the separation of the time scales is often violated. To get rid of this inconsistency, closed equations are derived in wave turbulence without assuming the weak nonlinearity according to Direct-Interaction Approximation (DIA), which has been successful in Navier-Stokes turbulence. The DIA equations is a natural extension of the conventional kinetic equation to not-necessarily-weak wave turbulence.     

Dynamics in a multicomponent plasma near the low-frequency cutoff

The distinctive feature of a multicomponent plasma is found to be a first order rotation of the light ion fluid with a characteristic frequency, Omega(r). A resonance at omega=Omega(r) results in plasma oscillation with unique properties. Coupling of the fast rotation time scale with the slow magnetosonic time scale leads to a nonlinear Schro dinger equation for the system and suggests the possibility of strong structural turbulence at MHD scales.     

Prediction of significant tokamak turbulence at electron gyroradius scales

The experimental conditions under which tokamak turbulence at hyperfine (electron gyroradius) scales is predicted to be significant and observable are described. The first quantitative predictions of fluctuation amplitudes, spectral features, and the associated electron energy transport are presented. A novel theoretical model which quantitatively describes the boundaries of the high-amplitude streamer transport regime is presented and shown to explain the gyrokinetic simulation results. This model uniquely includes consideration of two distinct secondary instabilities.     

Intermittent transport associated with the geodesic acoustic mode near the critical gradient regime

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.     

Electrostatic instability in magnetically confined inhomogeneous plasma driven by nonlinear force

A theoretical investigation on amplification of electrostatic ion acoustic wave in magnetically confined plasma has been presented in this paper. This investigation considers nonlinear wave–particle interaction process, called plasma maser effect, in presence of drift wave turbulence supported by magnetically confined inhomogeneous plasma. The role of associated nonlinear dissipative force in this effect in a confined plasma has been analyzed. The nonlinear force, which arises as a result of resonant interaction between electrons and modulated fields, is shown to drive the instability. Using the ion fluid equation and the ion equation of continuity, the nonlinear dispersion relation of a test ion acoustic wave has been derived, and the growth rate of ion acoustic wave in presence of low frequency drift wave turbulence has been estimated using Helimak data.     

Generation and amplification of magnetic islands by drift interchange turbulence

We investigate the multiscale nonlinear dynamics of a linearly stable or unstable tearing mode with small-scale interchange turbulence using 2D MHD numerical simulations. For a stable tearing mode, the nonlinear beating of the fastest growing small-scale interchange modes drives a magnetic island with an enhanced growth rate to a saturated size that is proportional to the turbulence generated anomalous diffusion. For a linearly unstable tearing mode the island saturation size scales inversely as one-fourth power of the linear tearing growth rate in accordance with weak turbulence theory predictions. Turbulence is also seen to introduce significant modifications in the flow patterns surrounding the magnetic island.     

Effects of intermittency and stratification on the evaluation of optical propagation

To evaluate the effects of intermittency and stratification on propagation of laser beams, we analyze experimental radiosonde profiles using complete ensemble empirical mode decomposition. First, we extract intermittency and stratification at different scales from the original C_n~2 profiles. Second, we establish an intrinsic optical turbulence model(IOTM), which includes fine structures. Third, we calculate several turbulence moments using IOTM. The results quantitatively evaluate the effects of variations at different scales on typical parameters relevant to optical turbulence effects.     

Imprint of large-scale flows on turbulence

We investigate the locality of interactions in hydrodynamic turbulence using data from a direct numerical simulation on a grid of 1024(3) points; the flow is forced with the Taylor-Green vortex. An inertial range for the energy is obtained in which the flux is constant and the spectrum follows an approximate Kolmogorov law. Nonlinear triadic interactions are dominated by their nonlocal components, involving widely separated scales. The resulting nonlinear transfer itself is local at each scale but the step in the energy cascade is independent of that scale and directly related to the integral scale of the flow. Interactions with large scales represent 20% of the total energy flux. Possible explanations for the deviation from self-similar models, the link between these findings and intermittency, and their consequences for modeling of turbulent flows are briefly discussed.     

A nonlinear wave coupling algorithm and its programing and application in plasma turbulences

The fully developed turbulence can be regarded as a nonlinear system, with wave coupling inside, which causes the nonlinear energy to transfer, and drives the turbulence to develop further or be suppressed. Spectral analysis is one of the most effective methods to study turbulence system. In order to apply it to the study of the nonlinear wave coupling process of edge plasma turbulence, an efficient algorithm based on spectral analysis technology is proposed to solve the nonlinear wave coupling equation. The algorithm is based on a mandatory temporal static condition with the nonideal spectra separated from the ideal spectra. The realization idea and programing flow are given. According to the characteristics of plasma turbulence, the simulation data are constructed and used to verify the algorithm and its implementation program. The simulation results and experimental results show the accuracy of the algorithm and the corresponding program, which can play a great role in the studying the energy transfer in edge plasma turbulences. As an application, the energy cascade analysis of typical edge plasma turbulence is carried out by using the results of a case calculation. Consequently, a physical picture of the energy transfer in a kind of fully developed turbulence is constructed, which confirms that the energy transfer in this turbulent system develops from lower-frequency region to higher-frequency region and from linear growing wave to damping wave.     

Density Fluctuation Spectrum of Solar Wind Turbulence between Ion and Electron Scales

We present a measurement of the spectral index of density fluctuations between ion and electron scales in solar wind turbulence using the EFI instrument on the ARTEMIS spacecraft. The mean spectral index at 1?AU was found to be -2.75±0.06, steeper than predictions for pure whistler or kinetic Alfvén wave turbulence but consistent with previous magnetic field measurements. The steep spectra are also consistent with expectations of increased intermittency or damping of some of the turbulent energy over this range of scales. Neither the spectral index nor the flattening of the density spectra before ion scales were found to depend on the proximity to the pressure anisotropy instability thresholds, suggesting that they are features inherent to the turbulent cascade.     

Suppression of turbulence by self-generated and imposed mean flows

The first direct experimental evidence of the suppression of quasi-two-dimensional turbulence by mean flows is presented. The flow either is induced externally or appears in the process of spectral condensation due to an inverse cascade in bounded turbulence. The observed suppression of large scales is consistent with an expected reduction in the correlation time of turbulent eddies due to shearing. At high flow velocities, sweeping of the forcing-scale vortices reduces the energy input, leading to a reduction in the turbulence level.     

Energy spectra stemming from interactions of Alfvén waves and turbulent eddies

We present a numerical analysis of an incompressible decaying magnetohydrodynamic turbulence run on a grid of 1536{3} points. The Taylor Reynolds number at the maximum of dissipation is approximately 1100, and the initial condition is a superposition of large-scale Arn'old-Beltrami-Childress flows and random noise at small scales, with no uniform magnetic field. The initial kinetic and magnetic energies are equal, with negligible correlation. The resulting energy spectrum is a combination of two components, each moderately resolved. Isotropy obtains in the large scales, with a spectral law compatible with the Iroshnikov-Kraichnan theory stemming from the weakening of nonlinear interactions due to Alfvén waves; scaling of structure functions confirms the non-Kolmogorovian nature of the flow in this range. At small scales, weak turbulence emerges with a k{perpendicular}{-2} spectrum, the perpendicular direction referring to the local quasiuniform magnetic field.     

Free energy cascade in gyrokinetic turbulence

In gyrokinetic theory, the quadratic nonlinearity is known to play an important role in the dynamics by redistributing (in a conservative fashion) the free energy between the various active scales. In the present study, the free energy transfer is analyzed for the case of ion temperature gradient driven turbulence. It is shown that it shares many properties with the energy transfer in fluid turbulence. In particular, one finds a (strongly) local, forward (from large to small scales) cascade of free energy in the plane perpendicular to the background magnetic field. These findings shed light on some fundamental properties of plasma turbulence, and encourage the development of large-eddy-simulation techniques for gyrokinetics.     

On the chaotic nature of turbulence observed in benchmark analysis of nonlinear plasma simulation

Simulational results of two dissipative interchange turbulence (Rayleigh-Taylor-type instability with dissipation) models with the same physics are compared. The convective nonlinearity is the nonlinear mechanism in the models. They are shown to have different time evolutions in the nonlinear phase due to the different initial value which is attributed to the initial noise. In the first model (A), a single pressure representing the sum of the ion and electron components is used (one-fluid model). In the second model (B) the ion and electron components of the pressure fields are independently solved (two-fluid model). Both models become physically identical if we set ion and electron pressure fields to be equal in the model (B). The initial conditions only differ by the infinitesimally small initial noise due to the roundoff errors which comes from the finite difference but not the differentiation. This noise grows in accordance with the nonlinear development of the turbulence mode. Interaction with an intrinsic nonlinearity of the system makes the noise grow, whose contribution becomes almost the same magnitude of the fluctuation itself in the results. The instantaneous deviation shows the chaotic characteristics of the turbulence. (c) 1997 American Institute of Physics.     

A speculative study of 23-order fractional Laplacian modeling of turbulence: some thoughts and conjectures

This study makes the first attempt to use the 23-order fractional Laplacian modeling of Kolmogorov -53 scaling of fully developed turbulence and enhanced diffusing movements of random turbulent particles. Nonlinear inertial interactions and molecular Brownian diffusivity are considered to be the bifractal mechanism behind multifractal scaling of moderate Reynolds number turbulence. Accordingly, a stochastic equation is proposed to describe turbulence intermittency. The 23-order fractional Laplacian representation is also used to model nonlinear interactions of fluctuating velocity components, and then we conjecture a fractional Reynolds equation, underlying fractal spacetime structures of Levy 23 stable distribution and the Kolmogorov scaling at inertial scales. The new perspective of this study is that the fractional calculus is an effective approach to modeling the chaotic fractal phenomena induced by nonlinear interactions.     

Identifying coherent structures in nonlinear wave propagation

Nonlinear wave phenomena are often characterized by the appearance of "solitary wave coherent structures" traveling at speeds determined by their amplitudes and morphologies. Assuming that time intervals exist in which these structures are essentially noninteracting, a method for identifying the number of independent features and their respective speeds is proposed and developed. The method is illustrated with a variety of increasingly realistic specific applications, beginning with a simple nonlinear but analytically tractable Gaussian model, continuing with (numerically generated) data describing multisoliton solutions to the Korteweg-de Vries equation, and concluding with (numerical) data from a realistic simulation of nonlinear wave interactions in plasma turbulence. These studies reveal both strengths and limitations of the method in its present incarnation and suggest topics for future investigations.     

Magnetosheath plasma turbulence and its spatiotemporal evolution as observed by the cluster spacecraft

We study the plasma turbulence, at scales larger than the ion inertial length scale, downstream of a quasiparallel bow shock using Cluster multispacecraft measurements. We show that turbulence is intermittent and well described by the extended structure function model, which takes into account the spatial inhomogeneity of the cascade rate. For the first time we use multispacecraft observations to characterize the evolution of magnetosheath turbulence, particularly its intermittency, as a function of the distance from the bow shock. The intermittency significantly changes over the distance of the order of 100 ion inertial lengths, being increasingly stronger and anisotropic away from the bow shock.