Three dimensional particle-in-cell simulation of particle acceleration by circularly polarised inertial Alfven waves in a transversely inhomogeneous plasma

The process of particle acceleration by left-hand, circularly polarised inertial Alfven waves (IAW) in a transversely inhomogeneous plasma is studied using 3D particle-in-cell simulation. A cylindrical tube with, transverse to the background magnetic field, inhomogeneity scale of the order of ion inertial length is considered on which IAWs with frequency $0.3 \omega_{ci}$ are launched that are allowed to develop three wavelength. As a result time-varying parallel electric fields are generated in the density gradient regions which accelerate electrons in the parallel to magnetic field direction. Driven perpendicular electric field of IAWs also heats ions in the transverse direction. Such numerical setup is relevant for solar flaring loops and earth auroral zone. This first, 3D, fully-kinetic simulation demonstrates electron acceleration efficiency in the density inhomogeneity regions, along the magnetic field, of the order of 45% and ion heating, in the transverse to the magnetic field direction, of 75%. The latter is a factor of two times higher than the previous 2.5D analogous study and is in accordance with solar flare particle acceleration observations. We find that the generated parallel electric field is localised in the density inhomogeneity region and rotates in the same direction and with the same angular frequency as the initially launched IAW. Our numerical simulations seem also to suggest that the "knee" often found in the solar flare electron spectra can alternatively be interpreted as the Landau damping (Cerenkov resonance effect) of IAWs due to the wave-particle interactions.Comment: Physics of Plasmas, in-press, September 2012 issue, final accepted versio

Modelling the measured local time evolution of strongly nonlinear heat pulses in the Large Helical Device

In some magnetically confined plasmas, an applied pulse of rapid edge cooling can trigger either a positive or negative excursion in the core electron temperature from its steady state value. We present a new model which captures the time evolution of the transient, non-diffusive local dynamics in the core plasma. We show quantitative agreement between this model and recent spatially localized measurements (Inagaki et al 2010 Plasma Phys. Control. Fusion 52 075002) of the local time-evolving temperature pulse in cold pulse propagation experiments in the Large Helical Device

Particle acceleration by circularly and elliptically polarised dispersive Alfven waves in a transversely inhomogeneous plasma in the inertial and kinetic regimes

Dispersive Alfven waves (DAWs) offer, an alternative to magnetic reconnection, opportunity to accelerate solar flare particles. We study the effect of DAW polarisation, L-, R-, circular and elliptical, in different regimes inertial and kinetic on the efficiency of particle acceleration. We use 2.5D PIC simulations to study how particles are accelerated when DAW, triggered by a solar flare, propagates in transversely inhomogeneous plasma that mimics solar coronal loop. (i) In inertial regime, fraction of accelerated electrons (along the magnetic field), in density gradient regions is ~20% by the time when DAW develops 3 wavelengths and is increasing to ~30% by the time DAW develops 13 wavelengths. In all considered cases ions are heated in transverse to the magnetic field direction and fraction of the heated particles is ~35%. (ii) The case of R-circular, L- and R- elliptical polarisation DAWs, with the electric field in the non-ignorable transverse direction exceeding several times that of in the ignorable direction, produce more pronounced parallel electron beams and transverse ion beams in the ignorable direction. In the inertial regime such polarisations yield the fraction of accelerated electrons ~20%. In the kinetic regime this increases to ~35%. (iii) The parallel electric field that is generated in the density inhomogeneity regions is independent of m_i/m_e and exceeds the Dreicer value by 8 orders of magnitude. (iv) Electron beam velocity has the phase velocity of the DAW. Thus electron acceleration is via Landau damping of DAWs. For the Alfven speeds of 0.3c the considered mechanism can accelerate electrons to energies circa 20 keV. (v) The increase of mass ratio from m_i/m_e=16 to 73.44 increases the fraction of accelerated electrons from 20% to 30-35% (depending on DAW polarisation). For the mass ratio m_i/m_e=1836 the fraction of accelerated electrons would be >35%.Comment: Final accepted version. To appear in Physics of Plasmas, volume 18, issue 9 (September 2011

Robustness of predator-prey models for confinement regime transitions in fusion plasmas

Energy transport and confinement in tokamak fusion plasmas is usually determined by the coupled nonlinear interactions of small-scale drift turbulence and larger scale coherent nonlinear structures, such as zonal flows, together with free energy sources such as temperature gradients. Zero-dimensional models, designed to embody plausible physical narratives for these interactions, can help to identify the origin of enhanced energy confinement and of transitions between confinement regimes. A prime zero-dimensional paradigm is predator-prey or Lotka-Volterra. Here, we extend a successful three-variable (temperature gradient; microturbulence level; one class of coherent structure) model in this genre [M. A. Malkov and P. H. Diamond, Phys. Plasmas 16, 012504 (2009)], by adding a fourth variable representing a second class of coherent structure. This requires a fourth coupled nonlinear ordinary differential equation. We investigate the degree of invariance of the phenomenology generated by the model of Malkov and Diamond, given this additional physics. We study and compare the long-time behaviour of the three-equation and four-equation systems, their evolution towards the final state, and their attractive fixed points and limit cycles. We explore the sensitivity of paths to attractors. It is found that, for example, an attractive fixed point of the three-equation system can become a limit cycle of the four-equation system. Addressing these questions which we together refer to as “robustness” for convenience is particularly important for models which, as here, generate sharp transitions in the values of system variables which may replicate some key features of confinement transitions. Our results help to establish the robustness of the zero-dimensional model approach to capturing observed confinement phenomenology in tokamak fusion plasmas

Quantifying fusion born ion populations in magnetically confined plasmas using ion cyclotron emission

Ion cyclotron emission (ICE) offers unique promise as a diagnostic of the fusion born alpha-particle population in magnetically confined plasmas. Pioneering observations from JET and TFTR found that ICE intensity $P_{ICE}$ scales approximately linearly with the measured neutron flux from fusion reactions, and with the inferred concentration, $n_\alpha/n_i$, of fusion-born alpha-particles confined within the plasma. We present fully nonlinear self-consistent kinetic simulations that reproduce this scaling for the first time. This resolves a longstanding question in the physics of fusion alpha-particle confinement and stability in MCF plasmas. It confirms the magnetoacoustic cyclotron instability (MCI) as the likely emission mechanism and greatly strengthens the basis for diagnostic exploitation of ICE in future burning plasmas

Transitions to improved confinement regimes induced by changes in heating in zero-dimensional models for tokamak plasmas

It is shown that rapid substantial changes in heating rate can induce transitions to improved energy confinement regimes in zero-dimensional models for tokamak plasma phenomenology. We examine for the first time the effect of step changes in heating rate in the models of E-J.Kim and P.H.Diamond, Phys.Rev.Lett. 90, 185006 (2003) and M.A.Malkov and P.H.Diamond, Phys.Plasmas 16, 012504 (2009) which nonlinearly couple the evolving temperature gradient, micro-turbulence and a mesoscale flow; and in the extension of H.Zhu, S.C.Chapman and R.O.Dendy, Phys.Plasmas 20, 042302 (2013), which couples to a second mesoscale flow component. The temperature gradient rises, as does the confinement time defined by analogy with the fusion context, while micro-turbulence is suppressed. This outcome is robust against variation of heating rise time and against introduction of an additional variable into the model. It is also demonstrated that oscillating changes in heating rate can drive the level of micro-turbulence through a period-doubling path to chaos, where the amplitude of the oscillatory component of the heating rate is the control parameter.Comment: 8 pages, 14 figure