“…Over the last 10 years, various numerical methods have been developed to solve the viscous coupled Burgers' equations. For example, Kapoor and Joshi developed a differential quadrature method with uniform algebraic trigonometric tension B-spline [11], Başhan applied a mixed method with the finite difference and differential quadrature method based on third order modified cubic B-spline functions [4], Hussein and Kashkool used a weak Galerkin finite element method [9], Zhang et al introduced an improved backward substitution method [19], Abdullah et al developed a numerical procedure based on the cubic B-spline and the Hermite formula [1], Mohanty and Sharma presented a high accuracy two-level implicit method based on cubic spline approximation [16], Bhatt and Khaliq modified exponential time-differencing Runge-Kutta method based on Padé approximation [6], Kumar and Pandit proposed a composite scheme based on finite difference and Haar wavelets [12], Jiwari and Alshomrani developed a new collocation method based on modified cubic trigonometric B-spline function [10], Mohanty et al proposed a two-level implicit compact operator method [13]. Among the abovementioned methods, there is a backward semi-Lagrangian (BSLM) to solve the model problem.…”