We consider extension of some established techniques of study of tensor fields on Lorentzian manifolds of arbitrary dimension to non-Abelian gauge covariant fields. These are then applied to study of gauge fields with vanishing scalar curvature invariants and universal Yang-Mills fields, for which all possible higher-order corrections to the Yang-Mills equation identically vanish.