Neumann solution of the Stefan problem (source code)

= Neumann solution of the Stefan problem
{c}

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.