High-order accurate methods based on difference potentials for 2D parabolic interface models

Highly-accurate numerical methods that can efficiently handle problems with interfaces and/or problems in domains with complex geometry are essential for the resolution of a wide range of temporal and spatial scales in many partial differential equations based models from Biology, Materials Science and Physics. In this paper we continue our work started in 1D, and we develop high-order accurate methods based on the Difference Potentials for 2D parabolic interface/composite domain problems. Extensive numerical experiments are provided to illustrate high-order accuracy and efficiency of the developed schemes.