Search

The triangular mesh-based gyrokinetic scheme enables comprehensive axis-to-edge studies across the entire plasma volume. Our approach employs triangular finite elements with first-derivative continuity (C1), building on previous work to facilitate gyrokinetic simulations. Additionally, we have adopted the mixed variable/pullback scheme for gyrokinetic electromagnetic particle simulations. The filter-free treatment in the poloidal cross-section with triangular meshes introduces unique features and challenges compared to previous treatments using structured meshes. Our implementation has been validated through benchmarks using ITPA-TAE (Toroidicity-induced Alfv\'en Eigenmode) parameters, showing its capability in moderate to small electron skin depth regimes. Additional examinations using experimental parameters confirm its applicability to realistic plasma conditions., Comment: 31 pages, 9 figures. Joint Varenna-Lausanne International Workshop on Theory of Fusion Plasmas. Submitted to Plasma Physics and Controlled Fusion (2024)

A simple Hamiltonian modeling framework for general models in nonlinear optics is given. This framework is specialized to describe the Hamiltonian structure of electromagnetic phenomena in cubicly nonlinear optical media. The model has a simple Poisson bracket structure with the Hamiltonian encoding all of the nonlinear coupling of the fields. The field-independence of the Poisson bracket facilitates a straightforward Hamiltonian structure-preserving discretization using finite element exterior calculus. The generality and relative simplicity of this Hamiltonian framework makes it amenable for simulating a broad class of time-domain nonlinear optical problems. The main contribution of this work is a finite element discretization of Maxwell's equations in cubicly nonlinear media which is energy-stable and exactly conserves Gauss's laws. Moreover, this approach may be readily adapted to consider more general nonlinear media in subsequent work.

In this manuscript we present a novel and efficient numerical method for the compressible viscous and resistive MHD equations for all Mach number regimes. The time-integration strategy is a semi-implicit splitting, combined with a hybrid finite-volume and finite-element (FE) discretization in space. The non-linear convection is solved by a robust explicit FV scheme, while the magneto-acoustic terms are treated implicitly in time. The resulting CFL stability condition depends only on the fluid velocity, and not on the Alfv\'enic and acoustic modes. The magneto-acoustic terms are discretized by compatible FE based on a continuous and a discrete de Rham complexes designed using Finite Element Exterior Calculus (FEEC). Thanks to the use of FEEC, energy stability, magnetic-helicity conservation and the divergence-free conditions can be preserved also at the discrete level. A very efficient splitting approach is used to separate the acoustic and the Alfv\'enic modes in such a fashion that the original symmetries of the PDE governing equations are preserved. In this way, the algorithm relies on the solution of linear, symmetric and positive-definite algebraic systems, that are very efficiently handled by the simple matrix-free conjugate-gradient method. The resulting algorithm showed to be robust and accurate in low and high Mach regimes even at large Courant numbers. Non-trivial tests are solved in one-, two- and three- space dimensions to confirm the robustness, accuracy, and the low-dissipative and conserving properties of the final algorithm. While the formulation of the method is very general, numerical results for a second-order accurate FV-FE scheme will be presented., Comment: 46 pages, 18 figures, 3 tables

The use of gyrokinetics, wherein phase-space coordinate transformations result in a phase-space dimensionality reduction as well as removal of fast time scales, has enabled the simulation of microturbulence in fusion devices. The state-of-the-art gyrokinetic models used in practice are parallel-only models wherein the perpendicular part of the vector potential is neglected. Such models are inherently not gauge invariant. We generalize the work of [Burby, Brizard. Physics Letters A, 383(18):2172-2175] by deriving a sufficient condition on the gyrocentre coordinate transformation which ensures gauge invariance. This leads to a parametrized family of gyrokinetic models for which we motivate a specific choice of parameters that results in a physically meaningful gauge invariant model. Due to gauge invariance this model can be expressed directly in terms of the electromagnetic fields, rather than the potentials, and the gyrokinetic model thereby results in the macroscopic Maxwell's equations. For the linearized model, it is demonstrated that the shear and compressional Alfv\'en waves are present with the correct frequencies. The fast compressional Alfv\'en wave can be removed by making use of a Darwin-like approximation. This approximation breaks gauge invariance, but we find that the field equations are still compatible without the need of a Lagrange multiplier.

In the presence of an inhomogeneous oscillatory electric field, charged particles experience a net force, averaged over the oscillatory timescale, known as the ponderomotive force. We derive a one-dimensional Hamiltonian model which self-consistently couples the electromagnetic field to a plasma which experiences the ponderomotive force. We derive a family of structure preserving discretizations of the model of varying order in space and time using conforming and broken finite element exterior calculus spectral element methods. In all variants of our discretization framework, the method is found to conserve the Casimir invariants of the continuous model to machine precision and the energy to the order of the splitting method used.

The gyro-kinetic model is an approximation of the Vlasov-Maxwell system in a strongly magnetized magnetic field. We propose a new algorithm for solving it combining the Semi-Lagrangian (SL) method and the Arakawa (AKW) scheme with a time-integrator. Both methods are successfully used in practice for different kinds of applications, in our case, we combine them by first decomposing the problem into a fast (parallel) and a slow (perpendicular) dynamical system. The SL approach and the AKW scheme will be used to solve respectively the fast and the slow subsystems. Compared to the scheme in [1], where the entire model is solved using only the SL method, our goal is to replace the method used in the slow subsystem by the AKW scheme, in order to improve the conservation of the physical constants.

We explore the possibilities of applying structure-preserving numerical methods to a plasma hybrid model with kinetic ions and mass-less fluid electrons satisfying the quasi-neutrality relation. The numerical schemes are derived by finite element methods in the framework of finite element exterior calculus (FEEC) for field variables, particle-in-cell (PIC) methods for the Vlasov equation, and splitting methods in time based on an anti-symmetric bracket proposed. Conservation properties of energy, quasi-neutrality relation, positivity of density, and divergence-free property of the magnetic field are given irrespective of the used resolution and metric. Local quasi-interpolation is used for dealing with the current terms in order to make the proposed methods more efficient. The implementation has been done in the framework of the Python package STRUPHY [1], and has been verified by extensive numerical experiments., Comment: 30 pages, 8 figures

We study the canonical variables based numerical schemes of a hybrid model with kinetic ions and mass-less electrons. Two equivalent formulations of the hybrid model are presented with the vector potentials in different gauges and the distribution functions depending on canonical momentum (not velocity), which constitutes a pair of canonical variables with the position variable. Particle-in-cell methods are used for the distribution functions, and the vector potentials are discretized by the finite element methods in the framework of finite element exterior calculus. Splitting methods are used for the time discretizations. It is illustrated that the second formulation is numerically superior and the schemes constructed based on the anti-symmetric bracket proposed have better conservation properties and lower noise, although the filters can be used to improve the schemes of the first formulation., Comment: 25 pages, 8 figures

Three geometric formulations of the Hamiltonian structure of the macroscopic Maxwell equations are given: one in terms of the double de Rham complex, one in terms of L2 duality, and one utilizing an abstract notion of duality. The final of these is used to express the geometric and Hamiltonian structure of kinetic theories in general media. The Poisson bracket so stated is explicitly metric free. Finally, as a special case, the Lorentz covariance of such kinetic theories is investigated. We obtain a Lorentz covariant kinetic theory coupled to nonlinear electrodynamics such as Born-Infeld or Euler-Heisenberg electrodynamics.

Structure-preserving methods can be derived for the Vlasov-Maxwell system from a discretisation of the Poisson bracket with compatible finite-elements for the fields and a particle representation of the distribution function. These geometric electromagnetic particle-in-cell (GEMPIC) discretisations feature excellent conservation properties and long-time numerical stability. This paper extends the GEMPIC formulation in curvilinear coordinates to realistic boundary conditions. We build a de Rham sequence based on spline functions with clamped boundaries and apply perfect conductor boundary conditions for the fields and reflecting boundary conditions for the particles. The spatial semi-discretisation forms a discrete Poisson system. Time discretisation is either done by Hamiltonian splitting yielding a semi-explicit Gauss conserving scheme or by a discrete gradient scheme applied to a Poisson splitting yielding a semi-implicit energy-conserving scheme. Our system requires the inversion of the spline finite element mass matrices, which we precondition with the combination of a Jacobi preconditioner and the spectrum of the mass matrices on a periodic tensor product grid.

A mimetic spectral element discretization, utilizing a novel Galerkin projection Hodge star operator, of the macroscopic Maxwell equations in Hamiltonian form is presented. The idea of splitting purely topological and metric dependent quantities is natural in the Hamiltonian modeling framework as the Poisson bracket is metric free with the Hamiltonian containing all metric information. This idea may be incorporated into the mimetic spectral element method by directly discretizing the Poincar\'e duality structure. This "split exterior calculus mimetic spectral element method" yields spatially discretized Maxwell's equations which are Hamiltonian and exactly and strongly conserve Gauss's laws. Moreover, the new discrete Hodge star operator is itself of interest as a partition of the purely topological and metric dependent portions of the Hodge star operator. As a simple test case, the numerical results of applying this method to a one-dimensional version of Maxwell's equations are given., Comment: 25 pages, 4 figures

A Hamiltonian field theory for the macroscopic Maxwell equations with fully general polarization and magnetization is stated in the language of differential forms. The precise procedure for translating the vector calculus formulation into differential forms is discussed in detail. We choose to distinguish between straight and twisted differential forms so that all integrals be taken over densities (i.e. twisted top forms). This ensures that the duality pairings, which are stated as integrals over densities, are orientation independent. The relationship between functional differentiation with respect to vector fields and with respect to differential forms is established using the chain rule. The theory is developed such that the Poisson bracket is metric and orientation independent with all metric dependence contained in the Hamiltonian. As is typically seen in the exterior calculus formulation of Maxwell's equations, the Hodge star operator plays a key role in modeling the constitutive relations. As a demonstration of the kind of constitutive models this theory accommodates, the paper concludes with several examples., Comment: 24 pages

In this article we describe a unifying framework for variational electromagnetic particle schemes of spectral type, and we propose a novel spectral Particle-In-Cell (PIC) scheme that preserves a discrete Hamiltonian structure. Our work is based on a new abstract variational derivation of particle schemes which builds on a de Rham complex where Low's Lagrangian is discretized using a particle approximation of the distribution function. In this framework, which extends the recent Finite Element based Geometric Electromagnetic PIC (GEMPIC) method to a variety of field solvers, the discretization of the electromagnetic potentials and fields is represented by a de Rham sequence of compatible spaces, and the particle-field coupling procedure is described by approximation operators that commute with the differential operators in the sequence. In particular, for spectral Maxwell solvers the choice of truncated $L^2$ projections using continuous Fourier transform coefficients for the commuting approximation operators yields the gridless Particle-in-Fourier method, whereas spectral Particle-in-Cell methods are obtained by using discrete Fourier transform coefficients computed from a grid. By introducing a new sequence of spectral pseudo-differential approximation operators, we then obtain a novel variational spectral PIC method with discrete Hamiltonian structure that we call Fourier-GEMPIC. Fully discrete schemes are then derived using a Hamiltonian splitting procedure, leading to explicit time steps that preserve the Gauss laws and the discrete Poisson bracket associated with the Hamiltonian structure. These explicit steps share many similarities with standard spectral PIC methods. As arbitrary filters are allowed in our framework, we also discuss aliasing errors and study a natural back-filtering procedure to mitigate the damping caused by anti-aliasing smoothing particle shapes., Comment: 40 pages, 34 figures

In this article we apply a discrete action principle for the Vlasov--Maxwell equations in a structure-preserving particle-field discretization framework. In this framework the finite-dimensional electromagnetic potentials and fields are represented in a discrete de Rham sequence involving general finite element spaces, and the particle-field coupling is represented by a set of projection operators that commute with the differential operators. With a minimal number of assumptions which allow for a variety of finite elements and shape functions for the particles, we show that the resulting variational scheme has a general discrete Poisson structure and thus leads to a semi-discrete Hamiltonian system. By introducing discrete interior products we derive a second type of space discretization which is momentum preserving, based on the same finite elements and shape functions. We illustrate our method by applying it to spline finite elements, and to a new spectral discretization where the particle-field coupling relies on discrete Fourier transforms., Comment: 37 pages, 18 figures

Numerical schemes that preserve the structure of the kinetic equations can provide stable simulation results over a long time. An electromagnetic particle-in-cell solver for the Vlasov-Maxwell equations that preserves at the discrete level the non-canonical Hamiltonian structure of the Vlasov-Maxwell equations has been presented in [Kraus et al. 2017]. Whereas the original formulation has been obtained for Cartesian coordinates, we extend the formulation to curvilinear coordinates in this paper. For the discretisation in time, we discuss several (semi-)implicit methods either based on a Hamiltonian splitting or a discrete gradient method combined with an antisymmetric splitting of the Poisson matrix and discuss their conservation properties and computational efficiency.

This paper discusses energy-conserving time-discretizations for finite element particle-in-cell discretizations of the Vlasov--Maxwell system. A geometric spatially discrete system can be obtained using a standard particle-in-cell discretization of the particle distribution and compatible finite element spaces for the fields to discretize the Poisson bracket of the Vlasov--Maxwell model (see Kraus et al., J Plasma Phys 83, 2017). In this paper, we derive energy-conserving time-discretizations based on the discrete gradient method applied to an antisymmetric splitting of the Poisson matrix. Firstly, we propose a semi-implicit method based on the average-vector-field discretization of the subsystems. Moreover, we devise an alternative discrete gradient that yields a time discretization that can additionally conserve Gauss' law. Finally, we explain how substepping for fast species dynamics can be incorporated.

A two-fold analysis of electromagnetic core tokamak instabilities in the framework of the gyrokinetic theory is presented. First principle theoretical foundations of the gyrokinetic theory are used to explain and justify the numerical results obtained with the global electromagnetic particle-in-cell code Orb5 whose model is derived from the Lagrangian formalism. The energy conservation law corresponding to the Orb5 model is derived from the Noether theorem and implemented in the code as a diagnostics for energy balance and conservation verification. An additional Noether theorem based diagnostics is implemented in order to analyse destabilising mechanisms for the electrostatic and the electromagnetic Ion Temperature Gradient (ITG) instabilities in the core region of the tokamak. The transition towards the Kinetic Ballooning Modes (KBM) at high electromagnetic $\beta$ is also investigated., Comment: 22 pages, 10 Figures, material form the ICPP conference 2018, invited

We present a novel framework for Finite Element Particle-in-Cell methods based on the discretization of the underlying Hamiltonian structure of the Vlasov-Maxwell system. We derive a semi-discrete Poisson bracket, which retains the defining properties of a bracket, anti-symmetry and the Jacobi identity, as well as conservation of its Casimir invariants, implying that the semi-discrete system is still a Hamiltonian system. In order to obtain a fully discrete Poisson integrator, the semi-discrete bracket is used in conjunction with Hamiltonian splitting methods for integration in time. Techniques from Finite Element Exterior Calculus ensure conservation of the divergence of the magnetic field and Gauss' law as well as stability of the field solver. The resulting methods are gauge invariant, feature exact charge conservation and show excellent long-time energy and momentum behaviour. Due to the generality of our framework, these conservation properties are guaranteed independently of a particular choice of the Finite Element basis, as long as the corresponding Finite Element spaces satisfy certain compatibility conditions., Comment: 57 Pages

The main idea of Gyrokinetic dynamical reduction consists in systematical removing of fastest scale of motion (the gyro motion) from plasma's dynamics, resulting in a considerable model simplification and gain of computing time. Gyrokinetic Maxwell-Vlasov system is broadly implemented in nowadays numerical experiments for modeling strongly magnetized plasma (both laboratory and astrophysical). Different versions of reduced set of equations exist depending on the construction of the Gyrokinetic reduction procedure and approximations assumed while their derivation. The purpose of this paper is to explicitly show the connection between the general second order gyrokinetic Maxwell-Vlasov system issued from the Modern Gyrokinetic theory derivation and the model currently implemented in global electromagnetic Particle in Cell code ORB5. Strictly necessary information about the Modern Gyrokinetic formalism is given together with the consistent derivation of the gyrokinetic Maxwell-Vlasov equations from the first principle of dynamics. Variational formulation of dynamics is also used for simultaneously obtaining the corresponding exact energy conservation law. The result of that explicit derivation is used for verification of the energy conservation diagnostics currently implemented in ORB5 code. This work subscribes into the context of codes verification project VeriGyro currently run in IPP Max Planck Institut in collaboration with others European institutions., Comment: 53 pages, 6 figures