Variational approximation for a non-isothermal coupled phase-field system: Structure-preservation & Nonlinear stability.