Poiseuille flows at two Reynolds numbers (Re) 2.5 × 10−2 and 5.0 are simulated by two different smoothed particle hydrodynamics (SPH) schemes on regular and irregular initial particles' distributions. In the first scheme, the viscous stress is calculated directly by the basic SPH particle approximation, while in the second scheme, the viscous stress is calculated by the combination of SPH particle approximation and finite difference method (FDM). The main aims of this paper are (a) investigating the influences of two different schemes on simulations and reducing the numerical instability in simulating Poiseuille flows discovered by other researchers and (b) investigating whether the similar instability exists in other cases and comparing results with the two viscous stress approximations. For Re = 2.5 × 10−2, the simulation with the first scheme becomes instable after the flow approaches to steady-state. However, this instability could be reduced by the second scheme. For Re = 5.0, no instability for two schemes is found.