The effect of latent heat on numerical solutions of the heat flow equation

The stability and accuracy of finite difference solutions of the heat flow equation which involve latent heat evolution are discussed. A comparison of true and finite difference solutions is used to show how a numerical solution must be interpreted, and a dimensionless group is developed which governs the appearance of inaccuracies peculiar to numerical solutions involving a latent heat term.