An adaptive refinement scheme is presented to reduce the geometry discretization error and provide higherorder enrichment functions for the interface-enriched generalized FEM. The proposed method relies on the h-adaptive and p-adaptive refinement techniques to reduce the discrepancy between the exact and discretized geometries of curved material interfaces. A thorough discussion is provided on identifying the appropriate level of the refinement for curved interfaces based on the size of the elements of the background mesh. Varied techniques are then studied for selecting the quasi-optimal location of interface nodes to obtain a more accurate approximation of the interface geometry. We also discuss different approaches for creating the integration sub-elements and evaluating the corresponding enrichment functions together with their impact on the performance and computational cost of higher-order enrichments. Several examples are presented to demonstrate the application of the adaptive interface-enriched generalized FEM for modeling thermomechanical problems with intricate geometries. The accuracy and convergence rate of the method are also studied in these example problems.ADAPTIVE INTERFACE-ENRICHED GENERALIZED FEM 1353 approximation from the finite element mesh. Meshfree methods tackle this issue by either limiting the discretization to the domain boundaries [12,13] or replacing the surface/volume mesh with a domain (cloud) of influence [14,15]. In the realm of mesh-independent FEM (methods that eliminate the requirement of using conforming meshes), the generalized/extended FEM (GFEM/XFEM) [16][17][18][19][20] is one of the most successful techniques. This method employs a partition of unity combined with appropriate enrichment functions to capture the weak and/or strong discontinuities in a domain discretized with non-conforming FE meshes [21].Recently, Soghrati et al. [22,23] have developed an interface-enriched generalized FEM (IGFEM) for the mesh-independent treatment of interface problems. In this method, the nonconforming elements intersected by the material/phase interface are divided into smaller integration sub-elements, which also serve as means for constructing enrichment functions. Unlike the conventional GFEM/XFEM, the method relies on enriching the nodes created at the intersection of the interface with the non-conforming elements edges and does not use a partition of unity. The IGFEM yields the same precision and convergence rate as the standard FEM with conforming meshes, while providing advantages such as a lower computational cost and straightforward imposing of Dirichlet boundary conditions compared with the GFEM/XFEM [22]. In [24,25], the IGFEM was applied to the simulation of actively cooled microvascular materials and the multi-scale failure of heterogeneous adhesives. A detailed comparison of the IGFEM with standard FEM and its advantages over GFEM/XFEM are presented in [22,26].In this manuscript, we study the application of h-adaptive and p-adaptive refinement techniques for discretizing...