 Neumann solution of the Stefan problem

The Neumann solution places a planar phase boundary at $h(t)=2\lambda \sqrt{\kappa t}$. Similarity reduction of the heat equation gives temperature profiles in terms of the error function, while the Stefan condition determines $\lambda$ implicitly.