A polyhedral discrete de Rham numerical scheme for the Yang-Mills equations

Abstract

We present a discretisation of the 3+1 formulation of the Yang-Mills equations in the temporal gauge, using a Lie algebra-valued extension of the discrete de Rham (DDR) sequence, that preserves the non-linear constraint exactly. In contrast to Maxwell's equations, where the preservation of the analogous constraint only depends on reproducing some complex properties of the continuous de Rham sequence, the preservation of the non-linear constraint relies for the Yang-Mills equations on a constrained formulation, previously proposed in [10]. The fully discrete nature of the DDR method requires to devise appropriate constructions of the non-linear terms, adapted to the discrete spaces and to the need for replicating the crucial Ad-invariance property of the $L^2$-product. We then prove some energy estimates, and provide results of 3D numerical simulations based on this scheme.

References in corpus (6)

• Arbitrary-order pressure-robust DDR and VEM methods for the Stokes problem on polyhedral meshes
• Fully discrete polynomial de Rham sequences of arbitrary degree on polygons and polyhedra
• An arbitrary-order method for magnetostatics on polyhedral meshes based on a discrete de Rham sequence
• Algebraic stability analysis of constraint propagation
• An arbitrary-order discrete de Rham complex on polyhedral meshes: Exactness, Poincaré inequalities, and consistency
• Charge-conserving hybrid methods for the Yang-Mills equations

Cited by in corpus (3)

• An exterior calculus framework for polytopal methods
• A polytopal discrete de Rham scheme for the exterior calculus Einstein's equations
• 2D discrete Yang-Mills equations on the torus