This complete introduction to the use of modern ray-tracing techniques in plasma physics describes the powerful mathematical methods generally applicable to vector wave equations in nonuniform media, and clearly demonstrates the application of these methods to simplify and solve important problems in plasma wave theory. Key analytical concepts are carefully introduced as needed, encouraging the development of a visual intuition for the underlying methodology, with more advanced mathematical concepts succinctly explained in the appendices, and supporting MATLAB code available online. Covering variational principles, covariant formulations, caustics, tunneling, mode conversion, weak dissipation, wave emission from coherent sources, incoherent wave fields, and collective wave absorption and emission, all within an accessible framework using standard plasma physics notation, this is an invaluable resource for graduate students and researchers in plasma physics.

Numerical studies of magnetohydrodynamic (MHD) instabilities with feedback control in reversed field pinches (RFPs) are presented. Specifically, investigations are performed of the stability of m=1 modes in RFPs with control based on sensing the normal and tangential magnetic fields at the resistive wall and applying two-parameter feedback proportional to these fields. The control scheme is based on that of [J. M. Finn, Phys. Plasmas 13, 082504 (2006)], which is here modified to use a more realistic plasma model. The plasma model now uses full resistive MHD rather than reduced MHD, and it uses three realistic classes of equilibrium parallel current density profiles appropriate to RFPs. Results with these modifications are in qualitative agreement with [J. M. Finn, Phys. Plasmas 13, 082504 (2006)]: the feedback can stabilize tearing modes (with resistive or ideal-wall) and resistive wall ideal modes. The limit for stabilization is again found to be near the threshold for ideal modes with an ideal-wall. In addition to confirming these predictions, the nature of the instabilities limiting the range of feedback stabilization near the ideal-wall ideal-plasma threshold are studied, and the effects of viscosity, resistive wall time, and plasma resistivity are reported.

No abstract

Abstract. In recent decades, there have been many attempts to construct symplectic integrators with variable time steps, with rather disappointing results. In this paper we identify the causes for this lack of performance, and find that they fall into two categories. In the first, the time step is considered a function of time alone, ∆ = ∆(t). In this case, backwards error analysis shows that while the algorithms remain symplectic, parametric instabilities arise because of resonance between oscillations of ∆(t) and the orbital motion. In the second category the time step is a function of phase space variables ∆ = ∆(q, p). In this case, the system of equations to be solved is analyzed by introducing a new time variable τ with dt = ∆(q, p)dτ. The transformed equations are no longer in Hamiltonian form, and thus are not guaranteed to be stable even when integrated using a method which is symplectic for constant ∆. We analyze two methods for integrating the transformed equations which do, however, preserve the structure of the original equations. The first is an extended phase space method, which has been successfully used in previous studies of adaptive time step symplectic integrators. The second, novel, method is based on a non-canonical mixed-variable generating function. Numerical trials for both of these methods show good results, without parametric instabilities or spurious growth or damping. It is then shown how to adapt the time step to an error estimate found by backward error analysis, in order to optimize the time-stepping scheme. Numerical results are obtained using this formulation and compared with other time-stepping schemes for the extended phase space symplectic method.

In three-dimensional magnetic configurations for a plasma in which no closed field line or magnetic null exists, no magnetic reconnection can occur, by the strictest definition of reconnection. A finitely long pinch with line-tied boundary conditions, in which all the magnetic field lines start at one end of the system and proceed to the opposite end, is an example of such a system. Nevertheless, for a long system of this type, the physical behavior in resistive magnetohydrodynamics (MHD) essentially involves reconnection. This has been explained in terms comparing the geometric and tearing widths [1,2]. The concept of a quasi-separatrix layer [3,4] was developed for such systems. In this paper we study a model for a line-tied system in which the corresponding periodic system has an unstable tearing mode. We analyze this system in terms of two magnetic field line diagnostics, the squashing factor [3][4][5] and the electrostatic potential difference used in kinematic reconnection studies [6,7]. We discuss the physical and geometric significance of these two diagnostics and compare them in the context of discerning tearing-like behavior in line-tied modes.