In　this　article,　we　propose　a　novel　meshless　method　for　solving　two-dimensional　stationary　heat　conduction　problems　in　layered　materials.　The　proposed　method　is　a　recently　developed　boundary-type　meshless　method　which　combines　the　collocation　scheme　from　the　method　of　fundamental　solutions　(MFS)　with　the　collocation　Trefftz　method　(CTM)　to　improve　the　applicability　of　the　method　for　solving　boundary　value　problems.　Particular　non-singular　basis　functions　from　cylindrical　harmonics　are　adopted　in　which　the　numerical　approximation　is　based　on　the　superposition　principle　using　the　non-singular　basis　functions　expressed　in　terms　of　many　source　points.　For　the　modeling　of　multi-layer　composite　materials,　we　adopted　the　domain　decomposition　method　(DDM),　which　splits　the　domain　into　smaller　subdomains.　The　continuity　of　the　flux　and　the　temperature　has　to　be　satisfied　at　the　interface　of　subdomains　for　the　problem.　The　validity　of　the　proposed　method　is　investigated　for　several　test　problems.　Numerical　applications　were　also　carried　out.　Comparison　of　the　proposed　method　with　other　meshless　methods　showed　that　it　is　highly　accurate　and　computationally　efficient　for　modeling　heat　conduction　problems,　especially　in　heterogeneous　multi-layer　composite　materials.