In this study, numerical simulation of Newtonian fluids in contraction channels is performed based on the Taylor Galerkin-pressure correction finite element method (TGPC-FEM). Here, three different contraction channel geometries (4:1, 6:1, and 8:1) were simulated in this investigation. Usually, the continuity equation and momentum equation are governed this type of fluid. These equations are presented in our study in cylindrical coordinate system. This research examined the impact of varying Reynolds number and geometry area difference on the rate of convergence for solution components, as well as the influence of geometry area difference on level of velocity and pressure.