We model the latent heats of crystallization and fusion in phase change materials with a unified latent heat of phase change, ensuring energy conservation by coupling the heat of phase change with amorphous and crystalline specific heats. We demonstrate the model with 2-D finite element simulations of Ge2Sb2Te5 and find that the heat of phase change increases local temperature up to 180 K in 300 nm × 300 nm structures during crystallization, significantly impacting grain distributions. We also show in electrothermal simulations of 45 nm confined and 10 nm mushroom cells that the higher amorphous specific heat predicted by this model increases nucleation probability at the end of reset operations. These nuclei can decrease set time, leading to variability, as demonstrated for the mushroom cell.