Invariant-Domain Preserving implicit-explicit finite element methods for Gray Radiation Hydrodynamics