Se ha desarrollado un código tridimensional, empleando esquemas de alto orden para evaluar el fenómeno de la turbulencia mediante la implementación del modelo k-épsilon , en una tubería con flujo totalmente desarrollado. La metodología establece la resolución del conjunto de ecuaciones diferenciales usando los esquemas SOU-ULTIMATE y QUICKES-ULTIMATE, a fin de lograr un buen balance entre precisión y aplicabilidad. Los resultados obtenidos con el modelo TURBUL, se compararon contra tres bancos de datos experimentales y contra otras 10 metodologías que emplean el mismo modelo k-épsilon pero ajustado en alguno de sus parámetros. Se concluye que las evaluaciones realizadas sobre la energía cinética turbulenta, la tasa de disipación de energía y la vorticidad a diferentes números de Reynolds han resultado satisfactorias, ofreciendo la ventaja de un menor tiempo computacional