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

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 nonlinear 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 nonlinear 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 nonlinear 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.