LAUSR

Abstract A weakly coupled convection dominated system of m-equations is analyzed. A higher order accurate asymptotic-numerical method is presented. The solutions of convection dominated problem are known to exhibit multi-scale character. There exist narrow region across the boundary of the domain where the solution exhibit steep gradient. This region is termed as boundary layer region and the solution of problem is said to have a boundary layer. Outside of this region, the solution of system behaves smoothly. To capture this multi-scale nature given system is factorized into two explicit systems. The degenerate system of initial value problems (IVPs), obtained by setting ϵ = 0, corresponds to the smooth solution, which lies outside of boundary layers. For solution inside boundary layers, a system of boundary value problems (BVPs) is obtained using stretching transformation. Regardless of this simple factorization, solutions of these systems preserve the key features of the given coupled system. Runge–Kutta method is used to solve the degenerate system of IVPs, whereas the system of BVPs is solved analytically. Stability and consistency of the proposed method is established. A uniform convergence of higher order is obtained. Possible extension to differential difference equations are also brought to attention. A comparative study of the present method with some state of art existing numerical schemes is carried out by means of several test problems. The results so obtained demonstrate the effectiveness and potential of present approach.