This paper introduces new fully implicit numerical schemes for solving 1D and 2D unsteady Burgers’ equation. The non-linear Burgers’ equation is discretized in the spatial direction by using second order Finite difference method which converts the Burgers’ equation to nonlinear system of ODEs. Then, the backward differentiation formula of order two (BDF-2) is employed to march the solution in the time direction. The non-linear term in the obtained system is linearized without any transformation; and it forms a system of linear algebraic equations that is solved by using Thomas’s algorithm. Accuracy and performance of the proposed schemes are studied by using four test problems with Dirichlet and Neumann boundary conditions. Comparing the numerical results with exact solutions and the solutions of other schemes shows that the proposed schemes are simple, efficient and accurate even for the cases with high Reynolds numbers.