During the process of numerical approximations, the difference formula of time transient terms, and the central/upwind difference expressions for acquiring boundary value are chosen.