After analyzing the characteristics of different low-order and high-order schemes, it was pointed out that the schemes of second-order numerical accuracy were the optimum way in numerical simulation. Taking, the general two-dimensional steady convection-d