“…There is extensive existing work addressing numerical approximation of PDE posed on composite domains with interfaces or irregular domains, for example, the boundary integral method [11,56], difference potentials method [3,6,26,27,58,67], immersed boundary method [30,42,61,74], immersed interface method [2,46,47,49,69], ghost fluid method [31,32,50,51], the matched interface and boundary method [82,84,85,86], Cartesian grid embedded boundary method [19,41,57,83], multigrid method for elliptic problems with discontinuous coefficients on an arbitrary interface [18], virtual node method [9,39], Voronoi interface method [35,36], the finite difference method [8,10,24,75,78,79] and finite volume method [22,34] based on mapped grids, or cut finite element method [13,14,15,…”