Showing 1–18 of 18 results for author: Pfefferlé, D

Global realisation of magnetic fields as 1$\frac{1}{2}$D Hamiltonian systems

Authors: Nathan Duignan, David Perrella, David Pfefferlé

Abstract: The paper reviews the notion of $n+\frac{1}{2}$D non-autonomous Hamiltonian systems, portraying their dynamics as the flow of the Reeb field related to a closed two-form of maximal rank on a cosymplectic manifold, and naturally decomposing into time-like and Hamiltonian components. The paper then investigates the conditions under which the field-line dynamics of a (tangential) divergence-free vect… ▽ More The paper reviews the notion of $n+\frac{1}{2}$D non-autonomous Hamiltonian systems, portraying their dynamics as the flow of the Reeb field related to a closed two-form of maximal rank on a cosymplectic manifold, and naturally decomposing into time-like and Hamiltonian components. The paper then investigates the conditions under which the field-line dynamics of a (tangential) divergence-free vector field on a connected compact three-manifold (possibly with boundary) diffeomorphic to a trivial fibre bundle over the circle can be conversely identified as a non-autonomous $1\frac{1}{2}$D Hamiltonian system. Under the assumption that the field is transverse to a global compact Poincaré section, an adaptation of Moser's trick shows that all such fields are locally-Hamiltonian. A full identification is established upon further assuming that the Poincaré sections are planar, which crucially implies (together with Dirichlet boundary conditions) that the cohomology class of the generating closed one-forms on each section is constant. By reviewing the classification of fibre bundles over the circle using the monodromy representation, it is remarked that as soon as the Poincaré section of an alleged field is diffeomorphic to a disk or an annulus, the domain is necessarily diffeomorphic to a solid or hollow torus, and thus its field-line dynamics can always be identified as a non-autonomous Hamiltonian system. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: 21 pages (excluding appendix), 1 figure

Structure-preserving particle methods for the Landau collision operator using the metriplectic framework

Authors: Sandra Jeyakumar, Michael Kraus, Matthew J. Hole, David Pfefferlé

Abstract: We present a novel family of particle discretisation methods for the nonlinear Landau collision operator. We exploit the metriplectic structure underlying the Vlasov-Maxwell-Landau system in order to obtain disretisation schemes that automatically preserve mass, momentum, and energy, warrant monotonic dissipation of entropy, and are thus guaranteed to respect the laws of thermodynamics. In contras… ▽ More We present a novel family of particle discretisation methods for the nonlinear Landau collision operator. We exploit the metriplectic structure underlying the Vlasov-Maxwell-Landau system in order to obtain disretisation schemes that automatically preserve mass, momentum, and energy, warrant monotonic dissipation of entropy, and are thus guaranteed to respect the laws of thermodynamics. In contrast to recent works that used radial basis functions and similar methods for regularisation, here we use an auxiliary spline or finite element representation of the distribution function to this end. Discrete gradient methods are employed to guarantee the aforementioned properties in the time discrete domain as well. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 17 pages

Characterization of admissible quasisymmetries

Authors: J. W. Burby, N. Kallinikos, R. S. MacKay, D. Perrella, D. Pfefferlé

Abstract: We solve "half" the problem of finding three-dimensional quasisymmetric magnetic fields that do not necessarily satisfy force balance. This involves determining which hidden symmetries are admissible as quasisymmetries, and then showing explicitly how to construct quasisymmetric magnetic fields given an admissible symmetry. The admissibility conditions take the form of a system of overdetermined n… ▽ More We solve "half" the problem of finding three-dimensional quasisymmetric magnetic fields that do not necessarily satisfy force balance. This involves determining which hidden symmetries are admissible as quasisymmetries, and then showing explicitly how to construct quasisymmetric magnetic fields given an admissible symmetry. The admissibility conditions take the form of a system of overdetermined nonlinear partial differential equations involving second derivatives of the symmetry's infinitesimal generator. △ Less

Submitted 3 July, 2024; v1 submitted 5 March, 2024; originally announced March 2024.

Comments: 18 pages

A structure-preserving particle discretisation for the Lenard-Bernstein collision operator

Authors: Sandra Jeyakumar, Michael Kraus, Matthew Hole, David Pfefferlé

Abstract: Collisions are an important dissipation mechanism in plasmas. In one-dimensional modelling, a commonly used collision operator is the Lenard-Bernstein operator, or its modified energy- and momentum-conserving counterpart. When approximating such operators numerically, it is important to respect their structure in order to satisfy the laws of thermodynamics. It is, however, challenging to discretis… ▽ More Collisions are an important dissipation mechanism in plasmas. In one-dimensional modelling, a commonly used collision operator is the Lenard-Bernstein operator, or its modified energy- and momentum-conserving counterpart. When approximating such operators numerically, it is important to respect their structure in order to satisfy the laws of thermodynamics. It is, however, challenging to discretise such operators in a structure-preserving way when considering particle methods. In this work, we present a macro-particle discretisation of the Lenard-Bernstein collision operator that is energy and momentum preserving. △ Less

Submitted 1 February, 2024; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: 17 pages, 9 figures

Existence of global symmetries of divergence-free fields with first integrals

Authors: David Perrella, Nathan Duignan, David Pfefferlé

Abstract: The relationship between symmetry fields and first integrals of divergence-free vector fields is explored in three dimensions in light of its relevance to plasma physics and magnetic confinement fusion. A Noether-type Theorem is known: for each such symmetry, there corresponds a first integral. The extent to which the converse is true is investigated. In doing so, a reformulation of this Noether-t… ▽ More The relationship between symmetry fields and first integrals of divergence-free vector fields is explored in three dimensions in light of its relevance to plasma physics and magnetic confinement fusion. A Noether-type Theorem is known: for each such symmetry, there corresponds a first integral. The extent to which the converse is true is investigated. In doing so, a reformulation of this Noether-type Theorem is found for which the converse holds on what is called the toroidal region. Some consequences of the methods presented are quick proofs of the existence of flux coordinates for magnetic fields in high generality; without needing to assume a symmetry such as in the cases of magneto-hydrostatics (MHS) or quasi-symmetry. △ Less

Submitted 31 March, 2023; v1 submitted 6 March, 2023; originally announced March 2023.

Comments: 31 pages, 3 figures. This version of the article features an example involving Reeb cylinders whose idea was suggested by Daniel Peralta-Salas

MSC Class: 53Z05; 37C86; 58A30

Quasisymmetric magnetic fields in asymmetric toroidal domains

Authors: Naoki Sato, Zhisong Qu, David Pfefferlé, Robert L. Dewar

Abstract: We explore the existence of quasisymmetric magnetic fields in asymmetric toroidal domains. These vector fields can be identified with a class of magnetohydrodynamic equilibria in the presence of pressure anisotropy. First, using Clebsch potentials, we derive a system of two coupled nonlinear first order partial differential equations expressing a family of quasisymmetric magnetic fields in bounded… ▽ More We explore the existence of quasisymmetric magnetic fields in asymmetric toroidal domains. These vector fields can be identified with a class of magnetohydrodynamic equilibria in the presence of pressure anisotropy. First, using Clebsch potentials, we derive a system of two coupled nonlinear first order partial differential equations expressing a family of quasisymmetric magnetic fields in bounded domains. In regions where flux surfaces and surfaces of constant field strength are not tangential, this system can be further reduced to a single degenerate nonlinear second order partial differential equation with externally assigned initial data. Then, we exhibit regular quasisymmetric vector fields which correspond to local solutions of anisotropic magnetohydrodynamics in asymmetric toroidal domains such that tangential boundary conditions are fulfilled on a portion of the bounding surface. The problems of boundary shape and locality are also discussed. We find that symmetric magnetic fields can be fitted into asymmetric domains, and that the mathematical difficulty encountered in the derivation of global quasisymmetric magnetic fields lies in the topological obstruction toward global extension affecting local solutions of the governing nonlinear first order partial differential equations. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Comments: 19 pages, 5 figures

Journal ref: Physics of Plasmas 28, 112507 (2021)

Combined plasma-coil optimization algorithms

Authors: Sophia A. Henneberg, Stuart R. Hudson, David Pfefferlé, Per Helander

Abstract: Combined plasma-coil optimization approaches for designing stellarators are discussed and a new method for calculating free-boundary equilibria is proposed. Four distinct categories of stellarator optimization, two of which are novel approaches, are the fixed-boundary optimization, the generalized fixed-boundary optimization, the quasi free-boundary optimization, and the free-boundary (coil) optim… ▽ More Combined plasma-coil optimization approaches for designing stellarators are discussed and a new method for calculating free-boundary equilibria is proposed. Four distinct categories of stellarator optimization, two of which are novel approaches, are the fixed-boundary optimization, the generalized fixed-boundary optimization, the quasi free-boundary optimization, and the free-boundary (coil) optimization. These are described using the multi-region relaxed magnetohydrodynmics (MRxMHD) energy functional, the Biot-Savart integral, the coil-penalty functional and the virtual casing integral, and their derivatives. The proposed free-boundary equilibrium calculation differs from existing methods in how the boundary-value problem is posed, and for the new approach it seems that there is not an associated energy minimization principle because a non-symmetric functional arises. We propose to solve the weak formulation of this problem using a spectral-Galerkin method, and this will reduce the free-boundary equilibrium calculation to something comparable to a fixed-boundary calculation. In our discussion of combined plasma-coil optimization algorithms, we emphasize the importance of the stability matrix. △ Less

Submitted 16 December, 2020; originally announced December 2020.

Comments: to be submitted to JPP

Energy and momentum conservation in the Euler-Poincaré formulation of local Vlasov-Maxwell-type systems

Authors: Eero Hirvijoki, Joshua W. Burby, David Pfefferlé, Alain J. Brizard

Abstract: The action principle by Low [Proc. R. Soc. Lond. A 248, 282--287] for the classic Vlasov-Maxwell system contains a mix of Eulerian and Lagrangian variables. This renders the Noether analysis of reparametrization symmetries inconvenient, especially since the well-known energy- and momentum-conservation laws for the system are expressed in terms of Eulerian variables only. While an Euler-Poincaré fo… ▽ More The action principle by Low [Proc. R. Soc. Lond. A 248, 282--287] for the classic Vlasov-Maxwell system contains a mix of Eulerian and Lagrangian variables. This renders the Noether analysis of reparametrization symmetries inconvenient, especially since the well-known energy- and momentum-conservation laws for the system are expressed in terms of Eulerian variables only. While an Euler-Poincaré formulation of Vlasov-Maxwell-type systems, effectively starting with Low's action and using constrained variations for the Eulerian description of particle motion, has been known for a while [J. Math. Phys., 39, 6, pp. 3138-3157], it is hard to come by a documented derivation of the related energy- and momentum-conservation laws in the spirit of the Euler-Poincaré machinery. To our knowledge only one such derivation exists in the literature so far, dealing with the so-called guiding-center Vlasov-Darwin system [Phys. Plasmas 25, 102506]. The present exposition discusses a generic class of local Vlasov-Maxwell-type systems, with a conscious choice of adopting the language of differential geometry to exploit the Euler-Poincaré framework to its full extent. After reviewing the transition from a Lagrangian picture to an Eulerian one, we demonstrate how symmetries generated by isometries in space lead to conservation laws for linear- and angular-momentum density and how symmetry by time translation produces a conservation law for energy density. We also discuss what happens if no symmetries exist. Finally, two explicit examples will be given -- the classic Vlasov-Maxwell and the drift-kinetic Vlasov-Maxwell -- and the results expressed in the language of regular vector calculus for familiarity. △ Less

Submitted 29 April, 2020; v1 submitted 13 December, 2019; originally announced December 2019.

Comments: 20 pages

Rigidity of MHD equilibria to smooth incompressible ideal motion near resonant surfaces

Authors: David Pfefferlé, Lyle Noakes, Yao Zhou

Abstract: In ideal MHD, the magnetic flux is advected by the plasma motion, freezing flux-surfaces into the flow. An MHD equilibrium is reached when the flow relaxes and force balance is achieved. We ask what classes of MHD equilibria can be accessed from a given initial state via smooth incompressible ideal motion. It is found that certain boundary displacements are formally not supported. This follows fro… ▽ More In ideal MHD, the magnetic flux is advected by the plasma motion, freezing flux-surfaces into the flow. An MHD equilibrium is reached when the flow relaxes and force balance is achieved. We ask what classes of MHD equilibria can be accessed from a given initial state via smooth incompressible ideal motion. It is found that certain boundary displacements are formally not supported. This follows from yet another investigation of the Hahm--Kulsrud--Taylor (HKT) problem, which highlights the resonant behaviour near a rational layer formed by a set of degenerate critical points in the flux-function. When trying to retain the mirror symmetry of the flux-function with respect to the resonant layer, the vector field that generates the volume-preserving diffeomorphism vanishes at the identity to all order in the time-like path parameter. △ Less

Submitted 3 October, 2019; originally announced October 2019.

Comments: 9 pages, 3 figures

Gauge freedom in magnetostatics and the effect on helicity in toroidal volumes

Authors: David Pfefferlé, Lyle Noakes

Abstract: Magnetostatics defines a class of boundary value problems in which the topology of the domain plays a subtle role. For example, representability of a divergence-free field as the curl of a vector potential comes about because of homological considerations. With this in mind, we study gauge-freedom in magnetostatics and its effect on the comparison between magnetic configurations through key quanti… ▽ More Magnetostatics defines a class of boundary value problems in which the topology of the domain plays a subtle role. For example, representability of a divergence-free field as the curl of a vector potential comes about because of homological considerations. With this in mind, we study gauge-freedom in magnetostatics and its effect on the comparison between magnetic configurations through key quantities such as the magnetic helicity. For this, we apply the Hodge decomposition of $k$-forms on compact orientable Riemaniann manifolds with smooth boundary, as well as de Rham cohomology, to the representation of magnetic fields through potential $1$-forms in toroidal volumes. An advantage of the homological approach is the recovery of classical results without explicit coordinates and assumptions about the fields on the exterior of the domain. In particular, a detailed construction of a minimal gauge and a formal proof of relative helicity formulae are presented. △ Less

Submitted 6 September, 2019; originally announced September 2019.

Comments: 15 pages

WWhittle Maximum Likelihood Estimate of spectral properties of Rayleigh-Taylor interfacial mixing using hot-wire anemometry experimental data

Authors: David Pfefferlé, Snezhana I Abarzhi

Abstract: Investigating the power density spectrum of fluctuations in Rayleigh-Taylor (RT) interfacial mixing is a means of studying characteristic length- and time-scales, anisotropies and anomalous processes. Guided by group theory, analysing the invariance-based properties of the fluctuations, our work examines raw time series from hot-wire anemometry measurements in the experiment by Akula et al., JFM 8… ▽ More Investigating the power density spectrum of fluctuations in Rayleigh-Taylor (RT) interfacial mixing is a means of studying characteristic length- and time-scales, anisotropies and anomalous processes. Guided by group theory, analysing the invariance-based properties of the fluctuations, our work examines raw time series from hot-wire anemometry measurements in the experiment by Akula et al., JFM 816, 619-660 (2017). The results suggest that the power density spectrum can be modelled as a compound function presented as the product of a power law and an exponential. The data analysis is based on Whittle's approximation of the power density spectrum for independent zero-mean near-Gaussian signals to construct a Maximum likelihood Estimator (MLE) of the parameters. Those that maximise the log-likelihood are computed numerically through Newton-Raphson iteration. The Hessian of the log-likelihood is used to evaluate the Fisher information matrix and provide an estimate of the statistical error on the obtained parameters. The Kolmogorov-Smirnov test is applied to analyse the goodness-of-fit, by verifying the hypothesis that the ratio between the observed periodogram and the estimated power density spectrum follows a chi-squared probability distribution. The dependence of the parameters of the compound function is investigated on the range of mode numbers over which the fit is performed. In the domain where the relative errors of the power law exponent and the exponential decay rate are small and the goodness-of-fit is excellent, the parameters of the compound function are clearly defined, in agreement with the theory. The study of the power-law spectra in RT mixing data suggests that rigorous physics-based statistical methods can help researchers to see beyond visual inspection. △ Less

Submitted 7 October, 2020; v1 submitted 11 August, 2019; originally announced August 2019.

Comments: 14 pages, 13 figures

Non-planar elasticae as optimal curves for the magnetic axis of stellarators

Authors: David Pfefferlé, Lee Gunderson, Stuart R. Hudson, Lyle Noakes

Abstract: The problem of finding an optimal curve for the target magnetic axis of a stellarator is addressed. Euler-Lagrange equations are derived for finite length three-dimensional curves that extremise their bending energy while yielding fixed integrated torsion. The obvious translational and rotational symmetry is exploited to express solutions in a preferred cylindrical coordinate system in terms of el… ▽ More The problem of finding an optimal curve for the target magnetic axis of a stellarator is addressed. Euler-Lagrange equations are derived for finite length three-dimensional curves that extremise their bending energy while yielding fixed integrated torsion. The obvious translational and rotational symmetry is exploited to express solutions in a preferred cylindrical coordinate system in terms of elliptic Jacobi functions. These solution curves, which, up to similarity transformations, depend on three dimensionless parameters, do not necessarily close. Two closure conditions are obtained for the vertical and toroidal displacement (the radial coordinate being trivially periodic) to yield a countably infinite set of one-parameter families of closed non-planar curves. The behaviour of the integrated torsion (Twist of the Frenet frame), the Linking of the Frenet frame and the Writhe of the solution curves is studied in light of the \Calugareanu theorem. A refreshed interpretation of Mercier's formula for the on-axis rotational transform of stellarator magnetic field-lines is proposed. △ Less

Submitted 6 September, 2018; v1 submitted 6 May, 2018; originally announced May 2018.

Algebraic motion of vertically displacing plasmas

Authors: David Pfefferlé, Amitava Bhattacharjee

Abstract: The vertical motion of a tokamak plasma is analytically modelled during its non-linear phase by a free-moving current-carrying rod inductively coupled to a set of fixed conducting wires or a cylindrical conducting shell. The solutions capture the leading term in a Taylor expansion of the Green's function for the interaction between the plasma column and the surrounding vacuum vessel. The plasma sh… ▽ More The vertical motion of a tokamak plasma is analytically modelled during its non-linear phase by a free-moving current-carrying rod inductively coupled to a set of fixed conducting wires or a cylindrical conducting shell. The solutions capture the leading term in a Taylor expansion of the Green's function for the interaction between the plasma column and the surrounding vacuum vessel. The plasma shape and profiles are assumed not to vary during the vertical drifting phase such that the plasma column behaves as a rigid body. In the limit of perfectly conducting structures, the plasma is prevented to come in contact with the wall due to steep effective potential barriers created by the induced Eddy currents. Resistivity in the wall allows the equilibrium point to drift towards the vessel on the slow timescale of flux penetration. The initial exponential motion of the plasma, understood as a resistive vertical instability, is succeeded by a non-linear "sinking" behaviour shown to be algebraic and decelerating. The acceleration of the plasma column often observed in experiments is thus concluded to originate from an early sharing of toroidal current between the core, the halo plasma and the wall or from the thermal quench dynamics precipitating loss of plasma current. △ Less

Submitted 6 May, 2018; v1 submitted 18 October, 2017; originally announced October 2017.

Journal ref: Physics of Plasmas 25, 022516 (2018)

Differential formulation of the gyrokinetic Landau operator

Authors: Eero Hirvijoki, Alain J. Brizard, David Pfefferlé

Abstract: Subsequent to the recent rigorous derivation of an energetically consistent gyrokinetic collision operator in the so-called Landau representation, this paper investigates the possibility of finding a differential formulation of the gyrokinetic Landau collision operator. It is observed that, while a differential formulation is possible in the gyrokinetic phase space, reduction of the resulting syst… ▽ More Subsequent to the recent rigorous derivation of an energetically consistent gyrokinetic collision operator in the so-called Landau representation, this paper investigates the possibility of finding a differential formulation of the gyrokinetic Landau collision operator. It is observed that, while a differential formulation is possible in the gyrokinetic phase space, reduction of the resulting system of partial differential equations to five dimensions via gyroaveraging poses a challenge. Based on the present work, it is likely that the gyrocentre analogues of the Rosenbluth-MacDonald-Judd potential functions must be kept gyroangle dependent. △ Less

Submitted 18 July, 2017; originally announced July 2017.

Comments: 11 pages

Journal ref: Journal of Plasma Physics, Volume 83, Issue 1, February 2017

Exact collisional moments for plasma fluid theories

Authors: David Pfefferlé, Eero Hirvijoki, Manasvi Lingam

Abstract: The velocity-space moments of the often troublesome nonlinear Landau collision operator are expressed exactly in terms of multi-index Hermite-polynomial moments of the distribution functions. The collisional moments are shown to be generated by derivatives of two well-known functions, namely the Rosenbluth-MacDonald-Judd-Trubnikov potentials for a Gaussian distribution. The resulting formula has a… ▽ More The velocity-space moments of the often troublesome nonlinear Landau collision operator are expressed exactly in terms of multi-index Hermite-polynomial moments of the distribution functions. The collisional moments are shown to be generated by derivatives of two well-known functions, namely the Rosenbluth-MacDonald-Judd-Trubnikov potentials for a Gaussian distribution. The resulting formula has a nonlinear dependency on the relative mean flow of the colliding species normalised to the root-mean-square of the corresponding thermal velocities, and a bilinear dependency on densities and higher-order velocity moments of the distribution functions, with no restriction on temperature, flow or mass ratio of the species. The result can be applied to both the classic transport theory of plasmas, that relies on the Chapman-Enskog method, as well as to deriving collisional fluid equations that follow Grad's moment approach. As an illustrative example, we provide the collisional ten-moment equations with exact conservation laws for momentum- and energy-transfer rate. △ Less

Submitted 14 March, 2017; v1 submitted 27 January, 2017; originally announced January 2017.

Nonlinear Resistivity for Magnetohydrodynamical Models

Authors: Manasvi Lingam, Eero Hirvijoki, David Pfefferlé, Luca Comisso, Amitava Bhattacharjee

Abstract: A new formulation of the plasma resistivity that arises from the collisional momentum-transfer rate between electrons and ions is presented. The resistivity computed herein is shown to depend not only on the temperature and density but also on all other polynomial velocity-space moments of the distribution function, such as the pressure tensor and heat flux vector. The exact expression for the col… ▽ More A new formulation of the plasma resistivity that arises from the collisional momentum-transfer rate between electrons and ions is presented. The resistivity computed herein is shown to depend not only on the temperature and density but also on all other polynomial velocity-space moments of the distribution function, such as the pressure tensor and heat flux vector. The exact expression for the collisional momentum-transfer rate is determined, and is used to formulate the nonlinear anisotropic resistivity. The new formalism recovers the Spitzer resistivity, as well as the concept of thermal force if the heat flux is assumed to be proportional to a temperature gradient. Furthermore, if the pressure tensor is related to viscous stress, the latter enters the expression for the resistivity. The relative importance of the nonlinear term(s) with respect to the well-established electron inertia and Hall terms is also examined. The subtle implications of the nonlinear resistivity, and its dependence on the fluid variables, are discussed in the context of magnetized plasma environments and phenomena such as magnetic reconnection. △ Less

Submitted 15 March, 2017; v1 submitted 17 October, 2016; originally announced October 2016.

Comments: 9 pages; 0 figures

Journal ref: Phys. Plasmas 24, 042120 (2017)

Fluid moments of the Landau collision operator

Authors: Eero Hirvijoki, Manasvi Lingam, David Pfefferlé, Luca Comisso, Jeff Candy, Amitava Bhattacharjee

Abstract: One important problem in plasma physics is the lack of an accurate and complete description of Coulomb collisions in associated fluid models. To shed light on the problem, this Letter introduces an integral identity involving the multi-dimensional Hermite tensor polynomials and presents a method for computing exact expressions for the fluid moments of the nonlinear Landau collision operator. The p… ▽ More One important problem in plasma physics is the lack of an accurate and complete description of Coulomb collisions in associated fluid models. To shed light on the problem, this Letter introduces an integral identity involving the multi-dimensional Hermite tensor polynomials and presents a method for computing exact expressions for the fluid moments of the nonlinear Landau collision operator. The proposed methodology provides a systematic and rigorous means of extending the validity of fluid models that have an underlying inverse-square force particle dynamics to weakly collisional and strong flow regimes. △ Less

Submitted 24 May, 2016; originally announced May 2016.

Comments: 10 pages

Journal ref: Phys. Plasmas 23, 080701 (2016)

Hybrid guiding-centre/full-orbit simulations in non-axisymmetric magnetic geometry exploiting general criterion for guiding-centre accuracy

Authors: David Pfefferlé, Jonathan P. Graves, Wilfred A. Cooper

Abstract: To identify under what conditions guiding-centre or full-orbit tracing should be used, an estimation of the spatial variation of the magnetic field is proposed, not only taking into account gradient and curvature terms but also parallel currents and the local shearing of field-lines. The criterion is derived for general three-dimensional magnetic equilibria including stellarator plasmas. Details a… ▽ More To identify under what conditions guiding-centre or full-orbit tracing should be used, an estimation of the spatial variation of the magnetic field is proposed, not only taking into account gradient and curvature terms but also parallel currents and the local shearing of field-lines. The criterion is derived for general three-dimensional magnetic equilibria including stellarator plasmas. Details are provided on how to implement it in cylindrical coordinates, and in flux coordinates that rely on the geometric toroidal angle. A means of switching between guiding-centre and full-orbit equations at first order in Larmor radius with minimal discrepancy is shown. Techniques are applied to a MAST (Mega Amp Spherical Tokamak) helical core equilibrium in which the inner kinked flux-surfaces are tightly compressed against the outer axisymmetric mantle and where the parallel current peaks at the nearly rational surface. This is put in relation with the simpler situation $\vec{B}(x,y,z) = B_0 [\sin(kx) \vec{e_y} + \cos(kx)\vec{e_z}]$, for which full orbits and lowest order drifts are obtained analytically. In the kinked equilibrium, the full orbits of NBI fast ions are solved numerically and shown to follow helical drift surfaces. This result partially explains the off-axis redistribution of NBI fast particles in the presence of MAST Long-Lived Modes (LLM). △ Less

Submitted 26 February, 2015; v1 submitted 17 December, 2014; originally announced December 2014.

Comments: 12 pages, 12 figures