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.

Energy-consistent Petrov-Galerkin time discretization of port-Hamiltonian systems

Energy-Preserving and Passivity-Consistent Numerical Discretization of Port-Hamiltonian Systems

On a Generalized Energy Conservation/dissipation Time Finite Element Method for Hamiltonian Mechanics

Goal-oriented time adaptivity for port-Hamiltonian systems

Structure-Preserving Time Discretization of Port-Hamiltonian Systems via Discrete Gradient Pairs

Port-Hamiltonian discontinuous Galerkin finite element methods

Discrete-time port-Hamiltonian systems: A definition based on symplectic integration

General Local Energy-Preserving Integrators for Solving Multi-Symplectic Hamiltonian PDEs.

On Multiscale Methods in Petrov-Galerkin formulation

Efficient and Parallel Solution of High-order Continuous Time Galerkin for Dissipative and Wave Propagation Problems

Discrete nonlinear elastodynamics in a port-Hamiltonian framework

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

A partitioned finite element method for power-preserving discretization of open systems of conservation laws

Stable and efficient Petrov-Galerkin methods for a kinetic Fokker-Planck equation

Port-Hamiltonian systems with energy and power ports

A geometric framework for discrete time port-Hamiltonian systems

Efficient High-Order Space-Angle-Energy Polytopic Discontinuous Galerkin Finite Element Methods for Linear Boltzmann Transport

Port-Hamiltonian formulation and structure-preserving discretization of hyperelastic strings

Hamiltonian-Preserving Discontinuous Galerkin Methods for the Liouville Equation with Discontinuous Potential

The design of conservative finite element discretisations for the vectorial modified KdV equation

Convergence of a continuous Galerkin method for hyperbolic-parabolic systems