Discontinuous Galerkin schemes for a class of Hamiltonian evolution equations with applications to plasma fluid and kinetic problems

A. Hakim,G. Hammett,E. Shi,N. Mandell
DOI: https://doi.org/10.48550/arXiv.1908.01814
2019-08-06
Abstract:In this paper we present energy-conserving, mixed discontinuous Galerkin (DG) and continuous Galerkin (CG) schemes for the solution of a broad class of physical systems described by Hamiltonian evolution equations. These systems often arise in fluid mechanics (incompressible Euler equations) and plasma physics (Vlasov--Poisson equations and gyrokinetic equations), for example. The dynamics is described by a distribution function that evolves given a Hamiltonian and a corresponding Poisson bracket operator, with the Hamiltonian itself computed from field equations. Hamiltonian systems have several conserved quantities, including the quadratic invariants of total energy and the $L_2$ norm of the distribution function. For accurate simulations one must ensure that these quadratic invariants are conserved by the discrete scheme. We show that using a discontinuous Galerkin scheme to evolve the distribution function and ensuring that the Hamiltonian lies in its continuous subspace leads to an energy-conserving scheme in the continuous-time limit. Further, the $L_2$ norm is conserved if central fluxes are used to update the distribution function, but decays monotonically when using upwind fluxes. The conservation of density and $L_2$ norm is then used to show that the entropy is a non-decreasing function of time. The proofs shown here apply to any Hamiltonian system, including ones in which the Poisson bracket operator is non-canonical (for example, the gyrokinetic equations). We demonstrate the ability of the scheme to solve the Vlasov--Poisson and incompressible Euler equations in 2D and provide references where we have applied these schemes to solve the much more complex 5D electrostatic and electromagnetic gyrokinetic equations.
Computational Physics,Plasma Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a numerical method that can accurately preserve energy conservation for solving a class of physical systems described by Hamiltonian evolution equations. These systems usually appear in fluid mechanics (such as the incompressible Euler equations) and plasma physics (such as the Vlasov - Poisson equations and gyrokinetic equations). Specifically, the paper proposes a hybrid discontinuous Galerkin (DG) and continuous Galerkin (CG) scheme, aiming to solve the following problems: 1. **Energy Conservation**: For Hamiltonian systems, energy is an important conserved quantity. The method proposed in the paper achieves energy conservation by ensuring that the Hamiltonian lies in the continuous subspace in the time - continuous limit. 2. **L2 - norm Conservation**: The L2 - norm of the distribution function is also an important conserved quantity. The paper proves that when the central flux is used to update the distribution function, the L2 - norm can be accurately conserved; while when the upwind flux is used, the L2 - norm will monotonically decay. 3. **Density Conservation**: The paper proves that the proposed scheme can ensure the conservation of the total number of particles (i.e., the integral of the distribution function). 4. **Entropy Growth**: If the discrete distribution function remains positive - definite, the discrete scheme can make the discrete entropy grow monotonically. 5. **Momentum Conservation**: Although the main objective of the paper is energy conservation, the author also discusses the problem of momentum conservation and points out that the proposed scheme may not be able to accurately preserve momentum conservation in some cases, but the error converges rapidly as the spatial grid is refined. In summary, the main objective of this paper is to develop an efficient numerical method that can accurately preserve the important conservation properties of Hamiltonian systems in numerical simulations, especially for complex problems in plasma physics, such as the Vlasov - Poisson equations and the incompressible Euler equations.