  This paper systematically explains how to apply the invariant subspace method
using variable transformation for finding the exact solutions of the
(k+1)-dimensional nonlinear time-fractional PDEs in detail. More precisely, we
have shown how to transform the given (k+1)-dimensional nonlinear
time-fractional PDEs into (1+1)-dimensional nonlinear time-fractional PDEs
using the variable transformation procedure. Also, we explain how to derive the
exact solutions for the reduced equations using the invariant subspace method.
Additionally, in this careful and systematic study, we will investigate how to
find the various types of exact solutions of the (3+1)-dimensional nonlinear
time-fractional convection-diffusion-reaction equation along with appropriate
initial and boundary conditions for the first time. Moreover, the obtained
exact solutions of the equation as mentioned above can be written in terms of
polynomial, exponential, trigonometric, hyperbolic, and Mittag-Leffler
functions. Finally, the discussed method is extended for the (k+1)-dimensional
nonlinear time-fractional PDEs with several linear time delays, and the exact
solution of the (3+1)-dimensional nonlinear time-fractional delay
convection-diffusion-reaction equation is derived.
