“…In the past, many heat conduction problems in layered composites have been routinely solved by numerical methods such as the finite difference method [4][5][6], finite element method [7,8], finite volume method [9][10][11] and the boundary element method [12,13]. In contrast to mesh-based numerical methods, the meshless methods, which do not need the mesh generation and boundary integral have been proposed such as the analytical method [14], method of fundamental solutions (MFS) [15][16][17], boundary knot method (BKM) [18], collocation Trefftz method (CTM) [19][20][21], radial basis function collocation method (RBFCM) [22][23][24][25], element-free Galerkin method (EFG) [26], reproducing kernel particle method (RKPM) [27,28], modified polynomial expansion method [29], meshless local boundary integral equation method (LBIE) [30,31], and so on. Among these, boundary-type meshless methods have attracted considerable attention because of their simplicity.…”