-
Two-fluid simulations of driven reconnection in the Mega-Ampere Spherical Tokamak
Authors:
A. Stanier,
P. Browning,
M. Gordovskyy,
K. G. McClements,
M. P. Gryaznevich,
V. S. Lukin
Abstract:
In the merging-compression method of plasma start-up, two flux-ropes with parallel toroidal current are formed around in-vessel poloidal field coils, before merging to form a spherical tokamak plasma. This start-up method, used in the Mega-Ampere Spherical Tokamak (MAST), is studied as a high Lundquist number and low plasma-beta magnetic reconnection experiment.
In this paper, 2D fluid simulatio…
▽ More
In the merging-compression method of plasma start-up, two flux-ropes with parallel toroidal current are formed around in-vessel poloidal field coils, before merging to form a spherical tokamak plasma. This start-up method, used in the Mega-Ampere Spherical Tokamak (MAST), is studied as a high Lundquist number and low plasma-beta magnetic reconnection experiment.
In this paper, 2D fluid simulations are presented of this merging process in order to understand the underlying physics, and better interpret the experimental data. These simulations examine the individual and combined effects of tight-aspect ratio geometry and two-fluid physics on the merging. The ideal self-driven flux-rope dynamics are coupled to the diffusion layer physics, resulting in a large range of phenomena. For resistive MHD simulations, the flux-ropes enter the sloshing regime for normalised resistivity eta < 1E-5. In Hall-MHD three regimes are found for the qualitative behaviour of the current sheet, depending on the ratio of the current sheet width to the ion-sound radius. These are a stable collisional regime, an open X-point regime, and an intermediate regime that is highly unstable to tearing-type instabilities.
In toroidal axisymmetric geometry, the final state after merging is a MAST-like spherical tokamak with nested flux-surfaces. It is also shown that the evolution of simulated 1D radial density profiles closely resembles the Thomson scattering electron density measurements in MAST. An intuitive explanation for the origin of the measured density structures is proposed, based upon the results of the toroidal Hall-MHD simulations.
△ Less
Submitted 13 August, 2013;
originally announced August 2013.
-
Nonlinear wave propagation and reconnection at magnetic X-points in the Hall MHD regime
Authors:
J. W. Threlfall,
C. E. Parnell,
I. De Moortel,
K. G. McClements,
T. D. Arber
Abstract:
The highly dynamical, complex nature of the solar atmosphere naturally implies the presence of waves in a topologically varied magnetic environment. Here, the interaction of waves with topological features such as null points is inevitable and potentially important for energetics. The low resistivity of the solar coronal plasma implies that non-MHD effects should be considered in studies of magnet…
▽ More
The highly dynamical, complex nature of the solar atmosphere naturally implies the presence of waves in a topologically varied magnetic environment. Here, the interaction of waves with topological features such as null points is inevitable and potentially important for energetics. The low resistivity of the solar coronal plasma implies that non-MHD effects should be considered in studies of magnetic energy release in this environment. This paper investigates the role of the Hall term in the propagation and dissipation of waves, their interaction with 2D magnetic X-points and the nature of the resulting reconnection. A Lagrangian remap shock-capturing code (Lare2d) is used to study the evolution of an initial fast magnetoacoustic wave annulus for a range of values of the ion skin depth in resistive Hall MHD. A magnetic null-point finding algorithm is also used to locate and track the evolution of the multiple null-points that are formed in the system. Depending on the ratio of ion skin depth to system size, our model demonstrates that Hall effects can play a key role in the wave-null interaction. In particular, the initial fast-wave pulse now consists of whistler and ion-cyclotron components; the dispersive nature of the whistler wave leads to (i) earlier interaction with the null, (ii) the creation of multiple additional, transient nulls and, hence, an increased number of energy release sites. In the Hall regime, the relevant timescales (such as the onset of reconnection and the period of the oscillatory relaxation) of the system are reduced significantly, and the reconnection rate is enhanced.
△ Less
Submitted 16 May, 2012; v1 submitted 16 February, 2012;
originally announced February 2012.
-
Alfvèn wave phase-mixing and damping in the ion cyclotron range of frequencies
Authors:
J. W. Threlfall,
K. G. McClements,
I. De Moortel
Abstract:
Aims. To determine the effect of the Hall term in the generalised Ohm's law on the damping and phase mixing of Alfven waves in the ion cyclotron range of frequencies in uniform and non-uniform equilibrium plasmas. Methods. Wave damping in a uniform plasma is treated analytically, whilst a Lagrangian remap code (Lare2d) is used to study Hall effects on damping and phase mixing in the presence of an…
▽ More
Aims. To determine the effect of the Hall term in the generalised Ohm's law on the damping and phase mixing of Alfven waves in the ion cyclotron range of frequencies in uniform and non-uniform equilibrium plasmas. Methods. Wave damping in a uniform plasma is treated analytically, whilst a Lagrangian remap code (Lare2d) is used to study Hall effects on damping and phase mixing in the presence of an equilibrium density gradient. Results. The magnetic energy associated with an initially Gaussian field perturbation in a uniform resistive plasma is shown to decay algebraically at a rate that is unaffected by the Hall term to leading order in k^2di^2 where k is wavenumber and di is ion skin depth. A similar algebraic decay law applies to whistler perturbations in the limit k^2di^2>>1. In a non-uniform plasma it is found that the spatially-integrated damping rate due to phase mixing is lower in Hall MHD than it is in MHD, but the reduction in the damping rate, which can be attributed to the effects of wave dispersion, tends to zero in both the weak and strong phase mixing limits.
△ Less
Submitted 8 October, 2010; v1 submitted 27 July, 2010;
originally announced July 2010.
-
Electron Inertial Effects on Rapid Energy Redistribution at Magnetic X-points
Authors:
K. G. McClements,
A. Thyagaraja,
N. Ben Ayed,
L. Fletcher
Abstract:
The evolution of non-potential perturbations to a current-free magnetic X-point configuration is studied, taking into account electron inertial effects as well as resistivity. Electron inertia is shown to have a negligible effect on the evolution of the system whenever the collisionless skin depth is less than the resistive scale length. Non-potential magnetic field energy in this resistive MHD…
▽ More
The evolution of non-potential perturbations to a current-free magnetic X-point configuration is studied, taking into account electron inertial effects as well as resistivity. Electron inertia is shown to have a negligible effect on the evolution of the system whenever the collisionless skin depth is less than the resistive scale length. Non-potential magnetic field energy in this resistive MHD limit initially reaches equipartition with flow energy, in accordance with ideal MHD, and is then dissipated extremely rapidly, on an Alfvenic timescale that is essentially independent of Lundquist number. In agreement with resistive MHD results obtained by previous authors, the magnetic field energy and kinetic energy are then observed to decay on a longer timescale and exhibit oscillatory behavior, reflecting the existence of discrete normal modes with finite real frequency. When the collisionless skin depth exceeds the resistive scale length, the system again evolves initially according to ideal MHD. At the end of this ideal phase, the field energy decays typically on an Alfvenic timescale, while the kinetic energy (which is equally partitioned between ions and electrons in this case) is dissipated on the electron collision timescale. The oscillatory decay in the energy observed in the resistive case is absent, but short wavelength structures appear in the field and velocity profiles, suggesting the possibility of particle acceleration in oppositely-directed current channels. The model provides a possible framework for interpreting observations of energy release and particle acceleration on timescales down to less than a second in the impulsive phase of solar flares.
△ Less
Submitted 27 January, 2004;
originally announced January 2004.
-
Field-guided proton acceleration at reconnecting X-points in flares
Authors:
B. Hamilton,
K. G. McClements,
L. Fletcher,
A. Thyagaraja
Abstract:
An explicitly energy-conserving full orbit code CUEBIT, developed originally to describe energetic particle effects in laboratory fusion experiments, has been applied to the problem of proton acceleration in solar flares. The model fields are obtained from solutions of the linearised MHD equations for reconnecting modes at an X-type neutral point, with the additional ingredient of a longitudinal…
▽ More
An explicitly energy-conserving full orbit code CUEBIT, developed originally to describe energetic particle effects in laboratory fusion experiments, has been applied to the problem of proton acceleration in solar flares. The model fields are obtained from solutions of the linearised MHD equations for reconnecting modes at an X-type neutral point, with the additional ingredient of a longitudinal magnetic field component. To accelerate protons to the highest observed energies on flare timescales, it is necessary to invoke anomalous resistivity in the MHD solution. It is shown that the addition of a longitudinal field component greatly increases the efficiency of ion acceleration, essentially because it greatly reduces the magnitude of drift motions away from the vicinity of the X-point, where the accelerating component of the electric field is largest. Using plasma parameters consistent with flare observations, we obtain proton distributions extending up to gamma-ray-emitting energies (>1MeV). In some cases the energy distributions exhibit a bump-on-tail in the MeV range. In general, the shape of the distribution is sensitive to the model parameters.
△ Less
Submitted 21 February, 2003; v1 submitted 17 January, 2003;
originally announced January 2003.
-
Azimuthally symmetric MHD and two-fluid equilibria with arbitrary flows
Authors:
K. G. McClements,
A. Thyagaraja
Abstract:
Magnetohydrodynamic (MHD) and two-fluid quasi-neutral equilibria with azimuthal symmetry, gravity and arbitrary ratios of (nonrelativistic) flow speed to acoustic and Alfven speeds are investigated. In the two-fluid case, the mass ratio of the two species is arbitrary, and the analysis is therefore applicable to electron-positron plasmas. The methods of derivation can be extended in an obvious m…
▽ More
Magnetohydrodynamic (MHD) and two-fluid quasi-neutral equilibria with azimuthal symmetry, gravity and arbitrary ratios of (nonrelativistic) flow speed to acoustic and Alfven speeds are investigated. In the two-fluid case, the mass ratio of the two species is arbitrary, and the analysis is therefore applicable to electron-positron plasmas. The methods of derivation can be extended in an obvious manner to several charged species. Generalized Grad-Shafranov equations, describing the equilibrium magnetic field, are derived. Flux function equations and Bernoulli relations for each species, together with Poisson's equation for the gravitational potential, complete the set of equations required to determine the equilibrium. These are straightforward to solve numerically. The two-fluid system, unlike the MHD system, is shown to be free of singularities. It is demonstrated analytically that there exists a class of incompressible MHD equilibria with magnetic field-aligned flow. A special sub--class first identified by S. Chandrasekhar, in which the flow speed is everywhere equal to the local Alfven speed, is compatible with virtually any azimuthally symmetric magnetic configuration. Potential applications of this analysis include extragalactic and stellar jets, and accretion disks.
△ Less
Submitted 6 November, 2000; v1 submitted 8 June, 2000;
originally announced June 2000.
-
Electron acceleration due to high frequency instabilities at supernova remnant shocks
Authors:
M. E. Dieckmann,
K. G. McClements,
S. C. Chapman,
R. O. Dendy,
L. O'C. Drury
Abstract:
Observations of synchrotron radiation across a wide range of wavelengths provide clear evidence that electrons are accelerated to relativistic energies in supernova remnants (SNRs). However, a viable mechanism for the pre-acceleration of such electrons to mildly relativistic energies has not yet been established. In this paper an electromagnetic particle-in-cell (PIC) code is used to simulate ac…
▽ More
Observations of synchrotron radiation across a wide range of wavelengths provide clear evidence that electrons are accelerated to relativistic energies in supernova remnants (SNRs). However, a viable mechanism for the pre-acceleration of such electrons to mildly relativistic energies has not yet been established. In this paper an electromagnetic particle-in-cell (PIC) code is used to simulate acceleration of electrons from background energies to tens of keV at perpendicular collisionless shocks associated with SNRs. Free energy for electron energization is provided by ions reflected from the shock front, with speeds greater than the upstream electron thermal speed. The PIC simulation results contain several new features, including: the acceleration, rather than heating, of electrons via the Buneman instability; the acceleration of electrons to speeds exceeding those of the shock-reflected ions producing the instability; and strong acceleration of electrons perpendicular to the magnetic field. Electron energization takes place through a variety of resonant and non-resonant processes, of which the strongest involves stochastic wave-particle interactions. In SNRs the diffusive shock process could then supply the final step required for the production of fully relativistic electrons. The mechanisms identified in this paper thus provide a possible solution to the electron pre-acceleration problem.
△ Less
Submitted 17 February, 2000;
originally announced February 2000.