Bifurcation in electrostatic resistive drift wave turbulence

The Hasegawa-Wakatani equations, coupling plasma density and electrostatic potential through an approximation to the physics of parallel electron motions, are a simple model that describes resistive drift wave turbulence. We present numerical analyses of bifurcation phenomena in the model that provide new insights into the interactions between turbulence and zonal flows in the tokamak plasma edge region. The simulation results show a regime where, after an initial transient, drift wave turbulence is suppressed through zonal flow generation. As a parameter controlling the strength of the turbulence is tuned, this zonal flow dominated state is rapidly destroyed and a turbulence-dominated state re-emerges. The transition is explained in terms of the Kelvin-Helmholtz stability of zonal flows. This is the first observation of an upshift of turbulence onset in the resistive drift wave system, which is analogous to the well-known Dimits shift in turbulence driven by ion temperature gradients.Comment: 21 pages, 11 figure

Strong "quantum" chaos in the global ballooning mode spectrum of three-dimensional plasmas

The spectrum of ideal magnetohydrodynamic (MHD) pressure-driven (ballooning) modes in strongly nonaxisymmetric toroidal systems is difficult to analyze numerically owing to the singular nature of ideal MHD caused by lack of an inherent scale length. In this paper, ideal MHD is regularized by using a $k$-space cutoff, making the ray tracing for the WKB ballooning formalism a chaotic Hamiltonian billiard problem. The minimum width of the toroidal Fourier spectrum needed for resolving toroidally localized ballooning modes with a global eigenvalue code is estimated from the Weyl formula. This phase-space-volume estimation method is applied to two stellarator cases.Comment: 4 pages typeset, including 2 figures. Paper accepted for publication in Phys. Rev. Letter

A comparison of incompressible limits for resistive plasmas

The constraint of incompressibility is often used to simplify the magnetohydrodynamic (MHD) description of linearized plasma dynamics because it does not affect the ideal MHD marginal stability point. In this paper two methods for introducing incompressibility are compared in a cylindrical plasma model: In the first method, the limit $\gamma \to \infty$ is taken, where $\gamma$ is the ratio of specific heats; in the second, an anisotropic mass tensor $\mathbf{\rho}$ is used, with the component parallel to the magnetic field taken to vanish, $\rho_{\parallel} \to 0$. Use of resistive MHD reveals the nature of these two limits because the Alfv\'en and slow magnetosonic continua of ideal MHD are converted to point spectra and moved into the complex plane. Both limits profoundly change the slow-magnetosonic spectrum, but only the second limit faithfully reproduces the resistive Alfv\'en spectrum and its wavemodes. In ideal MHD, the slow magnetosonic continuum degenerates to the Alfv\'en continuum in the first method, while it is moved to infinity by the second. The degeneracy in the first is broken by finite resistivity. For numerical and semi-analytical study of these models, we choose plasma equilibria which cast light on puzzling aspects of results found in earlier literature.Comment: 14 pages, 10 figure

Additivity of Bond Energies in the Light of the Maximum Overlap Approximation (MOA) and MIND0/3

Optimum parameters are determined for the use of the maximum overlap approximation (MOA) to calculate heats of atomization of a varied set of eleven hydrocarbons. The agreement with experiment is generally good. The reasons for this success are discussed in terms of an analysis by energy partitioning of MIND0/3 calculations for the same hydrocarbons

Dressed test particles, oscillation centres and pseudo-orbits

A general semi-analytical method for accurate and efficient numerical calculation of the dielectrically screened ("dressed") potential around a non-relativistic test particle moving in an isotropic, collisionless, unmagnetised plasma is presented. The method requires no approximations and is illustrated using results calculated for two cases taken from the MSc thesis of the first author: test particles with velocities above and below the ion sound speed in plasmas with Maxwellian ions and warm electrons. The idea that the fluctuation spectrum of a plasma can be described as a superposition of the fields around \emph{non-interacting} dressed test particles is an expression of the quasiparticle concept, which has also been expressed in the development of the oscillation-centre and pseudo-orbit formalisms.Comment: 14 pages to Plasma Physics and Controlled Fusion for publication with a cluster of papers associated with workshop Stability and Nonlinear Dynamics of Plasmas, October 31, 2009 Atlanta, GA on occasion of the 65th birthday of R.L. Dewar. Version 2: Reference [27] added in Sec. 5. Version 3: Revised in response to referee

Analysing and modelling train driver performance

Arguments for the importance of contextual factors in understanding human performance have been made extremely persuasive in the context of the process control industries. This paper puts these arguments into the context of the train driving task, drawing on an extensive analysis of driver performance with the Automatic Warning System (AWS). The paper summarises a number of constructs from applied psychological research which are thought to be important in understanding train driver performance. A “Situational Model” is offered as a framework for investigating driver performance. The model emphasises the importance of understanding the state of driver cognition at a specific time (“Now”) in a specific situation and a specific context

Singularity theory study of overdetermination in models for L-H transitions

Two dynamical models that have been proposed to describe transitions between low and high confinement states (L-H transitions) in confined plasmas are analysed using singularity theory and stability theory. It is shown that the stationary-state bifurcation sets have qualitative properties identical to standard normal forms for the pitchfork and transcritical bifurcations. The analysis yields the codimension of the highest-order singularities, from which we find that the unperturbed systems are overdetermined bifurcation problems and derive appropriate universal unfoldings. Questions of mutual equivalence and the character of the state transitions are addressed.Comment: Latex (Revtex) source + 13 small postscript figures. Revised versio