The eruption of multiple flux tubes in a magnetised plasma atmosphere is proposed as a mechanism for explosive release of energy in plasmas. Linearly stable isolated flux tubes are shown to be metastable in a box model magnetised atmosphere in which
ends of the field lines are embedded in conducting walls. The energy released by destabilising such field lines can be a significant fraction of the gravitational energy stored in the system. This energy can be released in a fast dynamical time.
Calculations of tearing mode stability in tokamaks split conveniently into an external region, where marginally stable ideal MHD is applicable, and a resonant layer around the rational surface where sophisticated kinetic physics is needed. These two
regions are coupled by the stability parameter. Pressure and current perturbations localized around the rational surface alter the stability of tearing modes. Equations governing the changes in the external solution and - are derived for arbitrary perturbations in axisymmetric toroidal geometry. The relationship of - with and without pressure flattening is obtained analytically for four pressure flattening functions. Resistive MHD codes do not contain the appropriate layer physics and therefore cannot predict stability directly. They can, however, be used to calculate -. Existing methods (Ham et al. 2012 Plasma Phys. Control. Fusion 54 025009) for extracting - from resistive codes are unsatisfactory when there is a finite pressure gradient at the rational surface and favourable average curvature because of the Glasser stabilizing effect. To overcome this difficulty we introduce a specific pressure flattening function that allows the earlier approach to be used. The technique is first tested numerically in cylindrical geometry with an artificial favourable curvature. Its application to toroidal geometry is then demonstrated using the toroidal tokamak tearing mode stability code T7 (Fitzpatrick et al. 1993 Nucl. Fusion 33 1533) which uses an approximate analytic equilibrium. The prospects for applying this approach to resistive MHD codes such as MARS-F (Liu et al. 2000 Phys. Plasmas 7 3681) which utilize a fully toroidal equilibrium are discussed.
We have extended our study of the competition between the drive and stabilization of plasma microinstabilities by sheared flow to include electromagnetic effects at low plasma $beta$ (the ratio of plasma to magnetic pressure). The extended system of
characteristic equations is formulated, for a dissipative fluid model developed from the gyrokinetic equation, using a twisting mode representation in sheared slab geometry and focusing on the ion temperature gradient mode. Perpendicular flow shear convects perturbations along the field at the speed we denote as $Mc_s$ (where $c_s$ is the sound speed). $M > 1/ sqrt{beta}$ is required to make the system characteristics unidirectional and inhibit eigenmode formation, leaving only transitory perturbations in the system. This typically represents a much larger flow shear than in the electrostatic case, which only needs $M>1$. Numerical investigation of the region $M < 1/sqrt{beta}$ shows the driving terms can conflict, as in the electrostatic case, giving low growth rates over a range of parameters. Also, at modest drive strengths and low $beta$ values typical of experiments, including electromagnetic effects does not significantly alter the growth rates. For stronger flow shear and higher $beta$, geometry characteristic of the spherical tokamak mitigates the effect of an instability of the shear Alfv{e}n wave, driven by the parallel flow shear.
We consider the problem of self-regulated heating and cooling in galaxy clusters and the implications for cluster magnetic fields and turbulence. Viscous heating of a weakly collisional magnetised plasma is regulated by the pressure anisotropy with r
espect to the local direction of the magnetic field. The intracluster medium is a high-beta plasma, where pressure anisotropies caused by the turbulent stresses and the consequent local changes in the magnetic field will trigger very fast microscale instabilities. We argue that the net effect of these instabilities will be to pin the pressure anisotropies at a marginal level, controlled by the plasma beta parameter. This gives rise to local heating rates that turn out to be comparable to the radiative cooling rates. Furthermore, we show that a balance between this heating and Bremsstrahlung cooling is thermally stable, unlike the often conjectured balance between cooling and thermal conduction. Given a sufficient (and probably self-regulating) supply of turbulent power, this provides a physical mechanism for mitigating cooling flows and preventing cluster core collapse. For observed density and temperature profiles, the assumed balance of viscous heating and radiative cooling allows us to predict magnetic-field strengths, turbulent velocities and turbulence scales as functions of distance from the centre. Specific predictions and comparisons with observations are given for several different clusters. Our predictions can be further tested by future observations of cluster magnetic fields and turbulent velocities.
A new analytically and numerically manageable model collision operator is developed specifically for turbulence simulations. The like-particle collision operator includes both pitch-angle scattering and energy diffusion and satisfies the physical con
straints required for collision operators: it conserves particles, momentum and energy, obeys Boltzmanns H-theorem (collisions cannot decrease entropy), vanishes on a Maxwellian, and efficiently dissipates small-scale structure in the velocity space. The process of transforming this collision operator into the gyroaveraged form for use in gyrokinetic simulations is detailed. The gyroaveraged model operator is shown to have more suitable behavior at small scales in phase space than previously suggested models. A model operator for electron-ion collisions is also presented.
Howes et al. Reply to Comment on Kinetic Simulations of Magnetized Turbulence in Astrophysical Plasmas arXiv:0711.4355
This paper describes a conceptual framework for understanding kinetic plasma turbulence as a generalized form of energy cascade in phase space. It is emphasized that conversion of turbulent energy into thermodynamic heat is only achievable in the pre
sence of some (however small) degree of collisionality. The smallness of the collision rate is compensated by the emergence of small-scale structure in the velocity space. For gyrokinetic turbulence, a nonlinear perpendicular phase mixing mechanism is identified and described as a turbulent cascade of entropy fluctuations simultaneously occurring at spatial scales smaller than the ion gyroscale and in velocity space. Scaling relations for the resulting fluctuation spectra are derived. An estimate for the collisional cutoff is provided. The importance of adequately modeling and resolving collisions in gyrokinetic simulations is biefly discussed, as well as the relevance of these results to understanding the dissipation-range turbulence in the solar wind and the electrostatic microturbulence in fusion plasmas.
In turbulent high-beta astrophysical plasmas (exemplified by the galaxy cluster plasmas), pressure-anisotropy-driven firehose and mirror fluctuations grow nonlinearly to large amplitudes, dB/B ~ 1, on a timescale comparable to the turnover time of th
e turbulent motions. The principle of their nonlinear evolution is to generate secularly growing small-scale magnetic fluctuations that on average cancel the temporal change in the large-scale magnetic field responsible for the pressure anisotropies. The presence of small-scale magnetic fluctuations may dramatically affect the transport properties and, thereby, the large-scale dynamics of the high-beta astrophysical plasmas.
We present a theoretical framework for plasma turbulence in astrophysical plasmas (solar wind, interstellar medium, galaxy clusters, accretion disks). The key assumptions are that the turbulence is anisotropic with respect to the mean magnetic field
and frequencies are low compared to the ion cyclotron frequency. The energy injected at the outer scale scale has to be converted into heat, which ultimately cannot be done without collisions. A KINETIC CASCADE develops that brings the energy to collisional scales both in space and velocity. Its nature depends on the physics of plasma fluctuations. In each of the physically distinct scale ranges, the kinetic problem is systematically reduced to a more tractable set of equations. In the inertial range above the ion gyroscale, the kinetic cascade splits into a cascade of Alfvenic fluctuations, which are governed by the RMHD equations at both the collisional and collisionless scales, and a passive cascade of compressive fluctuations, which obey a linear kinetic equation along the moving field lines associated with the Alfvenic component. In the dissipation range between the ion and electron gyroscales, there are again two cascades: the kinetic-Alfven-wave (KAW) cascade governed by two fluid-like Electron RMHD equations and a passive phase-space cascade of ion entropy fluctuations. The latter cascade brings the energy of the inertial-range fluctuations that was damped by collisionless wave-particle interaction at the ion gyroscale to collisional scales in the phase space and leads to ion heating. The KAW energy is similarly damped at the electron gyroscale and converted into electron heat. Kolmogorov-style scaling relations are derived for these cascades. Astrophysical and space-physical applications are discussed in detail.
