Most physical phenomena are formulated in the form of non-linear fractional partial differential equations to better understand the complexity of these phenomena. This article introduces a recent attractive analytic-numeric approach to investigate the approximate solutions for nonlinear time fractional partial differential equations by means of coupling the Laplace transform operator and the fractional Taylor’s formula. The validity and the applicability of the used method are illustrated via solving nonlinear time-fractional Kolmogorov and Rosenau–Hyman models with appropriate initial data. The approximate series solutions for both models are produced in a rapid convergence McLaurin series based upon the limit of the concept with fewer computations and more accuracy. Graphs in two and three dimensions are drawn to detect the effect of time-Caputo fractional derivatives on the behavior of the obtained results to the aforementioned models. Comparative results point out a more accurate approximation of the proposed method compared with existing methods such as the variational iteration method and the homotopy perturbation method. The obtained outcomes revealed that the proposed approach is a simple, applicable, and convenient scheme for solving and understanding a variety of non-linear physical models.