Heat flux manipulation has garnered extensive research attention due to its potential application in thermal management devices, such as heat sink, [1,2] thermal diode, [3-6] thermal transistor, [7] and thermal memory. [8] Yet until now, we were not able to control heat flux as we do with electric current. Considerable effort of researchers to investigate possible transport mechanism of microscopic particles, such as photons [9] and phonons, [10-12] has led to some significant conclusions. For instance, in 2006, Leonhardt and Pendry. [13,14] established the transformation optics theory and proposed the concept of invisible cloak, which have also been verified by subsequent experiments. Inspired by this, in 2008, Fan et al. [15] first applied the transformation theory in the field of heat transfer, and theoretically predicted the mode of thermal invisible cloak based on graded thermal metamaterial, in which the thermal conduction equation remains form invariant under coordinate transformation. However, the adverse establishing conditions of transformation thermotics theory limit the practical implementation of thermal invisible cloak. Until 2012, Narayana and Sato [16] successfully fabricated a thermal invisible cloak by constructing a concentric layered structure consisting of alternating layers of two different homogeneous isotropic thermal conductivities. Since then, a variety of thermal metamaterials, including thermal cloak and thermal concentrator, have been proposed, and some of them have been experimentally realized, [17-25] paving the way for follow-up research. It is worth noting that these novel thermal phenomena in above-mentioned metamaterials are mostly based on macroscopic compound structures such as metal/polymer. Since 2004, significant research has been directed toward two-dimensional (2D) materials. [26-31] Graphene, as a representative 2D material, is considered as a promising candidate for next-generation electronic and optoelectronic devices due to its exceptional charge carrier mobility and 2D nature. [32-36] In addition, ultrahigh phonon-dominated thermal conductivity of graphene facilitates its use as heat sinks in integrated circuits. [1,2,37] Moreover, the fast response speed in graphene induced by large phonon group velocity can be applied in thermal devices performed with phonons. For example, Wang et al. studied the thermal rectification (TR) effect in asymmetrically defected graphene nanoribbon (GNR) and showed that singlevacancies and Si defects are more preferable in thermal rectifier fabrication. [38] A similar result has been obtained by experimental studies on asymmetrically defected GNR. [39] Furthermore, Pal and Puri realized the function of thermal logic gate based on two asymmetric graphene thermal diodes.