 Error-function isobaric cooling front

For conductivity $\lambda ={\lambda }_{0}T$, no cooling in the hot phase and a zero-temperature boundary at $m=0$, the mass-coordinate equation is linear diffusion. Its similarity profile solves $f^{\prime \prime }+\xi f^{\prime }/2=0$, $f(0)=0$, $f(\infty )=1$, giving the displayed error function. The thermal layer broadens as the square root of time. Conducted energy enters the cold radiating reservoir at rate per area $L={\lambda }_{0}A{T}_m(0)=A{T}_0\sqrt{{\lambda }_0}/\sqrt{\pi \tau }$, hence $L\propto t^{-1/2}$.