Higher-order field theory with constraints

A differential-geometric setting for the dynamics of a higher-order field theory is proposed, based on the Skinner and Rusk formalism for mechanics. This approach incorporates aspects of both, the Lagrangian and the Hamiltonian description, since the field equations are formulated using the Lagrangian on a higher-order jet bundle and the canonical multisymplectic form on its affine dual. The result is that we obtain a unique and global intrinsic description of the dynamics. The case of vakonomic constraints is also studied within this formalism.