Eero Hirvijoki

Eero Hirvijoki, Michael Kraus, Joshua W. Burby

In this paper, we present a new framework for addressing the nonlinear Landau collision operator in terms of particle-in-cell methods. We employ the underlying metriplectic structure of the collision operator and, using a macro particle discretization for the distribution function, we transform the infinite-dimensional system into a finite-dimensional time-continuous metriplectic system for advancing the macro particle weights. Temporal discretization is accomplished using the concept of discrete gradients. The conservation of density, momentum, and energy, as well as the positive semi-definite production of entropy in both the time-continuous and the fully discrete system is demonstrated algebraically. The new algorithm is fully compatible with the existing particle-in-cell Poisson integrators for the Vlasov-Maxwell system.

Eero Hirvijoki, Joshua W. Burby, Michael Kraus

This paper explores energy-, momentum-, density-, and positivity-preserving spatio-temporal discretizations for the nonlinear Landau collision operator. We discuss two approaches, namely direct Galerkin formulations and discretizations of the underlying infinite-dimensional metriplectic structure of the collision integral. The spatial discretizations are chosen to reproduce the time-continuous conservation laws that correspond to Casimir invariants and to guarantee the positivity of the distribution function. Both the direct and the metriplectic discretization are demonstrated to have exact H-theorems and unique, physically exact equilibrium states. Most importantly, the two approaches are shown to coincide, given the chosen Galerkin method. A temporal discretization, preserving all of the mentioned properties, is achieved with so-called discrete gradients. Hence the proposed algorithm successfully translates all properties of the infinite-dimensional time-continuous Landau collision operator to time- and space-discrete sparse-matrix equations suitable for numerical simulation.

Michael Kraus, Eero Hirvijoki

We present a novel framework for addressing the nonlinear Landau collision integral in terms of finite element and other subspace projection methods. We employ the underlying metriplectic structure of the Landau collision integral and, using a Galerkin discretization for the velocity space, we transform the infinite-dimensional system into a finite-dimensional, time-continuous metriplectic system. Temporal discretization is accomplished using the concept of discrete gradients. The conservation of energy, momentum, and particle densities, as well as the production of entropy is demonstrated algebraically for the fully discrete system. Due to the generality of our approach, the conservation properties and the monotonic behavior of entropy are guaranteed for finite element discretizations in general, independently of the mesh configuration.