Variational modeling of parabolic free boundary problems by