Supplementary Variable Method for Developing Structure-Preserving Numerical Approximations to Thermodynamically Consistent Partial Differential Equations.