In the present study, the software of Portable, Extensible Toolkit for Scientific (PETSc) Computation associated with a parallel multilevel iterative solver is embedded in our LES program to enhance the capability of computational power for the three dimensional turbulent flows. The three-dimensional Poiseuille, Couette and lid-driven cavity are conducted for the examination of numerical accuracy at laminar flows. For the moderate Reynolds number, the turbulent Poiseuille flow inside the square duct is simulated for the verification of turbulent model. The results show that the both velocity profiles and turbulence intensities are in line with the benchmark. Moreover, the downstream flow of primary eddy would be enhanced by the fast moving wall. However, as Reynolds number increases higher than 1,000, the location of primary eddy is saturated. Also, the secondary eddies increasing with Reynolds number is confirmed by the present study.