The one dimensional coupled Burgers’ equation can be taken as a simple model of sedimentation and evolution of scaled volume concentrations of two kinds of particles in fluid suspensions and colloids under the effect of gravity. Various researchers have proposed analytical solution to one dimensional coupled Burgers’ equation, e.g. Kaya 4 used Adomian decomposition method, Soliman 5 applied a modified extended tanh-function method, whereas numerical solutions to this system of equation have been attempted by many researchers. Esipov 6 had given numerical solutions and compared the obtained results with those given by the experiment. Abdou 7 used variational iteration method to solve coupled Burgers’ equation, whereas Wei 8 used a conjugate filter approach, Khater 9 applied the Chebyshev spectral collocation method, Dehghan 10 gave numerical solutions of coupled viscous Burgers equations by applying the Adomian- Pade technique; Rashid 11 applied Fourier pseudo-spectral method. Mittal 12 has applied a cubic B-spline collocation scheme while Mokhtari 13 used a generalized differential quadrature method. Recently, Srivastava et al. 1,2 used a fully implicit scheme and Crank-Nicolson scheme for solving this system of coupled Burgers’ equation. Further, Srivastava et al. 14,15 proposed two new finite difference schemes, namely an implicit exponential finite-difference and an implicit logarithmic finite-difference method for solving the two dimensional coupled viscous Burgers’ equation. One can refer 16–21 for various numerical schemes for two dimensional coupled Burgers’ equations whereas the exact solution of two, three and (n + 1)-dimensional Burgers’ equation can be seen in. 22–24
In this article, an implicit logarithmic finite-difference method (I-LFDM) has been applied for the numerical solution of one dimensional coupled Burgers’ equation, proposed by Srivastava et al. 15 The obtained results are compared well with the exact solutions and those already available in the literature. The accuracy and computational reliability of the I-LFDM scheme are demonstrated in terms of error norms by considering the following three test cases.