The present manuscript gives an overview of how two-dimensional heat diffusion models underwent a fractional transformation, system coupling as well as solution treatment. The governing diffusion models, which are endowed with Caputo's fractional-order derivatives in time $t$, are suitably coupled using the (1) convection phenomenon, (2) interfacial coupling by considering the mechanism of a double-layered bar, and the (3) nonlinear coupling due to temperature-dependent thermal diffusivities. Semi-analytical and analytical methods are considered for the solution treatment. Moreover, we seek a computational environment to graphically illustrate the systems' response to different fractional orders in each case through the determined diffusional fields. Besides, we supply certain concluding notes at the end.