This study presents, a numerical method for the solutions of the generalized nonlinear Benjamin-Bona-Mahony-Burgers' equation, with variable order local time fractional derivative. This derivative is expressed as a product of two functions, the usual integer order time derivative, and a function of time having a fractional exponent. Then, forward difference approximation is used for time derivative. The unknown solution of the differential problem and corresponding derivatives are estimated using Haar wavelet approximations (HWA). The collocation procedure is then implemented in HWA, to transform the given model to the system of linear algebraic equations for the determination of unknown constant coefficient of the Haar wavelet series, which update the derivatives and the numerical solutions. The sufficient condition is established for the stability of the proposed technique, and then verified computationally. To check the performance of the scheme, few illustrative examples in one and two dimensions along with l∞ and l2 error norms are also given. Besides this, the computational convergence rate is calculated for both type equations. Additionally, computed solutions are compared with available results in literature. Simulations and graphical data discloses, that suggested scheme works well for such complex problems.