This article presents the implementation of the k-epsilon model in a code b
ased on the control-volume finite-element method, for the simulation of tur
bulent flow in arbitrary geometries. In order to apply the boundary conditi
ons at the solid wall, equations of the "wall function" are modified to con
sider that the contours are not parallel to the coordinate axis. To validat
e the model, a problem of turbulent flow through an axisymmetric constricte
d tube is solved. The results are compared with experimental and numerical
data reported in the technical literature. The present results are in good
agreement with those obtained using a finite-volume method with a second-or
der upwind scheme, and they show a significant improvement when compared wi
th a finite-volume method with the power law scheme.