In this study, two-dimensional thermal cloak for circular domain is discussed numerically by using transformed steady heat conduction equation and finite-element analysis. It was found that both of ideal cloak and Han's cloak work well for cloaking a circular domain within a rectangular conductive domain having two thermal boundary conditions.