A Stable FE-FD Method for Anisotropic Parabolic PDEs with Moving Interfaces