This paper investigates the effects of f (R) = R + α R 2 model on the dissipative anisotropic non-static rotating axial stellar system using structure scalars. By orthogonal splitting of the Riemann tensor, we obtain structure scalars and explore their roles for dissipative relativistic fluid configuration. We investigate the contribution of structure scalars in the evolution of various kinematical variables (e.g. expansion, shear and vorticity). Finally, we discuss thermodynamical aspects and the effect of dissipation on inertial mass density through transport equation. It is interesting to mention here that all our results reduce to the previous known results in the limiting case.