A discrete variational scheme for isentropic processes in polyconvex thermoelasticity