This paper develops a robust three-level time split high-order Leapfrog/Crank-Nicolson technique for solving the two-dimensional unsteady sobolev and regularized long wave equations arising in fluid mechanics. A deep analysis of the stability and error estimates of the proposed approach is considered using the $L^{\infty}(0,T;H^{2})$-norm. Under a suitable time step requirement, the theoretical studies indicate that the constructed numerical scheme is strongly stable (in the sense of $L^{\infty}(0,T;H^{2})$-norm), temporal second-order accurate and convergence of order $O(h^{\frac{8}{3}})$ in space, where $h$ denotes the grid step. This result suggests that the proposed algorithm is less time consuming, fast and more efficient than a broad range of numerical methods widely discussed in the literature for the considered problem. Numerical experiments confirm the theory and demonstrate the efficiency and utility of the three-level time split high-order formulation.