An unconditionally stable fast high order method for thermal phase change models. (15th April 2022)