“…Various specially designed numerical schemes with the jump conditions being incorporated into discretization have been introduced in the literature for solving parabolic interface problems. These include finite element methods and finite volume methods based on body-fitted interface treatments [3,8,32,33,34], finite difference methods based on Cartesian grids [2,5,6,18,19,22,17,30,40,24,36,20], and immersed finite element methods based on Cartesian meshes [25,39]. For instance, the immersed interface method (IIM) is one of the most successful finite difference methods in solving parabolic interface problems [5,6,18,19,20].…”