Turbulent flows through helical pipes with circular cross section are numerically investigated comparing with the experimental results obtained by our team. Numerical calculations are carried out for two helical circular pipes having different pitches and the same non-dimensional curvature δ (= 0.1) over a wide range of the Reynolds number, Re, from 3000 to 21000 for torsion parameter ß (= torsion / = 0.02 and 0.45). We numerically obtained the secondary flow, the axial flow and the intensity of the turbulent kinetic energy by use of three turbulence models incorporated in OpenFOAM. We found that RNG k-ε turbulence model can predict excellently the fully developed turbulent flow with comparison to the experimental data. It is found that the momentum transfer due to turbulence dominates the secondary flow pattern of the turbulent helical pipe flow. It is interesting that torsion effect is more remarkable for turbulent flows than laminar flows.